CN109164422A - 一种基于抖动重频的脉冲多普勒雷达解模糊测距方法 - Google Patents
一种基于抖动重频的脉冲多普勒雷达解模糊测距方法 Download PDFInfo
- Publication number
- CN109164422A CN109164422A CN201811125933.8A CN201811125933A CN109164422A CN 109164422 A CN109164422 A CN 109164422A CN 201811125933 A CN201811125933 A CN 201811125933A CN 109164422 A CN109164422 A CN 109164422A
- Authority
- CN
- China
- Prior art keywords
- data matrix
- distance
- doppler
- pulse
- target
- 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 27
- 230000000737 periodic effect Effects 0.000 claims abstract description 44
- 238000001514 detection method Methods 0.000 claims abstract description 37
- 238000001914 filtration Methods 0.000 claims abstract description 33
- 210000001367 artery Anatomy 0.000 claims abstract description 22
- 210000003462 vein Anatomy 0.000 claims abstract description 22
- 238000009825 accumulation Methods 0.000 claims abstract description 16
- 230000006835 compression Effects 0.000 claims abstract description 7
- 238000007906 compression Methods 0.000 claims abstract description 7
- 239000011159 matrix material Substances 0.000 claims description 134
- 241000208340 Araliaceae Species 0.000 claims description 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims description 2
- 235000003140 Panax quinquefolius Nutrition 0.000 claims description 2
- 235000008434 ginseng Nutrition 0.000 claims description 2
- 238000002592 echocardiography Methods 0.000 claims 1
- 238000007689 inspection Methods 0.000 claims 1
- 230000001427 coherent effect Effects 0.000 description 2
- 235000013399 edible fruits Nutrition 0.000 description 2
- 230000035485 pulse pressure Effects 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000002035 prolonged effect Effects 0.000 description 1
- 238000004088 simulation 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/28—Details of pulse systems
- G01S7/285—Receivers
- G01S7/292—Extracting wanted echo-signals
- G01S7/2923—Extracting wanted echo-signals based on data belonging to a number of consecutive radar periods
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于抖动重频的脉冲多普勒雷达解模糊测距方法,实现步骤为:(1)获取抖动重频下的雷达回波信号;(2)对两个脉组回波信号做脉冲压缩;(3)对两个脉组的脉冲压缩结果进行多普勒滤波;(4)按照最大距离模糊次数对两个脉组的多普勒滤波结果在距离维进行周期延拓;(5)对周期延拓结果进行相参积累;(6)进行门限检测,并进行目标判别,确定目标真实距离对应的距离单元;(7)计算目标真实距离,完成解模糊测距。本发明可以将两个抖动重频脉组的所有脉冲都积累起来,可以在解模糊的同时,提高检测概率,对于信噪比较小的目标也能够进行解模糊测距。
Description
技术领域
本发明属于雷达信号处理技术领域,特别涉及一种基于抖动重频的脉冲多普勒雷达解模糊测距方法,适用于PD雷达(脉冲多普勒雷达)对目标进行解模糊测距的场合。
背景技术
PD雷达通常采用中、高重复频率来工作,当目标距离超过雷达最大无模糊距离时,其回波出现在非本周期内,使得所测距离不是真实距离,因此需要进行解距离模糊。
在文献[谭贤四,李志淮,王红等.参差加抖动解模糊方法应用分析[J].现代雷达.2011,33(5)]中,作者提出了一种利用抖动重频解距离模糊的方法,其基本步骤如下:
(1)雷达发射两个脉组信号,这两个脉组信号的脉冲重复周期之间有很小的抖动;
(2)雷达接收回波信号,进行距离压缩,多普勒滤波,门限检测等处理,获得目标的视在距离信息;
(3)将目标的视在距离以不同的脉冲重复周期对应的最大不模糊距离进行周期延拓;
(4)选择一个适当宽度的窗口,按特定步长顺序移动窗口,当在同一窗口的目标数目等于2,该窗口出现的位置就是目标的真实距离。
上述方法是先进行目标信号检测,当检测到目标信号后再进行周期延拓,最后完成解模糊测距。该现有技术的缺陷在于,先进行目标信号检测,只能对两个脉组分别进行相参积累,没有充分利用所有的脉冲信息,因此信噪比会有一定的损失,当目标信号比较微弱时,有可能因为信噪比过低而检测不到目标,也就无法在后续过程中进行解距离模糊。
发明内容
针对上述现有技术的缺点,本发明的目的在于提供一种基于抖动重频的脉冲多普勒雷达解模糊测距方法,旨在提高信噪比,从而提高目标信号的检测概率以及解距离模糊的成功率。
为达到上述目的,本发明采用如下技术方案予以实现。
一种基于抖动重频的脉冲多普勒雷达解模糊测距方法,所述方法包括:
步骤1,获取来自被测目标的两组雷达回波信号,分别组成第一回波数据矩阵和第二回波数据矩阵;
步骤2,分别对所述第一回波数据矩阵和所述第二回波数据矩阵进行脉冲压缩,得到脉冲压缩后的第一数据矩阵和脉冲压缩后的第二数据矩阵;
步骤3,分别对所述脉冲压缩后的第一数据矩阵和脉冲压缩后的第二数据矩阵进行多普勒滤波,得到多普勒滤波后的第一数据矩阵和多普勒滤波后的第二数据矩阵;
步骤4,分别对所述多普勒滤波后的第一数据矩阵和所述多普勒滤波后的第二数据矩阵在距离维进行周期延拓,得到周期延拓后的第一数据矩阵和周期延拓后的第二数据矩阵;
步骤5,对所述周期延拓后的第一数据矩阵和所述周期延拓后的第二数据矩阵进行相参积累,得到相参积累后的数据矩阵,进而对所述相参积累后的数据矩阵进行检波,得到检波后的数据矩阵;
步骤6,对所述检波后的数据矩阵中的每一个多普勒通道的元素进行恒虚警检测,得到所有过门限的元素分别对应的多普勒通道及距离单元;根据所述所有过门限的元素分别对应的多普勒通道及距离单元,以及所述周期延拓后的第一数据矩阵和周期延拓后的第二数据矩阵,确定目标所在的多普勒通道及目标所在的距离单元;
步骤7,根据目标所在的距离单元,完成距离解模糊,得到目标的真实距离。
本发明与现有技术相比,具有以下优点:该方法先对滤波后的目标信号进行距离维拓展,然后进行目标信号的检测,最后完成解模糊测距。抖动重频意味着两个脉组的脉冲重复周期相差很小,所以两个脉组的多普勒通道是一致的,本发明在检测时将抖动重频的两个脉组中全部脉冲相参积累起来,避免了由于抖动重频将脉冲拆分而带来信噪比损失,可以提高目标信号的检测概率以及解距离模糊的成功率,对于微弱的目标信号,也可以完成解模糊测距,工程实现容易,有较强的实用性和推广应用价值。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明技术方案的实现流程框图;
图2为回波信噪比SNR=-7dB时采用现有技术得到的目标信号所在多普勒通道的距离图;
图3为对图2结果进行恒虚警检测和周期拓展的结果图;
图4为回波信噪比SNR=-7dB时采用本发明得到的目标信号所在多普勒通道的距离维拓展图;
图5为对图4结果进行恒虚警检测的结果图;
图6为回波信噪比SNR=-13dB时采用现有技术时得到的目标信号所在多普勒通道的距离图;
图7为对图6结果进行恒虚警检测和周期拓展的结果图;
图8为回波信噪比SNR=-13dB时采用本发明得到的目标信号所在多普勒通道的距离维拓展图;
图9为对图8结果进行恒虚警检测的结果图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明实施例提供一种基于抖动重频的脉冲多普勒雷达解模糊测距方法,为简单起见,本实施例为单目标解模糊情形,如图1所示,所述方法包括:
步骤1,获取来自被测目标的两组雷达回波信号,分别组成第一回波数据矩阵和第二回波数据矩阵。
步骤1具体包括:
(1a)雷达发射机先向被测目标发射第一组脉冲信号,所述第一组脉冲信号为以Tr1为脉冲重复周期的连续M个脉冲信号,再向被测目标发射第二组脉冲信号,所述第二组脉冲信号为以Tr2为脉冲重复周期的连续M个脉冲信号,其中,M≥2,Tr2=Tr1+Δt,Δt表示脉冲重复周期抖动差,一般取Δt=2us;
(1b)雷达接收机接收所述被测目标反射的两组雷达回波信号,并以采样频率fs对所述两组雷达回波信号进行采样,得到M×Q1维的第一回波数据矩阵X1和M×Q2维的第二回波数据矩阵X2,其中,M表示第一组脉冲信号和第二组脉冲信号分别包含的脉冲个数,Q1表示第一回波数据矩阵的距离单元总个数,Q2表示第二回波数据矩阵的距离单元总个数。
步骤2,分别对所述第一回波数据矩阵和所述第二回波数据矩阵进行脉冲压缩,得到脉冲压缩后的第一数据矩阵和脉冲压缩后的第二数据矩阵;
对第一回波数据矩阵X1中的每一行元素(即同一脉冲的数据)用相应的匹配滤波器进行滤波,并对半脉压点进行剔除,得到M×Q1维的脉冲压缩后的第一数据矩阵X1’,对第二回波数据矩阵X2中的每一行元素也用相应的匹配滤波器进行滤波,并对半脉压点进行剔除,得到M×Q2维的脉冲压缩后的第二数据据矩阵X2’。
步骤3,分别对所述脉冲压缩后的第一数据矩阵和脉冲压缩后的第二数据矩阵进行多普勒滤波,得到多普勒滤波后的第一数据矩阵和多普勒滤波后的第二数据矩阵;
对脉冲压缩后的第一数据矩阵X1’中的每一列元素(即同一距离单元的数据)做傅里叶变换,得到M×Q1维的数据矩阵Y1,对脉冲压缩后的第二数据矩阵X2’中的每一列元素(即同一距离单元的数据)做傅里叶变换,得到M×Q2维的数据矩阵Y2。
步骤4,分别对所述多普勒滤波后的第一数据矩阵和所述多普勒滤波后的第二数据矩阵在距离维进行周期延拓,得到周期延拓后的第一数据矩阵和周期延拓后的第二数据矩阵。
步骤4具体包括:
(4a)记分别对所述第一回波数据矩阵和所述第二回波数据矩阵进行脉冲压缩,得到脉冲压缩后的第一数据矩阵X1’和脉冲压缩后的第二数据矩阵X2’,所述脉冲压缩后的第一数据矩阵X1’和所述脉冲压缩后的第二数据矩阵X2’的维数分别为M×Q1维和M×Q2维;
记分别对所述脉冲压缩后的第一数据矩阵X1’和所述脉冲压缩后的第二数据矩阵X2’进行多普勒滤波,得到多普勒滤波后的第一数据矩阵Y1和多普勒滤波后的第二数据矩阵Y2,所述多普勒滤波后的第一数据矩阵Y1的维数为M×Q1维,所述多普勒滤波后的第二数据矩阵Y2的维数为M×Q2维;
(4b)将所述多普勒滤波后的第一数据矩阵Y1在距离域以Q1为周期延拓F次,得到周期延拓后的M×P1维的第一数据矩阵Z1,其中,P1=Q1×(F+1),F表示距离最大模糊次数,Rmax为雷达最大可探测距离,Ru=min[R1,R2],R1为第一组脉冲信号对应的最大不模糊距离,R2为第二组脉冲信号对应的最大不模糊距离,表示向上取整,min[·]表示取小操作;
(4c)将所述多普勒滤波后的第二数据矩阵Y2在距离域以Q2为周期延拓F次,得到周期延拓后的M×P2维的第二数据矩阵Z2,其中,P2=Q2×(F+1);
(4d)将所述周期延拓后的第一数据矩阵Z1保持不变,所述周期延拓后的第二数据矩阵Z2中第P1+1列至第P2列元素剔除,使得所述周期延拓后的第二数据矩阵Z2的维数与所述周期延拓后的第一数据矩阵Z1的维数均为M×P1维。
步骤5,对所述周期延拓后的第一数据矩阵和所述周期延拓后的第二数据矩阵进行相参积累,得到相参积累后的数据矩阵,进而对所述相参积累后的数据矩阵进行检波,得到检波后的数据矩阵。
步骤5具体包括:
(5a)获取M×P1维的相位因子矩阵φ:
其中,φm=2πfmMTr1,m=1,2,3....,M,π表示圆周率,fm表示第m个多普勒通道的中心频率,M表示第一组脉冲信号包含的脉冲个数,Tr1表示第一组脉冲信号的脉冲重复周期;
(5c)对周期延拓后的第一数据矩阵Z1和周期延拓后的第二数据矩阵Z2进行相参积累,获得相参积累后的数据矩阵Z’:
Z’=Z1+Z2·exp[-jφ]
其中,exp表示以e为底的指数,·表示两个矩阵中对应行对应列的元素进行相乘,j表示虚数单位;
(5d)对所述相参积累后的数据矩阵Z’进行检波,得到检波后的M×P1维的数据矩阵Z;
步骤6,对所述检波后的数据矩阵中的每一个多普勒通道的元素分别进行恒虚警检测,得到所有过门限的元素分别对应的多普勒通道及距离单元;根据所述所有过门限的元素分别对应的多普勒通道及距离单元,以及所述周期延拓后的第一数据矩阵和周期延拓后的第二数据矩阵,确定目标所在的多普勒通道及目标所在的距离单元。
步骤6具体包括:
(6a)设置中间变量m=0;
(6b)令m的值加1,对所述检波后的数据矩阵Z中的第m个多普勒通道(第m行)元素进行恒虚警检测,得到所有过门限元素对应的距离单元,记为其中im=0,1,2,....Im,im表示第m个多普勒通道中过门限的元素的序号,Im表示第m个多普勒通道中过门限的元素的总数,Im为非负整数;
(6c)如果第m个多普勒通道中没有过门限的元素,则执行步骤(6b),否则分别确定所述周期延拓后的第一数据矩阵Z1的第m个多普勒通道中第个距离单元对应的元素,记为并确定所述周期延拓后的第二数据矩阵Z2的第m个多普勒通道中第个距离单元对应的元素,记为计算得到一组比值其中im=0,1,2,....Im,im表示第m个多普勒通道中过门限的元素的序号,Im表示第m个多普勒通道中过门限的元素的总数,Im为非负整数,|·|表示求模操作;
(6d)在所述一组比值中确定处于预设范围[1-Δ,1+Δ]中的比值对应的多普勒通道及距离单元为目标的多普勒通道及目标所对应的距离单元,将目标对应的距离单元记为r,其中Δ为比值误差,可根据实际的两组回波脉冲信号具体选取,一般取0<Δ<0.3;
(6e)判断m=M是否成立,若不成立,则执行(6b),否则执行步骤7;
步骤7,根据目标所在的距离单元,完成距离解模糊,得到目标的真实距离。
步骤7具体为:
计算目标的真实距离Rtrue=rρr,其中,ρr表示距离采样率,其中,c为光速,fs为采样频率。
下面结合仿真实验,对本发明的技术效果作进一步说明。
仿真条件:
假设两个脉组的脉冲重复周期Tr1=1200us,Tr2=1202us,两个脉组中的脉冲个数均为M=32,目标相对雷达的径向运动速度v=100m/s,目标距离雷达的距离R0=300km,距离最大模糊次数为4,发射的线性调频信号带宽B=2MHz,时宽T=100us,雷达载频f0=600MHz,信噪比SNR=-7dB或-13dB,这里的信噪比指的是未经任何处理的原始回波信噪比,距离向采样率fs=B。
仿真内容:
在仿真条件下,利用现有技术和本发明分别将对目标信号进行解模糊测距。当信噪比为-7dB时,采用现有技术得到目标信号所在多普勒通道的距离图如图2所示,对图2结果进行恒虚警检测和周期拓展的结果如图3所示,采用本发明方法得到的目标信号所在多普勒通道的距离拓展图如图4所示,对图4结果进行恒虚警检测的结果如图5所示;当信噪比为-13dB时,采用现有技术得到目标信号所在多普勒通道的距离图如图6所示,对图6结果进行恒虚警检测和周期拓展的结果如图7所示,采用本发明方法得到的目标信号所在多普勒通道的距离拓展图如图8所示,对图8结果进行恒虚警检测的结果如图9所示。
仿真结果分析:
图2为利用现有技术得到的两个抖动重频脉组的目标信号所在多普勒通道的距离图,对图2结果进行恒虚警检测后再在距离维进行周期延拓,得到图3,由图3可知两个脉组的距离维拓展图过门限的距离单元均为5个,采用2个距离单元作为滑窗宽度,只有在4000距离单元处滑动窗口内的目标数目等于2,因此目标真实距离对应的距离单元为4000,计算可得目标真实距离为300km,完成解模糊测距。
图4为利用本发明方法得到的目标信号所在多普勒通道的距离维拓展图,对图4结果进行恒虚警检测,得到图5,由图5可以看出过门限的距离单元有3个,分别在第4000、11200、11212个距离单元处,分别找出脉组一目标信号所在多普勒通道和脉组二目标信号所在多普勒通道在这3个距离单元处的幅度值,得到在11200距离单元处的幅度比值为18.1698,11212距离单元处的幅度比值为0.0985,而在4000距离单元出两个脉组的幅度比值为1.0019,因此目标真实距离对应的距离单元为4000,计算可得目标真实距离为300km,完成解模糊测距。
由图7和图9结果可得,当信噪比很低时,现有技术已经无法检测到目标信息,更无法进行解模糊测距,而本发明依然可以检测到目标信息,此时只有第4000个距离单元处的元素过检测门限,且在第4000个距离单元处两个脉组的幅度比值为0.9843,可得目标真实距离为300km,因此本发明在信噪比较低时依然可以正确的完成解模糊测距。
本领域普通技术人员可以理解:实现上述方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完成,前述的程序可以存储于计算机可读取存储介质中,该程序在执行时,执行包括上述方法实施例的步骤;而前述的存储介质包括:ROM、RAM、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以所述权利要求的保护范围为准。
Claims (6)
1.一种基于抖动重频的脉冲多普勒雷达解模糊测距方法,其特征在于,所述方法包括:
步骤1,获取来自被测目标的两组雷达回波信号,分别组成第一回波数据矩阵和第二回波数据矩阵;
步骤2,分别对所述第一回波数据矩阵和所述第二回波数据矩阵进行脉冲压缩,得到脉冲压缩后的第一数据矩阵和脉冲压缩后的第二数据矩阵;
步骤3,分别对所述脉冲压缩后的第一数据矩阵和脉冲压缩后的第二数据矩阵进行多普勒滤波,分别得到多普勒滤波后的第一数据矩阵和多普勒滤波后的第二数据矩阵;
步骤4,分别对所述多普勒滤波后的第一数据矩阵和所述多普勒滤波后的第二数据矩阵在距离维进行周期延拓,得到周期延拓后的第一数据矩阵和周期延拓后的第二数据矩阵;
步骤5,对所述周期延拓后的第一数据矩阵和所述周期延拓后的第二数据矩阵进行相参积累,得到相参积累后的数据矩阵,进而对所述相参积累后的数据矩阵进行检波,得到检波后的数据矩阵;
步骤6,对所述检波后的数据矩阵中的每一个多普勒通道的元素进行恒虚警检测,得到所有过门限的元素分别对应的多普勒通道及距离单元;根据所述所有过门限的元素分别对应的多普勒通道及距离单元,以及所述周期延拓后的第一数据矩阵和周期延拓后的第二数据矩阵,确定目标所在的多普勒通道及目标所在的距离单元;
步骤7,根据目标所在的距离单元,完成距离解模糊,得到目标的真实距离。
2.根据权利要求1所述的一种基于抖动重频的脉冲多普勒雷达解模糊测距方法,其特征在于,步骤1具体包括:
(1a)雷达发射机先向被测目标发射第一组脉冲信号,所述第一组脉冲信号为以Tr1为脉冲重复周期的连续M个脉冲信号,再向被测目标发射第二组脉组信号,所述第二组脉冲信号为以Tr2为脉冲重复周期的连续M个脉冲信号,其中,M≥2,Tr2=Tr1+Δt,Δt表示脉冲重复周期抖动差,且Δt大于零;
(1b)雷达接收机接收所述被测目标反射的两组雷达回波信号,并以采样频率fs对所述两组雷达回波信号进行采样,得到M×Q1维的第一回波数据矩阵X1和M×Q2维的第二回波数据矩阵X2,其中,M表示第一组脉冲信号和第二组脉冲信号分别包含的脉冲个数,Q1表示第一回波数据矩阵的距离单元总个数,Q2表示第二回波数据矩阵的距离单元总个数。
3.根据权利要求1所述的一种基于抖动重频的脉冲多普勒雷达解模糊测距方法,其特征在于,步骤4具体包括:
(4a)记分别对所述第一回波数据矩阵和所述第二回波数据矩阵进行脉冲压缩,得到脉冲压缩后的第一数据矩阵X1’和脉冲压缩后的第二数据矩阵X2’,所述脉冲压缩后的第一数据矩阵X1’和所述脉冲压缩后的第二数据矩阵X2’的维数分别为M×Q1维和M×Q2维;
记分别对所述脉冲压缩后的第一数据矩阵X1’和所述脉冲压缩后的第二数据矩阵X2’进行多普勒滤波,得到多普勒滤波后的第一数据矩阵Y1和多普勒滤波后的第二数据矩阵Y2,所述多普勒滤波后的第一数据矩阵Y1的维数为M×Q1维,所述多普勒滤波后的第二数据矩阵Y2的维数为M×Q2维;
(4b)将所述多普勒滤波后的第一数据矩阵Y1在距离域以Q1为周期延拓F次,得到周期延拓后的M×P1维的第一数据矩阵Z1,其中,P1=Q1×(F+1),F表示距离最大模糊次数,Rmax为雷达最大可探测距离,Ru=min[R1,R2],R1为第一组脉冲信号对应的最大不模糊距离,R2为第二组脉冲信号对应的最大不模糊距离,表示向上取整,min[·]表示取小操作;
(4c)将所述多普勒滤波后的第二数据矩阵Y2在距离域以Q2为周期延拓F次,得到周期延拓后的M×P2维的第二数据矩阵Z2,其中,P2=Q2×(F+1);
(4d)将所述周期延拓后的第一数据矩阵Z1保持不变,所述周期延拓后的第二数据矩阵Z2中第P1+1列至第P2列元素剔除,使得所述周期延拓后的第二数据矩阵Z2的维数与所述周期延拓后的第一数据矩阵Z1的维数均为M×P1维。
4.根据权利要求1所述的一种基于抖动重频的脉冲多普勒雷达解模糊测距方法,其特征在于,步骤5具体包括:
(5a)获取M×P1维的相位因子矩阵φ:
其中,φm=2πfmMTr1,m=1,2,3....,M,π表示圆周率,fm表示第m个多普勒通道的中心频率,M表示第一组脉冲信号包含的脉冲个数,Tr1表示第一组脉冲信号的脉冲重复周期;
(5c)对周期延拓后的第一数据矩阵Z1和周期延拓后的第二数据矩阵Z2进行相参积累,获得相参积累后的数据矩阵Z’:
Z’=Z1+Z2·exp[-jφ]
其中,exp表示以e为底的指数,·表示两个矩阵中对应行对应列的元素进行相乘,j表示虚数单位;
(5d)对所述相参积累后的数据矩阵Z’进行检波,得到检波后的M×P1维的数据矩阵Z。
5.根据权利要求1所述的一种基于抖动重频的脉冲多普勒雷达解模糊测距方法,其特征在于,步骤6具体包括:
(6a)设置中间变量m=0;
(6b)令m的值加1,对所述检波后的数据矩阵Z中的第m个多普勒通道元素进行恒虚警检测,得到所有过门限元素对应的距离单元,记为其中im=0,1,2,....Im,im表示第m个多普勒通道中过门限的元素的序号,Im表示第m个多普勒通道中过门限的元素的总数,Im为非负整数;
(6c)如果第m个多普勒通道中没有过门限的元素,则返回子步骤(6b);否则分别确定所述周期延拓后的第一数据矩阵Z1的第m个多普勒通道中第个距离单元对应的元素,记为并确定所述周期延拓后的第二数据矩阵Z2的第m个多普勒通道中第个距离单元对应的元素,记为计算得到一组比值其中im=0,1,2,....Im,im表示第m个多普勒通道中过门限的元素的序号,Im表示第m个多普勒通道中过门限的元素的总数,Im为非负整数,|·|表示求模操作;
(6d)在所述一组比值中确定处于预设范围[1-Δ,1+Δ]中的比值对应的多普勒通道及距离单元,作为目标的多普勒通道及目标所对应的距离单元,将目标所对应的距离单元记为r,其中Δ为比值误差;
(6e)当m<M时返回子步骤(6b)。
6.根据权利要求1所述的一种基于抖动重频的脉冲多普勒雷达解模糊测距方法,其特征在于,步骤7具体为:
计算目标的真实距离Rtrue=rρr,其中,ρr表示距离采样率,其中,c为光速,fs为采样频率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811125933.8A CN109164422A (zh) | 2018-09-26 | 2018-09-26 | 一种基于抖动重频的脉冲多普勒雷达解模糊测距方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811125933.8A CN109164422A (zh) | 2018-09-26 | 2018-09-26 | 一种基于抖动重频的脉冲多普勒雷达解模糊测距方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109164422A true CN109164422A (zh) | 2019-01-08 |
Family
ID=64880414
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811125933.8A Pending CN109164422A (zh) | 2018-09-26 | 2018-09-26 | 一种基于抖动重频的脉冲多普勒雷达解模糊测距方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109164422A (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109814073A (zh) * | 2019-01-21 | 2019-05-28 | 西安电子科技大学 | 一种mtd雷达解模糊测速方法 |
CN110109078A (zh) * | 2019-05-05 | 2019-08-09 | 西安电子科技大学 | 基于相参频率捷变雷达的速度解模糊方法 |
CN110673105A (zh) * | 2019-09-25 | 2020-01-10 | 武汉滨湖电子有限责任公司 | 一种脉冲多普勒雷达解速度模糊的方法 |
CN111044987A (zh) * | 2019-12-26 | 2020-04-21 | 成都纳雷科技有限公司 | 基于汽车雷达的环境目标速度解模糊方法、***及介质 |
CN111142086A (zh) * | 2020-01-19 | 2020-05-12 | 山东省科学院自动化研究所 | 一种pd雷达幅度抖动抑制方法、时间抖动检测方法及*** |
CN111289953A (zh) * | 2020-01-14 | 2020-06-16 | 北京理工大学 | 一种基于模糊矩阵更新的天基雷达距离/速度解模糊方法 |
CN113050058A (zh) * | 2021-03-23 | 2021-06-29 | 北京环境特性研究所 | 一种基于周期延拓的运动目标多普勒频移获得方法及装置 |
CN113791404A (zh) * | 2021-09-15 | 2021-12-14 | 电子科技大学长三角研究院(衢州) | 一种基于正交频分信号的雷达解模糊及遮挡的方法 |
CN116930907A (zh) * | 2023-09-15 | 2023-10-24 | 厦门大学 | 一种ddma-mimo雷达的目标相参积累检测方法 |
CN117471449A (zh) * | 2023-12-27 | 2024-01-30 | 中国电子科技集团公司第十四研究所 | 一种适用于机动目标的单组pd跟踪方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070046526A1 (en) * | 2005-08-26 | 2007-03-01 | Vaisala Oyj | Method for using pulse compression in weather radar |
CN101470202A (zh) * | 2007-12-26 | 2009-07-01 | 清华大学 | 一种脉冲多普勒雷达***及其信号处理方法 |
CN101738606A (zh) * | 2008-11-21 | 2010-06-16 | 清华大学 | 基于广义多普勒滤波器组的雷达目标相参积累检测方法 |
CN107783093A (zh) * | 2017-10-25 | 2018-03-09 | 哈尔滨工业大学 | 一种基于单重复频率脉冲雷达的解距离模糊和距离遮挡的方法 |
CN108535730A (zh) * | 2018-04-16 | 2018-09-14 | 青海大学 | 一种多普勒气象雷达解速度模糊方法和*** |
-
2018
- 2018-09-26 CN CN201811125933.8A patent/CN109164422A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070046526A1 (en) * | 2005-08-26 | 2007-03-01 | Vaisala Oyj | Method for using pulse compression in weather radar |
CN101470202A (zh) * | 2007-12-26 | 2009-07-01 | 清华大学 | 一种脉冲多普勒雷达***及其信号处理方法 |
CN101738606A (zh) * | 2008-11-21 | 2010-06-16 | 清华大学 | 基于广义多普勒滤波器组的雷达目标相参积累检测方法 |
CN107783093A (zh) * | 2017-10-25 | 2018-03-09 | 哈尔滨工业大学 | 一种基于单重复频率脉冲雷达的解距离模糊和距离遮挡的方法 |
CN108535730A (zh) * | 2018-04-16 | 2018-09-14 | 青海大学 | 一种多普勒气象雷达解速度模糊方法和*** |
Non-Patent Citations (1)
Title |
---|
王玮: "毫米波步进频雷达一维成像及目标检测算法研究", 《中国优秀硕士论文全文数据库信息科技辑》 * |
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109814073A (zh) * | 2019-01-21 | 2019-05-28 | 西安电子科技大学 | 一种mtd雷达解模糊测速方法 |
CN109814073B (zh) * | 2019-01-21 | 2022-11-11 | 西安电子科技大学 | 一种mtd雷达解模糊测速方法 |
CN110109078A (zh) * | 2019-05-05 | 2019-08-09 | 西安电子科技大学 | 基于相参频率捷变雷达的速度解模糊方法 |
CN110109078B (zh) * | 2019-05-05 | 2022-04-08 | 西安电子科技大学 | 基于相参频率捷变雷达的速度解模糊方法 |
CN110673105B (zh) * | 2019-09-25 | 2021-12-10 | 武汉滨湖电子有限责任公司 | 一种脉冲多普勒雷达解速度模糊的方法 |
CN110673105A (zh) * | 2019-09-25 | 2020-01-10 | 武汉滨湖电子有限责任公司 | 一种脉冲多普勒雷达解速度模糊的方法 |
CN111044987A (zh) * | 2019-12-26 | 2020-04-21 | 成都纳雷科技有限公司 | 基于汽车雷达的环境目标速度解模糊方法、***及介质 |
CN111289953A (zh) * | 2020-01-14 | 2020-06-16 | 北京理工大学 | 一种基于模糊矩阵更新的天基雷达距离/速度解模糊方法 |
CN111142086B (zh) * | 2020-01-19 | 2021-10-08 | 山东省科学院自动化研究所 | 一种pd雷达幅度抖动抑制方法、时间抖动检测方法及*** |
CN111142086A (zh) * | 2020-01-19 | 2020-05-12 | 山东省科学院自动化研究所 | 一种pd雷达幅度抖动抑制方法、时间抖动检测方法及*** |
CN113050058B (zh) * | 2021-03-23 | 2022-08-30 | 北京环境特性研究所 | 一种基于周期延拓的运动目标多普勒频移获得方法及装置 |
CN113050058A (zh) * | 2021-03-23 | 2021-06-29 | 北京环境特性研究所 | 一种基于周期延拓的运动目标多普勒频移获得方法及装置 |
CN113791404A (zh) * | 2021-09-15 | 2021-12-14 | 电子科技大学长三角研究院(衢州) | 一种基于正交频分信号的雷达解模糊及遮挡的方法 |
WO2023040858A1 (zh) * | 2021-09-15 | 2023-03-23 | 电子科技大学长三角研究院(衢州) | 一种基于正交频分信号的雷达解模糊及遮挡的方法 |
CN113791404B (zh) * | 2021-09-15 | 2024-05-07 | 电子科技大学长三角研究院(衢州) | 一种基于正交频分信号的雷达解模糊及遮挡的方法 |
CN116930907A (zh) * | 2023-09-15 | 2023-10-24 | 厦门大学 | 一种ddma-mimo雷达的目标相参积累检测方法 |
CN116930907B (zh) * | 2023-09-15 | 2024-01-12 | 厦门大学 | 一种ddma-mimo雷达的目标相参积累检测方法 |
CN117471449A (zh) * | 2023-12-27 | 2024-01-30 | 中国电子科技集团公司第十四研究所 | 一种适用于机动目标的单组pd跟踪方法 |
CN117471449B (zh) * | 2023-12-27 | 2024-03-22 | 中国电子科技集团公司第十四研究所 | 一种适用于机动目标的单组pd跟踪方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109164422A (zh) | 一种基于抖动重频的脉冲多普勒雷达解模糊测距方法 | |
CN106646446B (zh) | 脉冲压缩捷变频雷达的动目标检测方法 | |
CN109814073A (zh) | 一种mtd雷达解模糊测速方法 | |
AU2006306650B2 (en) | Sonar system and method providing low probability of impact on marine mammals | |
CN204495996U (zh) | 宽带声纳接收器 | |
Thode et al. | Automated detection and localization of bowhead whale sounds in the presence of seismic airgun surveys | |
US7250900B2 (en) | System and method for sidelobe reduction using detect-and-subtract techniques | |
CN107450055B (zh) | 基于离散线性调频傅立叶变换的高速机动目标检测方法 | |
CN109814075A (zh) | 一种mti雷达运动目标检测方法 | |
CN109581313A (zh) | 一种脉组频率捷变雷达的目标检测方法 | |
CN106597411A (zh) | 雷达信号处理方法 | |
US8121222B2 (en) | Systems and methods for construction of time-frequency surfaces and detection of signals | |
WO2013087871A2 (en) | Method and device for separating seismic signals from seismic sources | |
US11988772B2 (en) | Remote recovery of acoustic signals from passive sources | |
CN106872985B (zh) | 基于改进短时傅里叶变换的火箭弹炮口速度测量方法 | |
Catherall et al. | High resolution spectrograms using a component optimized short-term fractional Fourier transform | |
CN110007282A (zh) | 一种连续波体制1-比特雷达目标重构问题降维方法 | |
US8358233B2 (en) | Radar target detection process | |
US20060235635A1 (en) | Apparatus and method for performing the delay estimation of signals propagating through an environment | |
EP3111252B1 (en) | Systems and methods for pre-correlation noise attentuation using random sweeps in seismic data processing | |
CN105092022B (zh) | 基于抹香鲸叫声的水下目标测距单脉冲筛选方法及装置 | |
CN110244290A (zh) | 一种距离扩展目标的检测方法 | |
CN111443333B (zh) | 一种基于频谱综合的多方位匹配信号生成方法 | |
CN108784744A (zh) | 用于剪切波成像的可变焦点 | |
Yu et al. | A double-threshold target detection method in detecting low slow small target |
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: 20190108 |