发明名称 一种基于背景噪声频谱起伏特性的宽带阵列信号处理方法
摘要 本发明涉及一种基于背景噪声频谱起伏特性的宽带阵列信号处理方法,该处理方法为:首先通过多次快拍累积得到较准确的频域协方差矩阵估计,同时这也是后续对角减载系数以及后置加权系数得到准确求解的基础;然后对频域协方差矩阵进行特征分解求取对角减载系数,将频域协方差矩阵主对角线元素乘以减载系数后得到新的频域协方差矩阵;使用减载后的频域协方差矩阵进行MVDR波束形成,得到各个频点的空间能量谱;最后使用最优频段加权系数对各个频点的空间能量谱进行加权求和,得到最终的空间能量谱;本发明在提高弱目标检测性能的同时,改善多目标分辨性能,该技术较常规波束形成器以及MVDR波束形成器具有更好的探测性能。
申请公布号 CN105445723A 申请公布日期 2016.03.30
申请号 CN201510975958.7 申请日期 2015.12.22
申请人 中国船舶重工集团公司第七一五研究所 发明人 周彬;马启明;孙微;杜栓平
分类号 G01S7/52(2006.01)I;G01S7/539(2006.01)I 主分类号 G01S7/52(2006.01)I
代理机构 杭州九洲专利事务所有限公司 33101 代理人 陈继亮
主权项 一种基于背景噪声频谱起伏特性的宽带阵列信号处理方法,其特征在于:该处理方法包括如下步骤:1)、通过M次时域采样得到阵列信号x(m),为N×M矩阵,N为阵元数;然后进行傅里叶变换得到N×M维阵元域频域信号X<sub>k</sub>(f),k表示第k次频域快拍,再计算处理频段内各个频点所对应的频域协方差矩阵<img file="FDA0000886253620000011.GIF" wi="487" he="80" />2)、对各个频域协方差矩阵进行K次频域快拍的累积,积分时间可设定,得到累加后的阵列频域协方差矩阵<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mover><mi>R</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mi>K</mi></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><msub><mover><mi>R</mi><mo>^</mo></mover><mi>k</mi></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000886253620000012.GIF" wi="423" he="135" /></maths>3)、利用特征值分解得到<img file="FDA0000886253620000013.GIF" wi="116" he="78" />的特征值μ<sub>n</sub>,n=1,2,…,N,然后求取特征值的最小值min(μ<sub>n</sub>)以及协方差矩阵主对角线上的最大元素<img file="FDA0000886253620000014.GIF" wi="367" he="102" />从而得到对角减载系数λ<sub>op</sub>为<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>&lambda;</mi><mrow><mi>o</mi><mi>p</mi></mrow></msub><mo>=</mo><mn>1</mn><mo>-</mo><mfrac><mrow><mi>min</mi><mrow><mo>(</mo><msub><mi>&mu;</mi><mi>n</mi></msub><mo>)</mo></mrow></mrow><mrow><mi>m</mi><mi>a</mi><mi>x</mi><mrow><mo>(</mo><mi>d</mi><mi>i</mi><mi>a</mi><mi>g</mi><mo>(</mo><mrow><mover><mi>R</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow></mrow><mo>)</mo><mo>)</mo></mrow></mrow></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000886253620000015.GIF" wi="550" he="159" /></maths>4)、将<img file="FDA0000886253620000016.GIF" wi="118" he="80" />主对角线上的元素均乘以λ<sub>op</sub>,得到对角减载后的协方差矩阵<img file="FDA0000886253620000017.GIF" wi="167" he="85" />5)、对<img file="FDA0000886253620000018.GIF" wi="147" he="79" />进行MVDR波束形成得到各个频点的空间能量谱P<sup>λ</sup>(f,θ),并使用P<sup>λ</sup>(f,θ)中的最小值代替噪声功率作为最优频段加权系数<img file="FDA0000886253620000019.GIF" wi="479" he="150" />6)、使用a(f)对处理频段[fL,fH]内各频点的空间能量谱进行最优后置频段加权处理,得到最终的输出空间能量谱为<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msup><mi>P</mi><mi>&lambda;</mi></msup><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>f</mi><mo>=</mo><mi>f</mi><mi>L</mi></mrow><mrow><mi>f</mi><mi>H</mi></mrow></munderover><mi>a</mi><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><msup><mi>P</mi><mi>&lambda;</mi></msup><mrow><mo>(</mo><mi>f</mi><mo>,</mo><mi>&theta;</mi><mo>)</mo></mrow><mo>.</mo></mrow>]]></math><img file="FDA00008862536200000110.GIF" wi="542" he="147" /></maths>
地址 311499 浙江省杭州市富阳区富春街道桂花西路82号