现代信号处理例题和matlab代码实现.docVIP

《现代信号处理》期末考核作业 1 MATLAB仿真均值为0，方差为1500、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; %1000次实现的平均功率谱密度 s=zeros(1,N); for i=1:1000 wn=randn(1,N);%产生随机白噪声 Pxx=10*log10(abs(fft(wn,Nfft).^2)/N);%Fourier振幅谱平方的平均值，并转换为db s=s+Px