[发明专利]面向输入饱和的迹向欠驱动航天器编队重构控制方法在审
申请号: | 201610186011.2 | 申请日: | 2016-03-29 |
公开(公告)号: | CN105786013A | 公开(公告)日: | 2016-07-20 |
发明(设计)人: | 黄煦;闫野;周洋;杨跃能 | 申请(专利权)人: | 中国人民解放军国防科学技术大学 |
主分类号: | G05D1/10 | 分类号: | G05D1/10 |
代理公司: | 北京中济纬天专利代理有限公司 11429 | 代理人: | 陈立新 |
地址: | 410073 湖*** | 国省代码: | 湖南;43 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | 本发明提供一种面向输入饱和的迹向欠驱动航天器编队重构控制方法,该方法针对存在输入饱和作用的圆轨道迹向欠驱动航天器编队构型重构控制问题,建立了其动力学模型。基于该动力学模型,分析了缺失迹向控制情况下的系统能控性以及编队重构可行性。以此模型为受控对象,构建了辅助系统以解决输入饱和问题,同时采用反步控制方法构建了迹向欠驱动情况下的闭环控制律。该方法能够完成存在输入饱和的圆轨道迹向欠驱动航天器编队构型重构控制,且闭环系统一致最终有界稳定,对外部摄动及模型误差具有良好的鲁棒性和动态性能。 | ||
搜索关键词: | 面向 输入 饱和 驱动 航天器 编队 控制 方法 | ||
【主权项】:
1.一种面向输入饱和的迹向欠驱动航天器编队重构控制方法,其特征在于,包括以下步骤:步骤S100:给定待重构的名义构型:根据待重构的名义构型,计算对应的名义相对运动状态X2d ,其中,X2d 的下标2代表缺失迹向控制加速度的欠驱动情况;步骤S200:误差量计算:对当前构型计算实际相对运动状态X2 ,由此计算当前构型与所述名义相对运动状态之间的误差量e2 ,e2 按公式(1)计算: 步骤S300:控制律构建:采用反步控制方法构建迹向欠驱动航天器编队重构控制律,计算实际控制量U2 ;其中,实际相对运动状态X2 , 式中x、y和z分别为径向、迹向和法向的实际相对位置, 和 分别为径向、迹向和法向的实际相对速度;名义相对运动状态X2d , 式中xd 、yd 和zd 分别为径向、迹向和法向的名义相对位置, 和 分别为径向、迹向和法向的名义相对速度;实际控制量U2 =[Ux Uz ]T ,其中Ux 和Uz 分别为径向和法向控制加速度;步骤S400:计算得到具体问题的控制量U2 ,将所得U2 代入公式(12)中,判断所得各项性能参数是否满足预设的性能指标,如果判断为满足则结束控制;如果判断为不满足则调整U2 中的各控制参数直至判断结果为满足所述性能指标时停止;其中,步骤S100包括以下步骤:建立迹向欠驱动航天器编队动力学模型:所述航天器包括主航天器和从航天器,所述欠驱动航天器编队动力学模型的坐标系定义:OE XI YI ZI 为地心惯性坐标系,其中OE 为地心,OC xyz为相对运动坐标系,其中OC 为主航天器质心,x轴沿主航天器径向,z轴与主航天器轨道面法向重合,y轴与x、z轴构成右手笛卡尔直角坐标系,OD 为从航天器质心,RC 和RD 分别为所述主航天器与所述从航天器的地心距矢量,所述迹向欠驱动航天器编队动力学模型在相对运动坐标系中的描述为 其中, 为外部摄动矢量,sat(U2 )为考虑输入饱和作用后由实际控制执行机构输入的控制加速度,F2 =[01×3 fy fx fz ]T (3) B=[02×4 I2×2 ]T (5)U2 =[Ux Uz ]T (6)其中,uC 为所述主航天器的纬度幅角, 和 分别为所述主航天器轨道角速度和轨道角加速度, 且 其中,μ=3.986×1014 m3 /s2 为地球引力常数,RC 和RD =[(RC +x)2 +y2 +z2 ]1/2 分别为主航天器和从航天器的地心距;0m×n 和Im×n 分别表示维数为m×n的零矩阵或单位矩阵,其中,Ux 和Uz 分别为径向和法向控制加速度;所述步骤S300包括以下步骤:建立误差动力学模型与构建控制律,得到缺失迹向控制加速度时,考虑外部摄动及模型线性化误差的欠驱动编队动力学模型中的可控部分为 其中 式中, 且 d2u =[dy 0 0]T 和d2a =[dx dz ]T 为不确定扰动矢量,sat(U2 )的表达式为sat(U2 )=[sat(Ux )sat(Uz )]T ,其中sat(·)为符号函数,即 式中,Ujm (j=x,z)为j方向可提供的最大控制加速度,sgn(·)为符号函数,其定义式为 构建的控制律为 其中 式中, 和 均为正定对角参数矩阵, 的定义式为 其中P21 为定常参数矩阵,其表达式为 式中,p21 <0为控制器参数,P22 =P21 A21 ,E21 =diag(ε211 ,ε212 )为正定增益矩阵,其中ε211 >ξ21m 且ε212 >ξ21m ,ξ21m 为矢量ξ21 =P21 d2u 的上界,即||ξ21 ||≤ξ21m =||P21 ||dm ,其中||P21 ||为矩阵P21 的诱导范数,同理,E22 =diag(ε221 ,ε222 )为正定增益矩阵,其中ε221 >ξ22m 且ε222 >ξ22m ,ξ22m 为矢量ξ22 =d2a +(K21 P21 +P22 )d2u 的上界,即ξ22m =(1+||K21 ||||P21 ||+||P22 ||)dm ,其中||K21 ||和||P22 ||分别为矩阵K21 和P22 的诱导范数,dm >0为矢量D2 的上界,D2 =[01×3 dy dx dz ]T ,sat(η22 ,δ21 ,δ22 )=[sat(η221 ,δ21 )sat(η222 ,δ22 )]T ,其中 式中,δ21 >0和δ22 >0为边界层的厚度;λ21 和λ22 的值由如下辅助系统积分得到,即 式中,ΔU2 =sat(U2 )-U2 。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于中国人民解放军国防科学技术大学,未经中国人民解放军国防科学技术大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/patent/201610186011.2/,转载请声明来源钻瓜专利网。