[ywf]

做一下频谱分析吧。

x1=hann(N).*x;
x1=x1-mean(x1);
y=fft(x1)/N;
y=abs(y).^2;
y(2:N/2)=2*y(2:N/2);
y(N/2+1:N)=0;

M=N/2-1;
n=1:1:M;
f0=n*(fs/N);
f0=f0/1000;
y=10*log10(y(1:M)/H1);
plot(f0,y);
clear y
grid on;
axis([-0.5 20.5 -120 10]);
set(gca,'XTick',[0:1:20]);
set(gca,'YTick',[-120:10:0]);
xlabel('f (KHz)');ylabel('A (dB)')
title('功率谱')

这个程序可能帮你的忙。