一种非仿射高超声速飞行器输出反馈控制方法

路 遥

高超声速飞行器(Hypersonic flight vehicle,HFV)一般是指飞行速度超过5 倍声速的有翼或无翼飞行器.它具有速度快、反应时间短、突防能力强等特点,能够有效提高远程作战效能,具有重要的军事意义,因此成为当今世界军事强国所广泛关注的焦点[1-5].

控制系统作为HFV 的重要分系统之一,是飞行器安全飞行、完成既定任务的重要保证.自上世纪90 年代以来,许多学者致力于HFV 控制器设计研究.然而,HFV 具有飞行包线大、飞行特性复杂、外部环境多变、动态特性存在强非线性和不确定性、气动与推进系统之间存在严重耦合等特点,为其设计控制器具有很大的挑战性.为了降低控制器设计工作的复杂度,部分学者选择首先针对HFV 纵向通道模型进行控制器设计研究.但即使如此,HFV纵向通道数学模型仍是复杂的.为此,学者们通过对HFV 纵向模型数据进行研究分析,通过一定的假设条件,将HFV 纵向模型简化为最小相位系统,或严反馈仿射非线性系统等相对容易进行控制器设计的形式,然后基于自适应控制[6]、backstepping[7-8]、滑模控制[9]、鲁棒控制[10]等较成熟的设计工具设计HFV 控制器,取得了很好的控制效果.

近年来,为提高控制器的实际应用前景,许多学者放弃部分模型简化假设,针对非仿射形式HFV模型进行控制器设计研究.部分学者吸取控制理论研究中已有的针对非仿射非线性系统的相关成果[11-15],将其应用于非仿射HFV 控制器设计中.如文献[16]基于文献[15]中的方法,针对HFV 纵向短周期姿态控制系统设计了非线性跟踪控制器,取得了一定的控制效果.文献[17]考虑执行机构中存在未知非线性动态,同样基于文献[15]中的方法,针对变几何进气道HFV 设计了容错控制器,实现了对速度和航迹倾角的稳定跟踪.但是,由于针对非仿射非线性系统的理论研究成果通常基于特定的非线性模型和一定的假设,而非仿射HFV 模型具有一些独特特性,因此在将理论成果直接应用于HFV 实际模型中时,会被迫对HFV 模型进行一定的限制以满足理论成果的要求,如文献[16]要求飞行器的飞行工况必须保持恒定的航迹倾角,文献[17]设计的控制器难以保证俯仰角速度变化较快情况下的控制效果.因此,在非仿射HFV 控制器设计研究中,不能简单地套用已有的先进理论控制方法,需针对HFV 模型的特点设计合适的、适用性更广的飞行控制器.其他针对非仿射HFV 模型控制器设计研究中,文献[18]针对HFV 轨迹跟踪问题,设计了一种预设性能backstepping 控制器,采用中值定理处理模型中的非仿射项,保证了闭环系统具有较好的跟踪精度.文献[19]提出一种模糊自适应姿态控制方法,将模型中的非仿射项视为总不确定项中的一部分,仿真结果表明该方法具有较好的鲁棒性.文献[20]提出一种基于动态逆的backstepping 姿态控制方法,设计扩张状态观测器估计含非仿射特性的不确定项,取得了良好的控制效果.文献[21]设计了一种预设性能控制器,通过一定的合理性假设将模型中的非仿射项转化为仿射形式设计控制器,仿真结果验证了所提方法的有效性.文献[22]设计了一种基于神经网络的预设性能控制器,利用神经网络估计模型中的非仿射项,取得了良好的跟踪控制效果.

文献[16-22]中所提方法均取得了良好的控制效果.然而,以上方法均需使用攻角、航迹倾角等气流角的测量值构建控制律,这使得上述方法难以应用于HFV 工程实际中.在HFV 实际飞行过程中,相较速度、高度、俯仰角和俯仰角速度,攻角、航迹倾角气流角是难以准确测量的[23-26].这是由于在高超声速飞行状态下,气动加热作用很强,气流角传感器的性能会显著下降.此时,如果控制律设计中仍直接采用气流角的测量信息会带来较大的控制误差,严重时会导致闭环系统不稳定[27].虽然采用新型的埋装式大气数据传感器可在一定程度上提高气流角的测量精度,但这种方案需要在飞行器机身上布置多个压力孔和传感器,代价较高[28].因此,在HFV 控制器设计中,应避免使用航迹倾角或攻角的状态量构建控制器.对于这一点,已有部分文献针对仿射化HFV 模型进行了研究[23-26].然而,文献[23-24]设计的针对仿射化HFV 非线性模型的控制器结构比较复杂,实际应用时会严重影响飞控软件计算耗时,且未考虑执行机构饱和特性;文献[25-26]均需对仿射化HFV 模型进行线性化,而HFV 模型具有强非线性和模型参数不确定的特点,不利于使用小扰动线性化或反馈线性化方法进行研究.根据工程实际情况,直接针对非仿射HFV非线性模型,在不使用航迹倾角和攻角模型输出值的情况下进行控制器设计研究是非常必要的.此外,受产品物理实现的限制,HFV 的控制舵面存在饱和特性,同时在飞行过程中执行机构可能因故障激发某些未知非线性动态,这些也是在控制器设计过程中必须予以考虑的.

以前述需求为牵引,本文针对非仿射HFV 纵向通道非线性动力学模型,在不使用航迹倾角和攻角模型输出值以及考虑执行机构存在未知非线性动态和饱和特性的情况下,研究其轨迹跟踪控制问题,提出一种基于backstepping 的输出反馈控制方法.基于飞行过程中可得到的相对精确的速度和高度测量值,在线估计航迹倾角状态用于控制律设计中.通过跟踪微分器和高阶微分器在线估计信号的一阶、二阶导数,用于估计模型干扰项和解决backstepping 应用中存在的 “微分项爆炸”问题.设计辅助变量降低控制量饱和对控制性能的影响.基于Lyapunov 理论证明闭环系统的稳定性.最后,通过对比仿真试验验证所设计控制器的有效性.

1 非仿射形式HFV 模型描述

1.1 HFV 纵向通道非线性模型

本文选取文献[29]提出的一类吸气式HFV 作为研究对象,其巡航段纵向通道非线性动力学模型可按下式描述:

式中,V、γ、θ、Q为状态量,分别表示飞行速度、航迹倾角、俯仰角和俯仰角速度;α表示攻角,在纵向通道内可按α=θ-γ计算;m、g、Iyy分别表示飞行器质量、重力加速度和俯仰通道转动惯量;dV、dγ、dQ表示由于模型参数不确定性、模型简化以及外部环境干扰引起的扰动项;T、D、L、M分别表示发动机推力、阻力、升力和俯仰力矩,可按以下表达式进行计算[29]:

式中,φ、δe为控制量,分别表示发动机燃料空气混合比和升降舵偏角;分别表示动压、机身参考面积、推力-俯仰力矩耦合系数和平均气动弦长;CT,φ、CT、CD、CL和CM表示推力、气动力和气动力矩系数.由式(2)可知,升降舵对气动阻力D和升力L均有耦合作用,这增加了控制器设计的复杂性.为此,文献[29]提出一种增加鸭翼舵,设定鸭翼舵摆角δc与升降舵摆角δe按固定比例联动,以消除升降舵对升力影响的方法;同时,通过对HFV 模型进行气动分析,认为巡航段中升降舵和鸭翼舵对气动阻力的影响可忽略不计.本文借鉴上述方法,则式(2)中推力、气动力和气动力矩系数可按下式进行计算[29]:

考虑执行机构实际物理特性,控制量需限制在表1 所示容许范围内.需要指出的是,以上鸭翼舵控制量的容许范围能够包络住的范围.

表1 控制量的容许范围Table 1 Admissible ranges for control inputs

1.2 升降舵执行机构模型

HFV 通常在比较恶劣的外部环境下飞行,这会对升降舵的执行机构动态带来一定影响.有时,恶劣的外部飞行环境会激发执行机构中存在的某些未知模态,使得升降舵的实际摆角不能精确地复现期望的理想指令.文献[17]研究了这一问题,并给出了一个考虑未知模态的升降舵执行机构的模型:

式中,u表示指令信号,t表示飞行时间,v表示模态激发后实际产生的指令信号.β(·) 为一函数,表示模态激发的过渡过程:当t＜T1时,β=0;当T1≤t ＜T2时,β=(t-T1)/(T2-T1);当t≥T2时,β=1.p(·) 表示模态实际作用效果,为一未知连续光滑非线性函数.v和p(·) 满足如下假设:

假设1[17].存在有界未知连续光滑函数pˆ(θ,γ,Q,u),使得

除可能存在的未知模态外,升降舵执行机构还存在饱和约束.令分别表示表1 中升降舵摆角上、下边界饱和值,则升降舵实际摆角δe可表示为:

由式(5) 和式(7) 所示的执行机构模型可知,以控制指令信号u作为输入量时,俯仰角速度动态表达式具有非仿射形式.

2 控制器设计

控制器设计目标为:设计φ和u的控制律,使得飞行速度V和航迹倾角γ分别跟踪已知的连续光滑参考轨迹Vref和γref.设计过程中不能使用航迹倾角和攻角的模型输出值构建控制律,同时需要考虑控制量受到饱和特性约束.由于飞行速度V仅受控制量φ的影响.因此,可将HFV 非线性模型分解为由状态V组成的速度子系统和由状态(γ,θ,Q)组成的航迹倾角子系统,然后分别设计φ和u的控制律.

2.1 航迹倾角估计方法

由于不能使用航迹倾角和攻角的模型输出值,因此需要考虑如何构造它们的估计值用于控制律设计.考虑到α=θ-γ,因此只需要估计γ和α其中一个状态即可.对于HFV 纵向通道模型,有以下关系成立:

式中,H表示飞行高度.考虑到H和V可通过机载装置获得相对较准确的测量值,因此,可基于式(8)对γ进行在线估计.令表示γ的估计值,其在线计算方法可设计为:

式中,χV1、χV2为微分器状态,RV ＞0,aV1＞0,aV2＞0为设计参数.定义eV1=χV1-V,eV2=χV2-,由文献[22]的注3 和文献[30]的定理2可知,误差eV1、eV2是有界的,满足:

对于飞行高度H的一阶、二阶导数,引入如下形式的高阶微分器进行在线估计:

式中,χH1、χH2、χH3为微分器状态,RH ＞0,aH1＞0,aH2＞0,aH3＞0为设计参数.定义eH1=χH1-H,同样由文献[22]的注3 和文献[30]的定理2 可知,估计误差eH1、eH2、eH3是有界的,满足:

注1.对于微分器的设计参数,一般来说,R(·)越大,估计误差越小,但选取过大的R(·)可能导致系统状态出现振荡.对此,可以在实际飞行前根据其输入量可能的变化情况,通过仿真选择合适的参数R(·)和a(·).以式(11)所示跟踪微分器为例,其输入量为速度V,则可根据Vref和典型参数偏差组合设置生成若干实际飞行过程中V潜在的变化情况,然后根据这些信息事先检查不同设计参数选择下跟踪微分器的估计效果,从而在实际飞行前确定式(11)中参数的具体取值.

2.2 速度子系统控制器设计

定义速度跟踪误差=V -Vref,根据式(1),的动态可描述为:

由于γ是未知的,因此fV、gV为未知项.为解决这一问题,将式(17)改写为如下形式:

式中

因此,fV、gV可通过机载测量数据计算得到.

考虑到φ需满足表1 所示约束,引入如下辅助系统:

式中,ξV为辅助状态,k1＞0为待设计增益,φc为不考虑饱和约束情况下的控制指令.定义新的误差量,由gV的形式可知gV0,则可设计φc的控制律为:

注2.式(19)和式(20)中增益k1和kV大小的选取会产生不同的控制效果.当k1＞kV时,控制器会更快速地消除饱和引起的跟踪偏差,但控制量φ的饱和时间会变长;当k1＜kV时,控制量φ会更快速地脱离饱和状态,但系统的跟踪误差会变大.比较两种选择的影响,考虑到HFV 飞行中对控制偏差的要求较高,实际中建议选择k1＞kV的设置.

2.3 航迹倾角子系统控制器设计

基于backstepping 方法设计航迹倾角子系统控制器.

步骤1.定义航迹倾角跟踪误差根据式(1),由于Tsinα/(mV) 项的存在,以θ作为虚拟控制量时,γ的动态具有非仿射形式.令,根据文献[29]的气动数据可得.对此,可将Tsinα/(mV) 项视为一类模型干扰项,则的动态可描述为:

式中,gγ ＞0.考虑到γ的模型输出值不能直接使用,将式(22)进一步改写为如下形式:

backstepping 方法设计过程中需要在线计算虚拟控制量的一阶导数.对于,采用以下跟踪微分器进行在线估计:

步骤2.定义俯仰角跟踪误差根据式(1),的动态可描述为:

以Q作为虚拟控制量,则虚拟控制律Qcmd可初步设计为:

式中,kθ＞0为待设计增益.对于,采用以下跟踪微分器进行在线估计:

式中,χQc1和χQc2为微分器状态.RQc ＞0,aQc1＞0,aQc2＞0为设计参数.定义eQc1=χQc1-Qcmd,eQc2=χQc2-cmd,eQc1、eQc2是有界的,满足:

步骤3.定义俯仰角速度跟踪误差=Q-Qcmd.根据式(1),的动态可描述为:

考虑到γ的模型输出值不能直接使用的问题,将式(33)进一步改写为如下形式:

对于式(34)中的控制量δe,受伺服系统物理特性限制,其大小受限于表1 所示饱和边界约束.根据式(7),升降舵实际摆角δe可重新表示为

对于v,根据式(5),其大小由控制指令u决定.限制控制指令u的大小范围为:

式中,uc为待设计的不受限控制指令.在式(37)的限制下,由假设1 可知dδ是有界的;令du=β(t-T)p(θ,γ,Q,u),由假设1可知du是有界的,且∂v/∂u ＞0.则式(34)可表达为:

注3.对于式(37)所设计的控制指令u的限制范围,需要综合考虑舵面控制能力和对干扰项ΔQ的影响.一方面,若u的边界值较小,则有可能出现u已达到边界值而实际控制舵面摆角δe未达到饱和值的情况,这相当于限制了升降舵的控制能力,在实际应用中应尽量避免;另一方面,若u的边界值较大,可能会出现u未达到边界值,而实际控制舵面摆角δe已达到饱和的情况,这会提高dδ可能的最大幅值,进而提高 ΔQ可能的最大幅值.对于以上两点,式(37)所规定的u的范围可以保证摆角δe在任何情况下可以达到其能够达到的最大摆角,同时最大限度地降低 ΔQ可能的最大幅值.

考虑到u需满足式(37)所示约束,引入如下辅助系统:

式中,ξu为辅助状态,k2＞0 为待设计增益.引入辅助状态ξu后,为提高控制精度,将式(30)所示初步设计的Qcmd控制律改进为:

为估计干扰项 ΔQ,首先引入如下跟踪微分器在线估计俯仰角速度的一阶导数:

式中,χQ1、χQ2为微分器状态,RQ ＞0,aQ1＞0,aQ2＞0为设计参数.定义eQ1=χQ1-Q,eQ2=χQ2-,误差eQ1、eQ2是有界的,满足:

式中,kQ ＞0为待设计增益,Q为 ΔQ的估计值,通过下式计算:

注4.对于干扰项 ΔV、Δγ和 ΔQ,现有方法一般设计估计策略对其进行在线估计,但这些方法大都需要使用气流角作为估计器的组成部分[7,17,25].此外,部分研究[7-8,16-17,25]中需要假设干扰项的导数存在,然而,由干扰项的形式可知,当控制量触发饱和时,其导数并不存在.本文设计的式(21)、式(26)和式(44)干扰估计方法能够适应气流角测量不准确以及干扰项不光滑的情况.

综合以上控制策略,图1 给出了本文设计的输出反馈控制器结构图.

图1 控制器结构图Fig.1 Structure diagram of the controller

3 稳定性分析

首先分析式(19)和式(39)中辅助状态ξV和ξu的稳定性,给出如下定理:

定理1.若控制量φ的饱和时间tφ是有界的,则式(19)所示辅助状态ξV是有界的;若控制量u的饱和时间tu是有界的,则式(39)所示辅助状态ξu是有界的.

证明.考虑Lyapunov 函数,对其求导可得:

在有限时间tφ内,gV(φ-φc) 是有界的.设|gV(φ-φc)|≤ϑ(tφ), 则当|ξV|＞ϑ/k1时,有W˙ξV ＜0,这表明状态ξV是有界的.而对于辅助状态ξu,基于同样的分析方式可知,当控制量u出现有限时间饱和时ξu是有界的. □

定理2.通过式(21)、式(26)和式(44)得到的干扰估计值导致的干扰估计误差是有界稳定的.

