主权项 |
1.一种太阳敏感器故障情况下的对日定向姿态控制方法,其特征在于包括如下步骤:(1)在航天器的-Z面上布局两个太阳敏感器,其中一个太阳敏感器的测量轴与卫星本体的X轴平行,测量角为φ,另一个太阳敏感器的测量轴与卫星本体的Y轴平行,测量角为θ;(2)周期性的采集两个太阳敏感器的数据,计算m<sub>x,n</sub>=tan(φ<sub>n</sub>),m<sub>y,n</sub>=tan(θ<sub>n</sub>),其中φ<sub>n</sub>,θ<sub>n</sub>分别为本周期X轴和Y轴的太阳敏感器测量值,m<sub>x,n</sub>,m<sub>y,n</sub>分别为本周期X轴和Y轴的太阳敏感器测量角度的正切值;(3)建立观测方程,并利用陀螺本周期的测量值,利用观测方程预报下一周期X轴和Y轴的太阳敏感器测量角度的正切估计值,<maths num="0001"><![CDATA[<math><mrow><msub><mover><mi>m</mi><mo>~</mo></mover><mrow><mi>x</mi><mo>,</mo><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>=</mo><msub><mover><mi>m</mi><mo>~</mo></mover><mrow><mi>x</mi><mo>,</mo><mi>n</mi></mrow></msub><mo>+</mo><mrow><mo>(</mo><msub><mi>ω</mi><mi>x</mi></msub><mo>-</mo><msub><mi>ω</mi><mi>z</mi></msub><msub><mover><mi>m</mi><mo>~</mo></mover><mrow><mi>y</mi><mo>,</mo><mi>n</mi></mrow></msub><mo>+</mo><msub><mi>ω</mi><mi>y</mi></msub><msub><mover><mi>m</mi><mo>~</mo></mover><mrow><mi>x</mi><mo>,</mo><mi>n</mi></mrow></msub><msub><mover><mi>m</mi><mo>~</mo></mover><mrow><mi>y</mi><mo>,</mo><mi>n</mi></mrow></msub><mo>+</mo><msub><mi>ω</mi><mi>x</mi></msub><msubsup><mover><mi>m</mi><mo>~</mo></mover><mrow><mi>x</mi><mo>,</mo><mi>n</mi></mrow><mn>2</mn></msubsup><mo>)</mo></mrow><mo>·</mo><mi>Δt</mi></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><msub><mover><mi>m</mi><mo>~</mo></mover><mrow><mi>y</mi><mo>,</mo><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>=</mo><msub><mover><mi>m</mi><mo>~</mo></mover><mrow><mi>y</mi><mo>,</mo><mi>n</mi></mrow></msub><mo>+</mo><mrow><mo>(</mo><msub><mi>ω</mi><mi>y</mi></msub><mo>-</mo><msub><mi>ω</mi><mi>z</mi></msub><msub><mover><mi>m</mi><mo>~</mo></mover><mrow><mi>x</mi><mo>,</mo><mi>n</mi></mrow></msub><mo>+</mo><msub><mi>ω</mi><mi>x</mi></msub><msub><mover><mi>m</mi><mo>~</mo></mover><mrow><mi>x</mi><mo>,</mo><mi>n</mi></mrow></msub><msub><mover><mi>m</mi><mo>~</mo></mover><mrow><mi>y</mi><mo>,</mo><mi>n</mi></mrow></msub><mo>+</mo><msub><mi>ω</mi><mi>y</mi></msub><msubsup><mover><mi>m</mi><mo>~</mo></mover><mrow><mi>y</mi><mo>,</mo><mi>n</mi></mrow><mn>2</mn></msubsup><mo>)</mo></mrow><mo>·</mo><mi>Δt</mi></mrow></math>]]></maths>其中<img file="FDA0000371084530000013.GIF" wi="124" he="76" /><img file="FDA0000371084530000014.GIF" wi="91" he="77" />分别为本周期X轴和Y轴的太阳敏感器测量角度的正切估计值,初值为0,<img file="FDA0000371084530000015.GIF" wi="154" he="76" /><img file="FDA0000371084530000016.GIF" wi="126" he="77" />为下一周期X轴和Y轴的太阳敏感器测量角度的正切估计值,Δt为控制周期,ω<sub>x</sub>、ω<sub>y</sub>、ω<sub>z</sub>为本周期陀螺测量的卫星本体三轴角速度;(4)判断太阳敏感器的健康状态,利用健康太阳敏感器的测量值对步骤(3)计算得到的<img file="FDA0000371084530000017.GIF" wi="152" he="76" /><img file="FDA0000371084530000018.GIF" wi="126" he="77" />进行修正,当X轴太阳敏感器健康,Y轴太阳敏感器故障时,对X轴太阳敏感器测量角度的正切估计值进行修正,修正方程为<img file="FDA0000371084530000019.GIF" wi="629" he="76" /><maths num="0003"><![CDATA[<math><mrow><msub><mi>K</mi><mi>x</mi></msub><mo><</mo><mn>0</mn><mo>,</mo><msub><mi>ω</mi><mi>z</mi></msub><msub><mi>K</mi><mi>y</mi></msub><mo>+</mo><msubsup><mi>ω</mi><mi>z</mi><mn>2</mn></msubsup><mo>></mo><mn>0</mn><mo>,</mo></mrow></math>]]></maths>当X轴太阳敏感器故障,Y轴太阳敏感器健康时,对Y轴太阳敏感器测量角度的正切估计值进行修正,修正方程为<img file="FDA00003710845300000111.GIF" wi="637" he="77" /><maths num="0004"><![CDATA[<math><mrow><msub><mi>K</mi><mi>y</mi></msub><mo><</mo><mn>0</mn><mo>,</mo><msub><mi>ω</mi><mi>z</mi></msub><msub><mi>K</mi><mi>x</mi></msub><mo>-</mo><msubsup><mi>ω</mi><mi>z</mi><mn>2</mn></msubsup><mo><</mo><mn>0</mn><mo>,</mo></mrow></math>]]></maths>其中K<sub>x</sub>,K<sub>y</sub>为修正系数;(5)更新本周期X轴和Y轴的太阳敏感器测量角度的正切估计值,<img file="FDA0000371084530000021.GIF" wi="477" he="76" />由此计算本周期X轴和Y轴的太阳敏感器测量估计值<img file="FDA0000371084530000022.GIF" wi="590" he="88" />其中<img file="FDA0000371084530000023.GIF" wi="48" he="81" />为X轴太阳敏感器测量估计值,<img file="FDA0000371084530000024.GIF" wi="49" he="70" />为Y轴太阳敏感器测量估计值;(6)依据太阳敏感器健康状态和故障太阳敏感器所在轴的测量估计值重构姿态,作为控制的输入量;当X轴的太阳敏感器健康而Y轴的太阳敏感器故障时,令<img file="FDA0000371084530000025.GIF" wi="266" he="82" />当X轴的太阳敏感器故障而Y轴的太阳敏感器健康时,令<img file="FDA0000371084530000026.GIF" wi="269" he="82" />当两轴的太阳敏感器均健康时,令<img file="FDA0000371084530000027.GIF" wi="269" he="71" />当两轴的太阳敏感器均故障时,令φ=0,θ=0;(7)利用步骤(6)得到的计算结果φ和θ,计算对日定向所需的姿态控制量,控制律为:<img file="FDA0000371084530000028.GIF" wi="424" he="75" />U<sub>y</sub>=K<sub>py</sub>θ+K<sub>dy</sub>ω<sub>y</sub>,U<sub>z</sub>=K<sub>dz</sub>(ω<sub>z</sub>-ω<sub>zb</sub>),式中U<sub>m</sub>,m=x,y,z为卫星的三轴控制量,K<sub>dm</sub>为角速度控制系数,ω<sub>m</sub>为卫星陀螺实测三轴角速度,ω<sub>zb</sub>为卫星Z轴目标角速度,K<sub>px</sub>,K<sub>py</sub>分别为滚动角和俯仰角控制系数。 |