CN104950899A - 一种固定时间收敛的飞行器姿态控制方法 - Google Patents

一种固定时间收敛的飞行器姿态控制方法 Download PDF

Info

Publication number
CN104950899A
CN104950899A CN201510316613.0A CN201510316613A CN104950899A CN 104950899 A CN104950899 A CN 104950899A CN 201510316613 A CN201510316613 A CN 201510316613A CN 104950899 A CN104950899 A CN 104950899A
Authority
CN
China
Prior art keywords
centerdot
omega
cos
sin
formula
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.)
Granted
Application number
CN201510316613.0A
Other languages
English (en)
Other versions
CN104950899B (zh
Inventor
盛永智
任小欢
刘向东
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN201510316613.0A priority Critical patent/CN104950899B/zh
Publication of CN104950899A publication Critical patent/CN104950899A/zh
Application granted granted Critical
Publication of CN104950899B publication Critical patent/CN104950899B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
  • Feedback Control In General (AREA)
  • Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)

Abstract

本发明公开的一种固定时间收敛的飞行器姿态控制方法,涉及一种飞行器姿态控制方法,属于飞行器控制技术领域。本发明包括如下步骤:步骤一、建立再入飞行器动态模型,提出有限时间姿态跟踪任务;步骤二、对步骤一所建立的模型进行反馈线性化处理;步骤三、给出有限时间控制律,实现系统状态从任意初值出发,在设定时间T跟踪上参考轨迹。本发明可针对部分初值信息未知的情况,实现系统误差固定时间收敛,此外,本发明可在系统参数不确定和外部干扰存在的情况下具有良好的鲁棒性。

Description

