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]);

總結

以上是生活随笔為你收集整理的matlab肌电信号hz低通滤波器,燕山大学 数字信号处理 课程设计 肌电信号 Matlab 程序汇总...的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。