发明名称 一种机载多主动雷达测距的地面物体高精度定位方法
摘要 本发明涉及一种机载多主动雷达测距的地面物体高精度定位方法,通过采用多主动雷达对物体交汇定位的方法,利用平台的位置信息和雷达的测距信息建立空间解析几何模型,通过求解方程组直接解算出物体点的大地坐标系位置信息,不涉及传统方法中利用时间差求解高次方程组和迭代初值的问题,大大降低了计算复杂度,并显著提高了定位精度;解决了目前机载主动雷达受载机姿态和雷达测角信息影响较大、基于辐射源的多雷达无源定位对物体的辐射源有要求等问题,实现对地面静止物体和慢速移动物体的高精度定位。
申请公布号 CN104535993A 申请公布日期 2015.04.22
申请号 CN201410683350.2 申请日期 2014.11.24
申请人 中国运载火箭技术研究院 发明人 顾鑫;王华;张尧;李潇;费智婷;姜鹏;邓志均;岑小锋
分类号 G01S13/42(2006.01)I 主分类号 G01S13/42(2006.01)I
代理机构 中国航天科技专利中心 11009 代理人 范晓毅
主权项 一种机载多主动雷达测距的地面物体高精度定位方法,其特征在于:包括如下步骤:步骤(一)、采用三架无人机对同一地面物体进行搜索,将首先搜索到所述地面物体的无人机定义为主机,主机根据在WGS‑84坐标系中自身的位置坐标、对物体的测距值和对物体的测角值解算出所述地面物体的位置坐标S<sub>0</sub>,并将位置坐标S<sub>0</sub>发送给其余两架无人机,其余两架无人机根据位置坐标S<sub>0</sub>对所述地面物体进行搜素,当搜索到地面物体后,其余两架无人机将对地面物体的测距值和无人机在WGS‑84坐标系中自身的位置坐标发送给主机;步骤(二)、设主机在WGS‑84坐标系中自身的位置坐标为P<sub>1</sub>(x<sub>1</sub>,y<sub>1</sub>,z<sub>1</sub>)、其余两架无人机自身的位置坐标分别为P<sub>2</sub>(x<sub>2</sub>,y<sub>2</sub>,z<sub>2</sub>),P<sub>3</sub>(x<sub>3</sub>,y<sub>3</sub>,z<sub>3</sub>),主机对地面物体的测距值为r<sub>1</sub>,其余两架无人机对地面物体的测距值分别为r<sub>2</sub>,r<sub>3</sub>,求解地面物体的位置坐标S(x,y,z),具体方法如下:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><msup><mrow><mo>(</mo><mi>x</mi><mo>-</mo><msub><mi>x</mi><mn>1</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mi>y</mi><mo>-</mo><msub><mi>y</mi><mn>1</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mi>z</mi><mo>-</mo><msub><mi>z</mi><mn>1</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>=</mo><msup><msub><mi>r</mi><mn>1</mn></msub><mn>2</mn></msup></mtd></mtr><mtr><mtd><msup><mrow><mo>(</mo><mi>x</mi><mo>-</mo><msub><mi>x</mi><mn>2</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mi>y</mi><mo>-</mo><msub><mi>y</mi><mn>2</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mi>z</mi><mo>-</mo><msub><mi>z</mi><mn>2</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>=</mo><msup><msub><mi>r</mi><mn>2</mn></msub><mn>2</mn></msup></mtd></mtr><mtr><mtd><msup><mrow><mo>(</mo><mi>x</mi><mo>-</mo><msub><mi>x</mi><mn>3</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mi>y</mi><mo>-</mo><msub><mi>y</mi><mn>3</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mi>z</mi><mo>-</mo><msub><mi>z</mi><mn>3</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>=</mo><msup><msub><mi>r</mi><mn>3</mn></msub><mn>2</mn></msup></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000615225110000011.GIF" wi="1413" he="241" /></maths>令:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mfenced open='' close=''><mtable><mtr><mtd><msub><mi>d</mi><mn>1</mn></msub><mo>=</mo><msup><msub><mi>r</mi><mn>1</mn></msub><mn>2</mn></msup><mo>,</mo><msub><mi>d</mi><mn>2</mn></msub><mo>=</mo><msup><msub><mi>r</mi><mn>2</mn></msub><mn>2</mn></msup><mo>,</mo><msub><mi>d</mi><mn>3</mn></msub><mo>=</mo><msup><msub><mi>r</mi><mn>3</mn></msub><mn>2</mn></msup><mo>;</mo></mtd></mtr><mtr><mtd><mi>X</mi><mo>=</mo><mi>x</mi><mo>-</mo><msub><mi>x</mi><mn>1</mn></msub><mo>,</mo><mi>Y</mi><mo>=</mo><mi>y</mi><mo>-</mo><msub><mi>y</mi><mn>1</mn></msub><mo>,</mo><mi>Z</mi><mo>=</mo><mi>z</mi><mo>-</mo><msub><mi>z</mi><mn>1</mn></msub><mo>;</mo></mtd></mtr><mtr><mtd><msub><mi>a</mi><mn>1</mn></msub><mo>=</mo><msub><mi>x</mi><mn>2</mn></msub><mo>-</mo><msub><mi>x</mi><mn>1</mn></msub><mo>,</mo><msub><mi>b</mi><mn>1</mn></msub><mo>=</mo><msub><mi>y</mi><mn>2</mn></msub><mo>-</mo><msub><mi>y</mi><mn>1</mn></msub><mo>,</mo><msub><mi>c</mi><mn>1</mn></msub><mo>=</mo><msub><mi>z</mi><mn>2</mn></msub><mo>-</mo><msub><mi>z</mi><mn>1</mn></msub><mo>;</mo></mtd></mtr><mtr><mtd><msub><mi>a</mi><mn>2</mn></msub><mo>=</mo><msub><mi>x</mi><mn>3</mn></msub><mo>-</mo><msub><mi>x</mi><mn>1</mn></msub><mo>,</mo><msub><mi>b</mi><mn>2</mn></msub><mo>=</mo><msub><mi>y</mi><mn>3</mn></msub><mo>-</mo><msub><mi>y</mi><mn>1</mn></msub><mo>,</mo><msub><mi>c</mi><mn>2</mn></msub><mo>=</mo><msub><mi>z</mi><mn>3</mn></msub><mo>-</mo><msub><mi>z</mi><mn>1</mn></msub><mo>;</mo></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000615225110000012.GIF" wi="1396" he="288" /></maths>其中d<sub>1</sub>,d<sub>2</sub>,d<sub>3</sub>均为中间变量,无含义;X,Y,Z均为中间变量,无含义;a<sub>1</sub>b<sub>1</sub>,c<sub>1</sub>,a<sub>2</sub>,b<sub>2</sub>,c<sub>2</sub>均为中间变量,无含义;结合式(1)和式(2)可以得到:d<sub>1</sub>=X<sup>2</sup>+Y<sup>2</sup>+Z<sup>2</sup>             (3)d<sub>2</sub>=(X‑a<sub>1</sub>)<sup>2</sup>+(Y‑b<sub>1</sub>)<sup>2</sup>+(Z‑c<sub>1</sub>)<sup>2</sup>                (4)d<sub>3</sub>=(X‑a<sub>2</sub>)<sup>2</sup>+(Y‑b<sub>2</sub>)<sup>2</sup>+(Z‑c<sub>2</sub>)<sup>2</sup>            (5)令式(3)减去式(4)得到式(7),式(3)减去式(5)得到式(8):d<sub>1</sub>=X<sup>2</sup>+Y<sup>2</sup>+Z<sup>2</sup>            (3)d<sub>1</sub>‑d<sub>2</sub>+a<sub>1</sub><sup>2</sup>+b<sub>1</sub><sup>2</sup>+c<sub>1</sub><sup>2</sup>=2a<sub>1</sub>X+2b<sub>1</sub>Y+2c<sub>1</sub>Z         (7)d<sub>1</sub>‑d<sub>3</sub>+a<sub>2</sub><sup>2</sup>+b<sub>2</sub><sup>2</sup>+c<sub>2</sub><sup>2</sup>=2a<sub>2</sub>X+2b<sub>2</sub>Y+2c<sub>2</sub>Z     (8)令D<sub>2</sub>=(d<sub>1</sub>‑d<sub>2</sub>+a<sub>1</sub><sup>2</sup>+b<sub>1</sub><sup>2</sup>+c<sub>1</sub><sup>2</sup>)/2,D<sub>3</sub>=(d<sub>1</sub>‑d<sub>3</sub>+a<sub>2</sub><sup>2</sup>+b<sub>2</sub><sup>2</sup>+c<sub>2</sub><sup>2</sup>)/2,则上式变为:d<sub>1</sub>=X<sup>2</sup>+Y<sup>2</sup>+Z<sup>2</sup>        (3)D<sub>2</sub>=a<sub>1</sub>X+b<sub>1</sub>Y+c<sub>1</sub>Z         (10)D<sub>3</sub>=a<sub>2</sub>X+b<sub>2</sub>Y+c<sub>2</sub>Z        (11)其中D<sub>2</sub>,D<sub>3</sub>均为中间变量,无含义;联立式(10)和式(11)可以得出如下公式:X=(b<sub>2</sub>D<sub>2</sub>‑b<sub>1</sub>D<sub>3</sub>)/(b<sub>2</sub>a<sub>1</sub>‑a<sub>2</sub>b<sub>1</sub>)+(b<sub>1</sub>c<sub>2</sub>‑c<sub>1</sub>b<sub>2</sub>)/(b<sub>2</sub>a<sub>1</sub>‑a<sub>2</sub>b<sub>1</sub>)Z   (12)Y=(a<sub>2</sub>D<sub>2</sub>‑a<sub>1</sub>D<sub>3</sub>)/(a<sub>2</sub>b<sub>1</sub>‑b<sub>2</sub>a<sub>1</sub>)+(a<sub>1</sub>c<sub>2</sub>‑c<sub>1</sub>a<sub>2</sub>)/(a<sub>2</sub>b<sub>1</sub>‑b<sub>2a1</sub>)Z   (13)令m<sub>1</sub>,m<sub>2</sub>,n<sub>1</sub>,n<sub>2</sub>满足:X=(b<sub>2</sub>D<sub>2</sub>‑b<sub>1</sub>D<sub>3</sub>)/(b<sub>2</sub>a<sub>1</sub>‑a<sub>2</sub>b<sub>1</sub>)+(b<sub>1</sub>c<sub>2</sub>‑c<sub>1</sub>b<sub>2</sub>)/(b<sub>2</sub>a<sub>1</sub>‑a<sub>2</sub>b<sub>1</sub>)Z=m<sub>1</sub>+n<sub>1</sub>Z   (14)Y=(a<sub>2</sub>D<sub>2</sub>‑a<sub>1</sub>D<sub>3</sub>)/(a<sub>2</sub>b<sub>1</sub>‑b<sub>2</sub>a<sub>1</sub>)+(a<sub>1</sub>c<sub>2</sub>‑c<sub>1</sub>a<sub>2</sub>)/(a<sub>2</sub>b<sub>1</sub>‑b<sub>2</sub>a<sub>1</sub>)Z=m<sub>2</sub>+n<sub>2</sub>Z   (15)其中m<sub>1</sub>,m<sub>2</sub>,n<sub>1</sub>,n<sub>2</sub>均为中间变量,无含义;将式(14)和式(15)代入式(3)可以得到:d<sub>1</sub>=(m<sub>1</sub>+n<sub>1</sub>Z)<sup>2</sup>+(m<sub>2</sub>+n<sub>2</sub>Z)<sup>2</sup>+Z<sup>2</sup>          (16)将上式整理得:(1+n<sub>1</sub><sup>2</sup>+n<sub>2</sub><sup>2</sup>)Z<sup>2</sup>+(2m<sub>1</sub>n<sub>1</sub>+2m<sub>2</sub>n<sub>2</sub>)Z+(m<sub>1</sub><sup>2</sup>+m<sub>2</sub><sup>2</sup>‑d<sub>1</sub>)=0   (17)令β=1+n<sub>1</sub><sup>2</sup>+n<sub>2</sub><sup>2</sup>,δ=2m<sub>1</sub>n<sub>1</sub>+2m<sub>2</sub>n<sub>2</sub>,u=m<sub>1</sub><sup>2</sup>+m<sub>2</sub><sup>2</sup>‑d<sub>1</sub>解得:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>Z</mi><mo>=</mo><mfrac><mrow><mo>-</mo><mi>&delta;</mi><mo>&PlusMinus;</mo><msqrt><msup><mi>&delta;</mi><mn>2</mn></msup><mo>-</mo><mn>4</mn><mi>&beta;u</mi></msqrt></mrow><mrow><mn>2</mn><mi>&beta;</mi></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>18</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000615225110000021.GIF" wi="1120" he="142" /></maths>其中β,δ,u均为中间变量,无含义;舍去Z为负数的值,将Z的值分别代入式(12)和式(13)即可解得X,Y的值,根据式(2)进一步解算出地面物体的坐标S(x,y,z)。
地址 100076 北京市丰台区北京9200信箱38分箱