[发明专利]一种基于CZT的频率估计方法

说明书

本发明涉及一种基于CZT的频率估计方法，包括步骤：S1：对待估计信号进行采样，并对采样后的信号进行FFT变换，获得最大谱线频率索引；S2：根据FFT的最大谱线频率索引确定频率细化区间，对采样后的信号进行CZT变换；S3：对CZT变换后的频谱进行搜索，获得最大谱线的幅值及其左右两根谱线的幅值；S4：根据最大谱线右侧谱线的幅值与最大谱线的幅值的比值以及最大谱线左侧谱线的幅值与最大谱线的幅值的比值确定待估计信号的瞬时频率与CZT变换后的频率的误差；S5：根据所述误差确定待估计信号的瞬时频率。本发明的基于CZT的频率估计方法，能够以较低的细化倍数和较小的计算量，获得较高精度的频率估计结果。

技术领域

本发明涉及频率估计领域，更具体地涉及一种CZT的频率估计方法。

背景技术

信号的频率估计是工程应用中常常存在的问题，许多场景需要对信号的频率进行精确估计。例如，在频率调制连续波(FMCW)雷达测距系统中，其可以根据发送波和反射波的差频信号来得到相关的距离信息，即通过估计差频信号的频率来得到估计距离，因此对于差频信号的频率估计准确程度直接影响到了估计距离的准确程度，差频信号频率的测量精度直接决定了测距的精度。

为了提升频率估计的精度，许多算法在FFT变换(傅里叶变换)的基础上进行了改进，例如选带傅里叶变换(Zoom-FFT)、Rife算法、线性调频Z变换(CZT)等。Zoom-FFT算法(赵强,侯孝民,廉昕.基于Zoom-FFT的改进Rife正弦波频率估计算法[J].数据采集与处理,2017,32(04):731-736.DOI:10.16337/j.1004-9037.2017.04.009)对于信号进行移频，再进行降采样，通过降采样提升了频谱分辨率，从而提升了频率估计的精度；Rife算法利用频谱幅值最大的谱线及其相邻谱线中幅值较大谱线来计算实际频率与最大谱线对应的频率的偏差值，从而提升频率估计的精度。但是，Rife算法存在估计精度受噪声影响大、偏差值较小时估计误差较大的问题。孙宏军等(孙宏军,王小威.基于幅值-相角判据的修正Rife正弦波频率估计算法[J].天津大学学报(自然科学与工程技术版),2018,51(08):810-816)通过设定频移门限值，对频移因子满足门限值的信号使用Rife算法，否则便使用相角来估计频移因子，提升了算法的性能。Quinn(Quinn B G.Estimation of frequency，amplitudeand phase from the DFT of a time series[J].IEEE Trans on Signal Processing，1997，45(3)：814-817)提出利用FFT变换主瓣内次大值与最大值之比的实部代替幅值之比的频率估计方法，减少了由于频移因子方向判断错误引起的误差。CZT算法(欧阳鑫信,贺青.基于CZT的多子带信号时差估计方法[J].现代雷达,2019,41(02):31-34.DOI:10.16592/j.cnki.1004-7859.2019.02.007)对于部分频段进行细化，相较于FFT算法具有更高的频率估计精度，但是需要大量的计算量。刘上(刘上,朱国富,王玲,等.LFMCW雷达高精度测距相位差改进算法[J].雷达科学与技术,2021,19(03):287-291+297)提出了基于CZT的相位差算法，对于信号进行CZT变换，求出信号相位，利用相位差来求解信号的频率，解决了相位差算法受噪声影响大缺陷，提升了频率估计的精度。

综上所述，现有的频率估计算法均存在受噪声影响大、计算量大等缺陷，无法满足目前工业领域对于频率估计的精度高、运算快的要求。

发明内容

本发明的目的在于提供一种基于CZT的频率估计方法，其计算量低、抗噪性能好。

本发明提供一种基于CZT的频率估计方法，包括以下步骤：

S1：对待估计信号进行采样，并对采样后的信号进行FFT变换，获得最大谱线频率索引；

S2：根据FFT的最大谱线频率索引确定频率细化区间，对采样后的信号进行CZT变换；