高超声速飞行器Terminal滑模控制器设计方法技术

本发明专利技术涉及高超声速飞行器Terminal滑模控制器设计方法，Terminal滑模控制的优点是可以使在滑模面上的系统可以在有限时间内收敛。为了保持模型本身的非线性，并且给控制系统设计提供方便，文中采用输入输出线性化方法对模型进行了处理，之后，根据Terminal滑模控制设计的要求重新选择状态变量并设计了非线性滑模面。之后为设计的Terminal非线性滑模面设计了滑模控制律，保证了系统能够到达滑模面，并且证明了系统的稳定性。最后对所给出的控制律进行了仿真验证，结果表明所设计的滑模控制器能够实现对高超声速飞行器的有效控制。

Design method of Terminal sliding mode controller for hypersonic vehicle

The invention relates to a method for designing the Terminal sliding mode controller of a hypersonic vehicle. The advantage of the Terminal sliding mode control is that the system on the sliding surface can converge in finite time. In order to maintain the nonlinear model itself, and to the design of control system to provide convenient, after the input output linearization method of the model were processed, according to the design requirements of Terminal sliding mode control to re select the state variables and design a nonlinear sliding surface. Then, the sliding mode control law is designed for the Terminal nonlinear sliding surface, which guarantees that the system can reach the sliding mode surface and proves the stability of the system. Finally, the simulation results show that the proposed sliding mode controller can effectively control the hypersonic vehicle.

【技术实现步骤摘要】

本专利技术涉及高超声速飞行器Terminal滑模控制器设计方法，属于航空控制器调节

