CN111007867B - 一种可预设调整时间的高超声速飞行器姿态控制设计方法 - Google Patents
一种可预设调整时间的高超声速飞行器姿态控制设计方法 Download PDFInfo
- Publication number
- CN111007867B CN111007867B CN201911399686.5A CN201911399686A CN111007867B CN 111007867 B CN111007867 B CN 111007867B CN 201911399686 A CN201911399686 A CN 201911399686A CN 111007867 B CN111007867 B CN 111007867B
- Authority
- CN
- China
- Prior art keywords
- aircraft
- control
- attitude
- angle
- 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.)
- Active
Links
- 238000013461 design Methods 0.000 title claims abstract description 33
- 238000000034 method Methods 0.000 title claims abstract description 21
- 238000013178 mathematical model Methods 0.000 claims abstract description 8
- 238000004088 simulation Methods 0.000 claims abstract description 6
- 238000007689 inspection Methods 0.000 claims abstract description 3
- 230000003044 adaptive effect Effects 0.000 claims description 15
- 230000006978 adaptation Effects 0.000 claims description 2
- 239000011159 matrix material Substances 0.000 claims description 2
- 238000005096 rolling process Methods 0.000 claims 2
- 230000033001 locomotion Effects 0.000 claims 1
- 230000001133 acceleration Effects 0.000 description 3
- 238000005215 recombination Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000011056 performance test Methods 0.000 description 2
- 230000006798 recombination Effects 0.000 description 2
- 239000000758 substrate Substances 0.000 description 2
- 241000287196 Asthenes Species 0.000 description 1
- 101100391182 Dictyostelium discoideum forI gene Proteins 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- PHTXVQQRWJXYPP-UHFFFAOYSA-N ethyltrifluoromethylaminoindane Chemical compound C1=C(C(F)(F)F)C=C2CC(NCC)CC2=C1 PHTXVQQRWJXYPP-UHFFFAOYSA-N 0.000 description 1
- 208000011580 syndromic disease Diseases 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Landscapes
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
- Feedback Control In General (AREA)
Abstract
本发明记载一种可预设调整时间的高超声速飞行器姿态控制设计方法,属于航空航天领域,具体方案包括以下步骤:步骤一、明确控制系统的设计任务;步骤二、建立高超声速飞行器姿态系统数学模型;步骤三、设计自适应滑模控制律;步骤四、闭环系统分析;步骤五、利用计算机数值仿真工具Matlab/Simulink进行闭环系统的性能检验。本设计方法不仅使得高超声速飞行器姿态角在有限时间内达到了精度要求,而且可以根据性能指标要求预先设定需要的调整时间。此外,该控制律还能使不确定项估计值按需增长,也减小了控制器增益,很大程度上克服了一般鲁棒控制和自适应控制设计的保守性。
Description
技术领域
本发明属于航空航天领域,具体涉及一种可预设调整时间的、可使高超声速飞行器姿态角误差在指定时间内收敛到指定精度的自适应滑模控制律设计方法。
背景技术
在高超声速飞行器的控制中,需要保证飞行器的姿态角跟踪参考指令信号,否则难以实现飞行器的飞行任务。目前对于高超声速飞行器的控制,大多只是保证姿态角在无限时间内收敛到指定参考值,少数可以保证在有限时间内收敛,但无法预先设定调整时间。为解决控制中的不确定性问题,目前应用的大多数鲁棒控制律一般需要选取较大的控制增益;而大多自适应控制律也存在估计值持续增长的问题。
发明内容
本发明的目的是提供一种可预设调整时间的高超声速飞行器姿态控制设计方法,该方法能使系统姿态角误差在预设调整时间内达到指定精度要求;而且该方法能解决多数鲁棒控制律控制增益较大和多数自适应律估计值持续增长的问题,使控制器的增益按需增长。
为了实现上述目的,本发明采取的技术方案如下:
一种可预设调整时间的高超声速飞行器姿态控制设计方法,包括以下步骤:
步骤一、明确控制系统的设计任务,使高超声速飞行器姿态角在预定时间内跟踪参考指令信号,并具有良好的鲁棒性和对系统扰动及外界干扰的自适应能力;
步骤二、建立高超声速飞行器姿态系统数学模型;
步骤三、设计自适应滑模控制律;
步骤四、闭环系统分析;
步骤五、利用计算机数值仿真工具Matlab/Simulink进行闭环系统的性能检验。
进一步的,步骤一中,控制系统的设计任务是:给定参考信号xd,预设的调整时间Tf和允许的跟踪误差∈,设计合适的控制律使得闭环系统的状态x有界,并且当t≥Tf时,系统状态的跟踪误差|x-xd|≤∈。
进一步的,步骤二中,建立高超声速飞行器姿态系统数学模型的具体步骤如下:描述高超声速飞行器三通道姿态的二阶非线性微分方程如下所示:
式中,分别为飞行器俯仰角、偏航角和滚转角,分别为飞行器俯仰角速度、偏航角速度和滚转角速度;分别为飞行器俯仰角加速度、偏航角加速度和滚转角加速度,Jx,Jy,Jz分别为飞行器对弹体坐标系各轴的转动惯量,Mx,My,Mz分别为作用在飞行器上的合外力对质心的力矩在弹体坐标系各轴上的分量,计算公式为:
Mx=qsllatmx
My=qsllatmy
Mz=qsllonmz
式中,q,s分别表示动压和参考面积,llat,llon分别为侧向和纵向参考长度;气动力矩系数mx,my,mz为空速V,迎角α,侧滑角β和姿态角速率ωx,ωy,ωz以及舵偏角δx,δy,δz的函数,即:
mx=mx(V,α,β,ωx,ωy,ωz,δx,δy,δz)
my=my(V,α,β,ωx,ωy,ωz,δx,δy,δz)
mz=mz(V,α,β,ωx,ωy,ωz,δx,δy,δz)
下面定义系统状态变量和控制量为:
根据式(1)-(3)即可得到飞行器姿态控制系统二阶非线性模型:
G(x)=qsLB(x)T
其中,δe、δa和δr分别为左右升降舵偏转角和方向舵偏转角;L=diag(llat,llat,llon)为飞行器参考长度构成的矩阵;
其中,
进一步的,步骤三中,设计自适应滑模控制律的具体步骤如下:
步骤1:定义滑模面
首先,定义系统实际跟踪误差轨迹与期望跟踪误差轨迹之差
z=e-η
式中,z=zi=[z1,z2,z3]T,z为时间t的函数,e为系统跟踪误差;η为跟踪误差的预期轨迹,其满足以下条件:
(1)η在区间[0,∞)上有二阶连续导数;
(4)当t>Tf时,η=0,Tf为设定的调整时间;
η依据下列公式拟合:
式中
其中,
定义滑模面为
容易得到,该滑模面具有下列性质:
根据选取的滑模面,系统自适应滑模控制律设计如下:
μi>0为待设计参数;
定义估计误差为
所以有
进一步的,步骤四中,闭环系统分析的具体步骤如下:
飞行器姿态控制系统二阶非线性模型为:
那么,当t≥Tf时,各通道状态实际值xi与状态期望值xid满足:|xi-xid|≤∈i(i=1,2,3)。
证明:将控制律(7)代入式(6)得
考虑第i个子系统
定义
当|si|>ε,有
下面讨论si(t)的变化情况,si(0)=0,从t=0开始,si(t)值的变化有两种情况;第一种理想情况对于都有|si(t)|<εi,因此,对所以,另一种情况为存在时间T1>0使得|si(t)|≤εi,t∈[0,T1),|si(T1)|=εi,并且当t>T1时si(t)超出区间[-εi,εi]。由|si(t)|≤εi,t∈[0,T1)得,因此再结合条件|si(T1)|=εi可知
由式(9)得,当si(t)在区间[-εi,εi]之外时,Vi(t)严格单调递减,所以,si(t)不可能一直保持在区间[-εi,εi]之外;假设成立,那么有在时成立;所以成立,即当时|si(t)|≤εi成立,与假设相矛盾,所以一定存在时间T2>T1使得|si(T2)|=εi并且当t>T2时si(t)回到区间[-εi,εi]内;又Vi(t)在[T1,T2]上单调递减,所以因此当t∈[T1,T2]时,一定有成立,并且Vi(T2)<Vi(T1)<εi 2,因此
再结合性质1可得
等价于
因为当t≥Tf时ηi(t)=0,所以有
|xi-xid|<∈i,t>Tf
并且,
综上,该控制律达到了设计目标,能够使系统姿态角误差在预设调整时间内达到指定精度要求;控制器增益不会过大;并且不存在估计值持续增长的问题。
本发明的优点:提供了一种可预设调整时间的高超声速飞行器姿态控制设计方法,不仅使得高超声速飞行器姿态角在有限时间内达到了精度要求,而且可以根据性能指标要求预先设定需要的调整时间。此外,该控制律还能使不确定项估计值按需增长,从而也减小了控制器增益,一定程度上克服了一般鲁棒控制和自适应控制设计的保守性。
附图说明
图1设计方法流程图;
图2控制系统结构框图;
图3俯仰角及其参考信号;
图4偏航角及其参考信号;
图5滚转角及其参考信号;
图6自适应律估计值。
具体实施方式
具体实施方式一
一种可预设调整时间的高超声速飞行器姿态控制设计方法,包括以下步骤:
步骤一:明确控制系统的设计任务,给出了控制系统设计的目的是使高超声速飞行器姿态角信号在预定时间内跟踪参考指令信号;步骤二:建立高超声速飞行器姿态系统二阶非线性数学模型;步骤三:定义了滑模面,基于滑模面设计了自适应滑模控制律;步骤四:借助计算机数值仿真工具Matlab/Simulink进行了闭环系统的性能检验。经过上述步骤,设计结束。
具体步骤如下:
步骤一、明确控制系统的设计任务
控制系统的设计任务是:给定参考信号xd,预设的调整时间Tf=2s和允许的跟踪误差∈=0.5°,设计合适的控制律使得闭环系统的状态x有界,并且当t≥Tf时,系统状态的跟踪误差|x-xd|≤∈。
步骤二、建立高超声速飞行器姿态系统数学模型
飞行器姿态系统二阶非线性数学模型为
其中
模型中涉及到的变量说明如下:di>0为未知有界常数,取其保守估计值飞行器的侧向和纵向参考长度llat,llon分别为24.384m和18.288m;参考面积s为334.73m2;详细的空气动力学参数计算以及气动力和力矩的公式参考文献(侯明哲,谭峰.高超声速飞行器的增益协调鲁棒参数化控制[M].北京:科学出版社,2018.)。飞行器的初始姿态设置为γ(0)=3°,ψ(0)=2°,ωx=10°/s,ωy=10°/s,ωz=10°/s;为了验证所设计控制律的鲁棒性,假设实际气动力矩系数mx、my和mz的值都比它们的标称值增加了30%,实际转动惯量Jx、Jy和Jz的值均较其标称值增加20%。
步骤三、设计自适应滑模控制律
设计过程分为两小步:
步骤1:定义滑模面
定义滑模面为
步骤2:设计自适应滑模控制律
根据选取的滑模面,系统自适应滑模控制律设计如下:
步骤四:闭环系统分析及检验
借助计算机数值仿真工具Matlab/Simulink进行闭环系统的性能检验。
系统的控制结构框图如图2所示,系统以期望的飞行器姿态角与实际姿态角的偏差作为控制输入,控制器采用设计的自适应滑模控制律,依据相应的状态偏差解算出需要的控制力,使执行机构产生相应的舵偏,调整飞行器姿态角,使得实际姿态角收敛于期望姿态角。
最终得到的控制效果如图3-5所示,(其中参考指令信号xd由制导信号生成),图3中,飞行器的俯仰角在预设调整时间2s内跟踪上参考指令信号,动态性能良好,跟踪误差小于0.5°,达到了设计目标;图4中飞行器的偏航角在预设调整时间2s内跟踪上参考指令信号,动态性能良好,跟踪误差小于0.5°,达到了设计目标;图5中,飞行器的滚转角在预设调整时间2s内跟踪上参考指令信号,动态性能良好,跟踪误差小于0.5°,达到了设计目标。
图6为系统不确定项估计值,由图可见系统不确定项估计值有界且值较小,避免了估计值持续增长的现象,达到了设计目标。
综上所述,由仿真结果可看出在系统控制增益k较小的情况下,飞行器姿态角可在预先设定的2s内实现对指令信号的跟踪,跟踪误差稳定在零附近,能够达到精度要求,并且不确定项估计值有界,没有持续增长现象,很好的实现了设计要求。
Claims (3)
1.一种可预设调整时间的高超声速飞行器姿态控制设计方法,其特征在于,包括以下步骤:
步骤一、明确控制系统的设计任务,使高超声速飞行器姿态角在预定时间内跟踪状态参考指令信号,并具有良好的鲁棒性和对系统扰动及外界干扰的自适应能力;
步骤二、建立高超声速飞行器姿态系统数学模型;定义系统状态变量x和控制量u为:
式中,θ,ψ,γ分别为飞行器俯仰角、偏航角和滚转角,δx,δy,δz分别为控制飞行器进行滚转、偏航和俯仰运动的舵偏角;
建立高超声速飞行器姿态系统数学模型
其中,
G(x)=qsLB(x)T
式中,q,s分别表示动压和参考面积,L=diag(llat,llat,llon)为飞行器参考长度构成的矩阵,llat,llon分别为侧向和纵向参考长度;δe、δa和δr分别为左右升降舵偏转角和方向舵偏转角;
其中,
mx=mx(V,α,β,ωx,ωy,ωz,δx,δy,δz)
my=my(V,α,β,ωx,ωy,ωz,δx,δy,δz)
mz=mz(V,α,β,ωx,ωy,ωz,δx,δy,δz)
式中,气动力矩系数mx,my,mz为空速V,迎角α,侧滑角β和姿态角速率ωx,ωy,ωz以及舵偏角δx,δy,δz的函数;
式中,Jx,Jy,Jz分别为飞行器对弹体坐标系各轴的转动惯量;
其中,
步骤三、定义滑模面,基于滑模面设计自适应滑模控制律;具体步骤如下:
步骤1:定义滑模面为
步骤2:设计自适应滑模控制律
根据选取的滑模面,系统自适应滑模控制律设计如下:
式中,x=xi=[x1,x2,x3]T=[γ,ψ,θ]T(i=1,2,3),为状态变量,k=[k1,k2,k3]T,其中ki>0,i=1,2,3为待设计参数,ε=[ε1,ε2,ε3]T为一常数向量,SATε(s)为饱和函数,其中εi满足不等式∈i为系统期望的跟踪误差;si代表第i个子系统的滑模面;
μi>0为待设计参数;
步骤四、闭环系统分析;
步骤五、利用计算机数值仿真工具Matlab/Simulink进行闭环系统的性能检验。
2.根据权利要求1所述的一种可预设调整时间的高超声速飞行器姿态控制设计方法,其特征在于,步骤一中,所述控制系统的设计任务是:给定状态参考信号xd,预设的调整时间Tf和允许的跟踪误差∈,设计合适的控制律使得闭环系统的状态变量x有界,并且当t≥Tf时,系统状态的跟踪误差|x-xd|≤∈。
3.根据权利要求1所述的一种可预设调整时间的高超声速飞行器姿态控制设计方法,其特征在于,步骤四中,闭环系统分析的具体步骤如下:
飞行器姿态控制系统二阶非线性模型为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911399686.5A CN111007867B (zh) | 2019-12-30 | 2019-12-30 | 一种可预设调整时间的高超声速飞行器姿态控制设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911399686.5A CN111007867B (zh) | 2019-12-30 | 2019-12-30 | 一种可预设调整时间的高超声速飞行器姿态控制设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111007867A CN111007867A (zh) | 2020-04-14 |
CN111007867B true CN111007867B (zh) | 2021-01-29 |
Family
ID=70119681
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911399686.5A Active CN111007867B (zh) | 2019-12-30 | 2019-12-30 | 一种可预设调整时间的高超声速飞行器姿态控制设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111007867B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112082549B (zh) * | 2020-09-10 | 2023-01-10 | 中国人民解放军海军航空大学 | 一种仅测量加速度的飞行器简单质心控制方法 |
CN113325861B (zh) * | 2021-06-02 | 2023-03-24 | 上海海事大学 | 一种非奇异预设时间四旋翼无人机姿态跟踪控制方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6634594B1 (en) * | 2002-05-03 | 2003-10-21 | The Boeing Company | Hypersonic waverider variable leading edge flaps |
CN102880060B (zh) * | 2012-10-25 | 2014-09-10 | 北京理工大学 | 再入飞行器自适应指数时变滑模姿态控制方法 |
CN105116905A (zh) * | 2015-05-26 | 2015-12-02 | 芜湖航飞科技股份有限公司 | 一种飞行器姿态控制方法 |
CN109977613B (zh) * | 2019-04-19 | 2021-01-01 | 哈尔滨工业大学 | 一种可预先设定调整时间的自适应滑模末制导律设计方法 |
-
2019
- 2019-12-30 CN CN201911399686.5A patent/CN111007867B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN111007867A (zh) | 2020-04-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110377045B (zh) | 一种基于抗干扰技术的飞行器全剖面控制方法 | |
CN106444799B (zh) | 基于模糊扩张状态观测器和自适应滑模的四旋翼无人机控制方法 | |
CN109782795B (zh) | 一种利用耦合的面对称高超声速飞行器横侧向控制方法及控制系统 | |
CN111290421A (zh) | 一种考虑输入饱和的高超声速飞行器姿态控制方法 | |
CN111367182A (zh) | 考虑输入受限的高超声速飞行器抗干扰反步控制方法 | |
CN105629734B (zh) | 一种近空间飞行器的轨迹跟踪控制方法 | |
CN103558857A (zh) | 一种btt飞行器的分布式复合抗干扰姿态控制方法 | |
CN105607473B (zh) | 小型无人直升机的姿态误差快速收敛自适应控制方法 | |
CN107807657B (zh) | 一种基于路径规划的挠性航天器姿态自适应控制方法 | |
CN111831002B (zh) | 一种基于预设性能的超高声速飞行器姿态控制方法 | |
CN102880060A (zh) | 再入飞行器自适应指数时变滑模姿态控制方法 | |
CN104331084B (zh) | 一种基于方向舵控滚转策略的气动舵偏范围计算方法 | |
CN112486193B (zh) | 一种基于自适应增广控制理论的飞翼无人机三轴全权限控制方法 | |
CN110162071B (zh) | 一种高超声速飞行器再入末段姿态控制方法及系统 | |
CN111007867B (zh) | 一种可预设调整时间的高超声速飞行器姿态控制设计方法 | |
CN110244751B (zh) | 一种高超声速飞行器姿态自适应递推控制方法及系统 | |
CN106444822A (zh) | 一种基于空间矢量场制导的平流层飞艇路径跟踪控制方法 | |
CN110362110B (zh) | 一种固定时自适应神经网络无人机航迹角控制方法 | |
CN114281092A (zh) | 一种基于滑模干扰观测器的高超声速飞行器协调姿态控制方法 | |
CN113467245A (zh) | 飞行器的分数阶滑模控制方法、装置及系统 | |
CN105182990A (zh) | 具有输出受限的三自由度模型直升机的鲁棒控制方法 | |
CN112947498B (zh) | 飞行器航迹角控制方法、系统及存储介质 | |
CN110347036A (zh) | 基于模糊滑模控制的无人机自主抗风智能控制方法 | |
CN116382332B (zh) | 一种基于ude的战斗机大机动鲁棒飞行控制方法 | |
CN112445235B (zh) | 一种应用于高动态飞行器的滚转稳定控制方法及系统 |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |