CN107390531A - 参数学习有限时间收敛的高超声速飞行器控制方法 - Google Patents

参数学习有限时间收敛的高超声速飞行器控制方法 Download PDF

Info

Publication number
CN107390531A
CN107390531A CN201710789197.5A CN201710789197A CN107390531A CN 107390531 A CN107390531 A CN 107390531A CN 201710789197 A CN201710789197 A CN 201710789197A CN 107390531 A CN107390531 A CN 107390531A
Authority
CN
China
Prior art keywords
mrow
msub
mover
mfrac
centerdot
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
Application number
CN201710789197.5A
Other languages
English (en)
Other versions
CN107390531B (zh
Inventor
许斌
郭雨岩
程怡新
张睿
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Northwestern Polytechnical University
Shenzhen Institute of Northwestern Polytechnical University
Original Assignee
Northwestern Polytechnical University
Shenzhen Institute of Northwestern Polytechnical University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Northwestern Polytechnical University, Shenzhen Institute of Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN201710789197.5A priority Critical patent/CN107390531B/zh
Publication of CN107390531A publication Critical patent/CN107390531A/zh
Application granted granted Critical
Publication of CN107390531B publication Critical patent/CN107390531B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/0265—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric the criterion being a learning criterion
    • G05B13/027—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric the criterion being a learning criterion using neural networks only
    • 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)
  • Artificial Intelligence (AREA)
  • Evolutionary Computation (AREA)
  • Health & Medical Sciences (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Feedback Control In General (AREA)

Abstract

本发明公开了一种参数学习有限时间收敛的高超声速飞行器控制方法，用于解决现有高超声速飞行器控制方法实用性差的技术问题。技术方案是通过构造一阶滤波器与辅助信号获取建模误差信息，结合跟踪误差与辅助信号设计神经网络自适应律，针对高超声速飞行器模型中的未知非线性函数估计设计了RBF神经网络方法，将有限时间学习的思想引入神经网络权重更新律设计中，相比传统的神经网络方法，本发明通过构造一阶滤波器与辅助信号将建模误差信息引入权重更新律，能保证参数学习误差有限时间收敛，从而保证学习的快速性。由于采用神经网络学习对不确定性进行估计，无需进行模型线性参数化表达，可实现不确定高超声速飞行器控制，便于实际工程应用。

Description

参数学习有限时间收敛的高超声速飞行器控制方法
技术领域
本发明涉及一种高超声速飞行器控制方法，特别涉及一种参数学习有限时间收敛的高超声速飞行器控制方法。
背景技术
高超声速飞行器指能够在大气层内以超过五倍声速飞行的飞行器，由于其突出的飞行能力使得全球实时打击成为可能，因此受到国内外的广泛关注；由于飞行环境复杂且存在气动参数摄动，飞行器模型具有强不确定性，且模型线性参数化形式在实际中难以获得。
《Novel auxiliary error compensation design for the adaptive neuralcontrol of a constrained flexible air-breathing hypersonic vehicle》(XiangweiBu,Xiaoyan Wu,Zhen Ma,Rui Zhang,Jiaqi Huang,《Neurocomputing》,2016,171:313-324.)一文针对高超声速飞行器的未知非线性动态提出了一类神经网络估计方法，根据跟踪误差给出神经网络权重自适应律。然而，这类神经网络学习方法只依据Lyapunov理论设计自适应律以保证系统稳定性，不能保证学习速度，在实际系统中难以取得满意的跟踪性能。
发明内容
为了克服现有高超声速飞行器控制方法实用性差的不足，本发明提供一种参数学习有限时间收敛的高超声速飞行器控制方法。该方法通过构造一阶滤波器与辅助信号获取建模误差信息，结合跟踪误差与辅助信号设计神经网络自适应律，针对高超声速飞行器模型中的未知非线性函数估计设计了RBF神经网络方法，将有限时间学习的思想引入神经网络权重更新律设计中，相比传统的神经网络方法，本发明通过构造一阶滤波器与辅助信号将建模误差信息引入权重更新律，能保证参数学习误差有限时间收敛，从而保证学习的快速性。由于采用神经网络学习对不确定性进行估计，无需进行模型线性参数化表达，可实现不确定高超声速飞行器控制，便于实际工程应用，实用性好。
本发明解决其技术问题所采用的技术方案：一种参数学习有限时间收敛的高超声速飞行器控制方法，其特点是包括以下步骤：
步骤一、建立高超声速飞行器纵向通道动力学模型：
其中， 表示动压，ρ表示空气密度，CL、CD、CT、CM(α)、CM(q)、CM(δe)均为气动参数，表示平均气动弦长，S表示气动参考面积。V表示速度，γ表示航迹倾角，h表示高度，α表示攻角，q表示俯仰角速度，δe是舵偏角，β为节流阀开度；T、D、L和Myy分别代表推力、阻力、升力和俯仰转动力矩；m、Iyy、μ和r代表质量、俯仰轴的转动惯量、引力系数以及距地心的距离。
步骤二、定义高度跟踪误差eh＝h-hd，设计航迹角指令γd：
式中，hd为高度指令，为高度指令的一阶微分，kh>0，ki>0。考虑巡航段航迹角变化小，航迹角指令的一阶微分取为零。
取x1＝γ，x2＝θp，x3＝q，θp＝α+γ代表俯仰角；将式(3)-(5)写成如下严格反馈形式：
其中，fi,i＝1,3为未知函数，gi,i＝1,3为根据高超声速飞行器模型得到的已知函数。
步骤三、定义航迹角跟踪误差：
e1＝x1-γd (8)
设计俯仰角虚拟控制量：
式中，k1>0，为由RBF神经网络得到的f1估计值，其中为神经网络最优权重向量的估计值，θ1为RBF函数向量。
定义设计一阶滤波器如下：
式中xf1、θf1为x1、θ1通过滤波器(10)后获得的信号，设计辅助矩阵P1、辅助向量Q1、W1如下：
式中，l1>0，神经网络权重通过如下自适应律获得：
其中，γ1>0，Γ1为正定常数对角矩阵。
设计一阶滤波器如下：
式中x2d为x2c通过滤波器(13)后获得的信号，α2>0。
定义俯仰角跟踪误差：
e2＝x2-x2d (14)
设计俯仰角速率虚拟控制量：
式中，k2>0。
设计一阶滤波器如下：
式中x3d为x3c通过滤波器(16)后获得的信号，α3>0。
定义俯仰角速率跟踪误差：
e3＝x3-x3d (17)
设计舵偏角控制律如下：
式中，k3>0，为由RBF神经网络得到的f3估计值，其中为神经网络最优权重向量的估计值，θ3为RBF函数向量。
令设计一阶滤波器如下：
式中xf3、θf3为x3、θ3通过滤波器(19)后获得的信号，设计辅助矩阵P3、辅助向量Q3、W3如下：
式中，l3>0，神经网络权重通过如下自适应律获得：
其中，γ3>0，Γ3为正定常数对角矩阵。
步骤四、定义速度跟踪误差：
式中，Vd为速度指令。设计节流阀开度如下：
式中，kpV>0、kiV>0、kdV>0。
步骤五、根据得到的舵偏角δe和节流阀开度β，返回到高超声速飞行器的动力学模型(1)-(5)，对高度和速度进行跟踪控制。
本发明的有益效果是：该方法通过构造一阶滤波器与辅助信号获取建模误差信息，结合跟踪误差与辅助信号设计神经网络自适应律，针对高超声速飞行器模型中的未知非线性函数估计设计了RBF神经网络方法，将有限时间学习的思想引入神经网络权重更新律设计中，相比传统的神经网络方法，本发明通过构造一阶滤波器与辅助信号将建模误差信息引入权重更新律，能保证参数学习误差有限时间收敛，从而保证学习的快速性。由于采用神经网络学习对不确定性进行估计，无需进行模型线性参数化表达，可实现不确定高超声速飞行器控制，便于实际工程应用，实用性好。
下面结合附图和具体实施方式对本发明作详细说明。
附图说明
图1是本发明参数学习有限时间收敛的高超声速飞行器控制方法的流程图。
具体实施方式
参照图1。本发明参数学习有限时间收敛的高超声速飞行器控制方法具体步骤如下：
步骤一、建立高超声速飞行器纵向通道动力学模型：
其中，V表示速度，γ表示航迹倾角，h表示高度，α表示攻角，q表示俯仰角速度，δe是舵偏角，β为节流阀开度；T、D、L和Myy分别代表推力、阻力、升力和俯仰转动力矩；m、Iyy、μ和r代表质量、俯仰轴的转动惯量、引力系数以及距地心的距离；
相关的力矩及参数定义如下： CL＝0.6203α，CD＝0.6450α2+0.0043378α+0.003772，CM(α)＝-0.035α2+0.036617α+5.3261×10-6，CM(δe)＝0.0292(δe-α)。其中表示动压，ρ表示空气密度，CL、CD、CT、CM(α)、CM(q)、CM(δe)均为气动参数，表示平均气动弦长，S表示气动参考面积。
步骤二、定义高度跟踪误差eh＝h-hd，设计航迹角指令γd：
式中，hd为高度指令，由设计者给出，为高度指令的一阶微分，kh＝0.5,ki＝0.05。考虑巡航段航迹角变化小，航迹角指令的一阶微分取为零。
取x1＝γ，x2＝θp，x3＝q，θp＝α+γ代表俯仰角；可将式(3)-(5)写成如下严格反馈形式：
其中，fi,gi,i＝1,3为根据高超声速飞行器模型得到的未知非线性函数。上式中，
步骤三、定义航迹角跟踪误差：
e1＝x1-γd (8)
设计俯仰角虚拟控制量：
式中，k1＝2，为由RBF神经网络得到的f1估计值，其中为神经网络最优权重向量的估计值，θ1为RBF函数向量。
定义令x1、θ1通过一阶滤波器：
其中设计辅助矩阵P1、辅助向量Q1、W1如下：
式中，l1＝2，神经网络权重通过如下自适应律获得：
其中，γ1＝1，Γ1为对角元素为0.05的对角矩阵。
令x2c通过一阶滤波器得到x2d：
式中α2＝0.05。
定义俯仰角跟踪误差：
e2＝x2-x2d (14)
设计俯仰角速率虚拟控制量：
式中，k2＝2。
令x3c通过一阶滤波器得到x3d：
式中α3＝0.05。
定义俯仰角速率跟踪误差：
e3＝x3-x3d (17)
设计舵偏角控制律如下：
式中，k3＝5，为由RBF神经网络得到的f3估计值，其中为神经网络最优权重向量的估计值，θ3为RBF函数向量。
令分别使x3、θ3通过一阶滤波器获得xf3、θf3：
其中设计辅助矩阵P3、辅助向量Q3、W3如下：
式中，l3＝1，神经网络权重通过如下自适应律获得：
其中，γ3＝2，Γ3为对角元素为0.05的对角矩阵。
步骤四、给定速度指令定义Vd，定义速度跟踪误差：
设计节流阀开度如下：
式中，kpV＝0.5、kiV＝0.001、kdV＝0.01。
步骤五、根据得到的舵偏角δe和节流阀开度β，返回到高超声速飞行器的动力学模型，对高度和速度进行跟踪控制。
本发明未详细说明部分属于领域技术人员公知常识。

Claims (1)

1.一种参数学习有限时间收敛的高超声速飞行器控制方法，其特征在于包括以下步骤：
步骤一、建立高超声速飞行器纵向通道动力学模型：
<mrow> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mo>=</mo> <mfrac> <mrow> <mi>T</mi> <mi> </mi> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mi>&amp;alpha;</mi> <mo>-</mo> <mi>D</mi> </mrow> <mi>m</mi> </mfrac> <mo>-</mo> <mfrac> <mrow> <mi>&amp;mu;</mi> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mi>&amp;gamma;</mi> </mrow> <msup> <mi>r</mi> <mn>2</mn> </msup> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <mover> <mi>h</mi> <mo>&amp;CenterDot;</mo> </mover> <mo>=</mo> <mi>V</mi> <mi> </mi> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mi>&amp;gamma;</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <mover> <mi>&amp;gamma;</mi> <mo>&amp;CenterDot;</mo> </mover> <mo>=</mo> <mfrac> <mrow> <mi>L</mi> <mo>+</mo> <mi>T</mi> <mi> </mi> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mi>&amp;alpha;</mi> </mrow> <mrow> <mi>m</mi> <mi>V</mi> </mrow> </mfrac> <mo>-</mo> <mfrac> <mrow> <mi>&amp;mu;</mi> <mo>-</mo> <msup> <mi>V</mi> <mn>2</mn> </msup> <mi>r</mi> <mi> </mi> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mi>&amp;gamma;</mi> </mrow> <mrow> <msup> <mi>Vr</mi> <mn>2</mn> </msup> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <mover> <mi>&amp;alpha;</mi> <mo>&amp;CenterDot;</mo> </mover> <mo>=</mo> <mi>q</mi> <mo>-</mo> <mover> <mi>&amp;gamma;</mi> <mo>&amp;CenterDot;</mo> </mover> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <mover> <mi>q</mi> <mo>&amp;CenterDot;</mo> </mover> <mo>=</mo> <mfrac> <msub> <mi>M</mi> <mrow> <mi>y</mi> <mi>y</mi> </mrow> </msub> <msub> <mi>I</mi> <mrow> <mi>y</mi> <mi>y</mi> </mrow> </msub> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>
其中， 表示动压，ρ表示空气密度，CL、CD、CT、CM(α)、CM(q)、CM(δe)均为气动参数，表示平均气动弦长，S表示气动参考面积；V表示速度，γ表示航迹倾角，h表示高度，α表示攻角，q表示俯仰角速度，δe是舵偏角，β为节流阀开度；T、D、L和Myy分别代表推力、阻力、升力和俯仰转动力矩；m、Iyy、μ和r代表质量、俯仰轴的转动惯量、引力系数以及距地心的距离；
步骤二、定义高度跟踪误差eh＝h-hd，设计航迹角指令γd：
<mrow> <msub> <mi>&amp;gamma;</mi> <mi>d</mi> </msub> <mo>=</mo> <mi>a</mi> <mi>r</mi> <mi>c</mi> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mrow> <mo>(</mo> <mfrac> <mrow> <mo>-</mo> <msub> <mi>k</mi> <mi>h</mi> </msub> <msub> <mi>e</mi> <mi>h</mi> </msub> <mo>-</mo> <msub> <mi>k</mi> <mi>i</mi> </msub> <mrow> <mo>&amp;Integral;</mo> <mrow> <msub> <mi>e</mi> <mi>h</mi> </msub> <mi>d</mi> <mi>t</mi> </mrow> </mrow> <mo>+</mo> <msub> <mover> <mi>h</mi> <mo>&amp;CenterDot;</mo> </mover> <mi>d</mi> </msub> </mrow> <mi>V</mi> </mfrac> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>
式中，hd为高度指令，为高度指令的一阶微分，kh>0，ki>0；考虑巡航段航迹角变化小，航迹角指令的一阶微分取为零；
取x1＝γ，x2＝θp，x3＝q，θp＝α+γ代表俯仰角；将式(3)-(5)写成如下严格反馈形式：
<mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mover> <mi>x</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>1</mn> </msub> <mo>=</mo> <msub> <mi>f</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mi>g</mi> <mn>1</mn> </msub> <msub> <mi>x</mi> <mn>2</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mover> <mi>x</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>2</mn> </msub> <mo>=</mo> <msub> <mi>x</mi> <mn>3</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mover> <mi>x</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>3</mn> </msub> <mo>=</mo> <msub> <mi>f</mi> <mn>3</mn> </msub> <mo>+</mo> <msub> <mi>g</mi> <mn>3</mn> </msub> <msub> <mi>&amp;delta;</mi> <mi>e</mi> </msub> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>
其中，fi,i＝1,3为未知函数，gi,i＝1,3为根据高超声速飞行器模型得到的已知函数；
步骤三、定义航迹角跟踪误差：
e1＝x1-γd (8)
设计俯仰角虚拟控制量：
<mrow> <msub> <mi>x</mi> <mrow> <mn>2</mn> <mi>c</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <msub> <mi>g</mi> <mn>1</mn> </msub> </mfrac> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mover> <mi>f</mi> <mo>^</mo> </mover> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>k</mi> <mn>1</mn> </msub> <msub> <mi>e</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mover> <mi>&amp;gamma;</mi> <mo>&amp;CenterDot;</mo> </mover> <mi>d</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>9</mn> <mo>)</mo> </mrow> </mrow>
式中，k1>0，为由RBF神经网络得到的f1估计值，其中为神经网络最优权重向量的估计值，θ1为RBF函数向量；
定义设计一阶滤波器如下：
式中xf1、θf1为x1、θ1通过滤波器(10)后获得的信号，设计辅助矩阵P1、辅助向量Q1、W1如下：
式中，l1>0，神经网络权重通过如下自适应律获得：
<mrow> <msub> <mover> <mover> <mi>&amp;omega;</mi> <mo>^</mo> </mover> <mo>&amp;CenterDot;</mo> </mover> <mn>1</mn> </msub> <mo>=</mo> <mo>-</mo> <msub> <mi>&amp;gamma;</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>e</mi> <mn>1</mn> </msub> <msub> <mi>&amp;theta;</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>&amp;Gamma;</mi> <mn>1</mn> </msub> <mfrac> <mrow> <msubsup> <mi>P</mi> <mn>1</mn> <mi>T</mi> </msubsup> <msub> <mi>W</mi> <mn>1</mn> </msub> </mrow> <mrow> <mo>|</mo> <mo>|</mo> <msub> <mi>W</mi> <mn>1</mn> </msub> <mo>|</mo> <mo>|</mo> </mrow> </mfrac> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>12</mn> <mo>)</mo> </mrow> </mrow>
其中，γ1>0，Γ1为正定常数对角矩阵；
设计一阶滤波器如下：
<mrow> <msub> <mi>&amp;alpha;</mi> <mn>2</mn> </msub> <msub> <mover> <mi>x</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mn>2</mn> <mi>d</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>x</mi> <mrow> <mn>2</mn> <mi>d</mi> </mrow> </msub> <mo>=</mo> <msub> <mi>x</mi> <mrow> <mn>2</mn> <mi>c</mi> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>13</mn> <mo>)</mo> </mrow> </mrow>
式中x2d为x2c通过滤波器(13)后获得的信号，α2>0；
定义俯仰角跟踪误差：
e2＝x2-x2d (14)
设计俯仰角速率虚拟控制量：
<mrow> <msub> <mi>x</mi> <mrow> <mn>3</mn> <mi>c</mi> </mrow> </msub> <mo>=</mo> <mo>-</mo> <msub> <mi>k</mi> <mn>2</mn> </msub> <msub> <mi>e</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>g</mi> <mn>1</mn> </msub> <msub> <mi>e</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mover> <mi>x</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mn>2</mn> <mi>d</mi> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>15</mn> <mo>)</mo> </mrow> </mrow>
式中，k2>0；
设计一阶滤波器如下：
<mrow> <msub> <mi>&amp;alpha;</mi> <mn>3</mn> </msub> <msub> <mover> <mi>x</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mn>3</mn> <mi>d</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>x</mi> <mrow> <mn>3</mn> <mi>d</mi> </mrow> </msub> <mo>=</mo> <msub> <mi>x</mi> <mrow> <mn>3</mn> <mi>c</mi> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>16</mn> <mo>)</mo> </mrow> </mrow>
式中x3d为x3c通过滤波器(16)后获得的信号，α3>0；
定义俯仰角速率跟踪误差：
e3＝x3-x3d (17)
设计舵偏角控制律如下：
<mrow> <msub> <mi>&amp;delta;</mi> <mi>e</mi> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <msub> <mi>g</mi> <mn>3</mn> </msub> </mfrac> <mrow> <mo>(</mo> <mo>-</mo> <msub> <mover> <mi>f</mi> <mo>^</mo> </mover> <mn>3</mn> </msub> <mo>-</mo> <msub> <mi>k</mi> <mn>3</mn> </msub> <msub> <mi>e</mi> <mn>3</mn> </msub> <mo>-</mo> <msub> <mi>e</mi> <mn>2</mn> </msub> <mo>+</mo> <msub> <mover> <mi>x</mi> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mn>3</mn> <mi>d</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>18</mn> <mo>)</mo> </mrow> </mrow>
式中，k3>0，为由RBF神经网络得到的f3估计值，其中为神经网络最优权重向量的估计值，θ3为RBF函数向量；
令设计一阶滤波器如下：
式中xf3、θf3为x3、θ3通过滤波器(19)后获得的信号，设计辅助矩阵P3、辅助向量Q3、W3如下：
式中，l3>0，神经网络权重通过如下自适应律获得：
<mrow> <msub> <mover> <mover> <mi>&amp;omega;</mi> <mo>^</mo> </mover> <mo>&amp;CenterDot;</mo> </mover> <mn>3</mn> </msub> <mo>=</mo> <mo>-</mo> <msub> <mi>&amp;gamma;</mi> <mn>3</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>e</mi> <mn>3</mn> </msub> <msub> <mi>&amp;theta;</mi> <mn>3</mn> </msub> <mo>-</mo> <msub> <mi>&amp;Gamma;</mi> <mn>3</mn> </msub> <mfrac> <mrow> <msubsup> <mi>P</mi> <mi>3</mi> <mi>T</mi> </msubsup> <msub> <mi>W</mi> <mn>3</mn> </msub> </mrow> <mrow> <mo>|</mo> <mo>|</mo> <msub> <mi>W</mi> <mn>3</mn> </msub> <mo>|</mo> <mo>|</mo> </mrow> </mfrac> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mi>2</mi> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
其中，γ3>0，Γ3为正定常数对角矩阵；
步骤四、定义速度跟踪误差：
<mrow> <mover> <mi>V</mi> <mo>~</mo> </mover> <mo>=</mo> <msub> <mi>V</mi> <mi>d</mi> </msub> <mo>-</mo> <mi>V</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>22</mn> <mo>)</mo> </mrow> </mrow>
式中，Vd为速度指令；设计节流阀开度如下：
<mrow> <mi>&amp;beta;</mi> <mo>=</mo> <msub> <mi>k</mi> <mrow> <mi>p</mi> <mi>V</mi> </mrow> </msub> <mover> <mi>V</mi> <mo>~</mo> </mover> <mo>+</mo> <msub> <mi>k</mi> <mrow> <mi>i</mi> <mi>V</mi> </mrow> </msub> <mo>&amp;Integral;</mo> <mover> <mi>V</mi> <mo>~</mo> </mover> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mi>d</mi> <mi>t</mi> <mo>+</mo> <msub> <mi>k</mi> <mrow> <mi>d</mi> <mi>V</mi> </mrow> </msub> <mfrac> <mrow> <mi>d</mi> <mover> <mi>V</mi> <mo>~</mo> </mover> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>23</mn> <mo>)</mo> </mrow> </mrow>
式中，kpV>0、kiV>0、kdV>0；
步骤五、根据得到的舵偏角δe和节流阀开度β，返回到高超声速飞行器的动力学模型(1)-(5)，对高度和速度进行跟踪控制。
CN201710789197.5A 2017-09-05 2017-09-05 参数学习有限时间收敛的高超声速飞行器控制方法 Active CN107390531B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710789197.5A CN107390531B (zh) 2017-09-05 2017-09-05 参数学习有限时间收敛的高超声速飞行器控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710789197.5A CN107390531B (zh) 2017-09-05 2017-09-05 参数学习有限时间收敛的高超声速飞行器控制方法