一种固定时间收敛的飞行器姿态控制方法
技术领域
本发明涉及一种飞行器姿态控制方法,尤其涉及一种固定时间收敛的飞行器姿态控制方法,属于飞行器控制技术领域。
背景技术
对于再入飞行器来讲,再入过程中飞行条件(空域、速域)大范围变化,各通道间耦合严重,呈现出强烈的非线性动态特性。另外,各种不确定性外部扰动的存在以及飞行器气动特性不能精确获知,导致其姿态控制变得异常复杂。再入飞行器控制系统的设计要解决的关键问题是抑制上述非线性、强耦合和不确定性对系统性能的影响。
目前应用较多的非线性控制方法有模糊控制、最优控制、动态逆控制以及滑模变结构控制等。其中滑模控制技术具备很多优点,例如:对参数变化不敏感、能抵抗外界扰动以及快速动态响应等,广泛应用于飞行器姿态控制中。然而,传统的滑模面是线性的,系统渐进收敛,跟踪误差在无穷时间收敛至零,响应特性较差。在实时控制操作中,无限时间收敛特性往往是不够的。
有限时间收敛能够提供更加优越的特性,例如:更快的收敛速率,更高的精度,对不确定性和外部扰动更好的鲁棒性等。为了实现系统动态的有限时间收敛,有学者提出了终端滑模控制方法。该方法能够使得系统动态到达滑模面后误差在有限时间内收敛到0。在该理论的基础上,学者们又提出了快速终端滑模控制方法,使得误差收敛速度进一步得到提升。然而,在终端滑模控制过程中可能会遇到奇异问题。为了克服这个缺陷,学者们提出了非奇异终端滑模控制技术。该方法能够在不添加额外过程的情况下使得奇异问题得到解决。
进一步,有些学者提出了固定时间收敛问题,通过将分段函数引入终端滑模面中实现系统误差在设定的时间收敛。但是该控制器设计时必须精确知道系统的初值信息。但在输出反馈中,控制器设计时需要的初值信息不能完全获知。尤其在飞行器控制中,由于姿态角测量噪声等原因,姿态角导数信息较难通过直接求导获取,因而设定时间收敛问题就失去了意义。
发明内容
本发明的目的是为了解决在无法获知初值信息的情况下实现固定时间收敛的问题,提供一种固定时间收敛的飞行器姿态控制方法。
本发明的目的是通过下述技术方案实现的。
本发明公开的一种固定时间收敛的飞行器姿态控制方法,具体包括如下步骤:
步骤一、建立再入飞行器动态模型,提出有限时间姿态跟踪任务。
基于无动力再入飞行器的姿态控制问题,姿态动力学方程如下:
ω · x = - I zz ( I zz - I yy ) + I xz 2 I * ω y ω z + I xz ( I zz + I xx - I yy ) I * ω x ω y + I zz I * M x + I xz I * M z ,
ω · y = - I xx - I zz I yy ω x ω z + I xz I yy ( ω z 2 ω x 2 ) + 1 I yy M y , - - - ( 1 )
ω · z = I xx ( I xx - I yy ) + I xz 2 I * ω x ω y + I xz ( I yy + I xx - I zz ) I * ω y ω z + I xz I * M x + I xx I * M z ,
其中,ωxy和ωz分别为滚转角速度、偏航角速度和俯仰角速度。Mx,My,Mz分别为滚转、偏航、俯仰转矩。Iij(i=x,y,z;j=x,y,z)是转动惯量和惯量积。对于几何外形相对于xz平面对称,且质量分布也对称的飞行器。Ixy=Iyz=0, I * = I xx I zz - I xz 2 .
运动学方程为:
α . = - tan β ( ω x cos α + ω z sin α ) + ω y + sin μ cos β [ χ . cos γ - φ . sin χ sin γ + ( θ . + Ω E ) × ( cos φ cos χ sin γ - sin φ cos γ ) ] - cos μ cos β [ γ . - φ . cos χ - ( θ . + Ω E ) cos φ sin χ ] ,
β . = - ω z cos α + ω x sin α + sin μ [ γ . - φ . cos χ + ( θ . + Ω E ) cos φ sin χ ] + cos μ [ χ . cos γ - φ . sin χ sin γ - ( θ . + Ω E ) ( cos φ cos χ sin γ - sin φ cos γ ) ] , - - - ( 2 )
μ . = - ω x cos α cos β - ω y sin β - ω z sin α cos β + α . sin β - χ . sin γ - φ . sin χ cos γ + ( θ . + Ω E ) ( cos φ cos χ cos γ + sin φ sin γ ) ,
其中,α,β,μ分别为攻角、侧滑角和倾侧角。χ,γ分别为航向角和航迹角,φ,θ分别为纬度和经度,ΩE为地球自转角速度。
由舵面产生的控制力矩为:
M x M y M z = 1 2 ρV 2 Sb C Mx ( α , Mz , δ ) C My ( α , Ma , δ ) C Mz ( α , Ma , δ ) - - - ( 3 )
其中,ρ是大气密度,Ma是马赫数,V为相对地面的飞行速度,S,b分别为飞行器的参考面积和参考长度。CMx,CMy,CMz,分别是与α,Ma和舵面相关的力矩系数。δear分别为升降舵,滚转舵和偏航舵。
再入姿态控制的目的是设计控制力矩u,并根据上式的表达式映射成舵面偏角指令δ,使得姿态角在参数不确定性和外部干扰存在的情况下,在有限时间T跟踪上制导指令的输出。即:
lim t > T ( y - y c ) = 0
其中y=[α,β,μ]T,yc=[αccc]T
步骤二、对步骤一所建立的模型进行反馈线性化处理。
将步骤一所得系统模型公式(1)、(2)改写成MIMO仿射非线性形式:
x . = f ( x ) + g ( x ) u y = h ( x )
应用反馈线性化理论,对输出变量进行求导,直到输出方程中显含控制量u。并引入辅助控制量v。将系统解耦成如下的不确定二阶系统
y . . = v + Δv - - - ( 4 )
Δv代表聚合扰动,假设该扰动有界。
步骤三、给出有限时间控制律,实现系统状态从任意初值出发,在设定时间T跟踪上参考轨迹。
步骤3.1,给出高阶滑模观测器。
解耦后不确定二阶系统,即公式(4)改写成如下形式:
z . i 1 = z i 2 z . i 2 = v i + Δυ i
其中i=1,2,3,zi1=yi是姿态角。
设计高阶滑模观测器如下:
z ^ . i 1 = u i
u i = u i 1 - λ i | z ^ i 1 - z i 1 | 1 / 2 sign ( z ^ i 1 - z i 1 ) - - - ( 5 )
u . i 1 = - κ i sign ( z ^ i 1 - z i 1 )
其中λii∈R+,ui是观测器输出。
通过上述给出的观测器,即公式(5),得到姿态角导数。
步骤3.2,给出滑模面。
s ( t ) = e . + ke , t < t 1 - - - ( 6 )
s ( t ) = e . + ke + f ( t ) , t > t 1 - - - ( 7 )
其中,k是滑模面参数,姿态角导数信息由公式(5)得到。t1大于滑模观测器的收敛时间。f(t)是如下定义的分段函数:
f ( t ) = f 1 ( t ) t 1 < t &le; T 0 t > T
T是设定的误差收敛时间,函数f(t)及其参数选择满足如下条件:
s ( t 1 ) = e . ( t 1 ) + ke ( t 1 ) + f ( t 1 ) = 0 - - - ( 8 )
v(t1-)=v(t1+)    (9)
lim t &RightArrow; T f 1 ( t ) = 0 - - - ( 10 )
lim t &RightArrow; T f . 1 ( t ) = 0 - - - ( 11 )
式(8)表明系统状态在t1时刻处于滑模面上,式(9)表明系统控制输出在t1时刻是连续的,式(10)(11)表明T是期望的收敛时间,同时滑模面在该时刻也是连续的。f(t)的具体表达式如下:
f ( t ) = A 4 ( t - t 1 ) 4 + A 3 ( t - t 1 ) 3 + A 2 ( t - t 1 ) 2 t 1 < t &le; T + A 1 ( t - t 1 ) + A 0 0 t > T - - - ( 12 )
Ai=[Ai1,Ai2,Ai3],(i=0,...4)是分段函数系数。
步骤3.3,求解得到滑模控制量。
根据步骤3.2所得t<t1时的线性滑模面,即公式(6),求解得到滑模控制量
v 1 = &Omega; . . c - k e . - &eta;sgin ( s ) , t < t 1 - - - ( 13 )
根据步骤3.2所得t>t1时的线性滑模面,即公式(7),求解得到滑模控制量
v 2 = &Omega; . . c - k e . - f . ( t ) - &eta;sgin ( s ) , t > t 1 - - - ( 14 )
η=diag(η123)是切换增益。满足η1≥||Δvi||1i,ε1i为任意正数。||Δvi||是扰动上界。
sign(si)定义如下:
sign ( s i ) = - 1 s i < 0 0 s i = 0 1 s i > 0
为了减小控制量抖振,采用了如下饱和函数代替切换函数sgn(S):
sat ( S ) = &epsiv; - 1 S , | S | &le; &epsiv; sgn ( S ) , otherwise
通过步骤3得到的控制量v1和v2(其中在0≤t≤t1时间内选择v1,在t>t1时间内选择v2)即可使得在初值信息未知的情况下跟踪误差在期望的时间T收敛到0。即y-yc=0,t≥T
所述的公式(12)中分段函数系数选择方法为:
分段函数中的系数由下述方程组给出:
根据公式(8)得:
s ( t 1 ) = e . ( t 1 ) + ke ( t 1 ) + f ( t 1 ) = A 0 = 0 - - - ( 15 )
根据公式(9)得
f . ( t 1 + ) = A 1 = k 1 s ( t 1 - ) + &eta;sgn ( s ( t 1 - ) ) - - - ( 16 )
根据公式(10)得
f(T)=A4(T-t1)4+A3(T-t1)3+A2(T-t1)2+A1(T-t1)+A0=0  (17)
根据公式(11)得
f . ( T ) = 4 A 4 ( T - t 1 ) 3 + 3 A 3 ( T - t 1 ) 2 + 2 A 2 ( T - t 1 ) + A 1 = 0 - - - ( 18 )
系统在t>t1时的跟踪误差曲线为:
e ( t ) = Ce - k ( t - t 1 ) + B 4 ( t - t 1 ) 4 + B 3 ( t - t 1 ) 3 + B 2 ( t - t 1 ) 2 t 1 < t &le; T + B 1 ( t - t 1 ) + B 0 e - k ( t - T ) e ( T ) t > T - - - ( 19 )
C=[C1,C2,C3]T,Bi=[Bi1,Bi2,Bi3],(i=0,...4)是待定系数。将公式(19)代入滑模面(7)中,参数C=[C1,C2,C3]T,Bi=[Bi1,Bi2,Bi3],(i=0,...4)可以由参数Ai=[Ai1,Ai2,Ai3],(i=0,...4)表示,由于期望的收敛时间是T,误差方程满足:
e ( T ) = Ce - k ( T - t 1 ) + B 4 ( T - t 1 ) 4 + B 3 ( T - t 1 ) 3 + B 2 ( T - t 1 ) 2 + B 1 ( T - t 1 ) + B 0 = 0 - - - ( 20 )
解方程组(15),(16),(17),(18),(20)即可求得待定系数Ai=[Ai1,Ai2,Ai3],(i=0,...4)。
由于e(T)=0,从式(20)看出e(t)=0,(t>T),系统响应具有下述形式:
e ( t ) = Ce - k ( t - t 1 ) + B 4 ( t - t 1 ) 4 + B 3 ( t - t 1 ) 3 + B 2 ( t - t 1 ) 2 t 1 < t &le; T + B 1 ( t - t 1 ) + B 0 0 t > T
可以看出采用本发明给出的滑模面和控制器可以在初值信息未知的情况下实现设定时间收敛。
有益效果
1、本发明的一种固定时间收敛的飞行器姿态控制方法,用于飞行器控制中,由于存在姿态角测量噪声等原因,姿态角导数信息较难通过直接求导获取,设计高阶滑模观测器实现对姿态角导数信息的估计
2、本发明的一种固定时间收敛的飞行器姿态控制方法,针对部分初值信息未知的情况,通过固定时间收敛的控制方法,实现系统误差设定时间收敛。
3、本发明的一种固定时间收敛的飞行器姿态控制方法,系统在参数不确定和外部干扰存在的情况下,采用本发明设计的控制方法,可以保证良好的鲁棒性。
附图说明
图1为本发明方法的流程图;
图2为具体实施方式中系统姿态角跟踪曲线图;(i)攻角跟踪曲线;(ii)侧滑角跟踪曲线;(iii)倾侧角跟踪曲线;
图3为具体实施方式中舵面偏转曲线图;(i)CSMC作用下的舵面偏转曲线;(ii)HONTSM作用下的舵面偏转曲线;
图4为具体实施方式中滑模面响应曲线图;(i)CSMC作用下的滑模面响应曲线;(ii)HONTSM作用下的滑模面响应曲线。
具体实施方式
为了更好的说明本发明的目的和优点,下面结合附图和实例对技术方案做进一步详细说明。
步骤一、建立再入飞行器动态模型,提出有限时间姿态跟踪任务。
基于无动力再入飞行器的姿态控制问题,姿态动力学方程如下:
&omega; &CenterDot; x = - I zz ( I zz - I yy ) + I xz 2 I * &omega; y &omega; z + I xz ( I zz + I xx - I yy ) I * &omega; x &omega; y + I zz I * M x + I xz I * M z ,
&omega; &CenterDot; y = - I xx - I zz I yy &omega; x &omega; z + I xz I yy ( &omega; z 2 &omega; x 2 ) + 1 I yy M y , - - - ( 21 )
&omega; &CenterDot; z = I xx ( I xx - I yy ) + I xz 2 I * &omega; x &omega; y + I xz ( I yy + I xx - I zz ) I * &omega; y &omega; z + I xz I * M x + I xx I * M z ,
其中,ωxy和ωz分别为滚转角速度、偏航角速度和俯仰角速度。Mx,My,Mz分别为滚转、偏航、俯仰转矩。Iij(i=x,y,z;j=x,y,z)是转动惯量和惯量积。对于几何外形相对于xz平面对称,且质量分布也对称的飞行器。Ixy=Iyz=0, I * = I xx I zz - I xx 2 .
运动学方程为:
&alpha; . = - tan &beta; ( &omega; x cos &alpha; + &omega; z sin &alpha; ) + &omega; y + sin &mu; cos &beta; [ &chi; . cos &gamma; - &phi; . sin &chi; sin &gamma; + ( &theta; . + &Omega; E ) &times; ( cos &phi; cos &chi; sin &gamma; - sin &phi; cos &gamma; ) ] - cos &mu; cos &beta; [ &gamma; . - &phi; . cos &chi; - ( &theta; . + &Omega; E ) cos &phi; sin &chi; ] ,
&beta; . = - &omega; z cos &alpha; + &omega; x sin &alpha; + sin &mu; [ &gamma; . - &phi; . cos &chi; + ( &theta; . + &Omega; E ) cos &phi; sin &chi; ] + cos &mu; [ &chi; . cos &gamma; - &phi; . sin &chi; sin &gamma; - ( &theta; . + &Omega; E ) ( cos &phi; cos &chi; sin &gamma; - sin &phi; cos &gamma; ) ] , - - - ( 2 )
&mu; . = - &omega; x cos &alpha; cos &beta; - &omega; y sin &beta; - &omega; z sin &alpha; cos &beta; + &alpha; . sin &beta; - &chi; . sin &gamma; - &phi; . sin &chi; cos &gamma; + ( &theta; . + &Omega; E ) ( cos &phi; cos &chi; cos &gamma; + sin &phi; sin &gamma; ) ,
其中,α,β,μ分别为攻角、侧滑角和倾侧角。χ,γ分别为航向角和航迹角,φ,θ分别为纬度和经度,ΩE为地球自转角速度。
由舵面产生的控制力矩为:
M x M y M z = 1 2 &rho;V 2 Sb C Mx ( &alpha; , Mz , &delta; ) C My ( &alpha; , Ma , &delta; ) C Mz ( &alpha; , Ma , &delta; ) - - - ( 23 )
其中,ρ是大气密度,Ma是马赫数,V为相对地面的飞行速度,S,b分别为飞行器的参考面积和参考长度。CMx,CMy,CMz,分别是与α,Ma和舵面相关的力矩系数。δear分别为升降舵,滚转舵和偏航舵。
再入姿态控制的目的是设计控制力矩u,并根据上式的表达式映射成舵面偏角指令δ,使得姿态角在参数不确定性和外部干扰存在的情况下,在有限时间T跟踪上制导指令的输出。即:
lim t > T ( y - y c ) = 0
其中y=[α,β,μ]T,yc=[αccc]T
考虑初始高度为30km,初始速度2800m/s,初始角速度为ωx(0)=0deg,ωy(0)=1deg,ωz(0)=0.9deg,初始姿态角α0=0deg,β0=-1deg,μ0=0deg。姿态角的给定指令为αc=3deg,βc=0deg,μc=3deg的跟踪状况。
步骤二、对步骤一所建立的模型进行反馈线性化处理。
在如下情况下:
(1)不考虑地球自转的影响,ΩE=0
(2)飞行器旋转运动大于平移运动,因而忽略平移运动产生的角速度。 &chi; . = &phi; . = &theta; . = &gamma; . = 0 .
将步骤一所得系统模型公式(21)(22)改写成如下MIMO仿射非线性形式:
x . = f ( x ) + g ( x ) u y = h ( x )
其中状态向量为输出向量y=[α,β,μ]T,控制向量u=[Mx,My,Mz]T,g(x)=[g1(x),g2(x),g3(x)]h(x)=[h1(x),h2(x),h3(x)]T具体表达式可由公式(21)(22)整理得到。
应用反馈线性化理论,对输出变量进行求导,直到输出方程中显含控制量u。得到下式:
y . . = F ( x ) + E ( x ) u ,
F ( x ) = L f 2 h 1 ( x ) L f 2 h 2 ( x ) L f 2 h 3 ( x ) ,
E ( x ) = L g 2 L f h 1 ( x ) L g 2 L f h 1 ( x ) L g 3 L f h 1 ( x ) L g 1 L f h 2 ( x ) L g 2 L f h 2 ( x ) L g 3 L f h 3 ( x ) L g 1 L f h 3 ( x ) L g 2 L f h 2 ( x ) L g 3 L f h 3 ( x )
det E ( x ) = 1 I * I yy cos &beta; .
由于采用倾斜转弯控制方式,cosβ≈1,因而detE(x)≠0,E(x)可逆,设计控制律:
u=E(x)-1[v-F(x)],    (24)
可以实现输入输出反馈线性化,v=[v1,v2,v3]T是引入的辅助控制量。
由于系统的相对阶是6,等于系统方程的维数。可以完全线性化,且不存在内动态。模型参数不确定性和外部扰动的存在,反馈线性化不精确,通过控制律(24)可将系统解耦成如下的不确定二阶系统
y . . = v + &Delta;v - - - ( 25 )
Δv代表聚合扰动,假设该扰动有界。验证时扰动设置为:
大气密度拉偏20%,力系数和力矩系数拉偏20%,质量和转动惯量拉偏10%。并施与如下形式的外部干扰:
&Delta;M = 10000 sin ( t ) 50000 sin ( 0.8 t + &pi; 4 ) 50000 sin ( 0.5 t + 3 &pi; 4 ) N &CenterDot; m
步骤三、给出有限时间控制律,实现系统状态从任意初值出发,在设定时间T跟踪上参考轨迹。
步骤3.1,高阶滑模观测器设计。
解耦后不确定二阶系统(4)改写成如下形式:
z . i 1 = z i 2 z . i 2 = v i + &Delta;&upsi; i
其中i=1,2,3,zi1=yi是姿态角。
设计高阶滑模观测器如下:
z ^ . i 1 = u i
u i = u i 1 - &lambda; i | z ^ i 1 - z i 1 | 1 / 2 sign ( z ^ i 1 - z i 1 ) - - - ( 26 )
u . i 1 = - &kappa; i sign ( z ^ i 1 - z i 1 )
其中λii∈R+,ui是观测器输出。取κi=8,λi=10。
如果状态变量z1有界并且Lebesgue可测,则通过合适的参数选择,状态观测值在有限时间收敛到其真实值。其中参数选择满足:
κi>Ci
&lambda; i 2 &GreaterEqual; 4 C i &kappa; i + C i &kappa; i - C i
其中Ci是Lebesgue常数,满足取Ci=1.5,收敛时间由下式给出:
T &le; 2 &kappa; i - C i 1 1 - &Phi; ( &kappa; i , &lambda; i , C i )
其中Φ(κii,Ci)=|Ψ(t*)|具体定义如下:
&Sigma; . i = - | &Sigma; | i 1 2 + &Psi; i
&Psi; . i = - 1 &lambda; i 2 ( &kappa; i - C i ) - | &Sigma; | i 1 2 + &Psi; i > 0 - 1 &lambda; i 2 ( &kappa; i + C i ) - | &Sigma; | i 1 2 + &Psi; i &le; 0
Σi(0)=0 Ψi(0)=1
其中t*=inf{t|t>0,Σi(t)=0,Ψi(t)<0}。计算得到Φ(κii,C)=1.167×10-6。相应的收敛时间为:
t converge &le; 1 &kappa; i - C i 1 1 - &Phi; ( &kappa; i , &lambda; i , C i ) < 0.3077
通过上述给出的观测器(26),可以实现对姿态角导数的估计。
步骤3.2,设计滑模面。
s ( t ) = e . + ke , t < t 1 - - - ( 27 )
s ( t ) = e . + ke + f ( t ) , t > t 1 - - - ( 28 )
其中,k是滑模面参数,取为k=diag(1,1,1),姿态角导数信息由高阶滑模观测器(26)估计得到。t1大于滑模观测器的收敛时间。t1=0.5s>tconverge,f(t)是如下定义的分段函数:
f ( t ) = f 1 ( t ) t 1 < t &le; T 0 t > T
T是设定的误差收敛时间,函数f(t)及其参数选择满足如下条件:
s ( t 1 ) = e . ( t 1 ) + ke ( t 1 ) + f ( t 1 ) = 0 - - - ( 29 )
v(t1-)=v(t1+)    (30)
lim t &RightArrow; T f 1 ( t ) = 0 - - - ( 31 )
lim t &RightArrow; T f . 1 ( t ) = 0 - - - ( 32 )
式(29)表明系统状态在t1时刻处于滑模面上,式(30)表明系统控制输出在t1时刻是连续的,式(31)(32)表明T是期望的收敛时间,取为T=3s,同时滑模面在该时刻也是连续的。f(t)的具体表达式如下:
f ( t ) = A 4 ( t - t 1 ) 4 + A 3 ( t - t 1 ) 3 + A 2 ( t - t 1 ) 2 t 1 < t &le; T + A 1 ( t - t 1 ) + A 0 0 t > T - - - ( 33 )
Ai=[Ai1,Ai2,Ai3],(i=0,...4)是分段函数系数。
步骤3.3,求解得到滑模控制量
根据步骤3.2所得t<t1时的线性滑模面,即公式(27),求解得到滑模控制量
v 1 = &Omega; . . c - k e . - &eta;sgin ( s ) , t < t 1 - - - ( 34 )
根据步骤3.2所得t>t1时的线性滑模面,即公式(28),求解得到滑模控制量
v 2 = &Omega; . . c - k e . - f . ( t ) - &eta;sgin ( s ) , t > t 1 - - - ( 35 )
η=diag(η123)是切换增益。满足η1≥||Δvi||1i,ε1i为任意正数。||Δvi||是扰动上界。
sign(si)定义如下:
sign ( s i ) = - 1 s i < 0 0 s i = 0 1 s i > 0
为了减小控制量抖振,采用了如下饱和函数代替切换函数sgn(S):
sat ( S ) = &epsiv; - 1 S , | S | &le; &epsiv; sgn ( S ) , otherwise
针对式(25)所表示的不确定二阶系统,在t<t1时选择式(34)所示的滑模控制律,闭环系统是渐进稳定的。
定义如下正定的Lyapunov函数:
V = 1 2 s T s
对上式进行微分得到
V . = s T s . = s T ( v + &Delta;v - &Omega; . . c + k e . ) = S T ( &Delta;v - &eta;sgn ( S ) ) &le; &Sigma; i = 1 3 | s i | | | &Delta;v i | | - &eta; i | s i | &le; - &epsiv; min ( | s 1 | + | s 2 | + | s 3 | ) &le; - &epsiv; min ( s T s ) 1 2 = - &epsiv; min 2 ( v ) 1 2
其中εmin=min(ε123),根据Lyapunov有限时间稳定原理可知滑模面是渐进收敛的,系统闭环稳定。
针对式(25)所表示的不确定二阶系统,选择式(35)所示的时变滑模控制律,系统状态在t>t1之后处于滑模段运动,即对于t∈[t1,+∞),有s≡0。
如下的Lyapunov函数
V = 1 2 s T s
沿控制律(35)作用下的闭环轨迹求导:
V . = s T s . = s T ( e . . + k e . + f . ( t ) ) = s T ( v + &Delta;v - &Omega; . . c + k e . + f . ( t ) ) = s T ( &Omega; . . c - k e . - f . ( t ) - &eta;sgin ( s ) + &Delta;v - &Omega; . . c + k e . + f . ( t ) ) = s T ( &Delta;v - &eta;sgn ( s ) ) &le; | s i | | | &Delta;v | | - &eta; | s i | &le; 0
由于选择的实变函数f(t)使得s(t1)=0,也即由于V是不增的,即V(t)≤V(t1)=0,又由V的表达式可知V≥0,所以由上述过程可以得到对于t∈[t1,+∞),有V(t)≡0,也即对于t∈[t1,+∞),s(t)≡0。因此说明系统状态在t>t1之后处于滑模段运动。
公式(33)中分段函数系数选择方法:
根据t∈[t1,+∞),有s≡0。可以知道系统误差响应在t∈[t1,+∞)时完全由滑模面决定。即在t≥t1时:
s ( t ) = e . + ke + f ( t ) = 0 - - - ( 36 )
求解上述微分方程可以得到姿态跟踪误差响应如下:
e ( t ) = Ce - k ( t - t 1 ) + B 4 ( t - t 1 ) 4 + B 3 ( t - t 1 ) 3 + B 2 ( t - t 1 ) 2 t 1 < t &le; T + B 1 ( t - t 1 ) + B 0 e - k ( t - T ) e ( T ) t > T - - - ( 37 )
C=[C1,C2,C3]T,Bi=[Bi1,Bi2,Bi3],(i=0,...4)是待定系数。将公式(37)带入到滑模面方程(36),参数Bi=[Bi1,Bi2,Bi3],(i=0,...4)可以由参数Ai=[Ai1,Ai2,Ai3],(i=0,...4)表示:
B 4 = - A 4 k
B 3 = 4 k 2 A 4 - A 3 k
B 2 = - 12 k 3 A 4 + 3 k 2 A 3 - A 2 k - - - ( 38 )
B 1 = - 1 k A 1 + 24 k 4 A 4 - 6 k 3 A 3 + 2 k 2 A 2
B 0 = - A 0 k + 1 k 2 A 1 - 24 k 5 A 4 + 6 k 4 A 3 - 2 k 3 A 2
由于e(t1)=C+B0,C=[C1,C2,C3]T可以表示为
C=-B0+e(t1)    (39)
从而Ci,Bi均可以由Ai表示出。
由于期望的收敛时间是T,误差方程满足:
e ( T ) = Ce - k ( T - t 1 ) + B 4 ( T - t 1 ) 4 + B 3 ( T - t 1 ) 3 + B 2 ( T - t 1 ) 2 + B 1 ( T - t 1 ) + B 0 = 0 - - - ( 40 )
根据公式(29)得:
s ( t 1 ) = e . ( t 1 ) + ke ( t 1 ) + f ( t 1 ) = A 0 = 0 - - - ( 41 )
根据公式(30)
v ( t 1 - ) = &Omega; . . c ( t 1 ) - k e . ( t 1 ) - k 1 s ( t 1 - ) - &eta;sgn ( s ( t 1 - ) ) = v ( t 1 + ) = &Omega; . . c ( t 1 ) - k e . ( t 1 ) - f . ( t 1 + ) - &eta;sgn ( s ( t 1 + ) )
上式进一步简化为:
f . ( t 1 + ) = A 1 = k 1 s ( t 1 - ) + &eta;sgn ( s ( t 1 - ) ) - - - ( 42 )
根据公式(31)得
f(T)=A4(T-t1)4+A3(T-t1)3+A2(T-t1)2+A1(T-t1)+A0=0    (43)
根据公式(32)得
f . ( T ) = 4 A 4 ( T - t 1 ) 3 + 3 A 3 ( T - t 1 ) 2 + 2 A 2 ( T - t 1 ) + A 1 = 0 - - - ( 44 )
解方程组(40),(41),(42),(43),(44)即可求得待定系数Ai=[Ai1,Ai2,Ai3],(i=0,...4)。
由于e(T)=0,从式(20)看出e(t)=0,(t>T),系统响应具有下述形式:
e ( t ) = Ce - k ( t - t 1 ) + B 4 ( t - t 1 ) 4 + B 3 ( t - t 1 ) 3 + B 2 ( t - t 1 ) 2 t 1 < t &le; T + B 1 ( t - t 1 ) + B 0 0 t > T
可以看出采用本发明给出的滑模面和控制器可以在初值信息未知的情况下实现设定时间收敛。
为了对比本发明给出的基于高阶滑模观测器的时变非奇异终端滑模控制律(High-order sliding mode observer based time-varying nonsingular terminal slidingmode attitude controller-HONTSM)的快速收敛特性,对比传统滑模控制器(conventional sliding mode control-CSMC):
s ( t ) = e . + ke
v = &Omega; . . c - k e . - &eta;sat ( s )
滑模面参数选择为k=diag(1,1,1)
如图2所示,是姿态角跟踪曲线,由曲线可以看出,在系统不确定性和外部扰动存在的情况下,姿态角在设定的时间T=3s跟踪上了给定指令,同时表明系统具有良好的鲁棒性。同样对比本发明与传统的滑模控制器作用下的响应曲线,可以看出本发明提出的控制器具有较快的响应速度。如图3所示,是舵面偏角指令,可以看出在本发明提出的滑模控制器作用下的控制量输出平滑。如图4所示,给出了两种控制器作用下的滑模面响应曲线,从图3(ii)可以看出自t=t1=0.5时刻之后,滑模面一直处于边界层内。也说明了系统响应在t>t1之后完全由滑模面决定。
综上所述,该发明提出的控制律鲁棒性强,能够使得误差在初值状态未知的情况下实现固定时间收敛,并且可以调节收敛速度,具有很高的工程应用价值。
以上所述的具体描述,对发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例,用于解释本发明,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (2)

