s=?load（‘09121052HF.CSV‘）;
fs=3000;
x1=s（1:163841）;
x2=s（1:163842）;
x3=s（1:163843）;
x4=s（1:163844）;
x5=s（1:163845）;
x6=s（1:163846）;
iemg1=sum（abs（x1））/length（x1）
iemg2=sum（abs（x2））/length（x2）
iemg3=sum（abs（x3））/length（x3）
iemg4=sum（abs（x4））/length（x4）
iemg5=sum（abs（x5））/length（x5）
iemg6=sum（abs（x6））/length（x6）???%求積分肌電值
rms1=sqrt（sum（x1.^2）/length（x1））
rms2=sqrt（sum（x2.^2）/length（x2））
rms3=sqrt（sum（x3.^2）/length（x3））
rms4=sqrt（sum（x4.^2）/length（x4））
rms5=sqrt（sum（x5.^2）/length（x5））
rms6=sqrt（sum（x6.^2）/length（x6））??%求均方根值
L1=length（x1）;
cx1=xcorr（x1‘unbiased‘）;
cxk1=fft（cx1L1）;
px1=abs（cxk1）;%求功率譜密度
pxx1=10*log10（px1）;
f1=（0:L1-1）*fs/L1;
subplot（321）plot（f1（1:L1/2）pxx1（1:L1/2））
xlabel（‘頻率/Hz‘）;ylabel（‘功率譜/dB‘）;
title（‘平均功率譜圖‘）;
grid?on??%做功率譜圖
df1=

?屬性????????????大小?????日期????時(shí)間???名稱(chēng)
-----------?---------??----------?-----??----

?????文件???????4317??2009-12-10?10:16??MFMPF.m

-----------?---------??----------?-----??----

?????????????????4317????????????????????1