毕业设计(论文)外文翻译-基于加窗傅里叶变化的测频计算

上海电力学院本科毕业设计（论文）专业译文原文题目：DFT-basedfrequencyestimationunderharmonicinterference译文题目：谐波干扰下基于DFT的频率估计学生姓名：学号：指导教师：院系：电气工程学院专业年级：电气工程及其自动化2016届2017年1月3日谐波干扰下基于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个或者更多的正弦声波，因此我们主要关注频率估计算法为：避免出现迭代程序，并且计算简单；避免计算数据向量长度更大的DFT；当不仅存在噪声，而存在干扰正弦声波时，将对噪声的估算精确度和稳定性最大化。为了评价其性能，对于不同的估算错误我们采用了Cramér-Rao下限作为参考，此估算错误获得时具有无偏的最大可能（ML）估计器。本文的结构如下。在第二部分，我们叙述了存在的问题，确定了CRLB和测试的条件。在第三部分中，我们对测试中使用的不同窗之间的关系进行了定义和分类。在第四和第五部分中，我们给出了研究中估算得到并在矩形或汉宁窗中使用的参考非迭代频率估计器。在第六部分中，我们给出了两个新的反正切窗函数估计器，并在单个文章的进行了叙述，这些估计器也包含在我们的对比研究中。第七部分讨论了不同频率估计器的相对性能，第八部分对本文进行了总结。2存在的问题为了简单但不失一般性，我们使用x代表任意频率的正弦信号，此信号受到了高斯噪声r(n)的破坏。通过公式=2π/N（l+）得到了正弦声波的频率，其中l和分别代表DFT元素大小中的整数部分（0<l<N/2）和小数（0.0<<1.0，或者取决于插入规则-0.5<<0.5）部分。对于真实的正弦声波：其中A代表正弦声波的振幅，Φ代表正弦声波的相位。可以通过简单的方法进行估算：对于输入的信号，根据公式（1），可以通过分析h(n)得到信号加倍之后l和的值，然后使用N-点DFT将其转化成频域：因此，存在问题的主要目标是寻找一个算法或公式，这个公式可以使用V(k)的值，并插入的值，因为l相当于|V(k)|的最小值，所以很容易确定。在复杂的正弦声波中，对于无偏的ML频率估计值，可以用公式[9]计算得到CRLB：其中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，的取值为[0.0，1.0[(或者对于一些估计器，l=20或l=20，的取值为[-0.5，0.5])]。这意味着，目标正弦声波和两个相邻声波的分离在10.3二进制到11.3二进制之间变动。因此，从现在开始，可以我们假设的新的测试场景的条件，通过设定标准化的CRLB的大小。3窗函数及其影响本文中考虑了不断降低的选择性和不断增加的主侧旁波瓣衰减的三个窗：矩形窗正弦窗和汉宁窗矩形窗、正弦窗和汉宁窗主要旁波瓣的宽度分别为4π/N、6π/N和8π/N。主要旁波瓣越模糊，选择性越好，因为较近空间的正弦声波能够更好的保留下来。对于此观点，矩形窗具有最好的选择性，汉宁窗具有最差的选择性。另外，主要旁波瓣越大，主旁波瓣和侧旁波瓣之间的衰减就越强，这个特征可以成为“泄漏”。因此，可以说矩形窗具有最大的泄漏，汉宁窗具有最小的泄漏。泄漏越小，DFT谱图中两个保留正弦声波之间的相互影响就越小。当信号中存在噪声和其他干扰正弦声波时，有一些重要的相位可能会影响估算过程的性能。为了尽量避免泄漏，当使用矩形、正弦或汉宁窗时，需要分别使用光谱峰值周围2、3或4个最大的DFT谱线来插入的值。4基于矩形窗函数的估计器使用四个非迭代、基于DFT的频率估计器来假定我们的估算中选择的矩形窗，本研究以简化、性能和以前的仿真为基础。Jain等。频率估计者意识到如果，那么两条在wl处并与本地波峰有关系的光谱线|V(l)|和|V(l+1)|是最大的。由于负轴上存在峰值（并假设l>20），因此可以通过忽略泄漏并使用下列公式得到频率的估计值：Jain等认为性能中的损失是由“谐波干涉”太小引起的。Quinn提出了一个频率估计器，这个估计器使用了本地最大峰值（k=l）每一个上的DFT波谱线，以此来提升对噪声的稳定性。令αL=，其中，以及，其中α=，频率估计可以通过下面的公式得到：其中，Quinn给出了单个真实正弦声波小于CRLB的倍时这些估计器的渐进方差。Macleod开发了一个三样本频率插入程序，这个程序涉及到DFT谱图中的一个峰值样本和两个相邻的样本。为了提升其性能，有必要在频率估计器中使用DFT相位和振幅的信息，并第一次计算了、和，从而得出：最后计算得到：据报道，相关的平均方差（在复杂的正弦声波的情况下）是CRLB的1.32倍左右。Jacobsen最近使用矩形窗提出了非常简单、高效的DFT频率估计器。在以前的两种情况下，使用谱图峰值中心的三个DFT样本，并根据下面的公式可以得到频率的估计值：研究者认为这个简单的估计值对于非常低的性噪比来说出人意料的精确，通过其有能力消除统计上偏差可以从在一定程度上解释这种现象。5以汉宁窗为基础的估计器我们选择了五个非迭代、基于DFT的频率估计器，并以报道的简单化与性能以及少量初步仿真为基础来假定汉宁窗。Grandke认为由于“谐波干涉”引起的泄漏是一个具有矩形基础频率插入值的问题，并且认为使用汉宁窗的频率估计器（不存在较长的泄漏）为：在使用相同的条件时（使用20DFT二进制将三个真实的正弦声波分离），Grandke给出的DFT频率估计器比Jain的结果要精确很多。对于谐波干涉，Grandke还预测比汉宁窗多很多的复杂窗可能会包围现有的限制条件，因为这些声调都被充分的隔开了。除了频率估计器使用了前面一节中提到的矩形窗，Macleod还提出了一种使用汉宁窗的频率估计器，这种估计器具有“充分的泄漏抵抗性”。通过计算第一、和，可以使用下面的公式估算得到频率：当考虑单个复杂正弦声波的估算时，本估计器的平均方差大约是CRLB的2.13倍。一个非常简单、流行并且经常使用的DFT频率估计就是抛物线插入。抛物线插入涉及到的抛物线要和对数刻度上分析窗频率反应振幅的主波瓣相匹配。由于主波瓣的顶部存在一个凸出的形状，X-Y平面上的模型可以由公式计算得到，其中未知参数分别为水平位移（x0）、垂直位移（y0）和凸面参数（m）。因此，需要三个方程才能计算。使用代表正弦声波的未知振幅，而中只有是未知的。使用可以得到三个方程、和，使用下面的式子计算出：抛物线插入需要所有可能的，以及适应分析窗频率反应主波瓣内部的三条DFT直线。因此，这表明如果对矩形或者正弦窗使用抛物线插入，那么抛物线插入是无效的。当对汉宁窗使用抛物线插入时，最大的绝对估计误差是二进制宽度的1.6%，这是相对较差的情况，并且是由抛物线插入的窗不可知性导致的。可靠地估计需要两个正弦声波之间的频率分离至少在4条DFT频率线（或者8π/N）上。如果在频率插入的初级阶段使用补零操作，那么二次插入的性能可能会有所提升。与Quinn在第四部分中提出的频率估计器类似，Quinn最近还提出了一种在汉宁窗中使用的新的估计器。定义，其中，以及，其中，频率的估算方程为：其中由于主要来自特定的性噪比，因此报道的真实正弦声波的性能结果没有太大的意义。除了第四部分中Jacobsen提出的频率插入程序外，研究者还在同一篇文章中提出了两个使用汉宁窗的频率插入法，包括：然而性能的结果却是单个音调以及较小的性噪比范围内（-2dB至10dB）。6基于反正切窗函数的插入法在对比文章中，我们给出了两个基于反正切窗函数的DFT频率估计器，这两个估计器适用于矩形窗和正弦窗。第一个估计器使用了和，两个最大的DFT谱线，因为两条DFT线和矩形窗频率回应的主波瓣是匹配的：对于复杂的正弦声波和不存在噪声的情况是非常精确的。然而，对于真是的正弦声波，此频率插入的精确度取决于l和，并独立于正弦声波的振幅和相位。对l的依赖性可以通过较小的l进行确定，在l为20或更小时，泄漏会非常明显，并且会引入系统系估计误差。第二个估计器使用了三个最大的DFT二进制，因为这三个DFT二进制和正弦窗频率回应的主波瓣是匹配的。考虑下面的公式：其中，，是本地最大值，，这表明精确的频率可以估算为；而当精确的频率被估算为时，。因此，的估算必须在0.0至0.5之间。在第一种情况下，当0.0<<γ时，其中γ为优化的参数，可以通过使用公式得到改进的估计精确值，进而得到：当，时，可以通过使用公式得到改进的估计精确值，进而得到：极小极大场景中优化的估计误差可以确定为r、G和F。在第二种情况下，可以在重新确定和之后获得同样的表达方式。由文献[18]可知，与二进制宽度相关的最大绝对估计误差通常小于0.1%。7仿真本节中我们评价了基于矩形窗的频率估计器的相对性能，此估计器在第四部分中有所描述，并且包含反正切窗函数估计器和第五部分中给出了基于汉宁窗的频率估计器的相对性能，还包含结合公式（20）和（21）得到的反正切窗函数估计器。A基于矩形窗的插入的性能得到的基于DFT的频率估计器的仿真分别记为Jain79、Quinn97、Macleod98、反正切窗函数和Jacobsen79，并且和公式（9）、（10）、（11）、（18）和（12）是相关的。图1给出了使用二进制宽度标准化之后每一个频率估计器引起的误差方差的平方根。可以得到的结论是可以确定三个区域，这三个区域中特定的估计器给出了获得CRLB的优势：1）对于性噪比处于-10dB至0dB之间的区域，Jain70和反正切窗函数的估计器相对于其余估计器来说具有很小的优势；2）对于性噪比处于0dB至10dB之间的区域，Macleod98的估计器相性能出了边缘相对优势；3）对于性噪比大于10dB的区域，Jacobsen07没有竞争对手。总之，尽管在中等（例如10dB附近）和较高的性噪比时，所有的频率估计器都性能在RMSE附近或者比二进制宽度的0.02%还要好，这对于研究的目的是可以接受的，由于在非常低的性噪比时性能很差，因此需要Quinn97应该避免，而在较高的性噪比时具有很好的性能，因此最好选择Jacobsen97估计器。结果表明Jacobsen97估计器不仅能够消除“统计偏差”，而且能够稳定谐波干扰。最后的评论考虑了Jain79和反正切窗函数估计器之间非常类似的性能，以及相对于剩余基于矩形窗的估计器来说较差的整体性能。一种可信的解释为这些估计器只是用了来自两个DFT二进制的振幅信息，而剩余的三个估计器使用了来自三个DFT二进制的振幅信息和相位信息。图1以性噪比为根据时五种基于矩形窗的频率估计器的RMSE（占标准化二进制宽度的%）。以CRLB为参考。图2除了目标声波之外，当干扰声波为10二进制左右时且以性噪比为根据时六种汉宁/正弦窗基频率估计器的RMSE（占标准化二进制宽度的%）。以CRLB为参考。B基于汉宁/正弦窗的插入的性能对于基于DFT的频率插入获得的仿真可以分别确定为Grandke83、Macleod98、Quadratic、Quinn06、反正切窗函数和Jacobsen79，并且和公式（13）、（14）、（15）、（16）、（20-21）和（17）是一致的。由每个频率估计器导致的与标准化误差方差有关的结果见图2，这些结果表明，对于性噪比为15dB或者更大时，Quadratic估计器的性能是第一个需要饱和的，其次是性噪比为30dB或者更大时Jacobsen07估计器的性能，然后是性噪比为40dB左右时反正切窗函数S估计器的性能。剩余三个估计器在与CRLB具有几乎相同的距离时具有类似的趋势。然而，非常有趣的是，与其余估计器相比，反正切窗函数S估计器性能出了边缘优势，这是因为当性噪比在0-30dB之间时，标准化的RMSE能够达到更好的CRLB。为了更好的理解不同正弦声波（可能是谐波正弦声波）之间的贴近程度如何影响其性能的，使用5二进制代替10二进制将目标正弦声波周围的两个正弦声波分开。最新测试环境的仿真见图3。图3反映了由于干扰声波对目标声波的贴近性而对RMSE性能的影响，具有很好的教育意义。实际上，所有频率估计器的性能饱和都处在10dB性噪比至30dB性噪比之间。然而，可以确定的是，当性噪比处于0dB至20dB之间时，反正切窗函数S估计器与剩余的频率估计器相比保持着边缘优势，这是由于其RMSE能够得到更好的CRLB导致的。这个事实可能与窗的选择性有关。图4代表经过汉宁或正弦窗处理后结果信号的短时间傅里叶谱图。可以看出，在这种情况下，不同的正弦声波之间具有更好的分离。另一方面，正弦窗比汉宁窗具有更好的选择性。图3除了目标声波之外，当干扰声波为5二进制左右时且以性噪比为根据时六种基于汉宁/正弦窗的频率估计器的RMSE（占标准化二进制宽度的%）。以CRLB为参考。图4目标声波和除了目标之外另两种频率为5.3二进制的干扰声波的短时间傅里叶谱图（l+Δl=20.5），N=512总之，图3表明，对于较宽的性噪比以及正弦声波具有非常贴近的空间时，Grandke83、反正切窗函数S、Macleod98和Quinn06频率估计器是可以接受的选择。8结论本文中对比了几种非迭代DFT频率估计器的性能，这些估计器使用了矩形或汉宁窗，其中还包括两个使用矩形和正弦窗的基于新型反正切窗函数的频率估计器。由于存在噪声或者感染声波，因此在测试环境中将其考虑为最可能的具有自然声调的场景，比如歌声。结果表明在非常低的性噪比的环境下，基于矩形窗函数的插入比基于汉宁窗函数的插入要性能得更好。这正好与中等和较高的性噪比相反，其性能出的优势超过了振幅的顺序。还可以总结出的是正弦声波的干扰会从三个方面降低所有估计器的性能：更加偏离CRLB、在较高的性噪比时通过渐近线趋势来限制其性能、当声波干扰更强时，不同估计器的性能变得更加和谐，特别是使用汉宁窗时。最后来讨论抛物线插入。当使用汉宁窗时，这是一种非常流行的频率估算方法，并且结果表明其性能在所有使用这种窗的测试估计器中时最差的，这可以用抛物线插入的窗口不可执行来简单的解释。然而，如果实施了频率插入的第一阶段，例如补零，那么抛物线的插入就会对结果有所改善。补零是一种通过将0加入到数据向量中并计算DFT中更大的N来插入信息的技术。然而，由于补零会增加计算的复杂性，因此我们在本文中没有考虑这种方法。基于C8051F单片机直流电动机反馈控制系统的设计与研究基于单片机的嵌入式Web服务器的研究MOTOROLA单片机MC68HC（8）05PV8/A内嵌EEPROM的工艺和制程方法及对良率的影响研究基于模糊控制的电阻钎焊单片机温度控制系统的研制基于MCS-51系列单片机的通用控制模块的研究基于单片机实现的供暖系统最佳启停自校正（STR）调节器单片机控制的二级倒立摆系统的研究基于增强型51系列单片机的TCP/IP协议栈的实现基于单片机的蓄电池自动监测系统基于32位嵌入式单片机系统的图像采集与处理技术的研究基于单片机的作物营养诊断专家系统的研究基于单片机的交流伺服电机运动控制系统研究与开发基于单片机的泵管内壁硬度测试仪的研制基于单片机的自动找平控制系统研究基于C8051F040单片机的嵌入式系统开发基于单片机的液压动力系统状态监测仪开发模糊Smith智能控制方法的研究及其单片机实现一种基于单片机的轴快流CO〈,2〉激光器的手持控制面板的研制基于双单片机冲床数控系统的研究基于CYGNAL单片机的在线间歇式浊度仪的研制基于单片机的喷油泵试验台控制器的研制基于单片机的软起动器的研究和设计基于单片机控制的高速快走丝电火花线切割机床短循环走丝方式研究基于单片机的机电产品控制系统开发基于PIC单片机的智能手机充电器基于单片机的实时内核设计及其应用研究基于单片机的远程抄表系统的设计与研究基于单片机的烟气二氧化硫浓度检测仪的研制基于微型光谱仪的单片机系统单片机系统软件构件开发的技术研究基于单片机的液体点滴速度自动检测仪的研制基于单片机系统的多功能温度测量仪的研制基于PIC单片机的电能采集终端的设计和应用基于单片机的光纤光栅解调仪的研制气压式线性摩擦焊机单片机控制系统的研制基于单片机的数字磁通门传感器基于单片机的旋转变压器-数字转换器的研究基于单片机的光纤Bragg光栅解调系统的研究单片机控制的便携式多功能乳腺治疗仪的研制基于C8051F020单片机的多生理信号检测仪基于单片机的电机运动控制系统设计Pico专用单片机核的可测性设计研究基于MCS-51单片机的热量计基于双单片机的智能遥测微型气象站MCS-51单片机构建机器人的实践研究基于单片机的轮轨力检测基于单片机的GPS定位仪的研究与实现基于单片机的电液伺服控制系统用于单片机系统的MMC卡文件系统研制基于单片机的时控和计数系统性能优化的研究基于单片机和CPLD的粗光栅位移测量系统研究单片机控制的后备式方波UPS提升高职学生单片机应用能力的探究基于单片机控制的自动低频减载装置研究基于单片机控制的水下焊接电源的研究基于单片机的多通道数据采集系统基于uPSD3234单片机的氚表面污染测量仪的研制基于单片机的红外测油仪的研究96系列单片机仿真器研究与设计基于单片机的单晶金刚石刀具刃磨设备的数控改造基于单片机的温度智能控制系统的设计与实现基于MSP430单片机的电梯门机控制器的研制基于单片机的气体测漏仪的研究基于三菱M16C/6N系列单片机的CAN/USB协议转换器基于单片机和DSP的变压器油色谱在线监测技术研究基于单片机的膛壁温度报警系统设计基于AVR单片机的低压无功补偿控制器的设计基于单片机船舶电力推进电机监测系统基于单片机网络的振动信号的采集系统基于单片机的大容量数据存储技术的应用研究基于单片机的叠图机研究与教学方法实践基于单片机嵌入式Web服务器技术的研究及实现基于AT89S52单片机的通用数据采集系统基于单片机的多道脉冲幅度分析仪研究机器人旋转电弧传感角焊缝跟踪单片机控制系统基于单片机的控制系统在PLC虚拟教学实验中的应用研究基于单片机系统的网络通信研究与应用基于PIC16F877单片机的莫尔斯码自动译码系统设计与研究基于单片机的模糊控制器在工业电阻炉上的应用研究基于双单片机冲床数控系统的研究与开发基于Cygnal单片机的μC/OS-Ⅱ的研究基于单片机的一体化智能差示扫描量热仪系统研究基于TCP/IP协议的单片机与Internet互联的研究与实现变频调速液压电梯单片机控制器的研究基于单片机γ-免疫计数器自动换样功能的研究与实现基于单片机的倒立摆控制系统设计与实现单片机嵌入式以太网防盗报警系统基于51单片机的嵌入式Internet系统的设计与实现单片机监测系统在挤压机上的应用MSP430单片机在智能水表系统上的研究与应用基于单片机的嵌入式系统中TCP/IP协议栈的实现与应用单片机在高楼恒压供水系统中的应用基于ATmega16单片机的流量控制器的开发基于MSP430单片机的远程抄表系统及智能网络水表的设计基于MSP430单片机具有数据存储与回放功能的嵌入式电子血压计的设计HYPERLINK"/detai