发明名称 一种适用于高分辨阵列的快速鲁棒的信源个数检测方法
摘要 本发明公开了一种适用于高分辨阵列的快速鲁棒的信源个数检测方法,属于信号处理领域。建立阵列天线的接收数据模型,用多级维纳滤波器对接收数据进行逐级滤波,得到各级的最小均方误差,然后用最小均方误差计算最大似然函数。同时,用计算最大似然函数所涉及的自由参数来确定模型中的未知自由参数的个数,最后得到一个改进的最小描述长度方法,从而实现对空域中的未知信源个数的检测。本发明解决了传统信源数检测技术中计算复杂度高、对非均匀噪声敏感等问题,对比已有技术,本发明更能快速鲁棒地得到待测空域内的信源个数,为应用于实际中的相控阵雷达和通信系统进行目标/用户跟踪、波达方向估计、自适应滤波处理等提供了必要的保障。
申请公布号 CN101644760A 申请公布日期 2010.02.10
申请号 CN200910091746.7 申请日期 2009.08.27
申请人 北京理工大学 发明人 黄磊;龙腾;梁正平;孔繁玲
分类号 G01S7/14(2006.01)I;G01S13/66(2006.01)I 主分类号 G01S7/14(2006.01)I
代理机构 北京理工大学专利中心 代理人 张利萍
主权项 1、一种适用于高分辨阵列的快速鲁棒的信源个数检测方法,其特征在于,包括如下步骤:(I)、阵元个数为N的阵列天线接收到L个快拍数据,根据多级维纳滤波器原理,求出多级维纳滤波器各级的最小均方误差,令多级维纳滤波器第j级的最小均方误差为ρ<sub>i</sub>,其中,j=1,2,…,M,M=N-1;(II)、把步骤(I)中求出的各级最小均方误差ρ<sub>j</sub>代入下式,得到最小描述长度mMDL(k)与假定的信源个数k的关系:<maths num="0001"><![CDATA[<math><mrow><mi>mMDL</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mi>L</mi><mrow><mo>(</mo><mi>M</mi><mo>-</mo><mi>k</mi><mo>)</mo></mrow><mi>ln</mi><mrow><mo>(</mo><mfrac><mrow><mn>1</mn><mo>/</mo><mrow><mo>(</mo><mi>M</mi><mo>-</mo><mi>k</mi><mo>)</mo></mrow><msubsup><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mi>M</mi></msubsup><msub><mi>&rho;</mi><mi>j</mi></msub></mrow><msup><mrow><mo>(</mo><msubsup><mi>&Pi;</mi><mrow><mi>j</mi><mo>=</mo><mi>k</mi><mo>+</mo><mn>1</mn></mrow><mi>M</mi></msubsup><msub><mi>&rho;</mi><mi>j</mi></msub><mo>)</mo></mrow><mrow><mn>1</mn><mo>/</mo><mrow><mo>(</mo><mi>M</mi><mo>-</mo><mi>k</mi><mo>)</mo></mrow></mrow></msup></mfrac><mo>)</mo></mrow><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mi>k</mi><mrow><mo>(</mo><mn>2</mn><mi>M</mi><mo>-</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mi>ln</mi><mi>L</mi></mrow></math>]]></maths>(III)、令整数k从0递增变化到M-1,搜索最小描述长度mMDL(k)的最小值,此时该最小值所对应的k值即为检测到的信源个数<img file="A2009100917460002C2.GIF" wi="25" he="50" />,即:<maths num="0002"><![CDATA[<math><mrow><mover><mi>q</mi><mo>^</mo></mover><mo>=</mo><mi>arg</mi><munder><mi>min</mi><mrow><mi>k</mi><mo>=</mo><mn>0,1</mn><mo>,</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>,</mo><mi>M</mi><mo>-</mo><mn>1</mn></mrow></munder><mi>mMDL</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow></math>]]></maths>
地址 100081北京市海淀区中关村南大街5号
您可能感兴趣的专利