技术介绍
高超声速飞行器因其飞行速度快、机动性强、有效载荷大等优点受到世界各国的重视，成为下一代飞行器研发的重点。但由于飞行器本身采用机身发动机一体化设计，导致飞行器发动机状态与飞行状态互相耦合；并且高超声速飞行器的飞行包线大，飞行情况复杂，其整体动力学模型存在严重的非线性，同时结构与参数的不确定性以及外部干扰也不可忽视，这使得经典的控制器设计方法难以设计出符合其强鲁棒性、快速响应性等性能要求的控制器。针对高超声速飞行器的控制器设计问题，国内外展开了广泛的研究，目前采用的方法有H∞最优控制、模型参考自适应、线性变参数、自适应滑模控制等，高超声速飞行器控制器的设计吸引了世界很多学者的目光。滑模变结构控制因其本身对匹配不确定性以及外部干扰等的完全鲁棒性而受到欢迎。针对高超声速飞行器的纵向模型设计了自适应滑模控制器，证明了滑模控制的有效性，但采用的滑模面为传统的线性滑模面，其状态随时间指数收敛，收敛速度难以满足飞行器在高速飞行的条件下对响应快速性的要求。
技术实现思路
本方法将针对高超声速飞行器的指令跟踪问题进行研究，首先利用输入输出线性化的理论将飞行器动态模型进行线性化，采用Terminal滑模设计方法设计控制器，所设计滑模面为非线性的，使得系统的响应速度快，可以保证跟踪误差系统在有限时间内到达零点，从而保证了对指令信号的快速跟踪。选择线性的滑模面能够确保系统轨迹在达到滑动模态阶段以后，滑动模态的运动是渐进稳定的或者跟踪误差渐进的收敛到零，并且收敛速度是通过改变滑模面的参数矩阵来调节的，但是，线性滑模面的状态跟踪误差不能够在有限的时间内收敛到零。Terminal滑模的控制思想，通过将非线性项引入到滑模面的设计中，突破了以往线性滑模面的局限，使得系统的滑动模态变量能够在“有限的时间内”收敛至平衡点。最初Terminal滑动模态是从Terminal吸引子而来。Venkatar-manS.T.等首先分析了Terminal滑模的设计问题，并将其应用到机器人系统当中，随后YuX.和man.Z.H.等针对SISO系统和MIMO系统，对Terminal滑模的控制问题进行了详尽的研究。Wu.Y.Q.等还利用递归的方法深入讨论了Terminal滑模控制中容易出现的奇异问题。在本文中，我们将借鉴Terminal滑模的思想设计控制器，以保证高超声速飞行器的状态跟踪误差可以在有限时间内收敛。首先介绍一下Terminal的基本思想。一阶Terminal滑模中，一阶Terminal滑模定义如下：s=x·+βxq/p=0---(1)]]>其中x为一标量，β＞0,p,q(p＞q)，且p,q为正奇数。无论x为任何实数，xp/q的解必须为实数。系统在滑模面上的动态性能为x·=-βxq/p---(2)]]>给定任一初始状态x(0)≠0，则系统将在有限时间内收敛到原点。求解方程(2)pp-q[x(ts)(p-q)/p-x(0)(p-q)/p]=-βts---(3)]]>可得系统从状态x(0)到原点经历的时间ts=pβ(p-q)|x(0)|(p-q)/p---(4)]]>下面考虑Terminal有限机理的解释，考虑在平衡点x＝0处的Jacobian矩阵为J=∂x·∂x=-βqpx(p-q)/p---(5)]]>把J看做一阶矩阵的特征值λ，那么当x→0+时，J→-∞，那么系统的轨迹在负无穷大的特征值的驱动下，自然会以无穷大的速度收敛至平衡点。因此系统将在有限时间内收敛到原点。对于非线性系统，采用基于微分几何理论的输入输出线性化方法，可以在保留系统非线性特点的基础上降低系统控制器设计的复杂度，因此我们考虑采用输入输出线性化方法对系统进行处理。下面简要介绍一下该方法的原理。针对仿射非线性系统x·=f(x)+G(x)u,G(x)=[g1(x),g2(x),...,gm(x)]y=h(x)---(6)]]>其中f(x),g(x)为光滑函数。首先通过分别对函数f(x),g(x)求取关于输出函数h(x)系统的李导数，求系统的相对阶，具体形式为：Lfhi=∂hi∂xfLfkhi=Lf(Lfk-1hi)Lgihi=∂hi∂xgi---(7)]]>如果则称系统的相对阶为r。如果系统的相对阶r＝n，(n为系统阶数)，则系统为可以完全输入输出线性化的。选择微分同胚变换ζ1=Lfh(x)ζ2=Lf2h(x)...ζn-1=Lfr-1h(x)---(8)]]>对(8)式求导，可得变换系统的状态方程为：ζ·1=Lf2h(x)=ζ2ζ·2=Lf3h(x)=ζ3...ζ·n-1=Lfr-1h(x)=ζnζ·n=Lfrh(x)+Lgr-1h(x)u---(9)]]>观察方程(4)，可以发现除了最后关于ζn的方程外，其余n-1个方程已经是线性形式，并且不含控制量。只有关于ζn的状态方程是非线性的，但对输入u，方程形式上是线性的。重写状态方程(9)ζ·1=Lf2h(x)=ζ2ζ·2=Lf3h(x)=ζ3...ζ·n-1=Lfr-1h(x)=ζnζ·n=v---(10)]]>其中v=Lfrh(x)+Lgr-1h(x)u---(11)]]>则状态方程可以变为：ζ·=Aζ+Bv---(12)]]>其中这时，系统形式上变为线性，并且保留了系统的非线性特性，使得系统更加易于处理。本方法采用NASA兰利实验室公布的高超声速飞行器纵向模型进行研究，其模型如下：模型假设(1)高超声速飞行器为理想的刚体，即不考虑机翼等的弹性自由度；(2)质心位置，转动惯量是质量的函数，质心位置始终在机体轴纵轴变动；(3)飞行器中心和参考力矩中心在机体X轴上；(4)假设飞行器布局是对称的，也即惯性积Ixy，Ixz，Iyz恒为零；(5)忽略操纵面的转动惯量和发动机推力安装角。V&Cent本文档来自技高网...

