CN106093905B - 基于apFFT相位差的精密测量雷达频偏分类测频方法 - Google Patents

基于apFFT相位差的精密测量雷达频偏分类测频方法 Download PDF

Info

Publication number
CN106093905B
CN106093905B CN201610454526.6A CN201610454526A CN106093905B CN 106093905 B CN106093905 B CN 106093905B CN 201610454526 A CN201610454526 A CN 201610454526A CN 106093905 B CN106093905 B CN 106093905B
Authority
CN
China
Prior art keywords
frequency
apfft
phase
value
iffirst
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.)
Active
Application number
CN201610454526.6A
Other languages
English (en)
Other versions
CN106093905A (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.)
Chengdu Rainier Technology Co Ltd
Original Assignee
Chengdu Rainier Technology Co Ltd
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 Chengdu Rainier Technology Co Ltd filed Critical Chengdu Rainier Technology Co Ltd
Priority to CN201610454526.6A priority Critical patent/CN106093905B/zh
Publication of CN106093905A publication Critical patent/CN106093905A/zh
Application granted granted Critical
Publication of CN106093905B publication Critical patent/CN106093905B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Measuring Frequencies, Analyzing Spectra (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种基于apFFT相位差的精密测量雷达频偏分类测频方法，包括以下步骤：(1)将差频信号序列X0(m)做nΔt的延迟得到XN(m)，X0(m)频率为fIF，nΔt表示n个采样周期；(2)对X0(m)、XN(m)分别加双汉宁窗，做apFFT并求得其相位值；(3)对二者相位值求差并对其以2π求Mod得到相位差Δφ；(4)找出次大值位置，得到第一次估计的频率；(5)根据情况确定是否修正频率；(6)将新信号重复步骤(2)、(3)、(4)，得到二次测量结果fIFSECOND，将fIFSECOND减去rfs/3N得到最终频率。本发明利用apFFT“相位不变性”以及相位与频率的线性关系，在任意可测量的频率都有较好的测量精度，而且计算量小、抗噪性好，利于硬件实现。

Description

基于apFFT相位差的精密测量雷达频偏分类测频方法
技术领域
本发明涉及一种精密测量雷达测频方法，尤其涉及一种基于apFFT相位差的精密测量雷达频偏分类测频方法。
背景技术
目前，精密测量雷达测频测距方法已经很成熟了，测量方法越来越多。以前，测量雷达测频方法多基于FFT即快速傅氏变换算法，但该方法因相位会发生变化而逐渐被淘汰。apFFT是全相位FFT的简称，是对传统FFT进行改进而得到的，它具有FFT不具备的“相位不变性”等优良性能，所以现在测量雷达测频方法多基于apFFT。但是，目前采用的apFFT算法存在相位模糊的问题，从而导致测量精度不高、尤其在某些频率测量精度较低的问题。
发明内容
本发明的目的就在于为了解决上述问题而提供一种在任意频率都具有较高测量精度的基于apFFT相位差的精密测量雷达频偏分类测频方法。
本发明通过以下技术方案来实现上述目的：
一种基于apFFT相位差的精密测量雷达频偏分类测频方法，包括以下步骤：
(1)将差频信号序列X0(m)做nΔt的延迟得到XN(m)，X0(m)频率为fIF，nΔt表示n个采样周期；
(2)对X0(m)、XN(m)分别加双汉宁窗，做apFFT并求得其相位值；
(3)对二者相位值求差并对其以2π求Mod得到相位差Δφ；
(4)找出次大值位置，若在峰值左侧，则第一次估计得到的频率为：
fIFFIRST＝(k+Δφ-1)fs/N
其中，N为FFT点数，fs为采样频率，k为谱线峰值最大值所对应的位置序列号；
若在峰值右侧，则第一次估计得到的频率为：
fIFFIRST＝(k+Δφ)fs/N；
(5)若0.5fs/N≥|fIFFIRST-kfs/N|＞0.33fs/N，此时频偏较大，测量准确，直接将fIFFIRST作为测频结果，不需要修正；第二种情况：若|fIFFIRST-kfs/N|≤0.33fs/N或|fIFFIRST-kfs/N|＞0.5fs/N，则将X0(m)、XN(m)分别乘以次大值在峰值右侧时r＝1，次大值在峰值左侧时r＝-1；
(6)将新信号重复步骤(2)、(3)、(4)，得到二次测量结果fIFSECOND，将fIFSECOND减去rfs/3N得到最终频率为：
fIFFINAL＝fIFSECOND-rfs/3N。
优选地，所述步骤(4)、(5)中，N＝512，fs＝512khz。
本发明的有益效果在于：
本发明利用apFFT“相位不变性”以及相位与频率的线性关系，提出基于apFFT相位差的精密测量雷达频偏分类测频算法，对于频偏较小的频点，都能够通过频谱搬移使得频偏变大，能够很准确分清粗次大值方向，使得相位模糊补偿值的索引值估测准确，解决了传统apFFT相位差值算法在频偏绝对值较小时算法的缺陷，测频精度甚至可以达到传统apFFT相位差值算法的几十倍；本方法解决了“相位模糊”问题带来了测量精度不高的缺点，在任意可测量的频率都有较好的测量精度，而且计算量小、抗噪性好，利于硬件实现。
附图说明
图1-1是信噪比SNR＝5dB时利用本发明所述基于apFFT相位差的精密测量雷达频偏分类测频方法获得的均值误差绝对值与频率的对照关系示意图；
图1-2是信噪比SNR＝5dB时利用本发明所述基于apFFT相位差的精密测量雷达频偏分类测频方法获得的均方根误差与频率的对照关系示意图；
图2-1是信噪比SNR＝5dB时利用传统apFFT相位差测频方法获得的均值误差绝对值与频率的对照关系示意图；
图2-2是信噪比SNR＝5dB时利用传统apFFT相位差测频方法获得的均方根误差与频率的对照关系示意图；
图3-1是信噪比SNR＝3dB时利用本发明所述基于apFFT相位差的精密测量雷达频偏分类测频方法获得的均值误差绝对值与频率的对照关系示意图；
图3-2是信噪比SNR＝3dB时利用本发明所述基于apFFT相位差的精密测量雷达频偏分类测频方法获得的均方根误差与频率的对照关系示意图。
具体实施方式
下面结合附图和实施例对本发明作进一步说明：
本发明所述基于apFFT相位差的精密测量雷达频偏分类测频方法，包括以下步骤：
(1)将差频信号序列X0(m)做nΔt的延迟得到XN(m)，X0(m)频率为fIF，nΔt表示n个采样周期；
(2)对X0(m)、XN(m)分别加双汉宁窗，做apFFT并求得其相位值；
(3)对二者相位值求差并对其以2π求Mod得到相位差Δφ；
(4)找出次大值位置，若在峰值左侧，则第一次估计得到的频率为：
fIFFIRST＝(k+Δφ-1)fs/N
其中，N为FFT点数，fs为采样频率，k为谱线峰值最大值所对应的位置序列号；
若在峰值右侧，则第一次估计得到的频率为：
fIFFIRST＝(k+Δφ)fs/N；
(5)若0.5fs/N≥|fIFFIRST-kfs/N|＞0.33fs/N，此时频偏较大，测量准确，直接将fIFFIRST作为测频结果，不需要修正；第二种情况：若|fIFFIRST-kfs/N|≤0.33fs/N或|fIFFIRST-kfs/N|＞0.5fs/N，则将X0(m)、XN(m)分别乘以次大值在峰值右侧时r＝1，次大值在峰值左侧时r＝-1；
(6)将新信号重复步骤(2)、(3)、(4)，得到二次测量结果fIFSECOND，将fIFSECOND减去rfs/3N得到最终频率为：
fIFFINAL＝fIFSECOND-rfs/3N。
实施例：
图1-1至图3-2的仿真条件同为FFT点数为N＝512，采样频率设置为fs＝512khz，差拍信号频率设置为fIF∈[30000hz,31000hz]，两两相差10hz，做1000次蒙特卡洛。
由图1-1、图1-2可知，利用本发明所述基于apFFT相位差的精密测量雷达频偏分类测频方法获得的均值误差绝对值最大值不超过10，均方根误差不超过12.2；对应地，由图2-1、图2-2可知，利用传统apFFT相位差测频方法获得的均值误差绝对值最高可达300以上，均方根误差最大可达600以上；可以看出本发明所述基于apFFT相位差的精密测量雷达频偏分类测频方法的测量性能卓越。
由图3-1、图3-2可知，本发明所述基于apFFT相位差的精密测量雷达频偏分类测频方法应用在信噪比SNR＝3的环境中时，其测量特性相较于SNR＝5时会增大，但增大并不明显，相比于传统apFFT相位差测频方法仍然非常具有较高的测量精度。综上可以证明，本发明所述基于apFFT相位差的精密测量雷达频偏分类测频方法有着良好的测频性能。
上述实施例只是本发明的较佳实施例，并不是对本发明技术方案的限制，只要是不经过创造性劳动即可在上述实施例的基础上实现的技术方案，均应视为落入本发明专利的权利保护范围内。

Claims (2)

1.一种基于apFFT相位差的精密测量雷达频偏分类测频方法，其特征在于：包括以下步骤：
(1)将差频信号序列X0(m)做nΔt的延迟得到XN(m)，X0(m)频率为fIF，nΔt表示n个采样周期；
(2)对X0(m)、XN(m)分别加双汉宁窗，做apFFT并求得其相位值；
(3)对二者相位值求差并对其以2π求Mod得到相位差Δφ；
(4)找出次大值位置，若在峰值左侧，则第一次估计得到的频率为：
fIFFIRST＝(k+Δφ-1)fs/N
其中，N为FFT点数，fs为采样频率，k为谱线峰值最大值所对应的位置序列号；
若在峰值右侧，则第一次估计得到的频率为：
fIFFIRST＝(k+Δφ)fs/N；
(5)若0.5fs/N≥|fIFFIRST-kfs/N|＞0.33fs/N，此时频偏较大，测量准确，直接将fIFFIRST作为测频结果，不需要修正；第二种情况：若|fIFFIRST-kfs/N|≤0.33fs/N或|fIFFIRST-kfs/N|＞0.5fs/N，则将X0(m)、XN(m)分别乘以次大值在峰值右侧时r＝1，次大值在峰值左侧时r＝-1；
(6)将新信号重复步骤(2)、(3)、(4)，得到二次测量结果fIFSECOND，将fIFSECOND减去rfs/3N得到最终频率为：
fIFFINAL＝fIFSECOND-rfs/3N。
2.根据权利要求1所述的基于apFFT相位差的精密测量雷达频偏分类测频方法，其特征在于：所述步骤(4)、(5)中，N＝512，fs＝512kH z。
CN201610454526.6A 2016-06-21 2016-06-21 基于apFFT相位差的精密测量雷达频偏分类测频方法 Active CN106093905B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610454526.6A CN106093905B (zh) 2016-06-21 2016-06-21 基于apFFT相位差的精密测量雷达频偏分类测频方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610454526.6A CN106093905B (zh) 2016-06-21 2016-06-21 基于apFFT相位差的精密测量雷达频偏分类测频方法

Publications (2)

Publication Number Publication Date
CN106093905A CN106093905A (zh) 2016-11-09
CN106093905B true CN106093905B (zh) 2018-07-31

Family

ID=57237866

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610454526.6A Active CN106093905B (zh) 2016-06-21 2016-06-21 基于apFFT相位差的精密测量雷达频偏分类测频方法

Country Status (1)

Country Link
CN (1) CN106093905B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106644075A (zh) * 2016-11-17 2017-05-10 天津津航技术物理研究所 一种傅立叶光谱仪的高效去噪方法
CN113281566B (zh) * 2021-05-11 2023-11-14 重庆矩子兴智能科技有限公司 一种基于组合复信号相位差的频率估计方法
CN113985116A (zh) * 2021-09-14 2022-01-28 合肥工业大学 氧化锌避雷器泄漏电流的全相位fft时移相位差校正方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101136892A (zh) * 2007-10-10 2008-03-05 天津大学 基于全相位fft的ofdm调制解调方法
CN101388001A (zh) * 2008-06-25 2009-03-18 天津大学 基于全相位fft的高精度瞬间相位估计方法
CN103983957B (zh) * 2014-05-12 2016-09-07 天津大学 一种多普勒偏移测量方法及其装置
CN105547174B (zh) * 2015-11-27 2018-08-17 上海无线电设备研究所 分布式高精度激光在线测量系统

Also Published As

Publication number Publication date
CN106093905A (zh) 2016-11-09

Similar Documents

Publication Publication Date Title
CN1996986B (zh) 全相位时移相位差频谱校正法
CN108955734B (zh) 一种光纤f-p温度/压力复合传感器的腔长解调方法
CN105841762B (zh) 超声波水表的流量计量方法和系统
CN106093905B (zh) 基于apFFT相位差的精密测量雷达频偏分类测频方法
CN102539944B (zh) 一种基于鉴相法测量相位噪声的方法
CN105512431A (zh) 一种基于相位噪声数学模型的相位噪声测量方法
CN108414833B (zh) 一种信号分量频率的精确估计方法
CN103529293B (zh) 基于边沿效应的并行的频率和周期性信号参数测量方法
Huibin et al. Energy based signal parameter estimation method and a comparative study of different frequency estimators
CN114460527B (zh) 关联度延拓Hilbert移相电子式互感器校验仪溯源方法和系统
CN103926463B (zh) 一种数据延拓式相关的相位差测量方法
CN105307095B (zh) 一种基于fft的高分辨率音频频率测量方法
CN103969614A (zh) 一种数字万用表的校准方法
CN106154257B (zh) 基于FFT与apFFT的精密测量雷达二次测频方法
Zhang et al. Power system dynamic frequency measurement based on novel interpolated STFT algorithm
CN109030944B (zh) 一种基于频移和dft的实正弦信号相位差估计方法
CN102193029B (zh) 非常规采样时间短期频率稳定度测量方法
CN105785123B (zh) 一种基于apFFT相位差的雷达信号频率计算方法
Sun et al. Quadratic correlation time delay estimation algorithm based on Kaiser window and Hilbert transform
CN104062498A (zh) 对称频谱信号的中心频率的质心估计方法
CN107104699B (zh) 一种分数域中低运算复杂度的最优阶数搜索的装置和方法
CN105629254B (zh) 一种目标微动特征相干激光探测效果量化评价方法
CN105738878A (zh) 一种基于频偏分类的精密测量雷达测频方法
CN112595889A (zh) 非理想多指数衰减正弦信号的欠Nyquist采样与参数测量方法
CN105372493A (zh) 基于三条dft复数谱线的信号幅值和相位测量方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
PE01 Entry into force of the registration of the contract for pledge of patent right

Denomination of invention: Precision measurement radar frequency deviation classification frequency measurement method based on apFFT phase difference

Effective date of registration: 20200414

Granted publication date: 20180731

Pledgee: Chengdu SME financing Company Limited by Guarantee

Pledgor: CHENGDU RAINIER TECHNOLOGY Co.,Ltd.

Registration number: Y2020510000030

PE01 Entry into force of the registration of the contract for pledge of patent right
PC01 Cancellation of the registration of the contract for pledge of patent right

Date of cancellation: 20210512

Granted publication date: 20180731

Pledgee: Chengdu SME financing Company Limited by Guarantee

Pledgor: CHENGDU RAINIER TECHNOLOGY Co.,Ltd.

Registration number: Y2020510000030

PC01 Cancellation of the registration of the contract for pledge of patent right
PE01 Entry into force of the registration of the contract for pledge of patent right

Denomination of invention: Frequency offset classification method for precision measurement radar based on apFFT phase difference

Effective date of registration: 20210520

Granted publication date: 20180731

Pledgee: Chengdu SME financing Company Limited by Guarantee

Pledgor: CHENGDU RAINIER TECHNOLOGY Co.,Ltd.

Registration number: Y2021510000066

PE01 Entry into force of the registration of the contract for pledge of patent right