基于智能网络模型的飞行器预设性能滑模控制方法及装置

申请号 CN202410239277.3 申请日 2024-03-04 公开(公告)号 CN117826617B 公开(公告)日 2024-05-10
申请人 西北工业大学; 发明人 丁一波; 张天辰; 岳晓奎; 代洪华; 李勇; 张顺家; 孙军; 李勰; 刘传凯;
摘要 本 发明 公开一种基于智能网络模型的 飞行器 预设性能滑模控制方法及装置,涉及飞行器控制领域,方法包括:根据当前时刻飞行器的速度、 攻 角 、升降 舵 偏转角及燃油当量比,基于智能网络模型确定当前误差动 力 学函数;根据当前误差动力学函数及飞行器传统误差动力学模型建立飞行器智能误差动力学模型;确定速度设定时间预设性能函数及攻角设定时间预设性能函数,进一步确定速度 控制器 及攻角转换误差动力学模型;基于攻角转换误差动力学模型,采用自适应超螺旋滑模 算法 ,建立飞行器预设性能超螺旋滑模控制器,以确定攻角控制器;基于速度控制器及攻角控制器对飞行器进行控制。本发明提高了飞行器控制的鲁棒性和控制 精度 。
权利要求

1.一种基于智能网络模型的飞行器预设性能滑模控制方法,其特征在于,所述基于智能网络模型的飞行器预设性能滑模控制方法包括:
根据当前时刻飞行器的速度、、升降偏转角及燃油当量比,基于智能网络模型确定当前误差动学函数;所述智能网络模型为预先采用飞行器的历史飞行量测数据集对神经网络进行训练得到的;所述历史飞行量测数据集中包括历史设定时段内各时刻飞行器的速度、攻角、升降舵偏转角、燃油当量比及对应的误差动力学函数;所述当前误差动力学函数包括:当前第一速度误差动力学函数、当前第二速度误差动力学函数、当前攻角误差动力学函数、当前第一俯仰角速度误差动力学函数及当前第二俯仰角速度误差动力学函数;
根据所述当前误差动力学函数及飞行器传统误差动力学模型,建立飞行器智能误差动力学模型;所述飞行器智能误差动力学模型包括速度智能误差动力学模型及攻角智能误差动力学模型;
所述速度智能误差动力学模型为:

所述攻角智能误差动力学模型为:

其中, 为 对时间的一阶导数, 为速度误差, ,V为飞行器的速度,
Vd为速度指令, 为当前第一速度误差动力学函数, 为当前第二速度误差动力学函数, 为燃油当量比, 为速度集总扰动, 为Vd对时间的一阶导数, 为 对时间的一阶导数, 为攻角误差, , 为飞行器的攻角, 为攻角指令, 为
当前攻角误差动力学函数,Q为飞行器的俯仰角速度, 为攻角集总扰动, 为 对时间的一阶导数, 为Q对时间的一阶导数, 为当前第一俯仰角速度误差动力学函数,为当前第二俯仰角速度误差动力学函数, 为升降舵偏转角, 为俯仰角速度集总扰动;
确定速度设定时间预设性能函数及攻角设定时间预设性能函数;
基于所述速度设定时间预设性能函数、所述攻角设定时间预设性能函数及所述飞行器智能误差动力学模型,确定速度控制器及攻角转换误差动力学模型,具体包括:
基于所述速度设定时间预设性能函数,对速度误差进行限幅,以确定速度预设性能条件;
基于所述攻角设定时间预设性能函数,对攻角误差进行限幅,以确定攻角预设性能条件;
确定所述速度智能误差动力学模型的速度预设性能控制器,并基于所述速度预设性能条件,采用齐次高阶滑模观测器对速度集总扰动进行补偿,以确定速度控制器;
所述速度控制器为:

其中, 为速度集总扰动的估计值, 为速度设定时间预设性能函数, 为
对时间的一阶导数, 为速度正待设计参数, 为速度中间参数,
, 为t时刻的速
度误差, 为速度设定时间预设性能函数的下界系数, 为速度设定时间预设性能函数的上界系数, 为速度转换误差;
确定所述攻角智能误差动力学模型的攻角预设性能控制器,并基于所述攻角预设性能条件,采用齐次高阶滑模观测器对速度集总扰动及俯仰角速度集总扰动进行补偿,以确定虚拟俯仰角速度控制量;
根据所述虚拟俯仰角速度控制量,采用一阶滤波器确定滤波俯仰角速度输出信号
根据所述滤波俯仰角速度输出信号,确定攻角转换误差动力学模型;
所述攻角转换误差动力学模型为:

其中, 为 对时间的一阶导数, 为攻角转换误差, 为攻角中间参数,
, 为攻角设定
时间预设性能函数, 为 对时间的一阶导数, 为t时刻的攻角误差,
为攻角设定时间预设性能函数的下界系数, 为攻角设定时间预设性能函数的上界系数,为 对时间的一阶导数, 为俯仰角速度误差, ,Qd为滤波俯仰角速度
输出信号, 为滤波俯仰角速度输出信号与虚拟俯仰角速度控制量的差值,,Qc为虚拟俯仰角速度控制量, 为攻角集总扰动误差,
, 为攻角集总扰动的估计值, 为攻角正待设计参数;
基于所述攻角转换误差动力学模型,采用自适应超螺旋滑模算法,建立飞行器预设性能超螺旋滑模控制器,以确定攻角控制器;
所述攻角控制器为:


其中, 为攻角实际控制律,表示升降舵偏转角,a1和a2为常数,a1>a2,1基于所述速度控制器及所述攻角控制器对飞行器进行控制。
2.根据权利要求1所述的基于智能网络模型的飞行器预设性能滑模控制方法,其特征在于,所述速度设定时间预设性能函数为:

其中, 为速度设定时间预设性能函数, 为速度设定时间预设性能函数的第一系数, 为速度设定时间预设性能函数的第二系数, 为速度设定时间预设性能函数的第三系数,t为时刻, 为速度误差的初值, 为速度设定时间预设性能函数的变化方向, 为速度设定时间预设性能函数的稳态收敛值, 为速度误差收敛到的设定时间;
所述攻角设定时间预设性能函数为:

其中, 为攻角设定时间预设性能函数, 为攻角设定时间预设性能函数的第一系数, 为攻角设定时间预设性能函数的第二系数, 为攻角设定时间预设性能函数的第三系数, 为攻角误差的初值, 为攻角设定时间预设性能函数的变化方向,为攻角设定时间预设性能函数的稳态收敛值, 为攻角误差收敛到 的设定
时间。
3.根据权利要求1所述的基于智能网络模型的飞行器预设性能滑模控制方法,其特征在于,所述速度预设性能条件为:

其中, 为速度设定时间预设性能函数的下界系数, 为速度设定时间预设性能函数的上界系数, 为速度设定时间预设性能函数, 为t时刻的速度误差;
所述攻角预设性能条件为:

其中, 为攻角设定时间预设性能函数的下界系数, 为攻角设定时间预设性能函数的上界系数, 为攻角设定时间预设性能函数, 为t时刻的攻角误差。
4.根据权利要求1所述的基于智能网络模型的飞行器预设性能滑模控制方法,其特征在于,所述自适应超螺旋滑模算法包括非奇异快速终端滑模面及改进的自适应超螺旋趋近律;
所述非奇异快速终端滑模面为:

所述改进的自适应超螺旋趋近律为:

其中,s为滑模变量, 为s的导数, 为攻角转换误差, 为俯仰角速度误差, 、、k1、k2、k3均为正待设计参数,a1和a2为常数,a1>a2,15.一种计算机装置,包括:存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序以实现权利要求1‑4中任一项所述基于智能网络模型的飞行器预设性能滑模控制方法的步骤。

说明书全文

基于智能网络模型的飞行器预设性能滑模控制方法及装置

技术领域

[0001] 本发明涉及飞行器控制领域,特别是涉及一种基于智能网络模型的飞行器预设性能滑模控制方法及装置。

背景技术

