发明名称 视场方向可调的卫星对地凝视姿态控制方法
摘要 本发明提供了一种视场方向可调的卫星对地凝视姿态控制方法,包括以下步骤:1)给定的卫星轨道状态、地面点空间坐标、相机参数以及指定视场方向对应的当地方向矢量;2)计算期望的凝视姿态四元数和角速度;3)计算误差四元数与误差角速度;4)设计飞轮控制律。通过本发明可实现视场方向可调卫星对地凝视姿态控制,既保证图像在相机视场中不发生旋转,又确保凝视成像过程中相机视场某一方向始终指向设计的当地方向矢量,便于进行图像观察与分析。该方法算法简单,运算量小,易于工程实现,具有较高的指向控制精度和较好的指向稳定度。
申请公布号 CN106054910A 申请公布日期 2016.10.26
申请号 CN201610532747.0 申请日期 2016.07.07
申请人 中国人民解放军国防科学技术大学 发明人 黄富强;连一君;曾国强;李志军;袁福;税海涛;高玉东;项军华;吴国福;韩大鹏
分类号 G05D1/08(2006.01)I 主分类号 G05D1/08(2006.01)I
代理机构 北京中济纬天专利代理有限公司 11429 代理人 陈立新
主权项 一种视场方向可调的卫星对地凝视姿态控制方法,其特征在于,包括以下步骤:步骤S1:获取卫星轨道状态、地面点空间坐标、相机参数,指定视场方向对应的当地方向矢量卫星轨道状态为X<sub>S</sub>=[R<sub>S</sub>,V<sub>S</sub>],地面点T的经纬高坐标为(L<sub>T</sub>,B<sub>T</sub>,H<sub>T</sub>),R<sub>S</sub>为卫星的地心惯性系位置,V<sub>S</sub>为卫星的地心惯性系速度,L<sub>T</sub>为地面点的地理经度,B<sub>T</sub>为地面点的地理纬度,H<sub>T</sub>为地面点的高程;设定相机投影中心在卫星体坐标下位置为<img file="FDA0001043173200000011.GIF" wi="91" he="62" />相机焦距为f,相机在星上的安装矩阵为<img file="FDA0001043173200000012.GIF" wi="85" he="66" />指定视场方向对应的当地方向矢量<img file="FDA0001043173200000013.GIF" wi="71" he="54" />在当地的地理方位角始终保持为指定的κ值,确定当地方向矢量<img file="FDA0001043173200000014.GIF" wi="66" he="55" />端点地面点P的经纬高坐标(L<sub>P</sub>,B<sub>P</sub>,H<sub>P</sub>);步骤S2:计算卫星的凝视期望四元数与期望角速度步骤S3:获取星体实际姿态四元数和星体实际姿态角速度,计算期望四元数与星体实际姿态四元数之间的误差四元数,计算期望角速度与星体实际姿态角速度之间的误差角速度按公式(11)计算误差四元数<maths num="0001"><math><![CDATA[<mrow><msub><mover><mi>q</mi><mo>&OverBar;</mo></mover><mi>e</mi></msub><mo>=</mo><msubsup><mover><mi>q</mi><mo>&OverBar;</mo></mover><mi>r</mi><mo>*</mo></msubsup><mo>&CircleTimes;</mo><mover><mi>q</mi><mo>&OverBar;</mo></mover><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001043173200000015.GIF" wi="1006" he="63" /></maths>其中,<img file="FDA0001043173200000016.GIF" wi="52" he="62" />为<img file="FDA0001043173200000017.GIF" wi="51" he="53" />的共轭四元数,<img file="FDA0001043173200000018.GIF" wi="461" he="79" />为星体实际姿态四元数,<img file="FDA0001043173200000019.GIF" wi="52" he="47" />表示四元数乘法,其中<img file="FDA00010431732000000110.GIF" wi="574" he="75" />q<sub>e0</sub>为误差四元数标部,q<sub>e</sub>=[q<sub>e1</sub> q<sub>e2</sub> q<sub>e3</sub>]<sup>T</sup>为误差四元数矢量部分;按公式(12)计算误差角速度ω<sub>e</sub>=ω‑ω<sup>*</sup>   (12)其中,ω为星体实际角速度;步骤S4:控制律设计:按公式(13)计算飞轮控制力矩<maths num="0002"><math><![CDATA[<mrow><msub><mi>L</mi><mi>c</mi></msub><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>L</mi><mi>x</mi></msub></mtd></mtr><mtr><mtd><msub><mi>L</mi><mi>y</mi></msub></mtd></mtr><mtr><mtd><msub><mi>L</mi><mi>z</mi></msub></mtd></mtr></mtable></mfenced><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><mo>-</mo><msubsup><mi>k</mi><mi>p</mi><mi>x</mi></msubsup><msub><mi>q</mi><mrow><mi>e</mi><mn>1</mn></mrow></msub><mo>-</mo><msubsup><mi>k</mi><mi>d</mi><mi>x</mi></msubsup><msub><mi>&omega;</mi><mrow><mi>e</mi><mi>x</mi></mrow></msub></mrow></mtd></mtr><mtr><mtd><mrow><mo>-</mo><msubsup><mi>k</mi><mi>p</mi><mi>y</mi></msubsup><msub><mi>q</mi><mrow><mi>e</mi><mn>2</mn></mrow></msub><mo>-</mo><msubsup><mi>k</mi><mi>d</mi><mi>y</mi></msubsup><msub><mi>&omega;</mi><mrow><mi>e</mi><mi>y</mi></mrow></msub></mrow></mtd></mtr><mtr><mtd><mrow><mo>-</mo><msubsup><mi>k</mi><mi>p</mi><mi>z</mi></msubsup><msub><mi>q</mi><mrow><mi>e</mi><mn>3</mn></mrow></msub><mo>-</mo><msubsup><mi>k</mi><mi>d</mi><mi>z</mi></msubsup><msub><mi>&omega;</mi><mrow><mi>e</mi><mi>z</mi></mrow></msub></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>13</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00010431732000000111.GIF" wi="1182" he="236" /></maths>其中,<img file="FDA00010431732000000112.GIF" wi="286" he="79" />为误差四元数控制系数,q<sub>ei</sub>(i=1,2,3)为误差四元数的矢量部分,<img file="FDA00010431732000000113.GIF" wi="286" he="63" />为角速度偏差系数,ω<sub>ei</sub>(i=x,y,z)为误差角速度,按公式(13)得到的控制力矩对卫星上的飞轮进行控制,并调整卫星姿态。
地址 410073 湖南省长沙市开福区德雅路109号