发明名称 基于认知框架的雷达动目标检测方法
摘要 本发明公开了一种基于认知框架的雷达动目标检测方法,主要解决现有方法在强杂波背景中对低速目标的检测性能差的缺陷。其实现过程是:从回波中分离出地杂波分量;根据地杂波分量得到当前帧的杂波图;根据杂波图感知环境将当前帧划分为噪声区和杂波区;根据地杂波分量计算相关系数,利用它校正天线指向误差,并将校正后的回波数据存入环境动态数据库EDDB中;根据当前帧的划分,利用EDDB中的回波数据估计干扰协方差矩阵;根据估计的干扰协方差矩阵得到自适应滤波器系数,并对输入数据做自适应滤波处理;对滤波器的输出做恒虚警检测。本发明具有对现有装备不做大的改动即可大幅提高其检测性能的优点,可用于强杂波背景中对低速目标的检测。
申请公布号 CN103197297B 申请公布日期 2014.08.06
申请号 CN201310091565.0 申请日期 2013.03.21
申请人 西安电子科技大学 发明人 刘宏伟;方明;戴奉周;赵永波;纠博
分类号 G01S7/41(2006.01)I 主分类号 G01S7/41(2006.01)I
代理机构 陕西电子工业专利中心 61205 代理人 王品华;朱红星
主权项 一种基于认知框架的雷达动目标检测方法,包括如下步骤:1)令x<sub>mip</sub>表示雷达第m帧第i个方位角第p个距离单元的回波数据,其中1≤m≤m<sub>0</sub>,1≤i≤I,1≤p≤D,m<sub>0</sub>表示当前帧序号,I表示每一帧雷达方位角的个数,D表示每一个方位角的距离单元数,h<sub>0</sub>表示零速滤波器,其系数全为1;通过零速滤波器得到回波数据中的地杂波分量y<sub>mip</sub>为:<img file="FDA0000521804520000011.GIF" wi="316" he="83" />其中T表示转置;2)根据地杂波分量y<sub>mip</sub>,获得当前帧m<sub>0</sub>第i个方位角第p个距离单元的干扰电平功率的估计值<img file="FDA0000521804520000012.GIF" wi="152" he="92" />3)根据步骤2)得到的干扰电平功率的估计值,将当前帧m<sub>0</sub>划分为杂波区和噪声区这两个区域;4)计算当前帧m<sub>0</sub>第i个方位角与先前帧m<sub>1</sub>第j个方位角的相关系数<img file="FDA0000521804520000014.GIF" wi="115" he="83" />:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>&rho;</mi><mrow><msub><mi>m</mi><mn>0</mn></msub><msub><mi>m</mi><mn>1</mn></msub><mi>ij</mi></mrow></msub><mo>=</mo><mfrac><mrow><mo>|</mo><mo>|</mo><msup><msub><mi>Y</mi><mrow><msub><mi>m</mi><mn>0</mn></msub><mi>i</mi></mrow></msub><mi>H</mi></msup><mo>&CenterDot;</mo><msub><mi>Y</mi><mrow><msub><mi>m</mi><mn>1</mn></msub><mi>j</mi></mrow></msub><mo>|</mo><mo>|</mo></mrow><mrow><mo>|</mo><mo>|</mo><msup><msub><mi>Y</mi><mrow><msub><mi>m</mi><mn>0</mn></msub><mi>i</mi></mrow></msub><mi>H</mi></msup><mo>|</mo><mo>|</mo><mo>&CenterDot;</mo><mo>|</mo><mo>|</mo><msub><mi>Y</mi><mrow><msub><mi>m</mi><mn>1</mn></msub><mi>j</mi></mrow></msub><mo>|</mo><mo>|</mo></mrow></mfrac><mo>,</mo><mn>1</mn><mo>&le;</mo><msub><mi>m</mi><mn>1</mn></msub><mo>&lt;</mo><msub><mi>m</mi><mn>0</mn></msub><mo>,</mo><mn>1</mn><mo>&le;</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>&le;</mo><mi>I</mi></mrow>]]></math><img file="FDA0000521804520000013.GIF" wi="1023" he="187" /></maths>其中||·||表示2范数,<img file="FDA0000521804520000015.GIF" wi="574" he="111" />表示当前帧m<sub>0</sub>第i个方位角的一段长度为C的地杂波数据,<img file="FDA0000521804520000016.GIF" wi="581" he="109" />表示第m<sub>1</sub>帧第j个方位角的一段长度为C的地杂波数据,H表示共轭转置;5)根据当前帧m<sub>0</sub>第i个方位角与先前帧m<sub>1</sub>第j个方位角的相关系数<img file="FDA0000521804520000017.GIF" wi="110" he="58" />校正天线的指向误差:若<img file="FDA0000521804520000018.GIF" wi="528" he="93" />则判定当前帧m<sub>0</sub>与先前帧m<sub>1</sub>是对准的,并将当前帧m<sub>0</sub>的第i<sub>0</sub>个方位角和先前帧m<sub>1</sub>的第j<sub>0</sub>个方位角分别作为当前帧m<sub>0</sub>和先前帧m<sub>1</sub>的基准,完成指向误差的校正,校正后的回波数据存入环境动态数据库EDDB中;反之,判定当前帧m<sub>0</sub>与先前帧m<sub>1</sub>是未对准的,舍去先前帧m<sub>1</sub>的回波数据;6)根据步骤3)的划分,利用环境动态数据库中的回波数据,计算当前帧m<sub>0</sub>第i个方位角第p个距离单元的干扰协方差矩阵的估计值<img file="FDA0000521804520000019.GIF" wi="101" he="80" />:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mover><mi>R</mi><mo>^</mo></mover><mrow><msub><mi>m</mi><mn>0</mn></msub><mi>ip</mi></mrow></msub><mo>=</mo><mrow><mo>(</mo><mfrac><mn>1</mn><mrow><mn>2</mn><mi>L</mi></mrow></mfrac><mo>&CenterDot;</mo><munderover><mi>&Sigma;</mi><mrow><mi>u</mi><mo>=</mo><mi>p</mi><mo>-</mo><mi>L</mi><mo>,</mo><mi>u</mi><mo>&NotEqual;</mo><mi>p</mi></mrow><mrow><mi>p</mi><mo>+</mo><mi>L</mi></mrow></munderover><msub><mi>t</mi><mrow><msub><mi>m</mi><mn>0</mn></msub><mi>iu</mi></mrow></msub><mo>&CenterDot;</mo><msubsup><mi>t</mi><mrow><msub><mi>m</mi><mn>0</mn></msub><mi>iu</mi></mrow><mi>H</mi></msubsup><mo>+</mo><munder><mi>&Sigma;</mi><mrow><msub><mi>m</mi><mn>1</mn></msub><mo>&lt;</mo><msub><mi>m</mi><mn>0</mn></msub></mrow></munder><msup><mi>&beta;</mi><mrow><msub><mi>m</mi><mn>0</mn></msub><mo>-</mo><msub><mi>m</mi><mn>1</mn></msub></mrow></msup><mo>&CenterDot;</mo><mrow><mo>(</mo><mfrac><mn>1</mn><mrow><mn>2</mn><mi>L</mi><mo>+</mo><mn>1</mn></mrow></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>v</mi><mo>=</mo><mi>p</mi><mo>-</mo><mi>L</mi></mrow><mrow><mi>p</mi><mo>+</mo><mi>L</mi></mrow></munderover><msub><mi>t</mi><mrow><msub><mi>m</mi><mn>1</mn></msub><mi>iv</mi></mrow></msub><mo>&CenterDot;</mo><msubsup><mi>t</mi><mrow><msub><mi>m</mi><mn>1</mn></msub><mi>iv</mi></mrow><mi>H</mi></msubsup><mo>)</mo></mrow><mo>)</mo></mrow><mo>/</mo><mrow><mo>(</mo><mn>1</mn><mo>+</mo><munder><mi>&Sigma;</mi><mrow><msub><mi>m</mi><mn>1</mn></msub><mo>&lt;</mo><msub><mi>m</mi><mn>0</mn></msub></mrow></munder><msup><mi>&beta;</mi><mrow><msub><mi>m</mi><mn>0</mn></msub><mo>-</mo><msub><mi>m</mi><mn>1</mn></msub></mrow></msup><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000521804520000021.GIF" wi="1605" he="174" /></maths>其中,L表示单侧参考窗的长度,β表示遗忘因子,这两个参数取决于估计单元所处的区域;<img file="FDA0000521804520000029.GIF" wi="92" he="70" />表示环境动态数据库中当前帧m<sub>0</sub>第i个方位角第u个距离单元的回波数据;<img file="FDA00005218045200000210.GIF" wi="87" he="70" />表示环境动态数据库中先前帧m<sub>1</sub>第i个方位角第v个距离单元的回波数据;7)利用上述估计值<img file="FDA0000521804520000022.GIF" wi="98" he="84" />在最大信噪比MSN准则下,计算当前帧m<sub>0</sub>第i个方位角第p个距离单元的滤波器系数<img file="FDA00005218045200000211.GIF" wi="118" he="68" />;8)利用步骤7)的滤波器系数<img file="FDA00005218045200000212.GIF" wi="114" he="68" />,得到当前帧m<sub>0</sub>第i个方位角第p个距离单元的输出信号<img file="FDA00005218045200000213.GIF" wi="100" he="53" />;9)设定的目标检测虚警概率为P<sub>fa</sub>,则与其对应的检测阈值为T=‑ln(1‑(1‑P<sub>fa</sub>)<sup>1/N</sup>),将自适应滤波后输出信号的功率<img file="FDA00005218045200000214.GIF" wi="148" he="109" />与检测阈值T和背景功率估计值<img file="FDA0000521804520000023.GIF" wi="355" he="96" />的乘积进行比较,若在输出的N个通道中只要有一个通道满足<img file="FDA0000521804520000024.GIF" wi="642" he="108" />则表示检测到目标;若在输出的N个通道中所有的通道都满足<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msup><mrow><mo>|</mo><mo>|</mo><msub><mi>z</mi><mrow><msub><mi>m</mi><mn>0</mn></msub><mi>ipk</mi></mrow></msub><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msup><mo>&lt;</mo><mi>T</mi><mo>&CenterDot;</mo><msup><msub><mi>w</mi><mrow><msub><mi>m</mi><mn>0</mn></msub><mi>ipk</mi></mrow></msub><mi>H</mi></msup><msub><mover><mi>R</mi><mo>^</mo></mover><mrow><msub><mi>m</mi><mn>0</mn></msub><mi>ip</mi></mrow></msub><msub><mi>w</mi><mrow><msub><mi>m</mi><mn>0</mn></msub><mi>ipk</mi></mrow></msub><mo>,</mo></mrow>]]></math><img file="FDA0000521804520000025.GIF" wi="641" he="108" /></maths>则表示没有检测到目标。
地址 710071 陕西省西安市太白南路2号
您可能感兴趣的专利