CN108120975A - 基于梯形连续波的雷达测速测距方法 - Google Patents
基于梯形连续波的雷达测速测距方法 Download PDFInfo
- Publication number
- CN108120975A CN108120975A CN201711356001.XA CN201711356001A CN108120975A CN 108120975 A CN108120975 A CN 108120975A CN 201711356001 A CN201711356001 A CN 201711356001A CN 108120975 A CN108120975 A CN 108120975A
- Authority
- CN
- China
- Prior art keywords
- radar
- target
- speed
- continuous wave
- measured
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/06—Systems determining position data of a target
- G01S13/42—Simultaneous measurement of distance and other co-ordinates
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/50—Systems of measurement based on relative movement of target
- G01S13/58—Velocity or trajectory determination systems; Sense-of-movement determination systems
- G01S13/583—Velocity or trajectory determination systems; Sense-of-movement determination systems using transmission of continuous unmodulated waves, amplitude-, frequency-, or phase-modulated waves and based upon the Doppler effect resulting from movement of targets
- G01S13/584—Velocity or trajectory determination systems; Sense-of-movement determination systems using transmission of continuous unmodulated waves, amplitude-, frequency-, or phase-modulated waves and based upon the Doppler effect resulting from movement of targets adapted for simultaneous range and velocity measurements
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Medicines That Contain Protein Lipid Enzymes And Other Medicines (AREA)
Abstract
本发明公开了一种基于梯形连续波的雷达测速测距方法,用于测量多个目标相对雷达的速度和距离,其步骤包括:雷达发射梯形连续波信号,所述梯形连续波包括正线性调频段波形、点频段波形和负线性调频段波形;雷达回波解调后,确定所有待测目标的实际速度值集合;确定每个待测目标的可能速度值;以及将所述每个待测目标的可能速度值与所述所有待测目标的实际速度值集合进行匹配得到待测目标的实际速度并且进一步得到待测目标的实际距离。所述基于梯形连续波的雷达测速测距方法可以有效实现雷达对密集多目标的测距测速,从而极大提高雷达对密集多目标的监视性能。
Description
技术领域
本发明涉及雷达检测领域,特别涉及一种基于梯形连续波的雷达测速测距方法。
背景技术
连续波雷达具有成本低、重量轻等突出优点,在近作用距离的民用领域具有重要的应用价值。目前连续波雷达发射波形主要包括点频信号、锯齿波信号及三角波信号等。其中点频信号主要用于测量目标的速度,锯齿波信号用于测量目标的距离,而三角波信号可同时进行目标测速和测距。
三角波信号需要对同一目标在正负调频斜率段信号中的不同位置进行正确配对,才能获得目标的正确速度和距离。当雷达波束内目标数目较少时,可以通过比较目标回波的功率信息进行配对,其基本思想是不同目标之间由于距离或雷达截面积的不同会导致回波功率强度不同,从而将功率接近的进行配对处理。但是当目标较为密集时,例如监控路面车辆的交通雷达或者行人的安防类雷达,其目标距离相差不大,并且雷达截面积也十分类似,此时利用回波功率难以区分,因此会导致配对错误。同时也有相关文献提出通过接收回波匹配滤波后的波形形状进行配对,但是实际上当目标类型接近时(例如车辆之间,行人之间),其雷达回波形状也基本相同,当目标密集时也难以适用。
公开于该背景技术部分的信息仅仅旨在增加对本发明的总体背景的理解,而不应当被视为承认或以任何形式暗示该信息构成已为本领域一般技术人员所公知的现有技术。
发明内容
本发明的目的在于提供一种基于梯形连续波的雷达测速测距方法,从而克服雷达对密集多目标无法准确测速测距的缺点。
为实现上述目的,本发明提供了一种基于梯形连续波的雷达测速测距方法,用于测量多个目标相对雷达的速度和距离。其步骤包括:雷达发射梯形连续波信号,所述梯形连续波包括正线性调频段波形、点频段波形和负线性调频段波形;雷达回波解调后,确定所有待测目标的实际速度值集合;确定每个待测目标的可能速度值;以及将所述每个待测目标的可能速度值与所述所有待测目标的实际速度值集合进行匹配得到待测目标的实际速度并且进一步得到待测目标的实际距离。
优选地,上述技术方案中,梯形连续波波形ST(t)公式如下:
其中k表示雷达发射梯形波的周期;f0为雷达发射信号中心频率;c为光速,波长γ为线性调频段信号的调频斜率;T1为正线性调频段或负线性调频段的持续时间;线性调频段信号的带宽B满足B=γT1;T2为点频段的持续时间。
优选地,上述技术方案中,所述雷达回波的波形分为正线性调频段,点频段波形,负线性调频段,分别近似表示为SR1(t),SR2(t),SR3(t),波形公式如下:
Ri表示第i个目标相对雷达的初始距离;Vi表示第i个目标相对雷达的初始速度;N为目标数目;ρi为与雷达发射功率、天线增益、目标特性及距离相关的复常数;N1(t),N2(t),N3(t)表示接收机噪声。
优选地,上述技术方案中,对所述SR1(t),SR2(t),SR3(t)分别进行FFT变换:
其中G1表示线性调频段的相干累积增益,G2表示线性调频段的相干累积增益。
优选地,上述技术方案中,对所述SR1(f1),SR2(f2),SR3(f3)进行单元平均恒虚警检测,检测到的目标峰值点在正线性调频段信号、点频段信号及负线性调频段信号的位置分别可表示为:
所述所有待测目标的实际速度值集合的集合为
优选地,上述技术方案中,将所述正线性调频段信号第m个检测的目标与所述负线性调频段信号的第n个目标进行配对,确定所述待测目标的可能速度值为
优选地,上述技术方案中,将Vmn与Vi2集合中所有数值逐一比较找到最接近数值Vp;如果其中K为常数,取值范围为[0,0.5],表示配对错误,令n=n+1,继续所述逐一比较过程;表示配对正确,输出该目标的正确速度Vmn和距离Rmn,所述如果m<N,令m=m+1,n=1,继续所述逐一比较过程;m≥N,密集目标测速测距结束。
与现有技术相比,本发明具有如下有益效果:
可以有效实现对密集多目标的雷达测距测速,从而极大提高连续波雷达系统对密集多目标的监视性能。
附图说明
图1是根据本发明一实施方式的雷达发射的梯形连续波的时频图。
图2是根据本发明一实施方式的雷达测距测速仿真结果。
具体实施方式
下面结合附图,对本发明的具体实施方式进行详细描述,但应当理解本发明的保护范围并不受具体实施方式的限制。
除非另有其它明确表示,否则在整个说明书和权利要求书中,术语“包括”或其变换如“包含”或“包括有”等等将被理解为包括所陈述的元件或组成部分,而并未排除其它元件或其它组成部分。
基于梯形连续波的雷达测速测距方法的原理是通过发射梯形连续波,利用梯形连续波中的点频段信号估计所有目标的速度信息,以此作为目标配对是否正确的检验准则,然后利用正负线调频段信号的回波进行配对并估计目标速度距离信息,如果配对正确则获得的速度信息必定与点频信号段获得的某个目标速度重合,否则就进行重新配对,通过该方法可以有效实现密集目标的正确配对,从而极大提高连续波雷达系统密集目标监视性能。
基于梯形连续波的雷达测速测距方法具体步骤如下:
1,雷达发射梯形连续波信号对密集目标检测,雷达发射梯形连续波包括正线性调频段波形,点频段波形,和负线性调频段波形,梯形连续波波形ST(t)公式如下:
其中k表示雷达发射梯形波的周期;f0为雷达发射信号中心频率;c为光速,对应波长为γ为线性调频段信号的调频斜率;T1为正线性调频段或负线性调频段的持续时间;因此线性调频段信号的带宽B满足B=γT1;T2为点频段的持续时间。
2,雷达回波解调后,正线性调频段、点频段,负线性调频段的雷达回波分别近似表示为SR1(t),SR2(t),SR3(t),波形公式如下:
Ri表示第i个目标相对雷达的初始距离;Vi表示第i个目标相对雷达的初始速度;N为目标数目;ρi为与雷达发射功率、天线增益、目标特性及距离相关的复常数;N1(t),N2(t),N3(t)表示接收机噪声。
3,对SR1(t),SR2(t),SR3(t)分别进行FFT变换:
其中G1表示线性调频段的相干累积增益,G2表示线性调频段的相干累积增益。
4,对SR1(f1),SR2(f2),SR3(f3),进行单元平均恒虚警检测,目标峰值点在正线性调频段信号、负线性调频段信号及点频段信号的位置分别可表示为:
5,通过可以确定所有目标的速度值的集合为
6,将正线性调频段信号第m个检测的目标与负线性调频段信号的第n个目标进行配对,得到可能的目标速度和距离分别为 将Vmn与Vi2集合中所有数值进行逐一比较,假设Vi2集合中与Vmn最接近数值为Vp。
7,如果其中K为常数,取值范围为[0,0.5],表示配对错误,令n=n+1,返回第6步;表示配对正确,输出该目标的正确速度Vmn和距离Rmn;如果m<N,令m=m+1,n=1,返回第6步;m≥N,密集目标测速测距结束。
下面以一个具体实施方式做验证:
假设待测的密集目标实际的情况是:N为30个。Ri实际值(i=1,2,...30)分别设置为(345.6253,344.5226,244.7528,171.5827,318.9544,98.4815,69.4933,276.6075,295.4367,157.4790,62.6288,262.8239,195.9677,201.8835,303.9411,98.0711,115.9553,114.5254,129.3818,292.8553,287.3594,258.2660,199.8911,204.4307,283.8437,229.5064,222.5229,178.0724,53.6040,59.2983)米。Vi实际值(i=1,2,...30)分别设置为(106.9501,118.6260,74.9322,115.1183,80.8346,121.1538,105.1768,79.7366,101.5018,123.1792,84.5725,81.7519,94.2987,84.9611,109.8040,125.9196,84.5647,100.9914,89.5191,77.1490,100.7747,123.0129,86.8843,117.8698,78.0492,74.9271,76.4098,118.3363,111.3115,108.6568)千米/小时。
在该实施方式中,k取值为1;f0取值为24×109赫兹;c取值为3×108米/秒;λ为0.0125米;γ取值为5×109赫兹/秒;T1取值为10毫秒;B为50兆赫兹;T2取值为20毫秒;G1取值为5×104;G2取值为1×105;便于计算,ρi均统一取值为1,N1(t),N2(t),N3(t)均满足均值为零,方差为100的高斯分布。
图1是该实施方式的雷达发射的梯形连续波的时频图。
图2是该实施方式的雷达测距测速的仿真结果。其中“*”表示目标真实距离和速度,“o”表示雷达测量获得的目标速度和距离,根据仿真结果显示,30个目标均获得精确的速度和距离信息,测速精度均优于1千米/小时,测距精度优于1米。
上述基于梯形连续波的雷达测速测距方法可以有效实现雷达对密集多目标的测距测速,从而极大提高雷达对密集多目标的监视性能。
前述对本发明的具体示例性实施方案的描述是为了说明和例证的目的。这些描述并非想将本发明限定为所公开的精确形式,并且很显然,根据上述教导,可以进行很多改变和变化。对示例性实施例进行选择和描述的目的在于解释本发明的特定原理及其实际应用,从而使得本领域的技术人员能够实现并利用本发明的各种不同的示例性实施方案以及各种不同的选择和改变。本发明的范围意在由权利要求书及其等同形式所限定。
Claims (7)
1.一种基于梯形连续波的雷达测速测距方法,用于测量多个目标相对雷达的速度和距离,其特征在于,步骤包括:
雷达发射梯形连续波信号,所述梯形连续波包括正线性调频段波形、点频段波形和负线性调频段波形;
雷达回波解调后,确定所有待测目标的实际速度值集合;
确定每个待测目标的可能速度值;以及
将所述每个待测目标的可能速度值与所述所有待测目标的实际速度值集合进行匹配得到待测目标的实际速度并且进一步得到待测目标的实际距离。
2.根据权利要求1所述的基于梯形连续波的雷达测速测距方法,其特征在于,梯形连续波波形ST(t)公式如下:
<mrow>
<msub>
<mi>S</mi>
<mi>T</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<msup>
<mi>e</mi>
<mrow>
<mi>j</mi>
<mn>2</mn>
<msub>
<mi>&pi;f</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mi>k</mi>
<mo>(</mo>
<mrow>
<mn>2</mn>
<msub>
<mi>T</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>T</mi>
<mn>2</mn>
</msub>
</mrow>
<mo>)</mo>
<mo>-</mo>
<mfrac>
<msub>
<mi>T</mi>
<mn>1</mn>
</msub>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>j</mi>
<mi>&pi;</mi>
<mi>&gamma;</mi>
<msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mi>k</mi>
<mo>(</mo>
<mrow>
<mn>2</mn>
<msub>
<mi>T</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>T</mi>
<mn>2</mn>
</msub>
</mrow>
<mo>)</mo>
<mo>-</mo>
<mfrac>
<msub>
<mi>T</mi>
<mn>1</mn>
</msub>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msup>
<mo>,</mo>
<mi>t</mi>
<mo>&Element;</mo>
<mo>&lsqb;</mo>
<mi>k</mi>
<mo>(</mo>
<mn>2</mn>
<msub>
<mi>T</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>T</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
<mo>,</mo>
<mi>k</mi>
<mrow>
<mo>(</mo>
<mrow>
<mn>2</mn>
<msub>
<mi>T</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>T</mi>
<mn>2</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>T</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msup>
<mi>e</mi>
<mrow>
<mi>j</mi>
<mn>2</mn>
<msub>
<mi>&pi;f</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mi>k</mi>
<mo>(</mo>
<mrow>
<mn>2</mn>
<msub>
<mi>T</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>T</mi>
<mn>2</mn>
</msub>
</mrow>
<mo>)</mo>
<mo>-</mo>
<msub>
<mi>T</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>T</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mn>1</mn>
</msub>
</mrow>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<mo>,</mo>
<mi>t</mi>
<mo>&Element;</mo>
<mo>&lsqb;</mo>
<mi>k</mi>
<mo>(</mo>
<mn>2</mn>
<msub>
<mi>T</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>T</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
<mo>+</mo>
<msub>
<mi>T</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<mi>k</mi>
<mrow>
<mo>(</mo>
<mrow>
<mn>2</mn>
<msub>
<mi>T</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>T</mi>
<mn>2</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>T</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>T</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msup>
<mi>e</mi>
<mrow>
<mi>j</mi>
<mn>2</mn>
<msub>
<mi>&pi;f</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mi>k</mi>
<mo>(</mo>
<mrow>
<mn>2</mn>
<msub>
<mi>T</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>T</mi>
<mn>2</mn>
</msub>
</mrow>
<mo>)</mo>
<mo>-</mo>
<msub>
<mi>T</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<mfrac>
<msub>
<mi>T</mi>
<mn>1</mn>
</msub>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>j</mi>
<mi>&pi;</mi>
<mi>&gamma;</mi>
<msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mi>k</mi>
<mo>(</mo>
<mrow>
<mn>2</mn>
<msub>
<mi>T</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>T</mi>
<mn>2</mn>
</msub>
</mrow>
<mo>)</mo>
<mo>-</mo>
<msub>
<mi>T</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<mfrac>
<msub>
<mi>T</mi>
<mn>1</mn>
</msub>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msup>
<mo>,</mo>
<mi>t</mi>
<mo>&Element;</mo>
<mo>&lsqb;</mo>
<mi>k</mi>
<mo>(</mo>
<mn>2</mn>
<msub>
<mi>T</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>T</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
<mo>+</mo>
<msub>
<mi>T</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>T</mi>
<mn>2</mn>
</msub>
<mo>,</mo>
<mi>k</mi>
<mrow>
<mo>(</mo>
<mrow>
<mn>2</mn>
<msub>
<mi>T</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>T</mi>
<mn>2</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mn>2</mn>
<msub>
<mi>T</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>T</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中k表示雷达发射梯形波的周期;f0为雷达发射信号中心频率;c为光速,波长γ为线性调频段信号的调频斜率;T1为正线性调频段或负线性调频段的持续时间;线性调频段信号的带宽B满足B=γT1;T2为点频段的持续时间。
3.根据权利要求2所述的基于梯形连续波的雷达测速测距方法,其特征在于,所述雷达回波的波形分为正线性调频段,点频段波形,负线性调频段,分别近似表示为SR1(t),SR2(t),SR3(t),波形公式如下:
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>R</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msub>
<mi>&rho;</mi>
<mi>i</mi>
</msub>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>j</mi>
<mfrac>
<mrow>
<mn>4</mn>
<mi>&pi;</mi>
</mrow>
<mi>&lambda;</mi>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<msub>
<mi>V</mi>
<mi>i</mi>
</msub>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>j</mi>
<mn>2</mn>
<mi>&pi;</mi>
<mi>&gamma;</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<msub>
<mi>V</mi>
<mi>i</mi>
</msub>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mi>c</mi>
</mfrac>
<mi>t</mi>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<msub>
<mi>V</mi>
<mi>i</mi>
</msub>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mi>c</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<mo>+</mo>
<msub>
<mi>N</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mi>t</mi>
<mo>&Element;</mo>
<mo>&lsqb;</mo>
<mo>-</mo>
<mfrac>
<msub>
<mi>T</mi>
<mn>1</mn>
</msub>
<mn>2</mn>
</mfrac>
<mo>,</mo>
<mfrac>
<msub>
<mi>T</mi>
<mn>1</mn>
</msub>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>R</mi>
<mn>2</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msub>
<mi>&rho;</mi>
<mi>i</mi>
</msub>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>j</mi>
<mfrac>
<mrow>
<mn>4</mn>
<mi>&pi;</mi>
</mrow>
<mi>&lambda;</mi>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<msub>
<mi>V</mi>
<mi>i</mi>
</msub>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<mo>+</mo>
<msub>
<mi>N</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mi>t</mi>
<mo>&Element;</mo>
<mo>&lsqb;</mo>
<mo>-</mo>
<mfrac>
<msub>
<mi>T</mi>
<mn>2</mn>
</msub>
<mn>2</mn>
</mfrac>
<mo>,</mo>
<mfrac>
<msub>
<mi>T</mi>
<mn>2</mn>
</msub>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>R</mi>
<mn>3</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msub>
<mi>&rho;</mi>
<mi>i</mi>
</msub>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>j</mi>
<mfrac>
<mrow>
<mn>4</mn>
<mi>&pi;</mi>
</mrow>
<mi>&lambda;</mi>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<msub>
<mi>V</mi>
<mi>i</mi>
</msub>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>j</mi>
<mn>2</mn>
<mi>&pi;</mi>
<mi>&gamma;</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<msub>
<mi>V</mi>
<mi>i</mi>
</msub>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mi>c</mi>
</mfrac>
<mi>t</mi>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>R</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<msub>
<mi>V</mi>
<mi>i</mi>
</msub>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mi>c</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
</mrow>
</msup>
<mo>+</mo>
<msub>
<mi>N</mi>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mi>t</mi>
<mo>&Element;</mo>
<mo>&lsqb;</mo>
<mo>-</mo>
<mfrac>
<msub>
<mi>T</mi>
<mn>1</mn>
</msub>
<mn>2</mn>
</mfrac>
<mo>,</mo>
<mfrac>
<msub>
<mi>T</mi>
<mn>1</mn>
</msub>
<mn>2</mn>
</mfrac>
<mo>)</mo>
</mrow>
Ri表示第i个目标相对雷达的初始距离;Vi表示第i个目标相对雷达的初始速度;N为目标数目;ρi为与雷达发射功率、天线增益、目标特性及距离相关的复常数;N1(t),N2(t),N3(t)表示接收机噪声。
4.根据权利要求3所述的基于梯形连续波的雷达测速测距方法,其特征在于,对所述SR1(t),SR2(t),SR3(t)分别进行FFT变换:
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>R</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msub>
<mi>G</mi>
<mn>1</mn>
</msub>
<msub>
<mi>&rho;</mi>
<mi>i</mi>
</msub>
<mi>sin</mi>
<mi> </mi>
<mi>c</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<mrow>
<msub>
<mi>f</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msub>
<mi>V</mi>
<mi>i</mi>
</msub>
</mrow>
<mi>&lambda;</mi>
</mfrac>
<mo>-</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msub>
<mi>&gamma;R</mi>
<mi>i</mi>
</msub>
</mrow>
<mi>c</mi>
</mfrac>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>N</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>f</mi>
<mn>1</mn>
</msub>
<mo>&Element;</mo>
<mo>&lsqb;</mo>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>T</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>,</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>T</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>R</mi>
<mn>2</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msub>
<mi>G</mi>
<mn>2</mn>
</msub>
<msub>
<mi>&rho;</mi>
<mi>i</mi>
</msub>
<mi>sin</mi>
<mi> </mi>
<mi>c</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mn>2</mn>
</msub>
<mo>(</mo>
<mrow>
<msub>
<mi>f</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msub>
<mi>V</mi>
<mi>i</mi>
</msub>
</mrow>
<mi>&lambda;</mi>
</mfrac>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>N</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>f</mi>
<mn>2</mn>
</msub>
<mo>&Element;</mo>
<mo>&lsqb;</mo>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>T</mi>
<mn>2</mn>
</msub>
</mfrac>
<mo>,</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>T</mi>
<mn>2</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>R</mi>
<mn>3</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mn>3</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msub>
<mi>G</mi>
<mn>1</mn>
</msub>
<msub>
<mi>&rho;</mi>
<mi>i</mi>
</msub>
<mi>sin</mi>
<mi> </mi>
<mi>c</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<mrow>
<msub>
<mi>f</mi>
<mn>3</mn>
</msub>
<mo>-</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msub>
<mi>V</mi>
<mi>i</mi>
</msub>
</mrow>
<mi>&lambda;</mi>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msub>
<mi>&gamma;R</mi>
<mi>i</mi>
</msub>
</mrow>
<mi>c</mi>
</mfrac>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>N</mi>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mn>3</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mi>f</mi>
<mn>3</mn>
</msub>
<mo>&Element;</mo>
<mo>&lsqb;</mo>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>T</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>,</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>T</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
其中G1表示线性调频段的相干累积增益,G2表示线性调频段的相干累积增益。
5.根据权利要求4所述的基于梯形连续波的雷达测速测距方法,其特征在于,对所述SR1(f1),SR2(f2),SR3(f3)进行单元平均恒虚警检测,检测到的目标峰值点在正线性调频段信号、点频段信号及负线性调频段信号的位置分别可表示为:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<msub>
<mi>f</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<msub>
<mi>i</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
<mo>=</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msub>
<mi>V</mi>
<msub>
<mi>i</mi>
<mn>1</mn>
</msub>
</msub>
</mrow>
<mi>&lambda;</mi>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msub>
<mi>&gamma;R</mi>
<msub>
<mi>i</mi>
<mn>1</mn>
</msub>
</msub>
</mrow>
<mi>c</mi>
</mfrac>
<mo>,</mo>
<msub>
<mi>i</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mi>N</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>f</mi>
<mn>3</mn>
</msub>
<mo>(</mo>
<msub>
<mi>i</mi>
<mn>3</mn>
</msub>
<mo>)</mo>
<mo>=</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msub>
<mi>V</mi>
<msub>
<mi>i</mi>
<mn>3</mn>
</msub>
</msub>
</mrow>
<mi>&lambda;</mi>
</mfrac>
<mo>-</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msub>
<mi>&gamma;R</mi>
<msub>
<mi>i</mi>
<mn>3</mn>
</msub>
</msub>
</mrow>
<mi>c</mi>
</mfrac>
<mo>,</mo>
<msub>
<mi>i</mi>
<mn>3</mn>
</msub>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mi>N</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>f</mi>
<mn>2</mn>
</msub>
<mo>(</mo>
<msub>
<mi>i</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
<mo>=</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msub>
<mi>V</mi>
<msub>
<mi>i</mi>
<mn>2</mn>
</msub>
</msub>
</mrow>
<mi>&lambda;</mi>
</mfrac>
<mo>,</mo>
<msub>
<mi>i</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mi>N</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
所述所有待测目标的实际速度值集合的集合为
6.根据权利要求4所述的基于梯形连续波的雷达测速测距方法,其特征在于,将所述正线性调频段信号第m个检测的目标与所述负线性调频段信号的第n个目标进行配对,确定所述待测目标的可能速度值为
7.根据权利要求5所述的基于梯形连续波的雷达测速测距方法,其特征在于,将Vmn与Vi2集合中所有数值逐一比较找到最接近数值Vp;如果其中K为常数,取值范围为[0,0.5],表示配对错误,令n=n+1,继续所述逐一比较过程;表示配对正确,输出该目标的正确速度Vmn和距离Rmn,所述如果m<N,令m=m+1,n=1,继续所述逐一比较过程;m≥N,密集目标测速测距结束。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711356001.XA CN108120975A (zh) | 2017-12-16 | 2017-12-16 | 基于梯形连续波的雷达测速测距方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711356001.XA CN108120975A (zh) | 2017-12-16 | 2017-12-16 | 基于梯形连续波的雷达测速测距方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108120975A true CN108120975A (zh) | 2018-06-05 |
Family
ID=62230028
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711356001.XA Pending CN108120975A (zh) | 2017-12-16 | 2017-12-16 | 基于梯形连续波的雷达测速测距方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108120975A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110361726A (zh) * | 2019-06-21 | 2019-10-22 | 广东工业大学 | 一种毫米波雷达测速方法 |
CN110361735A (zh) * | 2019-07-22 | 2019-10-22 | 成都纳雷科技有限公司 | 一种基于测速雷达的车辆测速方法及装置 |
WO2022142231A1 (zh) * | 2020-12-31 | 2022-07-07 | 上海禾赛科技有限公司 | 利用调频连续波进行探测的方法和激光雷达 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1457436A (zh) * | 2001-02-06 | 2003-11-19 | 富士通天株式会社 | Fm-cw雷达装置 |
CN101147082A (zh) * | 2005-03-21 | 2008-03-19 | 罗伯特·博世有限公司 | 用于多个目标的间距和相对速度测量的方法和装置 |
CN102998671A (zh) * | 2012-11-23 | 2013-03-27 | 北京理工大学 | 一种利用梯形调制波判断调频连续波测距时有无强信号干扰的方法 |
CN105487060A (zh) * | 2015-11-25 | 2016-04-13 | 上海无线电设备研究所 | 一种双通道四斜率调制的多目标提取方法 |
WO2017177972A1 (en) * | 2016-04-15 | 2017-10-19 | Mediatek Inc. | Radar interference mitigation method and apparatus |
-
2017
- 2017-12-16 CN CN201711356001.XA patent/CN108120975A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1457436A (zh) * | 2001-02-06 | 2003-11-19 | 富士通天株式会社 | Fm-cw雷达装置 |
CN101147082A (zh) * | 2005-03-21 | 2008-03-19 | 罗伯特·博世有限公司 | 用于多个目标的间距和相对速度测量的方法和装置 |
CN102998671A (zh) * | 2012-11-23 | 2013-03-27 | 北京理工大学 | 一种利用梯形调制波判断调频连续波测距时有无强信号干扰的方法 |
CN105487060A (zh) * | 2015-11-25 | 2016-04-13 | 上海无线电设备研究所 | 一种双通道四斜率调制的多目标提取方法 |
WO2017177972A1 (en) * | 2016-04-15 | 2017-10-19 | Mediatek Inc. | Radar interference mitigation method and apparatus |
Non-Patent Citations (1)
Title |
---|
薛效龙: "基于梯形波FMCW雷达差拍信号配对的研究", 《2013年全国微波毫米波会议论文集》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110361726A (zh) * | 2019-06-21 | 2019-10-22 | 广东工业大学 | 一种毫米波雷达测速方法 |
CN110361726B (zh) * | 2019-06-21 | 2022-12-16 | 广东工业大学 | 一种毫米波雷达测速方法 |
CN110361735A (zh) * | 2019-07-22 | 2019-10-22 | 成都纳雷科技有限公司 | 一种基于测速雷达的车辆测速方法及装置 |
WO2022142231A1 (zh) * | 2020-12-31 | 2022-07-07 | 上海禾赛科技有限公司 | 利用调频连续波进行探测的方法和激光雷达 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108415010B (zh) | 一种基于梯形lfmcw调制的雷达多目标检测方法 | |
US10379204B2 (en) | Method for calibrating a MIMO radar sensor for motor vehicles | |
CN104246534A (zh) | 用于借助雷达信号确定对象的间距和径向速度的方法和设备 | |
CN112834992B (zh) | 一种脉冲多普勒雷达的信号处理方法、装置及存储介质 | |
CN107255814B (zh) | 一种基于lfmsk波形的雷达目标检测方法 | |
CN105738889B (zh) | 一种调频连续波测速测距方法 | |
CN105807267A (zh) | 一种mimo雷达扩展目标的检测方法 | |
US11255956B2 (en) | Method for measuring the flow velocity of a medium | |
CN107037432A (zh) | 用于车辆的雷达装置及其确定目标的方法 | |
CN108120975A (zh) | 基于梯形连续波的雷达测速测距方法 | |
US11117518B2 (en) | Method for detecting an obstacle by means of reflected ultrasonic waves | |
CN108693523A (zh) | 基于锯齿波线性调频雷达多目标测距测速的方法及系统 | |
CN106842182A (zh) | 基于对称三角lfmcw雷达的多目标测速测距方法 | |
CN105676212A (zh) | 一种近程测距雷达系统和基于该系统的目标测量方法 | |
CN105022036A (zh) | 风廓线雷达风速测定方法 | |
CN106443660A (zh) | 一种线性调频连续波雷达多运动目标匹配方法 | |
CN105957355A (zh) | 一种车辆测速方法 | |
RU2649073C1 (ru) | Способ определения координат подводного объекта гидроакустической системой подводной навигации с юстировочным маяком | |
JPH0527018A (ja) | レーダ信号処理装置 | |
CN105158754A (zh) | 一种利用多输入单输出无线电系统进行目标定位的方法 | |
CN103901432B (zh) | 一种多观测节点下非合作目标的轨迹跟踪方法及系统 | |
RU2334244C1 (ru) | Способ определения местоположения источника радиоизлучения | |
CN104360333A (zh) | 一种同时校正一阶和二阶距离徙动的相参积累检测方法 | |
CN105988104B (zh) | 测距方法、测距装置、定位装置与定位方法 | |
CN105988105B (zh) | 测距方法、测距装置、定位装置与定位方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20180605 |