CN107065554A - 近空间可变翼飞行器自适应小翼切换控制***及工作方法 - Google Patents
近空间可变翼飞行器自适应小翼切换控制***及工作方法 Download PDFInfo
- Publication number
- CN107065554A CN107065554A CN201710285460.7A CN201710285460A CN107065554A CN 107065554 A CN107065554 A CN 107065554A CN 201710285460 A CN201710285460 A CN 201710285460A CN 107065554 A CN107065554 A CN 107065554A
- Authority
- CN
- China
- Prior art keywords
- mrow
- mtd
- mtr
- msub
- delta
- 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
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
技术领域
本发明属于航空航宇推进控制技术领域,具体指代一种近空间可变翼飞行器自适应小翼切换控制***及工作方法。
背景技术
近空间飞行器所采用的推进动力是无需自带氧化剂的超声速燃烧冲压式发动机,气动布局为机体/发动机一体化设计。近空间飞行器的弹性机体与飞行推进***之间存在很强的耦合性,飞行器模型的非线性特性十分严重,飞行过程具有快速时变性,近空间飞行器的气动特性变化剧烈,各种不确定性严重,传统的经典控制方法无法很好地满足飞行控制***稳定性和强鲁棒性的性能要求。因此要保证飞行控制***的实时性,鲁棒性和稳定性,对响应速度和控制精度提出了更高的要求,这极大地推动了先进的控制方法和控制理论关键技术的发展。
近空间可变翼飞行器是具有可伸缩小翼的近空间飞行器,是为了解决起飞爬升段升力不足、升阻比过小,满足飞行速度和飞行高度包络范围大的特点而设计的一种可变翼飞行器。近空间可变翼飞行器兼具近空间飞行器与可变翼飞行器的特点,可以用于高空高速飞行,可变翼的特点可以使飞行器根据不同的飞行环境和飞行状态选择伸出或是收回小翼,一般在低速飞行时伸出小翼,提高升力,在高速飞行时,收回小翼以减小阻力,减小燃油消耗。
近空间可变翼飞行器所处的飞行环境、自身复杂多变的气动特性对飞行控制***的设计带来了很多技术上的难点。第一,飞行控制***必须满足稳定性要求。大跨度飞行包络,严重的外界干扰,弹性形变、高温和低密度流效应等因素会严重影响***的稳定性。第二,飞行控制***必须满足鲁棒性要求。在高动压、高速环境下,各种外界干扰和内部参数变化要求飞行控制***必须具有较强的鲁棒性。飞行器的所处大气环境复杂,使得飞行器异常敏感,时变性强,很容易产生基础结构失真和参数不确定性;第三,强耦合和非线性特点要求***的协调控制。近空间可变翼飞行器多采用乘波体或升力体气动布局以保证机动飞行的大迎角姿态,采用机体/发动机一体化设计可以保证高速飞行时不解体。第四,实时性要求。在高速飞行时,飞行参数具有激烈快时变的特征,而气动舵面的控制效果反而锐减,***反应时间加长,会出现控制延时问题。在控制器设计时要充分考虑机翼变形的实时性,控制算法的复杂度,要避免控制参数过多,提高算法的运行速度。第五,要满足约束条件。飞行控制***在保证控制精度的同时,还要满足一些约束条件,例如执行机构的饱和约束、迎角和侧滑角约束,在爬升和再入返回段的热流约束、为保证机体结构强度而设置的动压约束和过载约束等。
滑模控制方法通过设计不连续的控制器,迫使***产生滑动运动模态,一旦***进行滑动模态,***将对不确定性和干扰具有完全不变性。然而,实际***由于切换装置不可避免地存在惯性,变结构***在不同的控制逻辑中来回切换将导致实际滑动模态不是准确地发生在切换面上,容易引起***的剧烈抖动,不能保证***的鲁棒性。
增益预置控制方法在控制***设计中的应用比较成熟并取得了一定的成果,当非线性***变化范围较大时,采用这种方法需要设计多个平衡点,整个控制器的稳定性难以得到保证。在高超声速飞行器大迎角和高机动状态下,飞行状态呈现强非线性和高耦合性,增益预置方法无法满足性能指标的要求。
动态逆控制方法通过被控对象非线性耦合特性的准确建模,在线构成非线性耦合时变控制器,以抵消对象的非线性耦合时变特性,使***成为伪线性***。但是动态逆方法对建模误差敏感,且通常情况下,非线性***精确建模非常困难,一旦建模与实际***有差别,非线性耦合特性的对消就会有影响,导致控制性能的恶化,不能保证鲁棒性。
发明内容
针对于上述现有技术的不足,本发明的目的在于提供一种可变翼高超声速飞行器自适应控制***及工作方法,可有效解决近空间可变翼飞行器在小翼伸出与收回切换过程中的不确定性的飞行控制问题。
为达到上述目的,本发明的一种可变翼高超声速飞行器自适应控制***,其用于生成近空间可变翼飞行器的控制输入量u(t),以保证近空间可变翼飞行器飞行速度V、飞行俯仰角θ能快速跟踪到输入的飞行速度参考值Vr、飞行俯仰角参考值θr;该控制***包括:参考模型、标称控制器和自适应更新模块;其中,
参考模型,根据模型阶次,设计一阶、二阶参考模型:
式中,Δym(t)=[ΔVm,Δθm]T为参考模型输出,Δr(t)为参考模型输入,Wm(s)为参考模型的传递函数,为关联矩阵,p1,p2,p3为配置的稳定极点;
标称控制器,其形式如下:
其中,是标称控制器参数,
为高频增益矩阵;
标称控制器的作用下得到飞行器的输出Δy(t)=G(s)Δu(t),此时的输出Δy(t)能跟踪上述参考模型输出Δym(t);
自适应更新模块,其形式如下:
式中,K2为自适应更新参数,其为标称控制器参数的趋近律,K2选取的值为更新初值,自适应更新函数表达式如下:
误差更新函数为:
其中,定义一些辅助变量:
为估计误差,
为自适应更新函数中间变量,
为误差更新函数中间变量,
为归一化函数,
ω(t)=[ΔV,Δγ,Δα,Δq,Δθ,ΔVc,Δθc]T为新的状态量函数,
SP∈RM×M,为自适应增益矩阵;
自适应更新模块自动更新控制量Δu(t),得到飞行器的输出Δy(t)=G(s)Δu(t),此时的输出Δy(t)能精确跟踪参考模型输出Δym(t)。
优选地,所述的参考模型的传递函数
为一阶和二阶环节组成,p1,p2,p3为配置的稳定极点;速度的相对阶次是一阶,所以选择一阶惯性环节;俯仰角的相对阶次是二阶,所以选择二阶惯性环节。
优选地,所述的标称控制器和自适应更新模块的参数:
其中,Δδe(t)为升降舵增量,Δβc(t)为油门百分比增量,ΔV(t)为速度增量,Δγ(t)为航迹角增量,Δα(t)为迎角增量,Δq(t)为俯仰角速率增量,Δθ(t)为俯仰角增量,ΔVc(t)为速度指令信号,Δθc(t)为俯仰角指令信号,βc为可变翼高超声速飞行器的发动机节流阀调定值,δe为可变翼高超声速飞行器的升降舵偏转值。
优选地,所述的近空间可变翼飞行器的控制输入量u(t)=Δu(t)+u*,其中(x*,u*)为线性化的平衡点。
本发明的一种可变翼高超声速飞行器自适应控制***的工作方法,包括步骤如下:
1)将近空间飞行器非线性模型线性化得到状态方程表达式如下:
Δy=CΔx
式中,Δu(t)=[Δδe(t),Δηe(t)]T,Δy(t)=[ΔV(t),Δθ(t)]T,
Δx(t)=[ΔV(t),Δγ(t),Δα(t),Δq(t),Δθ(t)]T;
2)根据速度与俯仰角的相对阶次,选择参考模型如下:
Δym(t)=diag[1/(s+p1),1/(s+p2)(s+p3)]Δr(t)
选择速度与俯仰角指令作为参考模型输入Δr(t)=[ΔVc,Δθc]T,得到参考模型输出Δym(t)=[ΔVm,Δθm]T;
3)将标称控制器输出加上线性化平衡状态值u*,即u=Δu+u*输入到飞行器非线性模型,得到飞行器输出Δy(t),使其跟踪Δym(t);
4)在标称控制器作为初始输入基础之上加入自适应更新模块,得到自适应控制律即将u=Δu+u*输入到飞行器非线性模型;飞行器的小翼在25秒之后缓慢伸出,30秒时完全伸出,自适应更新模块在小翼伸出时更新参数,保证对速度与俯仰角的稳定跟踪。
本发明的有益效果:
1、本发明具有良好的全局稳定性,能有效利用基于状态方程设计控制器不依赖于非线性***本身的特性,基于该方法设计出的控制***在处理近空间飞行器问题时拥有更大的灵活性,可有效地保证其全局稳定性。
2、当控制***小翼由收缩到伸出或由伸出到收缩的情况下,自适应更新律能有效抵消参数不确定对***的影响,使得近空间可变翼飞行器能快速有效的响应跟踪信号,具有较好的跟踪性能和鲁棒性能。
3、所设计的自适应控制器的自适应更新律是实时、动态更新的过程,可以提高***的实时性。
附图说明
图1为本发明控制***的结构原理示意图;
图2为本发明控制***的速度跟踪响应曲线图;
图3为本发明控制***的俯仰角跟踪响应曲线图。
具体实施方式
为了便于本领域技术人员的理解,下面结合实施例与附图对本发明作进一步的说明,实施方式提及的内容并非对本发明的限定。
本发明基于状态反馈输出跟踪的模型参考自适应控制***设计方法是将建立的参考模型作为控制器设计的基础,因此它不依赖与飞行器的非线性模型,在实际应用中较为简单,控制器分为标称控制器和自适应更新律,根据状态方程设计标称控制器参数,并将标称控制器参数作为自适应更新律的初值,并借助李亚普诺夫函数设计自适应更新律,以保证***的全局稳定性。该方法适用于线性控制,不依赖于数学模型,具有非常独特的优势,可以提高控制***的稳定性和鲁棒性能。
参照图1所示,本发明的一种可变翼高超声速飞行器自适应控制***,其应用于可变翼高超声速飞行器中,该可变翼高超声速飞行器的输入输出形式:y(t)=G(s)u(t),增量形式为Δy(t)=G(s)Δu(t),状态方程表达式为:Δy(t)=CΔx(t),其中G(s)为飞行器传递函数,y(t)为飞行器的输出,u(t)为飞行器的控制输入,x(t)为飞行器状态量,A,B,C为***矩阵;该控制***用于生成近空间可变翼飞行器的控制输入量u(t),以保证近空间可变翼飞行器飞行速度V、飞行俯仰角θ能快速跟踪到输入的飞行速度参考值Vr、飞行俯仰角参考值θr;该控制***包括:参考模型、标称控制器和自适应更新模块;
以某型近空间可变翼飞行器为例,其采用翼身融合布局,机体轮廓为三角形,大后掠机翼与机身采用翼身融合方式,升降舵布置在机翼后缘,机翼为可变形翼面,采用在主翼两侧伸缩小翼的方式提高升力和飞行效率。该飞行器在高超声速巡航飞行条件下的纵向运动模型描述为:
式中,状态量V,γ,α,q,θ分别表示飞行器速度,飞行航迹倾斜角,飞行迎角,俯仰角速率和飞行俯仰角;L,D,T分别表示飞行器升力、阻力和发动机推力,M表示俯仰力矩,Iyy表示纵向转动惯量,∫(x2+z2)dm=Iyy;
气动力和力矩表示为:
L=0.5ρV2sCL
D=0.5ρV2sCD
Myy=0.5ρV2sc[CM(α)+CM(δe)+CM(q)]
式中,ρ为大气密度,s为机翼浸润面积,CL,CD,CM分别为升力系数,阻力系数和俯仰力矩系数;
小翼伸缩模型如下:
式中,s为机翼面积,小翼在25s时刻开始伸出,伸出前机翼面积s=369m/s,在30s时刻完成小翼的伸出,小翼完全伸出机翼面积s=389m/s;小翼伸出过程中机翼面积变化影响飞行器升力系数CL、阻力系数CD和俯仰力矩系数CM发生变化。
将上述近空间可变翼飞行器非线性模型线性化得到状态方程表达式为:
标称控制器和自适应更新模块的设计目标是当小翼切换时出现不确定性扰动的情况下,使得近空间可变翼飞行器能稳定跟踪给定速度参考信号和俯仰角参考信号yr=[Vr(t) θr(t)]T,t≥0。
在标称控制器和自适应更新模块设计之前,作如下假设:
假设1:G(s)=C(sI-A)-1B的所有零点稳定,并且***(A,B,C)可控且可观测。
假设2:G(s)严格有理,满秩,并且有一个对角关联矩阵
ξm(s)=diag{d1(s),d2(s),…,dM(s)}
式中,为首一稳定多项式,也就是说高频增益矩阵定义为:
该矩阵有界,并且非奇异。
假设3:存在已知矩阵SP∈RM×M满足
参考模型,根据模型阶次,设计一阶、二阶参考模型:
式中,Δym(t)=[ΔVm,Δθm]T为参考模型输出,Δr(t)为参考模型输入,Wm(s)为参考模型的传递函数,为关联矩阵,参考模型由一阶和二阶环节组成,p1,p2,p3为配置的稳定极点,速度的相对阶次是一阶,所以选择一阶惯性环节,俯仰角的相对阶次是二阶,所以选择二阶惯性环节。
标称控制器,其形式如下:
其中,是标称控制器参数, 为高频增益矩阵;
标称控制器的作用下得到飞行器的输出Δy(t)=G(s)Δu(t),此时的输出Δy(t)能跟踪上述参考模型输出Δym(t);
自适应更新模块,将标称控制器应用到***中,得到闭环***:
并且输出跟踪误差为:
式中,ω(t)=[xT(t),rT(t)]T,
为得到自适应更新模块的参数向量ΘT(t)的更新律,令Ψ*=Kp,h(s)=1/(s+p2)(s+p3),相对阶ρM的稳定多项式为1/(s+p2)(s+p3),定义估计误差为:
ε(t)=ξm(s)h(s)[Δy-Δym](t)+Ψ(t)ξ(t)
式中,Ψ(t)是Ψ*的估计误差,并且
ζ(t)=h(s)[ω](t)
ξ(t)=ΘT(t)ζ(t)-h(s)[Δu](t)
估计误差还可以写成
式中,
选择以下自适应参数更新律
式中,Γ=ΓT>0,且
得到自适应控制器如下:
式中,K2为标称控制参数的趋近律,且更新律为
定义跟踪误差e(t)=Δy(t)-Δym(t)满足
证明:可以找到一个正定的方程:
其导数为:
根据李雅普诺夫稳定性定理可知,闭环***信号有界且输出能够实现渐进跟踪参考模型。
为了验证本实例的效果,对其进行数值仿真验证。为了描述飞行器小翼伸缩的不确定性,仿真中小翼在25s时刻开始伸出,伸出前机翼面积s=369m/s,在30s时刻完成小翼的伸出,小翼完全伸出机翼面积s=389m/s。小翼伸出过程中机翼面积变化影响飞行器升力系数CL、阻力系数CD和俯仰力矩系数CM发生变化。利用本实例控制***所得到的近空间可变翼飞行器高度和俯仰角的跟踪响应曲线如图2、3所示。
本发明的一种可变翼高超声速飞行器自适应控制***的工作方法,包括步骤如下:
1)将近空间飞行器非线性模型线性化得到状态方程表达式如下:
Δy=CΔx
式中,Δu(t)=[Δδe(t),Δβc(t)]T,Δy(t)=[ΔV(t),Δθ(t)]T,
Δx(t)=[ΔV(t),Δγ(t),Δα(t),Δq(t),Δθ(t)]T;
2)根据速度与俯仰角的相对阶次,选择参考模型如下:
Δym(t)=diag[1/(s+p1),1/(s+p2)(s+p3)]Δr(t)
选择速度与俯仰角指令作为参考模型输入Δr(t)=[ΔVc,Δθc]T,得到参考模型输出Δym(t)=[ΔVm,Δθm]T;
3)将标称控制器输出加上线性化平衡状态值u*,即u=Δu+u*输入到飞行器非线性模型,得到飞行器输出Δy(t),使其跟踪Δym(t);
4)在标称控制器作为初始输入基础之上加入自适应更新模块,得到自适应控制律即将u=Δu+u*输入到飞行器非线性模型;飞行器的小翼在25秒之后缓慢伸出,30秒时完全伸出,自适应更新模块在小翼伸出时更新参数,保证对速度与俯仰角的稳定跟踪。
本发明具体应用途径很多,以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以作出若干改进,这些改进也应视为本发明的保护范围。
Claims (5)
1.一种可变翼高超声速飞行器自适应控制***,其用于生成近空间可变翼飞行器的控制输入量u(t),以保证近空间可变翼飞行器飞行速度V、飞行俯仰角θ能快速跟踪到输入的飞行速度参考值Vr、飞行俯仰角参考值θr;其特征在于,该控制***包括:参考模型、标称控制器和自适应更新模块;其中,
参考模型,根据模型阶次,设计一阶、二阶参考模型:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;y</mi>
<mi>m</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>W</mi>
<mi>m</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mi>&Delta;</mi>
<mi>r</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<msubsup>
<mi>&xi;</mi>
<mi>m</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mi>&Delta;</mi>
<mi>r</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mfrac>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>+</mo>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
</mfrac>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mfrac>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>+</mo>
<msub>
<mi>p</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
<mo>(</mo>
<mi>s</mi>
<mo>+</mo>
<msub>
<mi>p</mi>
<mn>3</mn>
</msub>
<mo>)</mo>
</mrow>
</mfrac>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>&Delta;</mi>
<mi>r</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
式中,Δym(t)=[ΔVm,Δθm]T为参考模型输出,Δr(t)为参考模型输入,Wm(s)为参考模型的传递函数,为关联矩阵,p1,p2,p3为配置的稳定极点;
标称控制器,其形式如下:
<mrow>
<mi>&Delta;</mi>
<mi>u</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mi>K</mi>
<mn>1</mn>
<mrow>
<mo>*</mo>
<mi>T</mi>
</mrow>
</msubsup>
<mi>&Delta;</mi>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msubsup>
<mi>K</mi>
<mn>2</mn>
<mo>*</mo>
</msubsup>
<mi>&Delta;</mi>
<mi>r</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
其中,为标称控制器参数, 为高频增益矩阵;
标称控制器的作用下得到飞行器的输出Δy(t)=G(s)Δu(t),此时的输出Δy(t)能跟踪上述参考模型输出Δym(t);
自适应更新模块,其形式如下:
<mrow>
<mi>&Delta;</mi>
<mi>u</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mi>K</mi>
<mn>1</mn>
<mi>T</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>&Delta;</mi>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>K</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>&Delta;</mi>
<mi>r</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
式中,K2为自适应更新参数,其为标称控制器参数的趋近律,K2选取的值为更新初值,自适应更新函数表达式如下:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msubsup>
<mi>K</mi>
<mn>1</mn>
<mi>T</mi>
</msubsup>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>K</mi>
<mn>2</mn>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<msup>
<mover>
<mi>&Theta;</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>T</mi>
</msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mo>-</mo>
<msub>
<mi>S</mi>
<mi>p</mi>
</msub>
<mi>&epsiv;</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>&zeta;</mi>
<mi>T</mi>
</msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msup>
<mi>m</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
误差更新函数为:
<mrow>
<mover>
<mi>&Psi;</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mo>-</mo>
<mi>&Gamma;</mi>
<mi>&epsiv;</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>&xi;</mi>
<mi>T</mi>
</msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msup>
<mi>m</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
其中,定义一些辅助变量:
为估计误差,
为自适应更新函数中间变量,
为误差更新函数中间变量,
为归一化函数,
ω(t)=[ΔV,Δγ,Δα,Δq,Δθ,ΔVc,Δθc]T为新的状态量函数,
SP∈RM×M,为自适应增益矩阵;
自适应更新模块自动更新控制量Δu(t),得到飞行器的输出Δy(t)=G(s)Δu(t),此时的输出Δy(t)能精确跟踪参考模型输出Δym(t)。
2.根据权利要求1所述的可变翼高超声速飞行器自适应控制***,其特征在于,所述的参考模型的传递函数
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mfrac>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>+</mo>
<msub>
<mi>p</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
</mfrac>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mfrac>
<mn>1</mn>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>+</mo>
<msub>
<mi>p</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
<mo>(</mo>
<mi>s</mi>
<mo>+</mo>
<msub>
<mi>p</mi>
<mn>3</mn>
</msub>
<mo>)</mo>
</mrow>
</mfrac>
</mtd>
</mtr>
</mtable>
</mfenced>
为一阶和二阶环节组成,p1,p2,p3为配置的稳定极点;速度的相对阶次是一阶,所以选择一阶惯性环节;俯仰角的相对阶次是二阶,所以选择二阶惯性环节。
3.根据权利要求1所述的可变翼高超声速飞行器自适应控制***,其特征在于,所述的标称控制器和自适应更新模块的参数:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<mi>u</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;&delta;</mi>
<mi>e</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;&beta;</mi>
<mi>c</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>&Delta;</mi>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<mi>V</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<mi>&gamma;</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<mi>&alpha;</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<mi>q</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<mi>&theta;</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>&Delta;</mi>
<mi>r</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;V</mi>
<mi>c</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;&theta;</mi>
<mi>c</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,Δδe(t)为升降舵增量,Δβc(t)为油门百分比增量,ΔV(t)为速度增量,Δγ(t)为航迹角增量,Δα(t)为迎角增量,Δq(t)为俯仰角速率增量,Δθ(t)为俯仰角增量,ΔVc(t)为速度指令信号,Δθc(t)为俯仰角指令信号,βc为可变翼高超声速飞行器的发动机节流阀调定值,δe为可变翼高超声速飞行器的升降舵偏转值。
4.根据权利要求1所述的可变翼高超声速飞行器自适应控制***,其特征在于,所述的近空间可变翼飞行器的控制输入量u(t)=Δu(t)+u*,其中(x*,u*)为线性化的平衡点。
5.一种可变翼高超声速飞行器自适应控制***的工作方法,其特征在于,包括步骤如下:
1)将近空间飞行器非线性模型线性化得到状态方程表达式如下:
<mrow>
<mi>&Delta;</mi>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<mi>A</mi>
<mi>&Delta;</mi>
<mi>x</mi>
<mo>+</mo>
<mi>B</mi>
<mi>&Delta;</mi>
<mi>u</mi>
</mrow>
Δy=CΔx
式中,Δu(t)=[Δδe(t),Δηe(t)]T,Δy(t)=[ΔV(t),Δθ(t)]T,
Δx(t)=[ΔV(t),Δγ(t),Δα(t),Δq(t),Δθ(t)]T;
2)根据速度与俯仰角的相对阶次,选择参考模型如下:
Δym(t)=diag[1/(s+p1),1/(s+p2)(s+p3)]Δr(t)
选择速度与俯仰角指令作为参考模型输入Δr(t)=[ΔVc,Δθc]T,得到参考模型输出Δym(t)=[ΔVm,Δθm]T;
3)将标称控制器输出加上线性化平衡状态值u*,即u=Δu+u*输入到飞行器非线性模型,得到飞行器输出Δy(t),使其跟踪Δym(t);
4)在标称控制器作为初始输入基础之上加入自适应更新模块,得到自适应控制器即将u=Δu+u*输入到飞行器非线性模型;飞行器的小翼在25秒之后缓慢伸出,30秒时完全伸出,自适应更新模块在小翼伸出时更新参数,保证对速度与俯仰角的稳定跟踪。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710285460.7A CN107065554B (zh) | 2017-04-27 | 2017-04-27 | 近空间可变翼飞行器自适应小翼切换控制***及工作方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710285460.7A CN107065554B (zh) | 2017-04-27 | 2017-04-27 | 近空间可变翼飞行器自适应小翼切换控制***及工作方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107065554A true CN107065554A (zh) | 2017-08-18 |
CN107065554B CN107065554B (zh) | 2019-04-09 |
Family
ID=59604693
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710285460.7A Expired - Fee Related CN107065554B (zh) | 2017-04-27 | 2017-04-27 | 近空间可变翼飞行器自适应小翼切换控制***及工作方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107065554B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107608385A (zh) * | 2017-10-25 | 2018-01-19 | 南京航空航天大学 | 无人直升机编队控制方法及控制*** |
CN107807663A (zh) * | 2017-11-29 | 2018-03-16 | 南京航空航天大学 | 基于自适应控制的无人机编队保持控制方法 |
CN107942651A (zh) * | 2017-10-20 | 2018-04-20 | 南京航空航天大学 | 一种近空间飞行器控制*** |
CN108595755A (zh) * | 2018-03-21 | 2018-09-28 | 南京航空航天大学 | 一种新的火星探测飞行器面向控制的快速建模方法 |
CN110320927A (zh) * | 2019-07-26 | 2019-10-11 | 中国人民解放军国防科技大学 | 智能变形飞行器的飞行控制方法及*** |
CN112147896A (zh) * | 2020-09-28 | 2020-12-29 | 中国科学院数学与***科学研究院 | 一种非标准型mimo离散非线性***的自适应控制方法及*** |
CN113377122A (zh) * | 2021-07-06 | 2021-09-10 | 西北工业大学 | 一种可栖落机动的变体飞行器切换自适应控制方法 |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4609855A (en) * | 1984-11-09 | 1986-09-02 | Motorola, Inc. | Adaptive servomotor control |
US20050038638A1 (en) * | 2003-07-25 | 2005-02-17 | Dorin Comaniciu | Density morphing and mode propagation for Bayesian filtering |
CN101937233A (zh) * | 2010-08-10 | 2011-01-05 | 南京航空航天大学 | 近空间高超声速飞行器非线性自适应控制方法 |
CN103235504A (zh) * | 2013-03-25 | 2013-08-07 | 上海交通大学 | 基于直接自适应控制重构的大型民用飞机飞行控制方法 |
CN103425135A (zh) * | 2013-07-30 | 2013-12-04 | 南京航空航天大学 | 一种具有输入饱和的近空间飞行器鲁棒控制方法 |
CN104950672A (zh) * | 2015-06-10 | 2015-09-30 | 北京理工大学 | 一种最优积分滑模控制方法 |
CN104950898A (zh) * | 2015-06-10 | 2015-09-30 | 北京理工大学 | 一种再入飞行器全阶非奇异终端滑模姿态控制方法 |
CN105182742A (zh) * | 2015-07-23 | 2015-12-23 | 北京航空航天大学 | 一种弹性体飞行器自适应受限跟踪控制间接法 |
CN105807613A (zh) * | 2016-03-08 | 2016-07-27 | 清华大学 | 一种基于涡扇发动机多工作点的鲁棒补偿控制方法 |
CN105911864A (zh) * | 2016-06-02 | 2016-08-31 | 西北工业大学 | 基于虚拟平衡点不变集的不确定受限***鲁棒控制方法 |
CN105938368A (zh) * | 2016-05-26 | 2016-09-14 | 北京航空航天大学 | 一种基于有限时间摩擦估计的挠性航天器姿态控制方法 |
CN106444368A (zh) * | 2015-11-18 | 2017-02-22 | 南京航空航天大学 | 具有输入非线性的近空间飞行器预设性能姿态跟踪控制方法 |
-
2017
- 2017-04-27 CN CN201710285460.7A patent/CN107065554B/zh not_active Expired - Fee Related
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4609855A (en) * | 1984-11-09 | 1986-09-02 | Motorola, Inc. | Adaptive servomotor control |
US20050038638A1 (en) * | 2003-07-25 | 2005-02-17 | Dorin Comaniciu | Density morphing and mode propagation for Bayesian filtering |
CN101937233A (zh) * | 2010-08-10 | 2011-01-05 | 南京航空航天大学 | 近空间高超声速飞行器非线性自适应控制方法 |
CN103235504A (zh) * | 2013-03-25 | 2013-08-07 | 上海交通大学 | 基于直接自适应控制重构的大型民用飞机飞行控制方法 |
CN103425135A (zh) * | 2013-07-30 | 2013-12-04 | 南京航空航天大学 | 一种具有输入饱和的近空间飞行器鲁棒控制方法 |
CN104950672A (zh) * | 2015-06-10 | 2015-09-30 | 北京理工大学 | 一种最优积分滑模控制方法 |
CN104950898A (zh) * | 2015-06-10 | 2015-09-30 | 北京理工大学 | 一种再入飞行器全阶非奇异终端滑模姿态控制方法 |
CN105182742A (zh) * | 2015-07-23 | 2015-12-23 | 北京航空航天大学 | 一种弹性体飞行器自适应受限跟踪控制间接法 |
CN106444368A (zh) * | 2015-11-18 | 2017-02-22 | 南京航空航天大学 | 具有输入非线性的近空间飞行器预设性能姿态跟踪控制方法 |
CN105807613A (zh) * | 2016-03-08 | 2016-07-27 | 清华大学 | 一种基于涡扇发动机多工作点的鲁棒补偿控制方法 |
CN105938368A (zh) * | 2016-05-26 | 2016-09-14 | 北京航空航天大学 | 一种基于有限时间摩擦估计的挠性航天器姿态控制方法 |
CN105911864A (zh) * | 2016-06-02 | 2016-08-31 | 西北工业大学 | 基于虚拟平衡点不变集的不确定受限***鲁棒控制方法 |
Non-Patent Citations (1)
Title |
---|
吴雨珊 等: "基于回馈递推的可变翼高超声速飞行器智能非线性控制", 《哈尔滨工程大学学报》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107942651A (zh) * | 2017-10-20 | 2018-04-20 | 南京航空航天大学 | 一种近空间飞行器控制*** |
CN107942651B (zh) * | 2017-10-20 | 2020-11-03 | 南京航空航天大学 | 一种近空间飞行器控制*** |
CN107608385A (zh) * | 2017-10-25 | 2018-01-19 | 南京航空航天大学 | 无人直升机编队控制方法及控制*** |
CN107608385B (zh) * | 2017-10-25 | 2019-10-11 | 南京航空航天大学 | 无人直升机编队控制方法及控制*** |
CN107807663A (zh) * | 2017-11-29 | 2018-03-16 | 南京航空航天大学 | 基于自适应控制的无人机编队保持控制方法 |
CN107807663B (zh) * | 2017-11-29 | 2019-12-24 | 南京航空航天大学 | 基于自适应控制的无人机编队保持控制方法 |
CN108595755A (zh) * | 2018-03-21 | 2018-09-28 | 南京航空航天大学 | 一种新的火星探测飞行器面向控制的快速建模方法 |
CN110320927A (zh) * | 2019-07-26 | 2019-10-11 | 中国人民解放军国防科技大学 | 智能变形飞行器的飞行控制方法及*** |
CN112147896A (zh) * | 2020-09-28 | 2020-12-29 | 中国科学院数学与***科学研究院 | 一种非标准型mimo离散非线性***的自适应控制方法及*** |
CN113377122A (zh) * | 2021-07-06 | 2021-09-10 | 西北工业大学 | 一种可栖落机动的变体飞行器切换自适应控制方法 |
CN113377122B (zh) * | 2021-07-06 | 2022-10-28 | 西北工业大学 | 一种可栖落机动的变体飞行器切换自适应控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107065554B (zh) | 2019-04-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107065554A (zh) | 近空间可变翼飞行器自适应小翼切换控制***及工作方法 | |
CN104950898B (zh) | 一种再入飞行器全阶非奇异终端滑模姿态控制方法 | |
CN104155990B (zh) | 考虑攻角约束的高超声速飞行器俯仰通道姿态控制方法 | |
CN107942651B (zh) | 一种近空间飞行器控制*** | |
Zhang et al. | An analysis of the stability and chattering reduction of high-order sliding mode tracking control for a hypersonic vehicle | |
CN104950901B (zh) | 无人直升机姿态误差有限时间收敛非线性鲁棒控制方法 | |
CN105138003B (zh) | 多操纵面无人机直接升力控制方法 | |
CN103558857A (zh) | 一种btt飞行器的分布式复合抗干扰姿态控制方法 | |
CN106444807A (zh) | 一种栅格舵与侧喷流的复合姿态控制方法 | |
CN105182742A (zh) | 一种弹性体飞行器自适应受限跟踪控制间接法 | |
CN101708780A (zh) | 用于目标姿态跟踪的刚性航天器的控制方法 | |
CN109062055A (zh) | 一种基于Back-stepping鲁棒自适应动态面的近空间飞行器控制*** | |
CN102880052A (zh) | 基于时标功能分解的高超声速飞行器执行器饱和控制方法 | |
CN109470252A (zh) | 一种基于凸优化的垂直起降重复使用运载器快速轨迹优化方法 | |
CN107450323B (zh) | 高超声速飞行器再入段神经网络复合学习控制方法 | |
CN104252133A (zh) | 一种无人机纵向控制律平滑切换方法 | |
CN107085435A (zh) | 基于耦合分析的高超声速飞行器姿态协调控制方法 | |
CN104932531A (zh) | 一种基于滑模控制的四旋翼飞行器的最优抗输入饱和控制方法 | |
CN106707759A (zh) | 一种飞机Herbst机动控制方法 | |
CN108828957A (zh) | 基于切换机制的飞行器全局有限时间神经网络控制方法 | |
CN106527128A (zh) | 兼顾瞬态响应与鲁棒稳定性的飞行控制律设计新方法 | |
CN110320927A (zh) | 智能变形飞行器的飞行控制方法及*** | |
CN106843261A (zh) | 一种变体飞行器过渡段的张量积插值建模与控制方法 | |
CN106970633A (zh) | 抑制控制输入饱和的飞行控制方法 | |
CN106383523A (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190409 Termination date: 20210427 |
|
CF01 | Termination of patent right due to non-payment of annual fee |