发明名称 一种空间绳系机器人中远距离逼近位置估计方法
摘要 本发明提出一种空间绳系机器人中远距离逼近位置估计方法,利用相机进行目标方位角测量,得到目标相对空间绳系机器人的方位角,结合从空间绳系机器人的空间系绳放绳机构获得的空间绳系机器人的系绳长度信息,利用无损卡尔曼滤波方法实现对空间绳系机器人的位置估计。本发明提出的空间绳系机器人中远距离逼近位置估计方法,充分利用了空间绳系机器人的系绳长度信息,结合了单目相机的视场角信息,实现对空间绳系机器人位置进行实时估计。
申请公布号 CN103900561A 申请公布日期 2014.07.02
申请号 CN201410145305.1 申请日期 2014.04.11
申请人 西北工业大学 发明人 黄攀峰;王东科;蔡佳;孟中杰
分类号 G01C21/00(2006.01)I;G01C11/00(2006.01)I 主分类号 G01C21/00(2006.01)I
代理机构 西北工业大学专利中心 61204 代理人 陈星
主权项 1.一种空间绳系机器人中远距离逼近位置估计方法,其特征在于:循环进行以下步骤,对空间绳系机器人中远距离逼近位置进行实时估计:步骤1:利用单目视觉相机获得第k时刻目标相对空间绳系机器人的方位角η<sub>1k</sub>和η<sub>2k</sub>:<maths num="0001"><![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>&eta;</mi><mrow><mn>1</mn><mi>k</mi></mrow></msub><mo>=</mo><mi>arctan</mi><mrow><mo>(</mo><mfrac><mrow><mrow><mo>(</mo><mi>u</mi><mo>-</mo><msub><mi>u</mi><mn>0</mn></msub><mo>)</mo></mrow><msub><mi>d</mi><mi>x</mi></msub></mrow><mi>f</mi></mfrac><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>&eta;</mi><mrow><mn>2</mn><mi>k</mi></mrow></msub><mo>=</mo><mi>arctan</mi><mrow><mo>(</mo><mfrac><mrow><mrow><mo>(</mo><mi>v</mi><mo>-</mo><msub><mi>v</mi><mn>0</mn></msub><mo>)</mo></mrow><msub><mi>d</mi><mi>y</mi></msub></mrow><mi>f</mi></mfrac><mo>)</mo></mrow></mtd></mtr></mtable></mfenced></math>]]></maths>其中u、v分别为目标在相机的图像坐标系中位置,u<sub>0</sub>、v<sub>0</sub>分别为图像坐标系原点位置,d<sub>x</sub>、d<sub>y</sub>为相机单个像素的物理尺寸,f为相机的焦距;步骤2:从空间绳系机器人的空间系绳放绳机构获得第k时刻空间系绳的长度L<sub>k</sub>;步骤3:根据空间绳系机器人的状态方程和量测方程,采用无损卡尔曼滤波方法对空间绳系机器人的位置进行实时估计;其中空间绳系机器人的状态方程和量测方程为:<maths num="0002"><![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>X</mi><mi>k</mi></msub><mo>=</mo><msub><mi>&Phi;</mi><mrow><mi>k</mi><mo>,</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><msub><mi>X</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>+</mo><msub><mi>W</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>Z</mi><mi>k</mi></msub><mo>=</mo><msub><mi>H</mi><mi>k</mi></msub><mrow><mo>(</mo><msub><mi>X</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>V</mi><mi>k</mi></msub></mtd></mtr></mtable></mfenced></math>]]></maths>其中<img file="FDA0000489981390000012.GIF" wi="576" he="68" />为第k时刻空间绳系机器人的状态向量,x<sub>k</sub>、y<sub>k</sub>和z<sub>k</sub>为第k时刻空间绳系机器人在目标轨道坐标系下的位置;<img file="FDA0000489981390000013.GIF" wi="154" he="72" />和<img file="FDA0000489981390000014.GIF" wi="46" he="62" />为第k时刻空间绳系机器人在目标第k坐标系下的三轴速度;所述目标轨道坐标系中x轴为目标轨道切线方向,z轴沿地球半径方向指向地心,y轴垂直与x轴和z轴,并且满足右手定则;Φ<sub>k,k-1</sub>为状态转移矩阵,<maths num="0003"><![CDATA[<math><mrow><msub><mi>&Phi;</mi><mrow><mi>k</mi><mo>,</mo><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>=</mo><mfenced open='(' close=')'><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mi>T</mi></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mi>T</mi></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mi>T</mi></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>2</mn><mi>nT</mi></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mo>-</mo><msup><mi>n</mi><mn>2</mn></msup><mi>T</mi></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>3</mn><msup><mi>n</mi><mn>2</mn></msup><mi>T</mi></mtd><mtd><mo>-</mo><mn>2</mn><mi>nT</mi></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>T为时间间隔,n为轨道平均角速度;Z<sub>k</sub>为观测向量,<maths num="0004"><![CDATA[<math><mrow><msub><mi>Z</mi><mi>k</mi></msub><mo>=</mo><mfenced open='(' close=')'><mtable><mtr><mtd><msub><mi>&eta;</mi><mrow><mn>1</mn><mi>k</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>&eta;</mi><mrow><mn>2</mn><mi>k</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>L</mi><mi>k</mi></msub></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>由步骤1和步骤2得到;H<sub>k</sub>(X<sub>k</sub>)为观测非线性函数,<maths num="0005"><![CDATA[<math><mrow><msub><mi>H</mi><mi>k</mi></msub><mrow><mo>(</mo><msub><mi>X</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>=</mo><mfenced open='(' close=')'><mtable><mtr><mtd><mi>arctan</mi><mfrac><mrow><msub><mi>R</mi><mn>21</mn></msub><msub><mi>x</mi><mi>k</mi></msub><mo>+</mo><msub><mi>R</mi><mn>22</mn></msub><msub><mi>y</mi><mi>k</mi></msub><mo>+</mo><msub><mi>R</mi><mn>23</mn></msub><msub><mi>z</mi><mi>k</mi></msub></mrow><mrow><msub><mi>R</mi><mn>11</mn></msub><msub><mi>x</mi><mi>k</mi></msub><mo>+</mo><msub><mi>R</mi><mn>12</mn></msub><msub><mi>y</mi><mi>k</mi></msub><mo>+</mo><msub><mi>R</mi><mn>13</mn></msub><msub><mi>z</mi><mi>k</mi></msub></mrow></mfrac></mtd></mtr><mtr><mtd><mi>arctan</mi><mfrac><mrow><msub><mi>R</mi><mn>31</mn></msub><msub><mi>x</mi><mi>k</mi></msub><mo>+</mo><msub><mi>R</mi><mn>32</mn></msub><msub><mi>y</mi><mi>k</mi></msub><mo>+</mo><msub><mi>R</mi><mn>33</mn></msub><msub><mi>z</mi><mi>k</mi></msub></mrow><mrow><msub><mi>R</mi><mn>11</mn></msub><msub><mi>x</mi><mi>k</mi></msub><mo>+</mo><msub><mi>R</mi><mn>12</mn></msub><msub><mi>y</mi><mi>k</mi></msub><mo>+</mo><msub><mi>R</mi><mn>13</mn></msub><msub><mi>z</mi><mi>k</mi></msub></mrow></mfrac></mtd></mtr><mtr><mtd><msqrt><msup><mrow><mo>(</mo><msub><mi>x</mi><mi>p</mi></msub><mo>-</mo><msub><mi>x</mi><mi>k</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>y</mi><mi>p</mi></msub><mo>-</mo><msub><mi>y</mi><mi>k</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>z</mi><mi>p</mi></msub><mo>-</mo><msub><mi>z</mi><mi>k</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>R<sub>11</sub>、R<sub>12</sub>、R<sub>13</sub>、R<sub>21</sub>、R<sub>22</sub>、R<sub>23</sub>、R<sub>31</sub>、R<sub>32</sub>和R<sub>33</sub>分别为目标轨道坐标系到空间绳系机器人本体坐标系转化矩阵R的9个分量,<maths num="0006"><![CDATA[<math><mrow><mi>R</mi><mo>=</mo><msub><mi>I</mi><mn>3</mn></msub><mo>-</mo><mfrac><mrow><mn>4</mn><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>&sigma;</mi><mn>2</mn></msup><mo>)</mo></mrow></mrow><msup><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msup><mi>&sigma;</mi><mn>2</mn></msup><mo>)</mo></mrow><mn>2</mn></msup></mfrac><mo>[</mo><msup><mi>&sigma;</mi><mo>&times;</mo></msup><mo>]</mo><mo>+</mo><mfrac><mn>8</mn><msup><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msup><mi>&sigma;</mi><mn>2</mn></msup><mo>)</mo></mrow><mn>2</mn></msup></mfrac><msup><mrow><mo>[</mo><msup><mi>&sigma;</mi><mo>&times;</mo></msup><mo>]</mo></mrow><mn>2</mn></msup></mrow></math>]]></maths>I<sub>3</sub>为3×3单位矩阵,σ为利用修正罗德里格斯参数表示的空间绳系机器人的姿态角,通过空间绳系机器人自身的姿态测量装置获得,<maths num="0007"><![CDATA[<math><mrow><mo>[</mo><msup><mi>&sigma;</mi><mo>&times;</mo></msup><mo>]</mo><mo>=</mo><mfenced open='(' close=')'><mtable><mtr><mtd><mn>0</mn></mtd><mtd><mo>-</mo><msub><mi>&sigma;</mi><mn>3</mn></msub></mtd><mtd><msub><mi>&sigma;</mi><mn>2</mn></msub></mtd></mtr><mtr><mtd><msub><mi>&sigma;</mi><mn>3</mn></msub></mtd><mtd><mn>0</mn></mtd><mtd><mo>-</mo><msub><mi>&sigma;</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><mo>-</mo><msub><mi>&sigma;</mi><mn>2</mn></msub></mtd><mtd><msub><mi>&sigma;</mi><mn>1</mn></msub></mtd><mtd><mn>0</mn></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>σ<sub>1</sub>、σ<sub>2</sub>和σ<sub>3</sub>分别为σ的三个分量;(x<sub>p</sub>,y<sub>p</sub>,z<sub>p</sub>)<sup>T</sup>为释放空间绳系机器人的机动平台在目标轨道坐标系下的位置坐标;W<sub>k-1</sub>为系统白噪声向量,V<sub>k</sub>观测白噪声向量。
地址 710072 陕西省西安市友谊西路127号
您可能感兴趣的专利