CN111103795A - 基于智能自适应优化算法的高超声速飞行器再入段轨迹优化控制器 - Google Patents

基于智能自适应优化算法的高超声速飞行器再入段轨迹优化控制器 Download PDF

Info

Publication number
CN111103795A
CN111103795A CN201911149108.6A CN201911149108A CN111103795A CN 111103795 A CN111103795 A CN 111103795A CN 201911149108 A CN201911149108 A CN 201911149108A CN 111103795 A CN111103795 A CN 111103795A
Authority
CN
China
Prior art keywords
aircraft
class
hypersonic
flight
visibility
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
CN201911149108.6A
Other languages
English (en)
Other versions
CN111103795B (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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN201911149108.6A priority Critical patent/CN111103795B/zh
Publication of CN111103795A publication Critical patent/CN111103795A/zh
Application granted granted Critical
Publication of CN111103795B publication Critical patent/CN111103795B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive 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/042Adaptive 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

Abstract

本发明公开了一种基于智能自适应优化算法的高超声速飞行器再入段轨迹优化控制器,该控制器由飞行器海拔高度传感器、飞行器速度传感器、飞行器飞行航道倾角传感器、飞行器水平飞行距离传感器、飞行器微控制单元MCU、飞行器攻角控制器构成。高超声速飞行器到达再入段空域后,各传感器开启,获得飞行器海拔高度、速度、飞行航道倾角和飞行水平距离信息,MCU根据所获得的信息自动执行内部智能自适应优化算法,得到水平飞行距离最长的轨迹优化控制策略,并将其转换为控制指令发送给攻角控制器执行。本发明能够快速地得到轨迹优化控制策略,使高超声速飞行器获得更长的水平飞行距离。

Description

基于智能自适应优化算法的高超声速飞行器再入段轨迹优化 控制器
技术领域
本发明涉及高超声速飞行器再入段轨迹优化领域,主要是一种基于智能自适应优化算法的高超声速飞行器再入段轨迹优化控制器。在高超声速飞行器到达再入段后能够给出高超声速飞行器轨迹优化控制策略并转换为飞行器攻角控制指令,使高超声速飞行器获得更长的水平飞行距离。
背景技术
高超声速飞行器是实现远程快速精确打击和全球快速到达的新型飞行器,在未来的军事、政治和经济中具有十分重要的战略地位,已经成为世界航空航天领域一个极其重要的发展方向,研究和发展高超声速飞行器在开发太空和国家安全方面具有非常重要的意义。
在高超声速飞行器的研究中,轨迹优化不仅有利于提高飞行器飞行品质以满足既定任务要求,同时也是完成飞行任务的重要保证和实现机动飞行的必要条件,近些年来一直受到国内外各军事强国的重视,是当前国内外研究的热点和难点。
由于从大气从外缘进入大气层,高度和速度的变化范围很大,高超声速飞行器面临各种严峻的再入环境,再入段轨迹优化技术则是保证高超声速飞行器完成飞行任务的关键,对提高其打击范围和落点精度具有更重要的实用价值。因此,研究高效的高超声速飞行器再入段轨迹优化方法显得尤为重要。
发明内容
为了使高超声速飞行器获得更长的水平飞行距离,提高高超声速飞行器的打击范围,本发明提供了一种基于智能自适应优化算法的高超声速飞行器再入段轨迹优化控制器。
高超声速飞行器再入段轨迹优化问题可以描述为
Figure BDA0002283046190000011
其中t表示时间,x(t)表示高超声速飞行器的状态变量,x1(t)表示飞行器海拔高度、x2(t)表示飞行器速度、x3(t)表示飞行器飞行航道倾角、x4(t)表示飞行器水平飞行距离,u(t)表示高超声速飞行器的攻角控制量,为本问题的控制变量;
Figure BDA0002283046190000021
表示状态变量x(t)的一阶导数,F(x(t),u(t),t)是根据高超声速飞行器再入段三维空间运动方程建立的微分方程组数学模型;t0表示再入段轨迹优化开始的时间点,h0表示优化开始时刻飞行器的初始海拔高度,v0表示优化开始时刻飞行器的初始速度,γ0表示优化开始时刻飞行器的初始飞行航道角,r0表示优化开始时刻飞行器的初始水平飞行距离,tf表示再入段轨迹优化结束时间点,hf表示优化结束时刻飞行器的海拔高度,vf表示优化结束时刻飞行器的速度,γf表示优化结束时刻飞行器的飞行航道角;J[u(t)]表示高超声速飞行器轨迹优化的目标函数即优化结束时刻飞行器的水平飞行距离,G[u(t),x(t),t]是高超声速飞行器再入段过程的约束条件,umin和umax表示攻角控制范围的下限值和上限值。
本发明解决其技术问题所采用的技术方案是:一种基于智能自适应优化算法的高超声速飞行器再入段轨迹优化控制器,它主要由飞行器海拔高度传感器、飞行器速度传感器、飞行器飞行航道倾角传感器、飞行器水平飞行距离传感器、高超声速飞行器MCU、飞行器攻角控制器构成。所述飞行器海拔高度传感器、飞行器速度传感器、飞行器飞行航道倾角传感器、飞行器水平航程传感器和飞行器攻角控制器均通过数据总线与高超声速飞行器MCU相连。所述的高超声速飞行器MCU包括信息采集模块、初始化模块、系统状态求解模块、收敛性判断模块、智能自适应优化算法模块、控制指令输出模块。
所述基于智能自适应优化算法的高超声速飞行器再入段轨迹优化控制器的运行过程如下:
步骤1):在高超声速飞行器MCU中输入对应于该飞行器的气动系数模型、飞行器性能约束条件、指定优化目标;
步骤2):高超声速飞行器到达再入段后,开启飞行器海拔高度传感器、飞行器速度传感器、飞行器飞行航道倾角传感器和飞行器水平飞行距离传感器,得到高超声速飞行器当前的海拔高度、速度、飞行航道倾角和飞行水平距离状态信息;
步骤3):飞行器MCU根据设定的海拔高度、速度、飞行航道倾角要求执行内部智能自适应优化算法,得到使高超声速飞行器水平飞行距离最长的轨迹优化控制策略;
步骤4):高超声速飞行器MCU将获得的轨迹优化控制策略发送给控制策略输出模块,并转换为控制指令发送给飞行器攻角控制器执行。
所述步骤3包括以下子步骤:
步骤3.1):信息采集模块获取步骤2得到的高超声速飞行器当前的海拔高度、速度、飞行航道倾角和飞行水平距离状态信息;
步骤3.2):采用控制变量参数化方法对系统模型的控制变量进行离散。初始化各项参数,设置轨迹优化过程时间的分段数NE,在全局范围内随机生成数量为N的初始细胞种群,设置细胞伸缩间隔的常量r1和r2,细胞可见度的常量v1和v2,最大迭代次数max_iters,初始迭代次数k=0。
步骤3.3):通过系统状态求解模块获得第k次迭代种群个体i的目标函数值
Figure BDA0002283046190000031
记录种群的平均目标值
Figure BDA0002283046190000032
及最优目标值
Figure BDA0002283046190000033
更新第k次迭代的种群最优解
Figure BDA0002283046190000034
步骤3.4):执行收敛性判断模块,若当前迭代次数k达到最大迭代次数max_iters,则停止迭代,以当前种群最优解作为最优控制策略的离散解,经控制指令输出模块输出;否则,迭代次数k加1,执行步骤3.5。
步骤3.5):执行智能自适应优化算法模块,将种群按照模糊C均值算法分为C1、C2、C3和C4四类,各类中的个体通过生成反射reflection及可见性visibility来更新细胞状态并独立地进行搜索,其中C1、C4类在全局范围内搜索最优解,C2、C3类在最优解的邻域内寻优,同时各类根据个体的适应度自适应调整权重分布,执行完此步后跳转至步骤3.3)。
所述步骤3.3中,所述的系统状态求解模块,采用如下求解公式:
Figure BDA0002283046190000035
其中,tm表示选择的积分时刻,tm+1表示位于时刻tm后的积分时刻,并且tm+1=tm+h,h为积分步长,F(·)是描述状态微分方程的函数,K1、K2、K3、K4分别表示积分过程中的4个节点的函数值。
所述步骤3.5具体包括以下子步骤:
步骤3.5.1):利用模糊C均值算法将种群分为C1、C2、C3、C4四类,其中C1类及C4类在全局范围内进行搜索,C2类及C3类在最优解的邻域内搜索,各类通过生成反射及可见性得到更新的细胞Snew,由公式(2)确定
Snew=reflection+visibility (2)
各类更新细胞状态的过程如下:
步骤3.5.2):C1类中的细胞由式(3)和式(4)生成反射及可见性:
Figure BDA0002283046190000041
Figure BDA0002283046190000042
其中
Figure BDA0002283046190000043
表示C1中第i个细胞在第j维上的位置,Soptj表示最优解在第j维上的位置,j={1,2,3},分别表示x、y、z坐标,R和V表示细胞的反光程度与可见程度,可由式(5)和(6)得到:
R=random()×(r1-r2)+r2 (5)
V=random()×(v1-v2)+v2 (6)
式中r1和r2是表示细胞伸缩间隔的两个常数,v1和v2是表示细胞可见度的两个常数,初始化参数时已给出。
步骤3.5.3):C2类中细胞的反射与可见性由式(7)和式(8)确定:
reflectionj=R×Soptj (7)
Figure BDA0002283046190000044
其中
Figure BDA0002283046190000045
表示C2中第i个细胞在第j维上的位置,R和V仍由式(5)(6)得到。
步骤3.5.4):C3类中细胞反射过程由式(7)确定,用(9)式表示其可见性:
Figure BDA0002283046190000046
其中
Figure BDA0002283046190000047
表示Soptj平均值。
步骤3.5.5):C4类由式(10)随机更新个体
Snewj=random()×(Lj-lj)+lj (10)
其中Lj和lj表示在第j维上搜索的上、下界。
对通过上述细胞更新过程得到的Snew,若与当前最优解Sopt对应的目标函数值满足关系
Figure BDA0002283046190000048
则用Snew替换当前最优解。
步骤3.5.6):为提高各类细胞的搜索效率,增强C1、C4类个体的全局搜索能力,改善C2、C3类的局部寻优能力,根据个体适应度自适应调整权重分布:具体如下:
3.5.6.1)对于C1中适应度较差的个体,即目标值小于平均值Javg的个体,应适当的增大其权重,对于大于平均值的个体,采用余弦函数设置权重,以保持其较强的全局搜索能力,并且在后期能进行更细致的搜索。
Figure BDA0002283046190000051
其中wa、wb、和wc为限制权重取值范围的常数,且满足wa>wb>wc,k为当前迭代次数,max_iters为最大迭代次数。
3.5.6.2)C2类对目标值越小的个体权重设置越大,同时为改善种群迭代过程中的多样性,对适应度较好的即接近全局最优解的个体,设置较大的权重。
Figure BDA0002283046190000052
其中Jopt为当前最优解对应的目标值。
3.5.6.3)C3类通过余弦函数更新权重,为提高该类个体的局部搜索能力,随着迭代次数的增大不断缩小搜索范围。
Figure BDA0002283046190000053
其中wd为小于1的常数。
在新一轮的迭代过程中各类细胞生成反射及可见性时需要加入上述权重因子,更新细胞状态的公式变形为
Figure BDA0002283046190000054
其中wi∈{w1,w2,w3}。
本发明的有益效果主要表现在:所述的智能自适应优化算法通过对种群分类进行独立地搜索,并根据个体适应度自适应调整权重,保证种群在具有较强的全局搜索能力的同时,也能增强局部细致搜索能力,从而获得高寻优精度的轨迹优化攻角控制指令,提高高超声速飞行器打击范围。
附图说明
图1是一种基于智能自适应优化算法的高超声速飞行器再入轨迹优化控制器的结构示意图;
图2是一种基于智能自适应优化算法的高超声速飞行器再入轨迹优化控制器的高超声速飞行器MCU内部模块结构图。
具体实施方式
实施例1
基于智能自适应优化算法的高超声速飞行器再入轨迹优化控制器的结构如图1所示,假设高超声速飞行器到达再入段空域,高超声速飞行器海拔高度传感器3、飞行器速度传感器4、飞行航道倾角传感器5、水平飞行距离传感器6和MCU2均已开启。信息采集模块9立即采集飞行器进入再入段时的初始海拔高度、速度、飞行航道倾角和水平飞行距离,设当前初始时刻t0=0s,海拔高度传感器传入MCU的海拔高度为h0=80 000m,速度传感器传入MCU的速度为v0=6400m/s,飞行航道倾角传感器传入MCU的飞行航道倾角为γ0=-0.052rad,水平飞行距离传感器感器传入MCU的水平飞行距离为r0=0m;终值时刻tf高超声速飞行器需要满足的条件为海拔高度设定为hf=24000m,速度设定为vf=760m/s,飞行航道倾角设定为γf=-0.08rad;结合飞行器的三维空间运动方程、气动系数模型、飞行器性能约束条件和指定优化目标,得到该问题的数学模型如下:
Figure BDA0002283046190000071
其中L表示升力,D表示阻力,CL表示升力系数,CD表示阻力系数。
为了便于表述,采用F(x(t),u(t),t)表示高超声速飞行器再入段三维空间运动方程建立的微分方程组数学模型,即:
Figure BDA0002283046190000072
采用G[u(t),x(t),t]表示高超声速飞行器再入段过程的约束条件,为:
Figure BDA0002283046190000081
此外,J[u(t)]表示高超声速飞行器轨迹优化的目标函数即优化结束时刻飞行器的水平飞行距离。
高超声速飞行器MCU自动产生攻角控制指令的智能自适应优化算法如图2所示,其运行步骤如下:
步骤1):在高超声速飞行器MCU2中输入对应于该飞行器的气动系数模型、飞行器性能约束条件、指定优化目标1;
步骤2):高超声速飞行器到达再入段后,开启飞行器海拔高度传感器3、飞行器速度传感器4、飞行器飞行航道倾角传感器5和飞行器水平飞行距离传感器6,得到高超声速飞行器当前的海拔高度、速度、飞行航道倾角和飞行水平距离状态信息;
步骤3):飞行器MCU 2根据设定的海拔高度、速度、飞行航道倾角要求自动执行内部智能自适应优化算法模块13,得到使高超声速飞行器水平飞行距离最长的轨迹优化控制策略;
步骤4):高超声速飞行器MCU 2将获得的轨迹优化控制策略发送给控制策略输出模块,并转换为控制指令发送给飞行器攻角控制器8执行。
所述的高超声速飞行器MCU部分如图2所示,包括信息采集模块9、初始化模块10、系统状态求解模块11、收敛性判断模块12、智能自适应优化算法模块13、控制指令输出模块14。其中,信息采集模块9包括飞行器海拔高度和速度采集、飞行器飞行航道倾角和飞行水平距离采集、飞行器海拔高度和速度设定采集、飞行器飞行航道倾角设定采集、飞行器的气动系数模型和性能约束条件以及指定优化目标参数采集五个子模块。
所述的高超声速飞行器MCU自动产生攻角控制指令的智能自适应优化算法运行步骤如下:
步骤1):高超声速飞行器到达再入段后,飞行器海拔高度传感器、飞行器速度传感器、飞行器飞行航道倾角传感器和飞行器水平飞行距离传感器开启,信息采集模块9获取初始时刻t0=0s时高超声速飞行器海拔高度h0=80 000m,速度为v0=6400m/s,飞行航道倾角为γ0=-0.052rad,水平飞行距离传感器感器水平飞行距离设置为r0=0m;终值时刻tf高超声速飞行器海拔高度要求设定为hf=24000m,速度要求设定为vf=760m/s,飞行航道倾角要求设定为γf=-0.08rad;;
步骤2):初始化模块10开始运行,采用控制变量参数化方法对系统模型的控制变量进行离散。设置轨迹优化过程时间的分段数NE,在全局范围内随机生成数量为N的初始细胞种群,设置细胞伸缩间隔的常量r1和r2,细胞可见度的常量v1和v2,最大迭代次数max_iters,初始迭代次数k=0。
步骤3):通过系统状态求解模块11,获得第k次迭代种群个体i的目标函数值
Figure BDA0002283046190000091
记录种群的平均目标值
Figure BDA0002283046190000092
及最优目标值
Figure BDA0002283046190000093
更新第k次迭代的种群最优解
Figure BDA0002283046190000094
步骤4):执行收敛性判断模块12,若当前迭代次数k达到最大迭代次数max_iters,则停止迭代,以当前种群最优解作为最优控制策略的离散解,经控制指令输出模块输出;否则,迭代次数k加1。
步骤5):执行智能自适应优化算法模块13,将种群按照模糊C均值算法分为C1、C2、C3和C4四类,各类中的个体通过生成反射reflection及可见性visibility来更新细胞状态并独立地进行搜索,其中C1、C4类在全局范围内搜索最优解,C2、C3类在最优解的邻域内寻优,同时各类根据个体的适应度自适应调整权重分布,执行完此步后跳转至步骤3)。
所述的智能自适应优化算法模块13的运行过程如下:
步骤1):利用模糊C均值算法将种群分为C1、C2、C3、C4四类,其中C1类及C4类在全局范围内进行搜索,C2类及C3类在最优解的邻域内搜索,各类通过生成反射及可见性得到更新的细胞Snew,由公式(33)确定
Snew=reflection+visibility (33)
各类更新细胞状态的过程如下:
步骤2):C1类中的细胞由式(34)和式(35)生成反射及可见性:
Figure BDA0002283046190000095
Figure BDA0002283046190000096
其中
Figure BDA0002283046190000097
表示C1中第i个细胞在第j维上的位置,Soptj表示最优解在第j维上的位置,j={1,2,3},分别表示x、y、z坐标,R和V表示细胞的反光程度与可见程度,可由式(36)和(37)得到:
R=random()×(r1-r2)+r2 (36)
V=random()×(v1-v2)+v2 (37)
式中r1和r2是表示细胞伸缩间隔的两个常数,v1和v2是表示细胞可见度的两个常数,初始化参数时已给出。
步骤3):C2类中细胞的反射与可见性由式(38)和式(39)确定:
reflectionj=R×Soptj (38)
Figure BDA0002283046190000101
其中
Figure BDA0002283046190000102
表示C2中第i个细胞在第j维上的位置,R和V仍由式(36)、37)得到。
步骤4):C3类中细胞反射过程由式(38)确定,用(40)式表示其可见性:
Figure BDA0002283046190000103
其中
Figure BDA0002283046190000104
表示Soptj平均值。
步骤5):C4类由式(41)随机更新个体
Snewj=random()×(Lj-lj)+lj (41)
其中Lj和lj表示在第j维上搜索的上、下界。
对通过上述细胞更新过程得到的Snew,若与当前最优解Sopt对应的目标函数值满足关系
Figure BDA0002283046190000105
则用Snew替换当前最优解。
步骤6):为提高各类细胞的搜索效率,增强C1、C4类个体的全局搜索能力,改善C2、C3类的局部寻优能力,根据个体适应度自适应调整权重分布:
6.1)对于C1中适应度较差的个体,即目标值小于平均值Javg的个体,应适当的增大其权重,对于大于平均值的个体,采用余弦函数设置权重,以保持其较强的全局搜索能力,并且在后期能进行更细致的搜索。
Figure BDA0002283046190000106
其中wa、wb、和wc为限制权重取值范围的常数,且满足wa>wb>wc,k为当前迭代次数,max_iters为最大迭代次数。
6.2)C2类对目标值越小的个体权重设置越大,同时为改善种群迭代过程中的多样性,对适应度较好的即接近全局最优解的个体,设置较大的权重。
Figure BDA0002283046190000111
其中Jopt为当前最优解对应的目标值。
6.3)C3类通过余弦函数更新权重,为提高该类细胞的局部搜索能力,随着迭代次数的增大不断缩小搜索范围。
Figure BDA0002283046190000112
其中wd为小于1的常数。
在新一轮的迭代过程中各类细胞生成反射及可见性时需要加入上述权重因子,更新细胞状态的公式变形为
Figure BDA0002283046190000113
其中wi∈{w1,w2,w3}。
最后,飞行器MCU将获得的优化轨迹作为指令输出到控制策略输出模块,转换为控制指令发送给攻角控制器8,完成轨迹优化的执行。
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明,不能认定本发明的具体实施只限于这些说明。对于本发明所属技术领域的普通技术人员来说,在不脱离发明构思的前提下,还可以做出若干简单推演或替换,都应当视为属于本发明的保护范围。

