发明名称 一种阵列通道不一致性误差快速校正方法
摘要 本发明公开一种阵列通道不一致性误差快速校正方法,使用单个校正源对阵列通道不一致性误差进行校正;同时利用校正源的方位和波形信息对阵列通道不一致性误差参数进行估计;无需估计协方差矩阵和进行特征值分解,具有较小的运算量;本发明的优点是利用阵列源的方位和波形信息对阵列进行校正,无需估计协方差矩阵和进行特征值分解,运算量较小,且具有与特征分解方法基本相同的幅相参数估计性能。EACDM算法和特征分解法的计算复杂度分别为O(3MN)和O(M<sup>2</sup>N+4M<sup>3</sup>/3),本发明的运算量仅为O(MN)。
申请公布号 CN105044660A 申请公布日期 2015.11.11
申请号 CN201510340223.7 申请日期 2015.06.18
申请人 许昌学院 发明人 张柯;程菊明;孟雷;王闯;王奎甫;张向群
分类号 G01S3/02(2006.01)I 主分类号 G01S3/02(2006.01)I
代理机构 代理人
主权项 一种阵列通道不一致性误差快速校正方法,其特征是:(1)使用单个校正源对阵列通道不一致性误差进行校正;(2)同时利用校正源的方位和波形信息对阵列通道不一致性误差参数进行估计;(3)无需估计协方差矩阵和进行特征值分解,具有较小的运算量;在阵列的远场处放置一个窄带校正源S<sub>s</sub>(t),其相对于阵列的方位是θ<sub>s</sub>,则阵列输出可表示为:X(t)=Γa<sub>s</sub>(θ)S<sub>s</sub>(t)+N(t)其中,待校正阵列为阵元数目为M的半波长均匀线阵,X(t)=[x<sub>1</sub>(t),x<sub>2</sub>(t),…,x<sub>M</sub>(t)]<sup>T</sup>为M×1维观测的数据向量,N(t)=[n<sub>1</sub>(t),n<sub>2</sub>(t),…,n<sub>M</sub>(t)]<sup>T</sup>为M×1维零均值高斯白噪声向量。a(θ<sub>s</sub>)=[1,e<sup>‑jω</sup>,e<sup>‑j2ω</sup>,…,e<sup>‑j(M‑1)ω</sup>]<sup>T</sup>为理想的阵列流型导向矢量,ω=2πd sin(θ<sub>s</sub>)/λ,<img file="FDA0000741612850000011.GIF" wi="924" he="73" />为包含幅相信息的M×M维对角矩阵,g<sub>i</sub>和<img file="FDA0000741612850000012.GIF" wi="51" he="54" />(i=1,2,...,M)分别表示第i个阵元的增益和相位;使用如下公式求出存在幅相误差的导向矢量:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>h</mi><mo>=</mo><mi>E</mi><mo>[</mo><mi>X</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><msubsup><mi>S</mi><mi>s</mi><mi>H</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>]</mo><mo>/</mo><msub><mrow><mo>|</mo><mo>|</mo><mi>E</mi><mo>[</mo><mi>X</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><msubsup><mi>S</mi><mi>s</mi><mi>H</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>]</mo><mo>|</mo><mo>|</mo></mrow><mn>2</mn></msub></mrow>]]></math><img file="FDA0000741612850000013.GIF" wi="590" he="90" /></maths>式中,H表示共轭转置,||·||<sub>2</sub>表示2‑范数;(4)通过对上式的推导,可得h与Γa<sub>s</sub>(θ)的关系如下<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>h</mi><mo>=</mo><mfrac><mrow><mi>&Gamma;</mi><msub><mi>a</mi><mi>s</mi></msub><mrow><mo>(</mo><mi>&theta;</mi><mo>)</mo></mrow></mrow><msqrt><mi>M</mi><mrow><mo>(</mo><msubsup><mi>g</mi><mn>1</mn><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>g</mi><mn>2</mn><mn>2</mn></msubsup><mo>+</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>+</mo><msubsup><mi>g</mi><mi>M</mi><mn>2</mn></msubsup><mo>)</mo></mrow></msqrt></mfrac></mrow>]]></math><img file="FDA0000741612850000014.GIF" wi="499" he="156" /></maths>(5)设<img file="FDA0000741612850000015.GIF" wi="44" he="68" />为h归一化后的值,则由上式可得第i(i=1,2,...,M)个阵元Γ<sub>i</sub>的估计值:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mover><mi>&Gamma;</mi><mo>^</mo></mover><mi>i</mi></msub><mo>=</mo><mfrac><msub><mover><mi>h</mi><mo>&OverBar;</mo></mover><mi>i</mi></msub><mrow><msub><mi>a</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>&theta;</mi><mi>s</mi></msub><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0000741612850000016.GIF" wi="219" he="141" /></maths>其中<img file="FDA0000741612850000017.GIF" wi="48" he="77" />和a<sub>i</sub>(θ<sub>s</sub>)分别为向量<img file="FDA0000741612850000018.GIF" wi="42" he="69" />和a(θ<sub>s</sub>)中的第i个元素;(6)由上式可得增益和相位的估计值为:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mover><mi>g</mi><mo>^</mo></mover><mi>i</mi></msub><mo>=</mo><mo>|</mo><msub><mover><mi>&Gamma;</mi><mo>^</mo></mover><mi>i</mi></msub><mo>|</mo><mo>,</mo></mrow>]]></math><img file="FDA00007416128500000112.GIF" wi="183" he="121" /></maths><img file="FDA0000741612850000019.GIF" wi="230" he="123" />(7)由上式得到的增益和相位估计值<img file="FDA00007416128500000110.GIF" wi="52" he="82" />和<img file="FDA00007416128500000111.GIF" wi="50" he="82" />对阵列进行校正。
地址 461000 河南省许昌市魏都区八一路88号许昌学院