[0002] 随着航天技术的迅速发展,高速飞行器的结构和任务需求越来越复杂,对控制系统的精度要求越来越高。由于高速飞行器飞行环境跨大空域宽速域,气动与气动热环境变化剧烈,飞行器动学模型复杂呈现强非线性,且模型参数呈现强不确定性。因此,如果可以提高飞行器建模精度,既可以增强系统的鲁棒性与自适应性,又可以提升系统的控制精度。
[0003] 传统飞行器动力学模型的结构参数和气动参数均为地面测试仿真后得到的拟合结果,在控制系统设计过程中均保持标称值不变。然而高速飞行器在飞行过程中受环境因素等诸多影响,模型参数呈现快时变性与强不确定性。传统控制器仅能保证飞行器飞行状态收敛,无法保证系统状态在设定时间按照预设状态收敛。而在实际飞行过程中,由于气推耦合与复杂环境因素的影响,需要考虑高速飞行器飞行过程中的瞬态与稳态性能,否则可能导致系统跟踪误差瞬时超调较大,产生过大,从而无法维持吸气式发动机进气量,导致飞行过程发动机熄火。

发明内容

[0004] 本发明的目的是提供一种基于智能网络模型的飞行器预设性能滑模控制方法及装置,可提高飞行器控制的鲁棒性和控制精度。
[0005] 为实现上述目的,本发明提供了一种基于智能网络模型的飞行器预设性能滑模控制方法,包括如下步骤。
[0006] 根据当前时刻飞行器的速度、攻角、升降偏转角及燃油当量比,基于智能网络模型确定当前误差动力学函数;所述智能网络模型为预先采用飞行器的历史飞行量测数据集对神经网络进行训练得到的;所述历史飞行量测数据集中包括历史设定时段内各时刻飞行器的速度、攻角、升降舵偏转角、燃油当量比及对应的误差动力学函数。
[0007] 根据所述当前误差动力学函数及飞行器传统误差动力学模型,建立飞行器智能误差动力学模型。
[0008] 确定速度设定时间预设性能函数及攻角设定时间预设性能函数。
[0009] 基于所述速度设定时间预设性能函数、所述攻角设定时间预设性能函数及所述飞行器智能误差动力学模型,确定速度控制器及攻角转换误差动力学模型。
[0010] 基于所述攻角转换误差动力学模型,采用自适应超螺旋滑模算法,建立飞行器预设性能超螺旋滑模控制器,以确定攻角控制器。
[0011] 基于所述速度控制器及所述攻角控制器对飞行器进行控制。
[0012] 为实现上述目的,本发明还提供了一种计算机装置,包括:存储器、处理器以存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序以实现上述基于智能网络模型的飞行器预设性能滑模控制方法的步骤。
[0013] 根据本发明提供的具体实施例,本发明公开了以下技术效果:使用智能网络模型对飞行器进行精确动力学建模,提升了飞行器智能误差动力学模型的精度,通过速度设定时间预设性能函数及攻角设定时间预设性能函数保证了飞行状态跟踪误差在设定时间收敛,提升了收敛速度与收敛精度,避免了初始时刻控制量过大的问题,针对攻角转换误差动力学模型,采用自适应超螺旋滑模算法建立飞行器预设性能超螺旋滑模控制器,以确定攻角控制器,避免了传统终端滑模的奇异问题并加快收敛速度,同时避免了控制增益过小导致控制发散与控制增益过大导致飞行器剧烈抖振的问题,提高了飞行器控制的鲁棒性和控制精度。附图说明
[0014] 为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
[0015] 图1为本发明提供的基于智能网络模型的飞行器预设性能滑模控制方法的流程图
[0016] 图2为智能网络模型训练集误差曲线图。
[0017] 图3为智能网络模型训练集均方误差曲线图。
[0018] 图4为飞行器的速度跟踪曲线图。
[0019] 图5为飞行器的速度误差跟踪曲线图。
[0020] 图6为飞行器的攻角跟踪曲线图。
[0021] 图7为飞行器的攻角误差跟踪曲线图。
[0022] 图8为飞行器的燃油当量比变化曲线图。
[0023] 图9为飞行器的升降舵偏转角变化曲线图。
[0024] 图10为虚拟俯仰角速度控制量和滤波俯仰角速度输出信号的跟踪曲线图。
[0025] 图11为滤波俯仰角速度输出信号与虚拟俯仰角速度控制量的差值跟踪曲线图。
[0026] 图12为速度转换误差的跟踪曲线图。
[0027] 图13为速度误差转换函数的跟踪曲线图。
[0028] 图14为攻角转换误差的跟踪曲线图。
[0029] 图15为攻角误差转换函数的跟踪曲线图。
[0030] 图16为齐次高阶滑模观测器对速度误差的观测跟踪曲线图。
[0031] 图17为齐次高阶滑模观测器对速度集总扰动的观测跟踪曲线图。
[0032] 图18为齐次高阶滑模观测器对攻角误差的观测跟踪曲线图。
[0033] 图19为齐次高阶滑模观测器对攻角集总扰动的观测跟踪曲线图。
[0034] 图20为齐次高阶滑模观测器对俯仰角速度的观测跟踪曲线图。
[0035] 图21为齐次高阶滑模观测器对俯仰角速度集总扰动的观测跟踪曲线图。
[0036] 图22为飞行器预设性能滑模控制器自适应参数k1的变化曲线图。
[0037] 图23为飞行器预设性能滑模控制器自适应参数k2的变化曲线图。
[0038] 图24为飞行器预设性能滑模控制器滑模变量s的变化曲线图。

