[发明专利]一种唯星敏定姿定速的欠测量干扰抑制姿态控制方法有效
申请号: | 201910974453.7 | 申请日: | 2019-10-14 |
公开(公告)号: | CN110775299B | 公开(公告)日: | 2021-04-13 |
发明(设计)人: | 于萍;张洪华;陈尧;王泽国;李骥;关轶峰;程铭;张晓文;杨巍;于洁;王志文;王华强 | 申请(专利权)人: | 北京控制工程研究所 |
主分类号: | B64G1/24 | 分类号: | B64G1/24;B64G1/36 |
代理公司: | 中国航天科技专利中心 11009 | 代理人: | 张晓飞 |
地址: | 100080 *** | 国省代码: | 北京;11 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 一种 唯星敏定姿定速 测量 干扰 抑制 姿态 控制 方法 | ||
1.一种唯星敏定姿定速的欠测量干扰抑制姿态控制方法,其特征在于步骤如下:
(1)星敏感器接收控制计算机发出的同步信号,使得星敏感器测量数据的测量时刻与控制计算机的控制周期同步;
(2)如果当前控制周期有新的星敏感器测量数据发送回来,则执行步骤(3a)-(3h);如果当前控制周期没有新的星敏感器测量数据发送回来,则执行步骤(4a)-(4c);
(3a)采用星敏感器测量的姿态四元素数据直接计算出航天器当前时刻t1的惯性姿态角θi;
(3b)采用当前时刻t1的星敏感器测量的姿态四元素数据q1和当前时刻t1前的某一时刻t2的星敏感器测量的姿态四元素数据q2,计算当前姿态平均角速度ω1,并记
(3c)采用t2时刻的星敏感器测量的姿态四元素数据q2和t2时刻前的某一时刻t3的星敏感器测量的姿态四元素数据q3,计算次时刻t2的姿态角速度ω2;
(3d)计算t1时刻至t3时刻间的平均角速度ω;
(3e)计算t1时刻至t3时刻间的星体受到的力矩冲量dT;
(3f)计算t1时刻至t3时刻间的星体受到的干扰力矩冲量dMd;
(3g)计算得到干扰力矩Md;
(3h)采用α-β滤波对干扰力矩计算值进行滤波估计,获得当前时刻干扰力矩估计值Md,i;并转到步骤(5);
(4a)根据运动学方程,采用外推方法计算出航天器当前时刻t1惯性姿态角θi;
(4b)假设干扰力矩不变,得到当前时刻干扰力矩,即
Md,i=Md,i-1
其中Md,i-1是上一拍的干扰力矩估计值,Md,0=0;
(4c)根据动力学方程,采用外推方法计算出航天器当前惯性姿态角速度并转到步骤(5);
(5)计算控制器的输出,得到控制力矩,完成航天器的姿态控制。
2.根据权利要求1所述的一种唯星敏定姿定速的欠测量干扰抑制姿态控制方法,其特征在于:所述步骤(3b)中计算当前姿态平均角速度ω1的具体过程为:
计算得到姿态四元素数据q1和姿态四元素数据q2的误差四元数qe:
qe(1)=-q2(1)*q1(4)-q2(2)*q1(3)+q2(3)*q1(2)+q2(4)*q1(1)
qe(2)=q2(1)*q1(3)-q2(2)*q1(4)-q2(3)*q1(1)+q2(1)*q1(2)
qe(3)=-q2(1)*q1(2)+q2(2)*q1(1)-q2(3)*q1(4)+q2(4)*q1(3)
qe(4)=q2(1)*q1(1)+q2(2)*q1(2)+q2(3)*q1(3)+q2(4)*q1(4)
其中q(i),i=1,2,3,4,为四元素q的第i个元素,所述四元素q包括q1、q2和qe,则当前姿态平均角速度ω1
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于北京控制工程研究所,未经北京控制工程研究所许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201910974453.7/1.html,转载请声明来源钻瓜专利网。