%含有噪聲信號的功率譜
clf;Fs=1000;
% 第一種情況: N=256
N=256;Nfft=256;   %數(shù)據(jù)長度和FFT所用的數(shù)據(jù)長度
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);%Fourier振幅譜平方的平均值，并轉(zhuǎn)換為dB
f=(0:length(Pxx)-1)*Fs/length(Pxx); %給出頻率序列
subplot(2,1,1),plot(f,Pxx); %繪制功率譜曲線
xlabel('頻率/Hz');ylabel('功率譜/dB');
title('周期圖N=256');grid on;
%第二種情況: N=1024
Fs=1000;N=1024;Nfft=1024; %采樣頻率、數(shù)據(jù)長度和FFT所用數(shù)據(jù)長度
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);%Fourier振幅譜平方的平均值，并轉(zhuǎn)換為dB
f=(0:length(Pxx)-1)*Fs/length(Pxx); %頻率序列
subplot(2,1,2),plot(f,Pxx);  %繪制功率譜曲線
xlabel('頻率/Hz');ylabel('功率譜/dB');
title('周期圖N=1024');grid on