[发明专利]一种基于反馈线性化技术的航天器相对姿态控制方法有效
申请号: | 201310308568.5 | 申请日: | 2013-07-22 |
公开(公告)号: | CN103412571A | 公开(公告)日: | 2013-11-27 |
发明(设计)人: | 李鹏;岳晓奎;袁建平 | 申请(专利权)人: | 西北工业大学 |
主分类号: | G05D1/08 | 分类号: | G05D1/08 |
代理公司: | 西北工业大学专利中心 61204 | 代理人: | 顾潮琪 |
地址: | 710072 *** | 国省代码: | 陕西;61 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | 本发明提供了一种基于反馈线性化技术的航天器相对姿态控制方法,首先建立两个航天器的相对姿态动力学模型,然后基于反馈线性化技术进行相对姿态控制器设计。本发明直接采用非线性相对姿态模型进行控制器设计,并采用反馈线性化技术设计相对姿态控制器。在保证模型精度和控制精度的同时,所设计的控制器简便实用,对星载计算机没有特殊要求,特别适用于需要实时解算且计算能力有限的空间任务。 | ||
搜索关键词: | 一种 基于 反馈 线性化 技术 航天器 相对 姿态 控制 方法 | ||
【主权项】:
1.一种基于反馈线性化技术的航天器相对姿态控制方法,其特征在于包括下述步骤:步骤一、建立两个航天器的相对姿态动力学模型追踪航天器和目标航天器的姿态动力学方程为:J c ω · c + ω c × J c ω c = T cb + T cd - - - ( 1 ) ]]>J t ω · t + ω t × J t ω t = T td - - - ( 2 ) ]]> 式中,Jc、Jt分别为追踪航天器和目标航天器的转动惯量,ωc、ωt分别为追踪航天器本体系和目标航天器本体系相对于惯性系的旋转角速度,Tcd、Ttd分别为追踪航天器和目标航天器受到的姿态干扰力矩,Tcb为追踪航天器施加的姿态控制力矩;使用四元数来描述航天器的姿态运动,则追踪航天器和目标航天器的姿态以及两个航天器的相对姿态qe可表示为q c = q c 0 q ‾ c T T , q t = q t 0 q ‾ t T T , q e = q e 0 q ‾ e T T ; ]]> 其中,qc0,qt0和qe0为四元数的标量部分,和为四元数的矢量部分;定义为向量的反对称矩阵,则两个航天器的相对姿态四元数qe为式中,为qt的共轭四元数,q t * = q t 0 - q ‾ t T T , ]]> I3为三阶单位矩阵;定义和分别为追踪航天器本体系和目标航天器本体系相对惯性系的转换矩阵,则目标航天器本体系stb到追踪航天器本体系scb的转换矩阵为C cb tb = C cb i ( C tb i ) T - - - ( 4 ) ]]>C cb i = ( q c 0 2 - q ‾ c T q ‾ c ) I 3 + 2 q ‾ c q ‾ c T - 2 q c 0 S ( q ‾ c ) - - - ( 5 ) ]]>C tb i = ( q t 0 2 - q ‾ t T q ‾ t ) I 3 + 2 q ‾ t q ‾ t T - 2 q t 0 S ( q ‾ t ) - - - ( 6 ) ]]> 其中,I3为三阶单位矩阵,为向量的反对称矩阵;定义ωe为追踪航天器相对目标航天器的角速度,则在追踪航天器体坐标系scb中可表示为ω e = ω c - C cb tb ω t - - - ( 7 ) ]]> 对进行求导后,两端同左乘矩阵Jc可得J c ω · e = J c ω · c - J c C cb tb ω · t - J c ω e × ( ω e - ω c ) - - - ( 8 ) ]]> 将(1),(2),(7)代入上式,经整理,可以得到J c ω · e + ( ω e + C cb tb ω t ) × J c ( ω e + C cb tb ω t ) - J c ω e × C cb tb ω t - J c C cb tb J t - 1 ( ω t × J t ω t ) = T cb + T cd - J c C cb tb J t - 1 ( T td ) - - - ( 9 ) ]]> 将上式展开,则有J c ω · e - S [ J c ( ω e + C cb tb ω t ) ] ω e + S ( C cb tb ω t ) J c ω e + J c S ( C cb tb ω t ) ω e - J c C cb tb J t - 1 S ( ω t ) J t ω t + S ( C cb tb ω t ) J c ( C cb tb ω t ) = T cb + ( T cd - J c C cb tb J t - 1 T td ) - - - ( 10 ) ]]> 将相对姿态动力学方程化为如下级联形式:q · e = 1 2 Ω ( q e ) e ω = 1 2 - q e 1 - q e 2 - q e 3 q e 0 - q e 3 q e 2 q e 3 q e 0 - q e 1 - q e 2 q e 1 q e 0 ω e - - - ( 11 ) ]]>J c ω · e = - C ω e - N ( ω e ) + T c + T d - - - ( 12 ) ]]> 上式当中各项的具体表达如下,其中Ω ( q e ) = - q e 1 - q e 2 - q e 3 q e 0 - q e 3 q e 2 q e 3 q e 0 - q e 1 - q e 2 q e 1 q e 0 , ]]> Tc=Tcb,T d = T cd - J c C cb tb J t - 1 T td , N ( ω e ) = S ( C cb tb ω t ) J c ( C cb tb ω t ) - J c C cb tb J t - 1 S ( ω t ) J t ω t , ]]>C = - S [ J c ( ω e + C cb tb ω t ) ] + S ( C cb tb ω t ) J c + J c S ( C cb tb ω t ) ; ]]> 由于为反对称矩阵,且S ( C cb tb ω t ) J c + J c S ( C cb tb ω t ) = - [ S ( C cb tb ω t ) J c + J c S ( C cb tb ω t ) ] T , ]]> 因此矩阵C为反对称矩阵;步骤二、控制律U=Tc,设计控制律其中,eω=ωe-ωd,eq=qe-qd,定义qd和ωd分别为理想相对四元数,理想相对角速度,kp为比例系数,kd为微分系数,则有将(13)代入式J c ω · e = - C ω e - N ( ω e ) + T c + T d , ]]> 可以得到
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于西北工业大学,未经西北工业大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/patent/201310308568.5/,转载请声明来源钻瓜专利网。
- 上一篇:一种松弛因子的迭代计算方法
- 下一篇:一种番荔枝储存时间的检测方法