【技术保护点】
高超声速飞行器Terminal滑模控制器设计方法，其特征在于：一阶Terminal滑模中，一阶Terminal滑模定义如下：其中x为一标量，β＞0,p,q(p＞q)，且p,q为正奇数；无论x为任何实数，xp/q的解必须为实数；系统在滑模面上的动态性能为给定任一初始状态x(0)≠0，则系统将在有限时间内收敛到原点；求解方程(2)可得系统从状态x(0)到原点经历的时间下面考虑Terminal有限机理的解释，考虑在平衡点x＝0处的Jacobian矩阵为把J看做一阶矩阵的特征值λ，那么当x→0+时，J→‑∞，那么系统的轨迹在负无穷大的特征值的驱动下，自然会以无穷大的速度收敛至平衡点；因此系统将在有限时间内收敛到原点；对于非线性系统，采用基于微分几何理论的输入输出线性化方法，可以在保留系统非线性特点的基础上降低系统控制器设计的复杂度，因此我们考虑采用输入输出线性化方法对系统进行处理；下面简要介绍一下该方法的原理；针对仿射非线性系统其中f(x),g(x)为光滑函数；首先通过分别对函数f(x),g(x)求取关于输出函数h(x)系统的李导数，求系统的相对阶，具体形式为：如果则称系统的相对阶为r；如果系统的相对阶r＝n，(n为系统阶数)，则系统为可以完全输入输出线性化的；选择微分同胚变换对(8)式求导，可得变换系统的状态方程为：观察方程(4)，可以发现除了最后关于ζn的方程外，其余n‑1个方程已经是线性形式，并且不含控制量；只有关于ζn的状态方程是非线性的，但对输入u，方程形式上是线性的；重写状态方程(9)其中则状态方程可以变为：其中这时，系统形式上变为线性，并且保留了系统的非线性特性，使得系统更加易于处理；本方法采用NASA兰利实验室公布的高超声速飞行器纵向模型进行研究，其模型如下：模型假设(1)高超声速飞行器为理想的刚体，即不考虑机翼等的弹性自由度；(2)质心位置，转动惯量是质量的函数，质心位置始终在机体轴纵轴变动；(3)飞行器中心和参考力矩中心在机体X轴上；(4)假设飞行器布局是对称的，也即惯性积Ixy，Ixz，Iyz恒为零；(5)忽略操纵面的转动惯量和发动机推力安装角；式中，L为升力，D为阻力，T发动机推力，Myy为滚转力矩，Iyy为飞行器自身转动惯量，r为飞行器与地心的距离，各参数具体的表达式如下：式中ρ为空气密度，S为发动机有效横截面积，各系数的表达式如下：公式中δe代表升降舵偏角，β代表发动机模态，其表达式为：βc为发动机节流阀控制输入；根据高超声速飞行器的实际，将发动机节流阀控制输入βc和升降舵偏角δe作为控制输入，输出选为速度V和高度H；根据(13)‑(16)式可以发现，高超声速飞行器的动力学模型存在严重的非线性与强耦合，并且在方程中不显含输入，无法直接设计控制器，需要对模型进行变换；输入输出线性化方法是非线性控制系统设计与处理的重要方法，通过对系统求导数使系统形式上表现出线性性质，而本身还保留这系统的非线性特性，因此在这里采用输入输出线性化方法将模型进行变换，然后进行控制器设计；选择状态变量x＝[V γ α β h]T，控制输入u＝[βc δe]T，定义系统的输出为y＝[V,h]T，根据输入输出线性化方法，分别对速度V和高度进行3次和4次微分，可以得到：其中为了便于分离出控制量，选择其中微分后，控制输入已经出现在微分方程中，输出动力学方程可以写为：式中Ω2＝[ω21 ω22 ω23 ω24 ω25]Π2＝[π21 π22 π23 π24 π25]经过输入输出线性化变换，高超声速飞行器的非线性模型已经转化为形式上的线性模型，同时模型本身的非线性特性也得到了最大程度的保留；高超声速飞行器的控制目标是控制飞行器输出跟踪一个给定的指令信号，并且保证系统本身的稳定；给定的指令信号为ycom＝[Vd(t),hd(t)]T，则跟踪目标可以表示为：定义系统的跟踪误差为：则控制的目标为保证系统的跟踪误差但高超声速飞行器的飞行速度快，相对的对系统的响应时间要求也高，传统的指数收敛的控制器难以满足飞行器对控制系统有限时间收敛的要求，因此在这里考虑用Terminal滑模控制来设计有限时间收敛的控制器，保证系统的快速收敛；按照Terminal滑模控制的思路，设计跟踪误差的滑模面如下：这里同理可以推的和所设计的滑模面qi和pi(i＝0,1,2)为正奇数，并且pi>qi；β1，β2，α1，α2和α3为正实数；根据滑模控制的条件，当系统到达滑模面时，S＝0，即观察式(22)可以发现，是方程(22)的一个平衡点，并且在滑模面上，系统方程收敛速度为幂函数收敛，所以系统可以在有限时间内收敛到原点；而在的情况下，高超声速飞行器也完成了对指令信号的有效跟踪；所以设计的滑模面符合有限时间收敛的要求；下面将针对设计的滑模面，设计相应的滑模控制律，以保证滑模控制系统...

