第五章功率谱估计1-2节

第五章

功率谱估计1第一节引言21、信号与系统的分析研究方法（1）在时域进行（一维信号）（2）在频域进行（一维信号）（3）在时频域进行（二维信号）（4）在空域进行（多维信号）3（1）在时域进行如我们前面所学的维纳滤波、卡尔曼滤波和自适应滤波都属于这种方法。4（2）在频域进行对确定性信号傅里叶变换是在频率域分析研究的理论基础。对随机性信号其傅里叶变换并不存在，因此研究它的功率谱。5（3）在时频域进行（二维信号）时频分析是一种新的信号分析方法，对于非平稳信号、宽带信号，采用时频分析方法显得特别重要，近年来得到越来越广泛的重视。由于时频特征分离的固有性质，使局限时间或频谱的一维信号处理扩展到二维的时频平面，信号的时频特征得到了分离，对信号的分析可以具体到信号的特定时间和特定频率，全面反映观测信号的时频联合特征，相应的时频滤波也变得更为容易、更为直观。时频分布可分为线性时频表示（如短时傅立叶变换、小波分析等）和非线性时频表示。6（4）在空域进行空域信号处理器是雷达必不可少的最重要分系统。阵列信号处理。属于空域信号处理。空域相关去噪法。是一种新的自适应小波滤噪法。极化域——空域联合谱估计。分布式移动通信系统。充分利用空域来提高系统容量及频谱利用率。空时二维自适应信号处理是新一代机载雷达的一项关键技术。72、随机信号的频域分析随机信号是无始无终，能量无限的信号，其傅里叶变换不收敛，因而不能确定此类信号的频谱。随机信号能量无限，其功率未必无限，因而常用功率谱来描述其频率特性。随机信号自相关函数的傅里叶变换是信号的功率谱密度。83、谱估计定义谱估计或功率谱估值：根据有限个观测数据，估计平稳随机信号的功率谱。93、谱估计的应用谱估计的应用遍及雷达、声纳、通信、生物医学、地震勘探等诸多领域。在雷达信号处理中，由回波信号功率谱密度、谱峰的密度、高度和位置，可以确定运动目标的位置、幅射强度和运动速度等。在被动式声纳信号处理中，谱峰的位置可提供鱼雷的方向（方位角）。10在生物医学工程中，有关生理电信号的功率谱密度的谱峰可以指示癫痫病的发作周期。在电子战中，谱分析可用来对目标进行分类，识别等。根据信号、干扰与噪声的功率谱，可设计适当的滤波器，尽量不失真地重现信号，最大限度地抑制干扰与噪声。11功率谱估计分为两大类：（1）经典谱估计（非参数化方法）（线性谱估计）。（2）现代谱估计（参数化方法）（非线性谱估计）。4、谱估计的分类12实质上仍依赖于传统的傅里叶变换法。经典的谱估计法又分为两种：（a）相关图法（又称间接法（BT法））。是由布莱克曼(R.B.Blackman)和图基(J.W.Tukey)于1958年提出的，称为BT法。（b）直接法（又称周期图periodogram法）。是由舒斯特(Schuster)于1898年提出的。（1）经典谱估计

（非参数化方法）13（a）间接法（BT法）直到1965年快速傅里叶变换算法（FFT）问世以前，是最流行的谱估计方法。BT法又称为相关图法对信号序列估计求其自相关函数值对自相关函数的估计进行加权对加权的自相关函数做傅里叶变换获得功率谱估计。14（b）直接法（又称周期图（periodogram)法）对观测到的数据样本直接进行傅里叶变换取模的平方，再除以N得到功率谱估计。不用估计自相关函数，且可以用FFT进行计算，在FFT出现以后，周期图法才得到了广泛的应用。15（2）现代谱估计其基本思想是根据已有的观测数据，建立信号所服从的模型，从而在观测不到的区间上，信号的取值服从模型的分布情况，不再认为是零。主要讨论参数模型（AR、MA、ARMA）法。165、经典谱估计与现代谱估计

方法的比较(1)经典谱估计的优缺点(2)现代谱估计的优缺点17(1)经典谱估计的优缺点优点：原理简单，便于实现，并可采用FFT等技术而使计算量大为减小。缺点：它的估计方差大，谱分辨率差（分辨率大约为数据长度的倒数）。使得这种方法难以应用于短数据记录等情况。原因：一个假设，即除了得到N个数据外，序列的其它值均被认为是零。但实际上序列或其自相关函数在我们未能观测到或未估计的值，实际上并不全是零。18(2)现代谱估计的优缺点优点：现代谱估计质量比经典谱估计的质量有很大提高。缺点：没有理论能指导选择一个合适的模型，只能根据功率谱的一些先验知识，即重要的谱特性，选择模型。选型合适，可以得到偏差小分辨率高、方差小的谱估计质量。196、功率谱估计的发展功率谱估计的分辨能力用参量法可以改进。如自回归模型法、最大熵法和最大似然估计等。在这些方法中，不再认为观测到的N个数据以外的数据全为零，因此克服了经典法的这个缺点，提高了谱估计的分辨率。20然而这些方法在信噪比（SNR）较低时性能并不好，为此，1982年以来，人们陆续提出了多种基于矩阵奇异值分解或特征值分解的改进的谱估计方法，也叫做超分辨方法。21第二节

经典谱估计方法22一、相关图法2324二、相关图法功率谱估计步骤2526272829三、相关图法功率谱估计质量303132333435363738394041424344二、周期图法1、周期图法定义2、周期图与BT法的等价关系3、周期图法的谱估计质量4、周期图法的改进451、周期图法定义46474849周期图法计算功率谱框图502、周期图与BT法的等价关系51523、周期图法的谱估计质量5354555657585960616263646566例子67图(a)68程序N=64;Nfft=64;sigma2=1;wn=sigma2.*randn(1,N);W=fft(wn,Nfft);Pxx=(abs(W).^2)/Nfft;f=(0:length(Pxx)-1)/length(Pxx)*2*pi;plot(f(1:Nfft/2),Pxx(1:Nfft/2))69图(b)70图(c)71图(d)72结论733、周期图法的改进周期图法作为功率谱的估计不满足一致估计的条件，必须进行改进。改进的措施主要是将周期图进行平滑（平均是一种主要的平滑方法），使估计方差减小，从而得到一致谱估计。平滑的方法主要有两种：（1）窗口处理法。目的平滑周期图（2）周期图平均法（3）Welch法：修正的周期图平均法周期法从原理上讲，是基于相关函数值取傅里叶变换的一种方法。74（1）窗口处理法75（a）窗口处理法原理7677787980（b)几个谱窗（1）矩形滞后窗81（2）巴特利特窗（Bartlett’sWindow)82（3）布莱克曼－图基窗（Blackman-TukeyWindow)83（4）帕曾窗（ParzenWindow)84结论85（2）周期图平均法原理：对一随机变量进行观测，得到L组独立记录数据。用每一组数据求其平均值，然后将L个均值加起来求平均。或是对一随机变量进行观测，得到一组独立记录数据，先将数据分L段，对各段数据分别求其周期图，再将各段周期图进行平均作为功率谱的估计。又称为Bartlett方法，由于其采用FFT运算可提高运算速度，是当前用得最多的一种平滑方法。8687888990919293949596例子97图a98图b99图c100101(3)Welch法目前在工程实践中，经典谱估计获得更多有效应用的是1967年由Welch提出的修正周期图法，他将前