主权项 |
一种基于背景噪声频谱起伏特性的宽带阵列信号处理方法,其特征在于:该处理方法包括如下步骤: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>Σ</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>λ</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>μ</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>λ</mi></msup><mrow><mo>(</mo><mi>θ</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>Σ</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>λ</mi></msup><mrow><mo>(</mo><mi>f</mi><mo>,</mo><mi>θ</mi><mo>)</mo></mrow><mo>.</mo></mrow>]]></math><img file="FDA00008862536200000110.GIF" wi="542" he="147" /></maths> |