CN107065544A - 基于攻角幂函数的高超飞行器神经网络控制方法 - Google Patents
基于攻角幂函数的高超飞行器神经网络控制方法 Download PDFInfo
- Publication number
- CN107065544A CN107065544A CN201710199103.9A CN201710199103A CN107065544A CN 107065544 A CN107065544 A CN 107065544A CN 201710199103 A CN201710199103 A CN 201710199103A CN 107065544 A CN107065544 A CN 107065544A
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- angle
- attack
- mover
- 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
技术领域
本发明属于高超声速飞行器控制技术领域，涉及一种基于攻角幂函数的高超飞行器神经网络控制方法。
背景技术
作为新一代的天地往返航空航天飞行系统，高超声速飞行器是未来进入空间、控制空间，确保空间优势核心能力的关键支柱，也是进行大规模空间开发的前提，因此越来越受到世界各国的高度重视。
目前除了传统的控制方法应用高超声速控制外，神经网络等智能控制方法由于具有良好的不确定性处理能力而被应用于高超声速飞行器。攻角稳定跟踪控制是高超声速控制中最基本的任务之一，它的重要性在于攻角的稳定有利于发动机与速度的控制，避免发动机推力与攻角之间的自激振荡。对传统的低速飞行器来说，攻角的稳定控制是比较容易的。由于高超声速飞行器攻角具有静不稳定的特点，攻角跟踪与稳定对高超声速飞行器来说，并不是一件容易的事情。
高超声速飞行器模型具有快时变特点，而且高超声速飞行器气动参数具有强不确定性，因此本发明提出了一类变结构控制与神经网络相结合的复合控制方法，主要是因为变结构控制具有较好的快速性与鲁棒性。同时神经网络自适应控制能有效补偿模型不确定性与未建模因素带来的干扰。因此本发明方法不仅具有很好的理论创新性，而且具有一定的工程应用价值。
发明内容
为实现上述目的，本发明提供一种基于攻角幂函数的高超飞行器神经网络控制方法，得到的攻角响应具有较好的快速性特点，同时也具有很强的鲁棒性。
本发明所采用的技术方案是，基于攻角幂函数的高超飞行器神经网络控制方法，按照以下步骤进行：
步骤一，构造攻角测量、俯仰角速率的测量与滑模面信号；
步骤二，构造高超声速飞行器滑模控制律；
步骤三，构造自适应神经网络控制律；
步骤四，构造复合控制律经过控制分配，分别输入给高超声速飞行器的鸭翼与升降翼，以控制高超声速飞行器实现俯仰通道的攻角跟踪控制。
进一步的，所述步骤一，具体按照以下步骤进行：
首先采用攻角传感器，测量高超声速飞行器的攻角，记为α；采用速率陀螺测量高超声速飞行器的俯仰角速率，记为q；
其次，假定期望的攻角信号为αd，在飞行器控制计算机中生成攻角误差信号，记为eα，其满足eα＝α-αd；
再次，采用飞行器控制计算机生成攻角误差积分信号，记作Se，其满足Se＝∫(α-αd)dt；
最后，采用上述攻角误差信号与攻角误差积分信号组合生成滑模面信号，记作Sα，表达式为：Sα＝eα+c1Se；其中c1为正的控制参数，选取为c1＝0.1。
进一步的，所述步骤二，具体按照以下步骤进行：
构造高超声速飞行器滑模控制律u1：
u1＝-a1-c1eα-k1Sα-k2∫Sαdt-fs(Sα)，
其中，a1用于抵消飞行器的角速度与重力加速度项影响，按照下式计算：q为高超声速飞行器的俯仰角速率，g为重力加速度，γ为飞行器的航迹角，V为飞行器的速度；k1、k2为待设计的正参数，选取为k1＝5、k2＝0.2，
fs(Sα)为滑模控制项，定义如下：
其中k3、k4与k5为待设计的正参数，选取为k3＝0.5、k4＝0.3、k5＝0.02；
ε的含义是柔化因子，选取ε＝0.5；
τ的含义是柔化时间因子，选取τ＝0.5。
进一步的，所述步骤三，具体按照以下步骤进行：
根据测量的攻角信息，设计神经网络自适应控制律u2：
其中，φ为高超声速飞行器发动机的供油因子；为神经网络的权值因子，是通过自适应规律自动调节变化的；
其中为神经网络的权值， 为神经网络的权值，R1×n的表示是一个n维向量空间；B1(α)与B2(α)为基函数，B1(α)＝[1 α1α2 … αn-1]T，B2(α)＝[1 α1 α2 … αn-1]T，T的含义是矩阵的转置，表示是列向量；
神经网络的权值调节规律如下：
其中τ1、τ2为待调节的正参数，选取为τ1＝0.1，τ2＝0.2，的含义是 的导数，也就是两者的自适应调节规律的依据。
进一步的，所述步骤四，具体按照以下步骤进行：
根据上述高超声速飞行器滑模控制律u1与神经网络自适应控制规律u2组成如下复合控制律u：u＝u1+u2，然后按照均匀分配法进行控制律的分配，即定义：
从而有
其中m为飞行器的质量，S为飞行器的前向推力特征面积，为动压头数据，ρ为大气密度；δe为高超声速飞行器俯仰舵的舵偏值，δc为高超声速飞行器鸭翼舵的舵偏值。
本发明的有益效果是：通过测量高超声速飞行器的攻角以及俯仰角速度，通过飞行器控制计算机与给定的期望攻角指令进行对比形成攻角误差，一方面由攻角误差信号与攻角误差积分信号组成滑模面，再设计变结构控制律；另一方面由攻角误差信号设计FLNN型神经网络与其权值自适应调节规律；两个控制信号叠加后形成的复合信号经过控制分配，最终分别输入给两个鸭翼与升降翼，以控制高超声速飞行器实现俯仰通道的攻角跟踪控制。由于变结构控制与神经网络自适应策略的采用，本发明提供方法得到的攻角响应具有较好的快速性特点，同时也具有很强的鲁棒性，因此本发明方法不仅有较大的创新性，而且具有较大的理论价值和工程价值。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明提供的基于攻角幂函数的高超飞行器神经网络控制方法原理图。
图2是本发明实施例提供的2度攻角跟踪情况下的高超声速飞行器攻角曲线。
图3是本发明实施例提供的2度攻角跟踪情况下的高超声速飞行器高度曲线。
图4是本发明实施例提供的2度攻角跟踪情况下的高超声速飞行器速度曲线。
图5是本发明实施例提供的2度攻角跟踪情况下的高超声速飞行器航迹角曲线。
图6是本发明实施例提供的2度攻角跟踪情况下的高超声速飞行器供油因子曲线。
图7是本发明实施例提供的2度攻角跟踪情况下的高超声速飞行器鸭翼曲线。
图8是本发明实施例提供的2度攻角跟踪情况下的高超声速飞行器升降舵曲线。
图9是本发明实施例提供的2度攻角跟踪情况下的高超声速飞行器姿态角曲线。
具体实施方式
下面将结合本发明实施例中，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明一种基于攻角幂函数的高超飞行器神经网络控制方法，其思路是：通过测量高超声速飞行器的攻角、俯仰角速度信号，并与期望攻角信号组成变结构控制的滑模面，设计变结构控制器，利用变结构控制器的良好快速性与鲁棒性来处理高超声速飞行器快时变特点；针对高超声速飞行器气动参数的强不确定性，采用了一类以幂函数为基函数的神经网络结构，设计神经网络权值的自适应调节规律，最终组成高超声速飞行器神经网络与变结构的复合控制器，实现对期望攻角信号的跟踪。
具体来说，按照以下步骤进行：
步骤一：攻角测量、俯仰角速率的测量与滑模面的构造；
首先采用攻角传感器，测量高超声速飞行器的攻角，记为α；采用速率陀螺测量高超声速飞行器的俯仰角速率，记为q；
其次，假定期望的攻角信号为αd，在飞行器控制计算机中生成攻角误差信号，记为eα，其满足eα＝α-αd；
再次，采用飞行器控制计算机生成攻角误差积分信号，记作Se，其满足Se＝∫(α-αd)dt；
最后，采用上述攻角误差信号与攻角误差积分信号组合生成滑模面信号，记作Sα，其表达式如下Sα＝eα+c1Se；
其中c1为正的控制参数，其可初步选取为c1＝0.1；
步骤二：滑模控制律的构造；
按照下式构造高超声速飞行器滑模控制律u1：
u1＝-a1-c1eα-k1Sα-k2∫Sαdt-fs(Sα)
其中，a1没有物理含义，目的是为了抵消飞行器的角速度与重力加速度项影响，按照下面的公式计算即可：
q为高超声速飞行器的俯仰角速率，g为重力加速度，γ为飞行器的航迹角，V为飞行器的速度；k1、k2为待设计的正参数，初步选取为k1＝5、k2＝0.2，fs(Sα)定义如下：
其中k3、k4与k5为待设计的正参数，初步选取为k3＝0.5、k4＝0.3、k5＝0.02；
fs(Sα)为滑模控制项，由上述三小项组成，主要是为了提高具有切换性质以及柔化函数性质的控制量，使得滑模控制量具有很强的控制力，同时又能采用柔化函数的性质来减弱颤振；
ε的含义是柔化因子，可选取为小的正常数，如ε＝0.5；
τ的含义是柔化时间因子，可选取为小的正常数，如τ＝0.5；
步骤三：自适应神经网络控制律的构造；
根据测量的攻角信息，设计如下的神经网络自适应控制律u2：
其中，φ为高超声速飞行器发动机的供油因子，属于发动机的设计，现有技术即可得到；为神经网络的权值因子，其定义为下面矩阵，是通过后面自适应规律自动调节变化的；其中 为神经网络的权值，为神经网络的权值，R1×n的表示是一个n维向量空间；两者形式基本相同，调节规律不同；B1(α)与B2(α)为基函数，B1(α)＝[1 α1 α2 … αn-1]T，B2(α)＝[1 α1 α2 … αn-1]T，T的含义是矩阵的转置，表示是列向量；
神经网络的权值调节规律如下：
其中τ1、τ2为待调节的正参数，可初步选取为τ1＝0.1，τ2＝0.2。的含义是的导数，也就是两者的自适应调节规律的依据。
步骤四：复合控制律的构造与控制分配；
根据上述高超声速飞行器滑模控制律u1与神经网络自适应控制规律u2组成如下复合控制律u：
u＝u1+u2，
然后按照均匀分配法进行控制律的分配，即定义：
从而有
其中m为飞行器的质量，S为飞行器的前向推力特征面积，为动压头数据，ρ为大气密度；
δe为高超声速飞行器俯仰舵的舵偏值，δc为高超声速飞行器鸭翼舵的舵偏值。
步骤5，为了确保上述步骤一至步骤四中控制器的参数选取合理，可用通过计算机模拟仿真的手段进行编程，从而模拟被控对象高超声速飞行器的运动特性，从而方便进行参数调整。在此以某一类高超声速飞行器模型为例示范说明，其俯仰通道可以采用如下9个微分方程建模表示：
第一个微分方程为
其中V为飞行器的速度，为飞行器速度的导数，T为飞行器推力，α为飞行器的攻角，m为飞行器的质量，γ为飞行器航迹角，g为重力加速度，D为高超声速飞行器阻力。飞行器推力按照如下方式计算：
其中，φ为发动机供油因子，
表示两个矩阵相乘，其中 为实验测量值)，而η1、η2与η3为飞行器的弹性模态(其状态无需测量，需要测量的是系数，而且该测量是风洞实验所研究问题，本文不讨论该方法，用其结果得到的模型，作为本方法是否有效的验证工具)；为动压头数据，ρ为大气密度，S为飞行器的前向推力特征面积。CTφ也是风洞实验测量出来的，本文不讨论其测量方法，仅用其结果得到飞行器
模型，来验证本控制方法的有效性，其计算方式如下：
上式各参数主要是和模型相关，由实验所得，本专业人员非常熟悉其含义。
空气动力系数CT计算如下：
公式中各符号含义就是模型的气动参数；
CD是模型的气动参数，CD计算方式如下：
其中δe与δc分别为俯仰舵与鸭翼舵的舵偏值；其它符号含义是模型的气动参数；
第二个微分方程为：
其中，为飞行器的航迹角速度，γ为飞行器的航迹角，L为飞行器升力，其计算方式如下：
其中CL是模型中和升力相关的气动参数；
上式中各符号含义是模型中和升力相关的气动参数；
第三个微分方程为：
其中M为俯仰力矩，I为转动惯量，为俯仰角加速度，q为俯仰角。俯仰力矩的计算方法如下：
其中zT＝8.36*0.3048，(二者均为飞行器的尺寸参数，几何设计上有)，CM的计算方式如下：
所有符号含义是模型中和俯仰力矩相关的气动参数；
第四个微分方程为
是攻角的导数；
第五个微分方程为
其中，为高度的导数，h为飞行器高度。
第六个微分方程为
其中ωn＝5,(二者含义是发动机的供油设计，和本专利无关，是发动机设计所固有的方法)，φc为发动机的供油规律(为发动机的供油设计，和本专利无关，是发动机设计所固有的方法，以验证本控制方法的正确性)；
的含义是发动机的供油因子的导数；
第(7)(8)(9)个微分方程为弹性模态微分方程：
含义是弹性形态的二阶导数；
εm含义是发动机弹性形态建模所得到的振动弹性阻尼系数，属于建模，和本专利无关，用该建立好的模型，来验证控制方案的有效性；
ωmi含义是发动机弹性形态建模所得到的振动弹性频率系数，属于建模，和本专利无关，用该建立好的模型，来验证控制方案的有效性；
含义是弹性形态的一阶导数；
ηi含义是弹性形态的0阶导数，也就是弹性形态本身；
Ni的含义是弹性形态变化因子，按照下面公式计算；
其中i＝1,2,3。 含义是发动机弹性形态建模所得到的弹性相关系数，属于建模，和本专利无关，用该建立好的模型，来验证控制方案的有效性；
其计算方式如下：
其中
以上公式所有符号含义是发动机弹性形态建模所得到的弹性相关系数，属于建模，和本专利无关，用该建立好的模型，来验证控制方案的有效性；
其中，η＝[η1 η2 η3]，其它类似；
选取高超声速飞行器的初始条件为
V0＝7846*0.3048，α0＝0.0174，γ0＝0，q0＝0，h0＝85000*0.3048(初始高度)，
η10＝0.4588*0.3048*14.59，表示弹性形态η1在0时刻的初始取值η10，以下其它类似；
η20＝-0.08726*0.3048*14.59，
η30＝-0.03671*0.3048*14.59；
至此高超声速飞行器的模型描述完成。
步骤六：将步骤一至步骤四所得高超声速飞行器俯仰舵的舵偏值与高超声速飞行器鸭翼舵的舵偏值，代入步骤五所建立的高超声速飞行器模型，通过不断调整控制参数，并观察飞行器各状态的数据并画图，尤其是观测高超声速飞行器攻角变化的数据曲线，分析飞行器攻角跟踪的动态特性，从而确定最终高超声速飞行器攻角跟踪的控制方案参数，使得高超声速飞行器攻角响应具有满意的动态响应与稳态响应性能。
实施例及计算机仿真模拟结果分析
按照上述所设计控制律进行仿真，选取初始高度为h0＝85000*0.3048，初始速度为V0＝7846*0.3048，初始攻角α0＝0.0174，其它初始状态如下：γ0＝0，q0＝0，η10＝0.4588*0.3048*14.59，η20＝-0.08726*0.3048*14.59，η30＝-0.03671*0.3048*14.59(三个含义表示弹性形态η1在0时刻的初始取值η10，其它类似)，设定发动机的油门开度为PID控制律与柔化控制律的切换控制律，选取如果ev＞0，则设计如果ev＜0，则选取φc＝-300，其中ev＝Vd-V，ket、kset、kde、ketb与εev为可调正数。
选取期望速度为Vd＝2391，选取控制参数如下:c1＝0.1，k1＝5、k2＝0.2，k3＝0.5，k4＝0.3，k5＝0.02，τ1＝0.1，τ2＝0.2，将步骤四所得到的综合控制律代入模型五进行仿真，得到仿真结果图2至图9所示。
通过以上仿真结果与曲线可以看出，由于变结构的快速响应使得攻角输出在初始段具有较弱的阻尼而使得超调较大，这使得该控制策略比较适用于攻角控制精度要求不高但需要飞行器机动性反应速度较快的场合。而其中后面的波动是由于发动机的速度控制策略引起的，可以从发动机供油曲线比照来映证。从以上案例仿真结果可以看出，本发明提供的攻角跟踪控制方法能够有效地稳定高超声速飞行器攻角，因此本发明具有很好的实用价值，同时也能推广应用于其它类飞行器中。
以上所述仅为本发明的较佳实施例而已，并非用于限定本发明的保护范围。凡在本发明的精神和原则之内所作的任何修改、等同替换、改进等，均包含在本发明的保护范围内。
Claims (5)
1.基于攻角幂函数的高超飞行器神经网络控制方法，其特征在于，按照以下步骤进行：
步骤一，构造攻角测量、俯仰角速率的测量与滑模面信号；
步骤二，构造高超声速飞行器滑模控制律；
步骤三，构造自适应神经网络控制律；
步骤四，构造复合控制律经过控制分配，分别输入给高超声速飞行器的鸭翼与升降翼，以控制高超声速飞行器实现俯仰通道的攻角跟踪控制。
2.根据权利要求1所述的基于攻角幂函数的高超飞行器神经网络控制方法，其特征在于，所述步骤一，具体按照以下步骤进行：
首先采用攻角传感器，测量高超声速飞行器的攻角，记为α；采用速率陀螺测量高超声速飞行器的俯仰角速率，记为q；
其次，假定期望的攻角信号为αd，在飞行器控制计算机中生成攻角误差信号，记为eα，其满足eα＝α-αd；
再次，采用飞行器控制计算机生成攻角误差积分信号，记作Se，其满足Se＝∫(α-αd)dt；
最后，采用上述攻角误差信号与攻角误差积分信号组合生成滑模面信号，记作Sα，表达式为：Sα＝eα+c1Se；其中c1为正的控制参数，选取为c1＝0.1。
3.根据权利要求2所述的基于攻角幂函数的高超飞行器神经网络控制方法，其特征在于，所述步骤二，具体按照以下步骤进行：
构造高超声速飞行器滑模控制律u1：
u1＝-a1-c1eα-k1Sα-k2∫Sαdt-fs(Sα)，
其中，a1用于抵消飞行器的角速度与重力加速度项影响，按照下式计算：
q为高超声速飞行器的俯仰角速率，g为重力加速度，γ为飞行器的航迹角，V为飞行器的速度；k1、k2为待设计的正参数，选取为k1＝5、k2＝0.2，
fs(Sα)为滑模控制项，定义如下：
<mrow>
<msub>
<mi>f</mi>
<mi>s</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>S</mi>
<mi>&alpha;</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>k</mi>
<mn>3</mn>
</msub>
<mfrac>
<msub>
<mi>S</mi>
<mi>&alpha;</mi>
</msub>
<mrow>
<mo>|</mo>
<msub>
<mi>S</mi>
<mi>&alpha;</mi>
</msub>
<mo>|</mo>
<mo>+</mo>
<mi>&epsiv;</mi>
</mrow>
</mfrac>
<mo>+</mo>
<msub>
<mi>k</mi>
<mn>4</mn>
</msub>
<mfrac>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>&tau;s</mi>
<mi>&alpha;</mi>
</msub>
</mrow>
</msup>
</mrow>
<mrow>
<mn>1</mn>
<mo>+</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>&tau;s</mi>
<mi>&alpha;</mi>
</msub>
</mrow>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
<msub>
<mi>k</mi>
<mn>5</mn>
</msub>
<mi>sgn</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>S</mi>
<mi>&alpha;</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
其中k3、k4与k5为待设计的正参数，选取为k3＝0.5、k4＝0.3、k5＝0.02；
ε的含义是柔化因子，选取ε＝0.5；
τ的含义是柔化时间因子，选取τ＝0.5。
4.根据权利要求3所述的基于攻角幂函数的高超飞行器神经网络控制方法，其特征在于，所述步骤三，具体按照以下步骤进行：
根据测量的攻角信息，设计神经网络自适应控制律u2：
<mrow>
<msub>
<mi>u</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mo>-</mo>
<msub>
<mover>
<mi>W</mi>
<mo>^</mo>
</mover>
<mi>a</mi>
</msub>
<msub>
<mi>B</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mover>
<mi>W</mi>
<mo>^</mo>
</mover>
<mi>b</mi>
</msub>
<msub>
<mi>B</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
<mi>&phi;</mi>
</mrow>
其中，φ为高超声速飞行器发动机的供油因子；为神经网络的权值因子，是通过自适应规律自动调节变化的；
其中为神经网络的权值， 为神经网络的权值，R1×n的表示是一个n维向量空间；B1(α)与B2(α)为基函数，B1(α)＝[1 α1α2 … αn-1]T，
B2(α)＝[1 α1 α2 … αn-1]T，T的含义是矩阵的转置，表示是列向量；
神经网络的权值调节规律如下：
<mrow>
<msub>
<mover>
<mover>
<mi>W</mi>
<mo>^</mo>
</mover>
<mo>&CenterDot;</mo>
</mover>
<mn>1</mn>
</msub>
<mo>=</mo>
<msub>
<mi>&tau;</mi>
<mn>1</mn>
</msub>
<msub>
<mi>S</mi>
<mi>&alpha;</mi>
</msub>
<msub>
<mi>B</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msub>
<mover>
<mover>
<mi>W</mi>
<mo>^</mo>
</mover>
<mo>&CenterDot;</mo>
</mover>
<mn>2</mn>
</msub>
<mo>=</mo>
<msub>
<mi>&tau;</mi>
<mn>2</mn>
</msub>
<msub>
<mi>S</mi>
<mi>&alpha;</mi>
</msub>
<msub>
<mi>B</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>&alpha;</mi>
<mo>)</mo>
</mrow>
<mi>&phi;</mi>
<mo>,</mo>
</mrow>
其中τ1、τ2为待调节的正参数，选取为τ1＝0.1，τ2＝0.2，的含义是 的导数，也就是两者的自适应调节规律的依据。
5.根据权利要求4所述的基于攻角幂函数的高超飞行器神经网络控制方法，其特征在于，所述步骤四，具体按照以下步骤进行：
根据上述高超声速飞行器滑模控制律u1与神经网络自适应控制规律u2组成如下复合控制律u：u＝u1+u2，然后按照均匀分配法进行控制律的分配，即定义：
<mrow>
<mfrac>
<mi>u</mi>
<mn>2</mn>
</mfrac>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mrow>
<mover>
<mi>q</mi>
<mo>&OverBar;</mo>
</mover>
<mi>S</mi>
</mrow>
<mrow>
<mi>m</mi>
<mi>V</mi>
</mrow>
</mfrac>
<mn>0.76224</mn>
<msub>
<mi>&delta;</mi>
<mi>e</mi>
</msub>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mrow>
<mover>
<mi>q</mi>
<mo>&OverBar;</mo>
</mover>
<mi>S</mi>
</mrow>
<mrow>
<mi>m</mi>
<mi>V</mi>
</mrow>
</mfrac>
<mn>0.7408</mn>
<msub>
<mi>&delta;</mi>
<mi>c</mi>
</msub>
</mrow>
从而有
<mrow>
<msub>
<mi>&delta;</mi>
<mi>e</mi>
</msub>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mrow>
<mi>u</mi>
<mi>m</mi>
<mi>V</mi>
</mrow>
<mrow>
<mn>2</mn>
<mover>
<mi>q</mi>
<mo>&OverBar;</mo>
</mover>
<mi>S</mi>
<mn>0.76224</mn>
</mrow>
</mfrac>
</mrow>
<mrow>
<msub>
<mi>&delta;</mi>
<mi>c</mi>
</msub>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mrow>
<mi>u</mi>
<mi>m</mi>
<mi>V</mi>
</mrow>
<mrow>
<mn>2</mn>
<mover>
<mi>q</mi>
<mo>&OverBar;</mo>
</mover>
<mi>S</mi>
<mn>0.7408</mn>
</mrow>
</mfrac>
</mrow>
其中m为飞行器的质量，S为飞行器的前向推力特征面积，为动压头数据，ρ为大气密度；δe为高超声速飞行器俯仰舵的舵偏值，δc为高超声速飞行器鸭翼舵的舵偏值。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710199103.9A CN107065544B (zh) | 2017-03-29 | 2017-03-29 | 基于攻角幂函数的高超飞行器神经网络控制方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710199103.9A CN107065544B (zh) | 2017-03-29 | 2017-03-29 | 基于攻角幂函数的高超飞行器神经网络控制方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107065544A true CN107065544A (zh) | 2017-08-18 |
| CN107065544B CN107065544B (zh) | 2019-12-06 |
Family
ID=59620734
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710199103.9A Expired - Fee Related CN107065544B (zh) | 2017-03-29 | 2017-03-29 | 基于攻角幂函数的高超飞行器神经网络控制方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107065544B (zh) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106896722A (zh) * | 2017-03-29 | 2017-06-27 | 郑州轻工业学院 | 采用状态反馈与神经网络的高超飞行器复合控制方法 |
| CN108646557A (zh) * | 2018-05-08 | 2018-10-12 | 中国人民解放军海军航空大学 | 一种基于跟踪微分与柔化函数的飞行器攻角跟踪控制方法 |
| CN110135561A (zh) * | 2019-04-29 | 2019-08-16 | 北京航天自动控制研究所 | 一种实时在线飞行器ai神经网络系统 |
| CN110889170A (zh) * | 2019-11-29 | 2020-03-17 | 中国人民解放军海军工程大学 | 飞行器大角度攻击目标的落角与攻角的估计方法 |
| CN111122901A (zh) * | 2019-12-31 | 2020-05-08 | 中国航空工业集团公司西安飞机设计研究所 | 一种攻角误差自动检测方法及其修正系统 |
| CN111273681A (zh) * | 2020-04-09 | 2020-06-12 | 中北大学 | 一种考虑攻角受限的高超声速飞行器高安全抗干扰控制方法 |
| CN116594414A (zh) * | 2023-03-28 | 2023-08-15 | 西北工业大学 | 一种高超声速飞行器纵向控制方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103425135A (zh) * | 2013-07-30 | 2013-12-04 | 南京航空航天大学 | 一种具有输入饱和的近空间飞行器鲁棒控制方法 |
| CN103809446A (zh) * | 2014-02-28 | 2014-05-21 | 西安费斯达自动化工程有限公司 | 飞行器多回路模型簇颤振抑制复合频率鲁棒控制器设计方法 |
| CN109426146A (zh) * | 2017-08-23 | 2019-03-05 | 中国空气动力研究与发展中心计算空气动力研究所 | 高超声速飞行器的高阶非奇异Terminal滑模控制方法 |
-
2017
- 2017-03-29 CN CN201710199103.9A patent/CN107065544B/zh not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103425135A (zh) * | 2013-07-30 | 2013-12-04 | 南京航空航天大学 | 一种具有输入饱和的近空间飞行器鲁棒控制方法 |
| CN103809446A (zh) * | 2014-02-28 | 2014-05-21 | 西安费斯达自动化工程有限公司 | 飞行器多回路模型簇颤振抑制复合频率鲁棒控制器设计方法 |
| CN109426146A (zh) * | 2017-08-23 | 2019-03-05 | 中国空气动力研究与发展中心计算空气动力研究所 | 高超声速飞行器的高阶非奇异Terminal滑模控制方法 |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106896722A (zh) * | 2017-03-29 | 2017-06-27 | 郑州轻工业学院 | 采用状态反馈与神经网络的高超飞行器复合控制方法 |
| CN106896722B (zh) * | 2017-03-29 | 2019-11-29 | 郑州轻工业学院 | 采用状态反馈与神经网络的高超飞行器复合控制方法 |
| CN108646557A (zh) * | 2018-05-08 | 2018-10-12 | 中国人民解放军海军航空大学 | 一种基于跟踪微分与柔化函数的飞行器攻角跟踪控制方法 |
| CN110135561A (zh) * | 2019-04-29 | 2019-08-16 | 北京航天自动控制研究所 | 一种实时在线飞行器ai神经网络系统 |
| CN110889170A (zh) * | 2019-11-29 | 2020-03-17 | 中国人民解放军海军工程大学 | 飞行器大角度攻击目标的落角与攻角的估计方法 |
| CN110889170B (zh) * | 2019-11-29 | 2022-09-30 | 中国人民解放军海军工程大学 | 飞行器大角度攻击目标的落角与攻角的估计方法 |
| CN111122901A (zh) * | 2019-12-31 | 2020-05-08 | 中国航空工业集团公司西安飞机设计研究所 | 一种攻角误差自动检测方法及其修正系统 |
| CN111273681A (zh) * | 2020-04-09 | 2020-06-12 | 中北大学 | 一种考虑攻角受限的高超声速飞行器高安全抗干扰控制方法 |
| CN116594414A (zh) * | 2023-03-28 | 2023-08-15 | 西北工业大学 | 一种高超声速飞行器纵向控制方法 |
| CN116594414B (zh) * | 2023-03-28 | 2024-03-15 | 西北工业大学 | 一种高超声速飞行器纵向控制方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107065544B (zh) | 2019-12-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107065544A (zh) | 基于攻角幂函数的高超飞行器神经网络控制方法 | |
| CN103488814B (zh) | 一种适用于再入飞行器姿态控制的闭环仿真系统 | |
| CN106896722A (zh) | 采用状态反馈与神经网络的高超飞行器复合控制方法 | |
| CN107608367A (zh) | 多变量干扰补偿四旋翼无人机轨迹与姿态协同控制方法 | |
| CN106527122B (zh) | 平流层飞艇定高飞行非线性pid控制方法 | |
| CN107450324A (zh) | 考虑攻角约束的高超声速飞行器自适应容错控制方法 | |
| CN105676641A (zh) | 针对四旋翼无人机非线性模型的基于反步和滑模控制技术的非线性鲁棒控制器的设计方法 | |
| CN105159305A (zh) | 一种基于滑模变结构的四旋翼飞行控制方法 | |
| CN104122794B (zh) | 微陀螺仪的自适应模糊神经补偿非奇异终端滑模控制方法 | |
| CN109062055A (zh) | 一种基于Back-stepping鲁棒自适应动态面的近空间飞行器控制系统 | |
| CN109270947A (zh) | 倾转旋翼无人机飞行控制系统 | |
| CN106444807A (zh) | 一种栅格舵与侧喷流的复合姿态控制方法 | |
| CN105629734A (zh) | 一种近空间飞行器的轨迹跟踪控制方法 | |
| CN107942651A (zh) | 一种近空间飞行器控制系统 | |
| CN104597911A (zh) | 空中加油受油机自适应最优对接轨迹跟踪飞行控制方法 | |
| CN107390531A (zh) | 参数学习有限时间收敛的高超声速飞行器控制方法 | |
| CN107085435A (zh) | 基于耦合分析的高超声速飞行器姿态协调控制方法 | |
| CN105807789A (zh) | 基于t-s模糊观测器补偿的uuv控制方法 | |
| CN106707749A (zh) | 一种针对仿生扑翼飞行机器人的控制方法 | |
| CN106814626A (zh) | 基于模糊变结构的超声速飞行器过载指令快速跟踪方法 | |
| CN106802570A (zh) | 一种无人直升机位置跟踪的方法与装置 | |
| CN107367941A (zh) | 基于非线性增益的高超声速飞行器攻角观测方法 | |
| CN104462810B (zh) | 一种适用于轮控卫星姿态机动与跟踪控制的sdre参数调节方法 | |
| CN113268084B (zh) | 一种无人机编队智能容错控制方法 | |
| CN107894775A (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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20191206 Termination date: 20210329 |