tmp605.m

来自「信号分析与处理的一些源代码」· M 代码 · 共 35 行

%MATLAB PROGRAM 6-5
clf
Fs=1000;
%Case 1:  N=256
%==============
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);
f=(0:length(Pxx)-1)*Fs/length(Pxx);
subplot(211)
plot(f,Pxx)
xlabel('Frequency (Hz)');
ylabel('Power Spectrum (dB)');
title('Periodogram N=256')
grid

%Case 2: N=1024
%==============
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('Periodogram N=1024')
grid