matlab功率谱分析.pdfVIP

功率谱估计及其 MATLAB 仿真 王凤瑛，张丽丽 ( 山东科技大学 信息与电气工程学院，山东 青岛 266510) 摘要：从介绍功率谱的估计原理入手分析了经典谱估计和现代谱估计两类估计方法的原理、各自 特点及在 Matlab 中的实现方法。 关键词：功率谱估计；周期图法；AR 参数法；Matlab 中图分类号：TP391.9 文献标识码：A Power Spectrum Density Estimation and the simulation in Matlab WANG Feng-ying ,Zhang Lili (College of Info and Elect. Eng., SUST, Qingdao, Shandong 266510, China) Abstract:Mainly introduces the principles of classical PSD estimation and modern PSD estimation, discusses the characteristics of the methods of realization in Matlab.Moreover,It gives an example of each part in realization using Matlab functions. Keywords:PSD Pstimation, Periodogram method, AR Parameter method, Matlab 现代信号分析中，对于常见的具有各态历经的平稳随机信号，不可能用清楚的数学关系式来描 述，但可以利用给定的 N个样本数据估计一个平稳随机信号的功率谱密度叫做功率谱估计(PSD) 。它 是数字信号处理的重要研究内容之一。功率谱估计可以分为经典功率谱估计(非参数估计)和现代功率 谱估计(参数估计) 。 功率谱估计在实际工程中有重要应用价值，如在语音信号识别、雷达杂波分析、波达方向估计、 地震勘探信号处理、水声信号处理、系统辨识中非线性系统识别、物理光学中透镜干涉、流体力学 的内波分析、太阳黑子活动周期研究等许多领域，发挥了重要作用。 Matlab 是 MathWorks 公司于 1982 年推出的一套高性能的数值计算和可视化软件，人称矩阵实 验室，它集数值分析、矩阵运算、信号处理和图形显示于一体，构成了一个方便的、界面友好的用 户环境，成为目前极为流行的工程数学分析软件。也为数字信号处理进行理论学习、工程设计分析 提供了相当便捷的途径。本文的仿真实验中，全部在 Matlab6.5 环境下调试通过；随机序列由频率 不同的正弦信号加高斯白噪声组成。 1 经典功率谱估计 经典功率谱估计是将数据工作区外的未知数据假设为零，相当于数据加窗。经典功率谱估计方 法分为：相关函数法（BT 法）、周期图法以及两种改进的周期图估计法即平均周期图法和平滑平均 周期图法，其中周期图法应用较多，具有代表性。 1．1 相关函数法（BT 法） 该方法先由序列 x(n)估计出自相关函数 R(n)，然后对 R(n)进行傅立叶变换，便得到 x(n)的功 率谱估计。当延迟与数据长度相比很小时，可以有良好的估计精度。 Matlab 代码示例 1： Fs=500;%采样频率 n=0:1/Fs:1; %产生含有噪声的序列 xn=cos(2*pi*40*n)+3*cos(2*pi*90*n)+randn(size(n)); nfft=512; cxn=xcorr(xn,unbiased); %计算序列的自相关函数 CXk=fft(cxn,nfft); Pxx=abs(CXk); 1 index=0:round(nfft/2-1); k=index*Fs/nfft; plot_Pxx=10*log10(Pxx(index+1)); figure(1)% plot