Claims (3)

1.一种基于智能自适应优化算法的高超声速飞行器再入段轨迹优化控制器,其特征在于,它主要由飞行器海拔高度传感器、飞行器速度传感器、飞行器飞行航道倾角传感器、飞行器水平飞行距离传感器、高超声速飞行器MCU、飞行器攻角控制器构成。所述飞行器海拔高度传感器、飞行器速度传感器、飞行器飞行航道倾角传感器、飞行器水平航程传感器和飞行器攻角控制器均通过数据总线与高超声速飞行器MCU相连。所述的高超声速飞行器MCU包括信息采集模块、初始化模块、系统状态求解模块、收敛性判断模块、智能自适应优化算法模块、控制指令输出模块。
所述基于智能自适应优化算法的高超声速飞行器再入段轨迹优化控制器的运行过程如下:
步骤1):在高超声速飞行器MCU中输入对应于该飞行器的气动系数模型、飞行器性能约束条件、指定优化目标;
步骤2):高超声速飞行器到达再入段后,开启飞行器海拔高度传感器、飞行器速度传感器、飞行器飞行航道倾角传感器和飞行器水平飞行距离传感器,得到高超声速飞行器当前的海拔高度、速度、飞行航道倾角和飞行水平距离状态信息;
步骤3):飞行器MCU根据设定的海拔高度、速度、飞行航道倾角要求执行内部智能自适应优化算法,得到使高超声速飞行器水平飞行距离最长的轨迹优化控制策略;
步骤4):高超声速飞行器MCU将获得的轨迹优化控制策略发送给控制策略输出模块,并转换为控制指令发送给飞行器攻角控制器执行。
所述步骤3包括以下子步骤:
步骤3.1):信息采集模块获取步骤2得到的高超声速飞行器当前的海拔高度、速度、飞行航道倾角和飞行水平距离状态信息;
步骤3.2):采用控制变量参数化方法对系统模型的控制变量进行离散。初始化各项参数,设置轨迹优化过程时间的分段数NE,在全局范围内随机生成数量为N的初始细胞种群,设置细胞伸缩间隔的常量r1和r2,细胞可见度的常量v1和v2,最大迭代次数max_iters,初始迭代次数k=0。
步骤3.3):通过系统状态求解模块获得第k次迭代种群个体i的目标函数值
Figure FDA0002283046180000011
记录种群的平均目标值
Figure FDA0002283046180000012
及最优目标值
Figure FDA0002283046180000013
更新第k次迭代的种群最优解
Figure FDA0002283046180000014
步骤3.4):执行收敛性判断模块,若当前迭代次数k达到最大迭代次数max_iters,则停止迭代,以当前种群最优解作为最优控制策略的离散解,经控制指令输出模块输出;否则,迭代次数k加1,执行步骤3.5。
步骤3.5):执行智能自适应优化算法模块,将种群按照模糊C均值算法分为C1、C2、C3和C4四类,各类中的个体通过生成反射reflection及可见性visibility来更新细胞状态并独立地进行搜索,其中C1、C4类在全局范围内搜索最优解,C2、C3类在最优解的邻域内寻优,同时各类根据个体的适应度自适应调整权重分布,执行完此步后跳转至步骤3.3)。
2.根据权利要求1所示基于智能自适应优化算法的高超声速飞行器再入段轨迹优化控制器,其特征在于,所述步骤3.3中,所述的系统状态求解模块,采用如下求解公式:
Figure FDA0002283046180000021
其中,tm表示选择的积分时刻,tm+1表示位于时刻tm后的积分时刻,并且tm+1=tm+h,h为积分步长,F(·)是描述状态微分方程的函数,K1、K2、K3、K4分别表示积分过程中的4个节点的函数值。
3.根据权利要求1所示基于智能自适应优化算法的高超声速飞行器再入段轨迹优化控制器,其特征在于,所述步骤3.5具体包括以下子步骤:
步骤3.5.1):利用模糊C均值算法将种群分为C1、C2、C3、C4四类,其中C1类及C4类在全局范围内进行搜索,C2类及C3类在最优解的邻域内搜索,各类通过生成反射及可见性得到更新的细胞Snew,由公式(2)确定
Snew=reflection+visibility (2)
各类更新细胞状态的过程如下:
步骤3.5.2):C1类中的细胞由式(3)和式(4)生成反射及可见性:
Figure FDA0002283046180000022
Figure FDA0002283046180000023
其中
Figure FDA0002283046180000024
表示C1中第i个细胞在第j维上的位置,Soptj表示最优解在第j维上的位置,j={1,2,3},分别表示x、y、z坐标,R和V表示细胞的反光程度与可见程度,可由式(5)和(6)得到:
R=random()×(r1-r2)+r2 (5)
V=random()×(v1-v2)+v2 (6)
式中r1和r2是表示细胞伸缩间隔的两个常数,v1和v2是表示细胞可见度的两个常数,初始化参数时已给出。
步骤3.5.3):C2类中细胞的反射与可见性由式(7)和式(8)确定:
reflectionj=R×Soptj (7)
Figure FDA0002283046180000031
其中
Figure FDA0002283046180000032
表示C2中第i个细胞在第j维上的位置,R和V仍由式(5)(6)得到。
步骤3.5.4):C3类中细胞反射过程由式(7)确定,用(9)式表示其可见性:
Figure FDA0002283046180000033
其中
Figure FDA0002283046180000034
表示Soptj平均值。
步骤3.5.5):C4类由式(10)随机更新个体
Snewj=random()×(Lj-lj)+lj (10)
其中Lj和lj表示在第j维上搜索的上、下界。
对通过上述细胞更新过程得到的Snew,若与当前最优解Sopt对应的目标函数值满足关系
Figure FDA0002283046180000035
则用Snew替换当前最优解。
步骤3.5.6):为提高各类细胞的搜索效率,增强C1、C4类个体的全局搜索能力,改善C2、C3类的局部寻优能力,根据个体适应度自适应调整权重分布:具体如下:
3.5.6.1)对于C1中适应度较差的个体,即目标值小于平均值Javg的个体,应适当的增大其权重,对于大于平均值的个体,采用余弦函数设置权重,以保持其较强的全局搜索能力,并且在后期能进行更细致的搜索。
Figure FDA0002283046180000036
其中wa、wb、和wc为限制权重取值范围的常数,且满足wa>wb>wc,k为当前迭代次数,max_iters为最大迭代次数。
3.5.6.2)C2类对目标值越小的个体权重设置越大,同时为改善种群迭代过程中的多样性,对适应度较好的即接近全局最优解的个体,设置较大的权重。
Figure FDA0002283046180000041
其中Jopt为当前最优解对应的目标值。
3.5.6.3)C3类通过余弦函数更新权重,为提高该类个体的局部搜索能力,随着迭代次数的增大不断缩小搜索范围。
Figure FDA0002283046180000042
其中wd为小于1的常数。
在新一轮的迭代过程中各类细胞生成反射及可见性时需要加入上述权重因子,更新细胞状态的公式变形为
Figure FDA0002283046180000043
其中wi∈{w1,w2,w3}。
CN201911149108.6A 2019-11-21 2019-11-21 基于智能自适应优化算法的高超声速飞行器再入段轨迹优化控制器 Expired - Fee Related CN111103795B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911149108.6A CN111103795B (zh) 2019-11-21 2019-11-21 基于智能自适应优化算法的高超声速飞行器再入段轨迹优化控制器

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911149108.6A CN111103795B (zh) 2019-11-21 2019-11-21 基于智能自适应优化算法的高超声速飞行器再入段轨迹优化控制器