1.一种固定时间收敛的飞行器姿态控制方法,包括如下步骤,步骤一、建立再入飞行器动态模型,提出有限时间姿态跟踪任务;
基于无动力再入飞行器的姿态控制问题,姿态动力学方程如下:
&omega; &CenterDot; x = - I zz ( I zz - I yy ) + I xz 2 I * &omega; y &omega; z + I xz ( I zz + I xx - I yy ) I * &omega; x &omega; y + I zz I * M x + I xz I * M z ,
&omega; &CenterDot; y = - I xx - I zz I yy &omega; x &omega; z + I xz I yy ( &omega; z 2 - &omega; x 2 ) + 1 I yy M y , - - - ( 1 )
&omega; &CenterDot; z = I xx ( I xx - I yy ) + I xz 2 I * &omega; x &omega; y + I xz ( I yy - I xx - I zz ) I * &omega; y &omega; z + I xz I * M x + I xx I * M z ,
其中,ωxy和ωz分别为滚转角速度、偏航角速度和俯仰角速度;Mx,My,Mz分别为滚转、偏航、俯仰转矩;Iij(i=x,y,z;j=x,y,z)是转动惯量和惯量积;对于几何外形相对于xz平面对称,且质量分布也对称的飞行器;Ixy=Iyz=0, I * = I xx I zz - I xz 2 ;
运动学方程为:
&alpha; &CenterDot; = - tan &beta; ( &omega; x cos &alpha; + &omega; z sin &alpha; ) + &omega; y + sin &mu; cos &beta; [ &chi; &CenterDot; cos &gamma; - &phi; &CenterDot; sin &chi; sin &gamma; + ( &theta; &CenterDot; + &Omega; E ) &times; ( cos &phi; cos &chi; sin &gamma; - sin &phi; cos &gamma; ) ] - cos &mu; cos &beta; [ &gamma; &CenterDot; - &phi; &CenterDot; cos &chi; - ( &theta; &CenterDot; + &Omega; E ) cos &phi; sin &chi; ] ,
&beta; &CenterDot; = - &omega; z cos &alpha; + &omega; x sin &alpha; + sin &mu; [ &gamma; &CenterDot; - &phi; &CenterDot; cos &chi; + ( &theta; &CenterDot; + &Omega; E ) cos &phi; sin &chi; ] + cos &mu; [ &chi; &CenterDot; cos &gamma; - &phi; &CenterDot; sin &chi; sin &gamma; - ( &theta; &CenterDot; + &Omega; E ) ( cos &phi; cos &chi; sin &gamma; - sin &phi; cos &gamma; ) ] , - - - ( 2 )
&mu; &CenterDot; = - &omega; x cos &alpha; cos &beta; - &omega; y sin &beta; - &omega; z sin &alpha; cos &beta; + &alpha; &CenterDot; sin &beta; - &chi; &CenterDot; sin &gamma; - &phi; &CenterDot; sin &chi; cos &gamma; + ( &theta; &CenterDot; + &Omega; E ) ( cos &phi; cos &chi; cos &gamma; + sin &phi; sin &gamma; ) ,
其中,α,β,μ分别为攻角、侧滑角和倾侧角;χ,γ分别为航向角和航迹角,φ,θ分别为纬度和经度,ΩE为地球自转角速度;
由舵面产生的控制力矩为:
M x M y M z = 1 2 &rho;V 2 Sb C Mx ( &alpha; , Ma , &delta; ) C My ( &alpha; , Ma , &delta; ) C Mz ( &alpha; , Ma , &delta; ) - - - ( 3 )
其中,ρ是大气密度,Ma是马赫数,V为相对地面的飞行速度,S,b分别为飞行器的参考面积和参考长度;CMx,CMy,CMz,分别是与α,Ma和舵面相关的力矩系数;δear分别为升降舵,滚转舵和偏航舵;
再入姿态控制的目的是设计控制力矩u,并根据上式的表达式映射成舵面偏角指令δ,使得姿态角在参数不确定性和外部干扰存在的情况下,在有限时间T跟踪上制导指令的输出;即:
lim t > t f ( y - y c ) = 0
其中y=[α,β,μ]T,yc=[αccc]T
步骤二、对步骤一所建立的模型进行反馈线性化处理;
将步骤一所得系统模型公式(1)、(2)改写成MIMO仿射非线性形式:
x &CenterDot; = f ( x ) + g ( x ) u y = h ( x )
应用反馈线性化理论,对输出变量进行求导,直到输出方程中显含控制量u;并引入辅助控制量v;将系统解耦成如下的不确定二阶系统
y &CenterDot; &CenterDot; = v + &Delta;v - - - ( 4 )
Δv代表聚合扰动,假设该扰动有界;
其特征在于:还包括如下步骤:
步骤三、给出有限时间控制律,实现系统状态从任意初值出发,在设定时间T跟踪上参考轨迹;
步骤3.1,给出高阶滑模观测器
解耦后不确定二阶系统,即公式(4)改写成如下形式:
z &CenterDot; i 1 = z i 2
z &CenterDot; i 2 = v i + &Delta; &upsi; i
其中i=1,2,3,zi1=yi是姿态角;
给出高阶滑模观测器如下:
z ^ &CenterDot; i 1 = u i
u i = u i 1 - &lambda; i | z ^ i 1 - z i 1 | 1 / 2 sign ( z ^ i 1 - z i 1 ) - - - ( 5 )
u &CenterDot; i 1 = - &kappa; i sign ( z ^ i 1 - z i 1 )
其中λii∈R+,ui是观测器输出;
通过上述设计的观测器,即公式(5),得到姿态角导数;
步骤3.2,给出滑模面;
s ( t ) = e &CenterDot; + ke , t < t 1 - - - ( 6 )
s ( t ) = e &CenterDot; + ke + f ( t ) , t > t 1 - - - ( 7 )
其中,k是滑模面参数,姿态角导数信息由公式(5)得到;t1大于滑模观测器的收敛时间;f(t)是如下定义的分段函数:
f ( t ) = f 1 ( t ) t 1 < t &le; T 0 t > T
T是设定的误差收敛时间,函数f(t)及其参数选择满足如下条件:
s ( t 1 ) = e &CenterDot; ( t 1 ) + ke ( t 1 ) + f ( t 1 ) = 0 - - - ( 8 )
v(t1-)=v(t1+)               (9)
lim t &RightArrow; T f 1 ( t ) = 0 - - - ( 10 )
lim t &RightArrow; T f &CenterDot; 1 ( t ) = 0 - - - ( 11 )
式(8)表明系统状态在t1时刻处于滑模面上,式(9)表明系统控制输出在t1时刻是连续的,式(10)(11)表明T是期望的收敛时间,同时滑模面在该时刻也是连续的;f(t)的具体表达式如下:
f ( t ) = A 4 ( t - t 1 ) 4 + A 3 ( t - t 1 ) 3 + A 2 ( t - t 1 ) 2 t 1 < t &le; T + A 1 ( t - t 1 ) + A 0 0 t > T - - - ( 12 )
Ai=[Ai1,Ai2,Ai3],(i=0,...4)是分段函数系数;
步骤3.3,求解得到滑模控制量
根据步骤3.2所得t<t1时的线性滑模面,即公式(6),求解得到滑模控制量
v 1 = &Omega; &CenterDot; &CenterDot; c - k e &CenterDot; - &eta;sgin ( s ) , t < t 1 - - - ( 13 )
根据步骤3.2所得t>t1时的线性滑模面,即公式(7),求解得到滑模控制量
v 2 = &Omega; &CenterDot; &CenterDot; c - k e &CenterDot; - f &CenterDot; ( t ) - &eta;sgin ( s ) , t > t 1 - - - ( 14 )
η=diag(η123)是切换增益;满足η1≥||Δvi||1i,ε1i为任意正数;||Δvi||是扰动上界;
sign(si)定义如下:
sign ( s i ) = - 1 s i < 0 0 s i = 0 1 s i > 0
为了减小控制量抖振,采用如下饱和函数代替切换函数sgn(S):
sat ( s i ) = &epsiv; - 1 S , | S | &le; &epsiv; sgn ( S ) , otherwise
通过步骤三得到的控制量v1和v2(其中在0≤t≤t1时间内选择v1,在t>t1时间内选择v2)即可使得在初值信息未知的情况下跟踪误差在期望的时间T收敛到0,即y-yc=0,t≥T。
2.如权利要求1所述的一种固定时间收敛的飞行器姿态控制方法,其特征在于:所述公式(12)中分段函数系数选择方法为,
分段函数中的系数由下述方程组给出
根据公式(8)得:
s ( t 1 ) = e &CenterDot; ( t 1 ) + ke ( t 1 ) + f ( t 1 ) = A 0 = 0 - - - ( 15 )
根据公式(9)得:
f &CenterDot; ( t 1 + ) = A 1 = k 1 s ( t 1 - ) + &eta;sgn ( s ( t 1 - ) ) - - - ( 16 )
根据公式(10)得:
f(T)=A4(T-t1)4+A3(T-t1)3+A2(T-t1)2+A1(T-t1)+A0=0       (17)
根据公式(11)得:
f &CenterDot; = 4 A 4 ( T - t 1 ) 3 + 3 A 3 ( T - t 1 ) 2 + 2 A 2 ( T - t 1 ) + A 1 = 0 - - - ( 18 )
系统在t>t1时的跟踪误差曲线为:
e ( t ) = Ce - k ( t - t 1 ) + B 4 ( t - t 1 ) 4 + B 3 ( t - t 1 ) 3 + B 2 ( t - t 1 ) 2 t 1 < t &le; T + B 1 ( t - t 1 ) + B 0 e - k ( t - T ) e ( T ) t > T - - - ( 19 )
C=[C1,C2,C3]T,Bi=[Bi1,Bi2,Bi3],(i=0,...4)是待定系数;将公式(19)代入滑模面(7)中,参数C=[C1,C2,C3]T,Bi=[Bi1,Bi2,Bi3],(i=0,...4)由参数Ai=[Ai1,Ai2,Ai3],(i=0,...4)表示,由于期望的收敛时间是T,误差方程满足:
e ( T ) = Ce - k ( T - t 1 ) + B 4 ( T - t 1 ) 4 + B 3 ( T - t 1 ) 3 + B 2 ( T - t 1 ) 2 + B 1 ( T - t 1 ) + B 0 = 0 - - - ( 20 )
解方程组(15),(16),(17),(18),(20)即可求得待定系数Ai=[Ai1,Ai2,Ai3],(i=0,...4)。
CN201510316613.0A 2015-06-10 2015-06-10 一种固定时间收敛的飞行器姿态控制方法 Expired - Fee Related CN104950899B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510316613.0A CN104950899B (zh) 2015-06-10 2015-06-10 一种固定时间收敛的飞行器姿态控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510316613.0A CN104950899B (zh) 2015-06-10 2015-06-10 一种固定时间收敛的飞行器姿态控制方法

