dft法求周期信号功率谱.m

matlab,功率谱估计

clf;
Fs=1000;
N=256;
Nfft=256;
n=0:N-1;
t=n/Fs;
xn=sin(2*pi*50*t)+2*sin(2*pi*120*t)+randn(1,N);
Pxx=10*log10(abs(fft(xn,Nfft).^2)/(N+1));
f=(0:length(Pxx)-1)*Fs/length(Pxx);
subplot(211)
plot(f,Pxx)
xlabel('Frequency(Hz)')
ylabel('Power Spectrum(dB)')
title('N=256')
grid

Fs=1000;
N=1024;
Nfft=1024;
n=0:N-1;
t=n/Fs;
xn=sin(2*pi*50*t)+2*sin(2*pi*120*t)+randn(1,N);
Pxx=10*log10(abs(fft(xn,Nfft).^2)/N);
f=(0:length(Pxx)-1)*Fs/length(Pxx);
subplot(212)
plot(f,Pxx)
xlabel('Frequency(Hz)')
ylabel('Power Spectrum(dB)')
title('N=1024')
grid