CN115166681B - 调频连续波信号体制穿墙雷达目标检测快速方法及系统 - Google Patents

调频连续波信号体制穿墙雷达目标检测快速方法及系统 Download PDF

Info

Publication number
CN115166681B
CN115166681B CN202211091992.4A CN202211091992A CN115166681B CN 115166681 B CN115166681 B CN 115166681B CN 202211091992 A CN202211091992 A CN 202211091992A CN 115166681 B CN115166681 B CN 115166681B
Authority
CN
China
Prior art keywords
sequence
fourier transform
preset condition
distance
frequency
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202211091992.4A
Other languages
English (en)
Other versions
CN115166681A (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.)
Beijing Lsj Technology Development Co ltd
Huaxia Skyworth Beijing Security Technology Research Institute Co ltd
Wuhan Xinlang Photoelectric Technology Co ltd
Original Assignee
Beijing Lsj Technology Development Co ltd
Huaxia Skyworth Beijing Security Technology Research Institute Co ltd
Wuhan Xinlang Photoelectric Technology Co ltd
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 Beijing Lsj Technology Development Co ltd, Huaxia Skyworth Beijing Security Technology Research Institute Co ltd, Wuhan Xinlang Photoelectric Technology Co ltd filed Critical Beijing Lsj Technology Development Co ltd
Priority to CN202211091992.4A priority Critical patent/CN115166681B/zh
Publication of CN115166681A publication Critical patent/CN115166681A/zh
Application granted granted Critical
Publication of CN115166681B publication Critical patent/CN115166681B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • 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/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • 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/887—Radar or analogous systems specially adapted for specific applications for detection of concealed objects, e.g. contraband or weapons
    • G01S13/888—Radar or analogous systems specially adapted for specific applications for detection of concealed objects, e.g. contraband or weapons through wall detection

Landscapes

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

Abstract

本发明涉及一种调频连续波信号体制穿墙雷达目标检测快速方法及系统，方法包括：采用零中频去调频方式接收FMCW接收机按照设定频率所发射的每帧发射信息，获得对应的差拍信号；根据获得的差拍信号，获取每一差拍信号分别所对应的回波距离像；针对任意相邻两帧回波距离像，进行差分获得差分距离像的序列；筛选出满足第一预先设定条件的差分距离像；针对满足第一预先设定条件的差分距离像对应距离处的复散射值沿序列方向进行傅里叶变换，得到傅里叶变换后的复散射值的函数；针对傅里叶变换后复散射值的函数沿序列方向求取满足第二预设条件的傅里叶变换后复散射值的函数值点，并将满足第二预设条件的函数值点所对应的距离处作为有效目标。

Description

