CN109471095B - 一种基于快速迭代插值的fmcw雷达距离估计方法 - Google Patents

一种基于快速迭代插值的fmcw雷达距离估计方法 Download PDF

Info

Publication number
CN109471095B
CN109471095B CN201811310159.8A CN201811310159A CN109471095B CN 109471095 B CN109471095 B CN 109471095B CN 201811310159 A CN201811310159 A CN 201811310159A CN 109471095 B CN109471095 B CN 109471095B
Authority
CN
China
Prior art keywords
target
spectral
deviation
frequency
unit
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
CN201811310159.8A
Other languages
English (en)
Other versions
CN109471095A (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.)
Harbin Engineering University
Original Assignee
Harbin Engineering University
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 Harbin Engineering University filed Critical Harbin Engineering University
Priority to CN201811310159.8A priority Critical patent/CN109471095B/zh
Publication of CN109471095A publication Critical patent/CN109471095A/zh
Application granted granted Critical
Publication of CN109471095B publication Critical patent/CN109471095B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/06—Systems determining position data of a target
    • G01S13/08—Systems for measuring distance only
    • G01S13/32—Systems for measuring distance only using transmission of continuous waves, whether amplitude-, frequency-, or phase-modulated, or unmodulated

Landscapes

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

Abstract

本发明涉及一种基于快速迭代插值的FMCW雷达距离估计算法，本发明针对多目标环境中车载防撞雷达的测距精度易受到频谱间干扰的问题，首先利用快速傅里叶变换并结合恒虚警检测算法估计出目标的个数及目标的粗略位置，然后对目标的频谱进行非均匀迭代插值，并在每次迭代的过程中消除频谱间干扰，从而得到更精确的距离。仿真验证所提出的算法在多目标环境中和低信噪比条件下均具有更高的估计精度。最后对该算法进行DSP实现，证明该算法满足实时性的要求，易于工程实现。

Description

一种基于快速迭代插值的FMCW雷达距离估计方法
技术领域
本发明涉及一种FMCW雷达距离估计算法，特别是一种基于快速迭代插值的FMCW雷达距离估计算法，本发明属于毫米波雷达技术领域和信号处理领域，特别涉及多目标之间的干扰。
背景技术
调频连续波(Frequency Modulated Continuous Wave,FMCW)雷达具有良好的测速测距性能，因此被广泛应用于无人驾驶，医疗监护，信道估计等方面。然而传统的FMCW雷达存在栅栏效应和频谱泄露的问题，导致距离的估计偏差较大。为了更准确的估计目标距离，目前常用的方法有：采样序列直接补零法，基于复调制的局部放大傅里叶变换法(ZoomFast Fourier Transformation,ZFFT)方法、傅里叶变换-离散时间傅里叶变换联合校正法(Fast Fourier Transformation-Discrete-time Fourier Transform,FFT-DTFT)、Chirp-Z变换等。但以上四种方法存在计算量大，估计误差大等缺点。
近些年来，一些基于最大似然估计的频率估计方法被提出，这些方法同样可被用于距离估计。文献(提高LFMCW雷达测距精度的改进型ZFFT算法,2015,(2):1-4.)提出了一种单目标信号下的迭代插值算法，利用插值结果在迭代的过程中逐步修正偏差；文献(Accurate Frequency Estimators,2007,24(1):123-125.)提出了一种基于离散傅里叶变换(Discrete Fourier Transform，DFT)样本的频率偏差校正方法，该算法无需迭代，直接利用与目标的临近的谱峰对偏差进行估计；在此基础上，文献(A Method For FineResolution Frequency Estimation From Three DFT Samples,IEEE Signal ProcessingLetters,2011,18(6):351-354.)增加了偏差校正因子，提高了精度，但存在估计误差不均匀的问题；文献(基于细化频谱的频率迭代插值估计算法,电子与信息学报, 2017,39(9):2141-2147.)提出了利用细化频谱对信号进行插值的方法，从而解决了估计误差不均匀的问题。以上文献提出的方法对单频信号的偏差估计具有良好的效果，但是应用在多目标环境时，频谱间干扰会造成距离估计的偏差增大。为了解决此问题，文献(DFT InterpolationAlgorithm for Kaiser–Bessel and Dolph–Chebyshev Windows,IEEE Transactions onInstrumentation&Measurement,2011,60(3):784-790.)提出了基于高阶多项式插值补偿估计偏差的方法，该方法利用高阶多项式拟合出信号的频谱，获得每个目标的精确频率值，但该算法计算量大，不利于硬件实现。文献(Fine resolution frequency estimationfrom three DFT samples: Case of windowed data,Signal Processing,2015,114(1):245-250.)提出了加窗抑制频谱干扰的方法，利用窗函数抑制旁瓣的干扰,但该文献所提出的方法存在频谱展宽影响估计性能的问题。在DOA估计领域中，文献(Fast iterativeinterpolated beamforming for accurate single-snapshot DOA estimation,IEEEGeoscience and Remote Sensing Letters,2017,14(4):574-578.)提出利用迭代逐步校正角度偏差的方法，可以避免频谱展宽的问题，同时还可以解决多目标角度之间的干扰，但该方法应用在FMCW雷达***存在以下问题：在粗估的过程中需要把目标个数作为先验信息，且对每个目标谱峰估计时都需要对原始数据修正，增加了算法的复杂度。
受迭代校正思想解决多目标情况下频谱干扰的启发，本发明提出一种改进的基于快速迭代插值的FMCW雷达距离估计算法，通过恒虚警检测的方法粗估出目标的位置解决目标个数作为先验信息的问题，利用迭代过程中修正偏差因子替代对原始数据重复修正的过程。本发明用于解决在多目标环境中频谱间干扰的问题，具体说是一种可以消除多目标频谱间干扰并得到更精确目标反射信号频率及目标距离的方法。简化了处理过程，更利于硬件实现。理论仿真和实测数据验证了算法的可行性。
发明内容
针对上述现有技术，本发明要解决的技术问题是提供一种能够克服多目标环境中频谱之间的干扰对频率估计的影响，从而利用由目标本身的反射信号估计目标的频率偏差的基于快速迭代插值的FMCW雷达距离估计算法。
为解决上述技术问题，本发明一种基于快速迭代插值的FMCW雷达距离估计算法，包括以下步骤：
步骤1：对AD采样之后的中频信号进行傅里叶变换，然后进行恒虚警检测估计出目标个数P及目标所在的频谱单元mi,i＝1,2...P；
步骤2：计算第i个目标的频谱单元mi对应的频谱单元偏差
Figure BDA0001854722580000021
i＝1,2…P，步骤2包括以下步骤：
步骤2.1：设定目标迭代次数Q,并令当前迭代次数q＝0；对单次迭代之后第i个目标的偏差估计值
Figure BDA0001854722580000022
及幅值估计量
Figure BDA0001854722580000023
进行初始化，令
Figure BDA0001854722580000024
i＝1,2…P；
步骤2.2：估计第p个目标的反射信号强度
Figure BDA0001854722580000025
p＝1,2...P；
步骤2.3：利用
Figure BDA0001854722580000026
计算第p个目标对第i个目标产生频谱干扰
Figure BDA0001854722580000027
步骤2.4：计算去除干扰后
Figure BDA0001854722580000028
频谱单元的谱峰
Figure BDA0001854722580000029
步骤2.5：利用谱峰
Figure BDA00018547225800000210
计算偏差校正因子
Figure BDA00018547225800000211
和频谱单元偏差
Figure BDA00018547225800000212
步骤2.6：q←q+1，判断q是否小于Q，若q≤Q，则重复步骤2.2～2.5；若q＞Q，执行步骤3；
步骤3：利用频谱单元mi及频谱单元偏差
Figure BDA0001854722580000031
计算第i个目标反射信号的频率fi及目标的距离值Ri。
本发明还包括：
1.步骤1中对AD采样之后的中频信号进行傅里叶变换表达式如下：
Figure BDA0001854722580000032
其中sb(n)表示AD采样之后的中频信号，N表示傅里叶变换的点数，Sb(k)表示频域内第 k个单元的幅值，估计出目标个数P的表达式满足：
P＝Count{Sb(k)＞CFAR(Sb(k))}
其中CFAR(Sb(k))表示对Sb(k)恒虚警检测的结果,Count表示对Sb(k)大于CFAR(Sb(k))的幅值个数进行统计，统计的结果即为目标的个数；
第p个目标所在的频谱单元mp表达式如下：
Figure BDA0001854722580000033
2.步骤2.1中迭代次数Q的设置为2。
3.步骤2.2中第p个目标的反射信号强度
Figure BDA0001854722580000034
的估计值表达式满足：
Figure BDA0001854722580000035
其中
Figure BDA0001854722580000036
是第l个目标的反射信号强度估计值，其初始值为0，Δmpl表示第p个目标对应的频谱单元mp与第l个目标对应的频谱单元ml在频域内的距离，表达式满足：
Δmpl＝Δmp-Δml
Δδpl表示第p个目标对应的频谱单元偏差
Figure BDA0001854722580000037
与第l个目标对应的频谱单元偏差
Figure BDA0001854722580000038
在频域内的距离，表达式如下：
Δδpl＝Δδp-Δδl。
4.步骤2.3中第p个目标对第i个目标产生频谱干扰
Figure BDA0001854722580000039
表达式满足：
Figure BDA0001854722580000041
其中Δmip表示第i个目标对应的频谱单元mi与第p个目标对应的频谱单元的mp在频域内的距离，表达式满足：
Δmip＝Δmi-Δmp
Δδip表示第i个目标对应的频谱单元偏差
Figure BDA0001854722580000042
与第p个目标对应的频谱单元偏差
Figure BDA0001854722580000043
在频域内的距离，表达式满足：
Δδip＝Δδi-Δδp。
5.步骤2.4中去除干扰后
Figure BDA0001854722580000044
频谱单元的谱峰
Figure BDA0001854722580000045
表达式满足：
Figure BDA0001854722580000046
其中Si,±0.5表示对中频信号直接傅里叶变换得到
Figure BDA0001854722580000047
频谱单元的谱峰，Si,±0.5由两部分组成：第i个目标反射信号在
Figure BDA0001854722580000048
频谱单元做傅里叶变换得到的谱峰及第p个目标对第i个目标产生频谱干扰，表达式满足：
Figure BDA0001854722580000049
6.步骤2.5中偏差校正因子
Figure BDA00018547225800000410
和频谱单元偏差
Figure BDA00018547225800000411
表达式分别满足：
Figure BDA00018547225800000412
Figure BDA00018547225800000413
7.步骤3中第i个目标反射信号的频率fi及目标的距离值Ri表达式分别满足：
Figure BDA00018547225800000414
其中fs是AD采样频率，c是光速，B是调频带宽，T调频三角波的周期。
本发明有益效果：
本发明首先利用快速傅里叶变换并结合恒虚警检测算法估计出目标的个数及目标的粗略位置，然后对目标的频谱进行非均匀迭代插值，并在每次迭代的过程中消除频谱间干扰，从而得到更精确的距离。
附图说明
图1是FMCW雷达***框图；
图2是目标1在不同信噪比下的距离均方根误差；
图3是目标2在不同信噪比下的距离均方根误差；
图4是目标2采样点对距离均方根误差的影响；
图5是多目标下不同算法的距离均方根误差；
图6是目标3局部细化的结果；
图7是目标8局部细化的结果；
图8是实际环境测试图；
图9是FMCW雷达***；
图10是目标真实距离与FMCW雷达输出结果。
具体实施方式
为了便于描述，令Si＝Sb(mi)，Sb(mi)表示频域内第mi个频谱单元的谱峰，则
Figure BDA0001854722580000051
个频谱单元的谱峰可表示为：
Figure BDA0001854722580000052
实现本发明目的具体实施方式：
步骤1：对AD采样之后的中频信号进行傅里叶变换(Fast FourierTransformation,FFT)，然后进行恒虚警检测(Constant False-Alarm Rate,CFAR)估计出目标个数P及目标所在的频谱单元mi,i＝1,2...P。
傅里叶变换表达式如下
Figure BDA0001854722580000053
其中sb(n)表示AD采样之后的中频信号，N表示FFT变换的点数，Sb(k)表示频域内第k 个单元的幅值。通过对Sb(k)做CFAR处理估计出背景噪声的功率水平，然后将Sb(k)与背景噪声功率水平比较，估计出目标个数P，表达式如下
P＝Count{Sb(k)＞CFAR(Sb(k))}
其中CFAR(Sb(k))表示对Sb(k)恒虚警检测的结果，即背景噪声功率,Count表示对Sb(k)大于CFAR(Sb(k))的谱峰个数进行统计，统计的结果即为目标的个数。另外，还需要提取目标所在的频谱单元，并作为估计频谱单元偏差
Figure BDA0001854722580000061
基础。
第p个目标所在的频谱单元mp表达式如下
Figure BDA0001854722580000062
步骤2：计算第i个目标的频谱单元mi对应的频谱单元偏差
Figure BDA0001854722580000063
i＝1,2…P。
步骤2.1：设定目标迭代次数Q,并令当前迭代次数q＝0；对单次迭代之后第i个目标的偏差估计值
Figure BDA0001854722580000064
及第i个目标的反射信号强度
Figure BDA0001854722580000065
进行初始化，令
Figure BDA0001854722580000066
i＝1,2…P。
迭代次数Q的设置为2即可，实验表明当Q取2即可达到很高的距离精度，再进行多次迭代时对距离估计精度影响较小。再进行多次迭代时只会增加计算量对距离估计精度影响较小。
步骤2.2：估计第p个目标的反射信号强度
Figure BDA0001854722580000067
p＝1,2...P。
直接利用中频信号在
Figure BDA0001854722580000068
频谱单元做FFT变换得到的频谱峰值包含三部分：第p个目标反射信号进行FFT变换得到幅值、其他目标反射信号做FFT之后在
Figure BDA0001854722580000069
频谱单元的泄露和噪声对于幅值的干扰，其中噪声的干扰对谱峰影响较小可忽略。第l个目标在
Figure BDA00018547225800000610
频谱单元的泄露可通过该目标反射信号在
Figure BDA00018547225800000611
频谱单元的位置做FFT得到，表达式如下
Figure BDA00018547225800000612
其中
Figure BDA00018547225800000613
是第l个目标的反射信号强度估计值，其初始值为0，sbl(n)表示第l个目标反射得到的中频信号。由于sbl(n)无法得到，因此
Figure BDA00018547225800000614
也通过迭代逐次更新，并不断趋近真实值的过程。Δmpl表示第p个目标对应的频谱单元mp与第l个目标对应的频谱单元的ml在频域内的距离，表达式如下
Δmpl＝Δmp-Δml
Δδpl表示第p个目标对应的频谱单元偏差
Figure BDA00018547225800000615
与第l个目标对应的频谱单元偏差
Figure BDA00018547225800000616
在频域内的距离，表达式如下
Δδpl＝Δδp-Δδl
由此可得第p个目标的反射信号强度
Figure BDA0001854722580000071
表达式如下：
Figure BDA0001854722580000072
步骤2.3：利用
Figure BDA0001854722580000073
计算第p个目标对第i个目标产生频谱干扰
Figure BDA0001854722580000074
第p个目标对第i个目标产生频谱干扰
Figure BDA0001854722580000075
可通过第p个目标中频信号sbp(n)在
Figure BDA0001854722580000076
频谱单元做FFT变换得到，但sbp(n)无法直接得到，因此需要利用
Figure BDA0001854722580000077
求取
Figure BDA0001854722580000078
表达式如下
Figure BDA0001854722580000079
其中Δmip表示第i个目标对应的频谱单元mi与第p个目标对应的频谱单元的mp在频域内的距离，表达式如下
Δmip＝Δmi-Δmp
Δδip表示第i个目标对应的频谱单元偏差
Figure BDA00018547225800000710
与第p个目标对应的频谱单元偏差
Figure BDA00018547225800000711
在频域内的距离，表达式如下
Δδip＝Δδi-Δδp
步骤2.4：计算去除干扰后
Figure BDA00018547225800000712
频谱单元的谱峰
Figure BDA00018547225800000713
对中频信号直接做FFT可得到谱峰Si,±0.5，Si,±0.5由两部分组成：第i个目标反射信号在
Figure BDA00018547225800000714
频谱单元做FFT得到的谱峰及第p个目标反射信号做FFT在
Figure BDA00018547225800000715
频谱单元的增强，由此可知
Figure BDA00018547225800000716
的表达式
Figure BDA00018547225800000717
步骤2.5：利用谱峰
Figure BDA0001854722580000081
计算偏差校正因子
Figure BDA0001854722580000082
和频谱单元偏差
Figure BDA0001854722580000083
由步骤2.4得到的谱峰
Figure BDA0001854722580000084
可近似为反射信号中只存在第i个目标时在
Figure BDA0001854722580000085
频谱单元的谱峰。因此可利用单目标频率估计(Iterative frequencyestimation by interpolation on Fourier coefficients[J].IEEE Transactions onSignal Processing,2005,53(4):1237-1242.)的方法得到第i 个目标的偏差校正因子
Figure BDA0001854722580000086
和频谱单元偏差
Figure BDA0001854722580000087
表达式如下
Figure BDA0001854722580000088
Figure BDA0001854722580000089
步骤2.5：q←q+1，判断q是否小于Q，若q≤Q，则重复步骤2.2～2.5
当前迭代次数为q每次迭代结束后更新q并判断q与设定迭代次数Q的关系。
步骤3：利用频谱单元mi及频谱单元偏差
Figure BDA00018547225800000810
计算第i个目标反射信号的频率fi及目标的距离值Ri
当得到频谱单元偏差
Figure BDA00018547225800000811
之后，即可之后第i个目标反射信号对应的真实频谱单元
Figure BDA00018547225800000812
从而利用
Figure BDA00018547225800000813
计算第i个目标对应的频率值fi，由fi与距离的关系式得到第i个目标的距离值Ri，表达式如下
Figure BDA00018547225800000814
其中fs是AD采样频率，c是光速，B是调频带宽，T调频三角波的周期。
为了验证算法的性能，本文选择三种典型的方法进行对比，改进型局部放大傅里叶算法 (Zoom-FFT)(详见提高LFMCW雷达测距精度的改进型ZFFT算法,2015,(2):1-4.)，基于离散傅里叶变换的非迭代频率估计方法(Candan算法)(详见A Method For FineResolution Frequency Estimation From Three DFT Samples[J].IEEE SignalProcessing Letters,2011,18(6):351-354.)，基于汉明窗频谱泄露抑制算法(Hanming-window算法)(详见Fine resolution frequency estimation from three DFT samples:Case of windowed data[J].Signal Processing,2015,114:245-250.)。主要给出了不同信噪比(Signal-to-Noise Ratio,SNR)条件下的距离均方根误差(Root Mean SquareError,RMSE)及多目标下不同算法的RMSE等仿真效果。
本发明的效果可通过以下仿真说明：
(一)仿真条件与内容：
仿真实验采用fc为24.125GHz，带宽B为250MHz，信号周期T为8us，采样率为95KHz，FFT点数N为512。图2、3中两目标的信号强度满足A2＝0.5A1。图4中两目标的信号强度满足A2＝0.5A1，且SNR＝15dB。表1是信号强度变化范围是A∈[0.3,1],其中A＝1表示目标反射时无信号损失。仿真实验均采用5000次Monte-Carlo仿真。
从图2可看出，相比于其它三种算法，本文提出的算法在同样的SNR下具有更好的估计性能。此外，Zoom-FFT算法RMSE曲线趋于平坦，这是因为该算法的估计精度受局部细化倍数的限制，导致距离的偏差值不随SNR变化。
从图3可看出，本文提出的算法及Candan算法对目标2的偏差估计都具有较好的性能，这是因为，两种算法都是通过目标临近的谱峰对偏差估计，左右两个临近谱峰受目标1干扰引起的峰值的变化几乎相同，偏差校正因子
Figure BDA0001854722580000091
的求取受峰值变化较小，因此不会影响偏差
Figure BDA0001854722580000092
的估计，因此两种算法的估计性能较好。另外，Hanming-window算法在目标1的干扰下影响了偏差校正因子的估计，导致距离估计失败。
图4是目标2的信号采样点数对RMSE的影响。由图2可知Hanming-window算法和Zoom-FFT算法在该条件下对偏差估计误差较大，因此选择Candan算法作对比。从图4中可以看出，目标1和目标2的仿真结果相同，因此只给出了目标2的仿真结果。从图中可以看出，本文提出的算法在迭代次数Q较少时就可获得较小的RMSE，且在相同的采样点数下 RMSE更小。
表1是10个目标下的仿真结果，目标的信号强度变化范围是A∈[0.3,1],其中A＝1表示目标反射时无信号损失；从表中可看出信号幅值对本文提出的算法估计性能影响很小，另外从表中可看出目标6、7、8之间的距离差距不足0.5m，算法对每个目标的估计误差都很小，这是因为算法在迭代的过程中抑制了频谱干扰对偏差估计的影响。
表1多目标的仿真参数及结果
Figure BDA0001854722580000101
图5是多种算法对表1中目标进行仿真的RMSE。从图中可看出，本文提出的算法估计误差稳定在10cm以下，相比于其它两种算法，估计误差更均匀，RMSE更小，结合表1中信号强度发现本文提出的算法估计偏差受到信号强度的影响较小，这是因为，信号的干扰引起的频谱的变化均可以通过本文所提出的算法抑制掉，利用无干扰的频谱估计出的距离误差更小。从图中还可看出Zoom-FFT算法估计误差表现为随机性误差，这是因为Zoom-FFT的插值方式是均匀插值，当目标的谱峰离插值点越近误差越小，当目标的谱峰恰好落在插值点之间时误差最大。
图6和图7是对目标3、目标8局部进行放大的结果，所使用的算法是本文提出的算法及Zoom-FFT算法，从图中可知，通过改变迭代次数Q可以提高精确度，与Zoom-FFT算法相比，所提出的算法更接近真实频率，精确度更高。
图8和图9分别是实际环境测试图与FMCW雷达硬件***，目标1保持静止在14.6m处，目标2在[8,10]之间，每隔0.2m移动一次，FMCW雷达对两目标探测的测试结果及目标的真实值如图10所示，其中Rreal表示目标的真实距离值，Rest是使用本文提出的算法估计出的目标距离。
从图10中可看出，FMCW雷达输出的结果非常接近目标的真实距离且误差保持在10cm 以内。另外，FMCW雷达数据每秒可更新20次，证明该算法完全达到实时性的要求。
本发明具体实施方式还包括：
步骤1：对AD采样之后的中频信号进行傅里叶变换(Fast FourierTransformation,FFT)，然后进行恒虚警检测(Constant False-Alarm Rate,CFAR)估计出目标个数P及目标所在的频谱单元mi,i＝1,2...P。
步骤2：计算第i个目标的频谱单元mi对应的频谱单元偏差
Figure BDA0001854722580000111
i＝1,2…P。
步骤2.1：设定目标迭代次数Q,并令当前迭代次数q＝0；对单次迭代之后第i个目标的偏差估计值
Figure BDA0001854722580000112
及幅值估计量
Figure BDA0001854722580000113
进行初始化，令
Figure BDA0001854722580000114
i＝1,2…P。
步骤2.2：估计第p个目标的反射信号强度
Figure BDA0001854722580000115
p＝1,2...P；
步骤2.3：利用
Figure BDA0001854722580000116
计算第p个目标对第i个目标产生频谱干扰
Figure BDA0001854722580000117
步骤2.4：计算去除干扰后
Figure BDA0001854722580000118
频谱单元的谱峰
Figure BDA0001854722580000119
步骤2.5：利用谱峰
Figure BDA00018547225800001110
计算偏差校正因子
Figure BDA00018547225800001111
和频谱单元偏差
Figure BDA00018547225800001112
步骤2.6：q←q+1，判断q是否小于Q，若q≤Q，则重复步骤2.2～2.5。
步骤3：利用频谱单元mi及频谱单元偏差
Figure BDA00018547225800001113
计算第i个目标反射信号的频率fi及目标的距离值Ri
步骤1处理过中对AD采样之后的中频信号进行傅里叶变换表达式如下
Figure BDA00018547225800001114
其中sb(n)表示AD采样之后的中频信号，N表示FFT变换的点数，Sb(k)表示频域内第k 个单元的幅值。估计出目标个数P的表达式如下
P＝Count{Sb(k)＞CFAR(Sb(k))}
其中CFAR(Sb(k))表示对Sb(k)恒虚警检测的结果,Count表示对Sb(k)大于CFAR(Sb(k))的幅值个数进行统计，统计的结果即为目标的个数。
第p个目标所在的频谱单元mp表达式如下
Figure BDA00018547225800001115
步骤2.1中迭代次数Q的设置为2即可，实验表明当Q取2即可达到很高的距离精度，再进行多次迭代时对距离估计精度影响较小。
步骤2.2中第p个目标的反射信号强度
Figure BDA00018547225800001116
的估计值表达式如下
Figure BDA0001854722580000121
其中
Figure BDA0001854722580000122
是第l个目标的反射信号强度估计值，其初始值为0，Δmpl表示第p个目标对应的频谱单元mp与第l个目标对应的频谱单元ml在频域内的距离，表达式如下
Δmpl＝Δmp-Δml
Δδpl表示第p个目标对应的频谱单元偏差
Figure BDA0001854722580000123
与第l个目标对应的频谱单元偏差
Figure BDA0001854722580000124
在频域内的距离，表达式如下
Δδpl＝Δδp-Δδl
步骤2.3中第p个目标对第i个目标产生频谱干扰
Figure BDA0001854722580000125
表达式如下
Figure BDA0001854722580000126
其中Δmip表示第i个目标对应的频谱单元mi与第p个目标对应的频谱单元的mp在频域内的距离，表达式如下
Δmip＝Δmi-Δmp
Δδip表示第i个目标对应的频谱单元偏差
Figure BDA0001854722580000127
与第p个目标对应的频谱单元偏差
Figure BDA0001854722580000128
在频域内的距离，表达式如下
Δδip＝Δδi-Δδp
步骤2.4中去除干扰后
Figure BDA0001854722580000129
频谱单元的谱峰
Figure BDA00018547225800001210
表达式如下
Figure BDA00018547225800001211
其中Si,±0.5表示对中频信号直接FFT得到
Figure BDA00018547225800001212
频谱单元的谱峰，Si,±0.5由两部分组成：第i个目标反射信号在
Figure BDA00018547225800001213
频谱单元做FFT得到的谱峰及第p个目标对第i个目标产生频谱干扰，表达式如下
Figure BDA0001854722580000131
步骤2.5中偏差校正因子
Figure BDA0001854722580000132
和频谱单元偏差
Figure BDA0001854722580000133
表达式如下
Figure BDA0001854722580000134
Figure BDA0001854722580000135
步骤3中第i个目标反射信号的频率fi及目标的距离值Ri表达式如下
Figure BDA0001854722580000136
其中fs是AD采样频率，c是光速，B是调频带宽，T调频三角波的周期
本发明针对多目标环境中车载防撞雷达的测距精度易受到频谱间干扰的问题，本发明提出一种改进的基于迭代插值的多目标距离估计算法。该方法首先利用快速傅里叶变换并结合恒虚警检测算法估计出目标的个数及目标的粗略位置，然后对目标的频谱进行非均匀迭代插值，并在每次迭代的过程中消除频谱间干扰，从而得到更精确的距离。仿真验证所提出的算法在多目标环境中和低信噪比条件下均具有更高的估计精度。最后对该算法进行DSP实现，证明该算法满足实时性的要求，易于工程实现。

Claims (8)

1.一种基于快速迭代插值的FMCW雷达距离估计方法，其特征在于，包括以下步骤：
步骤1：对AD采样之后的中频信号进行傅里叶变换，然后进行恒虚警检测估计出目标个数P及目标所在的频谱单元mi，i＝1，2...P；
步骤2：计算第i个目标的频谱单元mi对应的频谱单元偏差
Figure FDA0003915250360000011
i＝1，2...P，步骤2包括以下步骤：
步骤2.1：设定目标迭代次数Q，并令当前迭代次数q＝0；对单次迭代之后第i个目标的偏差估计值
Figure FDA0003915250360000012
及幅值估计量
Figure FDA0003915250360000013
进行初始化，令
Figure FDA0003915250360000014
步骤2.2：估计第p个目标的反射信号强度
Figure FDA0003915250360000015
p＝1，2...P；
步骤2.3：利用
Figure FDA0003915250360000016
计算第p个目标对第i个目标产生频谱干扰
Figure FDA0003915250360000017
步骤2.4：计算去除干扰后
Figure FDA0003915250360000018
频谱单元的谱峰
Figure FDA0003915250360000019
步骤2.5：利用谱峰
Figure FDA00039152503600000110
计算偏差校正因子
Figure FDA00039152503600000111
和频谱单元偏差
Figure FDA00039152503600000112
步骤2.6：q←q+1，判断q是否小于Q，若q≤Q，则重复步骤2.2～2.5；若q＞Q，执行步骤3；
步骤3：利用频谱单元mi及频谱单元偏差
Figure FDA00039152503600000113
计算第i个目标反射信号的频率fi及目标的距离值Ri。
2.根据权利要求1所述的一种基于快速迭代插值的FMCW雷达距离估计方法，其特征在于：步骤1中对AD采样之后的中频信号进行傅里叶变换表达式如下：
Figure FDA00039152503600000114
其中sb(n)表示AD采样之后的中频信号，N表示傅里叶变换的点数，Sb(k)表示频域内第k个单元的幅值，估计出目标个数P的表达式满足：
P＝Count{Sb(k)＞CFAR(Sb(k))}
其中CFAR(Sb(k))表示对Sb(k)恒虚警检测的结果，Count表示对Sb(k)大于CFAR(Sb(k))的幅值个数进行统计，统计的结果即为目标的个数；
第p个目标所在的频谱单元mp表达式如下：
Figure FDA00039152503600000115
3.根据权利要求1所述的一种基于快速迭代插值的FMCW雷达距离估计方法，其特征在于：步骤2.1中迭代次数Q的设置为2。
4.根据权利要求1所述的一种基于快速迭代插值的FMCW雷达距离估计方法，其特征在于：步骤2.2中第p个目标的反射信号强度
Figure FDA0003915250360000021
的估计值表达式满足：
Figure FDA0003915250360000022
其中
Figure FDA0003915250360000023
是第l个目标的反射信号强度估计值，其初始值为0，N表示傅里叶变换的点数，Δmpl表示第p个目标对应的频谱单元mp与第l个目标对应的频谱单元ml在频域内的距离，表达式满足：
Δmpl＝Δmp-Δml
Δδpl表示第p个目标对应的频谱单元偏差
Figure FDA0003915250360000024
与第l个目标对应的频谱单元偏差
Figure FDA0003915250360000025
在频域内的距离，表达式如下：
Δδpl＝Δδp-Δδl。
5.根据权利要求1所述的一种基于快速迭代插值的FMCW雷达距离估计方法，其特征在于：步骤2.3中第p个目标对第i个目标产生频谱干扰
Figure FDA0003915250360000026
表达式满足：
Figure FDA0003915250360000027
其中，N表示傅里叶变换的点数，Δmip表示第i个目标对应的频谱单元mi与第p个目标对应的频谱单元的mp在频域内的距离，表达式满足：
Δmip＝Δmi-Δmp
Δδip表示第i个目标对应的频谱单元偏差
Figure FDA0003915250360000028
与第p个目标对应的频谱单元偏差
Figure FDA0003915250360000029
在频域内的距离，表达式满足：
Δδip＝Δδi-Δδp。
6.根据权利要求1所述的一种基于快速迭代插值的FMCW雷达距离估计方法，其特征在于：步骤2.4中去除干扰后
Figure FDA00039152503600000210
频谱单元的谱峰
Figure FDA00039152503600000211
表达式满足：
Figure FDA0003915250360000031
其中Si，±0.5表示对中频信号直接傅里叶变换得到
Figure FDA0003915250360000032
频谱单元的谱峰，Si，±0.5由两部分组成：第i个目标反射信号在
Figure FDA0003915250360000033
频谱单元做傅里叶变换得到的谱峰及第p个目标对第i个目标产生频谱干扰，表达式满足：
Figure FDA0003915250360000034
其中，sb(n)表示AD采样之后的中频信号，N表示傅里叶变换的点数。
7.根据权利要求1所述的一种基于快速迭代插值的FMCW雷达距离估计方法，其特征在于：步骤2.5中偏差校正因子
Figure FDA0003915250360000035
和频谱单元偏差
Figure FDA0003915250360000036
表达式分别满足：
Figure FDA0003915250360000037
Figure FDA0003915250360000038
其中，N表示傅里叶变换的点数。
8.根据权利要求1所述的一种基于快速迭代插值的FMCW雷达距离估计方法，其特征在于：步骤3中第i个目标反射信号的频率fi及目标的距离值Ri表达式分别满足：
Figure FDA0003915250360000039
其中fs是AD采样频率，c是光速，B是调频带宽，T调频三角波的周期，N表示傅里叶变换的点数。
CN201811310159.8A 2018-11-06 2018-11-06 一种基于快速迭代插值的fmcw雷达距离估计方法 Active CN109471095B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811310159.8A CN109471095B (zh) 2018-11-06 2018-11-06 一种基于快速迭代插值的fmcw雷达距离估计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811310159.8A CN109471095B (zh) 2018-11-06 2018-11-06 一种基于快速迭代插值的fmcw雷达距离估计方法

Publications (2)

Publication Number Publication Date
CN109471095A CN109471095A (zh) 2019-03-15
CN109471095B true CN109471095B (zh) 2023-02-14

Family

ID=65666960

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811310159.8A Active CN109471095B (zh) 2018-11-06 2018-11-06 一种基于快速迭代插值的fmcw雷达距离估计方法

Country Status (1)

Country Link
CN (1) CN109471095B (zh)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109581052A (zh) * 2018-11-10 2019-04-05 中国人民解放军陆军勤务学院 一种迭代插值的实复转换频率估计方法
CN111103588B (zh) * 2019-12-31 2022-12-13 哈尔滨工程大学 一种利用信号能量的三角波多目标识别方法
CN111337917B (zh) * 2020-03-17 2023-04-28 重庆邮电大学 基于变步长插值迭代的fmcw雷达高精度距离估计方法
CN111505618B (zh) * 2020-04-14 2022-06-24 东南大学 一种适用于车载毫米波雷达领域的基于频率估计的去耦合校正测距方法
CN111521968B (zh) * 2020-05-22 2022-05-20 南京理工大学 基于目标空间分集的欠定doa估计方法
CN112526474B (zh) * 2020-11-23 2022-05-27 哈尔滨工程大学 基于全相位傅里叶变换的fmcw雷达距离速度联合估计方法
CN112731307B (zh) * 2020-12-19 2023-01-03 湖北工业大学 基于距离-角度联合估计的ratm-cfar检测器及检测方法
CN113835077B (zh) * 2021-11-23 2022-02-11 中国空气动力研究与发展中心计算空气动力研究所 基于变脉冲重复频率的搜索雷达目标检测方法及***

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102608606A (zh) * 2012-03-22 2012-07-25 河海大学 汽车盲区监控雷达准确检测多目标的有效方法
CN103823215A (zh) * 2014-03-03 2014-05-28 中国科学院电子学研究所 线性调频连续波雷达测距方法
CN103823216A (zh) * 2014-03-10 2014-05-28 中国科学院电子学研究所 一种调频连续波雷达***测距方法
KR101446439B1 (ko) * 2013-07-30 2014-10-06 재단법인대구경북과학기술원 Fmcw 레이더의 고정밀 주파수 추정 방법 및 고정밀 주파수 추정 장치
DE102013210256A1 (de) * 2013-06-03 2014-12-04 Robert Bosch Gmbh Interferenzunterdrückung bei einem fmcw-radar
CN108415010A (zh) * 2018-04-28 2018-08-17 西安电子科技大学 一种基于梯形lfmcw调制的雷达多目标检测方法
KR20180113238A (ko) * 2017-04-05 2018-10-16 주식회사화신 Ｆｍｃｗ 레이더 시스템의 파라미터 추정 장치 및 방법

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102608606A (zh) * 2012-03-22 2012-07-25 河海大学 汽车盲区监控雷达准确检测多目标的有效方法
DE102013210256A1 (de) * 2013-06-03 2014-12-04 Robert Bosch Gmbh Interferenzunterdrückung bei einem fmcw-radar
KR101446439B1 (ko) * 2013-07-30 2014-10-06 재단법인대구경북과학기술원 Fmcw 레이더의 고정밀 주파수 추정 방법 및 고정밀 주파수 추정 장치
CN103823215A (zh) * 2014-03-03 2014-05-28 中国科学院电子学研究所 线性调频连续波雷达测距方法
CN103823216A (zh) * 2014-03-10 2014-05-28 中国科学院电子学研究所 一种调频连续波雷达***测距方法
KR20180113238A (ko) * 2017-04-05 2018-10-16 주식회사화신 Ｆｍｃｗ 레이더 시스템의 파라미터 추정 장치 및 방법
CN108415010A (zh) * 2018-04-28 2018-08-17 西安电子科技大学 一种基于梯形lfmcw调制的雷达多目标检测方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于快速迭代插值的FMCW雷达距离估计算法;段永昌等;《通信学报》;20190930;全文 *

Also Published As

Publication number Publication date
CN109471095A (zh) 2019-03-15

Similar Documents

Publication Publication Date Title
CN109471095B (zh) 一种基于快速迭代插值的fmcw雷达距离估计方法
CN108710103B (zh) 基于稀疏阵列的强弱多目标超分辨测向与信源数估计方法
CN107153178B (zh) 外辐射源雷达参考信号含有多径干扰时的目标检测方法
KR101294681B1 (ko) 기상 신호 처리장치 및 그 처리방법
CN108398669B (zh) 一种基于无需预延迟处理的空时宽带自适应单脉冲测角方法
CN103616687A (zh) 分段线性估计的多项式拟合isar包络对齐方法
Vizitiu et al. The synthesis of some NLFM laws using the stationary phase principle
CN103760534A (zh) 一种星载sar数据的电离层色散效应校正方法
CN111007473B (zh) 基于距离频域自相关函数的高速微弱目标检测方法
CN104155653B (zh) 一种基于特征距离子空间的sar后向投影成像方法
CN114545342A (zh) 利用多通道侦察接收机的雷达脉冲信号参数测量方法
CN114545351A (zh) 基于距离频率轴反转变换和二阶wvd的机动目标相参检测方法及***
CN112130138B (zh) 一种超视距雷达的电离层相位污染校正方法
CN110658517B (zh) 一种基于不确定先验知识的降维稀疏stap方法及装置
CN117471397A (zh) 一种基于图信号处理的圆形阵列二维doa估计方法
CN115826004B (zh) 一种基于二维角度及时差联合的三星协同直接定位方法
Torres et al. The impact of signal processing on the range-weighting function for weather radars
CN111551931B (zh) 基于极化高分辨距离像的目标长度特征提取方法
Chen et al. Efficient FFT based multi source DOA estimation for ULA
CN116996137B (zh) 一种基于加权叠加的低信噪比宽带线性调频信号检测方法
CN116449304B (zh) 一种基于频率测量的sar发射脉冲到达时间测量方法
CN113702969B (zh) 基于自适应stft方法的微多普勒信号参数估计方法
CN111948451B (zh) 一种基于互谱相位加权平均法的调频连续波信号测向方法
CN117949945B (zh) 一种双天线毫米波雷达超分辨目标定位方法及装置
CN114609606B (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