现代信号处理 matlab,现代信号处理例题及matlab代码实现

《现代信号处理》期末考核作业

1 MATLAB仿真均值为0，方差为1的白噪声信号，信号长度N=1024，并用周期图法分别求500、1000和1500次实现的平均功率谱密度，画图。

程序代码如下：

clear;

clear all;

N=1024;%数据长度

Nfft=1024;%FFT所采用的数据长度

n=0:N-1;

wn=randn(1,N);%产生随机白噪声

subplot(2,2,1);%绘出白噪声序列

plot(n,wn);

title('白噪声');

%500次实现的平均功率谱密度

s=zeros(1,N);

for i=1:500

wn=randn(1,N);%产生随机白噪声

Pxx=10*log10(abs(fft(wn,Nfft).^2)/N);%Fourier振幅谱平方的平均值，并转换为db

s=s+Pxx;

end

s=s/500;

f=(0:length(Pxx)-1)/length(Pxx);%绘出频率序列

subplot(222);

plot(f,s);

xlabel('频率/Hz');ylabel('功率谱/dB');

title('500次实现的平均功率谱密度');

grid on;