调频连续波信号体制穿墙雷达目标检测快速方法及系统
技术领域
本发明涉及穿墙雷达技术领域，尤其涉及一种调频连续波信号体制穿墙雷达目标检测快速方法及系统。
背景技术
穿墙雷达可穿透墙壁对墙后运动目标进行探测，是城市环境中抢险救援等应用中非常重要的探测感知手段。穿墙雷达从信号体制上分类，可以分为冲激脉冲体制和连续波体制等，其中连续波体制又可以分为线性调频连续波（FMCW）、步进频率（SFCW）、伪随机相位编码等若干种。由于连续波体制可以保证脉冲间的严格相参，有利于通过多帧积累提高目标探测性能，或通过多个接收天线提高测向能力，因此信号处理增益潜力巨大，是穿墙雷达重要的技术发展方向之一。
连续波穿墙雷达特别是FMCW穿墙雷达发射功率很低，特别是墙壁衰减很大造成目标回波微弱，对目标检测带来了很大挑战。通过前后两帧数据进行差分检测的方法，虽然可以消除场景中静止目标回波的影响，但穿透较厚墙壁后的微弱运动目标信噪比很低，特别是穿墙后室内场景中有植物、电扇等固定位置的动目标，仅仅通过运动信息难以区分，仍然会造成大量的虚警。同时，考虑到穿墙雷达往往采用便携或手持的方式工作，处理器处理速度有限，因此很多运算量需求较大的方法也无法适用。为了降低整机的体积、重量和功耗，检测算法的运算量也应该尽可能低，使其可以在单片机等运算资源受限的处理器上运行。
发明内容
（一）要解决的技术问题
鉴于现有技术的上述缺点、不足，本发明提供一种调频连续波信号体制穿墙雷达目标检测快速方法及系统，其解决了现有技术中连续波穿墙雷达特别是FMCW穿墙雷达发射功率很低，特别是墙壁衰减很大造成目标回波微弱，对目标检测不是特别精确的技术问题。
（二）技术方案
为了达到上述目的，本发明采用的主要技术方案包括：
第一方面，本发明实施例提供一种调频连续波信号体制穿墙雷达目标检测快速方法，包括：
S1、采用零中频去调频方式接收FMCW接收机按照设定频率所发射的每帧发射信息，并获得与每帧发射信息分别一一对应的差拍信号；
S2、根据获得的与每帧发射信息分别一一对应的差拍信号，获取每一差拍信号分别所对应的回波距离像；
S3、针对任意相邻两帧回波距离像，进行差分获得差分距离像的序列；
所述差分距离像的序列包括：依次相邻两帧回波距离像之间的差分值；
S4、在差分距离像的序列中，筛选出满足第一预先设定条件的差分距离像；
S5、针对满足第一预先设定条件的差分距离像对应距离处的复散射值沿序列方向进行傅里叶变换，得到傅里叶变换后的复散射值的函数；
S6、针对傅里叶变换后复散射值的函数沿序列方向求取满足第二预设条件的傅里叶变换后复散射值的函数值点，并将满足第二预设条件的函数值点所对应的距离处作为有效目标。
优选地，S1具体包括：
将FMCW接收机按照设定频率所发射的任一帧发射信息
Figure 613360DEST_PATH_IMAGE001
的耦合放大一路 作为本振信号，并将该本振信号与FMCW接收机接收到的与所述发射信息对应的回波信号
Figure 899984DEST_PATH_IMAGE002
进行混频和低通滤波，获得与该发射信息
Figure 56421DEST_PATH_IMAGE001
所对应的零中频的差拍信号
Figure 668668DEST_PATH_IMAGE003
。
优选地，
其中，
Figure 143512DEST_PATH_IMAGE001
为FMCW接收机在t时刻所发射的第m个发射信息；
Figure 735293DEST_PATH_IMAGE004
；
其中，FMCW接收机的发射信息幅度归一化为1；
exp是指数；
j是虚部符号；
f 0 是线性调频起始频率；
T为扫描时间长度；
K为调频斜率，K=B/T；
B为调频信号带宽；
t为时刻；
m为发射信息的序号；
Figure 244771DEST_PATH_IMAGE002
为FMCW接收机所接收到的发射信息为
Figure 762340DEST_PATH_IMAGE001
的回波信号；
其中，
Figure 724480DEST_PATH_IMAGE005
；
Figure 119952DEST_PATH_IMAGE006
为发射信息为
Figure 218358DEST_PATH_IMAGE001
的回波信号的延迟；
A为回波强度。
优选地，
与所述发射信息为
Figure 906828DEST_PATH_IMAGE001
对应的差拍信号为
Figure 857729DEST_PATH_IMAGE003
；
其中，
Figure 24268DEST_PATH_IMAGE007
；
其中，
Figure 773918DEST_PATH_IMAGE008
；
Figure 633290DEST_PATH_IMAGE009
；
Figure 71487DEST_PATH_IMAGE010
为目标起始距离；c表示光速；
Figure 572875DEST_PATH_IMAGE011
；
其中，
Figure 645873DEST_PATH_IMAGE012
；
Figure 443190DEST_PATH_IMAGE013
；
Figure 601639DEST_PATH_IMAGE014
为多普勒频率。
优选地，所述S2具体包括：
所述差拍信号为
Figure 109981DEST_PATH_IMAGE003
所对应的回波距离像为
Figure 568644DEST_PATH_IMAGE015
；
Figure 5704DEST_PATH_IMAGE016
；
其中，
Figure 917028DEST_PATH_IMAGE017
；
f是t做傅里叶变换后的频率变量表示；
Figure 494640DEST_PATH_IMAGE018
。
优选地，所述S3中，
所述差分距离像的序列包括：依次相邻两帧回波距离像之间的差分值；
其中，差分距离像的序列是：
Figure 320993DEST_PATH_IMAGE019
；
Figure 224227DEST_PATH_IMAGE020
是第m对相邻两帧回波距离像之间的差分值；
N是预先设定值；
其中，
Figure 357268DEST_PATH_IMAGE021
。
优选地，所述S4中的所述第一预先设定条件为：差分距离像大于Tp；
其中，
Figure 738571DEST_PATH_IMAGE022
；
p为预先设定的正实数。
优选地，所述S5具体包括：
S51、基于满足第一预先设定条件的差分距离像，获取满足第一预先设定条件的差分距离像对应距离处的差拍信号；
Figure 407712DEST_PATH_IMAGE023
；
Figure 950689DEST_PATH_IMAGE024
指的是满足第一预先设定条件的差分距离像对应距离处的差拍信号中的 与发射信息的序号m对应的差拍信号；
M指的是满足第一预先设定条件的差分距离像的数量；
S52、对满足第一预先设定条件的差分距离像对应距离处的差拍信号沿序列方向进行傅里叶变换可得第一公式：
Figure 305447DEST_PATH_IMAGE025
；
Figure 991905DEST_PATH_IMAGE026
具体表示指的是满足第一预先设定条件的差分距离像对应距离处的差拍 信号中的与发射信息的序号m对应的差拍信号沿序列方向进行傅里叶变换的结果；
Figure 748509DEST_PATH_IMAGE027
是序号m发射信息的对应频率域变量表示；
S53、通过补偿函数抑制掉第一公式中
Figure 462387DEST_PATH_IMAGE028
，得到沿序列方向的傅里叶变换后 的复散射值的函数；
其中，傅里叶变换后的复散射值的函数为：
Figure 570020DEST_PATH_IMAGE029
。
优选地，
所述第二预设条件为：傅里叶变换后复散射值的函数沿序列方向中满足高于预设门限的傅里叶变换后复散射值的函数值点。
所述预设门限为沿序列方向的傅里叶变换后的复散射值的函数的幅度的均值的预先设定倍数。
另一方面，本实施例还提供一种调频连续波信号体制穿墙雷达目标检测快速系统，包括：
至少一个处理器；以及
与所述处理器通信连接的至少一个存储器，其中，所述存储器存储有可被所述处理器执行的程序指令，所述处理器调用所述程序指令能够执行如上述任一所述的调频连续波信号体制穿墙雷达目标检测快速方法。
（三）有益效果
本发明的有益效果是：本发明的一种调频连续波信号体制穿墙雷达目标检测快速方法及系统，由于在差分距离像的序列中，筛选出满足第一预先设定条件的差分距离像；针对满足第一预先设定条件的差分距离像对应距离处的复散射值沿序列方向进行傅里叶变换，得到傅里叶变换后的复散射值的函数；针对傅里叶变换后复散射值的函数沿序列方向求取满足第二预设条件的傅里叶变换后复散射值的函数值点，并将满足第二预设条件的函数值点所对应的距离处作为有效目标。在此过程中进行了两次傅里叶变换，对目标的检测更为准确。
附图说明
图1为本发明的一种调频连续波信号体制穿墙雷达目标检测快速方法流程图；
图2为本发明实施例中的H条距离维响应非相参叠加结果示意图；
图3为本发明实施例中的图2中峰值点位置为13米处结果的包络示意图。
具体实施方式
为了更好的解释本发明，以便于理解，下面结合附图，通过具体实施方式，对本发明作详细描述。
为了更好的理解上述技术方案，下面将参照附图更详细地描述本发明的示例性实施例。虽然附图中显示了本发明的示例性实施例，然而应当理解，可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了能够更清楚、透彻地理解本发明，并且能够将本发明的范围完整的传达给本领域的技术人员。
参见图1，本实施例提供一种调频连续波信号体制穿墙雷达目标检测快速方法，包括：
S1、采用零中频去调频方式接收FMCW接收机按照设定频率所发射的每帧发射信息，并获得与每帧发射信息分别一一对应的差拍信号；
在本实施例的实际应用中，S1具体包括：
将FMCW接收机按照设定频率所发射的任一帧发射信息
Figure 529011DEST_PATH_IMAGE001
的耦合放大一路 作为本振信号，并将该本振信号与FMCW接收机接收到的与所述发射信息对应的回波信号
Figure 140121DEST_PATH_IMAGE002
进行混频和低通滤波，获得与该发射信息
Figure 556059DEST_PATH_IMAGE001
所对应的零中频的差拍信号
Figure 885409DEST_PATH_IMAGE003
。
其中，
Figure 648091DEST_PATH_IMAGE001
为FMCW接收机在t时刻所发射的第m个发射信息。
Figure 644866DEST_PATH_IMAGE004
。
其中，FMCW接收机的发射信息幅度归一化为1。
exp是指数。
j是虚部符号。
f 0 是线性调频起始频率。
T为扫描时间长度。
K为调频斜率，K=B/T。
B为调频信号带宽。
t为时刻。
m为发射信息的序号。
Figure 966126DEST_PATH_IMAGE002
为FMCW接收机所接收到的发射信息为
Figure 284237DEST_PATH_IMAGE001
的回波信号。
其中，
Figure 880303DEST_PATH_IMAGE005
。
A为回波强度。
Figure 731584DEST_PATH_IMAGE006
为发射信息为
Figure 459631DEST_PATH_IMAGE001
的回波信号的延迟。
在本实施例的具体应用中，
Figure 763574DEST_PATH_IMAGE030
。
其中，
Figure 366593DEST_PATH_IMAGE031
为目标起始距离，对应回波延迟为
Figure 806802DEST_PATH_IMAGE032
。
Figure 236909DEST_PATH_IMAGE033
实际表示的是由目标运动速度v引起的多普勒频移
Figure 762568DEST_PATH_IMAGE034
（考 虑到一般系统相对带宽不大，以
Figure 169278DEST_PATH_IMAGE035
近似整个频率范围）与
Figure 729573DEST_PATH_IMAGE035
之间的比值，即
Figure 53283DEST_PATH_IMAGE036
。
与所述发射信息为
Figure 331817DEST_PATH_IMAGE001
对应的差拍信号为
Figure 807798DEST_PATH_IMAGE003
。
其中，
Figure 724064DEST_PATH_IMAGE007
。
其中，
Figure 463349DEST_PATH_IMAGE008
。
Figure 229180DEST_PATH_IMAGE009
；
Figure 977693DEST_PATH_IMAGE010
为目标起始距离；c表示光速。
Figure 748465DEST_PATH_IMAGE011
。
其中，
Figure 455390DEST_PATH_IMAGE012
。
Figure 442938DEST_PATH_IMAGE013
。
Figure 496607DEST_PATH_IMAGE014
为多普勒频率。
在具体应用中，差拍信号是一个线性调频信号，起始频率为：
Figure 886000DEST_PATH_IMAGE037
。
M为接收回波信号的数量，积累总时间为MT。
具体应用中，假设人行走的最大速度为7.2km/h或2m/s，雷达工作于L波段，
Figure 232668DEST_PATH_IMAGE038
，此时行走速度引起的多普勒频移为
Figure 707511DEST_PATH_IMAGE039
，远远低于
Figure 299292DEST_PATH_IMAGE040
，即
Figure 543191DEST_PATH_IMAGE041
，同时假设带宽
Figure 60760DEST_PATH_IMAGE042
和扫描时间
Figure 554059DEST_PATH_IMAGE043
，调频斜率
Figure 683951DEST_PATH_IMAGE044
，因此DK项不可忽略，而
Figure 579095DEST_PATH_IMAGE045
项可忽略，且
Figure 1986DEST_PATH_IMAGE046
。此时各相位 项可化简为：
Figure 484045DEST_PATH_IMAGE047
。
Figure 447322DEST_PATH_IMAGE048
。
S2、根据获得的与每帧发射信息分别一一对应的差拍信号，获取每一差拍信号分别所对应的回波距离像。
在本实施例的实际应用中，所述S2具体包括：
所述差拍信号为
Figure 931393DEST_PATH_IMAGE003
所对应的回波距离像为
Figure 88967DEST_PATH_IMAGE015
。
Figure 25699DEST_PATH_IMAGE016
。
其中，
Figure 995929DEST_PATH_IMAGE017
。
f是t做傅里叶变换后的频率变量表示。
Figure 113269DEST_PATH_IMAGE018
。
S3、针对任意相邻两帧回波距离像，进行差分获得差分距离像的序列；
所述差分距离像的序列包括：依次相邻两帧回波距离像之间的差分值；
在本实施例的实际应用中，所述S3中，所述差分距离像的序列包括：依次相邻两帧回波距离像之间的差分值。
其中，差分距离像的序列是：
Figure 409121DEST_PATH_IMAGE019
。
Figure 600193DEST_PATH_IMAGE020
是第m对相邻两帧回波距离像之间的差分值。
N是预先设定值。
其中，
Figure 374114DEST_PATH_IMAGE021
。
S4、在差分距离像的序列中，筛选出满足第一预先设定条件的差分距离像。
所述第一预先设定条件为：差分距离像大于Tp。
其中，
Figure 832777DEST_PATH_IMAGE022
；
p为预先设定的正实数。
参见图2，具体应用中，由于距离像目标的信杂比相对较低，因此直接通过单帧进行判断可能会造成误判，因此可以将H条（H是某个正整数，对H条回波进行非相参叠加。H可以等于M）回波做傅里叶变换后的回波距离像非相参叠加（如图2），并选择其峰值点（如图2中的圆圈处）。实际为了减少虚假峰值点，可选择某一门限（如均值的若干倍），只选择高于该门限的峰值点。图2中峰值点位置分别为13，15.5，18.5， 21.5，26.5，28米处。
S5、针对满足第一预先设定条件的差分距离像对应距离处的复散射值沿序列方向进行傅里叶变换，得到傅里叶变换后的复散射值的函数。
所述S5具体包括：
S51、基于满足第一预先设定条件的差分距离像，获取满足第一预先设定条件的差分距离像对应距离处的差拍信号。
Figure 33952DEST_PATH_IMAGE023
；
Figure 181161DEST_PATH_IMAGE024
指的是满足第一预先设定条件的差分距离像对应距离处的差拍信号中的 与发射信息的序号m对应的差拍信号。
M指的是满足第一预先设定条件的差分距离像的数量；
S52、对满足第一预先设定条件的差分距离像对应距离处的差拍信号沿序列方向进行傅里叶变换可得第一公式：
Figure 289932DEST_PATH_IMAGE025
；
Figure 603101DEST_PATH_IMAGE026
具体表示指的是满足第一预先设定条件的差分距离像对应距离处的差拍 信号中的与发射信息的序号m对应的差拍信号沿序列方向进行傅里叶变换的结果。
Figure 476642DEST_PATH_IMAGE027
是序号m发射信息的对应频率域变量表示。
S53、通过补偿函数抑制掉第一公式中
Figure 875262DEST_PATH_IMAGE028
，得到沿序列方向的傅里叶变换后 的复散射值的函数；
其中，傅里叶变换后的复散射值的函数为：
Figure 990986DEST_PATH_IMAGE029
。
S6、针对傅里叶变换后复散射值的函数沿序列方向求取满足第二预设条件的傅里叶变换后复散射值的函数值点，并将满足第二预设条件的函数值点所对应的距离处作为有效目标。
所述第二预设条件为：傅里叶变换后复散射值的函数沿序列方向中满足高于预设门限的傅里叶变换后复散射值的函数值点。
所述预设门限为沿序列方向的傅里叶变换后的复散射值的函数的幅度的均值的预先设定倍数。
参见图3，针对上例中距离为13，15.5，18.5，21.5，26.5，28米处峰值点，沿序列方向做傅里叶变换后截取人体行走速度范围内的包络，其中距离为13米处结果的包络如图3。其速度峰值为0.5m/s，该距离上速度峰值点为所有距离处的速度峰值点的最大值，因此可判断该距离处为有效目标。
本实施例中的一种调频连续波信号体制穿墙雷达目标检测快速方法及系统，由于在差分距离像的序列中，筛选出满足第一预先设定条件的差分距离像；针对满足第一预先设定条件的差分距离像对应距离处的复散射值沿序列方向进行傅里叶变换，得到傅里叶变换后的复散射值的函数；针对傅里叶变换后复散射值的函数沿序列方向求取满足第二预设条件的傅里叶变换后复散射值的函数值点，并将满足第二预设条件的函数值点所对应的距离处作为有效目标。在此过程中进行了两次傅里叶变换，对目标的检测更为准确。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例，或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质（包括但不限于磁盘存储器、CD-ROM、光学存储器等）上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备（系统）和计算机程序产品的流程图和／或方框图来描述的。应理解可由计算机程序指令实现流程图和／或方框图中的每一流程和／或方框、以及流程图和／或方框图中的流程和／或方框的结合。
应当注意的是，在权利要求中，不应将位于括号之间的任何附图标记理解成对权利要求的限制。词语“包含”不排除存在未列在权利要求中的部件或步骤。位于部件之前的词语“一”或“一个”不排除存在多个这样的部件。本发明可以借助于包括有若干不同部件的硬件以及借助于适当编程的计算机来实现。在列举了若干装置的权利要求中，这些装置中的若干个可以是通过同一个硬件来具体体现。词语第一、第二、第三等的使用，仅是为了表述方便，而不表示任何顺序。可将这些词语理解为部件名称的一部分。
此外，需要说明的是，在本说明书的描述中，术语“一个实施例”、“一些实施例”、“实施例”、“示例”、“具体示例”或“一些示例”等的描述，是指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
尽管已描述了本发明的优选实施例，但本领域的技术人员在得知了基本创造性概念后，则可对这些实施例作出另外的变更和修改。所以，权利要求应该解释为包括优选实施例以及落入本发明范围的所有变更和修改。
显然，本领域的技术人员可以对本发明进行各种修改和变型而不脱离本发明的精神和范围。这样，倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内，则本发明也应该包含这些修改和变型在内。

