CN113219243A - 一种校正频率搜索间隔量化误差的高精度测频方法 - Google Patents

一种校正频率搜索间隔量化误差的高精度测频方法 Download PDF

Info

Publication number
CN113219243A
CN113219243A CN202110458637.5A CN202110458637A CN113219243A CN 113219243 A CN113219243 A CN 113219243A CN 202110458637 A CN202110458637 A CN 202110458637A CN 113219243 A CN113219243 A CN 113219243A
Authority
CN
China
Prior art keywords
frequency
search interval
quantization error
amplitude spectrum
peak position
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202110458637.5A
Other languages
English (en)
Other versions
CN113219243B (zh
Inventor
陈镜
倪丽花
万群
王宗全
朱玉军
郭贤生
庄杰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN202110458637.5A priority Critical patent/CN113219243B/zh
Publication of CN113219243A publication Critical patent/CN113219243A/zh
Application granted granted Critical
Publication of CN113219243B publication Critical patent/CN113219243B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R23/00Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
    • G01R23/02Arrangements for measuring frequency, e.g. pulse repetition rate; Arrangements for measuring period of current or voltage
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R23/00Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
    • G01R23/16Spectrum analysis; Fourier analysis

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Measuring Frequencies, Analyzing Spectra (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明属于电子信息技术领域,具体为一种校正频率搜索间隔量化误差的高精度测频方法。本发明首先利用信号采样的频谱序列和信号采样的尾首差确定幅度谱及其谱峰位置,然后由频谱序列在幅度谱谱峰位置的值和信号采样的尾首差确定频率搜索间隔的量化误差校正值,进而由幅度谱谱峰位置和频率搜索间隔的量化误差校正值确定频率测量。本发明的有益效果是:使用本发明提出的一种校正频率搜索间隔量化误差的高精度测频方法,不仅利用了幅度谱谱峰位置,还利用了频率搜索间隔的量化误差校正值,实现高精度测频。本发明提出的方法的计算量与快速傅里叶变换法的计算量是同一数量级的。

Description

一种校正频率搜索间隔量化误差的高精度测频方法
技术领域
本发明属于电子信息技术领域,具体为一种校正频率搜索间隔量化误差的高精度测频方法。
背景技术
对信号频率进行快速、准确的测量在雷达、通信、声纳等领域具有重要意义。测量频率的方法有很多,按照工作原理可以分为计数法、相位法、快速傅里叶变换法等。测量频率的计数法、相位法适用于单频信号,当信号中包含多个频率成分时,不再适用。快速傅里叶变换法先对信号进行傅里叶变换,再搜索幅度谱的峰值,最后利用幅度谱的峰值位置确定频率测量,适用于多频信号,但是受到幅度谱频率搜索间隔的量化误差的影响,频率测量的误差较大。
由于快速傅里叶变换的幅度谱频率搜索间隔与信号采样数或快速傅里叶变换的长度成反比,为了降低快速傅里叶变换法幅度谱频率搜索间隔的量化误差对频率测量精度的影响,常常采用填零的方法增加信号采样数或快速傅里叶变换的长度。但是,快速傅里叶变换的计算量与信号采样数或快速傅里叶变换长度成正比,这种提高测频精度的方法也会使频率测量的计算量成倍的增加,难以同时满足高精度、有实时性要求的应用场合。
发明内容
针对上述问题,本发明提出一种同时利用信号采样的快速傅里叶变换的幅度谱峰值位置和频率搜索间隔量化误差校正值确定信号频率的方法,能实现高精度快速测频的目的。
本发明的技术方案是:
一种校正频率搜索间隔量化误差的高精度测频方法,首先利用信号采样的频谱序列和信号采样的尾首差确定幅度谱及其谱峰位置,然后由频谱序列在幅度谱谱峰位置的值和信号采样的尾首差确定频率搜索间隔的量化误差校正值,进而由幅度谱谱峰位置和频率搜索间隔的量化误差校正值确定频率测量。记信号采样数为N,信号采样为x(n),n=0,1,2,…,N-1,包括以下步骤:
S1、由前N-1个信号采样x(n),n=0,1,2,…,N-2,采用快速傅里叶变换,确定频谱序列y(k),为:
Figure BDA0003041417220000021
其中,k=0,1,2,…,N-2;
S2、由信号采样的首尾样本,确定信号采样的尾首差v,为:
v=x(N-1)-x(0)
S3、由信号采样的尾首差v和频谱序列y(k),确定幅度谱z(k),为:
z(k)=|y(k)+v|+|y(k)|
其中,k=0,1,2,…,N-2;
S4、令k=0,1,2,…,N-2,搜索幅度谱的最大值对应的序列号,即幅度谱谱峰位置m,为:
Figure BDA0003041417220000022
S5、由频谱序列的第m个值y(m),即频谱序列在幅度谱谱峰位置的值,与信号采样的尾首差v,确定频率搜索间隔的量化误差校正值u,为:
Figure BDA0003041417220000023
其中,angle表示复数的相位角。
S6、由幅度谱谱峰位置m和频率搜索间隔的量化误差校正值u,确定频率测量f,为:
Figure BDA0003041417220000024
本发明的有益效果是:使用本发明提出的一种校正频率搜索间隔量化误差的高精度测频方法,不仅利用了幅度谱谱峰位置,还利用了频率搜索间隔的量化误差校正值,实现高精度测频。本发明提出的方法的计算量与快速傅里叶变换法的计算量是同一数量级的。
具体实施方式
下面结合实施例对本发明的实用性进行说明。
实施例
在本例中,信号采样数分四种情况,N=32、64、128、256,信号采样为x(n),n=0,1,2,…,N-1,信号采样中包括归一化(相对于采样频率)频率等于0.2143、信噪比为18dB的信号,以及归一化频率等于0.1162、干信比为11dB的干扰。
测量信号频率时,快速傅里叶变换法只利用幅度谱的谱峰位置测量信号频率,本发明方法不仅利用了幅度谱的谱峰位置,还利用了频率搜索间隔的量化误差校正值。在本例中,利用同样的采样信号,当N=32、64、128、256时,快速傅里叶变换法的信号归一化测频均方根误差分别为0.0091、0.0045、0.0023和0.0011;本发明方法的信号归一化测频均方根误差分别为0.0014、0.0007、0.0003和0.0002。可见,与快速傅里叶变换法相比,利用同样的采样信号,当N=32、64、128、256时,本发明方法的归一化测频均方根误差分别降低了84.12%、85.05%、84.77%和84.67%,本发明方法在信号采样数N=32时的测频精度接近快速傅里叶变换法在信号采样数N=256时的测频精度,实现了高精度测频的目的。

Claims (1)

1.一种校正频率搜索间隔量化误差的高精度测频方法,其特征在于,包括以下步骤:
S1、采用快速傅里叶变换,根据前N-1个采样信号x(n)确定频谱序列y(k):
Figure FDA0003041417210000011
其中,n=0,1,2,…,N-2,k=0,1,2,…,N-2;
S2、由采样信号的首尾样本,确定信号采样的尾首差v:
v=x(N-1)-x(0)
S3、由采样信号的尾首差v和频谱序列y(k),确定幅度谱z(k):
z(k)=|y(k)+v|+|y(k)|
S4、搜索幅度谱的最大值对应的序列号,即幅度谱谱峰位置m:
Figure FDA0003041417210000012
S5、根据步骤S4得到的频谱序列在幅度谱谱峰位置的值y(m),与采样信号的尾首差v,确定频率搜索间隔的量化误差校正值u:
Figure FDA0003041417210000013
其中,angle表示复数的相位角;
S6、由幅度谱谱峰位置m和频率搜索间隔的量化误差校正值u,获得频率f:
Figure FDA0003041417210000014
CN202110458637.5A 2021-04-27 2021-04-27 一种校正频率搜索间隔量化误差的高精度测频方法 Active CN113219243B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110458637.5A CN113219243B (zh) 2021-04-27 2021-04-27 一种校正频率搜索间隔量化误差的高精度测频方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110458637.5A CN113219243B (zh) 2021-04-27 2021-04-27 一种校正频率搜索间隔量化误差的高精度测频方法

Publications (2)

Publication Number Publication Date
CN113219243A true CN113219243A (zh) 2021-08-06
CN113219243B CN113219243B (zh) 2021-12-31

Family

ID=77089734

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110458637.5A Active CN113219243B (zh) 2021-04-27 2021-04-27 一种校正频率搜索间隔量化误差的高精度测频方法

Country Status (1)

Country Link
CN (1) CN113219243B (zh)

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1847837A1 (en) * 2006-04-21 2007-10-24 Fondazione Torino Wireless Method and apparatus for frequency estimation
US20130015890A1 (en) * 2011-07-15 2013-01-17 Askey Computer Corp. Method and system for calibrating frequency
CN105158740A (zh) * 2015-08-24 2015-12-16 西安电子科技大学 基于高精度频率估计的噪声调幅干扰抑制方法
CN105353215A (zh) * 2015-12-23 2016-02-24 合肥工业大学 基于Nuttall窗四谱线插值FFT的谐波检测方法
CN105675083A (zh) * 2014-11-20 2016-06-15 中国航空工业集团公司雷华电子技术研究所 一种带有频率细化插值的高精度液位测量方法
CN107064629A (zh) * 2017-06-07 2017-08-18 东南大学 一种基于频率相对偏差预估的分段综合单频信号频率估计方法
CN107305223A (zh) * 2016-04-19 2017-10-31 天津大学 一种改进的相位差频率估计方法
CN107632961A (zh) * 2017-07-12 2018-01-26 天津大学 基于全相位谱分析的多频内插迭代频率估计方法及估计器
CN108108692A (zh) * 2017-12-19 2018-06-01 湖北工业大学 一种正弦信号频率识别方法
CN108414833A (zh) * 2018-03-07 2018-08-17 南京大学 一种信号分量频率的精确估计方法
CN109490862A (zh) * 2018-11-23 2019-03-19 山东航天电子技术研究所 一种基于相位差分统计谱的载频估计方法
CN110333389A (zh) * 2019-05-07 2019-10-15 东南大学 基于插值dft的正弦信号频率估计方法

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1847837A1 (en) * 2006-04-21 2007-10-24 Fondazione Torino Wireless Method and apparatus for frequency estimation
US20130015890A1 (en) * 2011-07-15 2013-01-17 Askey Computer Corp. Method and system for calibrating frequency
CN105675083A (zh) * 2014-11-20 2016-06-15 中国航空工业集团公司雷华电子技术研究所 一种带有频率细化插值的高精度液位测量方法
CN105158740A (zh) * 2015-08-24 2015-12-16 西安电子科技大学 基于高精度频率估计的噪声调幅干扰抑制方法
CN105353215A (zh) * 2015-12-23 2016-02-24 合肥工业大学 基于Nuttall窗四谱线插值FFT的谐波检测方法
CN107305223A (zh) * 2016-04-19 2017-10-31 天津大学 一种改进的相位差频率估计方法
CN107064629A (zh) * 2017-06-07 2017-08-18 东南大学 一种基于频率相对偏差预估的分段综合单频信号频率估计方法
CN107632961A (zh) * 2017-07-12 2018-01-26 天津大学 基于全相位谱分析的多频内插迭代频率估计方法及估计器
CN108108692A (zh) * 2017-12-19 2018-06-01 湖北工业大学 一种正弦信号频率识别方法
CN108414833A (zh) * 2018-03-07 2018-08-17 南京大学 一种信号分量频率的精确估计方法
CN109490862A (zh) * 2018-11-23 2019-03-19 山东航天电子技术研究所 一种基于相位差分统计谱的载频估计方法
CN110333389A (zh) * 2019-05-07 2019-10-15 东南大学 基于插值dft的正弦信号频率估计方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
FENG-XIANG GE等: ""Frequency Estimation of Sinusoidal Signals in Multiplicative and Additive Noise"", 《IEEE JOURNAL OF OCEANIC ENGINEERING》 *
XUN HUANG等: ""Frequency estimation of cyclic spectrum carrier based on compressive sampling of BPSK signal"", 《IEEE》 *
罗蓬 等: ""高精度LFM信号参数估计的谱校正方法"", 《系统工程与电子技术》 *
葛凤翔 等: ""基于特征分析和二次规划的窄带信号超分辨率频率估计"", 《电子学报》 *

Also Published As

Publication number Publication date
CN113219243B (zh) 2021-12-31

Similar Documents

Publication Publication Date Title
CN109471095B (zh) 一种基于快速迭代插值的fmcw雷达距离估计方法
CN111610503B (zh) 基于改进的lvd的线性调频信号参数估计方法
CN112162152B (zh) 基于相位直线拟合的正弦波相参脉冲串信号频率估计方法
CN112881796A (zh) 频谱泄漏校正的多频实信号频率估计算法
CN109324309A (zh) 一种旋转单天线测角装置及其测量方法
CN110702042B (zh) 一种超声脉冲回波测厚中信号混叠条件下的厚度估算方法
CN109597137B (zh) 基于半导体磁传感器的Overhauser磁力仪跟踪配谐方法
CN113219243B (zh) 一种校正频率搜索间隔量化误差的高精度测频方法
CN114265050A (zh) 一种脉冲压缩和差测距方法
CN105652254A (zh) 室外场rcs测量方法及系统
CN114035149B (zh) 一种提升干涉仪测向系统灵敏度的方法
CN112541157B (zh) 一种信号频率精确估计方法
CN110441748A (zh) 一种基于幅度信息的α-β滤波方法
Belega et al. Estimation of the multifrequency signal parameters by interpolated DFT method with maximum sidelobe decay
Attivissimo et al. A study on nonlinear averagings to perform the characterization of power spectral density estimation algorithms
Chen et al. A frequency estimator for real valued sinusoidal signals using three dft samples
Belega et al. Multipoint interpolated DFT method for frequency estimation
CN112541156B (zh) 一种信号幅值精确估计方法
CN110133738A (zh) 基于IpDFT的质子磁力仪自由感应衰减信号的频率估计方法
CN112162153A (zh) 基于相位直线拟合的正弦波信号频率估计方法
Dao et al. Research on Improved Algorithm of Frequency Estimation Based on Complex Modulation.
CN110426610A (zh) 基于幅频响应和最小二乘法的谐振频率提取方法及系统
CN114252699B (zh) 一种离散傅里叶变换比幅测频方法
CN117347712A (zh) 一种基于线性调频z变换的信号频率估计方法
CN112485521B (zh) 一种基于对称离散傅里叶变换的频率估计方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant