[go: up one dir, main page]

CN110208804A - Transmitted waveform design method and target acquisition algorithm suitable for automobile collision avoidance radar - Google Patents

Transmitted waveform design method and target acquisition algorithm suitable for automobile collision avoidance radar Download PDF

Info

Publication number
CN110208804A
CN110208804A CN201910473989.0A CN201910473989A CN110208804A CN 110208804 A CN110208804 A CN 110208804A CN 201910473989 A CN201910473989 A CN 201910473989A CN 110208804 A CN110208804 A CN 110208804A
Authority
CN
China
Prior art keywords
signal
frequency
lfsk
sweep
distance
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.)
Pending
Application number
CN201910473989.0A
Other languages
Chinese (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.)
Nanjing University of Science and Technology
Original Assignee
Nanjing University of Science and Technology
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 Nanjing University of Science and Technology filed Critical Nanjing University of Science and Technology
Priority to CN201910473989.0A priority Critical patent/CN110208804A/en
Publication of CN110208804A publication Critical patent/CN110208804A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/88Radar or analogous systems specially adapted for specific applications
    • G01S13/93Radar or analogous systems specially adapted for specific applications for anti-collision purposes
    • G01S13/931Radar or analogous systems specially adapted for specific applications for anti-collision purposes of land vehicles
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/88Radar or analogous systems specially adapted for specific applications
    • G01S13/93Radar or analogous systems specially adapted for specific applications for anti-collision purposes
    • G01S13/931Radar or analogous systems specially adapted for specific applications for anti-collision purposes of land vehicles
    • G01S2013/9315Monitoring blind spots

Landscapes

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

Abstract

本发明公开了一种适用于汽车防撞雷达的发射波形设计方法及目标探测算法,包括以下步骤:确定雷达发射信号的体制为多阶线性频率键控体制;根据给定的指标确定发射波形的参数,指标有:雷达发射频段、最大探测距离、最大探测速度、距离分辨力、测距精度、测速精度和频域最小信可检测噪比,确定一个波形需要的参数有:调频带宽、相参处理间隔、LFSK信号的数量、每个LFSK信号的步进数以及LFSK信号之间的起始频差;对发射信号的回波数据进行处理,求解目标的距离和速度。本发明计算量适中,在测距不模糊的基础上提高了探测精度,并有效减小了收发泄漏造成的探测盲区。

The invention discloses a transmission waveform design method and a target detection algorithm suitable for automobile anti-collision radars, comprising the following steps: determining that the system of the radar transmission signal is a multi-order linear frequency keying system; determining the transmission waveform according to a given index Parameters and indicators include: radar transmission frequency band, maximum detection distance, maximum detection speed, distance resolution, ranging accuracy, speed measurement accuracy and minimum signal detectable noise ratio in the frequency domain. The parameters required to determine a waveform are: FM bandwidth, coherent Processing interval, the number of LFSK signals, the number of steps of each LFSK signal, and the initial frequency difference between LFSK signals; process the echo data of the transmitted signal, and calculate the distance and speed of the target. The invention has a moderate amount of calculation, improves the detection accuracy on the basis of unambiguous distance measurement, and effectively reduces the detection blind area caused by the transmission and reception leakage.

Description

适用于汽车防撞雷达的发射波形设计方法及目标探测算法Transmission Waveform Design Method and Target Detection Algorithm Applicable to Automotive Anti-collision Radar

技术领域technical field

本发明属于雷达信息技术领域,特别是一种适用于汽车防撞雷达的发射波形设计方法及目标探测算法。The invention belongs to the technical field of radar information, in particular to a transmission waveform design method and a target detection algorithm suitable for automobile anti-collision radars.

背景技术Background technique

为了提高汽车驾驶的安全性能,减少交通事故的发生,高级驾驶辅助系统(Advanced Driver Assistance System.ADAS)成为研究的热点,它能够自动分析采集到的信息,判断是否存在安全隐患,以及是否需要发出相应的制动信号给汽车,使汽车拥有碰撞预警或自主避撞的能力。In order to improve the safety performance of automobile driving and reduce the occurrence of traffic accidents, Advanced Driver Assistance System (ADAS) has become a research hotspot. The corresponding braking signal is given to the car, so that the car has the ability of collision warning or autonomous collision avoidance.

