samp9_5.m

MATLAB7.x数字信号处理 光盘内容

%Samp9_5
clf;Fs=1000;
% 第一种情况: N=256
N=256;Nfft=256;   %数据长度和FFT所用的数据长度
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振幅谱平方的平均值，并转换为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; %采样频率、数据长度和FFT所用数据长度
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振幅谱平方的平均值，并转换为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