matlab肌电信号hz低通滤波器 燕山大学 数字信号处理 课程设计 肌电信号 Matlab 程序汇总...

a=load('C:\Users\Administrator\Desktop\***\60\3\3.txt');

b=a(1,:);

c=a(2,:);

figure

subplot(2,1,1);

plot(b);

xlabel('时间 ');

ylabel('被测变量');

title('60度第三组原始信号(通道一)');

grid on;

subplot(2,1,2);

plot(c);

xlabel('时间 ');

ylabel('被测变量');

title('60度第三组原始信号(通道二)');

grid on;

fprintf('\n60度第二组n通道一时域参数：\n')

fprintf(' 均值 = %7.5f \n',mean(b)); fprintf(' 标准差 = %7.5f \n', sqrt(var(b)));

fprintf(' 方差 = %7.5f \n', var(b));

fprintf(' 积分肌电值 = %7.5f \n', mean(abs(b)));

fprintf(' 均方根有效值 = %7.5f \n', sqrt(mean(b.^2)) );

fprintf('\n60度第二组通道二时域参数：\n')

fprintf(' 均值 = %7.5f \n',mean(c)); fprintf(' 标准差 = %7.5f \n', sqrt(var(c)));

fprintf(' 方差 = %7.5f \n', var(c));

fprintf(' 积分肌电值 = %7.5f \n', mean(abs(c)));

fprintf(' 均方根有效值= %7.5f \n', sqrt(mean(c.^2)) );

a=load('C:\Users\***_\Desktop\Qu180.txt');

%b=a(1,:);

b=a(2,:);

N=16384;

fs=1000;

y=fft(b,N);%做16384点傅立叶变换

L=length(y);

mag=abs(y);

f=(0:L-1)/L*fs;

figure;

plot(f,mag);%幅频谱

axis([0 1000 0 1000]);