Publications (2)

Publication Number Publication Date
CN107390531A true CN107390531A (zh) 2017-11-24
CN107390531B CN107390531B (zh) 2019-10-18

Family

ID=60348028

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710789197.5A Active CN107390531B (zh) 2017-09-05 2017-09-05 参数学习有限时间收敛的高超声速飞行器控制方法

Country Status (1)

Country Link
CN (1) CN107390531B (zh)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108303889A (zh) * 2018-02-07 2018-07-20 中国航空工业集团公司西安飞机设计研究所 一种基于非线性信息的时标分离飞行器弹性体控制方法
CN108333939A (zh) * 2018-02-07 2018-07-27 中国航空工业集团公司西安飞机设计研究所 一种基于神经网络的时标分离飞行器弹性体智能控制方法
CN108415247A (zh) * 2018-02-07 2018-08-17 中国航空工业集团公司西安飞机设计研究所 一种基于标称信息的时标分离飞行器弹性体鲁棒控制方法
CN110320807A (zh) * 2019-07-24 2019-10-11 西北工业大学 基于奇异摄动分解的弹性飞行器数据筛选自适应控制方法
CN110456642A (zh) * 2019-07-24 2019-11-15 西北工业大学 基于奇异摄动分析的弹性飞行器鲁棒有限时间控制方法
CN111487866A (zh) * 2020-04-09 2020-08-04 中北大学 一种基于混合事件触发机制的高超声速飞行器神经抗干扰控制方法
CN113985901A (zh) * 2021-09-14 2022-01-28 中国人民解放军海军工程大学 基于扰动估计的高超声速飞行器预设性能控制方法及装置
CN113985732A (zh) * 2021-10-11 2022-01-28 广州大学 针对飞行器系统的自适应神经网络控制方法及装置
CN116594414A (zh) * 2023-03-28 2023-08-15 西北工业大学 一种高超声速飞行器纵向控制方法

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050137724A1 (en) * 2003-10-10 2005-06-23 Georgia Tech Research Corporation Adaptive observer and related method
GB2423377A (en) * 2002-12-09 2006-08-23 Georgia Tech Res Inst Adaptive output feedback apparatuses and methods capable of controlling a non-minimum phase system.
US8285659B1 (en) * 2009-08-18 2012-10-09 The United States of America as represented by the Administrator of the National Aeronautics & Space Administration (NASA) Aircraft system modeling error and control error
CN102880052A (zh) * 2012-09-29 2013-01-16 西北工业大学 基于时标功能分解的高超声速飞行器执行器饱和控制方法
CN102880055A (zh) * 2012-09-29 2013-01-16 西北工业大学 基于预测模型的高超声速飞行器神经网络控制方法
CN102929141A (zh) * 2012-10-10 2013-02-13 西北工业大学 飞行器时间滞后时变模型逼近及控制器设计方法
CN103365296A (zh) * 2013-06-29 2013-10-23 天津大学 一种四旋翼无人飞行器非线性输出反馈飞行控制方法
CN104022742A (zh) * 2014-06-19 2014-09-03 金陵科技学院 基于神经网络观测器的飞行器姿态鲁棒反演容错控制方法
CN104765272A (zh) * 2014-03-05 2015-07-08 北京航空航天大学 一种基于pid神经元网络控制(pidnn)的四旋翼飞行器控制方法
US9146557B1 (en) * 2014-04-23 2015-09-29 King Fahd University Of Petroleum And Minerals Adaptive control method for unmanned vehicle with slung load

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2423377A (en) * 2002-12-09 2006-08-23 Georgia Tech Res Inst Adaptive output feedback apparatuses and methods capable of controlling a non-minimum phase system.
US20050137724A1 (en) * 2003-10-10 2005-06-23 Georgia Tech Research Corporation Adaptive observer and related method
US8285659B1 (en) * 2009-08-18 2012-10-09 The United States of America as represented by the Administrator of the National Aeronautics & Space Administration (NASA) Aircraft system modeling error and control error
CN102880052A (zh) * 2012-09-29 2013-01-16 西北工业大学 基于时标功能分解的高超声速飞行器执行器饱和控制方法
CN102880055A (zh) * 2012-09-29 2013-01-16 西北工业大学 基于预测模型的高超声速飞行器神经网络控制方法
CN102929141A (zh) * 2012-10-10 2013-02-13 西北工业大学 飞行器时间滞后时变模型逼近及控制器设计方法
CN103365296A (zh) * 2013-06-29 2013-10-23 天津大学 一种四旋翼无人飞行器非线性输出反馈飞行控制方法
CN104765272A (zh) * 2014-03-05 2015-07-08 北京航空航天大学 一种基于pid神经元网络控制(pidnn)的四旋翼飞行器控制方法
US9146557B1 (en) * 2014-04-23 2015-09-29 King Fahd University Of Petroleum And Minerals Adaptive control method for unmanned vehicle with slung load
CN104022742A (zh) * 2014-06-19 2014-09-03 金陵科技学院 基于神经网络观测器的飞行器姿态鲁棒反演容错控制方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
BIN XU: "《Neural Control for Longitudinal Dynamics of Hypersonic Aircraft》", 《2013 INTERNATIONAL CONFERENCE ON UNMANNED AIRCRAFT SYSTEMS (ICUAS)》 *
JIANMIN WANG: "《RBF Neural Network based Adaptive Sliding Mode Control forHypersonic Flight Vehicles》", 《PROCEEDINGS OF 2016 IEEE CHINESE GUIDANCE, NAVIGATION AND CONTROL CONFERENCE》 *
李砚浓: "《基于RBF 神经网络自适应PID四旋翼飞行器控制》", 《控制工程》 *
王丽: "《基于混合神经网络的鲁棒自适应飞行控制器的设计》", 《电光与控制》 *

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108415247B (zh) * 2018-02-07 2019-12-20 中国航空工业集团公司西安飞机设计研究所 一种基于标称信息的时标分离飞行器弹性体鲁棒控制方法
CN108333939A (zh) * 2018-02-07 2018-07-27 中国航空工业集团公司西安飞机设计研究所 一种基于神经网络的时标分离飞行器弹性体智能控制方法
CN108415247A (zh) * 2018-02-07 2018-08-17 中国航空工业集团公司西安飞机设计研究所 一种基于标称信息的时标分离飞行器弹性体鲁棒控制方法
CN108303889A (zh) * 2018-02-07 2018-07-20 中国航空工业集团公司西安飞机设计研究所 一种基于非线性信息的时标分离飞行器弹性体控制方法
CN110320807B (zh) * 2019-07-24 2022-05-17 西北工业大学 基于奇异摄动分解的弹性飞行器数据筛选自适应控制方法
CN110456642A (zh) * 2019-07-24 2019-11-15 西北工业大学 基于奇异摄动分析的弹性飞行器鲁棒有限时间控制方法
CN110320807A (zh) * 2019-07-24 2019-10-11 西北工业大学 基于奇异摄动分解的弹性飞行器数据筛选自适应控制方法
CN111487866A (zh) * 2020-04-09 2020-08-04 中北大学 一种基于混合事件触发机制的高超声速飞行器神经抗干扰控制方法
CN113985901A (zh) * 2021-09-14 2022-01-28 中国人民解放军海军工程大学 基于扰动估计的高超声速飞行器预设性能控制方法及装置
CN113985901B (zh) * 2021-09-14 2023-06-27 中国人民解放军海军工程大学 基于扰动估计的高超声速飞行器预设性能控制方法及装置
CN113985732A (zh) * 2021-10-11 2022-01-28 广州大学 针对飞行器系统的自适应神经网络控制方法及装置
CN113985732B (zh) * 2021-10-11 2023-06-20 广州大学 针对飞行器系统的自适应神经网络控制方法及装置
CN116594414A (zh) * 2023-03-28 2023-08-15 西北工业大学 一种高超声速飞行器纵向控制方法
CN116594414B (zh) * 2023-03-28 2024-03-15 西北工业大学 一种高超声速飞行器纵向控制方法