Claims (10)

1.一种调频连续波信号体制穿墙雷达目标检测快速方法，其特征在于，包括：
S1、采用零中频去调频方式接收FMCW接收机按照设定频率所发射的每帧发射信息，并获得与每帧发射信息分别一一对应的差拍信号；
S2、根据获得的与每帧发射信息分别一一对应的差拍信号，获取每一差拍信号分别所对应的回波距离像；
S3、针对任意相邻两帧回波距离像，进行差分获得差分距离像的序列；
所述差分距离像的序列包括：依次相邻两帧回波距离像之间的差分值；
S4、在差分距离像的序列中，筛选出满足第一预先设定条件的差分距离像；
S5、针对满足第一预先设定条件的差分距离像对应距离处的差拍信号沿序列方向进行傅里叶变换，得到傅里叶变换后的复散射值的函数；
S6、针对傅里叶变换后复散射值的函数沿序列方向求取满足第二预设条件的傅里叶变换后复散射值的函数值点，并将满足第二预设条件的函数值点所对应的距离处作为有效目标。
2.根据权利要求1所述的方法，其特征在于，S1具体包括：
将FMCW接收机按照设定频率所发射的任一帧发射信息
Figure 81352DEST_PATH_IMAGE001
的耦合放大一路作为本 振信号，并将该本振信号与FMCW接收机接收到的与所述发射信息对应的回波信号
Figure 901540DEST_PATH_IMAGE002
进行混频和低通滤波，获得与该发射信息
Figure 693172DEST_PATH_IMAGE001
所对应的零中频的差拍信号
Figure 573403DEST_PATH_IMAGE003
。
3.根据权利要求2所述的方法，其特征在于，
其中，
Figure 949021DEST_PATH_IMAGE001
为FMCW接收机在t时刻所发射的第m个发射信息；
Figure 307321DEST_PATH_IMAGE004
；
其中，FMCW接收机的发射信息幅度归一化为1；
exp是指数；
j是虚部符号；
f 0 是线性调频起始频率；
T为扫描时间长度；
K为调频斜率，K=B/T；
B为调频信号带宽；
t为时刻；
m为发射信息的序号；
Figure 950530DEST_PATH_IMAGE002
为FMCW接收机所接收到的发射信息为
Figure 1663DEST_PATH_IMAGE001
的回波信号；
其中，
Figure 598997DEST_PATH_IMAGE005
；
Figure 620043DEST_PATH_IMAGE006
为发射信息为
Figure 619223DEST_PATH_IMAGE001
的回波信号的延迟；
A为回波强度。
4.根据权利要求3所述的方法，其特征在于，
与所述发射信息为
Figure 342722DEST_PATH_IMAGE001
对应的差拍信号为
Figure 692931DEST_PATH_IMAGE003
；
其中，
Figure 658613DEST_PATH_IMAGE007
；
其中，
Figure 777879DEST_PATH_IMAGE008
；
Figure 403770DEST_PATH_IMAGE009
；
Figure 241276DEST_PATH_IMAGE010
为目标起始距离；c表示光速；
Figure 886015DEST_PATH_IMAGE011
；
其中，
Figure 984421DEST_PATH_IMAGE012
；
Figure 913370DEST_PATH_IMAGE013
；
Figure 97227DEST_PATH_IMAGE014
为多普勒频率。
5.根据权利要求4所述的方法，其特征在于，所述S2具体包括：
所述差拍信号为
Figure 139132DEST_PATH_IMAGE003
所对应的回波距离像为
Figure 967411DEST_PATH_IMAGE015
；
Figure 935105DEST_PATH_IMAGE016
；
其中，
Figure 481624DEST_PATH_IMAGE017
；
f是t做傅里叶变换后的频率变量表示；
Figure 592799DEST_PATH_IMAGE018
。
6.根据权利要求5所述的方法，其特征在于，所述S3中，
所述差分距离像的序列包括：依次相邻两帧回波距离像之间的差分值；
其中，差分距离像的序列是：
Figure 275584DEST_PATH_IMAGE019
；
Figure 40278DEST_PATH_IMAGE020
是第m对相邻两帧回波距离像之间的差分值；
N是预先设定值；
其中，
Figure 575558DEST_PATH_IMAGE021
。
7.根据权利要求6所述的方法，其特征在于，所述S4中的所述第一预先设定条件为：差分距离像大于Tp；
其中，
Figure 224845DEST_PATH_IMAGE022
；
p为预先设定的正实数。
8.根据权利要求7所述的方法，其特征在于，所述S5具体包括：
S51、基于满足第一预先设定条件的差分距离像，获取满足第一预先设定条件的差分距离像对应距离处的差拍信号；
Figure 762137DEST_PATH_IMAGE023
；
Figure 838677DEST_PATH_IMAGE024
指的是满足第一预先设定条件的差分距离像对应距离处的差拍信号中的与发 射信息的序号m对应的差拍信号；
M指的是满足第一预先设定条件的差分距离像的数量；
S52、对满足第一预先设定条件的差分距离像对应距离处的差拍信号沿序列方向进行傅里叶变换可得第一公式：
Figure 858324DEST_PATH_IMAGE025
；
Figure 45722DEST_PATH_IMAGE026
具体表示指的是满足第一预先设定条件的差分距离像对应距离处的差拍信号 中的与发射信息的序号m对应的差拍信号沿序列方向进行傅里叶变换的结果；
Figure 703100DEST_PATH_IMAGE027
是序号m发射信息的对应频率域变量表示；
S53、通过补偿函数抑制掉第一公式中
Figure 684962DEST_PATH_IMAGE028
，得到沿序列方向的傅里叶变换后的复 散射值的函数；
其中，傅里叶变换后的复散射值的函数为：
Figure 194835DEST_PATH_IMAGE029
。
9.根据权利要求8所述的方法，其特征在于，
所述第二预设条件为：傅里叶变换后复散射值的函数沿序列方向中满足高于预设门限的傅里叶变换后复散射值的函数值点；
所述预设门限为沿序列方向的傅里叶变换后的复散射值的函数的幅度的均值的预先设定倍数。
10.一种调频连续波信号体制穿墙雷达目标检测快速系统，其特征在于，包括：
至少一个处理器；以及
与所述处理器通信连接的至少一个存储器，其中，所述存储器存储有可被所述处理器执行的程序指令，所述处理器调用所述程序指令能够执行如权利要求1-9中任一所述的调频连续波信号体制穿墙雷达目标检测快速方法。
CN202211091992.4A 2022-09-07 2022-09-07 调频连续波信号体制穿墙雷达目标检测快速方法及系统 Active CN115166681B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211091992.4A CN115166681B (zh) 2022-09-07 2022-09-07 调频连续波信号体制穿墙雷达目标检测快速方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211091992.4A CN115166681B (zh) 2022-09-07 2022-09-07 调频连续波信号体制穿墙雷达目标检测快速方法及系统

