经典功率谱估计课件.ppt

《经典功率谱估计课件.ppt》由会员分享，可在线阅读，更多相关《经典功率谱估计课件.ppt（74页珍藏版）》请在三一办公上搜索。

1、第十一章 经典谱估计,11.1 概述11.2 自相关函数的估计11.3 经典谱估计的基本方法11.4 经典谱估计的质量11.5 经典谱估计的改进11.6 经典谱估计算法比较11.7 短时傅里叶变换,第十一章 经典谱估计11.1 概述,请抓住并搞清楚如下四个问题：功率谱为什么要估计？如何估计？如何评价估计质量？如不理想，如何改进？,11.1 概述,请抓住并搞清楚如下四个问题：11.1 概述,平稳随机信号功率谱的两个定义：,随机信号的单个样本,集总平均,平稳随机信号功率谱的两个定义：随机信号的单个样本求均值运算求,平稳信号,单一样本,可将 看作能量信号，因此，可对它作傅立叶变换，并得到功率谱：,问

2、题：的功率谱 和单个样本的功率谱 有何关系？和整个随机信号的功率谱 有何关系？,平稳信号单一样本可将 看作能量信号，,1.求极限：,2.求均值：,1.求极限：2.求均值：单一样本的功率谱不能收敛到所有样,双求和变成单求和：,双求和变成单求和：证明了两个公式等效。所以自相关函数是集总自,功率谱的两个定义都要求：样本无穷多，时间无限长，即需要集总平均。,功率谱估计：古老而又年轻的话题！,实际工作中，我们往往能得到的是：1.单一的样本；2.单一样本的有限长数据；,功率谱的两个定义都要求：样本无穷多，时间无限长，即需要集总平,经典功率谱估计课件,11.2 自相关函数估计,目的：自身估计的需要；功率谱估

3、计的需要,定义：,11.2 自相关函数估计目的：自身估计的需要；集总自相关,近似质量如何,？,Estimation Estimate Estimator（估计子）,估计方法：,从估计方法上看，实际上是把随机信号“视为”单样本有限长的确定性信号。问题是：,实际求出的自相关函数近似质量如何？Estimation 估计,偏差,自相关函数估计的质量：,估计方法,单个样本,1.偏差,来自定义,所有样本,偏差自相关函数估计的质量：估计方法单个样本1.偏差来自定义,所以：,所以：含义？渐近无偏估计对固定的N，此结论给出了m的选取原则,在数据上加矩形窗，长度为 N，该矩形窗函数的自相关函数正是三角窗！注意矩形

4、窗加在数据上，三角窗加在相关函数上，体现在估计的自相关函数的均值上。,那儿来的三角窗？,在数据上加矩形窗，长度为 N，该矩形窗函数的自相关函数正是,方差,2.方差,来自定义,包含两项,前面结果,四阶统计量！,方差2.方差来自定义包含两项前面结果四阶统计量！,由：,最后导出：,零均值高斯分布,由：最后导出：有：渐近一致估计零均值高斯分布,3.自相关函数的计算,已知单个样本的 N 点数据,估计,两个方法:（1）直接按定义：,3.自相关函数的计算已知单个样本的 N 点数据估计两个方法:,（2）利用FFT：,Step1:将 补 个零得;,Step2:对 做FFT,得；,Step3:对 求幅平方，得；,

5、Step4:由 得，对其作IFFT,得。,（2）利用FFT：Step1:将 补 个,自相关函数的另一个估计方法（估计子）：,很容易证明：是 的无偏估计，但方差性能不好。在一些谱估计的方法中，有时用到该公式。,要求：很好掌握自相关函数的估计方法及估计性质。,自相关函数的另一个估计方法（估计子）：很容易证明：,11.3 经典谱估计,问题的提出：对随机信号，我们往往只能得到它的：1.单一的样本；并且仅是 2.单一样本的有限长数据；如何用这 N 数据去估计原随机信号真实的功率谱,？,11.3 经典谱估计问题的提出：对随机信号,1.周期图（Periodogram）法：,经典谱估计中有两个基本的方法：,思

6、路：对 做DTFT(DFT),得到频谱；对该频谱求幅平方，再除以N，即得到“周期图”功率谱，以此作为对真谱的估计。,1.周期图（Periodogram）法：经典谱估计中有两个基,2.自相关（Blackman-Tukey BT法）法:,Step1,Step2,因为先要估计自相关函数，所以又称间接法。与此相对应，周期图法又称直接法。,2.自相关（Blackman-Tukey BT法）法:Ste,3.直接法和间接法的关系：,需要考虑两种情况：,（一）,（二）,数据的范围,自相关函数的范围,3.直接法和间接法的关系：需要考虑两种情况：（一）（二）,（一）,比较用两种方法的估计出的离散谱：,2N 点的谱

