测试信号处理（第3版）程序代码 第9章.docx

例9.6

%例9.6中周期图法计算信号功率谱的Matlab程序

clf

Fs=2000;

%情况1:数据长度N1=256

N1=256;

N1fft=256;

n1=0:N1-1;

t1=n1/Fs;

f1=50;

f2=100;

xn1=sin(2*pi*f1*t1)+2*sin(2*pi*f2*t1)+randn(1,N1);

Pxx1=10*log10(abs(fft(xn1,N1fft).^2)/N1);

f1=(0:length(Pxx1)-1)*Fs/length(Pxx1);

subplot(2,1,1)

plot(f1,Pxx1)

ylabel(′功率谱(dB)′);

ti