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 PDFInfo
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 23
- 238000013461 design Methods 0.000 title claims abstract description 15
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 12
- 238000001514 detection method Methods 0.000 claims abstract description 55
- 230000005540 biological transmission Effects 0.000 claims abstract description 30
- 238000005259 measurement Methods 0.000 claims abstract description 23
- 238000012545 processing Methods 0.000 claims abstract description 20
- 230000001427 coherent effect Effects 0.000 claims abstract description 15
- 230000008569 process Effects 0.000 claims abstract description 9
- 238000004364 calculation method Methods 0.000 claims abstract description 6
- 230000003595 spectral effect Effects 0.000 claims description 17
- 238000005070 sampling Methods 0.000 claims description 10
- 238000009825 accumulation Methods 0.000 claims description 8
- 238000001228 spectrum Methods 0.000 claims description 8
- 230000035559 beat frequency Effects 0.000 claims description 5
- 230000007423 decrease Effects 0.000 claims description 5
- 238000006073 displacement reaction Methods 0.000 claims description 3
- 238000003672 processing method Methods 0.000 claims description 3
- 230000004044 response Effects 0.000 claims description 3
- 230000008707 rearrangement Effects 0.000 claims description 2
- 238000010009 beating Methods 0.000 claims 2
- 238000010586 diagram Methods 0.000 description 13
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 230000007812 deficiency Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 206010039203 Road traffic accident Diseases 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
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/88—Radar or analogous systems specially adapted for specific applications
- G01S13/93—Radar or analogous systems specially adapted for specific applications for anti-collision purposes
- G01S13/931—Radar or analogous systems specially adapted for specific applications for anti-collision purposes of land vehicles
-
- 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/88—Radar or analogous systems specially adapted for specific applications
- G01S13/93—Radar or analogous systems specially adapted for specific applications for anti-collision purposes
- G01S13/931—Radar or analogous systems specially adapted for specific applications for anti-collision purposes of land vehicles
- G01S2013/9315—Monitoring 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
技术领域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和每个步进持续的时间TStep:Therefore, 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,得到最终测量的目标距离R0。Iterate until k=M-1 to obtain the final measured target distance R 0 .
求解目标的速度。将R0代入下面公式得到目标的速度V0:The 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)
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)
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)
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 |
-
2019
- 2019-06-02 CN CN201910473989.0A patent/CN110208804A/en active Pending
Patent Citations (3)
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)
Title |
---|
何俊杰: "基于多阶线性频率键控的汽车防撞雷达设计与实现", 《中国优秀硕士学位论文全文数据库 (信息科技辑)》 * |
Cited By (12)
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 |