mpf与mf算法matlab实现.docVIP

mpf与mf算法matlab实现 s= loadF.CSV); fs=3000; x1=s(1:16384,1); x2=s(1:16384,2); x3=s(1:16384,3); x4=s(1:16384,4); x5=s(1:16384,5); x6=s(1:16384,6); 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(cx1,L1); px1=abs(cxk1);%求功率谱密度 pxx1=10*log10(px1); f1=(0:L1-1)*fs/L1; subplot(3,2,1),plot(f1(1:L1/2),pxx1(1:L1/2)) xlabel(频率/Hz);ylabel(功率谱/dB); title(平均功率谱图); grid on %做功率谱图 df1=fs/L1; p1=(sum(px1(1:L1/2-1))+sum(px1(1:L1/2)))/2.*df1; pf1=(sum(px1(1:L1/2-1).*[1:L1/2-1].*df1)+sum(px1(1:L1/2).*[1:L1/2].*df1))/2*df1; MPF1=pf1/p1 %求平均功率频率 N1=1;pp1=0; while abs(pp1-p1/2)(px1(N1)+px1(N1+1))/2*df1 pp1=pp1+(px1(N1)+px1(N1+1))/2*df1; N1=N1+1; end n_1=(N1+N1+1)/2; MF1=df1*n_1 %求中值频率 L2=length(x2); cx2=xcorr(x2,unbiased); cxk2=fft(cx2,L2); px2=abs(cxk2);%求功率谱密度 pxx2=10*log10(px2); f2=(0:L2-1)*fs/L2; subplot(3,2,2),plot(f2(1:L2/2),pxx2(1:L2/2)) xlabel(频率/Hz);ylabel(功率谱/dB); title(平均功率谱图); grid on %做功率谱图 df2=fs/L2; p2=(sum(px2(1:L2/2-1))+sum(px2(1:L2/2)))/2.*df2; pf2=(sum(px2(1:L2/2-1).*[1:L2/2-1].*df2)+sum(px2(1:L2/2).*[1:L2/2].*df2))/2*df2; MPF2=pf2/p2 %求平均功率频率 N2=1;pp2=0; while abs(pp2-p2/2)(px2(N2)+px2(N2+1))/2*df2 pp2=pp2+(px2(N2)+px2(N2+1))/2*df2; N2=N2+1; end n_2=(N2+N2+1)/2; MF2=df2*n_2 %求中值频率 L3=length(x3); cx3=xcorr(x3,unbiased); cxk3=fft(cx3,L3); px3=abs(cxk3);%求功率谱密度 pxx3=10*log10(px3); f3=(0:L3-1)*fs/L3; subplot(3,2,3),plot(f3(1:L3/2),pxx3(1:L3/2)) xlabel(频率/Hz);ylabel(功率谱/dB); title(平均功率谱图); grid on %做功率谱图 df3=fs/L3; p3=(sum(px3(1:L3/2-1))+sum(px3(1:L3/2)))/2.*df3; pf3=(sum(px3(1:L3/2-1).*[1:L3/2-1].*df3)+sum(px3(1:L3/2).*[1:L3/2].*df3))/2*df3; MPF3=pf3/p3 %求平均