发明名称 一种基于回波扩展熵的穿墙雷达成像方法
摘要 本发明提供一种基于回波扩展熵的穿墙雷达成像方法,具体步骤包括:建立穿墙成像模型,获得原始信号数据;对天线阵元所接受的回波信号进行扩展;计算熵值并设定门限;对天线阵元所接收的回波信号反复进行扩展,得到阵元不同时间采样点所对应的扩展熵值,利用该熵值设定门限,获得经过处理后的处理图。本发明能够消除墙体强杂波对穿墙成像的影响,提高成像精度,减少了天线阵元数目,对硬件设备的要求较低,大大简化了穿墙雷达硬件设备的设计。
申请公布号 CN104297749B 申请公布日期 2017.01.04
申请号 CN201410607213.0 申请日期 2014.10.31
申请人 南京信息工程大学 发明人 李家强;蔡洪渊;陈金立;李鹏
分类号 G01S13/89(2006.01)I 主分类号 G01S13/89(2006.01)I
代理机构 南京众联专利代理有限公司 32206 代理人 顾进
主权项 一种基于回波扩展熵的穿墙雷达成像方法,其特征在于:具体步骤包括:(1)建立穿墙成像模型,获得原始信号数据:利用均匀分布于测线的自发自收天线阵列,依次发射并接收信号,并且导出未经处理的成像图;(2)对天线阵元所接收的回波信号进行扩展:收发天线阵元数目为N,对接收天线接收到的原始回波信号做时间采样,采样点数为M,构成M×N的矩阵e,对矩阵e做如下处理:<maths num="0001"><math><![CDATA[<mrow><mi>q</mi><mrow><mo>(</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>e</mi><msup><mrow><mo>(</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mn>2</mn></msup></mrow><mrow><msubsup><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></msubsup><mi>e</mi><msup><mrow><mo>(</mo><mi>m</mi><mo>,</mo><mi>k</mi><mo>)</mo></mrow><mn>2</mn></msup></mrow></mfrac><mo>,</mo><mi>m</mi><mo>=</mo><mn>0</mn><mo>,</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>M</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>n</mi><mo>=</mo><mn>0</mn><mo>,</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>N</mi><mo>-</mo><mn>1</mn><mo>,</mo></mrow>]]></math><img file="FDA0001074259430000011.GIF" wi="1334" he="151" /></maths>其中<img file="FDA0001074259430000012.GIF" wi="299" he="126" />再根据离散无记忆的扩展信源熵值理论,对天线阵元所接收的回波信号做扩展处理,扩展次数p=2,3,...,P,P∈N<sup>*</sup>,N<sup>*</sup>为正整数,扩展后构成M×N<sup>p</sup>的矩阵E,并对矩阵E进行处理,处理公式为:<maths num="0002"><math><![CDATA[<mrow><mi>Q</mi><mrow><mo>(</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>E</mi><msup><mrow><mo>(</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mn>2</mn></msup></mrow><mrow><msubsup><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>0</mn></mrow><mrow><msup><mi>N</mi><mi>p</mi></msup><mo>-</mo><mn>1</mn></mrow></msubsup><mi>E</mi><msup><mrow><mo>(</mo><mi>m</mi><mo>,</mo><mi>k</mi><mo>)</mo></mrow><mn>2</mn></msup></mrow></mfrac><mo>,</mo><mi>m</mi><mo>=</mo><mn>0</mn><mo>,</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>M</mi><mo>-</mo><mn>1</mn><mo>;</mo></mrow>]]></math><img file="FDA0001074259430000013.GIF" wi="958" he="158" /></maths>(3)计算熵值并设定门限:利用如下公式计算熵值:<img file="FDA0001074259430000014.GIF" wi="708" he="135" />并根据最大熵值理论,设置门限;设置门限后对矩阵e继续处理,获得经过处理后的处理图,并且设定门限的公式为:<img file="FDA0001074259430000015.GIF" wi="686" he="151" />其中α为控制门限的参数;所述对矩阵e的处理公式为:e<sub>Z</sub>(m,n)=Z(m)e(m,n);所述e<sub>Z</sub>为矩阵e经过处理后的最终结果。
地址 210044 江苏省南京市宁六路219号