【技术特征摘要】
1.高超声速飞行器Terminal滑模控制器设计方法，其特征在于：
一阶Terminal滑模中，一阶Terminal滑模定义如下：
其中x为一标量，β＞0,p,q(p＞q)，且p,q为正奇数；无论x为任何实数，xp/q的解必须为
实数；
系统在滑模面上的动态性能为
给定任一初始状态x(0)≠0，则系统将在有限时间内收敛到原点；求解方程(2)
可得系统从状态x(0)到原点经历的时间
下面考虑Terminal有限机理的解释，考虑在平衡点x＝0处的Jacobian矩阵为
把J看做一阶矩阵的特征值λ，那么当x→0+时，J→-∞，那么系统的轨迹在负无穷大的特
征值的驱动下，自然会以无穷大的速度收敛至平衡点；因此系统将在有限时间内收敛到原
点；
对于非线性系统，采用基于微分几何理论的输入输出线性化方法，可以在保留系统非
线性特点的基础上降低系统控制器设计的复杂度，因此我们考虑采用输入输出线性化方法
对系统进行处理；下面简要介绍一下该方法的原理；
针对仿射非线性系统
其中f(x),g(x)为光滑函数；
首先通过分别对函数f(x),g(x)求取关于输出函数h(x)系统的李导数，求系统的相对
阶，具体形式为：
如果则称系统的相对阶为r；
如果系统的相对阶r＝n，(n为系统阶数)，则系统为可以完全输入输出线性化的；选择
微分同胚变换
对(8)式求导，可得变换系统的状态方程为：
观察方程(4)，可以发现除了最后关于ζn的方程外，其余n-1个方程已经是线性形式，并
且不含控制量；只有关于ζn的状态方程是非线性的，但对输入u，方程形式上是线性的；重写
状态方程(9)
其中
则状态方程可以变为：
其中这时，系统形式上变为线性，并且保留了系统的
非线性特性，使得系统更加易于处理；
本方法采用NASA兰利实验室公布的高超声速飞行器纵向模型进行研究，其模型如下：
模型假设
(1)高超声速飞行器为理想的刚体，即不考虑机翼等的弹性自由度；
(2)质心位置，转动惯量是质量的函数，质心位置始终在机体轴纵轴变动；
(3)飞行器中心和参考力矩中心在机体X轴上；
(4)假设飞行器布局是对称的，也即惯性积Ixy，Ixz，Iyz恒为零；
(5)忽略操纵面的转动惯量和发动机推力安装角；
式中，L为升力，D为阻力，T发动机推力，Myy为滚转力矩，Iyy为飞行器自身转动惯量，r
为飞行器与地心的距离，各参数具体的表达式如下：
式中ρ为空气密度，S为发动机有效横截面积，各系数的表达式如下：
公式中δe代表升降舵偏角，β代表发动机模态，其表达式为：
βc为发动机节流阀控制输入；
根据高超声速飞行器的实际，将发动机节流阀控制输入βc和升降舵偏角δe作为控制输
入，输出选为速度V和高度H；
根据(13)-(16)式可以发现，高超声速飞行器的动力学模型存在严重的非线性与强耦
合，并且在方程中不显含输入，无法直接设计控制器，需要对模型进行变换；输入输出线性
化方法是非线性控制系统设计与处理的重要方法，通过对系统求导数使系统形式上表现出
线性性质，...

【专利技术属性】
技术研发人员：姬庆庆，杨祎，陈楠，石莹，李晨宇，
申请(专利权)人：北京工业大学，
类型：发明
国别省市：北京;11