毫米波雷达探测距离远、精度高,是ADAS传感器的一个可靠选择。目前的车载毫米波雷达以24GHz与77GHz为主,分别用于汽车的后向辅助与前方探测。在波形的选择上,现有的防撞雷达普遍采用线性调频连续波(Linear Frequency Modulated ContinuousWaveform.LFMCW),也有的系统采用多重频率键控(Multiple Frequency ShiftKeying.MFSK)体制,不同的发射波形带来不同的探测性能及处理复杂度。Millimeter wave radar has a long detection range and high accuracy, making it a reliable choice for ADAS sensors. The current vehicle-mounted millimeter-wave radars are mainly 24GHz and 77GHz, which are used for rearward assistance and forward detection of vehicles respectively. In terms of waveform selection, existing anti-collision radars generally use Linear Frequency Modulated Continuous Waveform (LFMCW), and some systems use Multiple Frequency ShiftKeying (Multiple Frequency ShiftKeying. Different detection performance and processing complexity.

线性调频连续波雷达利用不同扫频段对应的差拍信号频率解距离和速度的耦合,结构简单、发射功率低,在探测近距离目标时比较有优势。但在排除虚假目标、选择正确配对方面,该体制存在一些不足。经典的对称三角线性调频连续波体制无法进行目标配对,多目标时一般不采用;变调频斜率的线性调频连续波体制在目标配对时需要求解多个方程组,算法较为复杂,而且当目标数量增多时,其计算量和内存占用量会呈指数式增长;引入点频信号的线性调频连续波体制降低了目标配对的复杂度,但配对的准确度也随之降低;若改用动目标检测(Moving Targets Detection.MTD)求解目标的距离和速度,尽管省略了目标配对的步骤,但对于高速目标的探测而言,它对硬件的要求较高,增加了产品的成本。Linear frequency modulation continuous wave radar uses the beat signal frequency corresponding to different sweep frequency bands to solve the coupling of distance and speed. It has a simple structure and low transmission power, and has advantages in detecting short-range targets. However, there are some deficiencies in this system in terms of excluding false targets and selecting correct pairs. The classic symmetrical triangular LFM-CW system cannot perform target matching, and it is generally not used when there are multiple targets; the LFM-CW system with variable FM slope needs to solve multiple equations when pairing targets, and the algorithm is more complicated, and when the number of targets increases , the amount of computation and memory usage will increase exponentially; the introduction of the linear frequency modulation continuous wave system of point frequency signals reduces the complexity of target pairing, but the accuracy of pairing also decreases; if moving target detection ( Moving Targets Detection.MTD) solves the distance and speed of the target. Although the step of target pairing is omitted, it requires high hardware requirements for the detection of high-speed targets, which increases the cost of the product.

多重频率键控体制不用进行目标配对,它利用每个LFSK信号对应的差拍信号频率和相位差解距离、速度的耦合,完全排除了虚假目标的干扰,而且计算量相对较小,在多目标探测时具有优势。但探测的结果很大程度上取决于相位差的估计,由于相位差本身的取值范围较小,所以在信噪比不高或相参积累时长较短的情况下,相位差估计的小误差也会带来距离速度测量精度的较为明显的下降。The multiple frequency keying system does not need to perform target pairing. It uses the beat signal frequency and phase difference corresponding to each LFSK signal to solve the coupling of distance and speed, completely eliminating the interference of false targets, and the calculation amount is relatively small. Advantages when probing. However, the result of the detection largely depends on the estimation of the phase difference. Since the value range of the phase difference itself is small, the small error of the phase difference estimation is small when the SNR is not high or the coherent accumulation time is short. It will also bring about a more obvious decline in the accuracy of distance velocity measurement.

发明内容Contents of the invention

本发明的目的在于提供一种适用于汽车防撞雷达的发射波形设计方法及目标探测算法,在计算量适中的前提下提高目标探测的精度,同时减小收发泄漏造成的探测盲区。The purpose of the present invention is to provide a transmission waveform design method and a target detection algorithm suitable for automobile anti-collision radars, improve the accuracy of target detection under the premise of moderate calculation amount, and reduce the detection blind area caused by transmission and reception leakage.

实现本发明目的的技术解决方案为:一种适用于汽车防撞雷达的发射波形设计方法,包括以下步骤:The technical solution that realizes the object of the present invention is: a kind of transmission waveform design method that is applicable to automobile anti-collision radar, comprises the following steps:

步骤1,确定雷达发射信号的体制为MS-LFSK体制,一个周期的MS-LFSK波形由上扫频段和下扫频段构成,每个扫频段由若干LFSK信号合成;Step 1, determine that the system of the radar transmitting signal is the MS-LFSK system, a cycle of MS-LFSK waveform is composed of an up-sweep frequency segment and a down-sweep frequency segment, and each frequency scan segment is synthesized by several LFSK signals;

步骤2,根据给定的指标确定发射波形的参数,指标有:雷达发射频段、最大探测距离、最大探测速度、距离分辨力、测距精度、测速精度和频域最小信可检测噪比,确定一个周期的波形需要的参数有:调频带宽、相参处理间隔、LFSK信号的数量、每个LFSK信号的步进数以及LFSK信号之间的起始频差。Step 2. Determine the parameters of the transmission waveform according to the given indicators. The indicators include: radar transmission frequency band, maximum detection distance, maximum detection speed, distance resolution, ranging accuracy, speed measurement accuracy and minimum signal detectable noise ratio in the frequency domain. Determine The parameters required for a waveform of one cycle include: FM bandwidth, coherent processing interval, number of LFSK signals, number of steps for each LFSK signal, and initial frequency difference between LFSK signals.

一种适用于汽车防撞雷达的目标探测算法,包括以下步骤:A target detection algorithm suitable for automobile collision avoidance radar, comprising the following steps:

步骤1,确定雷达发射信号的体制为MS-LFSK体制,一个周期的MS-LFSK波形由上扫频段和下扫频段构成,每个扫频段由若干LFSK信号合成;Step 1, determine that the system of the radar transmitting signal is the MS-LFSK system, a cycle of MS-LFSK waveform is composed of an up-sweep frequency segment and a down-sweep frequency segment, and each frequency scan segment is synthesized by several LFSK signals;

步骤2,根据给定的指标确定发射波形的参数,指标有:雷达发射频段、最大探测距离、最大探测速度、距离分辨力、测距精度、测速精度和频域最小信可检测噪比,确定一个波形需要的参数有:调频带宽、相参处理间隔、LFSK信号的数量、每个LFSK信号的步进数以及LFSK信号之间的起始频差;Step 2. Determine the parameters of the transmission waveform according to the given indicators. The indicators include: radar transmission frequency band, maximum detection distance, maximum detection speed, distance resolution, ranging accuracy, speed measurement accuracy and minimum signal detectable noise ratio in the frequency domain. Determine The parameters required for a waveform are: FM bandwidth, coherent processing interval, number of LFSK signals, number of steps for each LFSK signal, and initial frequency difference between LFSK signals;

步骤3,对发射信号的回波数据进行处理,求解目标的距离和速度:将接收到的回波信号与本振混频,得到差拍信号,对其采样并按照回波数据对应LFSK信号的情况对它重排,对重排后的数据按扫频周期做FFT和周期图积累,并对得到的二维频谱进行OS-CFAR,得到的频谱峰值点就对应一个目标,根据峰值坐标可以计算出目标对应的差拍信号频率,再由相频响应计算得到峰值点在相邻LFSK信号处的相位差,综合差拍信号频率和相位差,再经过解模糊的处理,就可以求解出目标的距离和速度。Step 3, process the echo data of the transmitted signal to solve the distance and speed of the target: mix the received echo signal with the local oscillator to obtain the beat signal, sample it and correspond to the LFSK signal according to the echo data Rearrange it according to the situation, do FFT and periodogram accumulation on the rearranged data according to the sweep cycle, and perform OS-CFAR on the obtained two-dimensional spectrum, the obtained spectrum peak point corresponds to a target, and can be calculated according to the peak coordinates The frequency of the beat signal corresponding to the target is obtained, and then the phase difference of the peak point at the adjacent LFSK signal is obtained by calculating the phase-frequency response, and the frequency and phase difference of the beat signal are integrated, and after defuzzification processing, the target's distance and speed.

本发明与现有技术相比,其显著优点为:(1)本发明结合了线性调频连续波体制与多重频率键控体制的优点,又弥补了它们的不足;(2)在探测目标时,利用差拍信号的频率和相位差解距离和速度的模糊,没有虚假目标的干扰,不需要做目标配对,因此计算量适中;(3)发射波形由不同频率间隔的多组LFSK信号组合而成,在距离不模糊的基础上显著提高了测距测速的精度,测量更加准确;(4)探测目标时,利用上下扫频段协同工作,有效减小了探测盲区。Compared with the prior art, the present invention has the remarkable advantages as follows: (1) the present invention combines the advantages of the linear frequency modulation continuous wave system and the multiple frequency keying system, and makes up for their deficiencies; (2) when detecting the target, Using the frequency and phase difference of the beat signal to solve the ambiguity of distance and speed, there is no interference from false targets, and no target pairing is required, so the amount of calculation is moderate; (3) The transmission waveform is composed of multiple groups of LFSK signals with different frequency intervals , On the basis of unambiguous distance, the accuracy of distance measurement and speed measurement is significantly improved, and the measurement is more accurate; (4) When detecting targets, the upper and lower scanning frequency bands are used to work together, effectively reducing the detection blind zone.

附图说明Description of drawings

图1为MS-LFSK发射波形的整体示意图。FIG. 1 is an overall schematic diagram of an MS-LFSK transmission waveform.

图2(a)为MS-LFSK的上扫频段示意图,图2(b)为MS-LFSK下扫频段示意图。Fig. 2(a) is a schematic diagram of MS-LFSK up-sweep frequency band, and Fig. 2(b) is a schematic diagram of MS-LFSK down-sweep frequency band.

图3为本发明中相位差估计误差和它引起的测距误差的关系图,其中(a)为相位差估计误差的概率密度示意图,(b)为测距误差的概率密度示意图。Fig. 3 is the relationship diagram of the phase difference estimation error and the ranging error caused by it in the present invention, wherein (a) is a schematic diagram of the probability density of the phase difference estimation error, and (b) is a schematic diagram of the probability density of the ranging error.

图4是本发明以上扫频段为例的信号处理流程图。Fig. 4 is a flow chart of signal processing of the above frequency sweep band as an example in the present invention.

图5是本发明距离解模糊的示意图。Fig. 5 is a schematic diagram of distance deblurring in the present invention.

图6是本发明的盲区示意图,其中(a)为上下扫频段各自的探测盲区示意图,(b)为MS-LFSK上下扫频联合探测时的盲区示意图。Fig. 6 is a schematic diagram of blind areas of the present invention, wherein (a) is a schematic diagram of detection blind areas in up and down frequency scanning bands, and (b) is a schematic diagram of blind areas in MS-LFSK up and down frequency scanning joint detection.

图7是本发明MS-LFSK与LFMCW、MFSK的探测精度随输入信噪比变化的对比图,其中(a)、(b)、(c)、(d)分别为测距均值、测速均值、测距标准差、测速标准差的对比示意图。Fig. 7 is a comparison diagram of the detection accuracy of MS-LFSK and LFMCW and MFSK of the present invention as the input signal-to-noise ratio changes, wherein (a), (b), (c), and (d) are respectively the average value of distance measurement, the average value of speed measurement, Schematic diagram of the comparison between the distance measurement standard deviation and the speed measurement standard deviation.

具体实施方式Detailed ways

本发明提出一种适用于汽车防撞雷达的发射波形设计方法,以及该发生波形对应的目标探测算法,主要包括波形结构设计、波形参数设计和信号处理设计三个步骤,其中前两个步骤为波形设计方法步骤,第三个步骤是利用设计好的发射波形进行目标探测的过程。MS-LFSK发射波形的结构设计如图1所示。The present invention proposes a transmission waveform design method suitable for automobile anti-collision radar, and the target detection algorithm corresponding to the generated waveform, which mainly includes three steps: waveform structure design, waveform parameter design and signal processing design, wherein the first two steps are In the steps of the waveform design method, the third step is the process of using the designed launch waveform for target detection. The structural design of the MS-LFSK launch waveform is shown in Figure 1.

步骤一:设计MS-LFSK发射波形的结构:确定雷达发射信号的体制为MS-LFSK体制,一个周期的MS-LFSK波形由上扫频段和下扫频段构成,每个扫频段又由若干LFSK信号合成;Step 1: Design the structure of the MS-LFSK transmission waveform: determine that the system of the radar transmission signal is the MS-LFSK system, and a cycle of MS-LFSK waveform is composed of an up-sweep frequency band and a down-sweep frequency band, and each sweep frequency band is composed of several LFSK signals synthesis;

本发明一个周期的MS-LFSK波形由上扫频段和下扫频段构成,如图1所示。一个上扫频段和一个下扫频段共同组成一个相参处理周期,每个扫频段各自由M个LFSK信号交替排列而成。The MS-LFSK waveform of one cycle of the present invention is composed of an up-sweep frequency segment and a down-sweep frequency segment, as shown in FIG. 1 . An up-sweep frequency segment and a down-sweep frequency segment together form a coherent processing cycle, and each frequency scan segment is formed by alternately arranging M LFSK signals.

上扫频段波形如图2(a)所示,它由M个LFSK信号Sm,n构成:每个频率步进信号,如S1,N,有N阶步进,其中每个步进的持续时间为TStep,两个步进之间的时间间隔为(M-1)TStep,频率差为fIncr;就这M个频率步进信号整体而言,每个信号的起始时间都比上一个信号滞后TStep,起始频率也比上一个信号的减小即有:The waveform of the up-sweep frequency band is shown in Figure 2(a), which is composed of M LFSK signals S m,n : each frequency step signal, such as S 1,N , has N steps, and the The duration is T Step , the time interval between two steps is (M-1)T Step , and the frequency difference is f Incr ; with respect to the M frequency step signals as a whole, the start time of each signal is T Step lags behind the previous signal, and the starting frequency is also reduced compared to the previous signal That is:

为步进Sm,n对应的频率。 is the frequency corresponding to the step S m,n .

下扫频的波形是上扫频波形的镜向对称形式,波形如图2(b)所示。The waveform of the down-sweep frequency is a mirror-symmetrical form of the up-sweep frequency waveform, and the waveform is shown in Fig. 2(b).

步骤二:设计MS-LFSK发射信号的波形参数:根据给定的指标确定发射波形的参数,指标有:雷达发射频段、最大探测距离、最大探测速度、距离分辨力、测距精度、测速精度和频域最小信可检测噪比,确定一个波形需要的参数有:调频带宽B,相参处周期TCPI,LFSK信号的数量M,每个LFSK信号的步进数N以及LFSK信号之间的起始频差 Step 2: Design the waveform parameters of the MS-LFSK transmission signal: Determine the parameters of the transmission waveform according to the given indicators. The indicators include: radar transmission frequency band, maximum detection distance, maximum detection speed, distance resolution, ranging accuracy, speed measurement accuracy and The minimum detectable noise ratio in the frequency domain, the parameters required to determine a waveform are: FM bandwidth B, phase coherent period T CPI , the number of LFSK signals M, the number of steps of each LFSK signal N and the starting point between LFSK signals initial frequency difference

以发射信号的上扫频为例,调频带宽由距离分辨力ΔR决定,即Taking the up-sweep frequency of the transmitted signal as an example, the FM bandwidth is determined by the distance resolution ΔR, namely

其中c为光速;where c is the speed of light;

确定相参处理周期TCPI时,为了防止在做FFT时出现距离门的走动,要求在一个扫频周期内,车辆的最大位移不得超过距离分辨力ΔR,即When determining the coherent processing period T CPI , in order to prevent the range gate from moving when doing FFT, it is required that the maximum displacement of the vehicle must not exceed the distance resolution ΔR within a frequency sweep period, that is,

Vmax为最大可探测速度;V max is the maximum detectable speed;

求得TCPI的值以后,根据公式还可以得到扫频斜率;After obtaining the value of T CPI , according to the formula The frequency sweep slope can also be obtained;

根据采样定理,对每个频率键控信号的采样率应不小于差拍信号频率的最大值与最小值之差,即According to the sampling theorem, the sampling rate of each frequency keyed signal It should not be less than the difference between the maximum value and the minimum value of the beat signal frequency, that is

其中Rmax是最大探测距离,|Vmax|是最大探测速度的绝对值,fd.max是多普勒频率的最大值,τmax是回波时延的最大值,λ是发射信号的波长。且为了能使用FFT来提高运算速度,在上式的基础上,每个LFSK信号的步进数N还应满足:Where R max is the maximum detection distance, |V max | is the absolute value of the maximum detection velocity, f d.max is the maximum value of the Doppler frequency, τ max is the maximum value of the echo delay, λ is the wavelength of the transmitted signal . And in order to use FFT to improve the operation speed, on the basis of the above formula, the step number N of each LFSK signal should also satisfy:

log2(N)∈Z (5)log 2 (N)∈Z (5)

Z为整数;Z is an integer;

于是,根据N就可以得到一个LFSK信号中相邻步进的频率差fIncr和每个步进持续的时间TStepTherefore, according to N, the frequency difference f Incr of adjacent steps in an LFSK signal and the duration of each step T Step can be obtained:

在MS-LFSK体制中,目标在不同的LFSK信号处对应的差拍信号频率相同、相位不同,利用目标谱线处的频率和相位差可以解距离和速度的模糊。对相邻两个LFSK信号对应的峰值谱线处的相位做差,得到相位差 In the MS-LFSK system, the frequency and phase of the beat signals corresponding to different LFSK signals of the target are the same, and the ambiguity of distance and velocity can be resolved by using the frequency and phase difference at the target spectral line. The phase at the peak spectral lines corresponding to two adjacent LFSK signals and do difference, get phase difference

与目标的距离R0满足: The distance R 0 from the target satisfies:

fb为目标谱线处的差拍信号频率;f b is the beat signal frequency at the target spectral line;

的估计误差为越小测距精度越高,它服从期望为0、方差为的正态分布,其中满足:remember The estimated error of The smaller the distance is, the higher the accuracy is, and it follows that the expectation is 0 and the variance is normal distribution of , where Satisfy:

SNR为时域信噪比,(SNR·N)就构成差拍信号匹配滤波后的频域信噪比,δ为差拍频率的真实值与测量值的误差,δ满足|δ|≤0.5;SNR is the signal-to-noise ratio in the time domain, (SNR·N) constitutes the signal-to-noise ratio in the frequency domain after the beat signal is matched and filtered, δ is the error between the real value and the measured value of the beat frequency, and δ satisfies |δ|≤0.5;

如图3(a)所示,这是目标在频域最小可检测信噪比处、δ=0.5时的相位差估计误差的概率密度图,由于相位差估计误差和测距误差服从的分布一样,所以测距误差的概率密度与相位差的类似,如图3(b)所示。As shown in Figure 3(a), this is the probability density map of the phase difference estimation error of the target at the minimum detectable signal-to-noise ratio in the frequency domain and when δ=0.5, because the distribution of the phase difference estimation error and the ranging error is the same , so the probability density of the ranging error is similar to that of the phase difference, as shown in Figure 3(b).

在最小可检测信噪比的条件下,选择合适的置信区间(1-α)进行正态分布的区间估计,得到:Under the condition of the minimum detectable signal-to-noise ratio, choose an appropriate confidence interval (1-α) for interval estimation of normal distribution, and get:

由此可以估计出在最小可检测信噪比条件下的相位差估计误差的最大值 From this, the maximum value of the phase difference estimation error under the condition of the minimum detectable signal-to-noise ratio can be estimated

根据和测距精度στ计算不模糊距离的最小值 according to and ranging accuracy σ τ to calculate the minimum value of the unambiguous distance

当发射信号的不模糊距离不小于时,它对应求解出来的距离就能够满足测距精度的要求。When the unambiguous distance of the transmitted signal is not less than When , it corresponds to the calculated distance and can meet the requirements of ranging accuracy.

不模糊距离的上限理论上可以选取最大探测距离,考虑到噪声的影响,将其放宽至:The upper limit of the unambiguous distance can theoretically select the maximum detection distance. Considering the influence of noise, it can be relaxed to:

本发明中利用不模糊距离递推的方法提高探测精度,需要用二等分法将不模糊距离从开始等分作为In the present invention, the method of unambiguous distance recursion is used to improve detection accuracy, and the unambiguous distance needs to be bisected from start equally as which is

直至不模糊距离小于 up to unambiguous distance less than

与各LFSK信号之间初始频率的差值满足如下关系: The difference from the initial frequency of each LFSK signal Satisfy the following relationship:

越小,越大,就越不容易产生测距的模糊。 smaller, The larger it is, the less likely it is to produce ranging blur.

二等分的次数决定了LFSK信号的个数M,有:The number of bisections determines the number M of LFSK signals, which are:

由此可以得到各个不模糊距离再将的数值代入式(14),就可以求得 From this we can get each unambiguous distance then Substituting the value of into formula (14), we can get

至此,根据上述各参数可以确定MS-LFSK信号的上扫频波形,而下扫频波形为上扫频的镜像,所以也可以推算出下扫频的发射波形参数。So far, the up-sweep waveform of the MS-LFSK signal can be determined according to the above parameters, and the down-sweep waveform is the mirror image of the up-sweep, so the transmit waveform parameters of the down-sweep can also be deduced.

步骤三:设计MS-LFSK的信号处理方式:对发射信号的回波数据进行处理,求解目标的距离和速度:将接收到的回波信号与本振混频,得到差拍信号,对其采样并按照回波数据对应LFSK信号的情况对它重排,对重排后的数据按扫频周期做FFT和周期图积累,并对得到的二维频谱进行OS-CFAR,得到的频谱峰值点就对应一个目标,根据峰值坐标可以计算出目标对应的差拍信号频率,再由相频响应计算得到峰值点在相邻LFSK信号处的相位差,综合差拍信号频率和相位差,再经过解模糊的处理,就可以求解出目标的距离和速度。Step 3: Design the signal processing method of MS-LFSK: process the echo data of the transmitted signal, and solve the distance and speed of the target: mix the received echo signal with the local oscillator to obtain the beat signal, and sample it And rearrange it according to the echo data corresponding to the LFSK signal, do FFT and periodogram accumulation for the rearranged data according to the frequency sweep cycle, and perform OS-CFAR on the obtained two-dimensional spectrum, and the peak point of the obtained spectrum is Corresponding to a target, the beat signal frequency corresponding to the target can be calculated according to the peak coordinates, and then the phase difference of the peak point at the adjacent LFSK signal is calculated from the phase-frequency response, and the beat signal frequency and phase difference are integrated, and then defuzzified After processing, the distance and speed of the target can be calculated.

在MS-LFSK的上扫频段中,信号处理的过程如图4所示。具体过程如下:In the up-sweep frequency band of MS-LFSK, the process of signal processing is shown in Fig. 4 . The specific process is as follows:

采样差拍信号。将接收到的回波信号与本振混频,得到差拍信号,并在每一个TStep的结尾对差拍信号采样,得到采样结果的数组 Sample the beat signal. Mix the received echo signal with the local oscillator to obtain the beat signal, and sample the beat signal at the end of each T Step to obtain an array of sampling results

数据重排。将数组按M的值重排成M×N的二维数组,即二维数组的每一行对应的是一个完整的频率步进信号按发射顺序排列的差拍结果。Data rearrangement. the array According to the value of M, it is rearranged into a two-dimensional array of M×N, that is, each row of the two-dimensional array corresponds to a beat result of a complete frequency step signal arranged in the order of transmission.

FFT与周期图积累。对二维数组的每一行单独做FFT,得到M个频谱Wk[N],k=1,2,...,M。将Wk[N]取模,做模平方之后把M行结果对应相加,完成周期图的积累,即FFT with periodogram accumulation. Perform FFT on each row of the two-dimensional array separately to obtain M frequency spectra W k [N], k=1,2,...,M. Take the modulus of W k [N], do the modulus square and then add up the results of the M rows correspondingly to complete the accumulation of the periodogram, that is

进行谱峰的恒虚警检测。由于实际环境中目标可能比较多,所以采用OS-CFAR。Perform constant false alarm detection for spectral peaks. Since there may be many targets in the actual environment, OS-CFAR is adopted.

求解谱峰对应的频率和相位。假设超过门限的谱峰所在的检测单元为p,那么对应的差拍频率fb(p)可由下式确定:Find the frequency and phase corresponding to the spectral peak. Assuming that the detection unit where the spectral peak exceeding the threshold is located is p, then the corresponding beat frequency f b (p) can be determined by the following formula:

该谱峰对应的相位为:The phase corresponding to the spectral peak for:

那么相邻两个谱峰的相位差为:Then the phase difference between two adjacent spectral peaks is:

求解模糊距离根据已知的参数fb(p)和由下式求解模糊距离 Solve for blur distance According to known parameters f b (p) and Solve the fuzzy distance by the following formula

k=1,2,...,M-1k=1,2,...,M-1

上式中,的值随着的减小而逐渐精确,但也越来越模糊,所以还需要对距离解模糊,示意图如图5。用R0表示最终检测到的目标的距离,首先把的值赋给R0,然后计算k从2到M-1时的寻找满足下式的qk(p):In the above formula, The value of The reduction of the distance gradually becomes more accurate, but it also becomes more and more blurred, so it is necessary to deblur the distance, as shown in Figure 5. Use R 0 to represent the distance of the final detected target, first put Assign the value to R 0 , and then calculate when k is from 2 to M-1 Find q k (p) satisfying:

并将新的结果迭代给给R0,即and iterate the new result to R 0 , ie

一直迭代到k=M-1,得到最终测量的目标距离R0Iterate until k=M-1 to obtain the final measured target distance R 0 .

求解目标的速度。将R0代入下面公式得到目标的速度V0The speed at which the target is solved. Substitute R 0 into the following formula to get the speed V 0 of the target:

至此,在上扫频段中完成了对目标距离和速度的计算,下扫频段的数据处理方式与上扫频一致。So far, the calculation of the target distance and speed has been completed in the up-sweep frequency band, and the data processing method in the down-sweep frequency band is consistent with the up-sweep frequency band.

另外在实际系统中,由于连续波雷达的天线是收发分置并且一直工作的,所以接收天线在接收目标回波信号的同时也会接收到发射信号,那么经混频后就会产生干扰的零频信号,这可能会淹没零频附近的目标回波信号,形成探测盲区,如图6(a)。由示意图可以看出,上下扫频段的探测盲区绝大部分并不重合,因此MS-LFSK中利用上下扫频段联合探测目标,可以在上扫频段检测到原本盲区中远离的目标,在下扫频段检测到原本盲区中驶近的目标,以减小探测盲区,如图6(b)所示。In addition, in the actual system, since the antenna of the continuous wave radar is separated from the transceiver and works all the time, the receiving antenna will also receive the transmitted signal while receiving the target echo signal, then the interference zero will be generated after mixing. frequency signal, which may overwhelm the target echo signal near zero frequency, forming a detection blind zone, as shown in Figure 6(a). It can be seen from the schematic diagram that most of the detection blind areas in the upper and lower scanning frequency bands do not overlap. Therefore, in MS-LFSK, using the upper and lower scanning frequency bands to jointly detect targets can detect targets that are far away from the original blind area in the upper scanning frequency band, and detect them in the lower scanning frequency band. Go to the approaching target in the original blind area to reduce the detection blind area, as shown in Figure 6(b).

下面结合实施例对本发明进行详细说明。The present invention will be described in detail below in conjunction with examples.

实施例Example

本实施例中,设定的雷达指标如表1所示:In the present embodiment, the radar index of setting is as shown in table 1:

表1Table 1

根据上述指标设计的MS-LFSK波形的参数如表2所示。The parameters of the MS-LFSK waveform designed according to the above indicators are shown in Table 2.

表2Table 2

此外,为了对比各种发射波形的性能,我们在相同的指标下设计了LFMCW和MFSK的波形参数。In addition, in order to compare the performance of various transmission waveforms, we designed the waveform parameters of LFMCW and MFSK under the same index.

以变调频斜率LFMCW为例,该体制下的发射波形参数如表3所示:Taking LFMCW as an example, the transmit waveform parameters under this system are shown in Table 3:

表3table 3

MFSK体制下的发射波形参数如表4所示:The transmit waveform parameters under the MFSK system are shown in Table 4:

表4Table 4

图7为matlab的仿真结果,设置一个距离为30m,速度为10m/s的单一目标,将频域信噪比在15dB到40dB之间取值,每隔1dB对3种波形的测距测速做1000次蒙特卡洛实验,并研究它们求得的距离和速度的均值与标准差。Figure 7 shows the simulation results of matlab. Set a single target with a distance of 30m and a speed of 10m/s, set the frequency domain signal-to-noise ratio between 15dB and 40dB, and measure the distance and speed of the three waveforms every 1dB. Take 1000 Monte Carlo experiments and study the mean and standard deviation of the distances and velocities obtained from them.

图中,*表示变调频斜率LFMCW的探测效果,+表示MFSK信号的探测结果,o刻画了MS-LFSK体制下的探测结果。In the figure, * indicates the detection effect of the frequency modulation slope LFMCW, + indicates the detection result of the MFSK signal, and o depicts the detection result under the MS-LFSK system.

由图7(a)、图7(b)可以看出,信噪比较高时,3种波形均能准确的测距测速;而当频域信噪比低于20dB时,变调频斜率LFMCW的测量均值发生较大偏差,MFSK信号已不能正确测量,相反,MS-LFSK波形的测量结果还是很准确。It can be seen from Fig. 7(a) and Fig. 7(b) that when the SNR is high, all three waveforms can accurately measure distance and speed; and when the SNR in the frequency domain is lower than 20dB, the FM slope LFMCW The measurement mean value of MFSK has a large deviation, and the MFSK signal can no longer be measured correctly. On the contrary, the measurement result of the MS-LFSK waveform is still very accurate.

由图7(c)、图7(d)可以看出,信噪比较高时,3种波形测距测速的精度都比较高;而当频域信噪比低于19dB时,变调频斜率LFMCW的测距测速精度出现显著的下降;当频域信噪比低于24dB时,MFSK信号的测量精度就已不达标。而在频域信噪比的变化过程中,MS-LFSK波形的测量精度一直很高。It can be seen from Fig. 7(c) and Fig. 7(d) that when the SNR is high, the accuracy of the distance measurement and speed measurement of the three waveforms is relatively high; and when the SNR in the frequency domain is lower than 19dB, the FM slope The accuracy of distance measurement and speed measurement of LFMCW has dropped significantly; when the signal-to-noise ratio in the frequency domain is lower than 24dB, the measurement accuracy of MFSK signals has not reached the standard. In the process of changing the SNR in the frequency domain, the measurement accuracy of the MS-LFSK waveform is always high.

综合上述,从探测的准确度和精度两个方面来衡量,当频域信噪比下降时,变调频斜率LFMCW的探测性能显著下降,MFSK体制在性能变差方面表现的更加明显,相反,在频域信噪比变化的过程中,MS-LFSK波形的探测效果一直很好。To sum up the above, measured from the two aspects of detection accuracy and precision, when the signal-to-noise ratio in the frequency domain decreases, the detection performance of the FM slope LFMCW decreases significantly, and the performance of the MFSK system is more obvious in terms of performance degradation. On the contrary, in In the process of frequency domain SNR change, the detection effect of MS-LFSK waveform has been very good.

Claims (7)

1.一种适用于汽车防撞雷达的发射波形设计方法,其特征在于,包括以下步骤:1. a method for designing a launch waveform suitable for automobile anti-collision radar, is characterized in that, comprises the following steps: 步骤1,确定雷达发射信号的体制为MS-LFSK体制,一个周期的MS-LFSK波形由上扫频段和下扫频段构成,每个扫频段由若干LFSK信号合成;Step 1, determine that the system of the radar transmitting signal is the MS-LFSK system, a cycle of MS-LFSK waveform is composed of an up-sweep frequency segment and a down-sweep frequency segment, and each frequency scan segment is synthesized by several LFSK signals; 步骤2,根据给定的指标确定发射波形的参数,指标有:雷达发射频段、最大探测距离、最大探测速度、距离分辨力、测距精度、测速精度和频域最小信可检测噪比,确定一个周期的波形需要的参数有:调频带宽、相参处理间隔、LFSK信号的数量、每个LFSK信号的步进数以及LFSK信号之间的起始频差。Step 2. Determine the parameters of the transmission waveform according to the given indicators. The indicators include: radar transmission frequency band, maximum detection distance, maximum detection speed, distance resolution, ranging accuracy, speed measurement accuracy and minimum signal detectable noise ratio in the frequency domain. Determine The parameters required for a waveform of one cycle include: FM bandwidth, coherent processing interval, number of LFSK signals, number of steps for each LFSK signal, and initial frequency difference between LFSK signals. 2.根据权利要求1所述适用于汽车防撞雷达的发射波形设计方法,其特征在于,步骤1具体为:2. according to claim 1, be applicable to the transmission waveform design method of automobile anti-collision radar, it is characterized in that, step 1 is specially: 一个周期的MS-LFSK波形由上扫频段和下扫频段构成,一个上扫频段和一个下扫频段共同组成一个相参处理周期TCPI,每个扫频段各自由M个LFSK信号交替排列而成;A cycle of MS-LFSK waveform consists of an up-sweep frequency band and a down-sweep frequency band, and an up-sweep frequency band and a down-sweep frequency band together form a coherent processing cycle T CPI , and each sweep frequency band is formed by alternate arrangement of M LFSK signals ; 上扫频段由M个LFSK信号Sm,n构成:每个频率步进信号,有N阶步进,其中每个步进的持续时间为TStep,两个步进之间的时间间隔为(M-1)TStep,频率差为fIncr;每个信号的起始时间都比上一个信号滞后TStep,起始频率也比上一个信号的减小即有:The up-sweep frequency band is composed of M LFSK signals S m,n : each frequency step signal has N steps, where the duration of each step is T Step , and the time interval between two steps is ( M-1) T Step , the frequency difference is f Incr ; the start time of each signal is T Step behind the previous signal, and the start frequency is also smaller than that of the previous signal That is: 为步进Sm,n对应的频率; is the frequency corresponding to the step S m,n ; 下扫频的波形是上扫频波形的镜向对称形式。The waveform of the frequency sweep down is the mirror symmetrical form of the frequency sweep up waveform. 3.根据权利要求1所述适用于汽车防撞雷达的发射波形和对应的目标探测算法,其特征在于,步骤2具体方法为:3. according to claim 1, be applicable to the launch waveform of automobile anti-collision radar and corresponding target detection algorithm, it is characterized in that, the specific method of step 2 is: 首先需要知道的发射波形参数如下:调频带宽B,相参处理周期TCPI,LFSK信号的数量M,每个LFSK信号的步进数N和各LFSK信号之间的起始频差 First, the transmit waveform parameters that need to be known are as follows: FM bandwidth B, coherent processing period T CPI , number of LFSK signals M, step number N of each LFSK signal and initial frequency difference between each LFSK signal 以发射信号的上扫频为例,调频带宽由距离分辨力ΔR决定,即Taking the up-sweep frequency of the transmitted signal as an example, the FM bandwidth is determined by the distance resolution ΔR, namely 其中c为光速;where c is the speed of light; 确定相参处理周期TCPI时,在一个扫频周期内,车辆的最大位移不得超过距离分辨力ΔR,即When determining the coherent processing period T CPI , the maximum displacement of the vehicle must not exceed the distance resolution ΔR within a frequency sweep period, that is, Vmax为最大可探测速度;V max is the maximum detectable speed; 求得TCPI的值以后,根据公式可得到扫频斜率;After obtaining the value of T CPI , according to the formula The frequency sweep slope can be obtained; 根据采样定理,对每个频率键控信号的采样率应不小于差拍信号频率的最大值与最小值之差,即According to the sampling theorem, the sampling rate of each frequency keyed signal It should not be less than the difference between the maximum value and the minimum value of the beat signal frequency, that is 其中Rmax是最大可探测距离,|Vmax|是最大可探测速度的绝对值,fd.max是多普勒频率的最大值,τmax是回波时延的最大值,λ是发射信号的波长;在上式的基础上,每个LFSK信号的步进数N还应满足:where R max is the maximum detectable distance, |V max | is the absolute value of the maximum detectable velocity, f d.max is the maximum value of Doppler frequency, τ max is the maximum value of echo delay, λ is the transmitted signal wavelength; on the basis of the above formula, the step number N of each LFSK signal should also satisfy: log2(N)∈Z (5)log 2 (N)∈Z (5) Z为整数;Z is an integer; 于是,根据N可得到一个LFSK信号中相邻步进的频率差fIncr和每个步进持续的时间TStepTherefore, according to N, the frequency difference f Incr of adjacent steps in an LFSK signal and the duration of each step T Step can be obtained: 对相邻两个LFSK信号对应的峰值谱线处相位做差,得到相位差 The phase at the peak spectral line corresponding to two adjacent LFSK signals and do difference, get phase difference 与目标的距离R0满足: The distance R 0 from the target satisfies: fb为目标谱线处对应的差拍信号频率;f b is the beat signal frequency corresponding to the target spectral line; 的估计误差为它服从期望为0、方差为的正态分布,其中满足: The estimated error of It has an expectation of 0 and a variance of normal distribution of , where Satisfy: SNR为时域信噪比,(SNR·N)构成差拍信号匹配滤波后的频域信噪比,δ为差拍频率真实值与测量值的误差,δ满足|δ|≤0.5;SNR is the signal-to-noise ratio in the time domain, (SNR·N) constitutes the signal-to-noise ratio in the frequency domain after the beat signal is matched and filtered, δ is the error between the real value of the beat frequency and the measured value, and δ satisfies |δ|≤0.5; 在最小可检测信噪比的条件下,选择置信区间(1-α)进行正态分布的区间估计,得到:Under the condition of the minimum detectable signal-to-noise ratio, choose the confidence interval (1-α) for the interval estimation of the normal distribution, and get: 由此估计出在最小可检测信噪比条件下的相位差估计误差的最大值 From this, the maximum value of the phase difference estimation error under the condition of the minimum detectable signal-to-noise ratio is estimated 根据和测距精度στ计算不模糊距离的最小值 according to and ranging accuracy σ τ to calculate the minimum value of the unambiguous distance 求解不模糊距离的上限 Find an upper bound on the unambiguous distance 采用二等分法将不模糊距离从开始等分作为Use the bisection method to divide the unambiguous distance from start equally as which is 直至不模糊距离小于 up to unambiguous distance less than 与各LFSK信号之间初始频率的差值满足如下关系: The difference from the initial frequency of each LFSK signal Satisfy the following relationship: 二等分的次数决定了LFSK信号的个数M,有:The number of bisections determines the number M of LFSK signals, which are: 由此可得到各个不模糊距离再将的数值代入式(14),求得 From this, each unambiguous distance can be obtained then Substituting the value of into formula (14), we get 根据上述各参数确定MS-LFSK信号的上扫频波形,下扫频波形为上扫频的镜像。Determine the up-sweep waveform of the MS-LFSK signal according to the above parameters, and the down-sweep waveform is the mirror image of the up-sweep. 4.一种适用于汽车防撞雷达的目标探测算法,其特征在于,包括以下步骤:4. A target detection algorithm applicable to automobile collision avoidance radar, is characterized in that, comprises the following steps: 步骤1,确定雷达发射信号的体制为MS-LFSK体制,一个周期的MS-LFSK波形由上扫频段和下扫频段构成,每个扫频段由若干LFSK信号合成;Step 1, determine that the system of the radar transmitting signal is the MS-LFSK system, a cycle of MS-LFSK waveform is composed of an up-sweep frequency segment and a down-sweep frequency segment, and each frequency scan segment is synthesized by several LFSK signals; 步骤2,根据给定的指标确定发射波形的参数,指标有:雷达发射频段、最大探测距离、最大探测速度、距离分辨力、测距精度、测速精度和频域最小信可检测噪比,确定一个周期的波形需要的参数有:调频带宽、相参处理间隔、LFSK信号的数量、每个LFSK信号的步进数以及LFSK信号之间的起始频差;Step 2. Determine the parameters of the transmission waveform according to the given indicators. The indicators include: radar transmission frequency band, maximum detection distance, maximum detection speed, distance resolution, ranging accuracy, speed measurement accuracy and minimum signal detectable noise ratio in the frequency domain. Determine The parameters required for a cycle waveform are: FM bandwidth, coherent processing interval, number of LFSK signals, number of steps for each LFSK signal, and initial frequency difference between LFSK signals; 步骤3,对发射信号的回波数据进行处理,求解目标的距离和速度:将接收到的回波信号与本振混频,得到差拍信号,对其采样并按照回波数据对应LFSK信号的情况对它重排,对重排后的数据按扫频周期做FFT和周期图积累,并对得到的二维频谱进行OS-CFAR,得到的频谱峰值点对应一个目标,根据峰值坐标计算出目标对应的差拍信号频率,再由相频响应计算得到峰值点在相邻LFSK信号处的相位差,综合差拍信号频率和相位差,再经过解模糊的处理,求解出目标的距离和速度。Step 3, process the echo data of the transmitted signal to solve the distance and speed of the target: mix the received echo signal with the local oscillator to obtain the beat signal, sample it and correspond to the LFSK signal according to the echo data Rearrange it according to the situation, do FFT and periodogram accumulation on the rearranged data according to the sweep cycle, and perform OS-CFAR on the obtained two-dimensional spectrum, the obtained spectrum peak point corresponds to a target, and calculate the target according to the peak coordinates The corresponding beating signal frequency is calculated from the phase-frequency response to obtain the phase difference of the peak point at the adjacent LFSK signal, and the frequency and phase difference of the beating signal are integrated, and then the distance and speed of the target are obtained through defuzzification processing. 5.根据权利要求4所述的适用于汽车防撞雷达的目标探测算法,其特征在于,步骤1具体为:5. the target detection algorithm applicable to automobile collision avoidance radar according to claim 4, is characterized in that, step 1 is specially: 一个周期的MS-LFSK波形由上扫频段和下扫频段构成,一个上扫频段和一个下扫频段共同组成一个相参处理周期TCPI,每个扫频段各自由M个LFSK信号交替排列而成;A cycle of MS-LFSK waveform consists of an up-sweep frequency band and a down-sweep frequency band, and an up-sweep frequency band and a down-sweep frequency band together form a coherent processing cycle T CPI , and each sweep frequency band is formed by alternate arrangement of M LFSK signals ; 上扫频段由M个LFSK信号Sm,n构成:每个频率步进信号,有N阶步进,其中每个步进的持续时间为TStep,两个步进之间的时间间隔为(M-1)TStep,频率差为fIncr;就这M个频率步进信号整体而言,每个信号的起始时间都比上一个信号滞后TStep,起始频率也比上一个信号的减小即有:The up-sweep frequency band is composed of M LFSK signals S m,n : each frequency step signal has N steps, where the duration of each step is T Step , and the time interval between two steps is ( M-1) T Step , the frequency difference is f Incr ; as far as the M frequency step signals are concerned as a whole, the start time of each signal is T Step behind the previous signal, and the start frequency is also higher than that of the previous signal. decrease That is: 为步进Sm,n对应的频率; is the frequency corresponding to the step S m,n ; 下扫频的波形是上扫频波形的镜向对称形式。The waveform of the frequency sweep down is the mirror symmetrical form of the frequency sweep up waveform. 6.根据权利要求4所述的适用于汽车防撞雷达的目标探测算法,其特征在于,步骤2具体为:6. the target detection algorithm applicable to automobile collision avoidance radar according to claim 4, characterized in that, step 2 is specifically: 首先需要知道的发射波形参数如下:调频带宽B,相参处理周期TCPI,LFSK信号的数量M,每个LFSK信号的步进数N和各LFSK信号之间的起始频差 First, the transmit waveform parameters that need to be known are as follows: FM bandwidth B, coherent processing period T CPI , number of LFSK signals M, step number N of each LFSK signal and initial frequency difference between each LFSK signal 以发射信号的上扫频为例,调频带宽由距离分辨力ΔR决定,即Taking the up-sweep frequency of the transmitted signal as an example, the FM bandwidth is determined by the distance resolution ΔR, namely 其中c为光速;where c is the speed of light; 确定相参处理周期TCPI时,在一个扫频周期内,车辆的最大位移不得超过距离分辨力ΔR,即When determining the coherent processing period T CPI , the maximum displacement of the vehicle must not exceed the distance resolution ΔR within a frequency sweep period, that is, Vmax为最大可探测速度;V max is the maximum detectable speed; 求得TCPI的值以后,根据公式可得到扫频斜率;After obtaining the value of T CPI , according to the formula The frequency sweep slope can be obtained; 根据采样定理,对每个频率键控信号的采样率应不小于差拍信号频率的最大值与最小值之差,即According to the sampling theorem, the sampling rate of each frequency keyed signal It should not be less than the difference between the maximum value and the minimum value of the beat signal frequency, that is 其中Rmax是最大可探测距离,|Vmax|是最大可探测速度的绝对值,fd.max是多普勒频率的最大值,τmax是回波时延的最大值,λ是发射信号的波长;在上式的基础上,每个LFSK信号的步进数N还应满足:where R max is the maximum detectable distance, |V max | is the absolute value of the maximum detectable velocity, f d.max is the maximum value of Doppler frequency, τ max is the maximum value of echo delay, λ is the transmitted signal wavelength; on the basis of the above formula, the step number N of each LFSK signal should also satisfy: log2(N)∈Z (5)log 2 (N)∈Z (5) Z为整数;Z is an integer; 于是,根据N可得到一个LFSK信号中相邻步进的频率差fIncr和每个步进持续的时间TStepTherefore, according to N, the frequency difference f Incr of adjacent steps in an LFSK signal and the duration of each step T Step can be obtained: 对相邻两个LFSK信号对应的峰值谱线处相位做差,得到相位差 The phase at the peak spectral line corresponding to two adjacent LFSK signals and do difference, get phase difference 与目标的距离R0满足: The distance R 0 from the target satisfies: fb为目标谱线处对应的差拍信号频率;f b is the beat signal frequency corresponding to the target spectral line; 的估计误差为它服从期望为0、方差为的正态分布,其中满足: The estimated error of It has an expectation of 0 and a variance of normal distribution of , where Satisfy: SNR为时域信噪比,(SNR·N)构成差拍信号匹配滤波后的频域信噪比,δ为差拍频率真实值与测量值的误差,δ满足|δ|≤0.5;SNR is the signal-to-noise ratio in the time domain, (SNR·N) constitutes the signal-to-noise ratio in the frequency domain after the beat signal is matched and filtered, δ is the error between the real value of the beat frequency and the measured value, and δ satisfies |δ|≤0.5; 在最小可检测信噪比的条件下,选择置信区间(1-α)进行正态分布的区间估计,得到:Under the condition of the minimum detectable signal-to-noise ratio, choose the confidence interval (1-α) for the interval estimation of the normal distribution, and get: 由此估计出在最小可检测信噪比条件下的相位差估计误差的最大值 From this, the maximum value of the phase difference estimation error under the condition of the minimum detectable signal-to-noise ratio is estimated 根据和测距精度στ计算不模糊距离的最小值 according to and ranging accuracy σ τ to calculate the minimum value of the unambiguous distance 求解不模糊距离的上限 Find an upper bound on the unambiguous distance 采用二等分法将不模糊距离从开始等分作为Use the bisection method to divide the unambiguous distance from start equally as which is 直至不模糊距离小于 up to unambiguous distance less than 与各LFSK信号之间初始频率的差值满足如下关系: The difference from the initial frequency of each LFSK signal Satisfy the following relationship: 二等分的次数决定了LFSK信号的个数M,有:The number of bisections determines the number M of LFSK signals, which are: 由此可得到各个不模糊距离再将的数值代入式(14),求得 From this, each unambiguous distance can be obtained then The value of is substituted into formula (14), to get 根据上述各参数确定MS-LFSK信号的上扫频波形,下扫频波形为上扫频的镜像。Determine the up-sweep waveform of the MS-LFSK signal according to the above parameters, and the down-sweep waveform is the mirror image of the up-sweep. 7.根据权利要求4所述的适用于汽车防撞雷达的目标探测算法,其特征在于,步骤3具体为:7. the target detection algorithm that is applicable to automobile anti-collision radar according to claim 4, is characterized in that, step 3 is specially: 采样差拍信号;将接收到的回波信号与本振信号混频,得到差拍信号,并在每一个TStep的结尾对差拍信号采样,得到采样结果的数组 Sampling the beat signal; mixing the received echo signal with the local oscillator signal to obtain the beat signal, and sampling the beat signal at the end of each T Step to obtain an array of sampling results 数据重排;将数组按M的值重排成M×N的二维数组,即二维数组的每一行对应的是一个完整的频率步进信号按发射顺序排列的差拍结果;data rearrangement; array According to the value of M, it is rearranged into a two-dimensional array of M×N, that is, each row of the two-dimensional array corresponds to the beat result of a complete frequency step signal arranged in the order of transmission; FFT与周期图积累;对二维数组的每一行单独做FFT,得到M个频谱Wk[N],k=1,2,...,M;将Wk[N]取模,做模平方之后把M行结果对应相加,完成周期图的积累,即FFT and periodogram accumulation; perform FFT on each row of the two-dimensional array separately to obtain M spectrum W k [N],k=1,2,...,M; take the modulus of W k [N] and do the modulus After squaring, the results of the M rows are added correspondingly to complete the accumulation of the periodogram, that is 进行谱峰的恒虚警检测;Carry out constant false alarm detection of spectral peaks; 求解谱峰对应的频率和相位;假设超过门限的谱峰所在的检测单元为p,那么对应的差拍频率fb(p)可由下式确定:Solve the frequency and phase corresponding to the spectral peak; assuming that the detection unit where the spectral peak exceeds the threshold is p, then the corresponding beat frequency f b (p) can be determined by the following formula: 该谱峰对应的相位为:The phase corresponding to this spectral peak for: 相邻两个谱峰的相位差为:The phase difference between two adjacent spectral peaks is: 求解模糊距离根据已知的参数fb(p)和由下式求解模糊距离 Solve for blur distance According to known parameters f b (p) and Solve the fuzzy distance by the following formula 求解不模糊距离;用R0表示最终检测到的目标的距离,首先把的值赋给R0,然后计算k从2到M-1时的寻找满足下式的qk(p):Solve the unambiguous distance; use R 0 to represent the distance of the final detected target, first put Assign the value of R 0 to R 0 , and then calculate when k is from 2 to M-1 Find q k (p) satisfying: 并将新的结果迭代给R0,即and iterate the new result to R 0 , ie 一直迭代到k=M-1,得到最终测量的目标距离R0Iterate until k=M-1 to obtain the final measured target distance R 0 . 求解目标的速度;将R0代入下面的公式得到目标的速度V0Find the speed of the target; Substitute R 0 into the following formula to get the speed V 0 of the target: 至此,在上扫频段中完成了对目标距离和速度的计算,下扫频段的数据处理方式与上扫频一致。So far, the calculation of the target distance and speed has been completed in the up-sweep frequency band, and the data processing method in the down-sweep frequency band is consistent with the up-sweep frequency band.
CN201910473989.0A 2019-06-02 2019-06-02 Transmitted waveform design method and target acquisition algorithm suitable for automobile collision avoidance radar Pending CN110208804A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910473989.0A CN110208804A (en) 2019-06-02 2019-06-02 Transmitted waveform design method and target acquisition algorithm suitable for automobile collision avoidance radar

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910473989.0A CN110208804A (en) 2019-06-02 2019-06-02 Transmitted waveform design method and target acquisition algorithm suitable for automobile collision avoidance radar

Publications (1)

Publication Number Publication Date
CN110208804A true CN110208804A (en) 2019-09-06

Family

ID=67790192

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910473989.0A Pending CN110208804A (en) 2019-06-02 2019-06-02 Transmitted waveform design method and target acquisition algorithm suitable for automobile collision avoidance radar

Country Status (1)

Country Link
CN (1) CN110208804A (en)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110837081A (en) * 2019-11-08 2020-02-25 南京理工大学 High-speed target detection method based on P/D band radar signal fusion processing
CN110850400A (en) * 2019-11-13 2020-02-28 南京理工大学 Multi-target detection method for LFMCW radar based on interferometer direction finding
CN112198365A (en) * 2020-09-29 2021-01-08 中国船舶重工集团公司第七二四研究所 Method for realizing real-time frequency spectrum monitoring equipment
CN112462336A (en) * 2020-10-10 2021-03-09 东华大学 Self-adaptive elimination method for FMCW foreign matter detection radar leakage signal
CN112782687A (en) * 2021-03-18 2021-05-11 江苏蛮酷科技有限公司 Radar waveform generation circuit, method, radar, and computer storage medium
CN113759327A (en) * 2021-09-29 2021-12-07 桂林长海发展有限责任公司 Interference method and system for linear frequency modulation continuous wave radar and electronic equipment
WO2022037101A1 (en) * 2020-08-21 2022-02-24 上海禾赛科技股份有限公司 Method for using frequency modulation continuous wave to perform detection, and radar and computer-readable storage medium
CN114217288A (en) * 2022-02-22 2022-03-22 湖南纳雷科技有限公司 Method and system for synchronizing high coherence between chips of echo signals of radar
WO2022247593A1 (en) * 2021-05-28 2022-12-01 华为技术有限公司 Distance measurement method and apparatus
WO2024197774A1 (en) * 2023-03-31 2024-10-03 华为技术有限公司 Signal processing method, and apparatus

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US2682035A (en) * 1950-10-26 1954-06-22 Collins Radio Co Linear frequency shift keying circuit
US20160291156A1 (en) * 2013-10-23 2016-10-06 Ladar Limited A laser detection and ranging device for detecting an object under a water surface
CN106338727A (en) * 2016-09-05 2017-01-18 南京理工大学 Target detection method of auxiliary vehicle driving radar

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US2682035A (en) * 1950-10-26 1954-06-22 Collins Radio Co Linear frequency shift keying circuit
US20160291156A1 (en) * 2013-10-23 2016-10-06 Ladar Limited A laser detection and ranging device for detecting an object under a water surface
CN106338727A (en) * 2016-09-05 2017-01-18 南京理工大学 Target detection method of auxiliary vehicle driving radar

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
何俊杰: "基于多阶线性频率键控的汽车防撞雷达设计与实现", 《中国优秀硕士学位论文全文数据库 (信息科技辑)》 *

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110837081A (en) * 2019-11-08 2020-02-25 南京理工大学 High-speed target detection method based on P/D band radar signal fusion processing
CN110850400A (en) * 2019-11-13 2020-02-28 南京理工大学 Multi-target detection method for LFMCW radar based on interferometer direction finding
WO2022037101A1 (en) * 2020-08-21 2022-02-24 上海禾赛科技股份有限公司 Method for using frequency modulation continuous wave to perform detection, and radar and computer-readable storage medium
CN112198365A (en) * 2020-09-29 2021-01-08 中国船舶重工集团公司第七二四研究所 Method for realizing real-time frequency spectrum monitoring equipment
CN112462336A (en) * 2020-10-10 2021-03-09 东华大学 Self-adaptive elimination method for FMCW foreign matter detection radar leakage signal
CN112462336B (en) * 2020-10-10 2022-11-29 东华大学 Self-adaptive elimination method for FMCW foreign matter detection radar leakage signal
CN112782687A (en) * 2021-03-18 2021-05-11 江苏蛮酷科技有限公司 Radar waveform generation circuit, method, radar, and computer storage medium
WO2022247593A1 (en) * 2021-05-28 2022-12-01 华为技术有限公司 Distance measurement method and apparatus
CN113759327A (en) * 2021-09-29 2021-12-07 桂林长海发展有限责任公司 Interference method and system for linear frequency modulation continuous wave radar and electronic equipment
CN113759327B (en) * 2021-09-29 2024-05-14 桂林长海发展有限责任公司 Interference method, system and electronic equipment for linear frequency modulation continuous wave radar
CN114217288A (en) * 2022-02-22 2022-03-22 湖南纳雷科技有限公司 Method and system for synchronizing high coherence between chips of echo signals of radar
WO2024197774A1 (en) * 2023-03-31 2024-10-03 华为技术有限公司 Signal processing method, and apparatus

Similar Documents

Publication Publication Date Title
CN110208804A (en) Transmitted waveform design method and target acquisition algorithm suitable for automobile collision avoidance radar
CN108008391B (en) A Multi-Target Recognition Method for Vehicle Radar Based on FMCW
CN106338727A (en) Target detection method of auxiliary vehicle driving radar
CN101470202B (en) Pulse Doppler radar system and its signal processing method
US6606052B1 (en) Method and apparatus for detecting multiple objects with frequency modulated continuous wave radar
US9244164B2 (en) Method for unambiguously determining a range and/or a relative speed of an object, driver assistance device and motor vehicle
KR101737590B1 (en) Method and apparatus for radar signal processing
CN114137522B (en) Millimeter-wave orthogonal waveform optimization method and vehicle-mounted radar system for autonomous driving
CN105093227B (en) A kind of traffic flux measurement device and vehicle operating information preparation method
CN115575922B (en) Moving target detection method and device based on vehicle-mounted FMCW millimeter wave radar
CN102798863A (en) Road central isolation belt detection method based on automobile anti-collision radar
JP7138970B2 (en) Method and apparatus for determining at least one parameter of an object
US20080129582A1 (en) Tracking target method of radar with frequency modulation continuous wave
CN108415010A (en) A radar multi-target detection method based on trapezoidal LFMCW modulation
JP2004163340A (en) Onboard radar system
CN102788980A (en) Automobile anticollision radar system based on frequency-modulated continuous wave
JP7480422B2 (en) Radar method and radar system having high range resolution at low signal processing load
CN107783121B (en) Unmanned automobile anti-collision radar system signal processing system and method based on combined waveform
CN110161472A (en) A kind of broadband vehicle-mounted millimeter wave radar solution velocity ambiguity method based on signal multiplexing
CN115657007A (en) River surface flow velocity measuring method based on millimeter wave radar
CN111929685A (en) Radar detection method and device based on linear frequency modulation continuous wave, radar device and motor vehicle
CN112639521A (en) Method and device for controlling radar frequency hopping, radar speed measuring method and radar
CN115128592A (en) Debris flow surface flow velocity monitoring method and system
Yang et al. Anti-interference waveform design for automotive radar
CN112859061A (en) Multi-target detection method based on frequency modulation continuous wave radar

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20190906