7、，把所能估计出的自相关函数都使用上了，而估计自相关函数时，把 N 点数据也全都使用上了。,（一）比较用两种方法的估计出的离散谱：2N 点的谱，把所能,对 补N 个零，做DFT,得到,结论：在 时，直接法和间接法估计的结果是一样的。,？,使用间接法时，往往取，这时二者是不一样的。因此，直接法可看作是间接法的特例。,对 补N 个零，做DFT,得到IFFT,不补零，,思考：,即：,如何和 相等？,不补零，思考：即：N点离散谱如何和,（二）,所以：,相当于只用了部分自相关函数,（二）所以：加在自相关函数上。目的是将其截短。第二次加窗。,直接法和间接法之间的关系,直接法和间接法之间的关系,11.4经典谱

8、估计的质量,11.4经典谱估计的质量也分两种情况讨论主要考察的是均值无偏,1.偏差,估计值的均值,1.偏差估计值的均值自相关函数估计的性质,于是有：,的真实功率谱；,的频谱；,的频谱；,三角窗；,注意：三角窗频谱恒为正,于是有：的真实功率谱；的频谱；的频谱；三角窗；注意：三角窗,最后有：,由于,最后有：由于如何理解这一结果？,所以：,因为：,所以：周期图和自相关法都是渐近无偏估计因为：,2.方差,又遇到四阶矩问题，直接求解困难。,（2）求 在 处的协方差：,定义：,2.方差又遇到四阶矩问题，直接求解困难。（1）假定,有关方差公式的推导不作要求。主要是掌握结论，并用来说明问题。,（3）令，则,求

9、解的关键,有关方差公式的推导不作要求。主要是掌握结论，并用来说明问题。,推导的结果：方差,（1）时,经典功率谱估计不是一致估计,推导的结果：方差（1）时经典功率谱估,解释：,解释：,推导的结果：协方差,推导的结果：协方差 假定 在主瓣外为零；,经典功率谱估计课件,在,处，,说明：随机变量 在 处不相关;,原因：功率谱的定义中即要求极限，又要求均 值；而实际的估计方法，仅靠单次实现 的有限长，无极限、又无均值运算，因 此产生上述问题。,后果：使估计出的谱曲线起伏加剧；,在 处，说明：随,增大，的主瓣（）将变窄，因此，引起不相关的区域进一步增多，从而引起谱曲线的更加起伏，实际上是方差变大。,分辨率

10、和方差（体现在曲线起伏上），是经典谱估计中的一对矛盾。,通常，增加，会提高谱的分辨率，对经典谱估计来说，增加 固然会有利于提高分辨率，但谱曲线的起伏令使用者难以接受，这是经典谱估计的一个致命缺点。,增大，的主瓣（,对白噪声在不同长度情况下估计出的谱曲线：,N16 N 32,N64 N 128,对白噪声在不同长度情况下估计出的谱曲线：N16,经典谱估计质量的讨论：,(二)、,：加在估计的自相关函数上，,周期图谱估计和自相关法的谱估计 不再一样！,经典谱估计质量的讨论：(二)、：加在估计的自相关函数上,1.偏差,1.偏差谁的主瓣比较宽?,假定1：是慢变谱，在 的主瓣内近似为一个常数,假定1：是慢变

11、谱，在,2.方差：,考虑特殊情况，为白噪序列，其功率谱应为常数，即,时对白噪声功率谱估计的方差,2.方差：考虑特殊情况，,：方差改进之比,两种情况下估计的方差之比：,?,取哈明窗：,：方差改进之比 两种情况下估计的方差之比：,1.在 加上 后，估计的谱 的偏差劣于 M=N1 时估计的谱，而方差优于 M=N1 时估计的谱；,（2）在 的范围上，因为B变大，不相关的点变少。,2.上加窗 以后，估计谱 方差的改进体现在两个方面：,（1）,结论：,1.在 加上 后，估计的,原主瓣宽，取决于,现主瓣宽，取决于,3.方差的减小是以牺牲分辨率为代价的！,原主瓣宽，取决于现主瓣宽，取决于3.方差的减小是以牺牲

12、分辨,11.5 直接法估计的改进,任务：改进 对 估计的性能；,目标：主要是改进方差的性能,方法：平滑与平均；,11.5 直接法估计的改进任务：改进,理论依据:L个独立同分布随即变量和的分布，方差减小 倍，即：,将一个较长的信号分成若干段，对每一段求功率谱，每一段的功率谱都是随机变量，然后平均之。类似相干平均，用以弥补经典谱估计中缺少的求均值运算。注意：信号应是平稳的，且每一段的统计特性基本一样。,2.平均（Average）,理论依据:L个独立同分布随即变量和的分布，,（1）.Bartlett平均,将 分成 段，每段 点，即,（1）.Bartlett平均将 分成,每一段谱平均后谱平均后估计出的