Publications (2)

Publication Number Publication Date
CN111103795A true CN111103795A (zh) 2020-05-05
CN111103795B CN111103795B (zh) 2021-07-13

Family

ID=70421461

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911149108.6A Expired - Fee Related CN111103795B (zh) 2019-11-21 2019-11-21 基于智能自适应优化算法的高超声速飞行器再入段轨迹优化控制器

Country Status (1)

Country Link
CN (1) CN111103795B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111856944A (zh) * 2020-08-05 2020-10-30 重庆大学 一种基于事件触发的高超声速飞行器模糊控制方法
CN113283164A (zh) * 2021-05-17 2021-08-20 北京航空航天大学 一种基于遗传算法的飞行器无动力段性能优化方法
CN114355770A (zh) * 2021-11-30 2022-04-15 浙江大学 入侵杂草种群分布优化的高超声速飞行器再入段控制器
CN114488788A (zh) * 2021-11-30 2022-05-13 浙江大学 概率自适应花授粉优化的高超声速飞行器再入段控制器
CN115167126A (zh) * 2022-06-24 2022-10-11 西北工业大学 一种两级入轨高超声速飞行器上升段纵向轨迹设计和优化方法
CN115308301A (zh) * 2022-08-16 2022-11-08 中山大学 一种可测量细胞及细胞核的弹性模量的测量装置

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103197543A (zh) * 2013-02-25 2013-07-10 西北工业大学 基于运动状态综合识别的高速飞行器自适应控制方法
CN106227972A (zh) * 2016-08-04 2016-12-14 北京航空航天大学 一种高超声速飞行器平稳滑翔弹道的优化方法
CN107491090A (zh) * 2017-08-25 2017-12-19 中国人民解放军火箭军工程大学 基于检测点自适应伪谱法的飞行器再入轨迹快速规划方法
CN107908109A (zh) * 2017-11-13 2018-04-13 浙江大学 一种基于正交配置优化的高超声速飞行器再入段轨迹优化控制器
CN109725644A (zh) * 2019-01-22 2019-05-07 湖南云顶智能科技有限公司 一种高超声速飞行器线性优化控制方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103197543A (zh) * 2013-02-25 2013-07-10 西北工业大学 基于运动状态综合识别的高速飞行器自适应控制方法
CN106227972A (zh) * 2016-08-04 2016-12-14 北京航空航天大学 一种高超声速飞行器平稳滑翔弹道的优化方法
CN107491090A (zh) * 2017-08-25 2017-12-19 中国人民解放军火箭军工程大学 基于检测点自适应伪谱法的飞行器再入轨迹快速规划方法
CN107908109A (zh) * 2017-11-13 2018-04-13 浙江大学 一种基于正交配置优化的高超声速飞行器再入段轨迹优化控制器
CN109725644A (zh) * 2019-01-22 2019-05-07 湖南云顶智能科技有限公司 一种高超声速飞行器线性优化控制方法

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111856944A (zh) * 2020-08-05 2020-10-30 重庆大学 一种基于事件触发的高超声速飞行器模糊控制方法
CN111856944B (zh) * 2020-08-05 2022-01-28 重庆大学 一种基于事件触发的高超声速飞行器模糊控制方法
CN113283164A (zh) * 2021-05-17 2021-08-20 北京航空航天大学 一种基于遗传算法的飞行器无动力段性能优化方法
CN113283164B (zh) * 2021-05-17 2022-07-22 北京航空航天大学 一种基于遗传算法的飞行器无动力段性能优化方法
CN114355770A (zh) * 2021-11-30 2022-04-15 浙江大学 入侵杂草种群分布优化的高超声速飞行器再入段控制器
CN114488788A (zh) * 2021-11-30 2022-05-13 浙江大学 概率自适应花授粉优化的高超声速飞行器再入段控制器
CN115167126A (zh) * 2022-06-24 2022-10-11 西北工业大学 一种两级入轨高超声速飞行器上升段纵向轨迹设计和优化方法
CN115167126B (zh) * 2022-06-24 2024-04-30 西北工业大学 一种两级入轨高超声速飞行器上升段纵向轨迹设计和优化方法
CN115308301A (zh) * 2022-08-16 2022-11-08 中山大学 一种可测量细胞及细胞核的弹性模量的测量装置

