CN113219243A - 一种校正频率搜索间隔量化误差的高精度测频方法 - Google Patents
一种校正频率搜索间隔量化误差的高精度测频方法 Download PDFInfo
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/02—Arrangements for measuring frequency, e.g. pulse repetition rate; Arrangements for measuring period of current or voltage
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/16—Spectrum 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),为:
其中,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,为:
S5、由频谱序列的第m个值y(m),即频谱序列在幅度谱谱峰位置的值,与信号采样的尾首差v,确定频率搜索间隔的量化误差校正值u,为:
其中,angle表示复数的相位角。
S6、由幅度谱谱峰位置m和频率搜索间隔的量化误差校正值u,确定频率测量f,为:
本发明的有益效果是:使用本发明提出的一种校正频率搜索间隔量化误差的高精度测频方法,不仅利用了幅度谱谱峰位置,还利用了频率搜索间隔的量化误差校正值,实现高精度测频。本发明提出的方法的计算量与快速傅里叶变换法的计算量是同一数量级的。
具体实施方式
下面结合实施例对本发明的实用性进行说明。
实施例
在本例中,信号采样数分四种情况,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):
其中,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:
S5、根据步骤S4得到的频谱序列在幅度谱谱峰位置的值y(m),与采样信号的尾首差v,确定频率搜索间隔的量化误差校正值u:
其中,angle表示复数的相位角;
S6、由幅度谱谱峰位置m和频率搜索间隔的量化误差校正值u,获得频率f:
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)
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的正弦信号频率估计方法 |
-
2021
- 2021-04-27 CN CN202110458637.5A patent/CN113219243B/zh active Active
Patent Citations (12)
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)
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 |