基于补零的主瓣拟合三谱线插值频谱校正算法研究
胡婧雯1 贾豫东1,2（通讯作者）

1 北京信息科技大学 仪器与光电学院，北京，100192；

2 生物医学成像北京实验室，北京，100871；

摘要：针对离散信号在低信噪比条件下频谱峰值易受噪声淹没、以及离散傅里叶变换栅栏效应导致频率估计精度不足的问题，本文提出一种基于补零的窗函数主瓣拟合三谱线插值算法，该算法利用汉宁窗主瓣的精确几何特征建立频偏与谱线幅值比的数值映射关系并推导至各种通用窗函数，通过构建频偏与不同窗函数谱线幅值比的映射关系，在消除算法系统误差的同时，避免复杂实时方程的求解。仿真结果表明，在低信噪比下，该算法基于各个窗函数的主瓣拟合误差在[-0.5，0.5]内均在10-4左右，系统误差均小于10-4，此外，采用本文算法在不同信噪比下的均方根误差RMSE约为Hz，相较于其他经典插值算法具有更低的RMSE和及更强的抗噪声性能。

关键词：离散信号；主瓣拟合；加窗补零；插值法

参考文献

[1]靳笑晗, 汪岳峰, 竹孝鹏, 等. 相干多普勒测风激光雷达低信噪比区域回波信号的估计方法[J]. 光学与光电技术, 2013, 11(3): 10-14. DOI: CNKI:SUN:GXGD.0.2013-03-004.

[2]孙宸, 胡嘉龙, 王昆, 等. 基于小波变换的激光雷达回波信号去噪方法研究[J]. 雷达与对抗, 2022, 42(4): 20-22, 57. DOI: 10.19341/j.cnki.1009-0401.2022.04.006.

[3]Agrež D. Weighted multipoint interpolated DFT to improve amplitude estimation of multifrequency signal[J]. IEEE Transactions on Instrumentation and Measurement, 2002, 51(2): 287-292.

[4]Belega D, Dallet D. Frequency estimation via weighted multipoint interpolated DFT[J]. IET Science, Measurement and Technology, 2008, 2(1): 1-8.

[5]Quinn B G. Estimating frequency by interpolation using Fourier coefficients[J]. IEEE Transactions on Signal Processing, 1994, 42(5): 1264-1268.

[6]Aboutanios E, Mulgrew B. Iterative frequency estimation by interpolation on Fourier coefficients[J]. IEEE Transactions on Signal Processing, 2005, 53(4): 1237-1242.

[7]Rife D, Boorstyn R R. Single tone parameter estimation from discrete-time observations[J]. IEEE Transactions on Information Theory, 1974, 20(5): 591-598.

[8]罗涛. 相干测风激光雷达多普勒频移估计技术研究[D]. 成都: 电子科技大学, 2022. DOI: 10.27005/d.cnki.gdzku.2022.002703.

[9]步志超. 相干测风激光雷达系统设计及数据处理算法研究[D]. 北京: 北京理工大学, 2014.

[10]Nickalls R W D. A new approach to solving the cubic: Cardan’s solution revealed[J]. The Mathematical Gazette, 1993, 77(480): 354-359.

[11]罗久飞. 离散频谱校正新方法及其抗干扰性能研究[D]. 重庆: 重庆大学, 2015.