证明.对于干扰误差,由式(11)和式(18)可得:

由式(16)可知σ2有界,因此误差量是有界稳定的,满足

由式(42)可知eQ2有界,因此误差量是有界稳定的,满足□

定理3.若控制量φ和u的饱和时间是有界的,则由式(1)、式(5)、式(7)、式(11)、式(13)、式(15)、式(19)～(21)、式(24)～(27)、式(31)、式(35)、式(39)～(41)、式(43)和式(44)组成的闭环系统中飞行状态跟踪误差是有界稳定的.

证明.考虑如下Lyapunov 函数:

综合式(51)、式(53)、式(55)和式(57)的结果,可得W的一阶导数为:

由定理2、式(16)、式(28)和式(32)可得:

式中,η为一选取的常数,满足0＜η ＜k.由式(60)可知,对∀‖x‖2≥ε/η,有:

这即证明误差量zV、μγ、和zQ是有界稳定的.而由式(16)已知σ1是有界的;由定理1 已知,当控制量φ和u的饱和时间有界时,ξV、ξu是有界的.因此,误差量zQ-ξu也是有界稳定的. □

4 仿真校验

本节通过对比仿真实验检验所设计方法的有效性.仿真模型基于式(1)～(4)建立,仿真模型中考虑控制舵摆角对CD的影响,即

具体飞行器总体和气动参数采用文献[29]附表中提供的数据.仿真初值设置为:

式(5)中函数p(·) 设置为

函数β(·) 中参数T1和T2分别设置为T1=7 s,T2=8 s, 假设1中参数设置为0.5.参考轨迹Vref和γref设置为:

仿真时,模型气动参数偏差取-20%,模型简化及外部环境干扰造成的影响设置为:

选取文献[17]方法作为对比方法与本文方法进行控制效果对比.两种方法的控制器参数设置见表2.

表2 控制器参数设置Table 2 Parameter settings of controllers

仿真过程中,为模拟航迹倾角测量不准确的实际情况,在航迹倾角的模型输出量上增加一个0.02°的模拟测量偏差.对比实验共进行3 组仿真:情形1 表示采用本文设计方法的仿真;情形2 表示采用对比方法的仿真;情形3 表示采用对比方法,但在航迹倾角的模型输出量上不加入模拟测量偏差的仿真.仿真结果如图2～13 所示:

图2 和图3 分别给出了速度和航迹倾角状态的跟踪对比情况.结合式(65)可以看出,速度轨迹在t=10 s和t=20 s 附近变化较快,航迹倾角轨迹在t=10 s和t=25 s 附近变化较快.

图2 速度跟踪曲线Fig.2 Tracking curves of velocity

图3 航迹倾角跟踪曲线Fig.3 Tracking curves of flight path angle

图4 和图5 分别给出了速度跟踪误差和速度子系统控制量φ的对比曲线.可以看出,当控制量φ未发生饱和时,3 种仿真状态飞行速度的跟踪误差均较小,且差异很小.在t=10 s和t=29 s 附近,受切变形外部干扰的影响,3 种仿真状态的控制量均出现了短时间的饱和情况,进而引起速度跟踪误差有明显增大现象;由于在t=10s 附近参考轨迹变化较快且控制量饱和时间更长,因此t=10 s附近速度跟踪误差较t=29 s 附近更大.放大t=10 s和t=29 s时刻前后的曲线可以看出,3 种仿真情形中控制量φ发生饱和的初始时刻和饱和时长相差不大,而情形1 因控制量饱和引起的速度跟踪误差明显较小;考虑到当控制量φ未发生饱和时3 种仿真状态的控制效果差异很小,说明式(19)所示抗饱和辅助系统在控制量达到饱和时起到了有益作用.

图4 速度跟踪误差曲线Fig.4 Curves of velocity tracking errors

图5 控制量燃料空气混合比曲线Fig.5 Curves of the control variable fuel-to-air ratio

图6 航迹倾角跟踪误差曲线Fig.6 Curves of flight path angle tracking errors

图7 俯仰角曲线Fig.7 Curves of pitch angle

图8 俯仰角速度曲线Fig.8 Curves of pitch rate

