一种航天器高精度控制的设计方法技术

技术编号:29790690 阅读:21 留言:0更新日期:2021-08-24 18:09
本发明专利技术涉及一种航天器高精度控制设计方法,属于航天器轨道控制领域;步骤一、根据牛顿运动理论建立目标航天器的非线性相对运动模型;步骤二、考虑系统的不确定特性,通过线性化得到动态特性方程;步骤三、将动态特性方程分解,获得x,y,z三个方向的动态特性;步骤四、根据z方向的动态特性,设计z方向的一阶特性的虚拟控制信号

【技术实现步骤摘要】
一种航天器高精度控制的设计方法
本专利技术属于航天器交会控制领域,涉及一种航天器高精度控制设计方法。
技术介绍
随着我国政治、经济、科技等的快速发展,我国综合国力不断增强,航天事业也取得了一个又一个进步。从2000年开始,我国的航天事业进入了一个快速发展的重要阶段,我国的航天事业已经走在了世界的前列。目前航天器信息预测方法通常基于传统以李雅普诺夫稳定性理论所构建,进而讨论被控对象在无穷时间条件下的渐进稳定性情况。然而,若航天器在连续时间域内每一特定时刻的指向任务要求,即使其在有限时间后最终进入稳态,但在由初始时刻开始的稳定调节时间段内依然有可能由于所受扰动和不确定特性的影响,造成系统失真,降低系统运动平顺性,严重影响控制性能。
技术实现思路
本专利技术解决的技术问题是:克服现有技术的不足,提出一种航天器高精度控制设计方法,实现在任意时刻内,轨道控制精度均能满足预先的设计要求。。本专利技术解决技术的方案是:一种航天器高精度控制设计方法,包括如下步骤:步骤一、根据牛顿运动理论建立目标航天器的非线性相对运动模型;步骤二、考虑系统的不确定特性,通过线性化得到动态特性方程;步骤三、将动态特性方程分解,获得x,y,z三个方向的动态特性;步骤四、根据z方向的动态特性,设计z方向的一阶特性的虚拟控制信号步骤五、基于z方向的一阶特性的虚拟控制信号设计z方向的二阶真实控制信号uz(t);步骤六、重复步骤四至步骤五,依次设计y方向的二阶真实控制信号uy(t)和x方向的二阶真实控制信号ux(t);完成针对航天器相对轨道特性的高精度控制率设计。在上述的一种航天器高精度控制设计方法,所述步骤一中,建立非线性相对运动模型的具体方法为:目标航天器运行在圆形轨道上,建立目标航天器旋转坐标系oxyz;其中,原点o位于目标航天器的质心;x方向为从原点o指向地心方向;y方向为目标航天器轨道的切线方向;z方向由右手定则确定;目标航天器旋转坐标系oxyz以地球质心o′为圆心,在预设轨道半径R,预设轨道角速度n条件下进行旋转;目标航天器的非线性相对运动模型为:式中,μ为地球引力常量;为对*一次求导;为对*二次求导;af=[ux,uy,uz]T,ux,uy,uz为目标航天器沿三个方向的推力加速度。在上述的一种航天器高精度控制设计方法,所述步骤二中,动态特性方程为:Y=CX式中,X为动态特性方程状态;Y为动态特性方程输出;*3×3为3×3的矩阵;I3×3为3×3的单位矩阵;n为轨道角速度;p为不确定性系数;B=[03×3,I3×3]T;C=[I3×3,03×3]。在上述的一种航天器高精度控制设计方法,所述步骤三中,获得x,y,z三个方向的动态特性的方法为:定义x1为x、x2为y1为y、y2为z1为z、z2为则x方向的动态特性为:式中,ux为目标航天器沿x方向的推力加速度;则y方向的动态特性为:式中,uy为目标航天器沿y方向的推力加速度;则z方向的动态特性为:式中,uz为目标航天器沿z方向的推力加速度。在上述的一种航天器高精度控制设计方法,所述步骤四中,一阶特性的虚拟控制信号的设计方法为:令z方向的给定信号为rz,且给定信号rz关于时间连续可微;令z2的合理规划为则误差ez2(t)的计算方法为:进一步得到:采用Funnel控制策略,定义z方向第一性能漏斗集合为其中的元素为(t,ez1(t)),且ez1(t)满足|ez1(t)|<ψz1(t);其中,ez1(t)=z1(t)-rz(t);设定λz1为待设计的漏斗参数,λz1>0;则第一漏斗函数ψz1(t)为:ψz1(t)=max{az1exp(-bz1t),λz1}式中,az1和bz1为待设计参数,且az1>0,bz1>0;基于z方向性能漏斗集合设计一阶特性的虚拟控制信号式中,kz1(t)为第一时变漏斗增益,kz1(t)=[ψz1(t)-|ez1(t)|]-1;为饱和函数,参数为z方向一阶特性的虚拟控制信号的最大容许值,记为Dz2(ψz1(t),ez1(t))。在上述的一种航天器高精度控制设计方法,所述步骤五中,基于一阶特性的虚拟控制信号设计二阶真实控制信号uz(t)的具体方法为:定义z方向第二性能漏斗集合为其中的元素为(t,ez2(t)),且ez2(t)满足|ez2(t)|<ψz2(t);其中,ez2(t)=z2(t)-Dz2(ψz1(t),ez1(t));设定λz1为待设计的漏斗参数,λz1>0;则第二漏斗函数ψz2(t)为:ψz2(t)=max{az2exp(-bz2t),λz2}式中,az2和bz2为待设计参数,且az2>0,bz2>0;则,二阶真实控制信号uz(t)为:式中,kz2(t)为z方向第二时变漏斗增益,kz2(t)=[ψz2(t)-|ez2(t)|]-1;为z方向第二饱和函数,参数为二阶真实控制信号uz(t)的最大容许值。在上述的一种航天器高精度控制设计方法,所述步骤六中,设计y方向的二阶真实控制信号uy(t)的方法为:根据步骤四获得y方向一阶特性的虚拟控制信号为Dy2(ψy1(t),ey1(t)),其中,[ay1,by1,λy1]为相对应的待设计参数;在确定虚拟信号Dy2(ψy1(t),ey1(t))后,根据步骤五设计二阶真实控制信号uy(t)为:式中,ky2(t)为y方向第二时变漏斗增益,ky2(t)=[ψy2(t)-|ey2(t)|]-1;为y方向第二饱和函数,设计x方向的二阶真实控制信号ux(t)的方法为:根据步骤四获得x方向一阶特性的虚拟控制信号为Dx2(ψx1(t),ex1(t)),其中,[ax1,bx1,λx1]为相对应的待设计参数;在确定虚拟信号Dx2(ψx1(t),ex1(t))后,根据步骤五设计二阶真实控制信号ux(t)为:式中,kx2(t)为x方向第二时变漏斗增益,kx2(t)=[ψx2(t)-|ex2(t)|]-1为x方向第二饱和函数,本专利技术与现有技术相比的有益效果是:(1)本专利技术采用的Funnel控制策略,是一种基于瞬时性能要求所设计的控制器设计方法,其核心特征为能够使得被控指标在任意时刻都能位于预先所设计的性能空间内,进而保证被控系统的瞬时动态性能;(2)本专利技术通过矩阵分解将所得总体动态特性分为x,y,z三个子系统。根据每个子系统所具有的积分链形式设计其一阶特性的本文档来自技高网...

【技术保护点】
1.一种航天器高精度控制设计方法,其特征在于:包括如下步骤:/n步骤一、根据牛顿运动理论建立目标航天器的非线性相对运动模型;/n步骤二、考虑系统的不确定特性,通过线性化得到动态特性方程;/n步骤三、将动态特性方程分解,获得x,y,z三个方向的动态特性;/n步骤四、根据z方向的动态特性,设计z方向的一阶特性的虚拟控制信号

【技术特征摘要】
1.一种航天器高精度控制设计方法,其特征在于:包括如下步骤:
步骤一、根据牛顿运动理论建立目标航天器的非线性相对运动模型;
步骤二、考虑系统的不确定特性,通过线性化得到动态特性方程;
步骤三、将动态特性方程分解,获得x,y,z三个方向的动态特性;
步骤四、根据z方向的动态特性,设计z方向的一阶特性的虚拟控制信号
步骤五、基于z方向的一阶特性的虚拟控制信号设计z方向的二阶真实控制信号uz(t);
步骤六、重复步骤四至步骤五,依次设计y方向的二阶真实控制信号uy(t)和x方向的二阶真实控制信号ux(t);完成针对航天器相对轨道特性的高精度控制率设计。


2.根据权利要求1所述的一种航天器高精度控制设计方法,其特征在于:所述步骤一中,建立非线性相对运动模型的具体方法为:
目标航天器运行在圆形轨道上,建立目标航天器旋转坐标系oxyz;其中,原点o位于目标航天器的质心;x方向为从原点o指向地心方向;y方向为目标航天器轨道的切线方向;z方向由右手定则确定;
目标航天器旋转坐标系oxyz以地球质心o′为圆心,在预设轨道半径R,预设轨道角速度n条件下进行旋转;
目标航天器的非线性相对运动模型为:



式中,μ为地球引力常量;

为对*一次求导;

为对*二次求导;
af=[ux,uy,uz]T,ux,uy,uz为目标航天器沿三个方向的推力加速度。


3.根据权利要求2所述的一种航天器高精度控制设计方法,其特征在于:所述步骤二中,动态特性方程为:



Y=CX
式中,X为动态特性方程状态;
Y为动态特性方程输出;



*3×3为3×3的矩阵;
I3×3为3×3的单位矩阵;






n为轨道角速度;
p为不确定性系数;
B=[03×3,I3×3]T;
C=[I3×3,03×3]。


4.根据权利要求3所述的一种航天器高精度控制设计方法,其特征在于:所述步骤三中,获得x,y,z三个方向的动态特性的方法为:
定义x1为x、x2为y1为y、y2为z1为z、z2为
则x方向的动态特性为:






式中,ux为目标航天器沿x方向的推力加速度;
则y方向的动态特性为:






式中,uy为目标航天器沿y方向的推力加速度;
则z方向的动态特性为:






式中,uz为目标航天器沿z方向的推力加速度。


5.根据权利要求4所述的一种航天器高精度控制设计方法,其特征在于:所述步骤四中,一阶特性的虚拟控制信号的设计方法为:
令z方向的给定信号为rz,且给定信号rz关于时间连续可微;令z2的合理规划为则误差ez2(t)的计算方法为:



进一步得到:...

【专利技术属性】
技术研发人员:黄元刘磊严晗郝仁剑魏春岭
申请(专利权)人:北京控制工程研究所
类型:发明
国别省市:北京;11

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1