Publications (2)

Publication Number Publication Date
CN115166681A CN115166681A (zh) 2022-10-11
CN115166681B true CN115166681B (zh) 2022-12-06

Family

ID=83481897

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211091992.4A Active CN115166681B (zh) 2022-09-07 2022-09-07 调频连续波信号体制穿墙雷达目标检测快速方法及系统

Country Status (1)

Country Link
CN (1) CN115166681B (zh)

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10451723B2 (en) * 2016-12-20 2019-10-22 National Chung-Shan Institute Of Science & Technology Signal processing apparatus of a continuous-wave (CW) radar sensing system
CN107153189B (zh) * 2017-04-18 2021-08-03 上海交通大学 线性调频连续波雷达测距的信号处理方法
CN107765240B (zh) * 2017-10-20 2020-11-06 中国人民解放军海军七〇一工厂 一种运动状态的判断方法、装置和电子设备
CN108387894B (zh) * 2018-04-13 2021-07-27 中南大学 穿墙雷达回波数据的处理方法
CN108776336B (zh) * 2018-06-11 2022-06-03 电子科技大学 一种基于emd的自适应穿墙雷达静止人体目标定位方法
CN112654888A (zh) * 2018-09-05 2021-04-13 京瓷株式会社 电子设备、电子设备的控制方法、以及电子设备的控制程序
CN110873877B (zh) * 2019-04-25 2021-04-23 北京航空航天大学 目标运动轨迹的确定方法及装置
CN110118966B (zh) * 2019-05-28 2020-10-13 长沙莫之比智能科技有限公司 基于毫米波雷达的人员检测与计数系统
CN110244296B (zh) * 2019-06-17 2020-07-14 中国人民解放军国防科技大学 基于距离像复包络相位匹配处理的空间进动目标成像方法
CN110161491B (zh) * 2019-06-28 2021-01-12 电子科技大学 一种针对微弱生命体的测距和呼吸频率估计方法
CN111352102B (zh) * 2020-03-18 2023-02-24 重庆邮电大学 一种基于调频连续波雷达的多目标个数检测方法及装置