Publications (2)

Publication Number Publication Date
CN104950899A true CN104950899A (zh) 2015-09-30
CN104950899B CN104950899B (zh) 2017-10-17

Family

ID=54165618

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510316613.0A Expired - Fee Related CN104950899B (zh) 2015-06-10 2015-06-10 一种固定时间收敛的飞行器姿态控制方法

Country Status (1)

Country Link
CN (1) CN104950899B (zh)

Cited By (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105739513A (zh) * 2016-02-05 2016-07-06 北京航空航天大学 一种四旋翼飞行机器人非线性轨迹跟踪控制器及其跟踪控制方法
CN106406102A (zh) * 2016-12-20 2017-02-15 哈尔滨工业大学 一种含干扰观测器的高超声速飞行器跟踪控制方法
CN107450584A (zh) * 2017-08-29 2017-12-08 浙江工业大学 一种基于固定时间滑模的飞行器自适应姿态控制方法
CN107608210A (zh) * 2017-08-25 2018-01-19 湘潭大学 输入饱和的航天器姿态终端滑模跟踪控制方法
CN107703952A (zh) * 2017-08-29 2018-02-16 浙江工业大学 一种刚性飞行器的非奇异固定时间自适应姿态控制方法
CN109213184A (zh) * 2018-11-06 2019-01-15 哈尔滨工业大学(深圳) 挠性航天器的有限时间多模态滑模姿态控制算法
CN109542112A (zh) * 2019-01-08 2019-03-29 哈尔滨工业大学 一种针对垂直起降可重复使用火箭返回飞行的固定时间收敛抗扰控制方法
CN109782795A (zh) * 2018-12-29 2019-05-21 南京航空航天大学 一种利用耦合的面对称高超声速飞行器横侧向控制方法及控制系统
CN109977613A (zh) * 2019-04-19 2019-07-05 哈尔滨工业大学 一种可预先设定调整时间的自适应滑模末制导律设计方法
CN109991848A (zh) * 2019-03-28 2019-07-09 鲁东大学 抗干扰控制方法和装置
CN110471439A (zh) * 2018-09-25 2019-11-19 浙江工业大学 一种基于神经网络估计的刚性飞行器固定时间姿态镇定方法
CN110488854A (zh) * 2018-09-25 2019-11-22 浙江工业大学 一种基于神经网络估计的刚性飞行器固定时间姿态跟踪控制方法
CN111506095A (zh) * 2020-04-24 2020-08-07 北京科技大学 一种双刚体特征点间饱和固定时间相对位姿跟踪控制方法
CN112526992A (zh) * 2020-11-27 2021-03-19 三峡大学 一种混凝土坝缆机入仓轨迹规划与控制方法
CN112650291A (zh) * 2020-12-21 2021-04-13 北京理工大学 基于rcs推力器的飞行器快速姿态控制系统
CN109946973B (zh) * 2019-04-12 2021-05-04 东北大学 一种结合快速终端滑模和线性滑模的组合滑模控制方法
CN112965371A (zh) * 2021-01-29 2021-06-15 哈尔滨工程大学 基于固定时间观测器的水面无人艇轨迹快速跟踪控制方法
CN117250854A (zh) * 2023-11-17 2023-12-19 北京中星时代科技有限公司 控制参数设计引入积分系数的飞行姿态控制方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102880060A (zh) * 2012-10-25 2013-01-16 北京理工大学 再入飞行器自适应指数时变滑模姿态控制方法
CN102929283A (zh) * 2012-11-07 2013-02-13 北京理工大学 基于sdre的再入飞行器自适应最优滑模姿态控制方法
CN103853157A (zh) * 2014-03-19 2014-06-11 湖北蔚蓝国际航空学校有限公司 一种基于自适应滑模的飞行器姿态控制方法
WO2014103030A1 (ja) * 2012-12-28 2014-07-03 日鍛バルブ株式会社 位相可変装置の制御方法及び制御装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102880060A (zh) * 2012-10-25 2013-01-16 北京理工大学 再入飞行器自适应指数时变滑模姿态控制方法
CN102929283A (zh) * 2012-11-07 2013-02-13 北京理工大学 基于sdre的再入飞行器自适应最优滑模姿态控制方法
WO2014103030A1 (ja) * 2012-12-28 2014-07-03 日鍛バルブ株式会社 位相可変装置の制御方法及び制御装置
CN103853157A (zh) * 2014-03-19 2014-06-11 湖北蔚蓝国际航空学校有限公司 一种基于自适应滑模的飞行器姿态控制方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
BAILING TIAN,ET: "Nonlinear robust control for reusable launch vehicles in reentry phase based on time-varying high order sliding mode", 《JOURNAL OF THE FRANKLIN INSTITUTE》 *
BAILING TIAN,ET: "Quasi-continuous high-order sliding mode controller design for reusable launch vehicles in reentry phase", 《AEROSPACE SCIENCE AND TECHNOLOGY》 *
JIANSHUANG SONG,ET: "An Integrated Sliding Model Dynamic Inversion Control Method for Reusable Launch Vehicle Reentry Attitude Control", 《IEEE》 *
王亮,等: "基于高阶滑模观测器的自适应时变滑模再入姿态控制", 《控制与决策》 *

Cited By (28)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105739513B (zh) * 2016-02-05 2018-06-12 北京航空航天大学 一种四旋翼飞行机器人非线性轨迹跟踪控制器及其跟踪控制方法
CN105739513A (zh) * 2016-02-05 2016-07-06 北京航空航天大学 一种四旋翼飞行机器人非线性轨迹跟踪控制器及其跟踪控制方法
CN106406102A (zh) * 2016-12-20 2017-02-15 哈尔滨工业大学 一种含干扰观测器的高超声速飞行器跟踪控制方法
CN106406102B (zh) * 2016-12-20 2019-06-11 哈尔滨工业大学 一种含干扰观测器的高超声速飞行器跟踪控制方法
CN107608210A (zh) * 2017-08-25 2018-01-19 湘潭大学 输入饱和的航天器姿态终端滑模跟踪控制方法
CN107608210B (zh) * 2017-08-25 2020-06-23 湘潭大学 输入饱和的航天器姿态终端滑模跟踪控制方法
CN107450584A (zh) * 2017-08-29 2017-12-08 浙江工业大学 一种基于固定时间滑模的飞行器自适应姿态控制方法
CN107703952A (zh) * 2017-08-29 2018-02-16 浙江工业大学 一种刚性飞行器的非奇异固定时间自适应姿态控制方法
CN107703952B (zh) * 2017-08-29 2020-10-30 浙江工业大学 一种刚性飞行器的非奇异固定时间自适应姿态控制方法
CN110471439A (zh) * 2018-09-25 2019-11-19 浙江工业大学 一种基于神经网络估计的刚性飞行器固定时间姿态镇定方法
CN110471439B (zh) * 2018-09-25 2022-06-17 浙江工业大学 一种基于神经网络估计的刚性飞行器固定时间姿态镇定方法
CN110488854A (zh) * 2018-09-25 2019-11-22 浙江工业大学 一种基于神经网络估计的刚性飞行器固定时间姿态跟踪控制方法
CN109213184B (zh) * 2018-11-06 2021-06-08 哈尔滨工业大学(深圳) 挠性航天器的有限时间多模态滑模姿态控制算法
CN109213184A (zh) * 2018-11-06 2019-01-15 哈尔滨工业大学(深圳) 挠性航天器的有限时间多模态滑模姿态控制算法
CN109782795A (zh) * 2018-12-29 2019-05-21 南京航空航天大学 一种利用耦合的面对称高超声速飞行器横侧向控制方法及控制系统
CN109542112A (zh) * 2019-01-08 2019-03-29 哈尔滨工业大学 一种针对垂直起降可重复使用火箭返回飞行的固定时间收敛抗扰控制方法
CN109991848B (zh) * 2019-03-28 2022-02-01 鲁东大学 抗干扰控制方法和装置
CN109991848A (zh) * 2019-03-28 2019-07-09 鲁东大学 抗干扰控制方法和装置
CN109946973B (zh) * 2019-04-12 2021-05-04 东北大学 一种结合快速终端滑模和线性滑模的组合滑模控制方法
CN109977613A (zh) * 2019-04-19 2019-07-05 哈尔滨工业大学 一种可预先设定调整时间的自适应滑模末制导律设计方法
CN111506095A (zh) * 2020-04-24 2020-08-07 北京科技大学 一种双刚体特征点间饱和固定时间相对位姿跟踪控制方法
CN111506095B (zh) * 2020-04-24 2021-04-06 北京科技大学 一种双刚体特征点间饱和固定时间相对位姿跟踪控制方法
CN112526992A (zh) * 2020-11-27 2021-03-19 三峡大学 一种混凝土坝缆机入仓轨迹规划与控制方法
CN112650291A (zh) * 2020-12-21 2021-04-13 北京理工大学 基于rcs推力器的飞行器快速姿态控制系统
CN112650291B (zh) * 2020-12-21 2022-04-15 北京理工大学 基于rcs推力器的飞行器快速姿态控制系统
CN112965371A (zh) * 2021-01-29 2021-06-15 哈尔滨工程大学 基于固定时间观测器的水面无人艇轨迹快速跟踪控制方法
CN117250854A (zh) * 2023-11-17 2023-12-19 北京中星时代科技有限公司 控制参数设计引入积分系数的飞行姿态控制方法
CN117250854B (zh) * 2023-11-17 2024-02-02 北京中星时代科技有限公司 控制参数设计引入积分系数的飞行姿态控制方法

Also Published As

Publication number Publication date
CN104950899B (zh) 2017-10-17

Similar Documents

Publication Publication Date Title
CN104950899A (zh) 一种固定时间收敛的飞行器姿态控制方法
CN104950898B (zh) 一种再入飞行器全阶非奇异终端滑模姿态控制方法
CN102880060B (zh) 再入飞行器自适应指数时变滑模姿态控制方法
CN105242676B (zh) 一种有限时间收敛时变滑模姿态控制方法
CN105700536B (zh) 基于绳系拖曳系统的主动星姿态和系绳摆振联合控制方法
CN106444799B (zh) 基于模糊扩张状态观测器和自适应滑模的四旋翼无人机控制方法
CN106325291B (zh) 基于滑模控制律和eso的四旋翼飞行器姿态控制方法及系统
CN102929151B (zh) 一种基于指数时变二阶滑模的再入飞行姿态控制方法
CN104238357A (zh) 一种近空间飞行器的容错滑模控制方法
CN107807657B (zh) 一种基于路径规划的挠性航天器姿态自适应控制方法
CN110347170B (zh) 可重复使用运载器再入段鲁棒容错制导控制系统及工作方法
CN104950671A (zh) 基于自适应模糊的再入飞行器pid型滑模姿态控制方法
CN105278545A (zh) 适用于高超声速机动飞行的自抗扰轨迹线性化控制方法
CN105912009A (zh) 基于极点配置和模糊自抗扰控制技术的四旋翼飞行器控制方法
CN104331084B (zh) 一种基于方向舵控滚转策略的气动舵偏范围计算方法
CN107272719B (zh) 基于协调因子的高超声速飞行器姿态运动协调控制方法
CN104460681A (zh) 倾转旋翼无人直升机过渡段的飞行控制方法
CN104317300A (zh) 一种基于模型预测控制的平流层飞艇平面路径跟踪控制方法
CN105138001A (zh) 一种四旋翼飞行器姿态控制方法
CN106406325A (zh) 基于模糊扩张状态观测器的四旋翼无人机反馈线性化控制方法
CN105182990B (zh) 具有输出受限的三自由度模型直升机的鲁棒控制方法
CN107203138A (zh) 一种输入输出饱和的飞行器鲁棒控制方法
CN105116914A (zh) 一种平流层飞艇解析模型预测路径跟踪控制方法
Yamasaki et al. Integrated guidance and autopilot for a path-following UAV via high-order sliding modes
CN103455035B (zh) 基于反步设计和非线性反馈的pd+姿态控制律设计方法

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: 20171017

Termination date: 20210610

CF01 Termination of patent right due to non-payment of annual fee