主权项 |
一种视场方向可调的卫星对地凝视姿态控制方法,其特征在于,包括以下步骤:步骤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>‾</mo></mover><mi>e</mi></msub><mo>=</mo><msubsup><mover><mi>q</mi><mo>‾</mo></mover><mi>r</mi><mo>*</mo></msubsup><mo>⊗</mo><mover><mi>q</mi><mo>‾</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>ω</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>ω</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>ω</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)得到的控制力矩对卫星上的飞轮进行控制,并调整卫星姿态。 |