CN103455035B - 基于反步设计和非线性反馈的pd+姿态控制律设计方法 - Google Patents
基于反步设计和非线性反馈的pd+姿态控制律设计方法 Download PDFInfo
- Publication number
- CN103455035B CN103455035B CN201310376138.7A CN201310376138A CN103455035B CN 103455035 B CN103455035 B CN 103455035B CN 201310376138 A CN201310376138 A CN 201310376138A CN 103455035 B CN103455035 B CN 103455035B
- Authority
- CN
- China
- Prior art keywords
- omega
- sigma
- centerdot
- attitude
- relative attitude
- 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
Landscapes
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明涉及一种基于反步设计和非线性反馈的PD+姿态控制律设计方法,属于航天器高性能姿态控制技术领域。本方法建立级联形式的相对姿态运动方程,根据反步设计思想将PD+姿态控制律的设计问题分解为两个相对姿态运动子***的稳定控制律设计问题。进一步根据相对姿态运动方程的级联关系设计PD+姿态控制律,以期望的动态响应稳定整个姿态控制***。最后通过对平衡点处的闭环***阻尼比进行约束来避免超调现象。本发明提高了现有基于反馈线性化方法设计PD+姿态控制律的灵活性,能够在不增加控制力矩幅值的前提下能够很大程度上地提高闭环***的响应速度。
Description
技术领域
本发明涉及一种基于反步设计和非线性反馈的PD+(ProportionalDerivativePlus)姿态控制律设计方法,属于航天器高性能姿态控制技术领域。
背景技术
为了完成交会对接、编队飞行、在轨服务等新兴空间应用任务,航天器需要频繁地进行大角度姿态机动/跟踪。在这种情况下,姿态控制***呈现出强非线性和强耦合的特点。此时,为了获得较高的控制性能,就需要用到非线性控制技术。其中,基于反馈线性化技术设计的PD+姿态控制律因其设计简单便于实现等特点被广泛采用。这种姿态控制律设计方法首先通过反馈线性化技术将相对姿态动力学方程线性化,然后针对线性化后的新***设计相对姿态变量(相对姿态参数和相对姿态角速度)的反馈项。WieB.[WieB.,WeissH.,ArapostathisA.Quaternionfeedbackregulatorforspacecrafteigenaxisrotations[J].JournalofGuidance,Control,andDynamics,12(3),1989:375-380.]针对刚性航天器的姿态稳定控制任务设计了PD+姿态调节律,并给出了控制参数的整定方法;WenJ.T.[WenJ.T.,KennethK.D.Theattitudecontrolproblem[J].IEEETransactionsonAutomaticControl,36(10),1991:1148-1162.]针对刚体姿态跟踪问题设计了PD+姿态控制律,并利用自适应控制技术提高了PD+姿态控制律的鲁棒性。随着滑模控制技术的发展,利用积分滑模控制技术可以保证PD+控制律的鲁棒性[UtkinV.,ShiJ.X.Integralslidingmodeinsystemsoperatingunderuncertaintyconditions[C].Proceedingsofthe35thIEEEConferenceonDecisionandControl,Japan,1996:4591-4596.]。在这种情况下,PD+姿态控制律作为积分滑模姿态控制律中的标称控制律并决定着闭环控制***的动态性能。
对于根据反馈线性化方法设计的PD+姿态控制律来说,当确定了需要抵消的非线性项后,设计者仅能通过调节控制增益来满足控制需求。而现有PD+姿态控制律中的PD反馈项是相对姿态变量的线性比例放大,若选择大控制增益来加快***响应,那么当相对姿态变量较大时,姿态控制律会产生很大的控制力矩幅值从而容易导致执行机构的饱和问题;反之,若选择小控制增益来降低执行机构发生饱和的几率,那么同时也会降低***的响应速度。为了解决这一问题,SchlanbuschR.[SchlanbuschR.,LoriaA.,KristiansenR.,NicklassonP.J.PD+attitudecontrolofrigidbodieswithimprovedperformance[C].Proceedingsofthe49thIEEEConferenceonDecisionandControl,Atlanta,2010:7069-7074.]利用非线性反馈技术使得PD反馈项具有指数增加的控制增益,增加了PD+姿态控制律的设计自由度。然而,这种改进的PD+姿态控制律在相对姿态变量较大时具有大控制增益,在加快***响应速度的同时更容易导致执行器饱和问题。另外,当相对姿态变量处于平衡点附近时,这种改进的PD+姿态控制律通过减少控制增益来避免超调,但同时也降低了***的响应速度。
发明内容
本发明的目的是为提高现有PD+姿态控制律的控制性能,提出一种PD+姿态控制律设计方法,结合反步设计方法和非线性反馈技术,在不增加控制力矩幅值的前提下提高***的响应速度,可用于刚性航天器进行姿态机动/跟踪时的姿态控制。
本发明的技术方案为:建立级联形式的相对姿态运动方程,根据反步设计思想将PD+姿态控制律的设计问题分解为两个相对姿态运动子***的稳定控制律设计问题。对于相对姿态运动学子***,利用非线性反馈技术设计虚拟控制律使闭环子***具有期望的动态响应;再针对相对姿态动力学子***设计辅助姿态控制律,使实际相对姿态角速度能够快速地跟踪上虚拟控制律的输出值;进一步根据相对姿态运动方程的级联关系设计PD+姿态控制律,以期望的动态响应稳定整个姿态控制***。最后通过对平衡点处的闭环***阻尼比进行约束来避免超调现象。
具体包括以下步骤:
步骤1,以进行姿态机动/跟踪的刚性航天器为对象,在姿态运动的构型空间内定义相对姿态变量,在航天器本体坐标系下建立级联形式的相对姿态动力学方程和相对姿态运动学方程。具体方法为:
在航天器本体坐标系下定义相对姿态变量如下:
ωe=ωb-R(σe)ωd(2)
式中,σb为航天器本体坐标系姿态对应的MRPs(ModifiedRodriguesParameters,修正罗德里格斯参数)矢量在本体坐标系下的向量表示,σd为惯性空间系姿态对应的MRPs矢量在惯性空间系下的向量表示,σe表示航天器本体坐标系与惯性空间系之间的相对姿态对应的MRPs矢量在本体坐标系下的向量表示,ωb表示航天器角速度矢量在本体坐标系下的向量表示,ωd表示参考角速度矢量在惯性空间系下的向量表示,ωe表示航天器本体坐标系与惯性空间系之间的相对姿态角速度矢量在本体坐标系下的向量表示;||·||表示向量的Euclidean范数,(·)×表示向量的反对称矩阵算子,(·)T表示向量或矩阵的转置算子,表示MRPs的乘法算子;航天器本体坐标系与惯性空间系之间的转移矩阵R为:
式中,I3表示3×3的单位矩阵。
建立相对姿态运动学方程为:
相对姿态动力学方程为:
式中,雅可比矩阵M为:
J为航天器惯量阵张量在本体坐标系下的矩阵表示,Tc为控制力矩矢量在本体坐标系下的向量表示。
所述惯性空间系以地心为坐标原点,xi轴在地球赤道平面内指向春分点,zi轴指向北极,与地球自旋轴重合,yi与xi和zi成右手正交系。
所述航天器本体坐标系,以航天器质心为坐标原点,本体系的x轴与航天器结构的主轴(对称轴)相重合,符合右手法则。航天器的惯量阵张量在本体系下的矩阵表示为对角阵。
步骤2,针对步骤1建立的相对姿态运动学方程,将相对姿态角速度看作虚拟控制输入,利用非线性反馈技术设计虚拟控制律保证相对姿态运动学子***的稳定性。基于非线性反馈技术设计的虚拟控制律在相对姿态变量较大时的输出为有界值,能够减少大误差情况下的控制力矩幅值;而在相对姿态变量较小时的输出大于线性反馈虚拟控制律,能够加快***的响应速度。具体方法为:
将相对姿态角速度看作相对姿态运动学子***的虚拟控制输入,并将相对姿态运动学方程改写为:
式中,是虚拟控制输入,表示相对姿态角速度ωe的期望值。
为达到限制虚拟控制输入在σe较大时的输出值,而在σe较小时使虚拟控制输入的输出大于线性反馈的目的,设计两种基于非线性反馈技术的虚拟控制律和(只使用一种即可实现PD+控制):
式中,ka>0,kt>0,向量函数arctan(aσe)和tanh(aσe)分别为:
arctan(aσe)=[arctan(aσe1),arctan(aσe2),arctan(aσe3)]T(10)
tanh(aσe)=[tanh(aσe1),tanh(aσe2),tanh(aσe3)]T(11)
arctan和tanh分别为反正切和双曲正切函数,参数a>0并满足
arctan(a|σei|)>|σei|,(i=1,2,3)(12)
tanh(a|σei|)>|σei|,(i=1,2,3)(13)
步骤3,引入辅助变量表示实际相对姿态角速度与步骤2设计的虚拟控制输入之间的误差,并以该辅助变量为状态变量来描述步骤1建立的相对姿态动力学方程,抵消相对姿态动力学方程中的非线性项并设计辅助变量的比例反馈项来建立辅助姿态控制律,保证实际相对姿态角速度对虚拟控制输入的跟踪。
引入辅助变量表示实际相对姿态角速度ωe与虚拟控制输入之间的误差,进一步以z为状态变量,将相对姿态动力学方程(5)改写为:
式中,表示辅助姿态控制律。对式(8)和(9)求导,得到虚拟控制率和的新形式分别为:
为稳定相对姿态动力学子***(18),将辅助姿态控制律设计为:
式中,kd>0。
步骤4,考虑相对姿态动力学方程和相对姿态运动学方程的级联关系,综合步骤2设计的虚拟控制律和步骤3设计的辅助姿态控制律,设计PD+姿态控制律以保证整个姿态控制***的稳定性。具体方法为:
将相对姿态运动学方程(4)和相对姿态动力学方程(5)改写为:
本发明提出的基于反步设计和非线性反馈技术的PD+姿态控制律包括两个部分。第一部分是相对姿态动力学方程(25)右端的抵消项,第二部分是状态变量σe和z的反馈项。具体形式为:
步骤5,将步骤4设计的PD+姿态控制律代入相对姿态动力学方程中,在平衡点附近利用小角度假设将闭环***方程近似为以欧拉主轴角为状态变量的二阶阻尼谐振***方程,限制平衡点处的闭环***阻尼比以避免超调现象并确定控制参数所需满足的条件。具体方法为:
姿态控制律(26)代入相对姿态动力学方程(25),整理有:
在相对姿态变量处于平衡点附近时使用小角度假设,即令:
(30)
M≈I3
arctan(aσe)≈aσe,tanh(aσe)≈aσe
式中,θe表示σe对应的欧拉主轴角,n表示σe对应的特征轴矢量在航天器本体系下的向量表示。
将式(30)代入式(29)中,对于虚拟控制律(8)有:
对于虚拟控制律(9)有:
根据二阶阻尼谐振***的阻尼比和自然振荡频率的概念,为了使闭环***在σe=0时为临界阻尼***,则参数ka,kt,kd,μ以及a需要满足:
有益效果
本发明方法能够有效地提高姿态控制***的控制性能,与现有技术相比的优点在于:
1)反步设计在很大程度上提高了现有基于反馈线性化方法设计PD+姿态控制律的灵活性。在反步设计框架内,设计者不仅可以通过调节控制增益来改变闭环***的动态响应,而且在两个子***的稳定控制律设计上可以采用多种线性或非线性控制方法来实现期望的控制性能;
2)利用非线性反馈技术能够有效地解决现有PD+姿态控制律在加快***响应速度和减少控制力矩幅值两种性能指标上的权衡问题,本发明提出的PD+姿态控制律通过对相对姿态变量进行非线性反馈,在不增加控制力矩幅值的前提下能够很大程度上地提高闭环***的响应速度。
附图说明
图1为本发明基于反步设计和非线性反馈技术的PD+姿态控制律设计流程图;
图2为具体实施中基于反步设计和非线性反馈技术的PD+姿态控制律设计框图;
图3为步骤2中虚拟控制律(8)和(9)使用的两种非线性函数在a=8时的输出曲线;
图4为具体实施中航天器姿态控制***跟踪给定姿态指令时,采用现有PD+控制以及本发明提出的PD+控制的σe响应对比图。其中,(a)为现有PD+控制律作用下的σe响应曲线;(b)为本发明在虚拟控制律(8)基础上设计的PD+控制律作用下的σe响应曲线;(c)为本发明在虚拟控制律(9)基础上设计的PD+控制律作用下的σe响应曲线;
图5为具体实施中航天器姿态控制***跟踪给定姿态指令时,采用现有PD+控制以及本发明提出的PD+控制的控制力矩对比图。其中,(a)为现有PD+控制律对应的控制力矩曲线;(b)为本发明在虚拟控制律(8)基础上设计的PD+控制律对应的控制力矩曲线;(c)为本发明在虚拟控制律(9)基础上设计的PD+控制律对应的控制力矩曲线。
具体实施方式
为了更好地说明本发明的目的和优点,下面结合附图和实施例加以进一步阐述。
本发明的设计流程如图1所示。依据本发明方法实施的PD+姿态控制律设计如图2所示,该控制律能够提高航天器姿态控制***对期望姿态变量的跟踪性能。
1)建立刚性航天器进行姿态机动/跟踪时的相对姿态运动模型。首先,在航天器本体坐标系下定义相对姿态变量如下:
ωe=ωb-R(σe)ωd(35)
式中,σb为航天器本体坐标系姿态对应的MRPs(ModifiedRodriguesParameters,修正罗德里格斯参数)矢量在本体坐标系下的向量表示,σd为惯性空间系姿态对应的MRPs矢量在惯性空间系下的向量表示,σe表示航天器本体坐标系与惯性空间系之间的相对姿态对应的MRPs矢量在本体坐标系下的向量表示,ωb表示航天器角速度矢量在本体坐标系下的向量表示,ωd表示参考角速度矢量在惯性空间系下的向量表示,ωe表示航天器本体坐标系与惯性空间系之间的相对姿态角速度矢量在本体坐标系下的向量表示;||·||表示向量的Euclidean范数,(·)×表示向量的反对称矩阵算子,(·)T表示向量或矩阵的转置算子,表示MRPs的乘法算子;航天器本体坐标系与惯性空间之间的转移矩阵R为:
式中,I3表示3×3的单位矩阵。
在相对姿态变量定义的基础上,相对姿态运动学方程为:
相对姿态动力学方程为:
式中,雅可比矩阵M为:
J为航天器惯量阵张量在本体坐标系下的矩阵表示,Tc为控制力矩矢量在本体坐标系下的向量表示。
2)相对姿态动力学方程(38)和相对姿态运动学方程(37)通过级联关系共同描述了航天器的相对姿态运动规律。根据反步设计思想,将相对姿态角速度看作相对姿态运动学子***的虚拟控制输入,并将相对姿态运动学方程(37)改写为:
式中,是虚拟控制输入,表示相对姿态角速度ωe的期望值。
为了在不增加控制力矩幅值的前提下增加***的响应速度,本发明采用非线性反馈技术设计虚拟控制律。具体实施是限制虚拟控制输入在σe较大时的输出值,而在σe较小时使虚拟控制输入的输出大于线性反馈的情况。为此,设计如下两种基于非线性反馈技术的虚拟控制律:
式中,ka>0,kt>0,向量函数arctan(aσe)和tanh(aσe)分别为:
arctan(aσe)=[arctan(aσe1),arctan(aσe2),arctan(aσe3)]T(43)
tanh(aσe)=[tanh(aσe1),tanh(aσe2),tanh(aσe3)]T(44)
arctan和tanh分别为反正切和双曲正切函数,参数a>0并满足
arctan(a|σei|)>|σei|,(i=1,2,3)(45)
tanh(a|σei|)>|σei|,(i=1,2,3)(46)
下面对所设计的非线性反馈虚拟控制律(41)和(42)进行稳定性分析。
考虑如下正定的Lyapunov函数:
式中,μ>0。
对式(47)求关于时间的导数:
将虚拟控制律(41)代入上式,可进一步表示为:
由于arctan(aσei)与σei同号,上式负定。
类似地,虚拟控制律(42)代入式(48),同样由于tanh(aσei)与σei同号有:
根据式(49)和(50)的推导结果可知,采用虚拟控制律(41)和(42)时,相对姿态运动学子***是全局一致渐近稳定的。
虚拟控制律(41)和(42)采用的非线性函数具有有界的输出值。其中,反正切函数arctan的值域为(-π/2,π/2),双曲正切函数tanh的值域为(-1,1)。因此,虚拟控制律(41)和(42)可以通过参数ka和kt来调节虚拟控制输入的幅值,从而减少相对姿态动力学子***需要跟踪的期望值并降低控制力矩幅值。另一方面,参数a能够改变非线性函数的变化规律。特别地,图3描述了arctan(aσe1)和tanh(aσe1)在a=8时的输出曲线。可以看出,当σe1<1.5时,arctan(aσe1)>σe1,当σe1<1时,tanh(aσe1)>σe1,即非线性虚拟控制律(41)和(42)在σe较小时的输出大于线性反馈的情况,从而能够加快***的响应速度。
由此证明所设计的虚拟控制律(41)和(42)能在不增加控制力矩幅值的前提下加快***的响应速度,从而提高现有PD+姿态控制律的控制性能。
3)在步骤2的基础上,引入辅助变量表示实际相对姿态角速度ωe与虚拟控制输入之间的误差,进一步以z为状态变量,将相对姿态动力学方程(38)改写为:
式中,表示辅助姿态控制律。对于虚拟控制律(41)和(42)来说,和分别为:
为稳定相对姿态动力学子***(51),将辅助姿态控制律设计为:
式中,kd>0。
下面对所设计的辅助姿态控制律(54)进行稳定性分析。
考虑如下正定的Lyapunov函数:
对式(55)求关于时间的导数并代入辅助姿态控制律(54):
根据式(56)的推导结果可知,采用辅助姿态控制律(54)时,相对姿态动力学子***是全局一致渐近稳定的。
由此证明,辅助姿态控制律(54)能够保证实际相对姿态角速度跟踪上虚拟控制输入,从而使姿态控制***具有虚拟控制输入确定的期望动态响应。
4)在步骤2和步骤3的基础上,将相对姿态运动学方程(37)和相对姿态动力学方程(38)改写为:
与现有基于反馈线性化方法设计的PD+姿态控制律一样,本发明提出的基于反步设计和非线性反馈技术的PD+姿态控制律也包括两个部分。第一部分是相对姿态动力学方程(58)右端的抵消项,第二部分是状态变量σe和z的反馈项。具体形式如下所示:
下面对所设计的PD+姿态控制律(59)进行稳定性分析:
考虑如下正定的Lyapunov函数:
对式(60)求关于时间的导数:
由于步骤2设计的虚拟控制律(41)和(42)都能够保证因此Lyapunov函数(60)的导数负定。根据式(61)的推导结果可知,采用PD+姿态控制律(59)时,姿态控制***是全局一致渐近稳定的。
5)将步骤4设计的姿态控制律(59)代入相对姿态动力学方程(58),整理有:
在相对姿态变量处于平衡点附近时使用小角度假设,即令:
(63)
M≈I3
arctan(aσe)≈aσe,tanh(aσe)≈aσe
式中,θe表示σe对应的欧拉主轴角,n表示σe对应的特征轴矢量在航天器本体系下的向量表示。
将式(63)代入式(62)中,对于虚拟控制律(41)有:
对于虚拟控制律(42)有:
根据二阶阻尼谐振***的阻尼比和自然振荡频率的概念,为了使闭环***在σe=0时为临界阻尼***,则参数ka(kt),kd,μ以及a需要满足:
实施例
本发明在Matlab2009a环境下进行仿真验证。航天器的惯量阵为
初始惯性姿态变量为:σb(t0)=[0,0,0]T,ωb(t0)=[0,0,0]T(rad/s)。期望姿态变量为目标轨道LVLH系对应的σd和ωd。其中,目标轨道为圆形轨道,轨道半径a=6899807(m),偏心率e=0,轨道倾角i=30(deg),升交点赤经Ω=60(deg),近地点幅角ω=0(deg),初始真近点角f(t0)=90(deg)。
控制器参数选择:μ=0.02,kd=0.05,ka=0.08,kt=0.12,a=12。
为了体现本发明提出的控制方法的优越性,与现有基于反馈线性化设计的PD+姿态控制律的姿态控制效果进行比较。
本发明基于反步设计和非线性反馈技术设计的PD+姿态控制律和现有基于反馈线性化设计的PD+姿态控制律作用下的σe响应曲线如图4所示。若以σe收敛至其初始值5%所需的时间作为调节时间,现有PD+姿态控制作用下的调节时间约为60秒,在本发明基于虚拟控制律(8)和(9)设计的PD+控制律作用下的调节时间分别约为33秒和32秒。可以看出,与现有PD+姿态控制相比,采用本发明提出的基于反步设计和非线性反馈技术设计的PD+姿态控制能够有效地减少姿态控制的调节时间,并且利用步骤5的参数整定方法不会引起超调现象。
图5分别给出了采用基于反步设计和非线性反馈技术设计的PD+姿态控制律和基于反馈线性化设计的PD+姿态控制律对应的控制力矩曲线。从图中可以看出本发明提出的PD+姿态控制律在减少调节时间的同时没有增加控制力矩的幅值,从而有效地解决了现有PD+姿态控制在加快***响应和减少控制力矩幅值两个性能指标上存在的权衡问题。
Claims (3)
1.基于反步设计和非线性反馈的PD+姿态控制律设计方法,其特征在于:具体包括以下步骤:
步骤1,以进行姿态机动/跟踪的刚性航天器为对象,在姿态运动的构型空间内定义相对姿态变量,在航天器本体坐标系下建立级联形式的相对姿态动力学方程和相对姿态运动学方程;具体方法为:
在航天器本体坐标系下定义相对姿态变量如下:
ωe=ωb-R(σe)ωd(2)
式中,σb为航天器本体坐标系姿态对应的MRPs矢量在本体坐标系下的向量表示,σd为惯性空间系姿态对应的MRPs矢量在惯性空间系下的向量表示,σe表示航天器本体坐标系与惯性空间系之间的相对姿态对应的MRPs矢量在本体坐标系下的向量表示,ωb表示航天器角速度矢量在本体坐标系下的向量表示,ωd表示参考角速度矢量在惯性空间系下的向量表示,ωe表示航天器本体坐标系与惯性空间系之间的相对姿态角速度矢量在本体坐标系下的向量表示;||·||表示向量的Euclidean范数,(·)×表示向量的反对称矩阵算子,(·)T表示向量或矩阵的转置算子,表示MRPs的乘法算子;航天器本体坐标系与惯性空间系之间的转移矩阵R为:
式中,I3表示3×3的单位矩阵;
建立相对姿态运动学方程为:
相对姿态动力学方程为:
式中,雅可比矩阵M为:
J为航天器惯量阵矢量在本体坐标系下的矩阵表示,Tc为控制力矩矢量在本体坐标系下的向量表示;
步骤2,针对步骤1建立的相对姿态运动学方程,将相对姿态角速度看作虚拟控制输入,利用非线性反馈技术设计虚拟控制律;具体方法为:
将相对姿态角速度看作相对姿态运动学子***的虚拟控制输入,将相对姿态运动学方程改写为:
式中,是虚拟控制输入,表示相对姿态角速度ωe的期望值;
设计基于非线性反馈技术的虚拟控制律
式中,ka>0,向量函数arctan(aσe)为:
arctan(aσe)=[arctan(aσe1),arctan(aσe2),arctan(aσe3)]T(10)arctan为反正切函数,参数a>0并满足
arctan(a|σei|)>|σei|,(i=1,2,3)(12)
步骤3,引入辅助变量表示实际相对姿态角速度ωe与虚拟控制输入之间的误差,以z为状态变量,相对姿态动力学方程改写为:
式中,表示辅助姿态控制律;对虚拟控制律求导:
设计辅助姿态控制律为:
式中,kd>0;
步骤4,综合步骤2设计的虚拟控制律和步骤3设计的辅助姿态控制律,设计PD+姿态控制律;具体方法为:
将相对姿态运动学方程和相对姿态动力学方程改写为:
所述PD+姿态控制律的具体形式为:
步骤5,将步骤4设计的PD+姿态控制律代入相对姿态动力学方程中,在平衡点附近利用小角度假设将闭环***方程近似为以欧拉主轴角为状态变量的二阶阻尼谐振***方程,限制平衡点处的闭环***阻尼比;具体方法为:
PD+姿态控制律代入改写后的相对姿态动力学方程,得到:
在相对姿态变量处于平衡点附近时使用小角度假设,令:
M≈I3(30)
arctan(aσe)≈aσe,tanh(aσe)≈aσe
式中,θe表示σe对应的欧拉主轴角,n表示σe对应的特征轴矢量在航天器本体系下的向量表示;
对于虚拟控制律有:
参数ka,kd,μ以及a满足:
2.根据权利要求1所述的基于反步设计和非线性反馈的PD+姿态控制律设计方法,其特征在于:所述惯性空间系以地心为坐标原点,xi轴在地球赤道平面内指向春分点,zi轴指向北极,与地球自旋轴重合,yi与xi和zi成右手正交系;所述航天器本体坐标系,以航天器质心为坐标原点,本体系的x轴与航天器结构的主轴,即对称轴相重合,符合右手法则;航天器的惯量阵矢量在本体系下的矩阵表示为对角阵。
3.根据权利要求1所述的基于反步设计和非线性反馈的PD+姿态控制律设计方法,其特征在于:步骤2所述的基于非线性反馈技术的虚拟控制律能够替换为另一种形式的虚拟控制律:
其中,kt>0,向量函数tanh(aσe)为:
tanh(aσe)=[tanh(aσe1),tanh(aσe2),tanh(aσe3)]T(11)tanh为双曲正切函数,参数a>0并满足
tanh(a|σei|)>|σei|,(i=1,2,3)(13)
求导得到:
对于虚拟控制律有:
参数kt,kd,μ以及a满足:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310376138.7A CN103455035B (zh) | 2013-08-26 | 2013-08-26 | 基于反步设计和非线性反馈的pd+姿态控制律设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310376138.7A CN103455035B (zh) | 2013-08-26 | 2013-08-26 | 基于反步设计和非线性反馈的pd+姿态控制律设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103455035A CN103455035A (zh) | 2013-12-18 |
CN103455035B true CN103455035B (zh) | 2015-12-23 |
Family
ID=49737524
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310376138.7A Active CN103455035B (zh) | 2013-08-26 | 2013-08-26 | 基于反步设计和非线性反馈的pd+姿态控制律设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103455035B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104635509B (zh) * | 2014-12-03 | 2016-09-21 | 陕西科技大学 | 多轴联动***精度控制的反步-滑模控制器及控制方法 |
CN104503241B (zh) * | 2014-12-23 | 2017-03-01 | 哈尔滨工业大学 | 卫星姿态控制***的转动惯量确定方法 |
CN108181916B (zh) * | 2017-12-29 | 2020-04-24 | 清华大学 | 小卫星相对姿态的控制方法及装置 |
CN108762070B (zh) * | 2018-05-10 | 2021-04-20 | 南京邮电大学 | 一种欠驱动无人机的分数阶控制方法 |
CN110647104B (zh) * | 2019-09-02 | 2021-05-11 | 华南理工大学 | 一种基于边界扰动观测器的柔性立管反步边界控制方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0493227A1 (fr) * | 1990-12-21 | 1992-07-01 | AEROSPATIALE Société Nationale Industrielle | Système de contrôle d'attitude pour satellite stabilisé 3-axes sur orbite à faible inclinaison |
EP1064591B1 (en) * | 1998-03-16 | 2003-09-10 | Honeywell Inc. | Continuous spacecraft attitude control that avoids cmg array singularities |
EP1840691A3 (en) * | 2006-03-29 | 2008-05-28 | Honeywell International Inc. | Satellite active motion damping |
CN102346486A (zh) * | 2011-05-27 | 2012-02-08 | 哈尔滨工业大学 | 一种基于pd控制律的姿态控制ip核及采用该ip核的卫星姿态控制力矩计算*** |
CN103076807A (zh) * | 2012-12-27 | 2013-05-01 | 北京航空航天大学 | 一种欠驱动柔性航天器姿态稳定的控制方法 |
-
2013
- 2013-08-26 CN CN201310376138.7A patent/CN103455035B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0493227A1 (fr) * | 1990-12-21 | 1992-07-01 | AEROSPATIALE Société Nationale Industrielle | Système de contrôle d'attitude pour satellite stabilisé 3-axes sur orbite à faible inclinaison |
EP1064591B1 (en) * | 1998-03-16 | 2003-09-10 | Honeywell Inc. | Continuous spacecraft attitude control that avoids cmg array singularities |
EP1840691A3 (en) * | 2006-03-29 | 2008-05-28 | Honeywell International Inc. | Satellite active motion damping |
CN102346486A (zh) * | 2011-05-27 | 2012-02-08 | 哈尔滨工业大学 | 一种基于pd控制律的姿态控制ip核及采用该ip核的卫星姿态控制力矩计算*** |
CN103076807A (zh) * | 2012-12-27 | 2013-05-01 | 北京航空航天大学 | 一种欠驱动柔性航天器姿态稳定的控制方法 |
Non-Patent Citations (2)
Title |
---|
Robust Backstepping Control for Slew Maneuver Using Nonlinear Tracking Function;Ki-Seok Kim,Youdan Kim;《IEEE TRANSACTIONS ON CONTROL SYSTEMS TECHNOLOGY》;20031130;第11卷(第6期);全文 * |
刚体航天器姿态跟踪***的自适应积分滑模控制;丛炳龙等;《航空学报》;20130325;第34卷(第3期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN103455035A (zh) | 2013-12-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Hu | Robust adaptive sliding mode attitude maneuvering and vibration damping of three-axis-stabilized flexible spacecraft with actuator saturation limits | |
CN103455035B (zh) | 基于反步设计和非线性反馈的pd+姿态控制律设计方法 | |
CN105242676A (zh) | 一种有限时间收敛时变滑模姿态控制方法 | |
CN102929283A (zh) | 基于sdre的再入飞行器自适应最优滑模姿态控制方法 | |
CN107807657B (zh) | 一种基于路径规划的挠性航天器姿态自适应控制方法 | |
Hu | Robust adaptive backstepping attitude and vibration control with L2-gain performance for flexible spacecraft under angular velocity constraint | |
CN104950671A (zh) | 基于自适应模糊的再入飞行器pid型滑模姿态控制方法 | |
Hu et al. | Robust adaptive backstepping attitude stabilization and vibration reduction of flexible spacecraft subject to actuator saturation | |
Hu et al. | Adaptive fuzzy integral sliding mode control for flexible air-breathing hypersonic vehicles subject to input nonlinearity | |
CN105629734A (zh) | 一种近空间飞行器的轨迹跟踪控制方法 | |
CN115639830B (zh) | 一种空地智能体协同编队控制***及其编队控制方法 | |
Zhang et al. | Nonlinear disturbance observer based adaptive super-twisting sliding mode control for generic hypersonic vehicles with coupled multisource disturbances | |
Xiao et al. | Adaptive sliding mode backstepping control for attitude tracking of flexible spacecraft under input saturation and singularity | |
CN115686048A (zh) | 执行器受限航天器交会***的动态触发有限时间控制方法 | |
Hu et al. | Control of three-axis stabilized flexible spacecraft using variable structure strategies subject to input nonlinearities | |
Tang et al. | Height and attitude active disturbance rejection controller design of a small-scale helicopter | |
CN113485396B (zh) | 一种航天器交会对接最后逼近段相对轨道与姿态跟踪控制方法 | |
CN104460678B (zh) | 一种基于Sigmoid型非线性滑模函数的航天器姿态控制方法 | |
Dong et al. | Adaptive nonsingular fixed‐time control for hypersonic flight vehicle considering angle of attack constraints | |
Ding et al. | Fault‐tolerant control of an aerial manipulator with guaranteed tracking performance | |
Hu | Robust Adaptive Attitude Tracking Control With L 2-Gain Performance and Vibration Reduction of an Orbiting Flexible Spacecraft | |
Castillo et al. | Small scale helicopter analysis and controller design for non-aggressive flights | |
Li et al. | Novel fuzzy approximation control scheme for flexible air-breathing hypersonic vehicles with non-affine dynamics and amplitude and rate constraints | |
Azza El-S et al. | Satellite attitude maneuver using sliding mode control under body angular velocity constraints | |
Gregory | Dynamic inversion to control large flexible transport aircraft |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |