一种DFT的频率校正算法.pdf

云南大学学报（自然科学版），２０１７，３９（２）：２００～２０６ ＤＯＩ：１０．７５４０／ｊ．ｙｎｕ Ｊｏｕｒｎａｌ ｏｆ Ｙｕｎｎａｎ Ｕｎｉｖｅｒｓｉｔｙ 一种ＤＦＴ 的频率校正算法∗ １，２ ２ 谢长贵 ，陈　 平 （１．重庆工程职业技术学院机械工程学院，重庆　 ４０２２６０，２．重庆大学 机械工程学院，重庆　 ４０００４４） 摘要：针对传统方法在工程实践中的频率估计精度不够的问题，提出了一种迭代频率估计方法．该方法具 有较强的傅里叶系数的加权样本插值功能，可适用于几乎所有窗函数；对于各种窗函数的系统误差及在高斯白 噪声下该算法表现出了比传统算法更优越的性能．实验结果表明该方法能显著减少错误引发的错误的谱线的 位置，高斯噪声条件下进行频率估计时该方法取得了比传统方法更小的估计误差，表现出较强的估计精度与抗 噪性能． 关键词：离散傅里叶变换；归一化频率；参数估计；补零 中图分类号：ＴＮ９１１　 　 文献标志码：Ａ　 　 － － － 文章编号：０２５８ ７９７１（２０１７）０２ ０２００ ０７ 频谱分析的关键在于根据信号的频域描述估计与分析信号的成分及特征量．实际信号中往往伴随多 个频率成分，为此需要对其进行频谱分析和谱估计．如何使谱估计值最大限度接近于真实值，已成为众多 领域专家急需解决的难题．工业生产中关于频谱分析的最典型的例子有电力系统谐波分析中电压与电流 波形失真问题，大型机械设备的噪声监测问题等，解决这些问题均需精确的频谱校正以得到各频率成分及 － ［１ ３］ 特征量 ．目前的谱分析方法中，大多数研究者主要通过对采样得到的离散信号运用ＦＦＴ 和ＤＦＴ算法以 获取各分量的幅值、频率和相位．由于存在功率泄露及栏栅效应，ＦＦＴ 和ＤＦＴ算法估计出的结果均存在较 大误差．对非相干采样信号，上述方法的频域校正误差尤其显著．当加窗时单谐波频率的最大误差可达 － ［４ ６］ ３６％ ．因此须采用各种分析方法对频谱进行校正以提高精度．离散傅里叶插值算法是众多误差校正方 法之一，该方法的思想是获得一个有一定数量的已知谱线加权平均值的信号的频率校正．该方法具有运算 速度较快，缺点是其频率校正精度还难以满足工程需要． 基于上述原因，本文提出了一种新的插值算法，该方法可通过补零技术以获得更多谱线．然后对新算 法进行扩展，并引入主瓣拟合技术以使它与其他典型窗函数相兼容．文章还研究了新算法对其它窗函数的 系统误差以及其在高斯白噪声下的性能，并与传统算法作了比较研究．通过计算机仿真表明，在加性噪声 条件下本文提出的新算法比传统算法表现更加稳定． １　 频率校正背景 设伴有高斯白噪声的单一频率的正弦信号的表达式为： ＝ ＋ ＋ 　 　 　 　 ｙ（ｔ） Ａｃｏｓ（２πｆ ｔ φ ） ｚ（ｔ）， （１） ０ ０ 式（１） 中ｔ为连续时间变量，０ ≤ｔ ≤Ｔ，Ａ、ｆ 和φ 分别为信号的幅度、频率及初始相位３个参数，ｚ（ｔ）为高 ０ ０ 斯白噪声．在工程应用，尤其在计算机处理中，仅能对离散的信号样点进行处理，在观察间隔ＮΔｔ 内对信号 以频率ｆ 进行采样，可得到Ｎ个离散正弦信号样本的表达式为：ｓ － － ∗ 收稿日期：２０１６ ０３ １５ 　 基金项目：国家自然科学基金委员会与中国工程物理研究院联合基金． － － 　 作者简介：谢长贵（１９８４ ），男