表面肌电信号时域频域分析.zip_时域 频域_肌电 Matlab_肌电 MATLAB_肌电信号 频域_肌电信号分析

clc; load('smd.mat', 's'); [hang,lie]=size(s); %% %均值 meanvalue=mean(s); %% %方差 varvalue=0; if lie>1 varvalue=var(s); end %% %绝对值积分平均值 IAV=1.0/lie*sum(abs(s)); %% %均方根值 s2=s'; rms=sqrt(1.0/(lie-1)*(s*s2)); %% %过零点数 num=0; for i=1:lie-1 temp=s(1,i)*s(1,i+1); if temp<0 num=num+1; end end %% %willison幅值 wamp=0; for i=1:lie-1 temp=s(1,i)*s(1,i+1); temp=abs(temp); threshold=50; if temp>threshold wamp=wamp+1; end end