发明名称 一种穿墙雷达对多运动目标实时探测、分离的方法
摘要 一种穿墙雷达对多运动目标实时探测、分离的方法,通过一种基于指数平滑滤波器的杂波抑制算法剔除复杂的背景信号,更快速、准确地保留运动人体目标的信息;通过一种基于希尔伯特变换的包络检测算法,从而实现将各运动目标置于不同包络的波峰;通过一种基于将固定阈值和自适应阈值相结合的算法,从而更快速地实现将各运动目标的分离在不同的区域;本发明方法简单、执行效率高,适用于便携式穿墙雷达隔墙对多运动人体目标的探测,为穿墙雷达对各运动目标准确定位提供了有力的保障,可有效提高穿墙雷达的工作效率。<b />
申请公布号 CN105137423A 申请公布日期 2015.12.09
申请号 CN201510642301.9 申请日期 2015.09.30
申请人 武汉大学 发明人 陈小莉;白迪;田茂
分类号 G01S13/88(2006.01)I 主分类号 G01S13/88(2006.01)I
代理机构 武汉科皓知识产权代理事务所(特殊普通合伙) 42222 代理人 张火春
主权项 一种穿墙雷达对多运动目标实时探测、分离的方法,其特征在于,包括如下步骤:步骤1、预采集前K次的一维回波信号,求出它们的平均值y<sub>0</sub>(n),n代表离散的时间,且n∈[1,N],其中N为每条扫描数据的采用点数;步骤2、从第K+1条扫描数据开始,将当前采集的扫描数据,假设为第k条扫描数据,经过指数平滑滤波处理,获得消除直达波及静止目标的信号x'<sub>k</sub>(n);步骤3、对当前消除直达波及静止目标的信号x'<sub>k</sub>(n)进行希尔伯特变换得到<img file="FDA0000815940450000014.GIF" wi="149" he="66" />进一步获得当前消除直达波及静止目标的信号x'<sub>k</sub>(n)的包络E(n);步骤4、利用固定阈值TH<sub>1</sub>对当前消除直达波及静止目标的信号x'<sub>k</sub>(n)的做粗分离,其实现步骤为,步骤4.1,计算包络E(n)所有点的幅度的平均值<img file="FDA0000815940450000011.GIF" wi="341" he="141" />步骤4.2,令阈值<img file="FDA0000815940450000012.GIF" wi="306" he="84" />将该阈值与包络E(n)幅度值进行比较,<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>X</mi><mrow><mi>o</mi><mi>u</mi><mi>t</mi><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mn>1</mn><mo>,</mo></mrow></mtd><mtd><mrow><mi>E</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>&GreaterEqual;</mo><msub><mi>TH</mi><mn>1</mn></msub></mrow></mtd></mtr><mtr><mtd><mrow><mn>0</mn><mo>,</mo></mrow></mtd><mtd><mrow><mi>E</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>&lt;</mo><msub><mi>TH</mi><mn>1</mn></msub></mrow></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000815940450000013.GIF" wi="623" he="161" /></maths>当X<sub>out1</sub>(n)为1就表明该位置有目标,X<sub>out1</sub>(n)为0就表明该位置没有目标,从而将包络E(n)的幅度值大于阈值TH<sub>1</sub>的点形成非零区域;步骤4.3,令x'<sub>k</sub>(n)=x'<sub>k</sub>(n)×X<sub>out1</sub>(n),将当前消除直达波及静止目标的信号x'<sub>k</sub>(n)的幅值不等于0的点构成若干有效区域;步骤5、利用自适应阈值TH<sub>2</sub>对步骤4中的某一个有效区域做细分离,假设有效区域为闭区间[a,b],其中<img file="FDA0000815940450000015.GIF" wi="308" he="65" />所述的步骤5包括如下步骤:步骤5.1,假设滑动窗的长度为2M‑1;步骤5.2,计算该闭区间第一个点,即a点,对应的自适应阈值TH<sub>2</sub>,其中<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>TH</mi><mn>2</mn></msub><mo>=</mo><mfrac><mrow><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mi>a</mi><mo>-</mo><mi>M</mi></mrow><mrow><mi>a</mi><mo>+</mo><mi>M</mi></mrow></munderover><mi>E</mi><mrow><mo>(</mo><mi>j</mi><mo>)</mo></mrow></mrow><mrow><mn>2</mn><mi>M</mi><mo>+</mo><mn>1</mn></mrow></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000815940450000021.GIF" wi="387" he="213" /></maths>步骤5.3,如果该点对应的包络幅度值E(a)≥C<sub>2</sub>×TH<sub>2</sub>,则令X<sub>out2</sub>(a)=1,否则X<sub>out2</sub>(a)=0,X<sub>out2</sub>(a)为1就表明该位置有目标,X<sub>out2</sub>(a)为0就表明该位置没有目标;步骤5.4,重复步骤5.2~5.3所述方法依次求出闭区间[a,b]其它点的输出结果X<sub>out2</sub>(m),再利用x'<sub>k</sub>(m)=x'<sub>k</sub>(m)×X<sub>out2</sub>(m),则对经过步骤4的粗分离后形成的有效区域[a,b]再进一步分离为若干不同的有效区域;步骤6、对步骤4中形成的其他有效区域重复步骤5.1‑5.4的方法分别完成对当前消除直达波及静止目标的信号x'<sub>k</sub>(n)在其他区域的细分离。
地址 430072 湖北省武汉市武昌区珞珈山武汉大学