您的当前位置:首页正文

基于LESO的高超声速飞行器动态面控制

来源:华佗健康网
2018年2月 西北工业大学学报 Feb.2018 第36卷第1期 Journal of Northwestern Polytechnical University V01.36 No.1 基于LESO的高超声速飞行器动态面控制 张科 ,杨文骏 ,张明环 ,王佩 '2 (1.航天飞行动力学技术重点实验室,陕西西安710072;2.西北工业大学航天学院,陕西西安710072) 摘要:针对含模型不确定性和外部扰动的高超声速飞行器纵向运动模型,提出了一种基于线性扩张 状态观测器的动态面控制方法。采用非线性动态逆技术实现了高超声速飞行器高度和速度通道的解 耦;结合传统反演设计方法,引入一阶低通滤波器求取虚拟控制量的微分,避免了“微分爆炸”的问题; 设计线性扩张状态观测器,实现了对模型不确定性和外部扰动等组成的“总和扰动”的精确估计,显著 提高了系统的扰动抑制能力;利用Lyapunov理论进行闭环稳定性分析。仿真结果表明,所提出控制 器参数配置简单,对参数不确定性和外部扰动等有较强的鲁棒性,具有良好的指令跟踪效果。 关键词:高超声速飞行器;动态面控制;线性扩张状态观测器;扰动抑制;Lyapunov方法 中图分类号:TJ765.2;V448.2 文献标志码:A 文章编号:1000—2758(2018)O1-0013-07 高超声速飞行器(hypersonic flight vehicle, 入LESO实现了对飞行器状态和总干扰的有效估 HFV)飞行速度快、高度高,具有高超声速巡航与突 计。注意到,上述部分研究成果存在算法复杂、参数 防、空天往返等广泛应用前景,成为了许多国家研究 设置繁琐和未考虑滤波误差等问题。 发展的重点¨J。也因其模型具有高非线性、强耦 本文基于文献[16],结合传统反演法和动态面 和参数不确定等特点,对HFv的控制系统设计提出 控制技术,引入一阶低通滤波器(1ow—pass filter, 了更高的挑战 J。 LPr)求取虚拟控制量的微分信号;设计了二阶 增益预置-4 J、滑模控制 和鲁棒控制 等方 LESO,实现对系统总和扰动的精确估计补偿;最后 法,都在HFV的控制系统设计中得到了广泛研究及 通过仿真实例验证了所提出控制方法的有效性。 应用。系统的非匹配不确定性和外部扰动是高超声 速飞行器飞行过程中的突出问题,反演法 副是处 1模型及问题描述 理非匹配不确定性的一种重要非线性控制方法,而 干扰观测器能实现对不确定性及各类扰动的有效估 1.1 lⅡ’V纵向模型 计和补偿 。 。扩张状态观测器是自抗扰控制技 术【I 的核心部分,不依赖精确模型,能实现对系统 HFV纵向模型如下 Tcosa-D不确定性、未建模动态和外部扰动等组成的“总和 :——/zsi一—ny —’ M r2 扰动”的有效估计。针对高超声速飞行器非线性和 L+Tsina易受扰动影响的问题,文献[12]提出了带有扩张状 :——一! 二 ! ‘ Mv vr 态观测器的连续滑模控制方法。文献[13]针对飞 疗: in (1) 行器复杂非线性和不确定性等问题,提出了一种含 扩张状态观测器的动态面控制方案。Gao【14]基于自 d q一 抗扰控制技术,提出了线性扩张状态观测器(1inear 口= extended state 0bserver,LESO),用带宽的概念确定 其参数,结构设计简便,需调参数少。文献[15]引 式中, ,y,H, ,q分别为HFV的速度、航迹角、高 收稿日期:2017-05—12 基金项目:国家自然科学基金(61174204、61101191、61502391)资助 作者简介:张科(1968一),西北工业大学教授,主要从事飞行器导航、制导与控制及深空探测等研究。 14 西北工业大学学报 第36卷 度、攻角和俯仰角速率; 为地心引力常数;r是HFV 到地心的距离且r=RE+日, 为地球半径; 和, 称控制矩阵;d (t)和d2(t)为速度和高度通道的综 合扰动。 系统对AB的敏感程度要远高于△F 和△ , 分别为HFV的质量及其沿Y轴的转动惯量。 ,D, 及 分别为HFV的升力、阻力、推力及俯仰力矩, 故AFy和AF 可以忽略不计。故(6)式可改为 FV: 其具体表达式可参考文献[17]。 简化的发动机模型采用一个二阶系统描述 +[Bo+AB [ 二 ㈩ =一2 ∞, 一∞: + :卢 (2) 式中, , 分别为二阶系统的阻尼比和自然频率, 卢,卢 分别表示发动机节流阀的实际开度及其指 令值。 取系统状态变量X=[V,y,H, ,q,卢, ]T,采用 非线性动态逆技术,对HFV的纵向非线性模型进行 精确反馈线性化,分别对速度 求导3次、高度日求 导4次后,出现控制量(发动机节流阀开度 和升降 舵偏转 ),系统的相对阶为3+4:7正好等于系统 的阶数,故不存在零动态。系统可描述为 ㈥H(4 叫 ㈩ 式中,, , 和曰的详细表达式见参考文 献[16一l7]。 假设1 矩阵 非奇异。 注1 考虑到’HFV在巡航飞行过程中,航迹角 很小,满足 ≠±÷,所以HFV的飞行航迹不是垂直 的,故矩阵曰非奇异,假设1成立。 1.2问题描述 考虑HFV的参数不确定性 r =眠(1+△ ),,"=1o(1+△,) { =S0(1+AS), =co(1+ ) (4) tp=P0(1+△!p),c =0.029 2(1+△c ) 式中,“ 。”表示参数标称值;“A ”表示参数偏 移量。其参数偏移量如下 rI△M I≤0.03,I A,f≤0.02, ?I△Ls l≤0.01,I△ l≤0.01, 【l△p j≤0.06,l△c l≤0.03,l ACM.1≤0.1 (5) 将系统的未知动态和外部扰动等视为系统的等 效综合扰动,考虑以上情况,系统动态可描述为 FH(4:v..+AFv 】+[Bo+AB][ 二 (6) 式中,△, ,△ 和AB为系统不确定项;曰。表示标 将速度和高度通道中的系统不确定性、等效综 合扰动视为“总和扰动” )卜 ㈤ 当AM和△,取最大负值,△p,AS,Ac。和 取最 大正值时,参数不确定性对系统的影响最大。 2基于LESO的动态面控制器设计 2.1动态面控制器设计 将(3)式写成紧凑形式 G=F +Bu (9) 式中 G= “= 可通过结合动态面控制技术的反演法设计得到 G,进而得到动态面控制器 0=Bo (G—F0) (10) 针对高度子系统,定义e =H一 ,E =[e , , , ,e: ’] ,其中 为高度指令信号,故 Ehl eh h1=E = E^2=E^3= (11) E^3=EJl4= ^ ^4=H ’一 ’=ei 定义4个误差变量 式中, (i=1,2,3)为需要设计的虚拟反馈控制 量。(12)式本质上为(11)式的微分同胚,因此镇 定(11)式,只需镇定(12)式即可。 第1步: 第1期 张科,等:基于LESO的高超声速飞行器动态面控制 15 对z 求导,得 : , 2^1=Z + ^l (13) (0): (0)(22) .r 式中, 为低通滤波器的时间常数 为虚拟控制量 针对 定义Lyapunov函数 。=÷z: ,并取虚 的估计值。 厶 拟控制量 ^。=一P Z 选取P^ >0,则有 L^l=一p^lz l+z^lz (14) 要使得 。即 收敛到零,需要满足厶 <0, 则Z 也要趋于零。 第2步: 对z 求导,得 Z =Z^3+ ^2一 hl (15) 1 对z…z^2定义Lyapunov函数L^2: 1+÷z 2, 二 并取虚拟控制量 ^2=一Z^1一Ph2Z^2+ 1,选取P^2> 0,则有 £^2=一phlz:1一Ph2z +z娩z^3 (16) 要使得 。和z :收敛到零,需满足£ :<0,则 。也需收敛到零。 第3步: 对z ,求导,得 Z^3=Z^4+ ^3一 (17) 对z z z^3定义Lyapunov函数L^3=L^2+ l z2 ,取虚拟控制量X^3=一Z^2一Ph3Z^3+ ,选取 P粘>O,则有 £^3:一p lz:l—Ph=z 一p船2 +zJl3zM(18) 。、Z :和z ,收敛到零,需满足£ 3<0,即z 需收敛到零。 第4步: 对 求导,得 Ⅵ=HH 一 一 3 (19) 对 I,Z z^3,z 定义Lyapunov函数LM:L^3 1 +—÷ ,取 Jl4=一pJl4z —z 3,其中P^4>0,贝Ⅱ有 L^4=一p^lz 1一p^2z:2一p^3zl3一Ph4Zh4<0 (20) 故 。,z , 和z 都能收敛到零,这也保证 了HFV的高度跟踪误差能稳定收敛到零。得到所 需的H( ’如下 H¨ = + ^3一pMz^4一Z^3(21) 设计如下一阶低通滤波器求取虚拟控制量的微 分信号,以避免“微分爆炸”的问题。 同理,针对速度子系统,采用相同的方法经3步 设计得到 ,最终可得动态面控制器如(10)式所示。 2.2 LESO设计 结合模型精确反馈线性化的部分信息,针对 HFV的速度和高度子系统,分别设计简单的二阶 LESO对系统的总和扰动进行估计补偿。 以速度子系统为例,令 。= , =d (t)为扩张 的状态(即总和扰动),并记 :=to( ),则将系统 (7)的速度子系统扩张成新的系统 fxl 2+F +bo,lllf +bo,128 {【 2=to(t) (23) y= 。 针对系统(23)设计二阶LESO如下 1 1一 1 { 1= 2一zle1+F +bo,11 ̄ +bo,12t ̄e (24) 【 :一f2e 式中,f (i=1,2)为LESO的可调参数;z 和 分别 是 。和 :的估计值。 取e = 一 ,i=1,2,则LESO相应的观测误 差动态为 P=Ae+Ed (25) el =㈥。 可调参数f 、Z 采用基于带宽的配置方法,则误 差动态的特征多项式可表示为 A(s)=s +f1s+f2=(s+∞。) (26) 只要满足 。>0,则式(26)具有负实部根,保 证了系数矩阵A是Hurwitz稳定的,其中z。=2oJ。,z =∞。2。 假设2 HFV速度和高度通道的总和扰动 d (t)和d (t)连续可微且一阶导数有界。 注2 高超声速飞行器巡航飞行过程中,其参 数不确定性和外部扰动均在一定范围内,故假设2 成立。 定理1[1 。 针对误差系统(25),若系数矩阵 A是Hurwitz稳定的,且假设2成立时,则误差e ,e2 是有界稳定的;选择合适的参数∞。,可使得e。和e 将在有限时间内收敛到包含原点、半径足够小的闭 16 西北工业大学学报 第36卷 球内。 同理,将该LESO应用于高度子系统,参数设置 耋 ÷+ 02' 2+2.5+Mv+Mh = II ,不变,可实现对总和扰动的有效估计。结合(10) 式,最终得到基于LESO的HFV动态面控制器 选取参数 ∑ ∑ 1 一 卜 ●一2 + H= cG—F。一【耋:】, c27 式中,a 和 分别是LESO对总和扰动d (t)和 d (t)的估计值。 2.3稳定性分析 Pq B,Phi ≥ >啪 max { 。, 5 -v-)+ } (34) 2 式中,占为待设计正数。则(33)式可转化为 ∑ ∑ ≤一2 +2.5+ +M^ 2△ 可 + 式中,W=∑÷z;+∑—, ;+∑ =1 厶 =1 厶 1 。一 + + 定义一阶低通滤波器估计误差 ' ∑—÷厶 。 (28) 选取 ≥ 一2 【△m= -X i:1,2,3 ,可使得£≤。,从而 对其求导可得 A :一鱼一 一 :12 ,保证了系统的半全局一致有界稳定性。 3仿真验证及分析 _1’2,3 易知虚拟控制量的一阶导数有界,且满足 以HFV在33 km高空以15马赫巡航飞行为 例,其主要初始参数为:高度日为33 000 m;速度 为4 565.309 m/s;航迹角 为0。;攻角 为1.955。; 俯仰角速率为0。/s;控制输入发动机节流阀开度 为0.181;升降舵偏角 为一0.427。。考虑现实中执 I 『≤D f I≤D ,其中D 和D巧为正常数。 定义总和扰动的估计误差 =ai—d ,i= ,h。 选取Lyapunov函数 L=L +己h (30) 式中 行机构的限制, ∈[0,1],6 ∈[-20。,20。]。 动态面控制器参数设置为P =1,P =1;一阶 低通滤波器参数设置为丁=0.05;LESO参数设置为 ∞ =30,其状态变量的初始状态设为0。 仿真过程中,给予速度通道60 m/s的阶跃信 Ao 一譬一 ≤一譬+Do△ 号,同时在高度通道中每隔100 s给予幅值为100 m 的阶跃信号,参考指令都通过一个二阶滤波器 ) s +2 +0 :, ≤一等+ i:一生一Zh = + 2丢 ≤一堡+D△ hi(35) 式中,阻尼比 r=0.9,自然频率 r=0.15 rad/s。 从仿真的第170 s开始引入等效综合扰动d (t) =△ 0.05cos(0.02-at),d2(t)=0.05sin(0.02,trt)。 ≤一生+等 =—q-—+— ——— —一= +I、 一—7— ,+ 2— I △ +-+丢 = ≤一为验证所提出方法的有效性,在仿真的170 s处 引入等效综合扰动d (£)=0.05cos(0.02,rrt),d2(£)= 0.05sin(0.02,rrt)。将本文所设计方法和滑模控制 方法进行对比,下标“ref”表示参考指令,下标“1” 表示滑模控制器+本文LESO,下标“2”表示本文方 法。仿真结果如图1~图5所示。 根据定理1,可知 a 和 是有界且收敛的, 分别取l a I和I a l的上界为M 和M , 和 M 均为正常数,则对L求导可得 £≤一 3 J = 1 p z;+耋(= 1 、 一 丁1+ 二譬1 i/ △ 2一 4=1 p z + 第1期 张科,等:基于LESO的高超声速飞行器动态面控制 17 3-36 : 3_30 4 l∞ 一自4  4 t/s b) 速度响应曲线 t/s b)速度跟踪误差 t/s b)航迹角响应曲线 图1高度、速度响应曲线 图2高度、速度跟踪误差 图3攻角、航迹角响应曲线 图1和图2分别为HFV的高度和速度指令跟 踪曲线和跟踪误差曲线。可以看出,采用本文提出 的控制方案,响应速度更快且跟踪误差小、收敛快, 图5中给出的是本文所设计控制器作用下的 LESO对总和扰动的估计曲线,可以看出在未引入 等效综合扰动前,LESO能够对模型不确定项进行 不到20 s就能稳定跟踪上指令信号;而滑模控制器 则需要50 s左右才能稳定跟踪上指令信号。在 LESO的估计补偿作用下,2种方案均对模型不确定 项、外部扰动等具有较强的鲁棒性。图3为攻角和 航迹角的响应曲线,可以看出,在170 s处引入等效 有效估计,引入等效综合扰动后,LESO能够实现对 模型不确定项和综合扰动组成的”总和扰动”的有 效估计。 扰动时,攻角和航迹角受到影响而偏离平衡状态,但 在控制器和LESO对扰动有效估计的作用下均能快 速收敛到平衡状态。 a) 速度予系统t/总和扰动估计 S 图4为控制量曲线,可以看出滑模控制器的控 制量全程伴随着高频抖振,这在工程实际中会造成 极大影响,而本文控制器则不存在这个问题。引入 等效综合扰动后,控制量也随之波动,这样的振荡是 为了使HFV能够抑制扰动作用的影响,保持稳定的 飞行状态。 t/s b) 高度子系统总和扰动估计 图5总和扰动估计 4结 论 针对HFV巡航飞行时存在模型不确定性和外 部扰动等问题,设计了一种基于线性扩张状态观测 器的动态面控制器。基于传统反演设计方法,引入 阶低通滤波器来求取虚拟控制量的微分,避免了 传统反演法中“微分爆炸”问题,简化了设计步骤; 设计了二阶LESO,实现了对总和扰动的有效估计及 一t/s b)Case1和Case 2的 变化曲线 图4控制量 补偿。仿真结果表明,所设计LESO能够实现对”总 和扰动”的有效估计,大幅提升了控制器的扰动抑 制能力和系统的鲁棒性;和滑模控制相比,本文提出 的动态面控制方案具有更好的高度和速度指令跟踪 西北工业大学学报 第36卷 性能,且该方法设计结构简单,参数设置简便,易于 实现,在工程应用上具有一定参考价值。 参考文献: [1]方洋旺,柴栋,毛东辉,等.吸气式高超声速飞行器制导与控制研究现状及发展趋势[J].航空学报,2014,35(7): 1776.1786 Fang Wangyang,Chai Dong,Mao Donghui,et a1.Status and Development Trend of the Guidance and Control for Air.Breathing Hypersonic Vehicle[J3.Acta Aeonautriea et Astornautica Sinica,2014,35(7):1776.1786(in Chinese) [2]Xu B,Shi Z K.An Ovevirew on Flight Dynamics and Control Approaches for Hypersonic Vehicles[J].Science China Informa. tion Sciences,2015,58(7):1-19 [3]孙长银,穆朝絮,余瑶.近空间高超声速飞行器控制的几个科学问题研究[J].自动化学报,2013,39(11):1901—1913 Sun Changyin,Mu Chaoxu,Yu Yao.Some Control Problems for Near Space Hypersonic Vehicles[J].Acta Automation Sinica。 2013,39(11):1901—1913(in Chinese) [4]Hou Y Z,Wang Q,Dong C Y.Gain Scheduled Control:Switched Polytopic System Approach[J].Journal of Guidance,Control and Dynamics,201l,34(2):623—629 [5]Zhang Y Y,Li R F,Xue T,et a1.An Analysis of the Stability and Chattering Reduction of High-Order Sliding Mode Tracking Control orf a Hypersonic Vehicle[J].Information Sciences,2016。348:25-48 [6]米鹏,罗建军,苏二龙.滑翔式高超声速飞行器H 回路成形控制[J].西北工业大学学报,2013,31(4):565.570 Mi Peng,Luo Jinjaun,Su Eflong.A Novel Gliding nd aEffective Hypersonic Vehicle Control Method Based on H Loop Shaping [J].Journal of Northwestern Polytechnieal University,2013,31(4):565—570(in Chinese) [7] Wang Q,Hou D L,Dong C Y,et a1.Reduced—Order Observer·Based Backstepping Tracking Control for a Class of Stochastic Nonlinear Systems[J].Nonlinear Dynamics,2014,78(2):939—956 [8] Bu X W,Wu X Y,Zhang R,et a1.Tracking Differentiator Design ofr the Robust Backstepping Control of a Flexible Air-Breath. ing Hypersonic Vehicle[J].Journal of the Franklin Institute,2015,352(4):1739—1765 [9]Wu H N,Feng S,Liu z Y,et a1.Disturbance Obseverr Based Robust Mixed H2/H Fuzzy Tracking Control ofr Hypesonirc Ve— hicles[J].Fuzzy Sets and Systems,2017,306:118—136 [10]Wang F,Zou Q,Hua C C,et a1.Disturbance Observer—Based Dynamic Surface Control Design ofr a Hypersonic Vehicle with In— put Constraints and Uncertainty[J].Proceedings of the Institution of Mechanical Engineers,Part I:Journal of Systems and Con— trol Engineering,2016,230(6):522-536 [11]韩京清.自抗扰控制技术——估计补偿不确定因素的控制技术[M].北京:国防工业出版社,2008 Han Jingqing.Active Disturbance R ̄eetion Control Technique一一the Technique for Estimating nd aCompensating the Uncertainties [M].Beijing,National Defense Industry Press,2008(in Chinese) [12]吴云沽,王建敏,刘晓东,等.带有干扰观测器的高超声速飞行器滑模控制[J].控制理论与应用,2015,32(6):717—724 Wu Yunjie,Wang Jianmin,Liu Xiaodong,et a1.Disturbance—Observer—based Sliding Mode Control for Hypersonic Flight Vehicle [J].Control Theoy&Applrications,2015,32(6):717-724(in Chinese) [13]刘晓东,黄万伟,禹春梅.含扩张状态观测器的高超声速飞行器动态面姿态控制[J].宇航学报,2015,36(8):916—922 Liu Xiaodong,Huang Wanwei,Yu Chunmei.Dyanamic Surface Attitude Control for Hypersonic Vehicle Contmning Extended State Observer[J].Journal ofAstronautics,2015,36(8):916—922(in Chinese) [14]Gao Z Q.Active Disturbance Rejection Control:A Paradigm Shift in Feedback Control System Design[C]∥Proceedings of the 2006 American Control Conference,Minneapolis,Minnesota,2006:2399—2405 [15]张银辉,杨华波,江振宇,等.基于干扰估计的高超声速飞行器鲁棒控制方法[J].系统工程与电子技术,2016,38(04): 875.881 Zhang Yinhui,Yang Huabo,Jiang Zhenyu,et a1.Roust Control of Hypersonic Flight Vehicle Based on Disturbance Estimation [J].System Engineering and Electonircs,2016,38(4):875—881(in Chinese) [16]杨文骏,张科,张明环,等.基于ESO的高超声速飞行器非线性动态逆控制[J].西北工业大学学报,2016,34(5): 805.81l Yang Wenjun,Zhang Ke,Zhang Minghuan,et 1.ESO aBased Nonlinear Dynamic Inversion Control ofr Hypersonic Flight Vehi— 第1期 张科,等:基于LESO的高超声速飞行器动态面控制 19 cle[J].Journal of Northwestern Polytechnieal University,2016,34(5):805-811(in Chinese) [17]xu H J,Mirmirani M D,Ioannou P A.Adaptive Sliding Mode Control Design for a Hypersonic Flight Vehicle[J].Joumal of Guidance。Control,and Dynamics。2004,27(5):829-838 [18]Guo B Z,Zhao Z L.On the Convergence of an Extended State Observer ofr Nonlinear Systems with Uncertainty[J].Systems& Control Letters,201 1,60(6):420-430 [19]邵星灵,王宏伦.线性扩张状态观测器及其高阶形式的性能分析[J].控制与决策,2015,30(5):815—822 Shao Xingling,Wang Honglun.Performance Analysis on Linear Extended State Observer and its Extension Case with Hihegr Ex- tended Order[J].Control and Decision,2015,30(5):815-822(in Chinese) LESO Based Dynamic Surface Control for Hypersonic Flight Vehicle Zhang Ke ,Yang Wenjun ,Zhang Minghuan ,Wang Pei , ,1.National Key Laboratory of Aerospace Flight Dynamics,Xi all 710072,China; 、 \2.School of Astronautics,Northwestern Polytechnieal University,Xi an 710072,China/ Abstract:In the presence of model parametirc uncertainties and external disturbances,a LESO—based dynamic sur— face control approach is designed for the longitudinal model of Hypersonic Flight Vehicle(HFV).Via Nonlinear Dynamic Inversion(NDI)technique,the decoupling of altitude and velocity is realized.Combining with conven- tional back-stepping technique,a low pass ifher(IJPF)is introduced to attain the derivation of vitrual control laws, which avoids the problem of“diferentiation explosion”.A linear extended state observer(LESO)is designed for the precise estimation and compensation of”lumped disturbance”containing parametric uncertainties and external disturbances.which tremendously improves the ability of disturbance rejection of the system.The stability of the proposed approach is analyzed by means of Lyapunov theory.11}le simulation results demonstrate that the proposed methodology has good command tracking performance,and the approach is robust in the presence of lumped dis— turbances. Keywords:hypersonic vehicles;dynamic surface control;linear extended state observer;disturbance rejection; Lyapunov methods 

因篇幅问题不能全部显示,请点此查看更多更全内容