发明名称 一种分布式阵列相参合成雷达发射相参参数校准方法
摘要 本发明涉及一种分布式阵列相参合成雷达发射相参参数校准方法,包括调整单元雷达阵面一致,单元雷达阵面参考点粗定位,计算满足方向增益要求的定位精度范围,根据粗定位精度选择合适的定位估计精度,通过优化算法估计各单元雷达定位误差,最后利用估计值校准加权形成发射波束,流程如图1所示。该方法适用于分布式阵列相参合成雷达在单元雷达密集分布而测量定位精度不够时,用软件方法进行阵列流形修正,校准发射相参参数,最终实现发射波束形成。
申请公布号 CN104375129A 申请公布日期 2015.02.25
申请号 CN201410554128.2 申请日期 2014.10.17
申请人 西安电子工程研究所 发明人 游俊;黄志忠;强勇;李万玉;肖开健;远望;肖秋
分类号 G01S7/40(2006.01)I 主分类号 G01S7/40(2006.01)I
代理机构 西北工业大学专利中心 61204 代理人 王鲜凯
主权项 一种分布式阵列相参合成雷达发射相参参数校准方法,其特征在于步骤如下:步骤1:调整N个雷达站的阵面一致,以其中一站作为参考站,空间坐标为(0,0,0),其余各站阵面参考点的实际空间坐标为(x<sub>n</sub>,y<sub>n</sub>,z<sub>n</sub>),构成一个坐标系,使各站在远场具有相同的发射波束<img file="FDA0000588678340000011.GIF" wi="314" he="83" />其中:θ是俯仰角,<img file="FDA0000588678340000012.GIF" wi="46" he="59" />是方位角,θ<sub>0</sub>和<img file="FDA0000588678340000013.GIF" wi="50" he="69" />分别表示扫描方向的俯仰角和方位角,<img file="FDA0000588678340000014.GIF" wi="421" he="95" />n=1,2,…,N‑1;步骤2:通过测量定位获得其余各站相对于参考站的坐标为(x′<sub>n</sub>,y′<sub>n</sub>,z′<sub>n</sub>),根据下式:(x′<sub>n</sub>,y′<sub>n</sub>,z′<sub>n</sub>)=(x<sub>n</sub>+Δx<sub>n</sub>,y<sub>n</sub>+Δy<sub>n</sub>,z<sub>n</sub>+Δz<sub>n</sub>)得到定位误差为(Δx<sub>n</sub>,Δy<sub>n</sub>,Δz<sub>n</sub>),Δx<sub>n</sub>∈[‑δ′<sub>x,n</sub>,δ′<sub>x,n</sub>],Δy<sub>n</sub>∈[‑δ′<sub>y,n</sub>,δ′<sub>y,n</sub>],Δz<sub>n</sub>∈[‑δ′<sub>z,n</sub>,δ′<sub>z,n</sub>],n=1,2,…,N‑1,(±δ′<sub>x,n</sub>,±δ′<sub>y,n</sub>,±δ′<sub>z,n</sub>)为粗定位精度;步骤3:按照测量定位值对各站发射信号加权:<img file="FDA0000588678340000015.GIF" wi="1638" he="171" />其中:n=1,2,…,N‑1,w<sub>0</sub>=1;加权后形成方向图为:<img file="FDA0000588678340000016.GIF" wi="1542" he="163" />在<img file="FDA0000588678340000017.GIF" wi="160" he="85" />方向增益损失:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>L</mi><mi>s</mi></msub><mo>=</mo><msub><mi>G</mi><mrow><mi>s</mi><mn>0</mn></mrow></msub><mo>-</mo><msubsup><mi>G</mi><mrow><mi>s</mi><mn>0</mn></mrow><mo>&prime;</mo></msubsup><mo>=</mo><mo>-</mo><mn>20</mn><mi>lg</mi><mo>[</mo><mo>|</mo><mn>1</mn><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover><mi>E</mi><mrow><mo>(</mo><msup><mi>e</mi><mrow><mi>j</mi><mn>2</mn><mi>&pi;</mi><mfrac><mrow><mi>&Delta;</mi><msub><mi>d</mi><mi>n</mi></msub></mrow><mi>&lambda;</mi></mfrac></mrow></msup><mo>)</mo></mrow><mo>|</mo><mo>/</mo><mi>N</mi><mo>]</mo></mrow>]]></math><img file="FDA0000588678340000018.GIF" wi="962" he="184" /></maths>其中:<img file="FDA0000588678340000019.GIF" wi="1822" he="83" /><img file="FDA00005886783400000110.GIF" wi="1140" he="87" />G<sub>s0</sub>=20lg(N|F<sub>0</sub>|),<img file="FDA00005886783400000113.GIF" wi="492" he="84" /><img file="FDA00005886783400000112.GIF" wi="1464" he="181" />Δx<sub>n</sub>服从[‑δ<sub>x</sub>,δ<sub>x</sub>]范围的均匀分布,Δy<sub>n</sub>服从[‑δ<sub>y</sub>,δ<sub>y</sub>]范围的均匀分布,Δz<sub>n</sub>服从[‑δ<sub>z</sub>,δ<sub>z</sub>]范围的均匀分布,且Δx<sub>n</sub>、Δy<sub>n</sub>、Δz<sub>n</sub>三者相互独立;保证在<img file="FDA0000588678340000021.GIF" wi="164" he="85" />方向增益损失不超过L<sub>0</sub>的前提下,定位估计精度δ<sub>x</sub>,δ<sub>y</sub>,δ<sub>z</sub>满足下式:<img file="FDA0000588678340000022.GIF" wi="1709" he="188" />步骤4:根据定位估计精度δ<sub>x</sub>,δ<sub>y</sub>,δ<sub>z</sub>计算定位误差(Δx<sub>n</sub>,Δy<sub>n</sub>,Δz<sub>n</sub>),其中n=1,2,…,N‑1:<maths num="0002" id="cmaths0002"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><mi>&Delta;</mi><msub><mi>x</mi><mi>n</mi></msub><mo>=</mo><mn>2</mn><msub><mi>k</mi><mrow><mi>x</mi><mo>,</mo><mi>n</mi></mrow></msub><msub><mi>&delta;</mi><mi>x</mi></msub><mo>+</mo><msub><mi>o</mi><mrow><mi>x</mi><mo>,</mo><mi>n</mi></mrow></msub><mo>,</mo><msub><mi>k</mi><mrow><mi>x</mi><mo>,</mo><mi>n</mi></mrow></msub><mo>=</mo><mn>0</mn><mo>,</mo><mo>&PlusMinus;</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mo>&PlusMinus;</mo><msub><mi>K</mi><mrow><mi>x</mi><mo>,</mo><mi>n</mi></mrow></msub></mtd></mtr><mtr><mtd><mi>&Delta;</mi><msub><mi>y</mi><mi>n</mi></msub><mo>=</mo><mn>2</mn><msub><mi>k</mi><mrow><mi>y</mi><mo>,</mo><mi>n</mi></mrow></msub><msub><mi>&delta;</mi><mi>y</mi></msub><mo>+</mo><msub><mi>o</mi><mrow><mi>y</mi><mo>,</mo><mi>n</mi></mrow></msub><mo>,</mo><msub><mi>k</mi><mrow><mi>y</mi><mo>,</mo><mi>n</mi></mrow></msub><mo>=</mo><mn>0</mn><mo>,</mo><mo>&PlusMinus;</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mo>&PlusMinus;</mo><msub><mi>K</mi><mrow><mi>y</mi><mo>,</mo><mi>n</mi></mrow></msub></mtd></mtr><mtr><mtd><mi>&Delta;</mi><msub><mi>z</mi><mi>n</mi></msub><mo>=</mo><mn>2</mn><msub><mi>k</mi><mrow><mi>z</mi><mo>,</mo><mi>n</mi></mrow></msub><msub><mi>&delta;</mi><mi>z</mi></msub><mo>+</mo><msub><mi>o</mi><mrow><mi>z</mi><mo>,</mo><mi>n</mi></mrow></msub><mo>,</mo><msub><mi>k</mi><mrow><mi>z</mi><mo>,</mo><mi>n</mi></mrow></msub><mo>=</mo><mn>0</mn><mo>,</mo><mo>&PlusMinus;</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mo>&PlusMinus;</mo><msub><mi>K</mi><mrow><mi>z</mi><mo>,</mo><mi>n</mi></mrow></msub></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000588678340000023.GIF" wi="889" he="199" /></maths>其中,o<sub>x,n</sub>∈[‑δ<sub>x</sub>,δ<sub>x</sub>];<img file="FDA0000588678340000024.GIF" wi="528" he="98" />o<sub>y,n</sub>∈[‑δ<sub>y</sub>,δ<sub>y</sub>];<img file="FDA0000588678340000025.GIF" wi="536" he="109" />o<sub>z,n</sub>∈[‑δ<sub>z</sub>,δ<sub>z</sub>];<img file="FDA0000588678340000026.GIF" wi="530" he="92" /><img file="FDA0000588678340000027.GIF" wi="89" he="83" />表示不小于·的最小整数;对定位误差进行量化估计:<maths num="0003" id="cmaths0003"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><mi>&Delta;</mi><msub><mover><mi>x</mi><mo>^</mo></mover><mi>n</mi></msub><mo>=</mo><mn>2</mn><msub><mi>k</mi><mrow><mi>x</mi><mo>,</mo><mi>n</mi></mrow></msub><msub><mi>&delta;</mi><mi>x</mi></msub><mo>,</mo><msub><mi>k</mi><mrow><mi>x</mi><mo>,</mo><mi>n</mi></mrow></msub><mo>=</mo><mn>0</mn><mo>,</mo><mo>&PlusMinus;</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mo>&PlusMinus;</mo><msub><mi>K</mi><mrow><mi>x</mi><mo>,</mo><mi>n</mi></mrow></msub></mtd></mtr><mtr><mtd><mi>&Delta;</mi><msub><mover><mi>y</mi><mo>^</mo></mover><mi>n</mi></msub><mo>=</mo><mn>2</mn><msub><mi>k</mi><mrow><mi>y</mi><mo>,</mo><mi>n</mi></mrow></msub><msub><mi>&delta;</mi><mi>y</mi></msub><mo>,</mo><msub><mi>k</mi><mrow><mi>y</mi><mo>,</mo><mi>n</mi></mrow></msub><mo>=</mo><mn>0</mn><mo>,</mo><mo>&PlusMinus;</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mo>&PlusMinus;</mo><msub><mi>K</mi><mrow><mi>y</mi><mo>,</mo><mi>n</mi></mrow></msub></mtd></mtr><mtr><mtd><mi>&Delta;</mi><msub><mover><mi>z</mi><mo>^</mo></mover><mi>n</mi></msub><mo>=</mo><mn>2</mn><msub><mi>k</mi><mrow><mi>z</mi><mo>,</mo><mi>n</mi></mrow></msub><msub><mi>&delta;</mi><mi>z</mi></msub><mo>,</mo><msub><mi>k</mi><mrow><mi>z</mi><mo>,</mo><mi>n</mi></mrow></msub><mo>=</mo><mn>0</mn><mo>,</mo><mo>&PlusMinus;</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mo>&PlusMinus;</mo><msub><mi>K</mi><mrow><mi>z</mi><mo>,</mo><mi>n</mi></mrow></msub></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000588678340000028.GIF" wi="775" he="209" /></maths>其中,<img file="FDA0000588678340000029.GIF" wi="1656" he="106" />定位估计精度δ<sub>x</sub>,δ<sub>y</sub>,δ<sub>z</sub>由下式确定:<img file="FDA00005886783400000210.GIF" wi="1262" he="190" />步骤5:按照定位误差量化估计值对各站发射信号进行加权:<img file="FDA00005886783400000211.GIF" wi="1173" he="150" />其中,<img file="FDA00005886783400000212.GIF" wi="862" he="84" />n=1,2,…,N‑1,w<sub>0</sub>=1;在信号通过天线辐射后,利用接收装置获得远场<img file="FDA00005886783400000213.GIF" wi="164" he="83" />方向的信号,记为S:<img file="FDA00005886783400000214.GIF" wi="804" he="146" />其中s表示接收到的(0,0,0)站参考单元发射的信号,采用穷举法或遗传算法搜索使接收信号S最大的定位误差的估计值<img file="FDA00005886783400000215.GIF" wi="394" he="81" />步骤6:利用定位误差的估计值<img file="FDA0000588678340000031.GIF" wi="366" he="77" />计算各站位置估计值:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mrow><mo>(</mo><msub><mover><mi>x</mi><mo>^</mo></mover><mrow><mi>n</mi><mo>,</mo><mn>0</mn></mrow></msub><mo>,</mo><msub><mover><mi>y</mi><mo>^</mo></mover><mrow><mi>n</mi><mo>,</mo><mn>0</mn></mrow></msub><mo>,</mo><msub><mover><mi>z</mi><mo>^</mo></mover><mrow><mi>n</mi><mo>,</mo><mn>0</mn></mrow></msub><mo>)</mo></mrow><mo>=</mo><mrow><mo>(</mo><msubsup><mi>x</mi><mi>n</mi><mo>&prime;</mo></msubsup><mo>-</mo><mi>&Delta;</mi><msub><mover><mi>x</mi><mo>^</mo></mover><mrow><mi>n</mi><mo>,</mo><mn>0</mn></mrow></msub><mo>,</mo><msubsup><mi>y</mi><mi>n</mi><mo>&prime;</mo></msubsup><mo>-</mo><mi>&Delta;</mi><msub><mover><mi>y</mi><mo>^</mo></mover><mrow><mi>n</mi><mo>,</mo><mn>0</mn></mrow></msub><mo>,</mo><msubsup><mi>z</mi><mi>n</mi><mo>&prime;</mo></msubsup><mo>-</mo><mi>&Delta;</mi><msub><mover><mi>z</mi><mo>^</mo></mover><mrow><mi>n</mi><mo>,</mo><mn>0</mn></mrow></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000588678340000032.GIF" wi="978" he="92" /></maths>利用估计位置值<img file="FDA0000588678340000033.GIF" wi="288" he="74" />进行相位校准加权:<img file="FDA0000588678340000034.GIF" wi="1303" he="155" />
地址 710100 陕西省西安市长安区凤栖东路