基于加窗傅里叶变化的测频计算

1上海电力学院本科毕业设计（论文）专业译文原文题目：DFT-basedfrequencyestimationunderharmonicinterference译文题目：谐波干扰下基于DFT的频率估计学生姓名：李诗雨学号：20133643指导教师：许仪勋院系：电气工程学院专业年级：电气工程及其自动化2016届2017年1月3日2谐波干扰下基于DFT的频率估计AníbalFerreira和RicardoSousa摘要：在本文中，我们介绍了自然信号正弦声波频率的精确估计，例如歌声、嗓音或乐曲。这些信号本质上都是谐波，并且通常混有一定的噪声。以Cramér-Rao的无偏频率估计器的下限为参考，我们对比了一些基于DFT的频率估计器的性能，这些估计器是非迭代的，并且使用了矩形窗或汉宁窗。测试环境模拟谐波干涉，在测试中还是用了两个新的基于反正切窗函数的频率估计器。结果表明，不同频率估计器的相对性能可以表示为性噪比的函数。1引言很多信号处理问题都需要估算正弦声波的频率、振幅和相位参数，特别是涉及讲话或者音频编码、用于MIDI转录的PCM和实时精确歌曲分析的情况。振幅和相位的估算通常取决于频率的估算，因此，本文仅关注频率参数。讲话或者唱歌的波谱通常由一些正弦声波组成，这些正弦声波大体上遵循谐波的形式。当存在干扰信号时，包括噪声或其他正弦声波，单个正弦声波的精确频率估算通常使用离散信号的离散傅里叶变换（DFT）和从DFT谱图中提取信号来实施。由于DFT谱图大部分在频率上都是离散的，并且受到DFT自然频率分辨率的限制（2π/N，其中N是DFT的大小），因此正弦声波频率的精确估算涉及到使用一些DFT谱图（或者是DFT谱图的片段）样本进行插入计算。在过去的40年中，研究者们提出了很多DFT插入算法。测试环境的较大差异使得相对性能的评估变得非常困难。例如，在DFT之前，可能使用不同的窗函数对信号进行了计算，这对DFT的频率选择性有很大的影响，并且会产生泄漏效应。测试正弦声波可能是复杂的正弦声波（例如顺向波）或者真实的正弦声波。DFT插入的程序可能是迭代的，也可能不是迭代的。例如对于歌声信号，一个更加真实的测试场景必然涉及到一些干扰正弦声波，因为这些信号具有谐波的特征。既然我们对歌声的谐波结构中单个正弦声波的实时和精确分析很感兴趣，这种歌声中可能包含100个或者更多的正弦声波，因此我们主要关注频率估计算法为：3避免出现迭代程序，并且计算简单；避免计算数据向量长度更大的DFT；当不仅存在噪声，而存在干扰正弦声波时，将对噪声的估算精确度和稳定性最大化。为了评价其性能，对于不同的估算错误我们采用了Cramér-Rao下限作为参考，此估算错误获得时具有无偏的最大可能（ML）估计器。本文的结构如下。在第二部分，我们叙述了存在的问题，确定了CRLB和测试的条件。在第三部分中，我们对测试中使用的不同窗之间的关系进行了定义和分类。在第四和第五部分中，我们给出了研究中估算得到并在矩形或汉宁窗中使用的参考非迭代频率估计器。在第六部分中，我们给出了两个新的反正切窗函数估计器，并在单个文章的进行了叙述，这些估计器也包含在我们的对比研究中。第七部分讨论了不同频率估计器的相对性能，第八部分对本文进行了总结。2存在的问题为了简单但不失一般性，我们使用x代表任意频率lw的正弦信号，此信号受到了高斯噪声r(n)的破坏。通过公式lw=2π/N（l+l）得到了正弦声波的频率，其中l和l分别代表DFT元素大小中的整数部分（0lN/2）和小数（0.0l1.0，或者取决于插入规则-0.5l0.5）部分。对于真实的正弦声波：其中A代表正弦声波的振幅，Φ代表正弦声波的相位。可以通过简单的方法进行估算：对于输入的信号，根据公式（1），可以通过分析h(n)得到信号加倍之后l和l的值，然后使用N-点DFT将其转化成频域：因此，存在问题的主要目标是寻找一个算法或公式，这个公式可以使用V(k)的值，并插入l的值，因为l相当于|V(k)|的最小值，所以很容易确定。在复杂的正弦声波中，对于无偏的ML频率估计值，可以用公式[9]计算得到CRLB：4其中A代表正弦声波的振幅，N代表DFT的长度，σ2代表不同噪声，在这里我们假设其为零均值、白色、复杂和高斯。为了简化不同频率插入值性能的相对评估，我们DFT（2π/N）的自然频率分解将方差{}标准化，将其平方根和结果作为正常元素宽度的百分比：在模拟中，我们考虑到A=1、N=512、Φ=π/6、l=20。Δl的取值为0.0到1.0，并且每次增加0.01，根据期望值，我们得到了辐照声音矢量r(n)的100种情况，RMSE根据每一个频率估计器确定，性噪比的值在-10dB到50dB之间，并且每次增加2.5dB。然而，在歌声中一个更加现实的测试场景必然在正弦声波的每一侧（真实、目标）都包含干扰声波，需要对此声波的频率进行估算。因此，我们将公式（1）进行拓展，使其成为更加合适的形式：在模拟中，设定A=1、Φ=π/6、l=20，l的取值为[0.0，1.0[(或者对于一些估计器，l=20或l=20，l的取值为[-0.5，0.5])]。这意味着，目标正弦声波和两个相邻声波的分离在10.3二进制到11.3二进制之间变动。因此，从现在开始，可以我们假设的新的测试场景的条件，通过设定标准化的CRLB的大小。3窗函数及其影响本文中考虑了不断降低的选择性和不断增加的主侧旁波瓣衰减的三个窗：矩形窗5正弦窗和汉宁窗矩形窗、正弦窗和汉宁窗主要旁波瓣的宽度分别为4π/N、6π/N和8π/N。主要旁波瓣越模糊，选择性越好，因为较近空间的正弦声波能够更好的保留下来。对于此观点，矩形窗具有最好的选择性，汉宁窗具有最差的选择性。另外，主要旁波瓣越大，主旁波瓣和侧旁波瓣之间的衰减就越强，这个特征可以成为“泄漏”。因此，可以说矩形窗具有最大的泄漏，汉宁窗具有最小的泄漏。泄漏越小，DFT谱图中两个保留正弦声波之间的相互影响就越小。当信号中存在噪声和其他干扰正弦声波时，有一些重要的相位可能会影响估算过程的性能。为了尽量避免泄漏，当使用矩形、正弦或汉宁窗时，需要分别使用光谱峰值周围2、3或4个最大的DFT谱线来插入l的值。4基于矩形窗函数的估计器使用四个非迭代、基于DFT的频率估计器来假定我们的估算中选择的矩形窗，本研究以简化、性能和以前的仿真为基础。Jain等。频率估计者意识到如果，那么两条在wl处并与本地波峰有关系的光谱线|V(l)|和|V(l+1)|是最大的。由于负轴上存在峰值（并假设l20），因此可以通过忽略泄漏并使用下列公式得到频率的估计值：Jain等认为性能中的损失是由“谐波干涉”太小引起的。Quinn提出了一个频率估计器，这个估计器使用了本地最大峰值（k=l）每一个上的DFT波谱线，以此来提升对噪声的稳定性。令αL=，其中6，以及，其中α=，频率估计可以通过下面的公式得到：其中，Quinn给出了单个真实正弦声波小于CRLB的倍时这些估计器的渐进方差。Macleod开发了一个三样本频率插入程序，这个程序涉及到DFT谱图中的一个峰值样本和两个相邻的样本。为了提升其性能，有必要在频率估计器中使用DFT相位和振幅的信息，并第一次计算了、和，从而得出：最后计算得到：据报道，相关的平均方差（在复杂的正弦声波的情况下）是CRLB的1.32倍左右。Jacobsen最近使用矩形窗提出了非常简单、高效的DFT频率估计器。在以前的两种情况下，使用谱图峰值中心的三个DFT样本，并根据下面的公式可以得到频率的估计值：研究者认为这个简单的估计值对于非常低的性噪比来说出人意料的精确，通过其有能力消除统计上偏差可以从在一定程度上解释这种现象。75以汉宁窗为基础的估计器我们选择了五个非迭代、基于DFT的频率估计器，并以报道的简单化与性能以及少量初步仿真为基础来假定汉宁窗。Grandke认为由于“谐波干涉”引起的泄漏是一个具有矩形基础频率插入值的问题，并且认为使用汉宁窗的频率估计器（不存在较长的泄漏）为：在使用相同的条件时（使用20DFT二进制将三个真实的正弦声波分离），Grandke给出的DFT频率估计器比Jain的结果要精确很多。对于谐波干涉，Grandke还预测比汉宁窗多很多的复杂窗可能会包围现有的限制条件，因为这些声调都被充分的隔开了。除了频率估计器使用了前面一节中提到的矩形窗，Macleod还提出了一种使用汉宁窗的频率估计器，这种估计器具有“充分的泄漏抵抗性”。通过计算第一、和，可以使用下面的公式估算得到频率：当考虑单个复杂正弦声波的估算时，本估计器的平均方差大约是CRLB的2.13倍。一个非常简单、流行并且经常使用的DFT频率估计就是抛物线插入。抛物线插入涉及到的抛物线要和对数刻度上分析窗频率反应振幅的主波瓣相匹配。由于主波瓣的顶部存在一个凸出的形状，X-Y平面上的模型可以由公式计算得到，其中未知参数分别为水平位移（x0）、垂直位移（y0）和凸面参数（m）。因此，需要三个方程才能计算。使用代表正弦声波的未知振幅，而中只有是未知的。使用可以得到三个方程、和8，使用下面的式子计算出l：抛物线插入需要所有可能的l，以及适应分析窗频率反应主波瓣内部的三条DFT直线。因此，这表明如果对矩形或者正弦窗使用抛物线插入，那么抛物线插入是无效的。当对汉宁窗使用抛物线插入时，最大的绝对估计误差是二进制宽度的1.6%，这是相对较差的情况，并且是由抛物线插入的窗不可知性导致的。可靠地估计需要两个正弦声波之间的频率分离至少在4条DFT频率线（或者8π/N）上。如果在频率插入的初级阶段使用补零操作，那么二次插入的性能可能会有所提升。与Quinn在第四部分中提出的频率估计器类似，Quinn最近还提出了一种在汉宁窗中使用的新的估计器。定义，其中，以及，其中，频率的估算方程为：其中由于主要来自特定的性噪比，因此报道的真实正弦声波的性能结果没有太大的意义。除了第四部分中Jacobsen提出的频率插入程序外，研究者还在同一篇文章中提出了两个使用汉宁窗的频率插入法，包括：然而性能的结果却是单个音调以及较小的性噪比范围内（-2dB至10dB）。96基于反正切窗函数的插入法在对比文章中，我们给出了两个基于反正切窗函数的DFT频率估计器，这两个估计器适用于矩形窗和正弦窗。第一个估计器使用了和，两个最大的DFT谱线，因为两条DFT线和矩形窗频率回应的主波瓣是匹配的：对于复杂的正弦声波和不存在噪声的情况是非常精确的。然而，对于真是的正弦声波，此频率插入的精确度取决于l和l，并独立于正弦声波的振幅和相位。对l的依赖性可以通过较小的l进行确定，在l为20或更小时，泄漏会非常明显，并且会引入系统系估计误差。第二个估计器使用了三个最大的DFT二进制，因为这三个DFT二进制和正弦窗频率回应的主波瓣是匹配的。考虑下面的公式：其中，，是本地最大值，，这表明精确的频率可以估算为；而当精确的频率被估算为时，。因此，l的估算必须在0.0至0.5之间。在第一种情况下，当0.0lγ时，其中γ为优化的参数，可以通过使用公式得到改进的估计精确值，进而得到：当，时，可以通过使用公式10得到改进的估计精确值，进而得到：极小极大场景中优化的估计误差可以确定为r、G和F。在第二种情况下，可以在重新确定和之后获得同样的表达方式。由文献[18]可知，与二进制宽度相关的最大绝对估计误差通常小于0.1%。7仿真本节中我们评价了基于矩形窗的频率估计器的相对性能，此估计器在第四部分中有所描述，并且包含反正切窗函数估计器和第五部分中给出了基于汉宁窗的频率估计器的相对性能，还包含结合公式（20）和（21）得到的反正切窗函数估计器。A基于矩形窗的插入的性能得到的基于DFT的频率估计器的仿真分别记为Jain79、Quinn97、Macleod98、反正切窗函数和Jacobsen79，并且和公式（9）、（10）、（11）、（18）和（12）是相关的。图1给出了使用二进制宽度标准化之后每一个频率估计器引起的误差方差的平方根。可以得到的结论是可以确定三个区域，这三个区域中特定的估计器给