图6～8 分别给出了航迹倾角跟踪误差、俯仰角和俯仰角速度的对比曲线,可以看出各状态均保持稳定;图9 给出了情形1 中航迹倾角估计误差σ1的曲线.由图6 可以看出,情形1 和情形3 的跟踪效果均较好,当航迹倾角变化平稳时两者的跟踪误差均较小.两者跟踪误差较明显时段均出现在航迹倾角变化较大的时段内,但两者跟踪误差的变化趋势存在明显不同.这是因为影响情形1 航迹倾角跟踪误差的主要因素是航迹倾角估计误差,由图3 和图9 可以看出,当航迹倾角变化较大时,航迹倾角的估计误差较大,进而导致图6 中航迹倾角的跟踪误差较大.而对于情形3 中采用的对比方法,引言中已述及,其在俯仰角速度变化较快的情况下控制效果较差,这是影响情形3 航迹倾角跟踪误差的主要因素.相较情形1 和情形3,图6 中情形2 的结果明显较差.由图6 可以看出,模拟的+0.02°航迹倾角测量误差贯穿全程,这是由于对比方法的控制方案中各环节使用的都是含测量误差的航迹倾角测量值,实际控制效果降低的也是含测量误差的航迹倾角与参考值之间的误差.比较图6 中的结果可知,气流角测量误差会严重影响控制品质,必须在控制器设计过程中予以考虑;而本文方法设计的控制器方案能够更好地适应工程实际应用的需要,取得更好的控制效果.

图9 航迹倾角估计误差曲线Fig.9 Curves of the estimated error of flight path angle

图10 给出了控制量升降舵偏角δe的结果曲线.可以看出,情形2 中初始阶段控制量存在明显的振荡,这是由于情形2 考虑了测量误差的影响,存在0.02°航迹倾角初始偏差所致.对于3 种仿真状态,在t=10 s附近,受切变形外部干扰的影响,控制量均出现了短暂的饱和情况;在t=29 s 附近,同样因切变形外部干扰的影响,控制量均出现了短暂的超调现象.综合图5 和图10 的结果可知,本文设计的方法能够在控制量出现有限时间饱和情况时保证飞行稳定性.

图10 控制量升降舵摆角曲线Fig.10 Curves of the control variable deflection of elevator

图11 给出了使用本文方法(情形1)时航迹倾角子系统各控制指令u、v和δe的结果曲线.可以看出,在t=7 s时刻前三者是相同的;在t=7 s 时刻激发执行机构未知模态后,受函数p(·)的影响,u和v、δe之间出现了一定的差距;在执行机构饱和情况发生时,δe=15°, 说明指令u没有限制升降舵的控制能力.在t=9.8 s 时刻附近控制量存在一个小的波动,对照图5 可知,是由于该时刻时控制量φ达到饱和,导致模型参数fQ、ΔQ变化不光滑,进而引起控制量δe出现超调所致.

图11 本文方法航迹倾角子系统控制指令曲线Fig.11 Curves of flight path angle subsystem control commands under proposed method

图12 给出了使用本文方法时系统干扰项ΔV、Δγ、ΔQ的估计曲线.可以看出,本文设计的式(21)、式(26)和式(44)的干扰估计方法,能够较好地跟踪干扰项的实际值,跟踪误差较小.

图12 本文方法总干扰项估计曲线Fig.12 Estimated curves of the lumped disturbances under proposed method

图13 给出了使用本文方法时抗饱和补偿信号ξV和ξu的结果曲线.可以看出,当控制量出现饱和情况时,ξV和ξu对控制量进行了补偿;又考虑到增益k1＞kV,k2＞kQ,由式(20)和式(43)可知,补偿效果为增大控制量的幅值以更快减小跟踪误差.控制量脱离饱和状态后,补偿信号能快速收敛至0附近,保持稳定.

图13 本文方法抗饱和补偿信号曲线Fig.13 Curves of anti-windup compensated signals under proposed method

5 结束语

本文针对非仿射HFV 轨迹跟踪控制问题,面向工程实际需求,提出了一种无需气流角测量值和考虑执行机构饱和特性的输出反馈控制方法.通过对比仿真实验验证了所设计方法在航迹倾角测量精度较差和控制量发生饱和的情况下仍具有良好的控制性能.与现有的针对非仿射HFV 轨迹跟踪控制问题的控制方案相比,所设计方法解决了工程实际中气流角测量精度差以及执行机构存在不确定非线性动态和饱和特性的问题,具有更好的工程应用前景.