一种穿墙雷达对多运动目标实时探测、分离的方法
【技术领域】
[0001] 本发明属于穿墙雷达技术领域，特别是涉及一种穿墙雷达对多运动目标实时探 测、分离的方法。
【背景技术】
[0002] 穿墙雷达自一面世就引起了社会各方的广泛关注，英国《新科学家》杂志曾经在 2009年将穿墙雷达列为未来30年最有前景的十项科技发明之一，并称它的作用将如同手 机、互联网等一样对现代人类生活产生巨大的影响。穿墙雷达作为一种正在发展中的新技 术，该技术可以广泛地应用在军事、反恐、救灾、医疗、警戒、安检等领域。
[0003] 对于穿墙雷达技术而言，国内外学者的研究主要集中在穿墙雷达的硬件设计以及 其信号处理技术两个方面。而目标探测、虚假目标剔除、成像及多运动目标实时搜索、跟踪 等是穿墙雷达的信号处理技术主要的研究方向。其中，对运动目标实时探测是穿墙雷达的 关键技术之一，它决定了穿墙雷达对各运动中的目标位置准确探测的策略，是目标的定位 和实时跟踪的基石，并直接影响了穿墙雷达的整体性能。目前，关于国内外的专家和学者对 单运动目标的实时探测取得了极大的成功，推出了各种有效地实时单运动目标探测、定位 技术。
[0004] 但是在多运动目标场景下，由于穿墙雷达的回波信号成分复杂，包含发射天线和 接收天线之间的直接耦合波、墙体界面的反射波、系统热噪声和周围环境的白噪声、射频干 扰信号和墙壁后面目标的反射回波等，特别是由于多运动目标场景下各目标的回波信号经 过多次反射后相互干扰，导致回波信号成分更加复杂，真实目标被淹没在其中从而难以探 测出真实目标，故在多运动目标自动探测、定位方面的进展却不尽人意。
[0005] 如果有一种新的方法能够对穿墙雷达的回波信号按照各运动目标到接收天线的 距离进行有效分离，即每区域里只含一个运动目标，那么再结合现有的单运动目标的探测 方法就能够实现对多运动目标的实时探测。因此，如何在多运动目标场景下实现对穿墙雷 达回波信号按照各运动目标实时分离的技术，成为本领域亟待解决的技术问题。
【发明内容】
[0006] 针对上述问题，本发明提出一种穿墙雷达对多运动目标实时探测、分离的方法。
[0007] 本发明提供的技术方案是：一种穿墙雷达对多运动目标实时探测、分离的方法，包 括如下步骤：
[0008] 步骤1、预采集前K次的一维回波信号，求出它们的平均值y?，η代表离散的时 间，且n e [1，Ν]，其中N为每条扫描数据的采用点数；
[0009] 步骤2、从第Κ+1条扫描数据开始，将当前采集的扫描数据，假设为第k条扫描数 据，经过指数平滑滤波处理，获得消除直达波及静止目标的信号X' k (η);
[0010] 步骤3、对当前消除直达波及静止目标的信号x'k(n)进行希尔伯特变换得到 之(《)，进一步获得当前消除直达波及静止目标的信号x' k(n)的包络E(n);
[0011] 步骤4、利用固定阈值1^对当前消除直达波及静止目标的信号x' k(n)的做粗分 离，其实现步骤为，
[0012] 步骤4. 1，计算包络E (η)所有点的幅度的平均僅
[0013] 步骤4. 2,令阈值
，将该阈值与包络E (η)幅度值进行比较，
>当Xciutl (η)为1就表明该位置有目标，Xciutl (η)为〇就表明该位 置没有目标，从而将包络E (η)的幅度值大于阈值TH1的点形成非零区域；
[0014] 步骤4. 3，令X' k (η) = X' k (η) X Xciutl (η)，将当前消除直达波及静止目标的信号 X' k (η)的幅值不等于0的点构成若干有效区域；
[0015] 步骤5、利用自适应阈值TH2对步骤4中的某一个有效区域做细分离，假设有效区 域为闭区间[a，b]，其中[a，b] ￡ [1，Ν]，所述的步骤5包括如下步骤：
[0016] 步骤5. 1，假设滑动窗的长度为2M-1 ;
[0017] 步骤5. 2,计算该闭区间第一个点，即a点，对应的自适应阈值TH2，其中
[0018] 步骤5. 3,如果该点对应的包络幅度值E(a)彡C2XTH2,则令Χ_2( &) = 1，否则 Xmn2 (a) = 〇, Xciut2 (a)为1就表明该位置有目标，Xciut2 (a)为0就表明该位置没有目标；
[0019] 步骤5. 4,重复步骤5. 2~5. 3所述方法依次求出闭区间[a，b]其它点的输出结 果Xciut2 (m)，再利用X' k (m) = X' k (m) X Xciut2 (m)，则对经过步骤4的粗分离后形成的有效区域 [a，b]再进一步分离为若干不同的有效区域；
[0020] 步骤6、对步骤4中形成的其他有效区域重复步骤5. 1-5. 4的方法分别完成对当前 消除直达波及静止目标的信号X' k (η)在其他区域的细分离。
[0021] 优选的，所述的步骤2包括如下步骤：
[0022] 步骤2. 1，利用前一条扫描数据获得的平滑序列yk Jn)和当前采集的扫描数据 Xk (η)计算出平滑参数α，其中\
.，min为求最小值运算，max为求最大 值运算；
[0023] 步骤2. 2,利用yk Jn)、xk(n)和α计算出当前的平滑序列yk(n)，其中yk(n)= α Yk ι(η) + (1-α )xk(n)；
[0024] 步骤2. 3,将当前采集的扫描数据xk(η)减去前一条扫描数据获得的平滑序列 yk I (η)，获得消除直达波及静止目标的信号X' k (η)，所述的X' k (n) = Xk (n) -yk i (η)。
[0025] 优选的，所述的步骤3包括如下步骤：
[0026] 步骤3. 1，对当前消除直达波及静止目标的信号x'k(n)做离散时间傅立叶变换得 到频谱Χ(ω); CN 105137423 A ~P 3/6 页
[0027] 步骤3· 2,将Χ(ω)的相位φ (ω)按照公式
做 90°或-90°移相得到相位Φ(ω)，如果ω>〇,做90°移相得到相位Φ(ω)，如果ω〈〇， 做-90°移相得到相位Φ (ω);再由Χ(ω)的模|Χ(ω) I和相位Φ (ω)构成新的频谱 Υ(ω) = |Χ(ω) |eJ〇(M);
[0028] 步骤3.3,对Υ(ω)做傅立叶逆变换后再乘以因子i得到对当前消除直达波及静 π 止目标的信号x'k(n)进行希尔伯特变换后的结果又⑶，其中η为圆周率；
[0029] 步骤3. 4,利用当前消除直达波及静止目标的信号X' k (η))和:? (?)构建 解析信号= + ⑷，j为虚部符号，求出解析信号Z( n)的瞬时幅值
即获得当前消除直达波及静止目标的信号X' k (η)的包络E (η)。
[0030] 优选的，所述的步骤4. 2中参数(^取0· 2。
[0031] 优选的，所述的步骤5. 3中其中(：2取0. 8。
[0032] 优选的，所述的步骤5中M取15。
[0033] 优选的，所述的步骤1和步骤2中K取5。
[0034] 本发明的有益效果是：一种穿墙雷达对多运动目标实时探测、分离的方法，通过一 种基于指数平滑滤波器的杂波抑制算法剔除复杂的背景信号，更快速、准确地保留运动人 体目标的信息；通过一种基于希尔伯特变换的包络检测算法，从而实现将各运动目标置于 不同包络的波峰；通过一种基于将固定阈值和自适应阈值相结合的算法，从而更快速地实 现将各运动目标的分离在不同的区域；本发明方法简单、执行效率高，适用于便携式穿墙雷 达隔墙对多运动人体目标的探测，为穿墙雷达对各运动目标准确定位提供了有力的保障， 可有效提尚穿墙雷达的工作效率。
【附图说明】
[0035] 图1为穿墙雷达回波示意图；
[0036] 图2为目标多次反射示意图；
[0037] 图3为对X ( ω )做的移相的程序框图；
【具体实施方式】
[0038] 1、理论基础
[0039] 如附图1所示，穿墙雷达的回波信号可表示为：
[0040] fr (t) = fd (t) +sr (t) +fc (t) +fn (t) +fw (t) (I)
[0041] 其中f；(t)是墙雷达的回波信号，sjt)是墙体后的运动目标信号，fd(t)是直接耦 合波，f w(t)是墙体的反射波，f Jt)是背景杂波信号，fn(t)是噪声，t表示连续的时间。其 中，墙体后的运动目标信号sjt)除了包含了目标的真实回波外还包含了目标的多次反射 回波，如附图2所示，图中别为对应的反射波达到接收天线的时间。
[0042] 下面从理论上讨论目标的多次反射回波对目标探测的影响。为了分析的简单，附 图2假设是在理想情况下的示意图，即忽略导电性（〇~0)且所有的物体都是非磁性的 (导磁系数μ =1)，那么折射率π *士，其中叫和η2分别代表空气和墙体的折射率，ε是 介电常数，反射回来的部分入射能量P为：
[0043]
樹
[0044] 按照能量守恒，透射系数τ = Ι-p，那么将反射回来的部分入射能量ρ和透射系 数τ分别用dB表示为L1^P L τ，则
[0045]
[0046] (4)
[0047] 考虑信号单程传输时能量的损失情况。假设入射信号的能量为P1，第一次从空气 透射进入墙体的能量P m= τ P i (W)，从墙体前侧反射回的能量Ph= P P i，从墙体后侧反射 回的能量Prt= P (TP1),则入射信号穿透墙体后的能量Pni= τ 2P1。假设入射信号的能量 P1Sl瓦特，即30毫分贝（dBm)，空气的介电常数ε 1= 1，墙体的介电常数ε 2= 10,那么 P = 0· 269，τ = (Ι-p ) = 〇· 73,则第一次从空气透射进入墙体的能量Pni= 0· 532W，即 27. 26dBm，故在信号单程传输时信号能量的损失为2. 7dBm。因此，通过上面的分析可得下面 结论：（一）、目标真实回波信号的幅度比多次反射后的回波信号幅度强；（二）、尽管多次 反射后的回波信号幅度比目标真实回波信号的幅度弱，但由于信号传输时信号能量损失较 小，所以在多目标场景下很难将多次反射后的回波信号和其它真实目标的回波信号分辨出 来。那么，如果能够按照不同的目标将每一次采集的回波信号进行分段，再利用结论一就能 够实现对对多运动人体目标的探测。2、对回波信号按照各运动目标的有效分离
[0048] 由于穿墙探测雷达回波信号中目标真实回波信号的幅度比多次反射后的回波信 号幅度强，则真实目标的位置就是回波信号的包络中的极大值，那么检测出回波信号包络 的波峰就可以将每个目标分隔在不同的区域。
[0049] 如果频谱Χ(ω) =XRf；(co)+jXIm(co) = |Χ(ω) |θ·]4>(ω)是去除直达波及静止目标后 的回波信号X'（t)的傅立叶变换，其中ω表示角频率，j为虚部符号，X rJco)表示Χ(ω) 的实部，ΧΙηι(ω)表示Χ(ω)的虚部，|Χ(ω) I表示Χ(ω)的模，φ (ω)表示Χ(ω)的相位， 则它的希尔伯特变换i⑴可以表示为：
[0050]
(5)
[0051] 将去除直达波及静止目标后的信号X'（t)和其希尔伯特变换W)构建