CN110231774A - 干扰观测变进气道高超声速飞行器模糊协调控制方法 - Google Patents
干扰观测变进气道高超声速飞行器模糊协调控制方法 Download PDFInfo
- Publication number
- CN110231774A CN110231774A CN201910276521.2A CN201910276521A CN110231774A CN 110231774 A CN110231774 A CN 110231774A CN 201910276521 A CN201910276521 A CN 201910276521A CN 110231774 A CN110231774 A CN 110231774A
- Authority
- CN
- China
- Prior art keywords
- designed
- disturbance observer
- air intake
- design
- intake duct
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Feedback Control In General (AREA)
Abstract
本发明涉及高超声速飞行器控制,为提出用于存在外界干扰的变进气道高超声速飞行器模糊***的协调控制方法,本发明,干扰观测变进气道高超声速飞行器模糊协调控制方法,步骤如下:第一部分,变进气道高超飞行器的严反馈形式的转化设计;第二部分,速度控制器‑非线性干扰观测器的设计;第三部分,高度控制器‑非线性干扰观测器的设计,使用非线性干扰观测器对所设计的子***控制器进行补偿,实时观测干扰,增强所设计的子***控制器的鲁棒性。本发明主要应用于高超声速飞行器控制场合。
Description
技术领域
本发明涉及高超声速飞行器控制***设计领域发明,尤其涉及变进气道高超声速飞行器的协调控制研究领域。具体讲,涉及带有干扰观测器的变进气道高超声速飞行器的模糊协调控制方法。具体涉及干扰观测变进气道高超声速飞行器模糊协调控制方法。
背景技术
高超声速飞行器是指飞行速度大于5马赫、在近空间飞行,集轨道战斗机、空间飞行器,甚至卫星空间站等各方面优点于一身的飞行器,具有大空域、超高速、长距离的特点,能够实现全球快速精确的战略打击,是确保国家战略安全的重要保障之一,具有飞行速度极高、运行周期长、运行性能稳定、生存能力强等特点,已经被广泛应用于军事战略、信息侦察、货物的快速运输等领域。超声速飞行器的技术研究有很强的前瞻性、战略性与带动性。高超声速技术的发展将会对未来军事发展战略、空间技术、武器体系构建乃至整个科学技术进步产生重大影响。
然而,高超声速飞行器采用超燃冲压发动机,其特有的机体/发动机一体化结构在减轻机体负担的同时,也增大了模型的非线性程度,因此,设计合适的控制器对高超声速飞行器***进行控制比较困难。另外,为了保证发动机在高马赫下发挥最佳性能,进气道的设计普遍采用定几何进气道,其结构设计只能够在其设计点处保证飞行器的运行性能,当速度偏离设计点较大时,控制性能急剧下降,特别在低马赫下,激波角增大,发动机进气罩不能捕获充足的气流,限制了发动机运行。变几何进气道的结构设计突破了这一局限,但是,关于变几何进气道在控制领域的研究仍处于起步阶段,存在的控制难点体现在以下几个方面:(1)引入变几何进气道结构后会引起气动力的不确定性,尤其对推力产生的影响较大;(2)变几何进气道结构的引入会增大飞行器控制模型的非线性程度;(3)变几何进气道控制与飞行控制的耦合影响,使得控制***的复杂度大大增加;(4)在近空间的复杂环境中飞行控制时,外界干扰变化剧烈,导致飞行器***的不确定性增强,飞行控制性能难以达到期望的控制效果。
因此,针对上述局限,本发明提出一种基于干扰观测器的变几何进气道高超声速飞行器的模糊协调控制方法,基于该算法设计的控制器考虑了变进气道控制与飞行控制的耦合影响,突破了定进气道结构下飞行控制不能满足机动性飞行的局限,实现飞行器的稳定条件下机动性的提高;此外,基于该算法设计的模糊2型控制***,能够对由可移动进气罩引入带来的参数不确定项进行逼近,简化了控制律,设计非线性干扰观测器对所设计的自适应模糊2型控制器进行补偿,实时的观测***的干扰,显著提高了变几何进气道高超声速飞行器在近空间复杂的飞行环境中稳定飞行的能力,具有很强的鲁棒性,提高了对于***不确定性和干扰的观测能力,有效增强***的控制性能。
发明内容
为克服现有技术的不足,本发明旨在提出一种用于存在外界干扰的变进气道高超声速飞行器模糊2型***的协调控制方法。为此,本发明采取的技术方案是,干扰观测变进气道高超声速飞行器模糊协调控制方法,步骤如下:
第一部分,变进气道高超飞行器的严反馈形式的转化设计:将可移动唇罩式变几何进气道飞行器的纵向模型转化为考虑干扰的严反馈形式,方便后续使用反步法进行控制器的设计;
第二部分,速度控制器-非线性干扰观测器的设计:基于非线性干扰观测器,采用模糊2型逻辑理论对速度控制器进行设计,其中,使用模糊2型逻辑***对控制律中的不确定复杂项进行在线逼近,达到简化控制律的目的,并使用非线性干扰观测器对所设计的自适应模糊2型速度控制器进行补偿,达到实时观测***干扰的目的;
第三部分,高度控制器-非线性干扰观测器的设计:使用反步法将高度控制***拆分为航迹角、攻角、俯仰角3个子***,选取每个子***中的状态量以及虚拟控制量,并在采用模糊2型逻辑***设计每个高度子***的控制律的同时,使用非线性干扰观测器对所设计的子***控制器进行补偿,实时观测干扰,增强所设计的子***控制器的鲁棒性。
具体实现过程如下:
第一步,考虑干扰的严反馈形式的转化:基于Bolender和Doman所建立的纵向刚性模型,并使用对变几何进气道飞行器进行受力分析得到的气动力的解析表达式,将可移动唇罩式变几何进气道飞行器的纵向模型转化为考虑干扰的严反馈形式,变几何进气道飞行器的纵向模型为:
式中:m、Iyy、g分别为质量、转动惯量和重力加速度,升力L、阻力D、推力T、以及俯仰力矩M的具体解析表达式如下:
式中:表示动压,ρ、S、分别表示大气密度、参考面积、平均气动弦长,zT表示推力系数:zT=Lftanτ1l+1/2hi,CL、CD、CT、CM表示在不考虑变几何进气道的情况下,飞行器的气动力系数;表示引入可移动唇罩对飞行器气动力造成的影响,其大小与唇罩伸长量有关;
转化的严反馈形式如下:
式中:
其中,是由于可移动唇罩的使用引入的非线性不确定项;干扰项Di包括:***中由于参数拟合造成的不确定性以及飞行器受到的外来干扰dii=V,γ,α,Q,表示如下:
第二步,基于非线性干扰观测器的自适应速度控制器设计:定义速度误差面SV=V-Vd,使用自适应模糊2型***对fV进行在线逼近:
定义速度子***干扰DV=ΩV+εV,则有:
设计如下非线性干扰观测器对DV进行在线观测:
式中:和分别是和的估计值,且定义为模糊2型***的估计误差;为干扰DV的估计值,且有zV为非线性干扰观测器的中间状态量;βV为待设计正数,则有:
则飞行器速度子***的控制量设计如下:
式中:kV0为待设计正数。
第三步,基于非线性干扰观测器的航迹角子***的控制器设计:首先,把飞行器的高度指令hd转化为航迹角指令γd:
式中:kh0为待设计正数,定义航迹角误差面Sγ=γ-γd,使用自适应模糊2型***对fγ进行在线逼近:
定义航迹角子***干扰Dγ=Ωγ+εγ,则有:
设计如下非线性干扰观测器对Dγ进行在线观测:
式中:和分别是和的估计值,且定义为模糊2型***的估计误差;为干扰Dγ的估计值,且有zγ为非线性干扰观测器的中间状态量;βγ为待设计正数,则有:
则虚拟控制输入可设计为:
式中:kγ0是待设计的正参数;
第四步,基于非线性干扰观测器的攻角子***的控制器设计:定义攻角误差面Sα=α-αd,使用自适应模糊2型***对fα进行在线逼近:
定义攻角子***干扰Dα=Ωα+εα,则有:
由于Dα未知,会对跟踪精度造成影响,因此设计如下非线性干扰观测器对Dα进行在线观测:
式中:和分别是和的估计值,且定义为模糊2型***的估计误差;为干扰Dα的估计值,且有zα为非线性干扰观测器的中间状态量;βα为待设计正数,则有:
则虚拟控制输入可设计为:
式中:kα0是待设计的正参数;
第五步,基于非线性干扰观测器的俯仰角速率子***的控制器设计:定义俯仰角速率误差面SQ=Q-Qd,使用自适应模糊2型***对fQ进行在线逼近:
定义俯仰角速率子***干扰DQ=ΩQ+εQ,则有:
设计如下非线性干扰观测器对DQ进行在线观测:
式中:和分别是和的估计值,且定义为模糊2型***的估计误差;为干扰DQ的估计值,且有zQ为非线性干扰观测器的中间状态量;βQ为待设计正数,则有:
则飞行器高度子***的控制量设计如下:
式中:kQ0为待设计正数。
本发明的特点及有益效果是:
基于该算法设计的控制器考虑了变进气道控制与飞行控制的耦合影响,突破了定进气道结构下飞行控制不能满足机动性飞行的局限,实现飞行器的稳定条件下机动性的提高;此外,基于该算法设计的模糊2型控制***,能够对由可移动进气罩引入带来的参数不确定项进行逼近,简化了控制律,设计非线性干扰观测器对所设计的自适应模糊2型控制器进行补偿,实时的观测***的干扰,显著提高了变几何进气道高超声速飞行器在近空间复杂的飞行环境中稳定飞行的能力,具有很强的鲁棒性,提高了对于***不确定性和干扰的观测能力,有效增强***的控制性能。
附图说明:
附图1基于干扰观测器的变进气道高超声速飞行器的模糊协调控制***结构图。
附图2速度及速度误差变化曲线(有无干扰观测器)。
附图3高度及高度误差变化曲线(有无干扰观测器)。
附图4燃料当量比变化曲线(有无干扰观测器)。
附图5升降舵偏角变化曲线(有无干扰观测器)。
附图6速度及速度误差变化曲线(有无可移动进气罩)。
附图7高度及高度误差变化曲线(有无可移动进气罩)。
附图8燃料当量比变化曲线(有无可移动进气罩)。
附图9升降舵偏角变化曲线(有无可移动进气罩)。
具体实施方式
本发明的目的在于提出一种用于存在外界干扰的变进气道高超声速飞行器模糊2型***的协调控制方法。一方面,变进气道高超声速飞行器中可移动唇罩的引入在提高飞行器发动机的机动性的同时,会使飞行器的气动力值发生改变,增加了变进气道高超声速飞行器的控制***设计的复杂性和难度。另一方面,飞行器在近空间飞行时,外界干扰变化剧烈,导致飞行器***的不确定性增强。因此,针对存在外来干扰时,变进气道高超声速飞行器***的不确定性和复杂性,本发明提出了一种基于干扰观测器的变进气道高超声速飞行器的模糊协调控制方法,首先,将可移动唇罩式变几何进气道的纵向模型转化为严反馈形式,其可由速度V、高度h、航迹角γ、攻角α、俯仰角速率Q来表示;其次,利用反步法分别对速度控制器和高度控制器进行设计;然后,采用区间模糊2型逻辑***对控制律中的不确定复杂性进行在线逼近;最后设计非线性干扰观测器对所设计的自适应模糊2型控制器进行补偿,实时的观测***的干扰。由此可知,本发明提出的基于干扰观测器的变进气道飞行器的模糊协调控制方法是一种实时控制和观测干扰的方法,可以实现对于***不确定性和干扰的观测能力,其控制精度高、响应速度快;该方法能够解决变几何进气道高超声速飞行器在近空间飞行中的不确定性和控制问题,具有很强的鲁棒性;该控制算法考虑变进气道控制与飞行控制的耦合影响,可以实现飞行器的稳定条件下机动性的提高,突破定进气道结构下飞行控制不能满足机动性飞行的不足,解决在外界干扰影响下变进气道高超声速飞行器宽速域、灵活机动的稳定飞行的控制问题,达到提高变进气道高超飞行器的控制性能,增强***鲁棒性和可靠性的目的。
本发明以控制理论方法和虚拟仿真技术相结合为主要研究手段,发明一种基于干扰观测器的变进气道飞行器的模糊协调控制方法,通过matlab软件进行虚拟仿真实验,验证了本方法的有效性。步骤如下:
第一部分,变进气道高超飞行器的严反馈形式的转化设计:考虑变进气道中可移动唇罩的引入造成的影响之后,变几何进气道高超声速飞行器的模型变得非常复杂,不适用于线性化的处理方法,将可移动唇罩式变几何进气道飞行器的纵向模型转化为考虑干扰的严反馈形式,方便后续使用反步法进行控制器的设计。
第二部分,速度控制器-非线性干扰观测器的设计:基于非线性干扰观测器,采用模糊2型逻辑理论对速度控制器进行设计,其中,使用模糊2型逻辑***对控制律中的不确定复杂项进行在线逼近,达到简化控制律的目的,并使用非线性干扰观测器对所设计的自适应模糊2型速度控制器进行补偿,达到实时观测***干扰的目的。
第三部分,高度控制器-非线性干扰观测器的设计:使用反步法将高度控制***拆分为航迹角、攻角、俯仰角3个子***,选取每个子***中的状态量以及虚拟控制量,并在采用模糊2型逻辑***设计每个高度子***的控制律的同时,使用非线性干扰观测器对所设计的子***控制器进行补偿,实时观测干扰,增强所设计的子***控制器的鲁棒性。
具体实现过程如下:
第一步,考虑干扰的严反馈形式的转化:基于Bolender和Doman所建立的纵向刚性模型,并使用对变几何进气道飞行器进行受力分析得到的气动力(升力L、阻力D、推力T、以及俯仰力矩M)的解析表达式,将可移动唇罩式变几何进气道飞行器的纵向模型转化为考虑干扰的严反馈形式:
变几何进气道飞行器的纵向模型为:
式中:m、Iyy、g分别为质量、转动惯量和重力加速度,升力L、阻力D、推力T、以及俯仰力矩M的具体解析表达式如下:
式中:表示动压,ρ、S、分别表示大气密度、参考面积、平均气动弦长。zT表示推力系数:zT=Lftanτ1l+1/2hi。值得注意到是,CL、CD、CT、CM表示在不考虑变几何进气道的情况下,飞行器的气动力系数;表示引入可移动唇罩对飞行器气动力造成的影响,其大小与唇罩伸长量有关。
转化的严反馈形式如下:
式中:
其中,是由于可移动唇罩的使用引入的非线性不确定项;干扰项Di(i=V,γ,α,Q)包括:***中由于参数拟合造成的不确定性以及飞行器受到的外来干扰di。表示如下:
第二步,基于非线性干扰观测器的自适应速度控制器设计:定义速度误差面SV=V-Vd,使用自适应模糊2型***对fV进行在线逼近:
定义速度子***干扰DV=ΩV+εV,则有:
设计如下非线性干扰观测器对DV进行在线观测:
式中:和分别是和的估计值,且定义为模糊2型***的估计误差;为干扰DV的估计值,且有zV为非线性干扰观测器的中间状态量;βV为待设计正数。则有:
则飞行器速度子***的控制量设计如下:
式中:kV0为待设计正数。
第三步,基于非线性干扰观测器的航迹角子***的控制器设计:首先,把飞行器的高度指令hd转化为航迹角指令γd:
式中:kh0为待设计正数。
定义航迹角误差面Sγ=γ-γd,使用自适应模糊2型***对fγ进行在线逼近:
定义航迹角子***干扰Dγ=Ωγ+εγ,则有:
设计如下非线性干扰观测器对Dγ进行在线观测:
式中:和分别是和的估计值,且定义为模糊2型***的估计误差;为干扰Dγ的估计值,且有zγ为非线性干扰观测器的中间状态量;βγ为待设计正数。则有:
则虚拟控制输入可设计为:
式中:kγ0是待设计的正参数。
第四步,基于非线性干扰观测器的攻角子***的控制器设计:定义攻角误差面Sα=α-αd,使用自适应模糊2型***对fα进行在线逼近:
定义攻角子***干扰Dα=Ωα+εα,则有:
由于Dα未知,会对跟踪精度造成影响,因此设计如下非线性干扰观测器对Dα进行在线观测:
式中:和分别是和的估计值,且定义为模糊2型***的估计误差;为干扰Dα的估计值,且有zα为非线性干扰观测器的中间状态量;βα为待设计正数。则有:
则虚拟控制输入可设计为:
式中:kα0是待设计的正参数。
第五步,基于非线性干扰观测器的俯仰角速率子***的控制器设计:定义俯仰角速率误差面SQ=Q-Qd,使用自适应模糊2型***对fQ进行在线逼近:
定义俯仰角速率子***干扰DQ=ΩQ+εQ,则有:
设计如下非线性干扰观测器对DQ进行在线观测:
式中:和分别是和的估计值,且定义为模糊2型***的估计误差;为干扰DQ的估计值,且有zQ为非线性干扰观测器的中间状态量;βQ为待设计正数。则有:
则飞行器高度子***的控制量设计如下:
式中:kQ0为待设计正数。
为了验证本发明提出的基于干扰观测器的变进气道高超声速飞行器的模糊协调控制方法的有效性,设计了变进气道高超声速飞行器控制综合虚拟仿真***,并在该***上进行了仿真实验。在飞行器虚拟仿真控制环境下,仿真参数设定如下:
1)变几何进气道高超飞行器物理参数:飞行器质量m=191slug,机体转动惯量Iyy=86723slug·ft2,重力加速度g=31.97ft/s2,大气密度ρ与飞行高度有关,参考面积S=100ft2,平均几何弦长
2)变几何进气道气动解析式中的气动系数:如表1-1:
表1-1进气伸长量以及气动力拟合系数
3)控制器的参数设计:如表1-2所示:
表1-2控制参数
分两种情况进行仿真测试验证,第一种情况将本发明提出的算法与未设计非线性干扰观测器的***的控制性能进行对比,如图2-图5所示;第二种情况将本发明提出的算法与定几何进气道高超飞行器***的控制性能进行对比,此种情况的***的控制性能的仿真结果如图6-图9所示。
第一种情况***的给定指令为:速度从7Ma变化到7.5Ma、高度从75000ft变化76000ft。控制***所受到的外部扰动为:在110<t≤150s时,dV=2ft/s、dγ=0.01deg、dα=0.02deg、dQ=0.3deg,在160<t≤250s时,dV=2sin(0.2t)ft/s、dγ=0.01sin(0.2t)deg、dα=0.02sin(0.2t)deg、dQ=0.3sin(0.3t)deg/s。拟合参数的参数摄动设定为-15%。
从图2-图5中,可以看到在未加入干扰的时候,两个控制***均能够实现对于给定指令的快速准确跟踪。相比之下,基于干扰观测器的控制***跟踪误差更小,响应更快,这是因为所设计非线性干扰观测器对***中的参数摄动更为敏感,从而对控制器进行补偿。在外来干扰作用阶段,很明显有干扰观测器的***抵御干扰的能力更强。从图2、图3中能够发现,有干扰观测器的***能够保证更高的速度和高度响应精度,而没有使用干扰观测器的控制***在干扰作用时,虽然也能够保证***的稳定性,但跟踪精度较差。从图4、图5中可以发现,基于干扰观测器的设计的执行机构的动作变化更小,这也减小了在面临干扰时,执行机构频繁动作造成的损耗。
第二种情况***的给定指令为:速度从7Ma变化到9.1Ma(速度指令变化较快)、高度从75000ft变化76000ft。在加速阶段(20<t≤100s),给***加入外部扰动:dV=2sin(0.2t)ft/s、dγ=0.01sin(0.2t)deg、dα=0.02deg、dQ=0.3deg。拟合参数的参数摄动设定为-20%。
从图6-图9中,可以看到,在加速阶段,定几何进气道飞行器由于执行机构饱和,***的运行性能较差,在外界干扰的影响下,直接导致控制失败。而变几何进气道飞行器在加速阶段控制性能较好,能够抵抗外界干扰,保证控制***的稳定性,完成对给定指令的跟踪。因而,使用变几何进气道也能够在一定程度上增强飞行器控制***的鲁棒性。
本发明提出的基于多变量干扰补偿的四旋翼无人机轨迹与姿态协同控制算法总体技术方案如图1所示,整个***主要包括三部分:严反馈形式的转化设计、速度控制器-非线性干扰观测器设计和高度控制器-非线性干扰观测器的设计,具体技术方案如下:
第一部分,变进气道高超飞行器的严反馈形式的转化设计:考虑变进气道中可移动唇罩的引入造成的影响之后,变几何进气道高超声速飞行器的模型变得非常复杂,不适用于线性化的处理方法,将可移动唇罩式变几何进气道飞行器的纵向模型转化为考虑干扰的严反馈形式,方便后续使用反步法进行控制器的设计。
第二部分,速度控制器-非线性干扰观测器的设计:基于非线性干扰观测器,采用模糊2型逻辑理论对速度控制器进行设计,其中,使用模糊2型逻辑***对控制律中的不确定复杂项进行在线逼近,达到简化控制律的目的,并使用非线性干扰观测器对所设计的自适应模糊2型速度控制器进行补偿,达到实时观测***干扰的目的。
第三部分,高度控制器-非线性干扰观测器的设计:使用反步法将高度控制***拆分为航迹角、攻角、俯仰角3个子***,选取每个子***中的状态量以及虚拟控制量,其中引入低阶滤波器,解决反步法中出现的“微分***”的问题,并在采用模糊2型逻辑***设计每个高度子***的控制律的同时,使用非线性干扰观测器对所设计的子***控制器进行补偿,实时观测干扰,增强所设计的子***控制器的鲁棒性。
最后为了验证本发明提出算法的有效性,搭建变几何进气道高超飞行器协同控制的MATLAB仿真***,并与情况一(有无干扰观测器)和情况二(有无可移动进气罩)的控制仿真结果进行对比,验证本发明提出的算法鲁棒性更强,对给定指令的跟踪效果更佳。
第一步,考虑干扰的严反馈形式的转化
基于Bolender和Doman所建立的纵向刚性模型,并使用对变几何进气道飞行器进行受力分析得到的气动力(升力L、阻力D、推力T、以及俯仰力矩M)的解析表达式,将可移动唇罩式变几何进气道飞行器的纵向模型转化为考虑干扰的严反馈形式:
变几何进气道飞行器的纵向模型为:
式中:m、Iyy、g分别为质量、转动惯量和重力加速度,升力L、阻力D、推力T、以及俯仰力矩M的具体解析表达式如下:
式中:表示动压,ρ、S、分别表示大气密度、参考面积、平均气动弦长。zT表示推力系数:zT=Lftanτ1l+1/2hi。值得注意到是,CL、CD、CT、CM表示在不考虑变几何进气道的情况下,飞行器的气动力系数;表示引入唇罩对飞行器气动力造成的影响,其大小与唇罩伸长量有关。
转化的严反馈形式如下:
式中:
其中,是由于可移动唇罩的使用引入的非线性不确定项。干扰项Di(i=V,γ,α,Q)包括:***中由于参数拟合造成的不确定性以及飞行器受到的外来干扰di。表示如下:
第二步,基于非线性干扰观测器的自适应速度控制器设计
定义速度误差面SV=V-Vd,使用自适应模糊2型***对fV进行在线逼近:
为了保证控制***的跟踪精度,将自适应模糊2型***的逼近误差εV视为***干扰的一部分,定义速度子***干扰DV=ΩV+εV,则有:
由于DV未知,会对跟踪精度造成影响,因此设计如下非线性干扰观测器对DV进行在线观测:
式中:和分别是和的估计值,且定义为模糊2型***的估计误差;为干扰DV的估计值,且有zV为非线性干扰观测器的中间状态量;βV为待设计正数。则有:
则飞行器速度子***的控制量设计如下:
式中:kV0为待设计正数。和的自适应更新律设计如下:
式中:μVl、χVl、μVr和χVr是待设计的正参数。
第三步,基于非线性干扰观测器的航迹角子***的控制器设计
首先,把飞行器的高度指令hd转化为航迹角指令γd:
式中:kh0为待设计正数。
定义航迹角误差面Sγ=γ-γd,使用自适应模糊2型***对fγ进行在线逼近:
为了保证控制***的跟踪精度,将自适应模糊2型***的逼近误差εγ视为***干扰的一部分,定义航迹角子***干扰Dγ=Ωγ+εγ,则有:
由于Dγ未知,会对跟踪精度造成影响,因此设计如下非线性干扰观测器对Dγ进行在线观测:
式中:和分别是和的估计值,且定义为模糊2型***的估计误差;为干扰Dγ的估计值,且有zγ为非线性干扰观测器的中间状态量;βγ为待设计正数。则有:
则虚拟控制输入可设计为:
式中:kγ0是待设计的正参数。将输入一阶滤波器,可得到控制量αd及其导数
定义滤波器的滤波误差为则有:
式中:且|Bα(·)|≤Mα,Mα是一个大于0的常数。和的自适应更新律设计如下:
式中:μγl、χγl、μγr和χγr是待设计的正参数。
第四步,基于非线性干扰观测器的攻角子***的控制器设计
定义攻角误差面Sα=α-αd,使用自适应模糊2型***对fα进行在线逼近:
为了保证控制***的跟踪精度,将自适应模糊2型***的逼近误差εα视为***干扰的一部分,定义攻角子***干扰Dα=Ωα+εα,则有:
由于Dα未知,会对跟踪精度造成影响,因此设计如下非线性干扰观测器对Dα进行在线观测:
式中:和分别是和的估计值,且定义为模糊2型***的估计误差;为干扰Dα的估计值,且有zα为非线性干扰观测器的中间状态量;βα为待设计正数。则有:
则虚拟控制输入可设计为:
式中:kα0是待设计的正参数。将输入一阶滤波器,可得到控制量Qd及其导数
定义滤波器的滤波误差为则有:
式中:且|BQ(·)|≤MQ,MQ是一个大于0的常数。和的自适应更新律设计如下:
式中:μγl、χγl、μγr和χγr是待设计的正参数。
第五步,基于非线性干扰观测器的俯仰角速率子***的控制器设计
定义俯仰角速率误差面SQ=Q-Qd,使用自适应模糊2型***对fQ进行在线逼近:
为了保证控制***的跟踪精度,将自适应模糊2型***的逼近误差εQ视为***干扰的一部分,定义俯仰角速率子***干扰DQ=ΩQ+εQ,则有:
由于DQ未知,会对跟踪精度造成影响,因此设计如下非线性干扰观测器对DQ进行在线观测:
式中:和分别是和的估计值,且定义为模糊2型***的估计误差;为干扰DQ的估计值,且有zQ为非线性干扰观测器的中间状态量;βQ为待设计正数。则有:
则飞行器高度子***的控制量设计如下:
式中:kQ0为待设计正数。和的自适应更新律设计如下:
式中:μQl、χQl、μQr和χQr是待设计的正参数。
基于以上5步,就完成了整个控制***设计过程,为了验证有效性具体步骤,首先将变几何进气道高超飞行器的协调控制***在Matlab进行设计,并进行了仿真实验,仿真过程如下:
(1)参数设置
1)变几何进气道高超飞行器物理参数:飞行器质量m=191slug,机体转动惯量Iyy=86723slug·ft2,重力加速度g=31.97ft/s2,大气密度ρ与飞行高度有关,参考面积S=100ft2,平均几何弦长
2)变几何进气道气动解析式中的气动系数:如表1-1:
表1-1进气伸长量以及气动力拟合系数
3)控制器的参数设计:如表1-2所示:
表1-2控制参数
分两种情况进行仿真测试验证,在第一种情况的仿真验证过程中,速度从7Ma变化到7.5Ma、高度从75000ft变化76000ft。控制***所受到的外部扰动为:在110<t≤150s时,dV=2ft/s、dγ=0.01deg、dα=0.02deg、dQ=0.3deg,在160<t≤250s时,dV=2sin(0.2t)ft/s、dγ=0.01sin(0.2t)deg、dα=0.02sin(0.2t)deg、dQ=0.3sin(0.3t)deg/s。拟合参数的参数摄动设定为-15%。第二种情况的仿真验证过程中,***的给定指令为:速度从7Ma变化到9.1Ma(速度指令变化较快)、高度从75000ft变化76000ft。在加速阶段(20<t≤100s),给***加入外部扰动:dV=2sin(0.2t)ft/s、dγ=0.01sin(0.2t)deg、dα=0.02deg、dQ=0.3deg。拟合参数的参数摄动设定为-20%。
(2)结果分析
在上述给定的条件下,将本发明提出的算法分两种情况进行对比:
第一种情况将本发明提出的算法与未设计非线性干扰观测器的***的控制性能进行对比,此种情况的***的控制性能的仿真结果如图2-图5所示。
第二种情况将本发明提出的算法与定几何进气道高超飞行器的***的控制性能进行对比,此种情况的***的控制性能的仿真结果如图6-图9所示。
情况一(有无干扰观测器)控制结果分析:从图2-图5中,可以看到在未加入干扰的时候,两个控制***均能够实现对于给定指令的快速准确跟踪。相比之下,基于干扰观测器的控制***跟踪误差更小,响应更快,这是因为所设计非线性干扰观测器对***中的参数摄动更为敏感,从而对控制器进行补偿。在外来干扰作用阶段,很明显有干扰观测器的***抵御干扰的能力更强。从图2、图3中能够发现,有干扰观测器的***能够保证更高的速度和高度响应精度,而没有使用干扰观测器的控制***在干扰作用时,虽然也能够保证***的稳定性,但跟踪精度较差。从图4、图5中可以发现,基于干扰观测器的设计的执行机构的动作变化更小,这也减小了在面临干扰时,执行机构频繁动作造成的损耗。
情况二(有无可移动进气罩)控制结果分析:从图6-图9中,可以看到,在加速阶段,定几何进气道飞行器由于执行机构饱和,***的运行性能较差,在外界干扰的影响下,直接导致控制失败。而变几何进气道飞行器在加速阶段控制性能较好,能够抵抗外界干扰,保证控制***的稳定性,完成对给定指令的跟踪。因而,使用变几何进气道也能够在一定程度上增强飞行器控制***的鲁棒性。
因此,本发明提出的基于干扰观测器的变进气道高超声速飞行器的模糊协调控制方法实用性更好、鲁棒性更强,也更适用于快速时变性、复杂性的可移动唇罩式变几何进气道高超声速飞行器控制***。并且,在考虑***干扰的情况下,变几何进气道飞行器控制***的鲁棒性更强,对给定指令的跟踪效果更佳。
Claims (2)
1.一种干扰观测变进气道高超声速飞行器模糊协调控制方法,其特征是,步骤如下:
第一部分,变进气道高超飞行器的严反馈形式的转化设计:将可移动唇罩式变几何进气道飞行器的纵向模型转化为考虑干扰的严反馈形式,方便后续使用反步法进行控制器的设计;
第二部分,速度控制器-非线性干扰观测器的设计:基于非线性干扰观测器,采用模糊2型逻辑理论对速度控制器进行设计,其中,使用模糊2型逻辑***对控制律中的不确定复杂项进行在线逼近,达到简化控制律的目的,并使用非线性干扰观测器对所设计的自适应模糊2型速度控制器进行补偿,达到实时观测***干扰的目的;
第三部分,高度控制器-非线性干扰观测器的设计:使用反步法将高度控制***拆分为航迹角、攻角、俯仰角3个子***,选取每个子***中的状态量以及虚拟控制量,并在采用模糊2型逻辑***设计每个高度子***的控制律的同时,使用非线性干扰观测器对所设计的子***控制器进行补偿,实时观测干扰,增强所设计的子***控制器的鲁棒性。
2.如权利要求1所述的干扰观测变进气道高超声速飞行器模糊协调控制方法,其特征是,具体实现过程如下:
第一步,考虑干扰的严反馈形式的转化:基于Bolender和Doman所建立的纵向刚性模型,并使用对变几何进气道飞行器进行受力分析得到的气动力的解析表达式,将可移动唇罩式变几何进气道飞行器的纵向模型转化为考虑干扰的严反馈形式,变几何进气道飞行器的纵向模型为:
式中:m、Iyy、g分别为质量、转动惯量和重力加速度,升力L、阻力D、推力T、以及俯仰力矩M的具体解析表达式如下:
式中:表示动压,ρ、S、分别表示大气密度、参考面积、平均气动弦长,zT表示推力系数:zT=Lftanτ1l+1/2hi,CL、CD、CT、CM表示在不考虑变几何进气道的情况下,飞行器的气动力系数;表示引入可移动唇罩对飞行器气动力造成的影响,其大小与唇罩伸长量有关;
转化的严反馈形式如下:
式中:
gα=1;
其中,是由于可移动唇罩的使用引入的非线性不确定项;干扰项Di包括:***中由于参数拟合造成的不确定性△fi,△gi,以及飞行器受到的外来干扰dii=V,γ,α,Q,表示如下:
第二步,基于非线性干扰观测器的自适应速度控制器设计:定义速度误差面SV=V-Vd,使用自适应模糊2型***对fV进行在线逼近:
定义速度子***干扰DV=ΩV+εV,则有:
设计如下非线性干扰观测器对DV进行在线观测:
式中:和分别是和的估计值,且定义为模糊2型***的估计误差;为干扰DV的估计值,且有zV为非线性干扰观测器的中间状态量;βV为待设计正数,则有:
则飞行器速度子***的控制量设计如下:
式中:kV0为待设计正数。
第三步,基于非线性干扰观测器的航迹角子***的控制器设计:首先,把飞行器的高度指令hd转化为航迹角指令γd:
式中:kh0为待设计正数,定义航迹角误差面Sγ=γ-γd,使用自适应模糊2型***对fγ进行在线逼近:
定义航迹角子***干扰Dγ=Ωγ+εγ,则有:
设计如下非线性干扰观测器对Dγ进行在线观测:
式中:和分别是和的估计值,且定义为模糊2型***的估计误差;为干扰Dγ的估计值,且有zγ为非线性干扰观测器的中间状态量;βγ为待设计正数,则有:
则虚拟控制输入可设计为:
式中:kγ0是待设计的正参数;
第四步,基于非线性干扰观测器的攻角子***的控制器设计:定义攻角误差面Sα=α-αd,使用自适应模糊2型***对fα进行在线逼近:
定义攻角子***干扰Dα=Ωα+εα,则有:
由于Dα未知,会对跟踪精度造成影响,因此设计如下非线性干扰观测器对Dα进行在线观测:
式中:和分别是和的估计值,且定义为模糊2型***的估计误差;为干扰Dα的估计值,且有zα为非线性干扰观测器的中间状态量;βα为待设计正数,则有:
则虚拟控制输入可设计为:
式中:kα0是待设计的正参数;
第五步,基于非线性干扰观测器的俯仰角速率子***的控制器设计:定义俯仰角速率误差面SQ=Q-Qd,使用自适应模糊2型***对fQ进行在线逼近:
定义俯仰角速率子***干扰DQ=ΩQ+εQ,则有:
设计如下非线性干扰观测器对DQ进行在线观测:
式中:和分别是和的估计值,且定义为模糊2型***的估计误差;为干扰DQ的估计值,且有zQ为非线性干扰观测器的中间状态量;βQ为待设计正数,则有:
则飞行器高度子***的控制量设计如下:
式中:kQ0为待设计正数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910276521.2A CN110231774A (zh) | 2019-04-08 | 2019-04-08 | 干扰观测变进气道高超声速飞行器模糊协调控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910276521.2A CN110231774A (zh) | 2019-04-08 | 2019-04-08 | 干扰观测变进气道高超声速飞行器模糊协调控制方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110231774A true CN110231774A (zh) | 2019-09-13 |
Family
ID=67860712
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910276521.2A Pending CN110231774A (zh) | 2019-04-08 | 2019-04-08 | 干扰观测变进气道高超声速飞行器模糊协调控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110231774A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116880198A (zh) * | 2023-07-25 | 2023-10-13 | 杭州牧星科技有限公司 | 一种用于超声速无人机的动力装备自适应控制***及其方法 |
CN117519257A (zh) * | 2024-01-04 | 2024-02-06 | 中国人民解放军国防科技大学 | 一种基于反步法的超声速巡航高度控制方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4243922A (en) * | 1971-12-31 | 1981-01-06 | Vereinigte Flugtechnische Werke-Fokker Gmbh | Aircraft control system |
CN103592847A (zh) * | 2013-10-30 | 2014-02-19 | 天津大学 | 一种基于高增益观测器的高超声速飞行器非线性控制方法 |
CN105278545A (zh) * | 2015-11-04 | 2016-01-27 | 北京航空航天大学 | 适用于高超声速机动飞行的自抗扰轨迹线性化控制方法 |
CN106406102A (zh) * | 2016-12-20 | 2017-02-15 | 哈尔滨工业大学 | 一种含干扰观测器的高超声速飞行器跟踪控制方法 |
CN107632518A (zh) * | 2017-09-05 | 2018-01-26 | 西北工业大学 | 基于高增益观测器的高超声速飞行器神经网络复合学习控制方法 |
CN108375907A (zh) * | 2018-03-28 | 2018-08-07 | 北京航空航天大学 | 基于神经网络的高超声速飞行器自适应补偿控制方法 |
-
2019
- 2019-04-08 CN CN201910276521.2A patent/CN110231774A/zh active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4243922A (en) * | 1971-12-31 | 1981-01-06 | Vereinigte Flugtechnische Werke-Fokker Gmbh | Aircraft control system |
CN103592847A (zh) * | 2013-10-30 | 2014-02-19 | 天津大学 | 一种基于高增益观测器的高超声速飞行器非线性控制方法 |
CN105278545A (zh) * | 2015-11-04 | 2016-01-27 | 北京航空航天大学 | 适用于高超声速机动飞行的自抗扰轨迹线性化控制方法 |
CN106406102A (zh) * | 2016-12-20 | 2017-02-15 | 哈尔滨工业大学 | 一种含干扰观测器的高超声速飞行器跟踪控制方法 |
CN107632518A (zh) * | 2017-09-05 | 2018-01-26 | 西北工业大学 | 基于高增益观测器的高超声速飞行器神经网络复合学习控制方法 |
CN108375907A (zh) * | 2018-03-28 | 2018-08-07 | 北京航空航天大学 | 基于神经网络的高超声速飞行器自适应补偿控制方法 |
Non-Patent Citations (1)
Title |
---|
苏沛华: "变几何进气道高超声速飞行器的建模与控制", 《中国优秀硕士学位论文全文数据库(电子期刊) 工程科技II辑》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116880198A (zh) * | 2023-07-25 | 2023-10-13 | 杭州牧星科技有限公司 | 一种用于超声速无人机的动力装备自适应控制***及其方法 |
CN116880198B (zh) * | 2023-07-25 | 2024-03-12 | 杭州牧星科技有限公司 | 一种用于超声速无人机的动力装备自适应控制***及其方法 |
CN117519257A (zh) * | 2024-01-04 | 2024-02-06 | 中国人民解放军国防科技大学 | 一种基于反步法的超声速巡航高度控制方法 |
CN117519257B (zh) * | 2024-01-04 | 2024-03-29 | 中国人民解放军国防科技大学 | 一种基于反步法的超声速巡航高度控制方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10787242B1 (en) | Variable camber continuous aerodynamic control surfaces and methods for active wing shaping control | |
Hull | Fundamentals of airplane flight mechanics | |
CN109614633A (zh) | 一种复合式旋翼飞行器非线性建模方法及配平方法 | |
CN109358634B (zh) | 一种高超声速飞行器鲁棒自适应控制方法 | |
CN107390531B (zh) | 参数学习有限时间收敛的高超声速飞行器控制方法 | |
CN107844128B (zh) | 一种基于复合比例导引的高超声速飞行器巡航段制导方法 | |
CN113778129B (zh) | 一种干扰补偿的高超声速变后掠翼飞行器跟踪控制方法 | |
CN109062055A (zh) | 一种基于Back-stepping鲁棒自适应动态面的近空间飞行器控制*** | |
CN109703769B (zh) | 一种基于预瞄策略的空中加油对接控制方法 | |
CN110568765A (zh) | 面向攻角跟踪的高超声速飞行器非对称输出受限控制方法 | |
CN113051662B (zh) | 一种基于cfd和datcom的折叠翼尖变体飞行器的气动建模及其性能评估方法 | |
CN111007720A (zh) | 高超声速飞行器多传感器故障的非线性自愈合控制方法 | |
CN114721266B (zh) | 飞机舵面结构性缺失故障情况下的自适应重构控制方法 | |
CN110231774A (zh) | 干扰观测变进气道高超声速飞行器模糊协调控制方法 | |
Anderson et al. | A comparison of aerodynamics models for optimizing the takeoff and transition of a bi-wing tailsitter | |
CN111240204A (zh) | 一种基于模型参考滑模变结构控制的巡飞弹控制方法 | |
Wilde et al. | Application of circulation control for three-axis control of a tailless flight vehicle | |
CN211685678U (zh) | 一种多旋翼无人机实时尾迹的仿真分析*** | |
CN115685764B (zh) | 变翼展飞行器任务自适应的抗干扰跟踪控制方法及*** | |
Luo et al. | Carrier-based aircraft precision landing using direct lift control based on incremental nonlinear dynamic inversion | |
CN110989397A (zh) | 一种航空器失事搜寻仿真方法及*** | |
Ricci et al. | Design and Experimental Validation of Gust Load Alleviation Systems based on Static Output Feedback | |
CN114237266A (zh) | 一种基于l1自适应的扑翼飞行姿态控制方法 | |
CN112441253A (zh) | 一种多旋翼无人机实时尾迹的仿真分析***及方法 | |
CN111061283A (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 | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20190913 |