发明名称 基于修正因子的最小方差谱估计方法
摘要 本发明公开了一种基于修正因子的最小方差谱估计方法,主要解决现有技术空间谱主瓣较宽,旁瓣较高的缺点。其实现过程是:1)提取均匀线阵的接收数据x(t);2)计算x(t)的协方差矩阵<img file="DDA0001018012320000011.GIF" wi="67" he="62" />并对其求逆,得到接收数据的协方差逆矩阵<img file="DDA0001018012320000012.GIF" wi="91" he="70" />3)构造N×N阶对角矩阵,根据<img file="DDA0001018012320000013.GIF" wi="72" he="71" />和对角矩阵Y,构造修正因子β,并计算Capon空间谱函数;4)根据Capon空间谱函数以及修正因子β,得到空间谱函数P<sub>p</sub>(θ);5)根据P<sub>p</sub>(θ),对波达方向进行最大似然估计,得到波达方向的估计值<img file="DDA0001018012320000014.GIF" wi="67" he="70" />本发明具有空间谱主瓣较窄,旁瓣较低的优点,且在分辨能力和测向稳健性上均优于现有MUSIC算法和Capon算法,可用于目标定位与跟踪。
申请公布号 CN106125056A 申请公布日期 2016.11.16
申请号 CN201610424037.6 申请日期 2016.06.13
申请人 西安电子科技大学 发明人 陈伯孝;余方伟;杨明磊
分类号 G01S7/41(2006.01)I 主分类号 G01S7/41(2006.01)I
代理机构 陕西电子工业专利中心 61205 代理人 王品华;朱红星
主权项 一种基于修正因子的最小方差谱估计方法,包括以下步骤:1)根据均匀线阵天线的阵列结构,获得整个均匀阵列的接收数据x(t);2)计算接收数据x(t)的协方差矩阵<img file="FDA0001018012290000011.GIF" wi="67" he="63" />并对其求逆,得到接收数据的协方差逆矩阵<img file="FDA0001018012290000012.GIF" wi="99" he="78" />3)构造N×N阶对角矩阵:<img file="FDA0001018012290000013.GIF" wi="543" he="303" />其中,ρ为常数,ρ∈[10<sup>‑8</sup>,10<sup>‑4</sup>];e<sup>‑jκd(i‑1)ρ</sup>为对角矩阵Y的第i个对角线元素,i=1,2…N,N表示阵元数,κ为波数,d为阵元间距,j为虚数单位,e为自然常数;4)根据协方差逆矩阵<img file="FDA0001018012290000014.GIF" wi="75" he="69" />和对角矩阵Y,构造修正因子β:4a)根据协方差逆矩阵<img file="FDA0001018012290000015.GIF" wi="74" he="71" />和对角矩阵Y,计算三个数值不同的扫描参数ξ<sub>1</sub>、ξ<sub>2</sub>和ξ<sub>3</sub>:<maths num="0001"><math><![CDATA[<mrow><msub><mi>&xi;</mi><mn>1</mn></msub><mo>=</mo><msup><mi>a</mi><mi>H</mi></msup><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow><msubsup><mover><mi>R</mi><mo>^</mo></mover><mi>x</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mi>a</mi><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow><mo>,</mo><msub><mi>&xi;</mi><mn>2</mn></msub><mo>=</mo><msup><mi>a</mi><mi>H</mi></msup><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow><msubsup><mover><mi>R</mi><mo>^</mo></mover><mi>x</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mi>Y</mi><mi>a</mi><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow><mo>,</mo><msub><mi>&xi;</mi><mn>3</mn></msub><mo>=</mo><msup><mi>a</mi><mi>H</mi></msup><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow><msup><mi>Y</mi><mi>H</mi></msup><msubsup><mover><mi>R</mi><mo>^</mo></mover><mi>x</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mi>Y</mi><mi>a</mi><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0001018012290000016.GIF" wi="1485" he="79" /></maths>其中:θ为扫描角度,a(θ)为阵列导向矢量,a(θ)=[1,e<sup>jκdsinθ</sup>,…,e<sup>jκ(N‑1)dsinθ</sup>]<sup>T</sup>,上标T表示转置,上标H表示共轭转置;4b)根据4a)中的第一扫描参数ξ<sub>1</sub>,得到Capon空间谱函数:<img file="FDA0001018012290000017.GIF" wi="294" he="127" />4c)根据4a)中的第一扫描参数ξ<sub>1</sub>、第二扫描参数ξ<sub>2</sub>和第三扫描参数ξ<sub>3</sub>,计算出Capon空间谱函数P<sub>capon</sub>(θ)的修正因子:<img file="FDA0001018012290000018.GIF" wi="270" he="135" />其中上标*表示共轭;5)根据4b)中的Capon空间谱函数P<sub>capon</sub>(θ)以及4c)中的修正因子β,得到空间谱函数:P<sub>p</sub>(θ)=βP<sub>capon</sub>(θ);6)根据空间谱函数P<sub>p</sub>(θ),对波达方向进行最大似然估计,得到波达方向的估计值<img file="FDA0001018012290000019.GIF" wi="58" he="71" />
地址 710071 陕西省西安市太白南路2号