13、功率谱的性能如何？在数据上加了,：的自相关；类似 引入的,：的自相关；统计性能分析：？,2.Welch平均,特点：交叠分段,2.Welch平均特点：交叠分段若重叠一半，段数,：不一定是矩形窗，如Hamming窗,Welch 平均是常用的经典谱估计方法，MATLAB中有相应的命令,：不一定是矩形窗，如Hamming窗归一化因子，保证无偏估计,Welch平均法的方差比Barttlett方法有明显的减小，而偏差几乎没有减小,3.Nottall 法：平滑与平均相结合,假定1：是慢变谱，在 的主瓣内近似为一个常数）,假定2,Welch平均法的方差比Barttlett方法有明显的减小，,经典功率谱估计课件

14、,11.6 总结与比较,请掌握如下的方法：,白噪声1,白噪声2,两个输出都是随机信号,由自己指定,H(z)H(z)11.6 总结与比较请掌握如下的方法：白噪声,令：,则：,构成一复信号,得到 的功率谱；,在 的基础上再加上四个复正弦，归一化频率分别是：,令：则：构成一复信号得到 的功率谱；在,调整，可以得到不同的信噪比，本例取,这样，的真实功率谱可得到，并可画出。我们可以此作为比较各种算法的依据。,实际工作中，对信号 总取有限长，如，由这128点去“求”功率谱，得到的当然是估计值。,调整，可以得到不同的信噪比，本例取这样，,（a）真实谱；（b）周期图；（c）Welch平均，四段，无迭合，Ham

15、ming窗；（d）同c,但迭合16点,（a）真实谱；（b）周期图；（c）Welch平均，四段，无迭,（e）BT法，M32；（f）BT法，M16,（e）BT法，M32；（f）BT法，M16,经典功率谱估计的特点：,1.物理概念明确，可用FFT快速算法。所以 是大众化的谱估计方法；,2.对周期图，分辨率受到 的限制；对自相关法，分辨率受到 的限制；,3.方差性能不好，不是一致估计，N 增 大时谱曲线反而起伏加剧；,经典功率谱估计的特点：1.物理概念明确，可用FFT快速算,4.改进方法是“平滑”与“平均”，改进的目 的是减小方差，但牺牲了分辨率；,5.注意窗函数的作用与影响：,4.改进方法是“平滑”

16、与“平均”，改进的目 5.注意窗,11.7 短时傅里叶变换,平稳信号：均值、方差及均方都不随时间变化，自相关函数仅和两个观察时间的差有关，和观察的具体位置无关；,非平稳信号：均值、方差都随时间变化，自相关函数也和观察的时间位置有关，信号的频率也随时间而变化，如语音、脑电及其他含有较多突变分量的信号。其一阶、二阶统计量和功率谱的估计显然不能简单地使用平稳信号的估计方法，必须考虑其时变因素。,方法：分段，每一小段可看作是平稳的。,11.7 短时傅里叶变换平稳信号：均值、方差及均方都不随时间,概念：,其STFT定义为:,并且窗函数应取对称函数。,式中,概念：其STFT定义为:并且窗函数应取对称函,x

17、()0FTFTFT0,谱图是恒正的，且是实的。,概念：,“谱图（spectrogram）”,由于,所以,考虑 是随机信号的一个样本，谱图可实现信号功率谱的估计。注意，它们是 的函数,谱图是恒正的，且是实的。概念：,将信号 变换为一个二维函数 的方法称为信号的联合时频分析：,STFT,谱图,Wigner分布,Cohen类分布,将信号 变换为一个二维函数,例1,信号x(n)由三个不同频率的正弦首尾相接所组成，即,例1信号x(n)由三个不同频率的正弦首尾相接所组成，即,例2,线性频率调制信号（chirp）,其频率与时间成正比非平稳,例2线性频率调制信号（chirp）其频率与时间成正比,例2,两个ch

18、irp信号，一个频率随时间增长，一个频率随时间减小，求它们和的谱图。,例2两个chirp信号，一个频率随时间增长，一个频率随时间减,与本章内容有关的MATLAB文件,pwelch.m 本文件用Welch平均法估计一个信号的功率谱，其基本调用格式是：Px,F=pwelch(x,Nfft,Fs,window,Noverlap)式中 x 是随机信号，Fs是抽样频率，Nfft是对x作FFT时的长度，window是选用的窗函数,Noverlap是估计x的功率谱时每一段叠合的长度。缺省时，Nfft=256,noverlap=0,window=Hanning(Nfft),Fs=2。输出的Px是估计出的功率谱，按上述调用格式给出的是幅平方值，F是频率轴坐标。,与本章内容有关的MATLAB文件 pwelch.m 本文件用,2.spectrum.m功能和pwelch.m类似，可用Welch平均法来估计一个信号的自功率谱，还可用于估计两个信号的互功率谱。,3specgram.m 估计信号的谱图，但实际上估计的是其短时傅里叶变换。该文件主要针对非平稳信号，当然也可用于平稳信号，甚至确定性信号。,2.spectrum.m3specgram.m,