Also Published As

Publication number Publication date
CN111103795B (zh) 2021-07-13

Similar Documents

Publication Publication Date Title
CN111103795B (zh) 基于智能自适应优化算法的高超声速飞行器再入段轨迹优化控制器
CN110703766B (zh) 一种基于迁移学习策略深度q网络的无人机路径规划方法
CN109725644A (zh) 一种高超声速飞行器线性优化控制方法
CN107908109A (zh) 一种基于正交配置优化的高超声速飞行器再入段轨迹优化控制器
CN110347181B (zh) 基于能耗的无人机分布式编队控制方法
CN108919818B (zh) 基于混沌种群变异pio的航天器姿态轨道协同规划方法
Liu et al. Novel docking controller for autonomous aerial refueling with probe direct control and learning-based preview method
CN106647264B (zh) 一种基于控制约束的扩展鲁棒h∞的无人机控制方法
CN108394429B (zh) 一种为城轨列车群生成自动驾驶曲线的方法
CN111506114B (zh) 一种飞行器编队控制方法
CN111045447B (zh) 高精度的高超声速飞行器轨迹优化多尺度最优控制系统
CN115903888A (zh) 一种基于天牛群算法的旋翼无人机自主路径规划方法
CN112327926A (zh) 一种无人机编队的自适应滑模控制方法
CN111324035A (zh) 一种高超声速飞行器轨迹优化自适应最优控制器
Xu et al. Coordinated intelligent control of the flight control system and shape change of variable sweep morphing aircraft based on dueling-DQN
CN113900440A (zh) 一种无人机控制律设计方法、装置及可读存储介质
CN110456790A (zh) 基于自适应权重的智能网联电动汽车队列优化控制方法
CN107422637A (zh) 一种集群智能控制系统和方法
Guo et al. PID parameter tuning algorithm of rotor UAV Based on Improved Particle Swarm Optimization
CN114326821B (zh) 基于深度强化学习的无人机自主避障系统及方法
CN114815875B (zh) 一种基于集合满射鸽群智能优化的无人机集群编队控制器调参方法
CN111061294B (zh) 非平稳高超声速飞行器轨迹优化自适应最优控制系统
CN114355770A (zh) 入侵杂草种群分布优化的高超声速飞行器再入段控制器
CN113359852B (zh) 一种仿原鸽个体属性智能行为的无人机集群协同控制方法
CN112698666B (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: 20210713

Termination date: 20211121