CN107861386B - 一种基于角速度观测器的抗干扰姿态控制地面验证系统及其控制方法 - Google Patents
一种基于角速度观测器的抗干扰姿态控制地面验证系统及其控制方法 Download PDFInfo
- Publication number
- CN107861386B CN107861386B CN201711094193.1A CN201711094193A CN107861386B CN 107861386 B CN107861386 B CN 107861386B CN 201711094193 A CN201711094193 A CN 201711094193A CN 107861386 B CN107861386 B CN 107861386B
- Authority
- CN
- China
- Prior art keywords
- attitude
- angular velocity
- control
- platform
- aircraft
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 54
- 238000004088 simulation Methods 0.000 claims abstract description 11
- 230000007246 mechanism Effects 0.000 claims description 27
- 238000012795 verification Methods 0.000 claims description 20
- 238000012544 monitoring process Methods 0.000 claims description 11
- 238000012360 testing method Methods 0.000 claims description 11
- 230000000694 effects Effects 0.000 claims description 9
- 230000014509 gene expression Effects 0.000 claims description 9
- 238000012546 transfer Methods 0.000 claims description 9
- 238000013461 design Methods 0.000 claims description 8
- 238000013329 compounding Methods 0.000 claims description 7
- 238000005457 optimization Methods 0.000 claims description 7
- 230000008859 change Effects 0.000 claims description 6
- 238000011161 development Methods 0.000 claims description 6
- 239000002131 composite material Substances 0.000 claims description 5
- 230000010354 integration Effects 0.000 claims description 5
- 206010034719 Personality change Diseases 0.000 claims description 3
- 230000005540 biological transmission Effects 0.000 claims description 3
- 230000007613 environmental effect Effects 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 238000010200 validation analysis Methods 0.000 abstract 3
- 238000004458 analytical method Methods 0.000 abstract 1
- 230000000052 comparative effect Effects 0.000 abstract 1
- 230000003245 working effect Effects 0.000 abstract 1
- 238000005259 measurement Methods 0.000 description 8
- 238000011160 research Methods 0.000 description 3
- 108091092878 Microsatellite Proteins 0.000 description 2
- 206010044565 Tremor Diseases 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 239000000428 dust Substances 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明涉及一种基于角速度观测器的抗干扰姿态控制地面验证系统及其控制方法,地面验证系统包括台下主控模块、台载控制计算机、姿态确定模块、执行机构模块以及三轴气浮台;所述地面验证系统是一类通用化验证系统,姿态确定模块选择该模块内不同敏感器组合验证敏感器失效等多种工况,通过台下主控模块选择台载控制计算机中姿态控制算法,完成不同姿态控制算法对比分析。所述基于角速度观测器的抗干扰姿态控制方法运行于台载控制计算机,能解决角速度信息缺失、载荷扰动情况下姿态控制系统可靠性差、抗干扰能力弱的问题。本发明的地面验证系统验证了基于角速度观测器的抗干扰姿态控制方法的有效性,提高了控制精度,适用于飞行器控制方法的地面仿真验证。
Description
技术领域
本发明涉及一种基于角速度观测器的抗干扰姿态控制方法及地面验证系统,可用于验证包括基于角速度观测器的抗干扰姿态控制方法在内的多种姿态控制方法;针对飞行器角速度信息的缺失和载荷扰动,提出了一种基于角速度观测器的抗干姿态扰控制方法,该方法能显著提高系统实时性、精度和稳定度,进一步改善姿控系统的控制性能,本发明属于飞行器的姿态控制领域。
背景技术
航天器长期在轨运行期间,控制回路中需要引入姿态和角速度测姿态测量信息进行反馈控制。然而航天工程中并不能保证所有的状态信息都一直高精度可测。角速率陀螺出故障或者测量精度大幅下降都会限制基于姿态和角速度测量的全状态反馈控制器的运用。因此,研究无角速度测量状况下的航天器姿态控制符合航天工程提高系统可靠性和容错性的需求,高效、实用的无角速度测量信息的姿态控制器设计越来越受到人们的关注。
在轨航天器含有多源干扰,既包括太阳光压、大气阻力、空间尘埃等外部环境干扰,卫星本身又有载荷转动、帆板振动、执行机构误差、敏感器测量噪声等内部扰动。多源干扰严重影响航天器的控制精度,尤其是在轨航天器硬件固定的情况下,难以在硬件上进一步挖掘控制精度提升的空间,因此对抗干扰姿态控制方法的研究及应用成为提高控制精度的重要新途径。随着科学技术发展和工程实际需要,航天器挠性化程度也越来越高。航天器有效载荷的扰动(如CCD立体相机的震颤、扫描镜转动等)对姿态带来的干扰是姿态控制的干扰源之一。经文献检索,王新升、韩建斌、梁斌在论文基于气浮台的微小卫星姿态控制实时仿真(见《北京航空航天大学学报》,2010年,第7期,第36卷,页码767-770)搭建了半物理仿真实验平台对单刚体微小卫星的姿态控制问题进行了实时仿真研究,其控制算法为传统PID控制方法,该方法不能有效抵消载荷扰动,未能达到精细抗干扰的目的,也无法验证角速度信息缺失情况下姿态控制算法的有效性。陆智俊、吴敬玉在论文基于MRP的航天器无角速度姿态控制算法(见《上海航天》,2013年,第2期,第30卷,页码13-21)用一阶无源滤波器对姿态参数进行差分,提出了基于修正罗德里格参数(MRP)的采用饱和函数和双曲余弦函数的无角速度控制律,实现了无角速度反馈情况下姿态控制,但不涉及地面验证系统。
姿控系统测试装置已经在飞行器研制过程中得到了广泛应用,国内航天五院502所、航天八院812所,还有哈工大、清华、北航等高校都曾研制出姿控系统测试装置,但现有的姿控系统测试装置通常针对特定型号的飞行器而研制,有的仅仅针对某特定航天任务,为验证某种特定方法搭建,忽略了载荷扰动的影响,缺乏敏感器信息缺失等不同工况下姿态控制算法验证,无法用于针对包括基于飞角速度观测器在内的多种抗干扰姿态控制方法研究。
发明内容
本发明的技术解决问题是:克服现有技术的不足,提供一种基于角速度观测器的抗干扰姿态控制地面验证系统及其控制方法,改善无角速度量测信息下航天器姿控系统的控制性能,提高姿态控制系统的响应速度,减小稳态偏差,同时验证了基于角速度观测器的抗干扰姿态控制方法的有效性。
本发明的技术解决方案是:一种基于角速度观测器的抗干扰姿态控制地面验证系统,所述地面验证系统包括:台下主控模块、台载控制计算机、姿态确定模块、执行机构模块以及三轴气浮台;所述姿态确定模块包括星敏感器、太阳敏感器、陀螺仪以及失效决策单元、无线发送单元;姿态确定模块中的星敏感器、太阳敏感器、陀螺仪用于实时获取飞行器三轴角度和三轴角速度信息,姿态确定模块中的失效决策单元可通过RS422总线向陀螺仪发送失效指令,同时通过无线发送单元向台载控制计算机发送敏感器失效指令;所述台载控制计算机包括姿态控制单元和无线接收单元,无线接收单元接收失效决策单元的敏感器失效指令以及台下主控模块的姿态控制方法的切换指令,姿态控制单元根据该切换指令完成姿态控制算法的选择与运行,姿态控制算法包括基于角速度观测器的抗干扰姿态控制方法、PID控制方法或鲁棒控制方法,为执行机构模块提供力矩控制指令;所述执行机构模块包括冷气推力器、力矩陀螺群、反作用飞轮组、偏置动量轮以及选择调度单元;执行机构模块的选择调度单元根据任务需求选择具体的执行机构类型;所述执行机构模块在接收台载控制计算机提供的力矩控制指令后,根据力矩控制指令输出力矩信号,输出力矩信号传给台载控制计算机;所述台下主控模块包括数据监测单元、无线发送单元、试验主控单元、集成优化单元;台下主控模块中的试验主控单元通过无线发送单元向台载控制计算机发送切换指令信号;数据监测单元存储不同姿态控制方法下的仿真实时运算数据,用于对比分析不同姿态控制方法下的控制效果;集成优化单元实现与ANSYS、STK、MATLAB的集成与二次开发;三轴气浮台作为仿真的支撑平台,姿态确定模块、台载控制计算机以及执行机构模块安装在气浮台面上,三轴气浮台的转动用来模拟飞行器在外层空间的姿态变化;地面验证系统中飞行器的期望姿态通过台载控制计算机输入,姿态确定模块中的失效决策单元向陀螺仪、台载控制计算机发送敏感器失效指令,姿态确定模块给出不含角速度信息的飞行器实际姿态,台载控制计算机接收敏感器失效指令后选择基于角速度观测器的抗干扰姿态控制方法,姿态控制算法解算出指令力矩信号并传送至选择调度单元确定的具体执行机构;执行机构在接收指令力矩信号后输出执行力矩信号,该力矩作用在三轴气浮台的台面引起台面转动,三轴气浮台模拟飞行器在外层空间的力学环境,姿态确定模块中的敏感器获取得到包括飞行器三轴转动角度和三轴转动角速度的姿态信息,飞行器姿态信息的实时数据传输给台下主控模块中的数据监测单元,该单元保存实时运算数据,同时该姿态信息传送至姿态确定模块;姿态确定模块得到的姿态信息与期望姿态作比较后获得新的偏差信号,形成了验证的数据流回路。
所述的基于角速度观测器的抗干扰姿态控制方法,包括以下步骤:
第一步,建立飞行器系统动力学模型;
其中,
I为飞行器惯量矩阵,q为飞行器姿态四元数,为飞行器姿态四元数的一阶导数,ω为飞行器角速度,u为复合控制器的输出,Td为干扰力矩,包括环境力矩和载荷扰动;ω1、ω2、ω3分别表示角速度ω在x、y、z轴上的分量,q0是姿态四元数标量部分;
第二步,针对飞行器姿态控制系统中的敏感器失效带来的角速度信息缺失,设计角速度观测器;
角速度观测器为:
其中,γ、λ为待设计参数,为角速度估计值的变化率,为角速度估计值,为四元数估计值的变化率,为四元数估计值,ω为角速度,u为复合控制器的输出,的表达式由下式给出:
其中,分别表示角速度估计值在x、y、z轴上的分量,的表达式由下式给出:
其中,q0是姿态四元数标量部分,qv=(q1,q2,q3)是姿态四元数矢量部分,为q0的估计值,为qv的估计值;的表达式分别由下式给出:
其中,I3×3为三阶单位阵;
第三步,针对飞行器系统中存在的载荷扰动设计载荷扰动估计器及PID控制器;
载荷扰动估计器为:
即采用Q(s)d(s)对载荷扰动d(s)进行估计,Q(s)与Q(s)G-1(s)构成了载荷扰动估计器;为载荷扰动的估计值,Y(s)为飞行器系统输出,表示为:
Y(s)=Guyu(s)+Gdyd(s);
其中,u(s)为控制输入,d(s)为载荷转动干扰,Guy为从输入到输出的闭环传递函数,Gdy为从干扰到输出的闭环传递函数,G(s)为飞行器系统模型,G0(s)为飞行器系统标称模型,Q(s)为滤波器,Ed(s)为扰动估计误差,表示为上述各式中的s代表载荷扰动估计器基于频域设计;
载荷扰动类型为慢时变低频扰动,Q(s)设计为低通滤波器,即干扰估计的效果由滤波器Q(s)的设计决定;为了达到最优扰动估计效果,使Q(s)接近1,即Ed(s)接近0,达到抵消干扰的效果;
PID控制器为:Gc(s)为PID控制器中实现的传递函数;
采用PID控制方法进行反馈补偿,PID控制律为:
△m=min-mout
其中,Kp、Ki、Kd分别为比例增益、积分增益、微分增益;Tc为PID控制器的输出,△m为姿态角偏差,△ω为姿态角速度偏差,min为期望姿态角,mout为输出姿态角,ωin为期望角速度,为角速度估计值;
第四步,将角速度观测器、载荷扰动观测器和PID控制器进行复合,实现基于角速度观测器的抗干扰姿态控制;角速度观测器、载荷扰动估计器和PID控制器进行复合如下: 为扰动的估计值,即载荷扰动估计器的输出,u(s)是角速度观测器、载荷扰动估计器和PID控制器复合后得到的控制输入。
本发明与现有技术相比的优点在于:
(1)基于角速度观测器的抗干扰控制方法在工程上易于实现,把基于干扰观测器控制的干扰抵消方法与PID控制方法有机结合,能够实现在角速度信息缺失情况下飞行器姿态高精度控制,不仅能显著提高系统可靠性,也能够减小稳态偏差,进一步改善姿控系统的控制性能,提高了姿态控制系统的控制精度。
(2)通过基于角速度观测器的抗干扰地面验证系统,提出了基于角速度观测器的抗干扰控制方法,改善了传统的单一鲁棒方法对于干扰抑制和抵消问题保守性大的缺陷。本发明的抗干扰地面验证系统通用性强、可扩展性高,不仅能够研究基于载体扰动的抗干扰姿态控制算法,还能够根据具体任务类型,通过选择调度单元选择不同类型执行机构、不同敏感器组合完成姿态控制算法验证,通过试验主控模块切换不同控制方法,姿态控制单元逐次验证多种控制方法。
附图说明
图1为本发明一种基于角速度观测器的抗干扰姿态控制地面验证系统的数据流回路;
图2为本发明一种基于角速度观测器的抗干扰姿态控制方法设计流程图。
具体实施方式
下面结合附图对本发明的具体实施方式做进一步详细说明。
以微纳三轴稳定卫星的姿控系统仿真试验为例,来说明地面验证系统以及基于角速度观测器的抗干扰姿态控制方法。
如图1所示,给出了本发明所述的地面验证系统,验证系统包括执行机构模块1、三轴气浮台2、姿态确定模块3、台载控制计算机4、台下主控模块5;所述姿态确定模块3包括太阳敏感器33、星敏感器34、陀螺仪35以及失效决策单元31、无线发送单元32;姿态确定模块3中的姿态敏感器实时获取微纳卫星三轴角度和三轴角速度信息,姿态确定模块3中的失效决策单元31可通过RS422总线向姿态敏感器发送失效指令,同时通过无线发送单元32向台载控制计算机4发送敏感器失效指令;所述台载控制计算机4包括姿态控制单元41和无线接收单元42,无线接收单元42接收失效决策单元的敏感器失效指令以及台下主控模块5的姿态控制方法的切换指令,姿态控制单元41根据该切换指令完成姿态控制算法的选择与运行,姿态控制算法包括基于角速度观测器的抗干扰姿态控制方法、PID控制方法或鲁棒控制方法,为执行机构模块1提供力矩控制指令;所述执行机构模块1包括选择调度单元11、冷气推力器12、力矩陀螺群13、反作用飞轮组14以及偏置动量轮15;执行机构模块1的选择调度单元11根据任务需求选择具体的执行机构类型;所述执行机构模块1在接收台载控制计算机4提供的力矩控制指令后,将输出力矩信号传给台载控制计算机4;所述台下主控模块5包括试验主控单元51、无线发送单元52、数据监测单元53、集成优化单元54;台下主控模块5中的试验主控单元51通过无线发送单元52向台载控制计算机4发送切换指令信号;数据监测单元53存储不同姿态控制方法下的仿真实时运算数据,用于对比分析不同姿态控制方法下的控制效果;集成优化单元54实现与ANSYS、STK、MATLAB的集成与二次开发;三轴气浮台2作为仿真的支撑平台,姿态确定模块3、台载控制计算机4以及执行机构模块1安装在气浮台面上,三轴气浮台面的转动用来模拟微纳卫星在外层空间的姿态变化;地面验证系统的数据流回路如下:微纳卫星的期望姿态通过台载控制计算机4输入,姿态确定模块3中的失效决策单元31向姿态敏感器、台载控制计算机4发送敏感器失效指令,姿态确定模块3给出不含角速度信息的微纳卫星实际姿态,台载控制计算机4接收敏感器失效指令后选择基于角速度观测器的抗干扰姿态控制方法,姿态控制算法解算出指令力矩信号并传送至选择调度单元11确定的具体执行机构;执行机构在接收指令力矩信号后输出执行力矩信号,该力矩作用在三轴气浮台2的台面引起台面转动,三轴气浮台2模拟微纳卫星在外层空间的力学环境,姿态确定模块3中的敏感器(太阳敏感器33、星敏感器、陀螺仪35)获取得到包括微纳卫星三轴转动角度和三轴转动角速度的姿态信息,微纳卫星姿态信息的实时数据传输给台下主控模块5中的数据监测单元53,该单元保存实时运算数据,集成优化单元54利用姿态信息可实现与ANSYS、STK、MATLAB的集成与二次开发,同时该姿态信息传送至姿态确定模块3;姿态确定模块3得到的姿态信息与期望姿态作比较后获得新的偏差信号,形成了验证的数据流回路。
如图2所示,基于角速度观测器的抗干扰姿态控制方法,包括以下步骤:
首先,建立微纳卫星的动力学及运动学模型;其次,针对微纳卫星姿态控制系统中的敏感器失效带来的角速度信息缺失,设计角速度观测器;再次,针对微纳卫星系统存在的载荷扰动设计载荷扰动估计器及PID控制器;最后,将角速度观测器、载荷扰动估计器和PID控制器进行复合,给出基于角速度观测器的抗干扰姿态控制方法;具体步骤如下:
1,建立微纳卫星动力学及运动学模型如下:
其中,
I为微纳卫星惯量矩阵,I=diag{5.50,6.14,2.28},ω=[ω1,ω2,ω3]为微纳卫星角速度,为角速度的导数,q为微纳卫星姿态四元数,为微纳卫星姿态四元数的导数,u为复合控制器的输出,Td为干扰力矩,包括环境力矩和载荷扰动;ω1、ω2、ω3分别表示角速度ω在x、y、z轴上的分量,q0是姿态四元数标量部分;微纳卫星初始姿态四元数为[0.9993,0.0246,0.0156,0.0254]T,初始角速度为[0.001,0.001,0.0015]rad/s。
2,针对微纳卫星姿态控制系统中的敏感器失效带来的角速度信息缺失,设计角速度观测器;
角速度观测器为:
其中,γ、λ为待设计参数,为角速度估计值的变化率,为角速度估计值,为四元数估计值的变化率,为四元数估计值,u为复合控制器的输出,的表达式由下式给出:
其中,分别表示角速度估计值在x、y、z轴上的分量,的表达式由下式给出:
其中,q0是姿态四元数标量部分,qv=(q1,q2,q3)是姿态四元数矢量部分,为q0的估计值,为qv的估计值;的表达式分别由下式给出:
其中,I3×3为三阶单位阵;
3,针对微纳卫星系统中存在的载荷扰动设计载荷扰动估计器及PID控制器;
载荷扰动估计器为:
即采用Q(s)d(s)对载荷扰动d(s)进行估计,Q(s)与Q(s)G-1(s)构成了载荷扰动估计器;为载荷扰动的估计值,Y(s)为微纳卫星系统输出,表示为:
Y(s)=Guyu(s)+Gdyd(s);
其中,u(s)为控制输入,d(s)为载荷转动干扰,Guy为从输入到输出的闭环传递函数,Gdy为从扰动到输出的闭环传递函数,G(s)为微纳卫星系统模型,G0(s)为微纳卫星系统标称模型,Q(s)为滤波器,Ed(s)为扰动估计误差,表示为上述各式中的s代表载荷扰动估计器基于频域设计;
载荷扰动类型为慢时变低频扰动,Q(s)设计为低通滤波器,即扰动估计的效果由滤波器Q(s)的设计决定;为了达到最优扰动估计效果,使Q(s)接近1,即Ed(s)接近0,达到抵消扰动的效果;
PID控制器为:Gc(s)为PID控制器中实现的传递函数;
采用PID控制方法进行反馈补偿,PID控制律为:
△m=min-mout
其中,Kp、Ki、Kd分别为比例增益、积分增益、微分增益,比例、积分、微分增益最优值分别是:Kp=[-5.6,-12.28,-2.18],Ki=[-0.55,-0.614,-0.218],Kd=[-12.44,-13.89,-4.93];Tc为PID控制器的输出,△m为姿态角偏差,△ω为姿态角速度偏差,min为期望姿态角,mout为输出姿态角,ωin为期望角速度,为角速度估计值;
4,将角速度观测器、载荷扰动估计器和PID控制器进行复合,给出基于角速度观测器的抗干扰姿态控制方法;
基于角速度观测器的抗干扰姿态控制方法为:为扰动的估计值,即载荷扰动估计器的输出,u(s)是角速度观测器、载荷扰动估计器和PID控制器复合后得到的控制输入。
本发明说明书中未作详细描述的内容属于本领域专业技术人员公知的现有技术。
Claims (2)
1.一种基于角速度观测器的抗干扰姿态控制地面验证系统,其特征在于:所述地面验证系统包括:台下主控模块、台载控制计算机、姿态确定模块、执行机构模块以及三轴气浮台;所述姿态确定模块包括星敏感器、太阳敏感器、陀螺仪以及失效决策单元、无线发送单元;姿态确定模块中的星敏感器、太阳敏感器、陀螺仪用于实时获取飞行器三轴角度和三轴角速度信息,姿态确定模块中的失效决策单元可通过RS422总线向陀螺仪发送失效指令,同时通过无线发送单元向台载控制计算机发送敏感器失效指令;所述台载控制计算机包括姿态控制单元和无线接收单元,无线接收单元接收失效决策单元的敏感器失效指令以及台下主控模块的姿态控制方法的切换指令,姿态控制单元根据该切换指令完成姿态控制算法的选择与运行,姿态控制算法包括基于角速度观测器的抗干扰姿态控制方法、PID控制方法或鲁棒控制方法,为执行机构模块提供力矩控制指令;所述执行机构模块包括冷气推力器、力矩陀螺群、反作用飞轮组、偏置动量轮以及选择调度单元;执行机构模块的选择调度单元根据任务需求选择具体的执行机构类型;所述执行机构模块在接收台载控制计算机提供的力矩控制指令后,根据力矩控制指令输出力矩信号,输出力矩信号传给台载控制计算机;所述台下主控模块包括数据监测单元、无线发送单元、试验主控单元、集成优化单元;台下主控模块中的试验主控单元通过无线发送单元向台载控制计算机发送切换指令信号;数据监测单元存储不同姿态控制方法下的仿真实时运算数据,用于对比分析不同姿态控制方法下的控制效果;集成优化单元实现与ANSYS、STK、MATLAB的集成与二次开发;三轴气浮台作为仿真的支撑平台,姿态确定模块、台载控制计算机以及执行机构模块安装在气浮台面上,三轴气浮台的转动用来模拟飞行器在外层空间的姿态变化;地面验证系统中飞行器的期望姿态通过台载控制计算机输入,姿态确定模块中的失效决策单元向陀螺仪、台载控制计算机发送敏感器失效指令,姿态确定模块给出不含角速度信息的飞行器实际姿态,台载控制计算机接收敏感器失效指令后选择基于角速度观测器的抗干扰姿态控制方法,姿态控制算法解算出指令力矩信号并传送至选择调度单元确定的具体执行机构;执行机构在接收指令力矩信号后输出执行力矩信号,该力矩作用在三轴气浮台的台面引起台面转动,三轴气浮台模拟飞行器在外层空间的力学环境,姿态确定模块中的敏感器获取得到包括飞行器三轴转动角度和三轴转动角速度的姿态信息,飞行器姿态信息的实时数据传输给台下主控模块中的数据监测单元,该单元保存实时运算数据,同时该姿态信息传送至姿态确定模块;姿态确定模块得到的姿态信息与期望姿态作比较后获得新的偏差信号,形成了验证的数据流回路。
2.根据权利要求1所述的基于角速度观测器的抗干扰姿态控制地面验证系统的抗干扰姿态控制方法,其特征在于:包括以下步骤:
第一步,建立飞行器系统动力学模型;
其中,
I为飞行器惯量矩阵,q为飞行器姿态四元数,为飞行器姿态四元数的一阶导数,ω为飞行器角速度,u为复合控制器的输出,Td为干扰力矩,包括环境力矩和载荷扰动;ω1、ω2、ω3分别表示角速度ω在x、y、z轴上的分量,q0是姿态四元数标量部分;
第二步,针对飞行器姿态控制系统中的敏感器失效带来的角速度信息缺失,设计角速度观测器;
角速度观测器为:
其中,γ、λ为待设计参数,为角速度估计值的变化率,为角速度估计值,为四元数估计值的变化率,为四元数估计值,ω为角速度,的表达式由下式给出:
其中,分别表示角速度估计值在x、y、z轴上的分量,的表达式由下式给出:
其中,q0是姿态四元数标量部分,qv=(q1,q2,q3)是姿态四元数矢量部分,为q0的估计值,为qv的估计值;的表达式分别由下式给出:
其中,I3×3为三阶单位阵;
第三步,针对飞行器系统中存在的载荷扰动设计载荷扰动估计器及PID控制器;
载荷扰动估计器为:
即采用Q(s)d(s)对载荷扰动d(s)进行估计,Q(s)与Q(s)G-1(s)构成了载荷扰动估计器;为载荷扰动的估计值,Y(s)为飞行器系统输出,表示为:
Y(s)=Guyu(s)+Gdyd(s);
其中,u(s)为控制输入,d(s)为载荷转动干扰,Guy为从输入到输出的闭环传递函数,Gdy为从干扰到输出的闭环传递函数,G(s)为飞行器系统模型,G0(s)为飞行器系统标称模型,Q(s)为滤波器,Ed(s)为扰动估计误差,表示为上述各式中的s代表载荷扰动估计器基于频域设计;
载荷扰动类型为慢时变低频扰动,Q(s)设计为低通滤波器,即干扰估计的效果由滤波器Q(s)的设计决定;为了达到最优扰动估计效果,使Q(s)接近1,即Ed(s)接近0,达到抵消干扰的效果;
PID控制器为:Gc(s)为PID控制器中实现的传递函数;
采用PID控制方法进行反馈补偿,PID控制律为:
Δm=min-mout
其中,Kp、Ki、Kd分别为比例增益、积分增益、微分增益;Tc为PID控制器的输出,Δm为姿态角偏差,Δω为姿态角速度偏差,min为期望姿态角,mout为输出姿态角,ωin为期望角速度,为角速度估计值;
第四步,将角速度观测器、载荷扰动观测器和PID控制器进行复合,实现基于角速度观测器的抗干扰姿态控制;角速度观测器、载荷扰动估计器和PID控制器进行复合如下: 为扰动的估计值,即载荷扰动估计器的输出,u(s)是角速度观测器、载荷扰动估计器和PID控制器复合后得到的控制输入。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711094193.1A CN107861386B (zh) | 2017-11-08 | 2017-11-08 | 一种基于角速度观测器的抗干扰姿态控制地面验证系统及其控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711094193.1A CN107861386B (zh) | 2017-11-08 | 2017-11-08 | 一种基于角速度观测器的抗干扰姿态控制地面验证系统及其控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107861386A CN107861386A (zh) | 2018-03-30 |
CN107861386B true CN107861386B (zh) | 2019-03-19 |
Family
ID=61701333
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711094193.1A Active CN107861386B (zh) | 2017-11-08 | 2017-11-08 | 一种基于角速度观测器的抗干扰姿态控制地面验证系统及其控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107861386B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110231831B (zh) * | 2019-07-23 | 2022-01-28 | 兰州理工大学 | 一种基于角度测量的航天器姿态解耦最优控制方法 |
CN111102887B (zh) * | 2019-12-31 | 2022-04-08 | 三峡大学 | 一种基于反作用飞轮的弹体滚转角控制方法 |
CN111158386B (zh) * | 2020-01-14 | 2021-02-02 | 北京航空航天大学 | 一种干扰自主免疫的无人机姿态控制系统及实现方法 |
CN112947116B (zh) * | 2021-02-09 | 2022-07-26 | 南京航空航天大学 | 一种三轴气浮台的联合执行机构控制方法 |
CN113401366B (zh) * | 2021-06-17 | 2022-07-05 | 北京控制工程研究所 | 一种克服周期性活动部件影响的强抗扰复合控制方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104932510A (zh) * | 2015-06-04 | 2015-09-23 | 北京航空航天大学 | 一种针对飞轮低速摩擦的挠性航天器姿态控制系统与方法 |
CN105807780A (zh) * | 2016-05-30 | 2016-07-27 | 北京航空航天大学 | 一种基于飞轮输出偏差的抗干扰姿态控制方法及验证装置 |
CN106020165A (zh) * | 2016-05-30 | 2016-10-12 | 北京航空航天大学 | 一种针对执行机构故障的航天器容错控制方法与验证装置 |
CN106017509A (zh) * | 2016-05-30 | 2016-10-12 | 北京航空航天大学 | 一种多源干扰环境下抗干扰姿态确定方法与测试平台 |
CN106020221A (zh) * | 2016-05-30 | 2016-10-12 | 北京航空航天大学 | 一种基于输出反馈的抗干扰姿态控制验证平台及验证方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101008176B1 (ko) * | 2008-12-30 | 2011-01-13 | 한국항공우주연구원 | 반작용휠과 추력기 기반 자세제어기를 동시에 이용한 자세기동 및 가제어성 향상 방법 |
-
2017
- 2017-11-08 CN CN201711094193.1A patent/CN107861386B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104932510A (zh) * | 2015-06-04 | 2015-09-23 | 北京航空航天大学 | 一种针对飞轮低速摩擦的挠性航天器姿态控制系统与方法 |
CN105807780A (zh) * | 2016-05-30 | 2016-07-27 | 北京航空航天大学 | 一种基于飞轮输出偏差的抗干扰姿态控制方法及验证装置 |
CN106020165A (zh) * | 2016-05-30 | 2016-10-12 | 北京航空航天大学 | 一种针对执行机构故障的航天器容错控制方法与验证装置 |
CN106017509A (zh) * | 2016-05-30 | 2016-10-12 | 北京航空航天大学 | 一种多源干扰环境下抗干扰姿态确定方法与测试平台 |
CN106020221A (zh) * | 2016-05-30 | 2016-10-12 | 北京航空航天大学 | 一种基于输出反馈的抗干扰姿态控制验证平台及验证方法 |
Non-Patent Citations (2)
Title |
---|
Improved satellite attitude control using a disturbance compensator;Toshiaki Yamashita 等;《Acta Astronautica》;20041231;15-25 |
基于气浮台的小卫星姿态控制全物理仿真实验系统;高华宇 等;《计算机测量与控制》;20090630;第17卷(第06期);1153-1156 |
Also Published As
Publication number | Publication date |
---|---|
CN107861386A (zh) | 2018-03-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107861386B (zh) | 一种基于角速度观测器的抗干扰姿态控制地面验证系统及其控制方法 | |
CN106020221B (zh) | 一种基于输出反馈的抗干扰姿态控制验证平台及验证方法 | |
CN104932510B (zh) | 一种针对飞轮低速摩擦的挠性航天器姿态控制系统与方法 | |
Zhu et al. | Adaptive sliding mode disturbance observer-based composite control with prescribed performance of space manipulators for target capturing | |
Ding et al. | Nonsmooth attitude stabilization of a flexible spacecraft | |
CN105843244A (zh) | 一种基于输出反馈的挠性航天器精细姿态控制方法 | |
Wang et al. | Tracking control for a velocity-sensorless VTOL aircraft with delayed outputs | |
CN110316402A (zh) | 一种编队控制模式下的卫星姿态控制方法 | |
CN105807780A (zh) | 一种基于飞轮输出偏差的抗干扰姿态控制方法及验证装置 | |
Xiao et al. | A general tracking control framework for uncertain systems with exponential convergence performance | |
Srivastava et al. | Attitude determination and control system for a leo debris chaser small satellite | |
Zhang et al. | Adaptive disturbance observer-based dual-loop integral-type fast terminal sliding mode control for micro spacecraft and its gimbal tracking device | |
CN109445448B (zh) | 一种轮控小卫星自适应积分滑模姿态控制器 | |
CN111258325A (zh) | 对地遥感卫星高精度高性能的姿态容错控制方法、装置及计算机存储介质 | |
Zhang et al. | Finite-time attitude optimization maneuver control for coupled spacecraft under attitude measurement errors and actuator faults | |
Yanyachi et al. | Low-cost Test System for 1U CubeSat Attitude Control with Reaction Wheels | |
Liao et al. | Finite-time stabilization for a novel non-contact spacecraft with a small air-clearance range constraint | |
CN103885451B (zh) | 一种新型抗干扰姿态控制电路解算系统 | |
Newton | Design, Development, and Experimental Validation of a Nanosatellite Attitude Control Simulator | |
طائی et al. | Optimal Control of a Tri-axial Spacecraft Simulator Test bed Actuated by Reaction Wheels | |
Chávez-Moreno et al. | Contracting angular velocity observer for small satellites | |
Hong et al. | Application of EKF for missile attitude estimation based on “SINS/CNS” integrated guidance system | |
Xiao et al. | Large-angle velocity-free attitude tracking control of satellites: An observer-free framework | |
Xuan et al. | From PID to L1 adaptive control for automatic balancing of a spacecraft three-axis simulator | |
Haotian et al. | Accurate attitude estimation of HB2 standard model based on QNCF in hypersonic wind tunnel test |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right |
Effective date of registration: 20230713 Address after: Room 313-2851, Floor 3, Shanshui Building, No. 3 Gulou Street, Miyun District, Beijing 101500 (Central Office Area of Yunchuanggu Economic Development Center) Patentee after: Beijing Qixing Hangyu Technology Co.,Ltd. Address before: 100191 No. 37, Haidian District, Beijing, Xueyuan Road Patentee before: BEIHANG University |
|
TR01 | Transfer of patent right |