资源描述：
《现代信号处理例题及matlab代码实现》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、《现代信号处理》期末考核作业1MATLAB仿真均值为0，方差为1的白噪声信号，信号长度N=1024，并用周期图法分别求500、1000和1500次实现的平均功率谱密度，画图。程序代码如下：clear;clearall;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);fori=1:50
2、0wn=randn(1,N);%产生随机白噪声Pxx=10*log10(abs(fft(wn,Nfft).^2)/N);%Fourier振幅谱平方的平均值，并转换为dbs=s+Pxx;ends=s/500;f=(0:length(Pxx)-1)/length(Pxx);%绘出频率序列subplot(222);plot(f,s);xlabel('频率/Hz');ylabel('功率谱/dB');title('500次实现的平均功率谱密度');gridon;%1000次实现的平均功率谱密度s=zeros(1,N)
3、;fori=1:1000wn=randn(1,N);%产生随机白噪声Pxx=10*log10(abs(fft(wn,Nfft).^2)/N);%Fourier振幅谱平方的平均值，并转换为dbs=s+Pxx;ends=s/1000;f=(0:length(Pxx)-1)/length(Pxx);%绘出频率序列subplot(223);plot(f,s);xlabel('频率/Hz');ylabel('功率谱/dB');title('1000次实现的平均功率谱密度');gridon;%500次实现的平均功率谱密度
4、s=zeros(1,N);fori=1:1500wn=randn(1,N);%产生随机白噪声Pxx=10*log10(abs(fft(wn,Nfft).^2)/N);%Fourier振幅谱平方的平均值，并转换为dbs=s+Pxx;ends=s/1500;f=(0:length(Pxx)-1)/length(Pxx);%绘出频率序列subplot(224);plot(f,s);xlabel('频率/Hz');ylabel('功率谱/dB');title('1500次实现的平均功率谱密度');gridon;实验结
5、果图如下：2仿真如下随机过程：其中：Vn是均值为0，方差为1的Gaussian白噪声过程，Φ为随机相位，在[0，2π]间服从均匀分布。试对其中的正弦波频率进行估计（在不同的数据长度下，N=16，64，128，1024，可使用经典谱估计中的任何一种方法），并讨论数据长度对估计分辨率和平滑特性的影响。解答：使用周期图法对不同数据长度的信号进行估计。程序代码如下：clearall;%********************N=16;第一种情况数据长度为16**********************N=16;Nff
6、t=16;n=0:N-1;xn=sin(0.5*pi*n+2*pi*rand)+sin(0.3333*pi*n+2*pi*rand)+randn(1,N);figure(1);subplot(2,1,1);plot(n,xn);ylabel('幅值(V)');xlabel('时间(s)');title('原始信号');Pxx=10*log10(abs(fft(xn,Nfft).^2)/N);Fourier振幅谱平方的平均值，并转换为dbf=0:length(Pxx)-1;绘出频率序列subplot(212);
7、plot(f,Pxx);xlabel('频率/Hz');ylabel('功率谱/dB');title('周期图N=16');gridon;%*******************N=64*********************N=64;Nfft=64;n=0:N-1;xn=sin(0.5*pi*n+2*pi*rand)+sin(0.3333*pi*n+2*pi*rand)+randn(1,N);figure(2);subplot(2,1,1);plot(n,xn);ylabel('幅值(V)');xlabe
8、l('时间(s)');title('原始信号');Pxx=10*log10(abs(fft(xn,Nfft).^2)/N);f=0:length(Pxx)-1;subplot(212);plot(f,Pxx);xlabel('频率/Hz');ylabel('功率谱/dB');title('周期图N=64');gridon;%*******************N=64***************