Also Published As

Publication number Publication date
CN115166681A (zh) 2022-10-11

Similar Documents

Publication Publication Date Title
CN107966688B (zh) 基于相位干涉技术的宽带雷达目标速度解模糊方法
CN110646773B (zh) 一种基于汽车毫米波雷达的弱目标检测方法、跟踪方法及系统
Li et al. Fast coherent integration for maneuvering target with high-order range migration via TRT-SKT-LVD
US8305257B2 (en) Method and apparatus for coherent marine radar measurements of properties of ocean waves and currents
US10054668B2 (en) Probabilistic signal, detection, and track processing architecture and system
CN107450055B (zh) 基于离散线性调频傅立叶变换的高速机动目标检测方法
CN107656255A (zh) 基于多径回波的超宽带雷达动目标二维定位方法
US20220179062A1 (en) Detection apparatus and method
CN108535730A (zh) 一种多普勒气象雷达解速度模糊方法和系统
CN110376563B (zh) 基于gps辐射源的雷达运动目标检测方法
CN109298417B (zh) 一种基于雷达信号处理的建筑物内部结构探测方法及装置
Heaney et al. Measurements of three-dimensional propagation in a continental shelf environment
CN105589061A (zh) 一种岸基雷达的信号处理算法
US7942821B2 (en) Doppler velocity detection device and ultrasonographic device using the same
CN115166681B (zh) 调频连续波信号体制穿墙雷达目标检测快速方法及系统
CN113419238A (zh) 基于毫米波雷达的山体滑坡监测方法、电子设备、存储介质
JP5423713B2 (ja) アクティブソーナー装置及びその信号処理方法
Kronauge Waveform design for continuous wave radars
US6721678B1 (en) Determining three dimensional velocity of a object
Ranney et al. Recent MTI experiments using ARL's synchronous impulse reconstruction (SIRE) radar
CN110907930B (zh) 一种基于角度估计的车载雷达目标检测估计方法及装置
CN116027288A (zh) 生成数据的方法、装置、电子设备及存储介质
CN108387891B (zh) 基于短基线一发两收的超宽带雷达干涉二维定位跟踪方法
CN108983192B (zh) 基于gps辐射源的雷达运动目标参数估计方法
CN112698291B (zh) 一种基于cnn的机载气象雷达气象目标检测方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant