CN107943062A - 带外部干扰力矩的挠性卫星姿态滑模控制方法 - Google Patents
带外部干扰力矩的挠性卫星姿态滑模控制方法 Download PDFInfo
- Publication number
- CN107943062A CN107943062A CN201710901110.9A CN201710901110A CN107943062A CN 107943062 A CN107943062 A CN 107943062A CN 201710901110 A CN201710901110 A CN 201710901110A CN 107943062 A CN107943062 A CN 107943062A
- Authority
- CN
- China
- Prior art keywords
- flexible
- attitude
- satellite
- mode control
- sliding mode
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 claims abstract description 30
- 239000011159 matrix material Substances 0.000 claims description 15
- 239000000126 substance Substances 0.000 claims description 6
- 238000013016 damping Methods 0.000 claims description 5
- 230000001808 coupling effect Effects 0.000 claims description 3
- 230000001629 suppression Effects 0.000 abstract description 3
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 238000009415 formwork Methods 0.000 abstract 3
- 238000004422 calculation algorithm Methods 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
- B64G1/244—Spacecraft control systems
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Automation & Control Theory (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Radar, Positioning & Navigation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明提供了一种带外部干扰力矩的挠性卫星姿态滑模控制方法,利用姿态四元数建立挠性卫星运动学方程,采用混合坐标法对中心刚体带有挠性附件、外部干扰匹配的挠性卫星建立动力学方程;采用等效控制思想,基于Lyapunov直接法设计以下两种滑模控制律:针对挠性模态可测量的情况,设计一种基于状态反馈的滑模控制律;针对挠性模态不可测量的情况设计一种基于动态观测器的滑模控制律。本发明的有益效果是:解决了挠性卫星在飞行过程中存在外界干扰时的姿态控制及挠性附件的振动控制问题,具有较强的振动抑制能力,完成了挠性卫星姿态机动的高鲁棒性控制。
Description
技术领域
本发明涉及挠性卫星姿态滑模控制方法,尤其涉及一种带外部干扰力矩的挠性卫星姿态滑模控制方法。
背景技术
传统的挠性卫星姿态滑模控制方法在运行时存在大量计算,使得控制时间长,导致控制系统无法及时抑制干扰,难以保证良好的鲁棒性。
发明内容
为了解决现有技术中的问题,本发明提供了一种带外部干扰力矩的挠性卫星姿态滑模控制方法。
本发明提供了一种带外部干扰力矩的挠性卫星姿态滑模控制方法,利用姿态四元数建立挠性卫星运动学方程,采用混合坐标法对中心刚体带有挠性附件、外部干扰匹配的挠性卫星建立动力学方程;采用等效控制思想,基于Lyapunov直接法设计以下两种滑模控制律:针对挠性模态可测量的情况,设计一种基于状态反馈的挠性卫星姿态滑模控制律;针对挠性模态不可测量的情况设计一种基于动态观测器的挠性卫星姿态滑模控制律。
作为本发明的进一步改进,以姿态四元数为基础的挠性卫星运动学方程和动力学方程如下:
其中,
其中,Jmb=J-δTδ,
其中,q0,qv分别为姿态四元数的标量部分与向量部分;ω表示挠性卫星的姿态角速度;δ为表示卫星本体与挠性部件的耦合作用矩阵;C,K分别表示阻尼矩阵和刚度矩阵,
考虑N个弹性模态,其对应的自然角频率为ωni,i=1,2,…,N,对应的阻尼为i=1,2,…,N;Jmb为刚体部分的转动惯量,u表示控制力矩,d表示有界外部干扰力矩,卫星带挠性附件,并且含有有界外部干扰力矩。
作为本发明的进一步改进,针对挠性模态可测的情形下,设计如下滑模面:
s=kqv+ω.
设计的基于状态反馈的挠性卫星姿态滑模控制律如下:
F(s)={sgn(s1),sgn(s2),sgn(s3)}T
min{D1(i)}>||d||
其中,K1,D1,k为正定对角矩阵。
作为本发明的进一步改进,针对挠性模态不可测的情形下,设计如下滑模面:
s=kqv+ω
设计的基于动态观测器的挠性卫星姿态滑模控制律如下:
F(s)={sgn(s1),sgn(s2),sgn(s3)}T
min{D1(i)}>||d||
其中,K1,D1,k为正定对角矩阵。正定对称矩阵P满足下述Lyapunov方程
本发明的有益效果是:通过上述方案,解决了挠性卫星在飞行过程中存在外界干扰时的姿态控制及挠性附件的振动控制问题,具有较强的振动抑制能力,完成了挠性卫星姿态机动的高鲁棒性控制;采用本发明设计的姿态控制器保证卫星姿态的良好稳定性和模态振动的衰减,当卫星系统存在外部干扰力矩时,卫星的姿态能很快趋于稳定,能够有效消除外界干扰带来的影响。
附图说明
图1是本发明一种带外部干扰力矩的挠性卫星姿态滑模控制方法的基于状态反馈的姿态滑模控制算法的simulink模型示意图。
图2是发明一种带外部干扰力矩的挠性卫星姿态滑模控制方法的基于动态观测器的姿态滑模控制算法的simulink模型示意图。
具体实施方式
下面结合附图说明及具体实施方式对本发明作进一步说明。
如图1至图2所示,一种带外部干扰力矩的挠性卫星姿态滑模控制方法,包括:
以姿态四元数为基础对存在有界外部干扰力矩的挠性卫星建立如下运动学方程和动力学方程:
其中,Jmb=J-δTδ,
其中,d是有界的外部干扰力矩,可为常值干扰或者正弦干扰。其特征在于:卫星带有挠性附件,并且在运行过程中受各种干扰力矩的影响。
考虑以下两种情况:
(1)假设卫星的状态η和ψ完全可测量,针对挠性卫星(1)-(2)设计基于状态反馈的滑模控制律。
Step1 选取滑模面
选取如下的滑模切换面函数:
s=kqv+ω.
并证明闭环系统在所选滑模面上稳定运行,保证该系统s=0时的全局稳定性得到保证。
证明:选取Lyapunov函数为:
V0=2(1-q0)
Step2 设计控制律
设计以下滑模控制律
u=-K1s-D1F(s)+ueq (3)
其中,
在滑模控制律(5)-(6)作用下,从任意初始状态出发挠性卫星系统能在有限时间内到达滑模面s=0,即ω=-kqv,并在滑模面上保持稳定运行。
以下根据等效控制思想设计等效控制项ueq,选取以下下的Lyapunov函数:
令动力学方程(7)中外部干扰力矩为0,并使得到
其中,
(2)在实际应用中模态η和ψ难以测量时,针对挠性卫星系统(1)-(2)设计基于动态观测器的滑模控制律。
采用下式的动态观测器
定义如下模态观测误差
其中,正定对称矩阵P满足以下Lyapunov方程
Step1 选取滑模面
选取如下的滑模切换面函数:
s=kqv+ω
并证明闭环系统在所选滑模面上稳定运行,保证该系统s=0时的全局稳定性得到保证。
证明:选取Lyapunov函数为:
V0=2(1-q0)
Step2 设计控制律
设计以下滑模控制律
u=-K1s-D1F(s)+ueq (10)
其中,
在滑模控制律(12)-(13)作用下,从任意初始状态出发挠性卫星系统能在有限时间内到达滑模面s=0,即ω=-kqv,并在滑模面上保持稳定运行。
以下根据等效控制思想设计等效控制项ueq,选取以下下的Lyapunov函数:
令
其中,
最后,通过计算证明根据LaSalle理论可得在该控制律下可保证系统的全局稳定性。
以下进行仿真实验验证提出的两种控制律。分别验证当模态完全可测时提出的基于状态反馈的滑模控制律和当模态难以测量时提出的基于动态观测器的滑模控制律。
挠性卫星主体惯量矩阵Jmb为:
卫星本体与挠性部件的耦合作用矩阵δ:
四模态挠性卫星的自然角频率为
ωn1=0.7681rad/s,ωn2=1.1038rad/s,
ωn3=1.8733rad/s,ωn4=2.5496rad/s,
挠性附件的阻尼系数:
初始角速度如下
ω(0)=[0 0 0]T.
此外,挠性附件四个模态的初始值为
ηi=0.001,ψi=0.001,i=1,2,3,4.
外部干扰力矩d为
基于状态反馈的滑模控制器参数为:
k=6I3×3;D1=2I3×3;K1=9I3×3
基于动态观测器的滑模控制器的参数为:
k=0.2I3×3;D1=1000I3×3;K1=500I3×3;Q=10I8×8
本发明提供的一种带外部干扰力矩的挠性卫星姿态滑模控制方法,针对挠性卫星在运行过程中受到各种力矩干扰的情况,设计了一种基于滑模控制的姿态控制算法。该发明算法的目的在于解决挠性卫星在飞行过程中存在有界干扰时的姿态控制及挠性附件的振动控制问题,完成挠性卫星姿态机动的高鲁棒性控制。该发明利用姿态四元数建立挠性卫星运动学方程,采用混合坐标法对中心刚体带有挠性附件、外部干扰匹配的挠性航天器建立姿态动力学方程。采用等效控制思想,基于Lyapunov直接法设计一种基于状态反馈的姿态滑模控制律。并构造挠性模态观测器测量挠性状态变量,设计带挠性模态观测器的姿态滑模控制律。最后,进行仿真实验验证所设计的控制算法的具有良好的鲁棒性且具有较强的振动抑制能力。
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明,不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干简单推演或替换,都应当视为属于本发明的保护范围。
Claims (4)
1.一种带外部干扰力矩的挠性卫星姿态滑模控制方法,其特征在于:利用姿态四元数建立挠性卫星的运动学方程,采用混合坐标法对中心刚体带有挠性附件、外部干扰匹配的挠性卫星建立动力学方程;采用等效控制思想,基于Lyapunov直接法设计以下两种滑模控制律:针对挠性模态可测量的情况,设计一种基于状态反馈的挠性卫星姿态滑模控制律;针对挠性模态不可测量的情况,设计一种基于动态观测器的挠性卫星姿态滑模控制律。
2.根据权利要求1所述的带外部干扰力矩的挠性卫星姿态滑模控制方法,其特征在于:以姿态四元数为基础的挠性卫星的运动学方程和动力学方程如下:
<mrow>
<mover>
<mi>q</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mover>
<mi>q</mi>
<mo>&CenterDot;</mo>
</mover>
<mn>0</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mover>
<mi>q</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>v</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msup>
<mi>E</mi>
<mi>T</mi>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>q</mi>
<mn>0</mn>
</msub>
<mo>,</mo>
<msub>
<mi>q</mi>
<mi>v</mi>
</msub>
<mo>)</mo>
</mrow>
<mi>&omega;</mi>
</mrow>
其中,
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mover>
<mi>&omega;</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<msubsup>
<mi>J</mi>
<mrow>
<mi>m</mi>
<mi>b</mi>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>&lsqb;</mo>
<mo>-</mo>
<mi>&omega;</mi>
<mo>&times;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>J</mi>
<mrow>
<mi>m</mi>
<mi>b</mi>
</mrow>
</msub>
<mi>&omega;</mi>
<mo>+</mo>
<msup>
<mi>&delta;</mi>
<mi>T</mi>
</msup>
<mi>&psi;</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msup>
<mi>&delta;</mi>
<mi>T</mi>
</msup>
<mrow>
<mo>(</mo>
<mi>C</mi>
<mi>&psi;</mi>
<mo>+</mo>
<mi>K</mi>
<mi>&eta;</mi>
<mo>-</mo>
<mi>C</mi>
<mi>&delta;</mi>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>u</mi>
<mo>+</mo>
<mi>d</mi>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mover>
<mi>&eta;</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<mi>&psi;</mi>
<mo>-</mo>
<mi>&delta;</mi>
<mi>&omega;</mi>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mover>
<mi>&psi;</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mi>C</mi>
<mi>&psi;</mi>
<mo>+</mo>
<mi>K</mi>
<mi>&eta;</mi>
<mo>-</mo>
<mi>C</mi>
<mi>&delta;</mi>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,Jmb=J-δTδ,
其中,q0,qv分别为姿态四元数的标量部分与向量部分;ω表示挠性卫星的姿态角速度;δ为表示卫星本体与挠性部件的耦合作用矩阵;C,K分别表示阻尼矩阵和刚度矩阵,
<mrow>
<mi>K</mi>
<mo>=</mo>
<mi>d</mi>
<mi>i</mi>
<mi>a</mi>
<mi>g</mi>
<mo>{</mo>
<msubsup>
<mi>&omega;</mi>
<mrow>
<mi>n</mi>
<mn>1</mn>
</mrow>
<mn>2</mn>
</msubsup>
<mo>,</mo>
<msubsup>
<mi>&omega;</mi>
<mrow>
<mi>n</mi>
<mn>2</mn>
</mrow>
<mn>2</mn>
</msubsup>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msubsup>
<mi>&omega;</mi>
<mrow>
<mi>n</mi>
<mi>N</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>}</mo>
<mo>;</mo>
</mrow>
考虑N个弹性模态,其对应的自然角频率为ωni,i=1,2,…,N,对应的阻尼为i=1,2,…,N;Jmb为刚体部分的转动惯量,u表示控制力矩,d表示有界外部干扰力矩,卫星带挠性附件,并且含有有界外部干扰力矩。
3.根据权利要求1所述的带外部干扰力矩的挠性卫星姿态滑模控制方法,其特征在于:针对挠性模态可测的情形下,设计如下滑模面:
s=kqv+ω.
设计的基于状态反馈的挠性卫星姿态滑模控制律如下:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>u</mi>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>K</mi>
<mn>1</mn>
</msub>
<mi>s</mi>
<mo>-</mo>
<msub>
<mi>D</mi>
<mn>1</mn>
</msub>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mi>&omega;</mi>
<mo>&times;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>J</mi>
<mrow>
<mi>m</mi>
<mi>b</mi>
</mrow>
</msub>
<mi>&omega;</mi>
<mo>+</mo>
<msup>
<mi>&delta;</mi>
<mi>T</mi>
</msup>
<mi>&psi;</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msup>
<mi>&delta;</mi>
<mi>T</mi>
</msup>
<mrow>
<mo>(</mo>
<mi>C</mi>
<mi>&psi;</mi>
<mo>+</mo>
<mi>K</mi>
<mi>&eta;</mi>
<mo>-</mo>
<mi>C</mi>
<mi>&delta;</mi>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>J</mi>
<mrow>
<mi>m</mi>
<mi>b</mi>
</mrow>
</msub>
<mi>k</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msup>
<mi>E</mi>
<mi>T</mi>
</msup>
<mo>(</mo>
<msub>
<mi>q</mi>
<mi>v</mi>
</msub>
<mo>)</mo>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
<mo>.</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
F(s)={sgn(s1),sgn(s2),sgn(s3)}T
<mrow>
<mi>sgn</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mrow>
<msub>
<mi>s</mi>
<mi>i</mi>
</msub>
<mo>></mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<msub>
<mi>s</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>s</mi>
<mi>i</mi>
</msub>
<mo><</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>3</mn>
</mrow>
<mrow>
<mi>E</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>q</mi>
<mi>v</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>q</mi>
<mn>0</mn>
</msub>
<mi>I</mi>
<mo>-</mo>
<msubsup>
<mi>q</mi>
<mi>v</mi>
<mo>&times;</mo>
</msubsup>
</mrow>
min{D1(i)}>||d||
其中,K1,D1,k为正定对角矩阵。
4.根据权利要求1所述的带外部干扰力矩的挠性卫星姿态滑模控制方法,其特征在于:针对挠性模态不可测的情形下,设计如下滑模面:
s=kqv+ω
设计的基于动态观测器的挠性卫星姿态滑模控制律如下:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>u</mi>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>K</mi>
<mn>1</mn>
</msub>
<mi>s</mi>
<mo>-</mo>
<msub>
<mi>D</mi>
<mn>1</mn>
</msub>
<mi>F</mi>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mi>&omega;</mi>
<mo>&times;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>J</mi>
<mrow>
<mi>m</mi>
<mi>b</mi>
</mrow>
</msub>
<mi>&omega;</mi>
<mo>+</mo>
<msup>
<mi>&delta;</mi>
<mi>T</mi>
</msup>
<mover>
<mi>&psi;</mi>
<mo>^</mo>
</mover>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msup>
<mi>&delta;</mi>
<mi>T</mi>
</msup>
<mrow>
<mo>(</mo>
<mi>C</mi>
<mover>
<mi>&psi;</mi>
<mo>^</mo>
</mover>
<mo>+</mo>
<mi>K</mi>
<mover>
<mi>&eta;</mi>
<mo>^</mo>
</mover>
<mo>-</mo>
<mi>C</mi>
<mi>&delta;</mi>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>J</mi>
<mrow>
<mi>m</mi>
<mi>b</mi>
</mrow>
</msub>
<mi>k</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msup>
<mi>E</mi>
<mi>T</mi>
</msup>
<mo>(</mo>
<msub>
<mi>q</mi>
<mi>v</mi>
</msub>
<mo>)</mo>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mover>
<mover>
<mi>&eta;</mi>
<mo>^</mo>
</mover>
<mo>&CenterDot;</mo>
</mover>
</mtd>
</mtr>
<mtr>
<mtd>
<mover>
<mover>
<mi>&psi;</mi>
<mo>^</mo>
</mover>
<mo>&CenterDot;</mo>
</mover>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mi>I</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mi>K</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mi>C</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mover>
<mi>&eta;</mi>
<mo>^</mo>
</mover>
</mtd>
</mtr>
<mtr>
<mtd>
<mover>
<mi>&psi;</mi>
<mo>^</mo>
</mover>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>I</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>C</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>&delta;</mi>
<mi>&omega;</mi>
<mo>+</mo>
<msup>
<mi>P</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>K</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>C</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>&lsqb;</mo>
<mi>&delta;</mi>
<mi>&omega;</mi>
<mo>+</mo>
<msub>
<mi>k&delta;q</mi>
<mi>v</mi>
</msub>
<mo>&rsqb;</mo>
<mo>+</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>I</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<msup>
<mi>k&delta;&omega;</mi>
<mo>&times;</mo>
</msup>
<msub>
<mi>q</mi>
<mi>v</mi>
</msub>
</mrow>
<mo>&rsqb;</mo>
</mrow>
</mrow>
F(s)={sgn(s1),sgn(s2),sgn(s3)}T
<mrow>
<mi>s</mi>
<mi>g</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mrow>
<msub>
<mi>s</mi>
<mi>i</mi>
</msub>
<mo>></mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<msub>
<mi>s</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>s</mi>
<mi>i</mi>
</msub>
<mo><</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>3</mn>
</mrow>
<mrow>
<mi>E</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>q</mi>
<mi>v</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>q</mi>
<mn>0</mn>
</msub>
<mi>I</mi>
<mo>-</mo>
<msubsup>
<mi>q</mi>
<mi>v</mi>
<mo>&times;</mo>
</msubsup>
</mrow>
min{D1(i)}>||d||
<mrow>
<msup>
<mi>&omega;</mi>
<mo>&times;</mo>
</msup>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>&omega;</mi>
<mn>3</mn>
</msub>
</mrow>
</mtd>
<mtd>
<msub>
<mi>&omega;</mi>
<mn>2</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>&omega;</mi>
<mn>3</mn>
</msub>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>&omega;</mi>
<mn>1</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>&omega;</mi>
<mn>2</mn>
</msub>
</mrow>
</mtd>
<mtd>
<msub>
<mi>&omega;</mi>
<mn>1</mn>
</msub>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中,K1,D1,k为正定对角矩阵。正定对称矩阵P满足下述Lyapunov方程
<mrow>
<mi>P</mi>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mi>I</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mi>K</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mi>C</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>+</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mi>I</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mi>K</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mi>C</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
<mi>P</mi>
<mo>=</mo>
<mo>-</mo>
<mn>2</mn>
<mi>Q</mi>
<mo>.</mo>
</mrow>
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710822183 | 2017-09-13 | ||
CN2017108221839 | 2017-09-13 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107943062A true CN107943062A (zh) | 2018-04-20 |
Family
ID=61935088
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710901110.9A Pending CN107943062A (zh) | 2017-09-13 | 2017-09-28 | 带外部干扰力矩的挠性卫星姿态滑模控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107943062A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109213184A (zh) * | 2018-11-06 | 2019-01-15 | 哈尔滨工业大学(深圳) | 挠性航天器的有限时间多模态滑模姿态控制算法 |
CN111258325A (zh) * | 2020-01-08 | 2020-06-09 | 哈尔滨工业大学 | 对地遥感卫星高精度高性能的姿态容错控制方法、装置及计算机存储介质 |
CN111498147A (zh) * | 2020-04-03 | 2020-08-07 | 哈尔滨工业大学(深圳)(哈尔滨工业大学深圳科技创新研究院) | 挠性航天器的有限时间分段滑模姿态跟踪控制算法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6341249B1 (en) * | 1999-02-11 | 2002-01-22 | Guang Qian Xing | Autonomous unified on-board orbit and attitude control system for satellites |
CN105843043A (zh) * | 2016-05-17 | 2016-08-10 | 南京理工大学 | 一种电液负载模拟装置自适应鲁棒力控制方法 |
CN105843238A (zh) * | 2016-03-31 | 2016-08-10 | 北京航空航天大学 | 一种未知频率谐波干扰下的火星着陆器抗干扰控制器 |
CN107065913A (zh) * | 2017-05-10 | 2017-08-18 | 哈尔滨工业大学深圳研究生院 | 挠性航天器的滑模变结构姿态控制算法 |
-
2017
- 2017-09-28 CN CN201710901110.9A patent/CN107943062A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6341249B1 (en) * | 1999-02-11 | 2002-01-22 | Guang Qian Xing | Autonomous unified on-board orbit and attitude control system for satellites |
CN105843238A (zh) * | 2016-03-31 | 2016-08-10 | 北京航空航天大学 | 一种未知频率谐波干扰下的火星着陆器抗干扰控制器 |
CN105843043A (zh) * | 2016-05-17 | 2016-08-10 | 南京理工大学 | 一种电液负载模拟装置自适应鲁棒力控制方法 |
CN107065913A (zh) * | 2017-05-10 | 2017-08-18 | 哈尔滨工业大学深圳研究生院 | 挠性航天器的滑模变结构姿态控制算法 |
Non-Patent Citations (2)
Title |
---|
H.FOURATI 等: ""A quaternion-based Complementary Sliding Mode Observer for attitude estimation: application in free-ranging animal motions"", 《49TH IEEE CONFERENCE ON DECISION AND CONTROL》 * |
靳永强 等: ""含有参数不确定性的挠性航天器姿态跟踪滑模控制"", 《控制理论与应用》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109213184A (zh) * | 2018-11-06 | 2019-01-15 | 哈尔滨工业大学(深圳) | 挠性航天器的有限时间多模态滑模姿态控制算法 |
CN109213184B (zh) * | 2018-11-06 | 2021-06-08 | 哈尔滨工业大学(深圳) | 挠性航天器的有限时间多模态滑模姿态控制算法 |
CN111258325A (zh) * | 2020-01-08 | 2020-06-09 | 哈尔滨工业大学 | 对地遥感卫星高精度高性能的姿态容错控制方法、装置及计算机存储介质 |
CN111258325B (zh) * | 2020-01-08 | 2022-08-02 | 哈尔滨工业大学 | 对地遥感卫星高精度高性能的姿态容错控制方法、装置及计算机存储介质 |
CN111498147A (zh) * | 2020-04-03 | 2020-08-07 | 哈尔滨工业大学(深圳)(哈尔滨工业大学深圳科技创新研究院) | 挠性航天器的有限时间分段滑模姿态跟踪控制算法 |
CN111498147B (zh) * | 2020-04-03 | 2021-09-21 | 哈尔滨工业大学(深圳)(哈尔滨工业大学深圳科技创新研究院) | 挠性航天器的有限时间分段滑模姿态跟踪控制算法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Yeh | Sliding-mode adaptive attitude controller design for spacecrafts with thrusters | |
CN103412491B (zh) | 一种挠性航天器特征轴姿态机动指数时变滑模控制方法 | |
CN105911866B (zh) | 四旋翼无人飞行器的有限时间全阶滑模控制方法 | |
CN111498147B (zh) | 挠性航天器的有限时间分段滑模姿态跟踪控制算法 | |
Meng et al. | Space robots with flexible appendages: dynamic modeling, coupling measurement, and vibration suppression | |
CN108241292B (zh) | 一种基于扩张状态观测器的水下机器人滑模控制方法 | |
CN107831655B (zh) | 微陀螺仪的分数阶自适应反演模糊滑模控制方法 | |
CN107943062A (zh) | 带外部干扰力矩的挠性卫星姿态滑模控制方法 | |
CN109213184B (zh) | 挠性航天器的有限时间多模态滑模姿态控制算法 | |
CN107807657B (zh) | 一种基于路径规划的挠性航天器姿态自适应控制方法 | |
Schlanbusch et al. | Hybrid attitude tracking of rigid bodies without angular velocity measurement | |
CN110543184B (zh) | 一种刚性飞行器的固定时间神经网络控制方法 | |
CN111874266B (zh) | 一种刚体航天器的抗退绕滑模姿态机动控制方法及系统 | |
CN111258216A (zh) | 一种适用于四旋翼飞行器的滑模重复控制器 | |
CN105223961A (zh) | 一种用于控制力矩陀螺奇异规避的航天器姿态控制方法 | |
CN108427272A (zh) | 基于挠性模态观测的挠性航天器姿态控制和振动抑制方法 | |
Zarafshan et al. | Adaptive hybrid suppression control of space free-flying robots with flexible appendages | |
Navabi et al. | Fuzzy control of fuel sloshing in a spacecraft | |
Meng et al. | Vibration suppression of a large flexible spacecraft for on-orbit operation | |
CN110333656A (zh) | 一种基于互联系统方法的挠性航天器容错控制方法 | |
Huang et al. | Robust H∞ control for spacecraft formation flying with coupled translational and rotation dynamics | |
Qi et al. | Constrained attitude tracking control and active sloshing suppression for liquid-filled spacecraft | |
CN111216928A (zh) | 一种携带太阳帆板的失效卫星绳系拖曳控制方法及系统 | |
CN109799704A (zh) | 一种挠性航天器姿态与结构振动的耦合控制方法 | |
CN105717791B (zh) | 一种自适应h无穷控制的悬臂梁振动控制方法 |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180420 |
|
RJ01 | Rejection of invention patent application after publication |