主权项 |
一种基于合成孔径天线阵列的雷达信号被动定位方法,包括如下步骤:(1)系统初始化:(1a)测量得到待定位雷达天线旋转角度速度及旋转角度周期;(1b)按照T=Π/Ω公式,计算待定位雷达天线的旋转时间周期,其中T表示待定位雷达天线的旋转时间周期,Π表示待定位雷达天线的旋转角度周期,Ω表示待定位雷达天线的旋转角度速度;(1c)实时记录两台飞行器所处位置;(2)对待定位雷达脉冲波形的复基带信号进行离散化:(2a)分别获得两台飞行器接收待定位雷达脉冲波形的复基带信号;(2b)对两台飞行器接收的待定位雷达脉冲波形的复基带信号,分别以采样速率进行离散化采样;(3)确定虚拟到达时间差:(3a)利用麦克劳林展开级数拟合法,拟合待定位雷达天线主瓣增益方向性;(3b)利用最大值查找法,估计第1个飞行器进入待定位雷达天线主瓣的第i个扫描周期增益方向性中心的对应时间,经过可调可预设飞行时间间隔后,估计第1个飞行器进入待定位雷达天线主瓣的第j个扫描周期的增益方向性中心的对应时间,估计第2个飞行器进入待定位雷达天线主瓣的第l个扫描周期增益方向性中心的对应时间,经过可调可预设飞行时间间隔后,估计第2个飞行器进入待定位雷达天线主瓣的第m个扫描周期的增益方向性中心的对应时间;(3c)按照下式,计算虚拟到达时间差:<maths num="0001"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>Ψ</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>=</mo><msub><mover><mi>t</mi><mo>^</mo></mover><mrow><mi>j</mi><mo>,</mo><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>-</mo><msub><mover><mi>t</mi><mo>^</mo></mover><mrow><mi>i</mi><mo>,</mo><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>Ψ</mi><mrow><mi>l</mi><mo>,</mo><mi>m</mi></mrow></msub><mo>=</mo><msub><mover><mi>t</mi><mo>^</mo></mover><mrow><mi>m</mi><mo>,</mo><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>-</mo><msub><mover><mi>t</mi><mo>^</mo></mover><mrow><mi>l</mi><mo>,</mo><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0001205170150000011.GIF" wi="397" he="151" /></maths>其中,Ψ<sub>i,j</sub>表示第1个飞行器的虚拟到达时间差,<img file="FDA0001205170150000012.GIF" wi="90" he="68" />表示待定位雷达天线主瓣的第i个扫描周期的增益方向性中心对应时间的估计值,<img file="FDA0001205170150000013.GIF" wi="98" he="76" />表示待定位雷达天线主瓣的第j个扫描周期的增益方向性中心对应时间的估计值,Ψ<sub>l,m</sub>表示第2个飞行器的虚拟到达时间差,<img file="FDA0001205170150000014.GIF" wi="88" he="71" />表示待定位雷达天线主瓣的第l个扫描周期的增益方向性中心对应时间的估计值,<img file="FDA0001205170150000021.GIF" wi="100" he="71" />表示待定位雷达天线主瓣的第m个扫描周期的增益方向性中心对应时间的估计值;(4)对待定位雷达的位置进行估计:(4a)按照下式,初步计算待定位雷达位置的经度、纬度;<maths num="0002"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><msubsup><mi>x</mi><mi>o</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>y</mi><mi>o</mi><mn>2</mn></msubsup><mo>+</mo><msub><mi>A</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><msub><mi>x</mi><mi>o</mi></msub><mo>+</mo><msub><mi>B</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><msub><mi>y</mi><mi>o</mi></msub><mo>+</mo><msub><mi>Θ</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>=</mo><mn>0</mn></mtd></mtr><mtr><mtd><msubsup><mi>x</mi><mi>o</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>y</mi><mi>o</mi><mn>2</mn></msubsup><mo>+</mo><msub><mi>A</mi><mrow><mi>l</mi><mo>,</mo><mi>m</mi></mrow></msub><msub><mi>x</mi><mi>o</mi></msub><mo>+</mo><msub><mi>B</mi><mrow><mi>l</mi><mo>,</mo><mi>m</mi></mrow></msub><msub><mi>y</mi><mi>o</mi></msub><mo>+</mo><msub><mi>Θ</mi><mrow><mi>l</mi><mo>,</mo><mi>m</mi></mrow></msub><mo>=</mo><mn>0</mn></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0001205170150000022.GIF" wi="1022" he="222" /></maths>其中,x<sub>o</sub>表示待定位雷达位置的经度,y<sub>o</sub>表示待定位雷达位置的纬度,A<sub>i,j</sub>表示第1个飞行器所构建待定位雷达位置的非线性方程中待定位雷达位置经度的系数矢量,A<sub>l,m</sub>表示第2个飞行器所构建待定位雷达位置的非线性方程中待定位雷达位置经度的系数矢量,Β<sub>i,j</sub>表示第1个飞行器所构建待定位雷达位置的非线性方程中待定位雷达位置纬度的系数矢量,Β<sub>l,m</sub>表示第2个飞行器所构建待定位雷达位置的非线性方程中待定位雷达位置纬度的系数矢量,Θ<sub>i,j</sub>表示第1个飞行器所构建待定位雷达位置的非线性方程中,仅与第1个飞行器位置有关而与待定位雷达位置无关的系数矢量,Θ<sub>l,m</sub>表示第2个飞行器所构建待定位雷达位置的非线性方程中,仅与第2个飞行器位置有关而与待定位雷达位置无关的系数矢量;(4b)确定待定位雷达的位置估计结果:对待定位雷达进行第二次定位,根据两次定位结果的收敛于待定位雷达的真实位置,确定待定位雷达的位置估计结果。 |