研究生课程论文封面.docVIP

研究生学位课程试卷 院（系、所） 信息与机电工程学院 专业 通信与信息系统 考试科目 现代数字信号处理及其应用 第 1 学期 研究生姓名 梁超平 学号 132201000 考试成绩 导师评语： 导师签字 年 月 日 随机信号功率谱估计及matlab仿真 姓名： 梁超平 学号：132201000 信息与机电工程学院 信息与通信系统 功率谱估计在实际工程中有重要应用价值，如在语音信号识别、雷达杂波分析、波达方向估计、地震勘探信号处理、水声信号处理、系统辨识中非线性系统识别、物理光学中透镜干涉、流体力学。本文主要讨论经典功率谱估计和现代功率谱估计方法的比较。 经典功率谱估计方法：主要是平均周期图法和分段周期图法，以及welch法。 现代功率谱估计这里主要利用最大熵法（MEM）。 时域波形 周期图法和分段周期图法 %Samp9_7 clf;Fs=1000; %采样频率 N=1024;Nsec=256;n=0:N-1;t=n/Fs; %数据长度，分段数据长度、时间序列 w=hanning(256); %采用的窗口数据 randn(state,0); %设置产生随机数的初始状态 xn=sin(2*pi*50*t)+2*sin(2*pi*120*t)+randn(1,N); %带噪声的信号 %采用不重叠加窗方法的功率谱估计 pxx1=abs(fft(w.*xn(1:256),Nsec).^2)/norm(w)^2; %第一段加窗振幅谱平方 pxx2=abs(fft(w.*xn(257:512),Nsec).^2)/norm(w)^2; %第二段加窗振幅谱平方 pxx3=abs(fft(w.*xn(513:768),Nsec).^2)/norm(w)^2; %第三段加窗振幅谱平方 pxx4=abs(fft(w.*xn(769:1024),Nsec).^2)/norm(w)^2; %第四段加窗振幅谱平方 Pxx=10*log10((pxx1+pxx2+pxx3+pxx4)/4); %求得平均功率谱，转换为dB f=(0:length(Pxx)-1)*Fs/length(Pxx); %求得频率序列 subplot(2,1,1),plot(f(1:Nsec/2),Pxx(1:Nsec/2)); %绘制功率谱曲线 xlabel(频率/Hz);ylabel(功率谱/dB); title(加窗平均周期图(无重叠) N=4*256); grid on %采用重叠加窗方法的功率谱估计 pxx1=abs(fft(w.*xn(1:256),Nsec).^2)/norm(w)^2; %第一段加窗振幅谱平方 pxx2=abs(fft(w.*xn(129:384),Nsec).^2)/norm(w)^2; %第二段加窗振幅谱平方 pxx3=abs(fft(w.*xn(257:512),Nsec).^2)/norm(w)^2; %第三段加窗振幅谱平方 pxx4=abs(fft(w.*xn(385:640),Nsec).^2)/norm(w)^2; %第四段加窗振幅谱平方 pxx5=abs(fft(w.*xn(513:768),Nsec).^2)/norm(w)^2; %第五段加窗振幅谱平方 pxx6=abs(fft(w.*xn(641:896),Nsec).^2)/norm(w)^2; %第六段加窗振幅谱平方