《基于matlab的功率谱分析方法分析》由会员分享，可在线阅读，更多相关《基于matlab的功率谱分析方法分析（35页珍藏版）》请在金锄头文库上搜索。
1、基于matlab的功率谱分析方法研究 作 者 姓 名 专 业 指导教师 专业技术职务 目 录摘 要1第一章 绪论3第二章 谱估计中的变量62.1随机信号简介62.1.1 随机变量62.1.2随机信号的特征72.2平稳随机信号82.2.1平稳随机信号的定义82.2.2 平稳随机信号的自相关函数92.2.3 平稳随机信号的功率谱92.3估计质量的评价标准10第三章 经典功率谱估计123.1谱估计与相关函数123.1.1 相关函数和功率谱123.1.2 相关函数的估计133.2 周期图法153.2.1周期图法的定义153.2.2 周期图的性能153.2.3 周期图法改进措施163.3自相关法173.
2、4 直接法和间接法的关系173.5谱估计仿真与比较173.6 本章小结23第四章 现代谱估计244.1平稳随机信号的参数模型244.2 AR模型的正则方程与参数计算254.2.1 正则方程的求导254.2.2 AR模型参数求解的典型算法264.3 AR模型谱估计的实现及性质274.3.1 谱估计的步骤274.3.2 AR模型谱估计的性质284.3.3 AR 模型阶次p的选择284.3.4 AR模型谱估计仿真294.4 MA模型谱估计304.5 ARMA模型谱估计314.6 小结32第五章 论文总结32参考文献34致 35 / 摘 要数字信号处理重要的应用领域之一,是建立在周期信号和随机信号基础
3、上的功率谱估计。在实际应用中往往不能获得具体信号的表达式,需要根据有限的数据样本来获得较好的谱估计效果,因而谱估计被广泛的应用于各种信号处理中。 本论文研究了功率谱估计的几种常用的方法,包括经典谱估计和现代谱估计的各种方法,且对每种方法的估计质量做了数学推导,并给出仿真程序及仿真图。经典法主要包括周期图法、自相关法,但这两种方法都存在缺陷,即认为观测数据之外的数据都为零,所以对经典法中的周期图法进行了加窗、平均等修正,因此提出了周期图法的改进方法；现代谱估计的方法分类比较多,AR模型法,MA模型法和ARMA模型法是现代功率谱估计中最主要的参数模型,本论文着重讨论了AR模型参数法。同时论文将通过
4、对经典谱估计和现代谱估计的实现方法及仿真图的比较,得出经典功率谱估计方法的方差性较差,分辨率较低,而现代谱估计的目标正是在于努力改善谱估计的分辨率,因此能得到较好的谱估计效果,为此应用更为广泛。关键字：数字信号处理；功率谱估计；周期图法；自相关法；AR模型法。ABSTRACTPerhaps one of the more important application areas of digital signal processing is builting on thePower Spectral Estimation of periodic and random signals. Actua
5、lly, wecant get the expression of a specific signal, so we need to estimate the power spectral of a signal according to some sample data sequences.so spectrum estimation which is widely used in various signal processing.In this thesis, some common methods of Power Spectral Estimation, such as classi
6、cal spectral estimation and modern spectral estimation,are studied. The quality of each estimation method isderived, simulation program and simulation figure is given. Classical methods of Power Spectral Estimation mainly include the Periodogram and the BT method. But both of them have a common draw
7、back: the data sequences, beyond the area of the observed sequences, are all presumed to zero. So the Windows and the average method are introduced to improve the quality of the Periodogram.Therefore the improvement of The Periodogram estimation method is proposed. The classification of modern spect
8、ral estimation methods are more ,AR,MA,andARMAis the most important parametersof modern spectral estimation.This thesis will focus on discussion of AR model parameters method. At the same time ,It can be seen from the comparison and realization of classical spectral estimation and modern spectral es
9、timation,classical power spectrum estimation variance is poor, low resolution .The goal of modern spectral estimation is woking to improve the resolution of spectral estimation, better results of the estimation of the power spectrum can be obtained, so it is applied more widely. Keywords: digital si
10、gnal processing;Power Spectrum Estimation;The Periodogram；the BT methods；AR model.第一章 绪论1.功率谱估计的发展 功率谱估计技术渊源流长,在过去的几十年获得了飞速的发展。功率谱估计涉及信号与系统、随机信号分析、概率统计、随机过程、矩阵代数等一系列的基础科学,广泛应用于雷达、声纳、通信、地址勘探、天文、生物医学工程等众多领域,其容、方法不断更新,是一个具有强大生命力的研究领域。 功率谱估计是用有限长的数据来估计信号的功率谱, 它对于认识一个随机信号或其它应用方面来讲都是极其重要的, 是数字信号处理的重要研究容之一
11、,在军事、生物医学、通信等领域得到了较为广泛的应用3。谱最早是由英国科学家牛顿提出来的,后来法国工程师傅里叶提出了著名的傅里叶谐波分析理论,该理论至今仍然是我们进行信号分析和处理的理论基础。傅里叶级数首先在观察自然界中的周期现象得到应用,但傅里叶的计算比较复杂,促使人们研制相应的机器来计算傅里叶级数。在19世纪末,Schuster提出傅里叶系数的平方,并命名为周期图,这是经典谱估计的最早提出法,至今仍被人们沿用。后来,鉴于周期图的起伏剧烈,提出了平均周期图的概念,并提出了在对有限长数据计算傅里叶系数时所存在的边瓣问题,这就是后来我们所熟悉的窗函数的影响。周期图较差的方差性能促使人们研究另外的分
12、析方法。Yule在1927年提出了用线性回归方程来模拟一个时间序列,从而发现隐含在该时间序列中的周期,从而发现了现代谱估计中最重要的方法参数模型法。Walker利用Yule的分析方法研究了衰减正弦时间序列,并得出了在对最小二乘分析中经常应用的Yule-Walker方程。Yule的工作使人们重新想起了早在1795年Prony提出的指数拟合法,从而Prony方法形成了现代谱估计的又一重要容。之后又陆续提出了Wiener-khintchine定理、谱估计自相关法BT法等。所有这些都为现代谱估计的发展打下了良好的基础1。2.功率谱估计的方法 功率谱估计可以分为经典谱估计 和现代谱估计。经典谱估计的方法
13、主要方法有自相关估计法和周期图法以及对周期图的改进方法; 现代谱估计的容极其丰富,涉及的学科及应用领域也相当广泛,方法大致可分为参数模型谱估计和非参数模型谱估计,前者有AR 模型法、MA模型,ARMA模型、Prony 指数模型等；后者有最小方差法,多分量的MUSIC方法等。其中周期图法和AR 模型法是用得较多且最具代表性的方法。从信号的来源分,又可分为一维谱估计、二维谱估计及多维谱估计。从使用的统计量来分,目前大部分工作是建立在二阶矩基础上的,但由于功率谱密度是频率的实函数,缺少相位信息,因此,建立在高阶矩基础上的谱估计方引起人们的注意。从信号的特征来分,在这之前所说的方法都是对平稳随机信号而
14、言,其谱分量不随时间变化,对非平稳随机信号,其谱是时变的,近20年来,以wigner分析为代表的时域分析引起了人们的广泛兴趣,形成了现代谱估计的一个新的研究领域1。3.功率谱估计的提出 在通信系统中,往往需要研究具有统计特性的随机信号。由于随机信号是一持续时间无限长,具有无限大能量的功率信号,它不满足傅里叶变换条件,而且也不存在解析表达式,因此就不能够应用确定信号的频谱计算方法去分析随机信号的频谱。然而,虽然随机信号的频谱不存在,但其相关函数是可以确定的。如果随机信号是平稳的,那么其相关函数的傅里叶变换就是它的功率谱密度函数,简称功率谱。功率谱反应了单位频带随机信号的一个样本信号来对该随机过程
15、的功率谱密度函数做出估计。4.经典谱估计 直接法：又称周期图法,它是把随机序列x的N个观测数据视为一能量有限的序列,直接计算x的离散傅立叶变换,得X,然后再取其幅值的平方,并除以N,作为序列x真实功率谱的估计。周期图这一概念早在1989年就提出了,但由于点数N一般比较大,该方法的计算量过大在当时无法使用,在1965年FFT出现后,此法才变成谱估计的一个常用方法。周期图法包含了下列二条假设：1. 认为随机序列是广义平稳且各态遍历的,可以用其一个样本x中的一段来估计该随机序列的功率谱,这当然必带来误差。2. 由于对有限序列采用DFT,就默认此有限序列时域是周期的,以及该有限序列在频域是周期的。这种方法把随机序列样本x看成是截得一段的有限序列的周期延拓,这也就是周期图法这个名字的来历。 间接法：也叫相关法。间接法先由序列x估计出自相关函数Rn
