数字信号处理课件胡广书第11章.pptVIP

第十一章 经典谱估计;请抓住并搞清楚如下四个问题： 功率谱为什么要估计？ 如何估计？ 如何评价估计质量？ 如不理想，如何改进？;平稳随机信号功率谱的两个定义：;平稳信号;1. 求极限：;双求和变成单求和：;功率谱的两个定义都要求：样本无穷多，时间无限长，即需要集总平均。;;11.2 自相关函数估计; 在实际应用中，我们所遇到的大都是实际物理信号，因此是因果性的，即当n0时，x(n)≡0，且x(n)是实信号，这样，其自相关函数r(m)可表示为 r(m)的估计方法通常有两种，一是利用公式直接计算，二是先计算出 的能量谱，然后对该能量谱作反变换。; 如果观察值的点数N为有限值，则求r(m)估计值的一种方法是 由于x(n)只有N个观察值，因此，对于每一个固定的延迟m，可以利用的数据只有N-1-|m|个，且在0~N-1的范围内， ＝x(n),所以在实际计算时，上式变为 它的长度为2N-1，是以m=0为偶对称的。;实际求出的自相关函数;偏差;所以：;在数据上加矩形窗，长度为 N ，该矩形窗函数的自相关函数正是三角窗！注意矩形窗加在数据上，三角窗加在相关函数上，体现在估计的自相关函数的均值上。;方差;由：;3.自相关函数的计算;（2） 利用FFT：;自相关函数的另一个估计方法（估计子）：;11.3 经典谱估计;;1.周期图（Periodogram）法：;2.自相关（Blackman-Tukey BT法）法:;3.直接法和间接法的关系：; （一）;对 补N 个零，做DFT,得到;不补零，; （二）;直接法和间接法之间的关系;直接法;11.4经典谱估计的质量;1. 偏差;于是有：;最后有：;所以：;;有关方差公式的推导不作要求。主要是掌握结论，并用来说明问题。;推导的结果：方差;解释：;推导的结果：协方差;;在 ; 增大， 的主瓣（ ）将变窄，因此，引起不相关的区域进一步增多，从而引起谱曲线的更加起伏，实际上是方差变大。;对白噪声在不同长度情况下估计出的谱曲线：;经典谱估计质量的讨论：;;假定1： 是慢变谱，在 的主瓣内近似为一个常数;2.方差; ：方差改进之比 ;1. 在 加上 后，估计的谱 的偏差劣于 M=N－1 时估计的谱，而方差优于 M=N－1 时估计的谱；;原主瓣宽，取决于;间接法与直接法缺点;直接法缺点： 周期图法应用比较广泛，主要是由于它与序列的频谱有直接的对应关系，并且可以采用FFT快速算法来计算。但是，这种方法需要对无限长的平稳随机序列进行截断，相当于对其加矩形窗，使之成为有限长数据。同时，这也意味着对自相关函数加三角窗，使功率谱与窗函数卷积，从而产生频谱泄漏，容易使弱信号的主瓣被强信号的旁瓣所淹没，造成频谱的模糊和失真，使得谱分辨率较低。;; 理论依据: L个独立同分布随即变量和的分布，方差减小 倍，即：;（1） Bartlett平均;每一段谱;： 的自相关； 类似 引入的; 的均值为 上式中 是矩形窗 的频谱， 是由 做自相关所得到的三角窗 的频谱， 的长度是2M-1。可见，不取平均的周期图和取平均后的都是有偏估计，且当N→∞时，二者都是渐近无偏的。但因为 主瓣的宽度远大于 ，所以取平均后，偏差加大，分辨率下降。;如果x(n)为一白噪声序列，有 因此，分的段数越多，方差越小。如若L能趋于∞，则 是P(ω)的一致估计。由上面的分析我们再一次看到，方差性能的改善是以牺牲偏差和分辨率为代价的。 每段数据长度M的选择主要取决于所需的分辨率。因为 主瓣的宽度是4π/M,若P(ω)中有两个相距为BW的谱峰，为了要分辨它们，需要4π/MBW,即M4π/BW。如果数据长度N已确定，根据所需的M，段数L也就自然被确定。如果N可以变化，则应根据方差要求确定L，然后再确定要记录的数据长度N。;Bartlett法;（2） Welch平均;：不一定是矩形窗，如Hamming窗;Welch平均法的方差比Barttlett方法有明显的减小，而偏差几乎没有减小;;;H(z);令：;调整 ，可以得到不同的信噪比，本例取;;（e）BT法，M＝32；（f）BT法，M＝