具体实施方式

[0039] 下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
[0040] 本发明的目的是提供一种基于智能网络模型的飞行器预设性能滑模控制方法及装置,解决现有控制算法无法保证系统状态瞬态与稳态性能,可能导致高速飞行器吸气式发动机熄火的问题,提高建模精度,抑制了系统抖振,提升的收敛速度、控制精度与鲁棒性。
[0041] 为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
[0042] 如图1所示,本发明提供的基于智能网络模型的飞行器预设性能滑模控制方法包括步骤100至步骤600。
[0043] 步骤100:根据当前时刻飞行器的速度、攻角、升降舵偏转角及燃油当量比,基于智能网络模型确定当前误差动力学函数。
[0044] 所述当前误差动力学函数包括:当前第一速度误差动力学函数、当前第二速度误差动力学函数、当前攻角误差动力学函数、当前第一俯仰角速度误差动力学函数及当前第二俯仰角速度误差动力学函数。
[0045] 所述智能网络模型为预先采用飞行器的历史飞行量测数据集对神经网络进行训练得到的。所述历史飞行量测数据集中包括历史设定时段内各时刻飞行器的速度、攻角、升降舵偏转角、燃油当量比及对应的误差动力学函数。
[0046] 本实施例中,神经网络采用长短时记忆网络(Long Short‑Term Memory,LSTM)。LSTM使用三个控制(输入门、输出门与遗忘门)和细胞单元状态交互建立模型。其存储单元模型如公式(1)。
[0047] (1)。
[0048] 其中,it为输入门在t时刻的状态,ft为遗忘门在t时刻的状态,ct为细胞单元在t时刻的状态,ct‑1为细胞单元在t‑1时刻的状态,ot为输出门在t时刻的状态,ht为LSTM存储单元在t时刻的隐藏层状态,ht‑1为LSTM存储单元在t‑1时刻的隐藏层状态,xt为LSTM存储单元在t时刻的输入, 为sigmoid激活函数,Wxi为输入数据对输入门的权重矩阵,Whi为隐藏层对输入门的权重矩阵,Wci为细胞单元对输入门的权重矩阵,Wxf为输入数据对遗忘门的权重矩阵,Whf为隐藏层对遗忘门的权重矩阵,Wcf为细胞单元对遗忘门的权重矩阵,Wxc为输入数据对细胞单元的权重矩阵,Whc为隐藏层对细胞单元的权重矩阵,Wxo为输入数据对输出门的权重矩阵,Who为隐藏层对输出门的权重矩阵,Wco为细胞单元对输出门的权重矩阵,bi为输入门的偏置向量,bf为遗忘门的偏置向量,bc为细胞单元的偏置向量,bo为输出门的偏置向量。
[0049] t时刻的隐藏层状态ht与输出层采用全连接的方式计算得到t时刻智能网络模型的输出yt,如公式(2)。
[0050] (2)。
[0051] 其中,Wd为全连接层的权重矩阵,bd为全连接层的偏置向量。
[0052] 步骤200:根据所述当前误差动力学函数及飞行器传统误差动力学模型,建立飞行器智能误差动力学模型。
[0053] 传统飞行器动力学模型如公式(3)。
[0054] (3)。
[0055] 其中,V为速度, 为攻角,Q为俯仰角速度,T为推力,D为阻力,m为质量,g为重力加速度, 为飞行路径角,L为升力,M为俯仰力矩,Iyy为转动惯量, 为V对时间的一阶导数, 为 对时间的一阶导数, 为Q对时间的一阶导数。
[0056] 各个力和力矩经由曲线拟合近似可得到公式(4)。
[0057] (4)。
[0058] 其中, 为燃油当量比, 、 、 、 为推力系数导数,为推力拟合参数,h为高度, 为动压,S为参考面积, 和
为升力系数导数, 为推力矩耦合系数, 为平均气动弦长, 、 、
、 为力矩系数导数, 为升降舵偏转角,d1和d2为外部扰动,在动力学方程中分别体现于 和 中, 为大气密度, 、h0、hs为大气密度的拟合参数。
[0059] 定义速度指令Vd由阶跃信号Vc通过二阶光滑滤波器得到,如公式(5)。
[0060] (5)。
[0061] 其中, 为角频率, 为阻尼比,s1为拉普拉斯变换对应于时域的符号或算子。
[0062] 定义攻角指令 由攻角最大限幅 和正弦信号组成,如公式(6)。
[0063] (6)。
[0064] 其中,t为时刻, 和 为正弦信号频率, 和 为正弦信号振幅,且, , 为常数,且 。
[0065] 分别定义速度误差 ,攻角误差 。
[0066] 根据公式(3)所示的传统飞行器动力学模型可建立飞行器传统误差动力学模型,如公式(7)。
[0067] (7)。
[0068] 其中, 为 对时间的一阶导数,FV为第一速度误差动力学函数,GV为第二速度误差动力学函数,dV为速度外部扰动, 为Vd对时间的一阶导数, 为 对时间的一阶导数, 为攻角误差动力学函数, 为 对时间的一阶导数,FQ为第一俯仰角速度误差动力学函数,GQ为第二俯仰角速度误差动力学函数,dQ为俯仰角速度外部扰动。
[0069] 上述各误差动力学函数和外部扰动如公式(8)。
[0070] (8)。
[0071] 由于飞行器气动系数和结构系数在飞行过程中均会受到参数摄动影响,定义Cr为气动系数(实际气动拟合参数和推力拟合参数),C0为其标称值;Pr为结构系数(包含m、Iyy和S),P0为其标称值。而在传统飞行器控制器设计中,仅有标称参数可以获得,FV0、GV0、 、FQ0与GQ0即为使用标称值C0与P0得到的误差动力学函数。
[0072] 本发明建立的飞行器智能误差动力学模型包括速度智能误差动力学模型及攻角智能误差动力学模型。
[0073] 速度智能误差动力学模型如公式(9),攻角智能误差动力学模型如公式(10)。
[0074] (9)。
[0075] (10)。
[0076] 其中, 为当前第一速度误差动力学函数, 为当前第二速度误差动力学函数, 为速度集总扰动, 为当前攻角误差动力学函数, 为攻角集总扰动, 为当前第一俯仰角速度误差动力学函数, 为当前第二俯仰角速度误差动力学函数, 为俯仰角速度集总扰动。 、 和 由外部扰动和参数摄动引起。
[0077] 步骤300:确定速度设定时间预设性能函数及攻角设定时间预设性能函数。
[0078] 速度设定时间预设性能函数如公式(11)和公式(12)。
[0079] (11)。
[0080] (12)。
[0081] 其中, 为速度设定时间预设性能函数, 为速度设定时间预设性能函数的第一系数, 为速度设定时间预设性能函数的第二系数, 为速度设定时间预设性能函数的第三系数,t为时刻, 为速度误差的初值, 为速度设定时间预设性能函数的变化方向, 为速度设定时间预设性能函数的稳态收敛值, 为速度误差收敛到 的设定时间。 、 、 与 为待设计参数。
[0082] 当 时, 及其二阶导数满足公式(13)的定义。
[0083] (13)。
[0084] 攻角设定时间预设性能函数如公式(14)。
[0085] (14)。
[0086] 其中, 为攻角设定时间预设性能函数, 为攻角设定时间预设性能函数的第一系数, 为攻角设定时间预设性能函数的第二系数, 为攻角设定时间预设性能函数的第三系数, 为攻角误差的初值, 为攻角设定时间预设性能函数的变化方向, 为攻角设定时间预设性能函数的稳态收敛值, 为攻角误差收敛到
的设定时间。 、 、 与 为待设计参数。
[0087] 当 时, 及其二阶导数满足公式(15)的定义。
[0088] (15)。
[0089] 步骤400:基于所述速度设定时间预设性能函数、所述攻角设定时间预设性能函数及所述飞行器智能误差动力学模型,确定速度控制器及攻角转换误差动力学模型。
[0090] 具体地,将飞行器智能误差动力学模型改写为转换误差动力学模型。步骤400包括如下步骤。
[0091] 1)为保证高速飞行器飞行过程系统状态瞬态与稳态性能要求与吸气式发动机进气约束条件,需对 和 进行限幅。即基于所述速度设定时间预设性能函数,对速度误差进行限幅,以确定速度预设性能条件。基于所述攻角设定时间预设性能函数,对攻角误差进行限幅,以确定攻角预设性能条件。
[0092] 速度预设性能条件如公式(16),攻角预设性能条件如公式(17)。
[0093] (16)。
[0094] (17)。
[0095] 其中, 为速度设定时间预设性能函数的下界系数, 为速度设定时间预设性能函数的上界系数, 为t时刻的速度误差, 为攻角设定时间预设性能函数的下界系数, 为攻角设定时间预设性能函数的上界系数, 为t时刻的攻角误差。 、、 、 为 之间的正数,一般均取值为1。
[0096] 2)确定所述速度智能误差动力学模型的速度预设性能控制器,并基于所述速度预设性能条件,采用齐次高阶滑模观测器对速度集总扰动进行补偿,以确定速度控制器。
[0097] 速度预设性能控制器如公式(18)。
[0098] (18)。
[0099] 其中, 为速度转换误差, 为速度误差转换函数,其表达式如公式(19)。
[0100] (19)。
[0101] 从而,可以得到 及其导数 的表达式如公式(20)。
[0102] (20)。
[0103] 其中, 为t时刻的速度转换误差, 为 对时间的一阶导数,为t时刻的速度中间参数, 为 对时间的一阶导数,其表达式如公式(21)。
[0104] (21)。
[0105] 加入齐次高阶滑模观测器对速度集总扰动 补偿,如公式(22)。
[0106] (22)。
[0107] 其中, 为 的估计值, 为 的估计值, 为 的估计值, 为的估计值, 为 的一阶导数, 为 的二阶导数, 为 的一阶导数,
为 的一阶导数, 为 的一阶导数, 为 的一阶导数, 、 、
、 为中间变量,LV为正的待设计参数, , , ,

[0108] 经过有限时间后, , 。
[0109] 因此,速度控制器设计为公式(23)的形式。
[0110] (23)。
[0111] 其中, 为速度正待设计参数, 为速度中间参数。
[0112] 进一步地,速度转换误差动力学模型如公式(24)。
[0113] (24)。
[0114] 其中, 为速度集总扰动误差, 。若速度转换误差 有界,则速度误差转换函数 ,因此速度子系统满足速度预设性能条件。
[0115] 3)确定所述攻角智能误差动力学模型的攻角预设性能控制器,并基于所述攻角预设性能条件,采用齐次高阶滑模观测器对速度集总扰动及俯仰角速度集总扰动进行补偿,以确定虚拟俯仰角速度控制量。
[0116] 针对公式(10)中的攻角智能误差动力学模型,首先使用攻角预设性能控制器对进行转换,如公式(25)。
[0117] (25)。
[0118] 其中, 为攻角转换误差, 为攻角误差转换函数,其表达式如公式(26)。
[0119] (26)。
[0120] 从而,可以得到 及其导数 的表达式如公式(27)。
[0121] (27)。
[0122] 为t时刻的攻角转换误差, 为 对时间的一阶导数, 为t时刻的攻角中间参数, 为 对时间的一阶导数,其表达式如公式(28)。
[0123] (28)。
[0124] 加入齐次高阶滑模观测器对速度集总扰动 补偿,如公式(29)。
[0125] (29)。
[0126] 其中, 为 的估计值, 为 的估计值, 为 的估计值, 为的估计值, 为 的一阶导数, 为 的二阶导数, 为 的一阶导数,
为 的一阶导数, 为 的一阶导数, 为 的一阶导数, 、 、
、 为中间变量,LV为正的待设计参数。
[0127] 经过有限时间后, , 。
[0128] 因此,虚拟俯仰角速度控制量Qc的表达式如公式(30)。
[0129] (30)。
[0130] 其中, 为攻角正待设计参数, 为攻角中间参数。
[0131] 4)根据所述虚拟俯仰角速度控制量,采用一阶滤波器确定滤波俯仰角速度输出信号Qd,如公式(31)。
[0132] (31)。
[0133] 其中, 为时间常量, 为Qd的一阶导数。
[0134] 5)根据所述滤波俯仰角速度输出信号,确定攻角转换误差动力学模型。
[0135] 具体地,定义俯仰角速度误差 与滤波俯仰角速度输出信号与虚拟俯仰角速度控制量的差值 ,则攻角转换误差动力学模型如公式(32)。
[0136] (32)。
[0137] 其中, 为 对时间的一阶导数, 为攻角集总扰动误差, 。
[0138] 若攻角转换误差 有界,则攻角误差转换函数 ,因此攻角子系统满足攻角预设性能条件。
[0139] 同样对俯仰角速度的动力学方程加入齐次高阶滑模观测器,对俯仰角速度集总扰动补偿,如公式(33)。
[0140] (33)。
[0141] 其中, 为 的估计值, 为 的估计值, 为 的估计值, 为的估计值, 为 的一阶导数, 为 的二阶导数, 为 的一阶导数,
为 的一阶导数, 为 的一阶导数, 为 的一阶导数, 、 、
、 为中间变量,LQ为正的待设计参数。
[0142] 经过有限时间后, , 。
[0143] 步骤500:基于所述攻角转换误差动力学模型,采用自适应超螺旋滑模算法,建立飞行器预设性能超螺旋滑模控制器,以确定攻角控制器。
[0144] 本实施例中,所述自适应超螺旋滑模算法包括非奇异快速终端滑模面及改进的自适应超螺旋趋近律。
[0145] 非奇异快速终端滑模面如公式(34)。
[0146] (34)。
[0147] 改进的自适应超螺旋趋近律如公式(35)。
[0148] (35)。
[0149] 其中,s为滑模变量, 为s的导数, 、 、k1、k2、k3均为正待设计参数,k1、k2满足公式(36)的自适应律,a1和a2为常数,a1>a2,1
[0150] (36)。
[0151] 其中, 为k1的一阶导数, 、 、 、 、 、 均为正参数,且 。本发明的自适应律可以实现控制增益灵活调整。具体表现形式为:当 时,
,k1增加,实现 减小;当 时, ,k1减小;当 ,
,k1不变。若系统稳定条件不满足,会导致 再次增大,自适应律再次作用,使系统趋于稳定。本自适应律避免了控制增益过小导致系统发散与控制增益过大导致系统剧烈抖振的问题。
[0152] 结合自适应超螺旋滑模算法与攻角转换误差动力学模型,建立飞行器预设性能超螺旋滑模控制器,得到攻角实际控制律,如公式(37)。
[0153] (37)。
[0154] 其中, 为攻角实际控制律,即升降舵偏转角, 为攻角控制器中间参数,kQ为俯仰角速度正待设计参数, 为 的一阶导数。 的定义如公式(38)。
[0155] (38)。
[0156] 的动力学方程如公式(39)。
[0157] (39)。
[0158] 为 的一阶导数, 为俯仰角速度集总扰动误差, 。
[0159] 公式(34)的非奇异快速终端滑模面由攻角转换误差 和俯仰角速度误差 组成,公式(32)给出攻角转换误差的导数。通过对公式(34)求导,代入公式(32)和公式(10)中的俯仰角速度误差动力学方程,从而得到一个含有升降舵偏转角 的表达式,令该表达式与公式(35)改进的自适应超螺旋趋近律相等,即可推导出公式(37)的攻角实际控制律。
[0160] 步骤600:基于所述速度控制器及所述攻角控制器对飞行器进行控制。
[0161] 本发明使用智能网络模型进行精确动力学建模,提升了飞行器智能误差动力学模型的精度。采用预设性能控制保证高速飞行器飞行过程中发动机进气条件约束,通过满足攻角幅值条件避免飞行器发动机出现壅塞而熄火。相比传统性能函数,本发明通过设定时间性能函数可以保证系统飞行状态跟踪误差在设定时间收敛,提升了收敛速度与收敛精度,避免了初始时刻控制量过大的问题。本发明设计了自适应超螺旋滑模算法,其中,非奇异快速终端滑模面避免了传统终端滑模的奇异问题并加快系统收敛速度;自适应超螺旋趋近律自动调整控制增益,避免了控制增益过小导致系统发散与控制增益过大导致系统剧烈抖振的问题。
[0162] 以下通过MATLAB仿真说明本发明提供的基于智能网络模型的飞行器预设性能滑模控制方法的具体实施过程。
[0163] 首先假设飞行器气动系数与结构系数受飞行环境参数摄动影响后,气动系数Cr与结构系数Pr的表达式如公式(40)。
[0164] (40)。
[0165] 然后假设外部扰动d1和d2分别发生在第70秒和第40秒,表达形式如公式(41)。
[0166] (41)。
[0167] 飞行器参数设置为: ,m=4378.17kg,Vc=2407.92m/s, , , , , ,
, 。
[0168] 初始飞行状态设置为:h=26212.8m,V=2347.59m/s, ,Q=0°/s,, 。
[0169] LSTM的超参数设置为:学习率为0.01,训练次数为180,隐藏层数量为5,隐藏层单元数为64,批样本数为800,输入维度为4。
[0170] 预设性能滑模控制器的参数的设定为: , ,, , , , , ,
, , , , , ,
, , , , , , ,
, , , , , , ,
, 。
[0171] 本发明基于上述设置的参数进行智能网络模型的仿真训练,从仿真训练结果可确定本发明构建的智能网络模型可以对未知数据进行预测,从而提高在外部扰动情形下高速飞行器的建模精度。图2和图3展示了智能网络模型对训练数据的拟合误差和均方误差。智能网络模型对预测集的预测结果显示,决定系数R2=0.9995、均方误差RMSE=0.0906,这表明了飞行器智能网络模型对时序飞行数据的拟合具有良好的外推性和泛化性能。
[0172] 为预设性能滑模控制器对速度指令的跟踪性能如图4、图5、图12和图13所示。预设性能滑模控制器对攻角指令的跟踪性能如图6、图7、图14和图15。在飞行器预定性能滑模控制器的作用下,速度跟踪精度达到 ,攻角跟踪精度达到 ,均可满足期望性能边界要求并快速收敛,而且攻角变化范围始终不会超过最大幅值约束,满足吸气式发动机正常工作条件。此外,转换误差函数始终可以限制在预定性能边界
与 内,进一步保证飞行器系统状态的瞬态与稳态性能要
求。
[0173] 图8和图9的曲线描述了控制器执行机构均在合理范围内变化。图10和图11给出通过一阶低通滤波器后Qd对虚拟俯仰角速度控制量Qc的跟踪情况。齐次高阶滑模观测器对、 、Q、 、 、 的观测效果如图16至图21所示,使用HSMO代指齐次高阶滑模观测器。图22和图23给出自适应参数变化曲线,展示其自动调节增益大小的攻能。图24给出滑模变量的变化曲线,进一步表明本发明的快速收敛性与强鲁棒性。
[0174] 在一个实施例中,提供了一种计算机装置,包括存储器和处理器以存储在存储器上并可在处理器上运行的计算机程序,该处理器执行计算机程序时实现上述基于智能网络模型的飞行器预设性能滑模控制方法的步骤。
[0175] 在一个实施例中,提供了一种计算机可读存储介质,其上存储有计算机程序,该计算机程序被处理器执行时实现上述基于智能网络模型的飞行器预设性能滑模控制方法的步骤。
[0176] 在一个实施例中,提供了一种计算机程序产品,包括计算机程序,该计算机程序被处理器执行时实现上述基于智能网络模型的飞行器预设性能滑模控制方法的步骤。
[0177] 以上实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。
[0178] 本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。
QQ群二维码
意见反馈