发明名称 用于月球及深空巡视器导航控制性能的地面室内验证方法
摘要 用于月球及深空巡视器导航控制性能的地面室内验证方法,通过坐标转换得到真实的位置、速度和姿态信息,将真实的位置、速度和姿态信息与巡视器导航计算得到的位置、速度和姿态信息做差,得到验证导航控制性能的参数,从而实现在地面验证导航控制性能。本发明提出的验证方案具有组成简便完备、验证项目全面、实时性好、精度高、跟踪能力强、试验连续性和自主性好、试验过程操作简便、试验设备性价比高等优点,能够同时验证姿态、速度、位置等导航信息的精度以及航向角控制精度;本发明所用的iGPS系统可实时高精度测量巡视器的真实位置、姿态,为评价导航控制性能提供基准,iGPS适用于室内测量,对试验场地无特殊要求,可增强对巡视器的试验验证能力;本发明所采用的公共标识能同时为iGPS和经纬仪提供测量目标,克服了iGPS测量和经纬仪测量之间转换的困难。
申请公布号 CN102322872B 申请公布日期 2013.03.13
申请号 CN201110120256.2 申请日期 2011.05.10
申请人 北京控制工程研究所 发明人 邢琰;贾永;毛晓艳;刘祥;滕宝毅;王保丰
分类号 G01C25/00(2006.01)I 主分类号 G01C25/00(2006.01)I
代理机构 中国航天科技专利中心 11009 代理人 杨虹
主权项 1.用于月球及深空巡视器导航控制性能的地面室内验证方法,其特征在于通过以下步骤实现:第一步,构建导航控制性能地面验证系统,导航控制性能地面验证系统包括iGPS系统、真北标识设备、3台经纬仪和若干公共标识,iGPS系统包括四台激光发射器和四个接收器,真北标识设备由安装在地基台上的立方镜组成,公共标识为安装接收器的同心圆环;第二步,建立北-东-地坐标系F<sub>n</sub>和真北指示立方镜坐标系F<sub>d</sub>,得到F<sub>n</sub>和F<sub>d</sub>的转换矩阵<img file="FSB00000981325200011.GIF" wi="81" he="57" />第三步,在室内试验场地四端固定安装四台激光发射器,建立iGPS系统测量坐标系F<sub>r</sub>;第四步,建立真北指示立方镜坐标系F<sub>d</sub>与iGPS系统测量坐标系F<sub>r</sub>的转换矩阵<img file="FSB00000981325200012.GIF" wi="80" he="65" />A4.1、在室内试验场地设置4个公共标识作为控制点,利用三台经纬仪联机互瞄得到控制点在真北指示立方镜坐标系F<sub>d</sub>下的坐标[x<sub>di0</sub> y<sub>di0</sub> z<sub>di0</sub>]<sup>T</sup>,i=1,2,3,4,即<maths num="0001"><![CDATA[<math><mrow><msub><mi>C</mi><mn>1</mn></msub><mo>=</mo><msup><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>x</mi><mrow><mi>di</mi><mn>0</mn></mrow></msub></mtd><mtd><msub><mi>y</mi><mrow><mi>di</mi><mn>0</mn></mrow></msub></mtd><mtd><msub><mi>z</mi><mrow><mi>di</mi><mn>0</mn></mrow></msub></mtd></mtr></mtable></mfenced><mi>T</mi></msup><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>x</mi><mrow><mi>d</mi><mn>10</mn></mrow></msub></mtd><mtd><msub><mi>x</mi><mrow><mi>d</mi><mn>20</mn></mrow></msub></mtd><mtd><msub><mi>x</mi><mrow><mi>d</mi><mn>30</mn></mrow></msub></mtd><mtd><msub><mi>x</mi><mrow><mi>d</mi><mn>40</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>y</mi><mrow><mi>d</mi><mn>10</mn></mrow></msub></mtd><mtd><msub><mi>y</mi><mrow><mi>d</mi><mn>20</mn></mrow></msub></mtd><mtd><msub><mi>y</mi><mrow><mi>d</mi><mn>30</mn></mrow></msub></mtd><mtd><msub><mi>y</mi><mrow><mi>d</mi><mn>40</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>z</mi><mrow><mi>d</mi><mn>10</mn></mrow></msub></mtd><mtd><msub><mi>z</mi><mrow><mi>d</mi><mn>20</mn></mrow></msub></mtd><mtd><msub><mi>z</mi><mrow><mi>d</mi><mn>30</mn></mrow></msub></mtd><mtd><msub><mi>z</mi><mrow><mi>d</mi><mn>40</mn></mrow></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>A4.2、将四个接收器安装在4个公共标识上,利用iGPS系统得到控制点在iGPS系统测量坐标系F<sub>r</sub>下的坐标[x<sub>ri0</sub> y<sub>ri0</sub> z<sub>ri0</sub>]<sup>T</sup>,i=1,2,3,4,即<maths num="0002"><![CDATA[<math><mrow><msub><mi>C</mi><mn>2</mn></msub><mo>=</mo><msup><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>x</mi><mrow><mi>ri</mi><mn>0</mn></mrow></msub></mtd><mtd><msub><mi>y</mi><mrow><mi>ri</mi><mn>0</mn></mrow></msub></mtd><mtd><msub><mi>z</mi><mrow><mi>ri</mi><mn>0</mn></mrow></msub></mtd></mtr></mtable></mfenced><mi>T</mi></msup><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>x</mi><mrow><mi>r</mi><mn>10</mn></mrow></msub></mtd><mtd><msub><mi>x</mi><mrow><mi>r</mi><mn>20</mn></mrow></msub></mtd><mtd><msub><mi>x</mi><mrow><mi>r</mi><mn>30</mn></mrow></msub></mtd><mtd><msub><mi>x</mi><mrow><mi>r</mi><mn>40</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>y</mi><mrow><mi>r</mi><mn>10</mn></mrow></msub></mtd><mtd><msub><mi>y</mi><mrow><mi>r</mi><mn>20</mn></mrow></msub></mtd><mtd><msub><mi>y</mi><mrow><mi>r</mi><mn>30</mn></mrow></msub></mtd><mtd><msub><mi>y</mi><mrow><mi>r</mi><mn>40</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>z</mi><mrow><mi>r</mi><mn>10</mn></mrow></msub></mtd><mtd><msub><mi>z</mi><mrow><mi>r</mi><mn>20</mn></mrow></msub></mtd><mtd><msub><mi>z</mi><mrow><mi>r</mi><mn>30</mn></mrow></msub></mtd><mtd><msub><mi>z</mi><mrow><mi>r</mi><mn>40</mn></mrow></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>A4.3、利用公式(1)和(2)通过公式(3)得到真北指示立方镜坐标系F<sub>d</sub>与iGPS系统测量坐标系F<sub>r</sub>的转换矩阵<img file="FSB00000981325200015.GIF" wi="82" he="50" /><maths num="0003"><![CDATA[<math><mrow><msubsup><mi>C</mi><mi>d</mi><mi>r</mi></msubsup><mo>=</mo><msub><mi>C</mi><mn>2</mn></msub><msup><msub><mi>C</mi><mn>1</mn></msub><mi>T</mi></msup><msup><mrow><mo>(</mo><msub><mi>C</mi><mn>1</mn></msub><msup><msub><mi>C</mi><mn>1</mn></msub><mi>T</mi></msup><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>第五步,利用公式(4)得到北-东-地坐标系F<sub>n</sub>与iGPS系统测量坐标系F<sub>r</sub>的转换矩阵<img file="FSB00000981325200022.GIF" wi="80" he="64" /><maths num="0004"><![CDATA[<math><mrow><msubsup><mi>C</mi><mi>n</mi><mi>r</mi></msubsup><mo>=</mo><msubsup><mi>C</mi><mi>n</mi><mi>d</mi></msubsup><msubsup><mi>C</mi><mi>d</mi><mi>r</mi></msubsup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>第六步,在巡视器上安装立方镜,建立巡视器本体坐标系F<sub>b</sub>;第七步,将4个公共标识安装在巡视器上表面,在4个公共标识安装四个接收器;第八步,利用三台经纬仪联机互瞄得到第七步中安装在巡视器上表面的四个接收器在巡视器本体坐标系Fb中的坐标C<sub>3</sub>,<maths num="0005"><![CDATA[<math><mrow><msub><mi>C</mi><mn>3</mn></msub><mo>=</mo><msup><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>x</mi><mi>bi</mi></msub></mtd><mtd><msub><mi>y</mi><mi>bi</mi></msub></mtd><mtd><msub><mi>z</mi><mi>bi</mi></msub></mtd></mtr></mtable></mfenced><mi>T</mi></msup><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>x</mi><mrow><mi>b</mi><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>x</mi><mrow><mi>b</mi><mn>2</mn></mrow></msub></mtd><mtd><msub><mi>x</mi><mrow><mi>b</mi><mn>3</mn></mrow></msub></mtd><mtd><msub><mi>x</mi><mrow><mi>b</mi><mn>4</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>y</mi><mrow><mi>b</mi><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>y</mi><mrow><mi>b</mi><mn>2</mn></mrow></msub></mtd><mtd><msub><mi>y</mi><mrow><mi>b</mi><mn>3</mn></mrow></msub></mtd><mtd><msub><mi>y</mi><mrow><mi>b</mi><mn>4</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>z</mi><mrow><mi>b</mi><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>z</mi><mrow><mi>b</mi><mn>2</mn></mrow></msub></mtd><mtd><msub><mi>z</mi><mrow><mi>b</mi><mn>3</mn></mrow></msub></mtd><mtd><msub><mi>z</mi><mrow><mi>b</mi><mn>4</mn></mrow></msub></mtd></mtr></mtable></mfenced><mo>,</mo><mi>i</mi><mo>=</mo><mn>1,2,3,4</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>第九步,实时获得iGPS测量坐标系F<sub>r</sub>与巡视器本体坐标系F<sub>b</sub>的转换矩阵<img file="FSB00000981325200025.GIF" wi="80" he="68" />A9.1、安装了四个接收器的巡视器在室内试验场地运动,利用iGPS系统获得四个接收器在iGPS测量坐标系Fr中的坐标C<sub>4</sub>,<maths num="0006"><![CDATA[<math><mrow><msub><mi>C</mi><mn>4</mn></msub><mo>=</mo><msup><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>x</mi><mi>ri</mi></msub></mtd><mtd><msub><mi>y</mi><mi>ri</mi></msub></mtd><mtd><msub><mi>z</mi><mi>ri</mi></msub></mtd></mtr></mtable></mfenced><mi>T</mi></msup><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>x</mi><mrow><mi>r</mi><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>x</mi><mrow><mi>r</mi><mn>2</mn></mrow></msub></mtd><mtd><msub><mi>x</mi><mrow><mi>r</mi><mn>3</mn></mrow></msub></mtd><mtd><msub><mi>x</mi><mrow><mi>r</mi><mn>4</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>y</mi><mrow><mi>r</mi><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>y</mi><mrow><mi>r</mi><mn>2</mn></mrow></msub></mtd><mtd><msub><mi>y</mi><mrow><mi>r</mi><mn>3</mn></mrow></msub></mtd><mtd><msub><mi>y</mi><mrow><mi>r</mi><mn>4</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>z</mi><mrow><mi>r</mi><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>z</mi><mrow><mi>r</mi><mn>2</mn></mrow></msub></mtd><mtd><msub><mi>z</mi><mrow><mi>r</mi><mn>3</mn></mrow></msub></mtd><mtd><msub><mi>z</mi><mrow><mi>r</mi><mn>4</mn></mrow></msub></mtd></mtr></mtable></mfenced><mo>,</mo><mi>i</mi><mo>=</mo><mn>1,2,3,4</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>A9.2、利用第八步得到的安装在巡视器上表面的四个接收器在巡视器本体坐标系Fb中的坐标C<sub>3</sub>和步骤A9.1得到的四个接收器在iGPS测量坐标系Fr中的坐标C<sub>4</sub>,根据公式(7)得到iGPS测量坐标系Fr与巡视器本体坐标系Fb的转换矩阵<img file="FSB00000981325200027.GIF" wi="81" he="68" /><maths num="0007"><![CDATA[<math><mrow><msubsup><mi>C</mi><mi>r</mi><mi>b</mi></msubsup><mo>=</mo><msub><mi>C</mi><mn>3</mn></msub><msup><msub><mi>C</mi><mn>4</mn></msub><mi>T</mi></msup><msup><mrow><mo>(</mo><msub><mi>C</mi><mn>4</mn></msub><msup><msub><mi>C</mi><mn>4</mn></msub><mi>T</mi></msup><mo>)</mo></mrow><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>第十步,实时获得巡视器本体坐标系Fb在北-东-地坐标系Fn的位置、姿态与速度信息,A10.1、利用公式(8)得到巡视器本体坐标系Fb在iGPS测量坐标系Fr中的实时位置信息P<sub>r</sub>,<maths num="0008"><![CDATA[<math><mrow><msub><mi>P</mi><mi>r</mi></msub><mo>=</mo><mfrac><mn>1</mn><mi>k</mi></mfrac><mfenced open='[' close=']'><mtable><mtr><mtd><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>k</mi></munderover><mrow><mo>(</mo><msub><mi>x</mi><mi>ri</mi></msub><mo>-</mo><msub><mi>x</mi><mi>bi</mi></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>k</mi></munderover><mrow><mo>(</mo><msub><mi>y</mi><mi>ri</mi></msub><mo>-</mo><msub><mi>y</mi><mi>bi</mi></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>k</mi></munderover><mrow><mo>(</mo><msub><mi>z</mi><mi>ri</mi></msub><mo>-</mo><msub><mi>z</mi><mi>bi</mi></msub><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中k为接收到激光发射器有效信号的接收器个数,且k≥3;A10.2、利用公式(9)得到巡视器本体坐标系Fb在北-东-地坐标系Fn中的实时位置信息P<sub>n</sub>,<maths num="0009"><![CDATA[<math><mrow><msub><mi>P</mi><mi>n</mi></msub><mo>=</mo><msup><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>x</mi><mi>n</mi></msub></mtd><mtd><msub><mi>y</mi><mi>n</mi></msub></mtd><mtd><msub><mi>z</mi><mi>n</mi></msub></mtd></mtr></mtable></mfenced><mi>T</mi></msup><mo>=</mo><msup><msubsup><mi>C</mi><mi>n</mi><mi>r</mi></msubsup><mi>T</mi></msup><msub><mi>P</mi><mi>r</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow><mo>;</mo></mrow></math>]]></maths>A10.3、利用公式(10)得到巡视器本体坐标系Fb相对iGPS测量坐标系Fr的姿态信息<img file="FSB00000981325200033.GIF" wi="73" he="50" /><img file="FSB00000981325200034.GIF" wi="1432" he="365" />其中r、p和ψ为巡视器滚动、俯仰和航向姿态角,-90°<r<sub>r</sub><90°,-90°<p<sub>r</sub><90°,-180°<ψ<sub>r</sub>≤180°,<img file="FSB00000981325200035.GIF" wi="170" he="58" />表示<img file="FSB00000981325200036.GIF" wi="53" he="58" />矩阵中的第i′行第j′列的元素;A10.4、利用公式(11)得到巡视器本体坐标系Fb相对北-东-地坐标系Fn中的实时姿态信息<img file="FSB00000981325200037.GIF" wi="77" he="54" /><img file="FSB00000981325200038.GIF" wi="1780" he="84" />A10.5利用公式(12)获得巡视器本体坐标系Fb相对北-东-地坐标系Fn的速度信息V<sub>n</sub>,<maths num="0010"><![CDATA[<math><mrow><msub><mi>V</mi><mi>n</mi></msub><mo>=</mo><msup><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>v</mi><mi>nx</mi></msub></mtd><mtd><msub><mi>v</mi><mi>ny</mi></msub></mtd><mtd><msub><mi>v</mi><mi>nz</mi></msub></mtd></mtr></mtable></mfenced><mi>T</mi></msup><mo>=</mo><mfrac><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>x</mi><mi>n</mi></msub><mo>-</mo><msub><mi>x</mi><mi>nlast</mi></msub></mtd></mtr><mtr><mtd><msub><mi>y</mi><mi>n</mi></msub><mo>-</mo><msub><mi>y</mi><mi>nlast</mi></msub></mtd></mtr><mtr><mtd><msub><mi>z</mi><mi>n</mi></msub><mo>-</mo><msub><mi>z</mi><mi>nlast</mi></msub></mtd></mtr></mtable></mfenced><mi>&Delta;T</mi></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>12</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,x<sub>n</sub>、y<sub>n</sub>和z<sub>n</sub>为步骤A10.2得到的本次位置获取时刻的巡视器本体坐标系Fb在北-东-地坐标系Fn中的位置信息,x<sub>nlast</sub>、y<sub>nlast</sub>和z<sub>nlast</sub>为上一位置获取时刻得到的实时位置信息,ΔT是本次位置确定时刻距上次位置获取时刻的时间差;第十一步,导航控制性能验证,A11.1、巡视器自主导航运算实时得到位置P<sub>m</sub>、姿态<img file="FSB00000981325200042.GIF" wi="51" he="37" />和速度V<sub>m</sub>;A11.2、根据步骤A11.1得到的位置P<sub>m</sub>、姿态<img file="FSB00000981325200043.GIF" wi="51" he="38" />和速度V<sub>m</sub>和第十步的得到的位置P<sub>n</sub>、姿态<img file="FSB00000981325200044.GIF" wi="48" he="39" />和速度V<sub>n</sub>,利用公式组(13)得到导航性能评价参数ΔP(t)、ΔV(t)、<img file="FSB00000981325200045.GIF" wi="107" he="48" />和ψ<sub>c</sub>(t),ΔP(t)=|P<sub>n</sub>(t)-P<sub>m</sub>(t)|ΔV(t)=|V<sub>n</sub>(t)-V<sub>m</sub>(t)|(13)<img file="FSB00000981325200046.GIF" wi="443" he="51" />Δψ(t)=|ψ<sub>n</sub>(t)-ψ<sub>c</sub>(t)|其中,t为时间,且t=0~T,T为试验时间,ψ<sub>c</sub>(t)是t时刻的期望航向角度,ψ<sub>n</sub>(t)是t时刻的巡视器本体坐标系Fb相对北-东-地坐标系Fn中的航向角度;A11.3、对步骤A11.2得到导航性能评价参数ΔP(t)、ΔV(t)、<img file="FSB00000981325200047.GIF" wi="109" he="48" />和ψ<sub>c</sub>(t)在试验时间T内的均值和均方差,利用均值和均方差来评价巡视器导航控制的精度和稳定度和巡视器航向角控制精度和稳定度。
地址 100080 北京市2729信箱