求肌电信号的功率谱分析程序

yatou502

在书上找了一个s周期图法功率谱，可是老师讲效果太差，请哪位大侠帮帮忙，给点提示吧

s=emg';
plot(s)
title('原始肌电信号');
fs=1000;
nfft=1024;
index=0:(nfft/2-1);
k=index*fs/nfft;
xk1=fft(s,nfft);
pxx1=abs(xk1).^2/length(s);
plot_pxx1=10*log10(pxx1(index+1));
cxn1=xcorr(s,'unbiased');
cxk1=fft(cxn1,nfft);
pxx11=abs(cxk1);
plot_pxx11=10*log10(pxx11(index+1));
figure;
plot(k,plot_pxx1);
title('s信号的周期图法功率谱');
figure;
plot(k,plot_pxx11*10);
title('间接法功率谱');

[ 本帖最后由 zhangnan3509 于 2009-6-5 09:18 编辑 ]