CN105353615B - 一种基于滑模观测器的四旋翼飞行器的主动容错控制方法 - Google Patents
一种基于滑模观测器的四旋翼飞行器的主动容错控制方法 Download PDFInfo
- Publication number
- CN105353615B CN105353615B CN201510776314.5A CN201510776314A CN105353615B CN 105353615 B CN105353615 B CN 105353615B CN 201510776314 A CN201510776314 A CN 201510776314A CN 105353615 B CN105353615 B CN 105353615B
- Authority
- CN
- China
- Prior art keywords
- mrow
- mtd
- msub
- msup
- mover
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 40
- 230000009466 transformation Effects 0.000 claims abstract description 8
- 238000002347 injection Methods 0.000 claims abstract description 7
- 239000007924 injection Substances 0.000 claims abstract description 7
- 239000011159 matrix material Substances 0.000 claims description 34
- RZVHIXYEVGDQDX-UHFFFAOYSA-N 9,10-anthraquinone Chemical compound C1=CC=C2C(=O)C3=CC=CC=C3C(=O)C2=C1 RZVHIXYEVGDQDX-UHFFFAOYSA-N 0.000 claims description 9
- 238000013461 design Methods 0.000 claims description 9
- 239000013256 coordination polymer Substances 0.000 claims 3
- 230000004044 response Effects 0.000 abstract description 7
- 238000009415 formwork Methods 0.000 abstract 2
- 238000004088 simulation Methods 0.000 description 7
- 230000008569 process Effects 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000006073 displacement reaction Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 238000005034 decoration Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000006798 recombination Effects 0.000 description 1
- 238000005215 recombination Methods 0.000 description 1
- 238000012360 testing method Methods 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
-
- 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/048—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 using a predictor
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)
- Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
- Feedback Control In General (AREA)
Abstract
本发明公开了一种基于滑模观测器的四旋翼飞行器的主动容错控制方法。考虑四旋翼飞行器在存在状态时滞情况下发生执行器故障,结合滑模观测器和滑模控制,提出一种主动容错控制方法。设计了滑模观测器,对系统进行线性变换,基于等效误差注入的思想对执行器故障进行重构,利用执行器故障的重构估计值,在滑模控制中加入补偿控制,最终构成完整的主动容错控制器。本发明方法通过设计滑模观测器,对故障进行重构和估计,可以实现控制器增益的在线调整,使得所提控制律为最优,有效地提高了四旋翼飞行器飞行的控制精度和响应速度,可为带有执行器故障的复杂四旋翼飞行器提供容错控制器设计依据。本发明用于带有时变时滞的四旋翼飞行器的主动容错控制。
Description
技术领域
本发明涉及一种基于滑模观测器的四旋翼飞行器的主动容错控制方法,属于飞行器故障诊断与容错控制领域。
背景技术
典型的四旋翼直升机,四个旋翼呈十字形分布在前、后、左、右四个方向。每个旋翼都连接一个直流电机,控制律通过调节直流电机的转速来控制旋翼的转速,从而调节四个方向顶点的升力。为了抵消旋翼旋转产生的反扭矩,前后为一组,左右为一组,两组旋翼按照相反的方向旋转,因此不需要像常规纵列式直升飞机借助尾桨来消除扭矩。四旋翼直升机是一个复杂的被控对象,具有多输入多输出以及非线性、强耦合、时滞等各种各样的复杂问题,且在飞行过程中会不可避免地遇到风扰、发动机振动等不确定因素,加之缺少人为实时操纵,直升机一旦发生故障,将会引起灾难性后果。因此容错控制器需要在系统存在时滞和不确定性的情况下仍然具备较强的容错能力。
目前对四旋翼飞行器的容错控制方法主要分为主动容错控制和被动容错控制,由于主动容错控制通过故障调节或信号重构,保证故障发生后系统的稳定性,该方法设计灵活,容错能力强,因此一般对四旋翼飞行器采用主动容错控制。主动容错控制方法一般有控制律重新调度方法、控制律重构方法和模型跟随重组方法等方法,然而这些方法却很难起到很好的容错效果,原因是真实的飞行器会存在时滞等多种非线性因素,这就使得现有方法很难很好地抑制故障的影响。
由于滑模控制的滑动模态对系统参数摄动和外加干扰有完全的自适应性,因此非常适合处理四旋翼直升机飞控系统的被动容错控制问题。它的控制是不连续的,控制过程中,闭环系统的结构不停的变化,迫使系统状态沿着预先设计好的滑模面运动,渐渐“滑”向状态平衡点,即渐近稳定。其最主要的优点是一旦系统状态量到达滑模面,系统便不受参数变化和外界扰动的影响。滑模控制广泛用于飞控系统中,为飞控系统的容错控制提供了新思路。
为了有效处理四旋翼飞行器飞控系统中存在的时滞和不确定性,叶思隽提出了一种鲁棒容错控制算法。孙新柱针对不满足匹配条件的不确定系统提出了一种可靠跟踪控制器,贾新春则研究了不确定线性时滞系统的可靠保性能问题。但现有方法多为结构固定单一的被动容错控制,对复杂的飞控系统很难有很好的控制效果,因此本发明有很好的实用性。
发明内容
发明目的:针对上述现有技术,提出一种基于滑模观测器的四旋翼飞行器的主动容错控制方法,设计滑模观测器对执行器故障进行重构和估计,构造线性变换,利用等效误差注入思想,得出故障的估计值,设计积分型滑模面,并利用故障估计值设计容错控制律,克服故障对系统的影响。
技术方案:一种基于滑模观测器的四旋翼飞行器的主动容错控制方法,其特征在于:考虑四旋翼飞行器存在时变时滞和执行器故障,结合滑模观测器和滑模控制,提出一种主动容错控制方法,使得飞行器在发生执行器故障后能够继续安全飞行。根据所获取的飞行器的模型参数,设计一种具有滑模补偿的观测器,对执行器故障进行重构和估计,进而设计相应的滑模面和滑模控制律,最终构成主动容错控制器。包括如下具体步骤:
步骤1)获取四旋翼飞行器的控制模型:
其中,x=[x1 x2 x3]T为系统状态变量,以四旋翼飞行器X轴方向的位置控制为例,分别表示X轴方向上的位置、速度和执行器动态,u(t)为控制输入,y为可测输出,h(t)为时变时滞,满足η是有界常数,f(x,t)为执行器故障项,满足||f(x,t)||≤M。
步骤2)针对以上具有时变时滞和执行器故障的四旋翼飞控系统,设计具有滑模补偿项的观测器:
其中,L∈Rn×m为待设计的观测器增益,v∈Rm是滑模切换项,矩阵CG是列满秩的,(A,G)可控。定义为状态估计误差,ey=Ce为输出误差,则由式(1)和式(2)可得误差系统的状态方程:
ey=Ce (4)
步骤2.1)构造观测器增益阵L使得A-LC稳定,则系统理想滑动模态ey=0,且
将(3)式代入得到等效控制为:
veq=(CG)-1C(Ace(t)+Ade(t-h(t))+Df(t)) (5)
将(5)式代入(3)式中,得到误差系统理想滑模的状态方程:
由于矩阵(I-(CG)-1C)具有m零特征值和n-m指定特征值,所以理想滑动模态是渐进稳定的。
步骤2.2)构造补偿控制器满足下式:
其中,ρ是标量函数,滑模切换项增益阵w满足:
如果存在常数ε1>0,ε2>0和正定矩阵P>0,使得如下线性矩阵不等式成立:
则误差系统(2.3)是渐进稳定的,即构造的时滞滑模观测器在系统出现故障的情况下实现了未知状态的估计。
步骤3)根据步骤2)中设计的观测器,定义如下线性变换,对故障进行重构。
T=[CT ⊥P C]T
其中,CT ⊥为CT的正交补矩阵,由误差方程(2.3)可以得到在新的坐标系下的误差方程和输出方程:
其中,
根据等效输出误差注入原理,执行器重构故障的估计值为:
步骤4)综合步骤2)和步骤3),设计完整的容错控制律:
步骤4.1)根据滑模控制器设计方法,首先设计滑模面:
其中,矩阵H满HB足非奇异,K是待定的常数矩阵。可以证明,在该滑模面上的系统滑动模态是渐进稳定的。
步骤4.2)设计如下所示的控制律:
u=ul+un (11)
其中,为滑模控制律的线性部分,用来维持系统在滑模面上的理想滑动模态运动。利用等效控制思想解出线性部分的控制律如下:
ul=-Kx(t) (12)
非线性部分需要知道故障项f(t)的上界信息,由于该信息是未知的,因此可以利用式(9)给出的故障项的估计值来进行控制律的设计:
式中,ε为某很小的正常数。结合式(12)和(13),可以得出完整的滑模容错控制律如下:
步骤5)根据四旋翼飞行器的飞行状态,选择合适的参数,完成对其的容错控制。有益效果:本发明提出的一种基于滑模观测器的四旋翼飞行器的主动容错控制方法,设计了滑模观测器,对系统进行线性变换,基于等效误差注入的思想对执行器故障进行重构,利用执行器故障的重构估计值,在滑模控制中加入补偿控制,最终构成完整的主动容错控制器。
具有如下优点:
(1)通过设计具有滑模切换项的滑模观测器,使得飞控系统具有更好的鲁棒性;
(2)利用线性矩阵不等式给出保证系统渐进稳定的时滞上界值,充分考虑到四旋翼飞行器在实际飞行过程中可能存在的时滞现象,使得控制器的设计具有更好的实用性;
(3)引入自适应边界估计的方法估计出四旋翼飞行器执行器故障的大小,容错控制律不断地改变参数,使得系统保守性更小,控制效果更佳。
本发明所用方法作为一种四旋翼飞行器的容错控制方法,具有一定的实际应用价值,易于实现,容错能力强,能够有效提高四旋翼飞行器的飞行安全性。该方法可操作性强,应用方便、可靠。
附图说明
图1是本发明方法的流程图;
图2是Quanser的四旋翼飞行器仿真实验系统;
图3是四旋翼飞行器姿态运动示意图;
图4是四旋翼飞行器控制系统原理框图;
图5是X轴位移响应曲线;
图6是X轴速度响应曲线;
图7是执行器动态响应曲线;
图8是控制输入响应曲线;
图9是simulink仿真图。
具体实施方式
下面结合附图对本发明做更进一步的解释。
如图1所示,考虑四旋翼飞行器在存在状态时滞情况下发生执行器故障,结合滑模观测器和滑模控制,提出一种主动容错控制方法,使得飞行器在发生执行器故障后能够继续安全飞行。设计了滑模观测器,对系统进行线性变换,基于等效误差注入的思想对执行器故障进行重构,利用执行器故障的重构估计值,在滑模控制中加入补偿控制,最终构成完整的主动容错控制器。包括如下具体步骤:
步骤1)获取四旋翼飞行器的控制模型:
其中,x=[x1 x2 x3]T为系统状态变量,以四旋翼飞行器X轴方向的位置控制为例,分别表示X轴方向上的位置、速度和执行器动态,u(t)为控制输入,y为可测输出,h(t)为时变时滞,满足η是有界常数,f(x,t)为执行器故障项,满足||f(x,t)||≤M。
步骤2)针对以上具有时变时滞和执行器故障的四旋翼飞控系统,设计具有滑模补偿项的观测器:
其中,L∈Rn×m为待设计的观测器增益,v∈Rm是滑模切换项,矩阵CG是列满秩的,(A,G)可控。定义为状态估计误差,ey=Ce为输出误差,则由式(1)和式(2)可得误差系统的状态方程:
ey=Ce (4)
步骤2.1)构造观测器增益阵L使得A-LC稳定,则系统理想滑动模态ey=0,且
将(3)式代入得到等效控制为:
veq=(CG)-1C(Ace(t)+Ade(t-h(t))+Df(t)) (5)
将(5)式代入(3)式中,得到误差系统理想滑模的状态方程:
由于矩阵(I-(CG)-1C)具有m零特征值和n-m指定特征值,所以理想滑动模态是渐进稳定的。
步骤2.2)构造补偿控制器满足下式:
G=P-1CTw-1
其中,ρ是标量函数,滑模切换项增益阵w满足:
如果存在常数ε1>0,ε2>0和正定矩阵P>0,使得如下线性矩阵不等式成立:
则误差系统(2.3)是渐进稳定的,即构造的时滞滑模观测器在系统出现故障的情况下实现了未知状态的估计。
步骤3)根据步骤2)中设计的观测器,定义如下线性变换,对故障进行重构。
T=[CT ⊥P C]T
其中,CT ⊥为CT的正交补矩阵,由误差方程(2.3)可以得到在新的坐标系下的误差方程和输出方程:
其中,
根据等效输出误差注入原理,执行器重构故障的估计值为:
步骤4)综合步骤2)和步骤3),设计完整的容错控制律:
步骤4.1)根据滑模控制器设计方法,首先设计滑模面:
其中,矩阵H满HB足非奇异,K是待定的常数矩阵。可以证明,在该滑模面上的系统滑动模态是渐进稳定的。
步骤4.2)设计如下所示的控制律:
u=ul+un (11)
其中,为滑模控制律的线性部分,用来维持系统在滑模面上的理想滑动模态运动。利用等效控制思想解出线性部分的控制律如下:
ul=-Kx(t) (12)
非线性部分需要知道故障项f(t)的上界信息,由于该信息是未知的,因此可以利用式(9)给出的故障项的估计值来进行控制律的设计:
式中,ε为某很小的正常数。结合式(12)和(13),可以得出完整的滑模容错控制律如下:
步骤5)根据四旋翼飞行器的飞行状态,选择合适的参数,完成对其的容错控制。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围,
下面以实际案例仿真说明实施方案的有效性。
采用加拿大quanser公司生产的Qball-X4四旋翼飞行器半物理仿真平台作为具体的算法实验仿真对象。图2是Quanser的四旋翼飞行器仿真实验系统,图3是四旋翼飞行器姿态运动示意图。该仿真实验系统由地面控制站、照相机定位系统和飞行器组成,主控制机通过无线局域网与各个无人工具进行通讯,主要对系统进行定位和任务规划。一旦整个控制系统的控制算法设计完成,可以使控制站仅仅起到定位作用,从而进行无人工具的自主控制及多个工具之间的协调控制研究。系统通过六个红外照相机实现空间三维定位,从而获取所需参数。
四旋翼飞行器的数学模型如下所示:
其中,各系数矩阵如下:
根据实验测量,本试验系统时滞波动范围为0~2s,实验中τ设为1s。实验中,四旋翼飞行器在飞行过程中遇到如下形式的突发故障:
取初始时刻系统的状态量矢量为:
x0=[x1 x2 x3]T=[1 1 0.25]T
根据本发明方法,对发生执行器故障的四旋翼飞行器进行容错控制。根据步骤1)-步骤5),其中待定的参数取值如下:观测器增益阵矩阵P由线性矩阵不等式(7)解出,线性变换矩阵T=[CT ⊥P C]T,其中C=[1 0 0],H=[0 0 1],连续部分反馈增益阵K=[2.0556 0.8722 0.8]。
图5-图8为容错控制结果。图5-图7分别是X轴方向位移、速度和执行器动态的响应曲线,图8是控制输入的曲线。
由图5-图7可知,当系统发生执行器故障后,在本发明的容错控制下,飞行器X轴位移和速度均能在较短的时间内趋于稳定,且响应速度快,超调小,也就是说当系统发生故障之后飞行器仍然能够维持原来的飞行状态,避免了事故的发生。由图8的控制曲线可知,即使系统存在时滞和执行器故障,滑模容错控制律变化幅度很小,即该容错控制律能够很好地保证飞行器的控制精度和安全性。
Claims (1)
1.一种基于滑模观测器的四旋翼飞行器的主动容错控制方法,其特征在于:考虑四旋翼飞行器存在时变时滞和执行器故障,结合滑模观测器和滑模控制,提出一种主动容错控制方法,使得飞行器在发生执行器故障后能够继续安全飞行,根据所获取的飞行器的模型参数,设计一种具有滑模补偿的观测器,对执行器故障进行重构和估计,进而设计相应的滑模面和滑模控制律,最终构成主动容错控制器,包括如下具体步骤:
步骤1)获取四旋翼飞行器的控制模型:
<mrow>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>A</mi>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>A</mi>
<mi>d</mi>
</msub>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mi>h</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>B</mi>
<mi>u</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>D</mi>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
y(t)=Cx(t) (1)
其中,x=[x1 x2 x3]T为系统状态变量,以四旋翼飞行器X轴方向的位置控制为例,分别表示X轴方向上的位置、速度和执行器动态,t为时间变量,A∈Rn×n,Ad∈Rn×n,B∈Rn×m,C∈Rp ×n为系统常值矩阵,D∈Rn×m为故障增益矩阵,u(t)为控制输入,y为可测输出,h(t)为时变时滞,满足η是有界常数,f(t)为执行器故障项,满足||f(t)||≤M,M是一个正的常数;
步骤2)针对以上具有时变时滞和执行器故障的四旋翼飞控系统,设计具有滑模补偿项的观测器:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>A</mi>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>A</mi>
<mi>d</mi>
</msub>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mi>h</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>B</mi>
<mi>u</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>L</mi>
<mrow>
<mo>(</mo>
<mi>y</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>-</mo>
<mi>C</mi>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>G</mi>
<mi>v</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>C</mi>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,L∈Rn×m为待设计的观测器增益,G∈Rn×m为待设计的滑模增益矩阵,v∈Rm是滑模切换项,矩阵CG是列满秩的,(A,G)可控,定义为状态估计误差,ey=Ce为输出误差,则由式(1)和式(2)可得误差系统的状态方程:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mover>
<mi>e</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mi>A</mi>
<mo>-</mo>
<mi>L</mi>
<mi>C</mi>
<mo>)</mo>
</mrow>
<mi>e</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>A</mi>
<mi>d</mi>
</msub>
<mi>e</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mi>h</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>D</mi>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>G</mi>
<mi>v</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<msub>
<mi>A</mi>
<mi>c</mi>
</msub>
<mi>e</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>A</mi>
<mi>d</mi>
</msub>
<mi>e</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mi>h</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>D</mi>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>G</mi>
<mi>v</mi>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
ey=Ce (4)
其中,Ac=A-LC,且Ac∈Rn×n为误差系统矩阵;
步骤2.1)构造观测器增益阵L使得A-LC稳定,则系统理想滑动模态ey=0,且
将(3)式代入得到等效控制为:
veq=(CG)-1C(Ace(t)+Ade(t-h(t))+Df(t)) (5)
将(5)式代入(3)式中,得到误差系统理想滑模的状态方程:
<mrow>
<mover>
<mi>e</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mi>I</mi>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mi>C</mi>
<mi>G</mi>
</mrow>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mi>C</mi>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>A</mi>
<mi>c</mi>
</msub>
<mi>e</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>+</mo>
<msub>
<mi>A</mi>
<mi>d</mi>
</msub>
<mi>e</mi>
<mo>(</mo>
<mrow>
<mi>t</mi>
<mo>-</mo>
<mi>h</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
<mo>+</mo>
<mi>D</mi>
<mi>f</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,I表示单位矩阵,由于矩阵(I-(CG)-1C)具有m个零特征值和n-m个指定特征值,所以理想滑动模态是渐进稳定的;
步骤2.2)构造补偿控制器满足下式:
<mrow>
<mi>v</mi>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>&rho;</mi>
<mi>w</mi>
<mfrac>
<mrow>
<mi>C</mi>
<mi>e</mi>
</mrow>
<mrow>
<mo>|</mo>
<mo>|</mo>
<mi>C</mi>
<mi>e</mi>
<mo>|</mo>
<mo>|</mo>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>y</mi>
<mo>&NotEqual;</mo>
<mover>
<mi>y</mi>
<mo>^</mo>
</mover>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>y</mi>
<mo>=</mo>
<mi>C</mi>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中,ρ是标量函数,滑模切换项增益阵w满足:
<mrow>
<msub>
<mi>&lambda;</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>w</mi>
<mo>)</mo>
</mrow>
<mo>&GreaterEqual;</mo>
<mi>M</mi>
<mfrac>
<mrow>
<msub>
<mi>&lambda;</mi>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msup>
<mi>CP</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msup>
<mi>C</mi>
<mi>T</mi>
</msup>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>&lambda;</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msup>
<mi>CP</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msup>
<mi>C</mi>
<mi>T</mi>
</msup>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
其中,λmax,λmin分别指矩阵的最大特征值和最小特征值,如果存在常数ε1>0,ε2>0和正定矩阵P>0,使得如下线性矩阵不等式成立:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<msub>
<mi>A</mi>
<mi>c</mi>
</msub>
<mi>T</mi>
</msup>
<mi>P</mi>
<mo>+</mo>
<msub>
<mi>PA</mi>
<mi>c</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<msup>
<msub>
<mi>A</mi>
<mi>d</mi>
</msub>
<mi>T</mi>
</msup>
</mrow>
</mtd>
<mtd>
<mi>P</mi>
</mtd>
<mtd>
<mi>K</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mfrac>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mi>&eta;</mi>
</mrow>
<msub>
<mi>&epsiv;</mi>
<mn>1</mn>
</msub>
</mfrac>
<mi>I</mi>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mrow>
<msub>
<mi>&epsiv;</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msup>
<msub>
<mi>&epsiv;</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
</mrow>
</mfrac>
<mi>I</mi>
</mrow>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<mo>*</mo>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>&epsiv;</mi>
<mn>2</mn>
</msub>
</mfrac>
<mi>I</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo><</mo>
<mn>0</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,K为待定的常数矩阵,则误差系统(2.3)是渐进稳定的,即构造的时滞滑模观测器在系统出现故障的情况下实现了未知状态的估计;
步骤3)根据步骤2)中设计的观测器,定义如下线性变换,对故障进行重构,
T=[CT 丄P C]T
其中,CT 丄为CT的正交补矩阵,由误差方程(3)可以得到在新的坐标系下的误差方程和输出方程:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mover>
<mover>
<mi>e</mi>
<mo>&OverBar;</mo>
</mover>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<msub>
<mover>
<mi>A</mi>
<mo>&OverBar;</mo>
</mover>
<mi>c</mi>
</msub>
<mo>+</mo>
<msub>
<mover>
<mi>A</mi>
<mo>&OverBar;</mo>
</mover>
<mi>d</mi>
</msub>
<mover>
<mi>e</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mi>h</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mover>
<mi>D</mi>
<mo>&OverBar;</mo>
</mover>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>T</mi>
<mi>G</mi>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>e</mi>
<mi>y</mi>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mi>I</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mover>
<mi>e</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,
<mrow>
<mover>
<mi>e</mi>
<mo>&OverBar;</mo>
</mover>
<mo>=</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>e</mi>
<mn>1</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>e</mi>
<mi>y</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
<mo>=</mo>
<mi>T</mi>
<mi>e</mi>
<mo>,</mo>
<msub>
<mover>
<mi>A</mi>
<mo>&OverBar;</mo>
</mover>
<mi>c</mi>
</msub>
<mo>=</mo>
<msub>
<mi>TA</mi>
<mi>c</mi>
</msub>
<msup>
<mi>T</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>A</mi>
<mrow>
<mi>c</mi>
<mn>11</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>A</mi>
<mrow>
<mi>c</mi>
<mn>12</mn>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>A</mi>
<mrow>
<mi>c</mi>
<mn>21</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>A</mi>
<mrow>
<mi>c</mi>
<mn>22</mn>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<msub>
<mover>
<mi>A</mi>
<mo>&OverBar;</mo>
</mover>
<mi>d</mi>
</msub>
<mo>=</mo>
<msub>
<mi>TA</mi>
<mi>d</mi>
</msub>
<msup>
<mi>T</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>A</mi>
<mrow>
<mi>d</mi>
<mn>11</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>A</mi>
<mrow>
<mi>d</mi>
<mn>12</mn>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>A</mi>
<mrow>
<mi>d</mi>
<mn>21</mn>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>A</mi>
<mrow>
<mi>d</mi>
<mn>22</mn>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
<mrow>
<mover>
<mi>D</mi>
<mo>&OverBar;</mo>
</mover>
<mo>=</mo>
<mi>T</mi>
<mi>D</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>C</mi>
<msup>
<mi>P</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msup>
<mi>C</mi>
<mi>T</mi>
</msup>
<msup>
<mi>&Lambda;</mi>
<mi>T</mi>
</msup>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mi>T</mi>
<mi>G</mi>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>&rho;</mi>
<mi>C</mi>
<msup>
<mi>P</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msup>
<mi>C</mi>
<mi>T</mi>
</msup>
<mfrac>
<mrow>
<mi>C</mi>
<mi>e</mi>
</mrow>
<mrow>
<mo>|</mo>
<mo>|</mo>
<mi>C</mi>
<mi>e</mi>
<mo>|</mo>
<mo>|</mo>
</mrow>
</mfrac>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>v</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中,Λ∈Rm×p为待定常数矩阵,根据等效输出误差注入原理,执行器重构故障的估计值为:
<mrow>
<mover>
<mi>f</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&ap;</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>C</mi>
<mi>D</mi>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<mi>C</mi>
<mi>D</mi>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msup>
<mrow>
<mo>(</mo>
<mi>C</mi>
<mi>D</mi>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<mi>C</mi>
<mi>v</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>C</mi>
<mi>D</mi>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<mi>C</mi>
<mi>D</mi>
<mo>&rsqb;</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msup>
<mrow>
<mo>(</mo>
<mi>C</mi>
<mi>D</mi>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<msup>
<mi>CP</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<msup>
<mi>C</mi>
<mi>T</mi>
</msup>
<mfrac>
<msub>
<mi>e</mi>
<mi>y</mi>
</msub>
<mrow>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>e</mi>
<mi>y</mi>
</msub>
<mo>|</mo>
<mo>|</mo>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤4)综合步骤2)和步骤3),设计完整的容错控制律:
步骤4.1)根据滑模控制器设计方法,首先设计滑模面:
<mrow>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>H</mi>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msubsup>
<mo>&Integral;</mo>
<mn>0</mn>
<mi>t</mi>
</msubsup>
<mi>H</mi>
<mrow>
<mo>(</mo>
<mi>A</mi>
<mo>-</mo>
<mi>B</mi>
<mi>K</mi>
<mo>)</mo>
</mrow>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>s</mi>
<mo>-</mo>
<msubsup>
<mo>&Integral;</mo>
<mn>0</mn>
<mi>t</mi>
</msubsup>
<mi>H</mi>
<msub>
<mi>A</mi>
<mi>d</mi>
</msub>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>-</mo>
<mi>h</mi>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>s</mi>
<mo>-</mo>
<mi>H</mi>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,矩阵H满足矩阵HB非奇异,K是待定的常数矩阵,可以证明,在该滑模面上的系统滑动模态是渐进稳定的;
步骤4.2)设计如下所示的控制律:
u=ul+un (11)
其中,ul为滑模控制律的线性部分,用来维持系统在滑模面上的理想滑动模态运动,利用等效控制思想解出线性部分的控制律如下:
ul=-Kx(t) (12)
非线性部分需要知道故障项f(t)的上界信息,由于该信息是未知的,因此可以利用式(9)给出的故障项的估计值来进行控制律的设计:
<mrow>
<msub>
<mi>u</mi>
<mi>n</mi>
</msub>
<mo>=</mo>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>H</mi>
<mi>B</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mo>|</mo>
<mo>|</mo>
<mover>
<mi>f</mi>
<mo>^</mo>
</mover>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>|</mo>
<mo>|</mo>
<mo>+</mo>
<mi>&epsiv;</mi>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<mo>|</mo>
<mo>|</mo>
<mi>H</mi>
<mi>B</mi>
<mo>|</mo>
<mo>|</mo>
<mfrac>
<mi>s</mi>
<mrow>
<mo>|</mo>
<mo>|</mo>
<mi>s</mi>
<mo>|</mo>
<mo>|</mo>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>13</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,ε为某正常数,结合式(12)和(13),可以得出完整的滑模容错控制律如下:
<mrow>
<mi>u</mi>
<mo>=</mo>
<msub>
<mi>u</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>u</mi>
<mi>n</mi>
</msub>
<mo>=</mo>
<mo>-</mo>
<mi>K</mi>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>H</mi>
<mi>B</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mo>|</mo>
<mo>|</mo>
<mover>
<mi>f</mi>
<mo>^</mo>
</mover>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>|</mo>
<mo>|</mo>
<mo>+</mo>
<mi>&epsiv;</mi>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<mo>|</mo>
<mo>|</mo>
<mi>H</mi>
<mi>B</mi>
<mo>|</mo>
<mo>|</mo>
<mfrac>
<mi>s</mi>
<mrow>
<mo>|</mo>
<mo>|</mo>
<mi>s</mi>
<mo>|</mo>
<mo>|</mo>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>14</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤5)根据四旋翼飞行器的飞行状态,选择合适的参数,完成对其的容错控制。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510776314.5A CN105353615B (zh) | 2015-11-10 | 2015-11-10 | 一种基于滑模观测器的四旋翼飞行器的主动容错控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510776314.5A CN105353615B (zh) | 2015-11-10 | 2015-11-10 | 一种基于滑模观测器的四旋翼飞行器的主动容错控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105353615A CN105353615A (zh) | 2016-02-24 |
CN105353615B true CN105353615B (zh) | 2018-02-09 |
Family
ID=55329604
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510776314.5A Expired - Fee Related CN105353615B (zh) | 2015-11-10 | 2015-11-10 | 一种基于滑模观测器的四旋翼飞行器的主动容错控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105353615B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107861383B (zh) * | 2017-10-23 | 2021-07-09 | 天津大学 | 基于自适应观测器的卫星故障诊断与容错控制方法 |
Families Citing this family (28)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106055770B (zh) * | 2016-05-26 | 2019-02-26 | 南京航空航天大学 | 一种基于滑模理论的航空发动机气路故障诊断方法 |
CN105911866B (zh) * | 2016-06-15 | 2019-01-22 | 浙江工业大学 | 四旋翼无人飞行器的有限时间全阶滑模控制方法 |
CN105912012B (zh) * | 2016-06-29 | 2018-09-11 | 哈尔滨工业大学 | 一种部分执行机构故障下航天器混沌姿态控制方法 |
CN106444799B (zh) * | 2016-07-15 | 2019-03-12 | 浙江工业大学 | 基于模糊扩张状态观测器和自适应滑模的四旋翼无人机控制方法 |
CN106597851B (zh) * | 2016-12-15 | 2019-04-30 | 南京航空航天大学 | 一种小型无人机飞控系统的鲁棒容错控制方法 |
CN106406102B (zh) * | 2016-12-20 | 2019-06-11 | 哈尔滨工业大学 | 一种含干扰观测器的高超声速飞行器跟踪控制方法 |
CN106842920B (zh) * | 2017-01-04 | 2019-04-30 | 南京航空航天大学 | 针对多时滞四旋翼直升机飞控系统的鲁棒容错控制方法 |
CN106774273B (zh) * | 2017-01-04 | 2019-01-01 | 南京航空航天大学 | 针对变时滞控制系统执行器故障的滑模预测容错控制方法 |
CN106647584A (zh) * | 2017-01-17 | 2017-05-10 | 南京航空航天大学 | 一种基于最优滑模的四旋翼飞行器的容错控制方法 |
CN108345212A (zh) * | 2017-01-24 | 2018-07-31 | 南京航空航天大学 | 一种基于滑模的三自由度直升机的鲁棒h∞控制方法 |
CN107256028B (zh) * | 2017-07-24 | 2020-05-12 | 大连理工大学 | 四旋翼飞行器对角动力损失状态下的失控保护控制算法 |
CN109491244B (zh) * | 2017-09-13 | 2021-10-22 | 南京航空航天大学 | 一种基于滑模观测器的无人机编队系统故障诊断方法 |
CN108180910B (zh) * | 2017-12-26 | 2019-01-08 | 北京航空航天大学 | 一种基于气动参数不确定的飞行器快速高精度制导方法 |
CN108398884A (zh) * | 2018-03-09 | 2018-08-14 | 南京航空航天大学 | 一种基于滑模的不确定时滞系统的自适应容错控制方法 |
CN110161847B (zh) * | 2018-09-18 | 2020-12-08 | 南京航空航天大学 | 一种基于分布式奇异观测器的无人机编队系统传感器故障估计方法 |
WO2020142984A1 (zh) * | 2019-01-10 | 2020-07-16 | 大连理工大学 | 一种基于误差区间观测器的航空发动机主动容错控制方法 |
CN109630281B (zh) * | 2019-01-10 | 2020-12-11 | 大连理工大学 | 一种基于误差区间观测器的航空发动机主动容错控制方法 |
CN110119089B (zh) * | 2019-03-29 | 2022-06-14 | 华东理工大学 | 一种基于积分滑模的浸入不变流型自适应四旋翼控制方法 |
CN110320799B (zh) * | 2019-06-13 | 2021-05-07 | 大连理工大学 | 一种面向航空发动机执行器故障的无扰切入容错控制方法 |
CN110262248B (zh) * | 2019-07-08 | 2021-05-11 | 南京航空航天大学 | 一种微型燃气轮机故障鲁棒自适应重构方法 |
CN110531778B (zh) * | 2019-07-17 | 2020-10-02 | 北京航空航天大学 | 一种多旋翼无人机自驾仪桨叶损伤的估计与自愈控制方法 |
CN111522241B (zh) * | 2020-05-08 | 2020-12-29 | 哈尔滨工业大学 | 一种基于固定时间观测器的主动容错控制方法及装置 |
CN111679579B (zh) * | 2020-06-10 | 2021-10-12 | 南京航空航天大学 | 针对传感器和执行器故障系统的滑模预测容错控制方法 |
CN111880407B (zh) * | 2020-07-16 | 2022-08-02 | 中国人民解放军海军航空大学岸防兵学院 | 一种非线性系统的稳定方法及装置 |
CN112147897B (zh) * | 2020-09-28 | 2021-10-12 | 南京航空航天大学 | 离散不确定多时滞四旋翼系统的容错控制方法 |
CN112180960B (zh) * | 2020-09-29 | 2021-09-14 | 西北工业大学 | 一种面向执行器故障的无人机容错飞行方法及飞行系统 |
CN115877717B (zh) * | 2022-12-26 | 2023-06-13 | 南通大学 | 一种基于自抗扰控制的飞行器容错控制结构与控制方法 |
CN117784616B (zh) * | 2024-02-23 | 2024-05-24 | 西北工业大学 | 一种基于智能观测器组的高速飞行器故障重构方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103105850A (zh) * | 2013-01-30 | 2013-05-15 | 南京航空航天大学 | 一种近空间飞行器故障诊断与容错控制方法 |
CN103135553A (zh) * | 2013-01-21 | 2013-06-05 | 南京航空航天大学 | 四旋翼飞行器容错控制方法 |
CN103838145A (zh) * | 2014-01-17 | 2014-06-04 | 南京航空航天大学 | 基于级联观测器的垂直起降飞机鲁棒容错控制系统及方法 |
CN104020774A (zh) * | 2014-06-19 | 2014-09-03 | 金陵科技学院 | 基于动态控制重新分配的飞行器姿态容错控制方法 |
CN104035447A (zh) * | 2014-06-27 | 2014-09-10 | 金陵科技学院 | 基于动态控制重新分配的无人飞行器姿态容错控制方法 |
CN104238357A (zh) * | 2014-08-21 | 2014-12-24 | 南京航空航天大学 | 一种近空间飞行器的容错滑模控制方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080133074A1 (en) * | 2006-12-01 | 2008-06-05 | Zyss Michael J | Autonomous rollout control of air vehicle |
US8346408B2 (en) * | 2008-11-10 | 2013-01-01 | The Boeing Company | Fault tolerant flight control system |
-
2015
- 2015-11-10 CN CN201510776314.5A patent/CN105353615B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103135553A (zh) * | 2013-01-21 | 2013-06-05 | 南京航空航天大学 | 四旋翼飞行器容错控制方法 |
CN103105850A (zh) * | 2013-01-30 | 2013-05-15 | 南京航空航天大学 | 一种近空间飞行器故障诊断与容错控制方法 |
CN103838145A (zh) * | 2014-01-17 | 2014-06-04 | 南京航空航天大学 | 基于级联观测器的垂直起降飞机鲁棒容错控制系统及方法 |
CN104020774A (zh) * | 2014-06-19 | 2014-09-03 | 金陵科技学院 | 基于动态控制重新分配的飞行器姿态容错控制方法 |
CN104035447A (zh) * | 2014-06-27 | 2014-09-10 | 金陵科技学院 | 基于动态控制重新分配的无人飞行器姿态容错控制方法 |
CN104238357A (zh) * | 2014-08-21 | 2014-12-24 | 南京航空航天大学 | 一种近空间飞行器的容错滑模控制方法 |
Non-Patent Citations (1)
Title |
---|
基于滑模与控制分配的六旋翼飞行器容错控制;陈阳 等;《电光与控制》;20140531;第21卷(第5期);第24-28页 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107861383B (zh) * | 2017-10-23 | 2021-07-09 | 天津大学 | 基于自适应观测器的卫星故障诊断与容错控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105353615A (zh) | 2016-02-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105353615B (zh) | 一种基于滑模观测器的四旋翼飞行器的主动容错控制方法 | |
CN106774373B (zh) | 一种四旋翼无人机有限时间姿态跟踪控制方法 | |
Wang et al. | A composite adaptive fault-tolerant attitude control for a quadrotor UAV with multiple uncertainties | |
Peng et al. | Modeling and robust backstepping sliding mode control with adaptive RBFNN for a novel coaxial eight-rotor UAV | |
CN106325291B (zh) | 基于滑模控制律和eso的四旋翼飞行器姿态控制方法及系统 | |
CN103838145B (zh) | 基于级联观测器的垂直起降飞机鲁棒容错控制系统及方法 | |
CN103135553B (zh) | 四旋翼飞行器容错控制方法 | |
CN105607473B (zh) | 小型无人直升机的姿态误差快速收敛自适应控制方法 | |
Gui et al. | Adaptive fault-tolerant spacecraft pose tracking with control allocation | |
CN107272639B (zh) | 刚性航天器反作用飞轮故障的检测、估计及其调节方法 | |
CN108490783B (zh) | 基于增强型双幂次趋近律和快速终端滑模面的刚性航天飞行器有限时间自适应容错控制方法 | |
CN108398884A (zh) | 一种基于滑模的不确定时滞系统的自适应容错控制方法 | |
Zou et al. | Adaptive trajectory tracking controller for quadrotor systems subject to parametric uncertainties | |
CN106647584A (zh) | 一种基于最优滑模的四旋翼飞行器的容错控制方法 | |
CN106707749B (zh) | 一种针对仿生扑翼飞行机器人的控制方法 | |
CN109634291B (zh) | 一种基于改进型障碍李雅普诺夫函数的刚性飞行器姿态约束跟踪控制方法 | |
CN109116737A (zh) | 基于观测器的多智能体系统故障上限评估以及容错一致性控制算法 | |
CN107422741A (zh) | 基于学习的保预设性能集群飞行分布式姿态追踪控制方法 | |
CN102183957A (zh) | 船舶航向变论域模糊与最小二乘支持向量机复合控制方法 | |
Sun et al. | Hybrid sensor-based backstepping control approach with its application to fault-tolerant flight control | |
CN105549607B (zh) | 一种卫星姿态控制系统故障可重构的执行器构型设计方法 | |
Sun et al. | Practical solution to attitude control within wide envelope | |
CN108762069A (zh) | 刚性航天器姿控系统飞轮故障辨识与调节方法 | |
Ortiz-Torres et al. | An actuator fault detection and isolation method design for planar vertical take-off and landing unmanned aerial vehicle modelled as a qLPV system | |
Liu et al. | Fault diagnosis and accommodation for multi-actuator faults of a fixed-wing unmanned aerial vehicle |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180209 Termination date: 20211110 |
|
CF01 | Termination of patent right due to non-payment of annual fee |