Also Published As

Publication number Publication date
CN107390531B (zh) 2019-10-18

Similar Documents

Publication Publication Date Title
CN107390531B (zh) 参数学习有限时间收敛的高超声速飞行器控制方法
CN107450324B (zh) 考虑攻角约束的高超声速飞行器自适应容错控制方法
CN109144084B (zh) 一种基于固定时间收敛观测器的垂直起降重复使用运载器姿态跟踪控制方法
CN107479383B (zh) 基于鲁棒设计的高超声速飞行器神经网络复合学习控制方法
CN105676641A (zh) 针对四旋翼无人机非线性模型的基于反步和滑模控制技术的非线性鲁棒控制器的设计方法
CN102880052B (zh) 基于时标功能分解的高超声速飞行器执行器饱和控制方法
CN108828957B (zh) 基于切换机制的飞行器全局有限时间神经网络控制方法
CN107942651A (zh) 一种近空间飞行器控制系统
CN102880055B (zh) 基于预测模型的高超声速飞行器神经网络控制方法
CN109426146A (zh) 高超声速飞行器的高阶非奇异Terminal滑模控制方法
CN102749851A (zh) 一种挠性高超声速飞行器的精细抗干扰跟踪控制器
CN109703768B (zh) 一种基于姿态/轨迹复合控制的软式空中加油对接方法
CN102866635B (zh) 基于等价模型的高超声速飞行器离散神经网络自适应控制方法
CN109828602B (zh) 一种基于观测补偿技术的航迹回路非线性模型变换方法
CN106896722B (zh) 采用状态反馈与神经网络的高超飞行器复合控制方法
CN103853157A (zh) 一种基于自适应滑模的飞行器姿态控制方法
CN113268084B (zh) 一种无人机编队智能容错控制方法
CN105759609B (zh) 一种基于显式非线性模型预测控制的舰载机自主着舰方法
CN105629734A (zh) 一种近空间飞行器的轨迹跟踪控制方法
CN102540882A (zh) 一种基于最小参数学习法的飞行器航迹倾角控制方法
CN110320927A (zh) 智能变形飞行器的飞行控制方法及系统
CN106527128B (zh) 兼顾瞬态响应与鲁棒稳定性的飞行控制律设计方法
CN110568765A (zh) 面向攻角跟踪的高超声速飞行器非对称输出受限控制方法
CN103592847A (zh) 一种基于高增益观测器的高超声速飞行器非线性控制方法
CN105843080A (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