CN107092987B - 一种中小型无人机自主着舰风速预测方法 - Google Patents
一种中小型无人机自主着舰风速预测方法 Download PDFInfo
- Publication number
- CN107092987B CN107092987B CN201710253166.8A CN201710253166A CN107092987B CN 107092987 B CN107092987 B CN 107092987B CN 201710253166 A CN201710253166 A CN 201710253166A CN 107092987 B CN107092987 B CN 107092987B
- Authority
- CN
- China
- Prior art keywords
- learning machine
- extreme learning
- network
- machine network
- parameter
- 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.)
- Expired - Fee Related
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N20/00—Machine learning
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Software Systems (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Physics (AREA)
- Strategic Management (AREA)
- Evolutionary Computation (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Artificial Intelligence (AREA)
- Development Economics (AREA)
- Entrepreneurship & Innovation (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Computational Linguistics (AREA)
- Molecular Biology (AREA)
- Game Theory and Decision Science (AREA)
- General Health & Medical Sciences (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Marketing (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Medical Informatics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供了一种基于智能蜂群算法与极限学习机网络相结合的中小型无人机自主着舰风速预测方法,其包括:对智能蜂群算法和极限学习机网络相关参数进行初始化设置,实时采样极限学习机网络输入节点变量然后存储;利用采样数据对极限学习机网络进行训练并采用智能蜂群算法对极限学习机网络参数进行优化;利用已训练优化好的极限学习机网络和当前时刻采样的网络输入节点变量对下一周期的风速进行预测。本发明为中小型无人机安全着舰提供所需的实时风速预测方法,有利于中小型无人机根据预测风速对着舰控制进行自适应调整,提高中小型无人机自主着舰成功率。
Description
技术领域
本发明涉及无人飞行器气象数据预测领域,尤其是一种中小型无人机自主着舰海面风速预测方法。
背景技术
中小型化是舰载无人作战飞机(Unmanned Combat Aerial Vehicle,UCAV)的发展趋势,未来UCAV将大量替代有人舰载飞机完成侦查、打击和搜救等任务。UCAV自主着舰是指UCAV在降落阶段,将机载设备和航母所采集的所有有用信息进行综合处理,以获得信息种类足够多、精度足够高的相对位置信息,使UCAV在没有外界或人的干预下自动控制完成着舰降落的过程。据统计在舰载飞行器的整个飞行阶段,着陆/着舰事故占到总事故的80%左右,对于无人机着舰过程,海面气流扰动成为影响其着舰成功的重要因素之一。中小型UCAV相对于大型UCAV,其体积小、重量轻,虽然机动灵活,但在着舰阶段更易受到海面气流扰动的影响。黄华等在文献《“大气扰动及其对无人机自动着陆影响仿真研究”、系统仿真学报、2009、Vol-21(21)、6821-6824》、《“风场扰动对无人机控制/导航影响效应仿真”、解放军理工大学学报、2012、Vol-13(5)、565-570》中对以上结论进行了仿真验证。因此,UCAV能够对风速进行准确的短期预测、并预先对自身的飞行状态进行调整,对于中小型UCAV的安全着舰具有重要的意义。
李艳晴在文献《“风速时间序列预测算法研究”、北京科技大学博士论文、2015》中指出:目前对于风速的预测方法大致可分为:物理模型预测法、统计预测法和智能预测法。物理模型预测法是考虑风速产生背景,建立气象预报模型对风速变化进行模拟计算,实现风速预测。物理模型预测法因考虑季节、地理等背景因素,其模型较为复杂,需借助大型计算机辅助。以卡尔曼滤波为代表的统计预测方法,虽计算速度较快、实时性强,但对于非线性较强海面气流扰动,其预测准确度相对较低。神经网络等智能预测方法,对于非线性较强系统的预测准确度较高,但需要大量的历史数据进行网络训练,因此其时间消耗过大。另外,上述方法也只针对于固定点位置的风速预测。
由此可见,上述预测方法针对运动中的中小型UCAV着舰风速预测显然存在不足,不适用于实时性较强的UCAV机载设备,研究新的预测准确度高、运算实时性强的风速预测方法成为发展中小型UCAV自主着舰的关键技术之一。
极限学习机(Extreme Learning Machine,ELM)是一种单隐层前馈神经网络,其隐层权值和阈值随机赋值,再通过广义逆矩阵计算求得网络输出权值,尹刚等在文献《“自适应集成极限学习机在故障诊断中的应用”、振动、测试与诊断、2013、Vol-13(5)、897-901》中指出:极限学习机比传统神经网络具有更快的计算速度,以及泛化能力更强。但是极限学习机隐层权值和阈值随机赋值,通常会导致某些节点参数对网络输出影响较小,从而使极限学习机网络性能下降。
智能蜂群算法是受蜂群采蜜时的个体分工、信息交流和相互协作启发,而提出的一种仿生优化方法,秦全德等在文献《“人工蜂群算法研究综述”、智能系统学报、2014、Vol-9(2)、127-135》中指出:智能蜂群算法与传统优化方法相比,其搜索精度高、鲁棒性较强。
因此,本发明结合上述两种算法的优点,设计一种基于智能蜂群算法和极限学习机相结合的动态中小型UCAV着舰风速预测方法。
发明内容
为了克服现有技术的不足,本发明提供一种中小型无人机自主着舰风速预测方法,具体包括:
步骤A:参数初始化及数据采样存储:对智能蜂群算法和极限学习机网络参数进行初始化设置,参数包括:智能蜂群算法的种群数量Csize、迭代次数I、最大迭代次数IMax、误差阈值EGoal、优化参数尺度Dparameter、迭代控制数ILimit、优化参数下限Rlow、优化参数上限Rup;极限学习机的网络输入节点数Ninput、网络隐层节点数Nhide、输出节点数Nout、训练样本数N、输入节点参数采集次数Z、输入节点参数采集周期Tr、网络参数优化条件阈值Num;
对风速预测所需的极限学习机网络输入节点数据进行采集,并按照采集时间的先后顺序编号将其进行存储,极限学习机网络输入节点数据G=[Vwind,Dwind,Vuav,T,P,H],其中:Vwind为UCAV所处空间点实时风速、Dwind为风向、Vuav为UCAV航速、T为大气温度、P为大气压强、H为大气湿度;
对采集到的极限学习机网络输入节点数据进行归一化处理;归一化处理公式为:
Y=2(G-Gmin)/(Gmax-Gmin)-1 (1)
式中:G为采集的极限学习机网络输入节点数据,Y为采集的极限学习机网络输入节点归一化数据,Gmax为G取值最大值,Gmin为G取值最小值;
步骤B:极限学习机网络参数优化:根据当前时刻的前N个采样周期的极限学习机网络输入节点数据作为极限学习机网络训练样本数据,采用智能蜂群算法对极限学习机网络隐层节点权值和阈值进行优化计算,以优化计算得到权值和阈值作为极限学习机最终网络参数;
其中包括以下步骤:
步骤B1:随机生成ELM网络参数初始解;
步骤B2:计算ELM网络隐层节点输出值;
步骤B3:计算ELM网络隐层节点权值;
步骤B4:计算样本风速预测值;
步骤B5:计算ELM网络参数解适应度值;
步骤B6:生成ELM网络参数候选解;
步骤B7:计算ELM网络参数候选解概率;
步骤B8:判断迭代次数是否达到门限,输出最优极限学习机网络参数解向量;
步骤C:风速预测:采用智能蜂群算法优化后的极限学习机网络,结合当前时刻采样的极限学习机网络输入节点数据,对下一采样周期的风速进行在线预测;同时判断当前时刻是否满足极限学习机网络参数优化条件,当满足条件时进行极限学习机网络参数优化计算;
其中包括以下步骤:
步骤C1:利用优化训练后的极限学习网络对风速进行预测;
步骤C2:判断当前输入节点参数采集次数Z是否为极限学习机网络参数优化条件阈值Num的整数倍,若当前输入节点参数采集次数Z是Num的整数倍,则转入步骤B利用智能蜂群算法对极限学习机网络参数进行在线优化;若Z不是Num的整数倍,则转入步骤C1继续对风速进行预测。
进一步的,步骤B实现极限学习机网络参数优化,包括:
步骤B1:随机生成ELM网络参数初始解:利用智能蜂群算法随机生成Csize个种群的网络隐层节点权值和阈值组合向量Xj作为极限学习机网络参数解向量,其中Xj=[Aj,Bj],为极限学习机网络输入节点到网络隐层节点的权值,为极限学习机网络隐层节点阈值,xji为第j个种群中、第i个极限学习机网络参数解,其中j={1,2,…Csize}为种群标号,i={1,2,…Ninput×Nhide+Nhide}为种群中参数标号;
步骤B2:计算ELM网络隐层节点输出值:利用前N个周期的极限学习机网络输入节点归一化数据YN作为训练样本数据,结合极限学习机网络参数解向量Xj,计算每个极限学习机网络参数解向量Xj所对应的网络隐层节点输出值,其计算表达式为:
步骤B3:计算ELM网络隐层节点权值:根据前N个周期的风速实际测量值Ok(k={1,2,…N}),计算极限学习机网络参数解向量Xj所对应的网络隐层节点到输出节点的权值其中βji表示第j个种群中、第i个网络隐层节点到输出节点的权值,网络隐层节点到输出节点的权值计算表达式为:
式中:H+为极限学习机网络隐层节点输出矩阵H的Moore-Penrose增广逆矩阵;
步骤B4:计算样本风速预测值:根据计算得到的极限学习机网络参数解向量Xj所对应的网络隐层节点到输出节点的权值和前N个周期的极限学习机网络输入节点归一化数据YN,计算极限学习机网络参数解向量Xj和第k个极限学习机网络输入节点归一化数据Yk所对应的极限学习机网络输出风速预测值其中k={1,2,…,N},其计算表达式为:
步骤B5:计算ELM网络参数解适应度值:利用适应度公式计算每个极限学习机网络参数解向量Xj的适应度值,并存储其中最大适应度值所对应的极限学习机网络参数解向量到Xbest变量中,Xbest为极限学习机网络参数最优解向量,适应度值计算公式表示为:
式中:fitj为第j个极限学习机网络参数解向量Xj所对应的适应度值;
步骤B6:生成ELM网络参数候选解:根据已有的极限学习机网络参数解向量Xj,应用候选解公式生成新的候选极限学习机网络参数解向量Vj,并应用适应度值计算公式(5)计算新解Vj的适应度值,其中并在Vj和Xj之间择优选择,如果Vj的适应度值大于Xj,则用Vj替换Xj,候选解公式表示为:
vji=xji+φji(xji-xjm) (6)
式中:vji为新的极限学习机网络参数解,φji为[-1,1]之间的随机数,m为不等于i的正整数;
步骤B7:计算ELM网络参数候选解概率:应用选择概率公式计算与Xj相关的选择概率Pj,并以选择概率Pj选择已有的极限学习机网络参数解向量Xj,根据候选解公式(6)进行相邻域搜索产生新解,并应用适应度值计算公式(5)计算新解适应度值,对Vj和Xj进行择优选择,其中选择概率公式为:
式中:Pj为极限学习机网络参数解向量Xj所对应的选择概率;
步骤B8:判断迭代次数是否达到门限,输出最优极限学习机网络参数解向量:根据迭代控制数ILimit,判断智能蜂群算法在ILimit次迭代后,是否发现新的适应度值更优的极限学习机网络参数解向量Xj,如果没有发现新的适应度值更优的极限学习机网络参数解向量Xj,则采用随机搜索公式随机搜索产生一个新解替换旧解;记录已搜索到的极限学习机网络参数解向量Xj中的最优解向量到Xbest变量中,智能蜂群算法迭代次数I=I+1,如果I<IMax,则转至步骤B2;否则,输出最优极限学习机网络参数解向量Xbest;随机搜索公式为:
wji=xjmin+rand(0,1)(xjmax-xjmin) (8)
式中:wji为新产生的第j个种群中、第i个极限学习机网络参数解,xjmin表示第j个种群中值最小的极限学习机网络参数解,xjmax表示第j个种群中值最大的极限学习机网络参数解;
步骤C:实现利用已优化的极限学习机网络对风速进行实时预测,包括:
步骤C1:根据智能蜂群算法迭代输出的最优极限学习机网络参数解向量Xbest和当前时刻采集的极限学习机网络输入节点归一化数据Y,并利用公式(4)对下一周期的风速进行在线预测。
进一步的,步骤A1中智能蜂群算法的种群数量取值范围为5≤Csize≤20、最大迭代次数取值范围为30≤IMax≤200、误差阈值取值范围为0≤EGoal≤5、迭代控制数10≤ILimit<IMax;极限学习机的网络隐层节点数取值范围为5≤Nhide≤50、训练样本数取值范围为N≥20、输入节点参数采集周期取值范围为Tr≥0、网络参数优化条件阈值取值范围为Num>0。
进一步的,步骤A1中智能蜂群算法的种群数量Csize=10、最大迭代次数IMax=50、误差阈值EGoal=5、优化参数尺度Dparameter=70、迭代控制数ILimit=20、优化参数下限Rlow=0、优化参数上限Rup=1;极限学习机的网络输入节点数Ninput=6、输出节点数Nout=1、网络隐层节点数Nhide=10、训练样本数N=100、输入节点参数采集周期Tr=1s。
进一步的,步骤B7中以选择概率Pj选择已有的极限学习机网络参数解向量Xj的方法为:对应每个Pj随机生成一个取值区间为[0,1]的随机数rand(j);如果Pj>randj,则Pj所对应的Xj为选中的极限学习机网络参数解向量。
进一步的,极限学习机网络参数优化条件阈值Num=100。
本发明方法采用智能蜂群算法和极限学习机对UCAV着舰时所处空间风速进行预测,可以根据风速的变化对着舰时飞行姿态进行实时调整,从而减小大气扰动对自主着舰的影响、提高自主着舰的成功率;所采用的极限学习机网络对于风速的非线性模型具有较强的适应性,同时也避免了传统神经网络对大样本数据进行训练耗时的缺点,从而提高了风速预测的实时性要求。
下面结合附图和实施例对本发明作详细说明。
附图说明
图1是实施例中所需预测的UCAV风速变化曲线图;
图2是本发明方法流程图;
图3是实施例中中采用本发明方法极限学习机网络参数优化适应度值变化曲线图;
图4是实施例中采用本发明方法对风速进行预测曲线图;
图5是实施例中采用本发明方法的风速预测误差曲线图;
图6是实施例中采用传统神经网络50次网络训练耗时曲线图;
图7是实施例中采用本发明方法50次网络训练耗时曲线图;
图8是实施例中未采用智能蜂群算法优化的极限学习机网络预测误差曲线图。
具体实施方式
以下将结合具体实施例对本发明提供的技术方案进行详细说明,应理解下述具体实施方式仅用于说明本发明而不用于限制本发明的范围。
本发明一个实施例中UCAV所需预测的风速变化曲线图如图1所示。
本发明在现有技术和需要解决的关键技术问题基础上,提出一种基于智能蜂群算法与极限学习机网络相结合的动态中小型UCAV着舰风速预测方法,其预测方法步骤如图2所示,包括:
步骤A:参数初始化及数据采样存储。对智能蜂群算法和极限学习机网络参数进行初始化设置;对风速预测所需的极限学习机网络输入节点数据进行采集,并按照采集时间的先后顺序编号将其进行存储;对采集到的极限学习机网络输入节点数据进行归一化处理。
步骤B:极限学习机网络参数优化。根据当前时刻的前N个采样周期的极限学习机网络输入节点数据作为极限学习机网络训练样本数据,采用智能蜂群算法对极限学习机网络隐层节点权值和阈值进行优化计算,以优化计算得到权值和阈值作为极限学习机最终网络参数。
步骤C:风速预测。采用智能蜂群算法优化后的极限学习机网络,结合当前时刻采样的极限学习机网络输入节点数据,对下一采样周期的风速进行在线预测;同时判断当前时刻是否满足极限学习机网络参数优化条件,当满足条件时进行极限学习机网络参数优化计算。
其中,步骤A实现参数初始化及数据采样存储。具体实施过程如下:
步骤A1:智能蜂群算法与极限学习机网络参数初始化。对智能蜂群算法和极限学习机的相关参数进行初始化人工设置。其参数包括:智能蜂群算法的种群数量Csize、迭代次数I、最大迭代次数IMax、误差阈值EGoal、优化参数尺度Dparameter、迭代控制数ILimit、优化参数下限Rlow、优化参数上限Rup;极限学习机的网络输入节点数Ninput、网络隐层节点数Nhide、输出节点数Nout、训练样本数N、输入节点参数采集次数Z、输入节点参数采集周期Tr、网络参数优化条件阈值Num。
进一步的,上述步骤中智能蜂群算法的种群数量取值范围为5≤Csize≤20、最大迭代次数取值范围为30≤IMax≤200、误差阈值取值范围为0≤EGoal≤5、迭代控制数10≤ILimit<IMax;极限学习机的网络隐层节点数取值范围为5≤Nhide≤50、训练样本数取值范围为N≥20、输入节点参数采集周期取值范围为Tr≥0、网络参数优化条件阈值取值范围为Num>0。
具体的,在本实施例中以上参数具体取值设定为:智能蜂群算法的种群数量Csize=10、最大迭代次数IMax=50、误差阈值EGoal=5、优化参数尺度Dparameter=70、迭代控制数ILimit=20、优化参数下限Rlow=0、优化参数上限Rup=1;极限学习机的网络输入节点数Ninput=6、输出节点数Nout=1、网络隐层节点数Nhide=10、训练样本数N=100、输入节点参数采集周期Tr=1s。
步骤A2:采集极限学习机网络输入节点数据并按时序存储。采集极限学习机网络输入节点数据(即风速预测所需的相关数据)G,并将其按时间序列保存到存储器中。极限学习机网络输入节点数据G=[Vwind,Dwind,Vuav,T,P,H],其中:Vwind为UCAV所处空间点实时风速、Dwind为风向、Vuav为UCAV航速、T为大气温度、P为大气压强、H为大气湿度。
步骤A3:采集数据归一化处理。对采集的极限学习机网络输入节点数据进行归一化处理,归一化处理公式表示为:
Y=2(G-Gmin)/(Gmax-Gmin)-1 (1)
式中:G为采集的极限学习机网络输入节点数据;Y为采集的极限学习机网络输入节点归一化数据;Gmax为G取值最大值,Gmin为G取值最小值。
步骤B实现极限学习机网络参数优化。具体实施过程如下:
步骤B1:随机生成ELM网络参数初始解。利用智能蜂群算法随机生成Csize个种群的网络隐层节点权值和阈值组合向量Xj作为极限学习机网络参数解向量,其中Xj=[Aj,Bj],为极限学习机网络输入节点到网络隐层节点的权值,为极限学习机网络隐层节点阈值;xji为第j个种群中、第i个极限学习机网络参数解,其中j={1,2,…Csize}为种群标号,i={1,2,…Ninput×Nhide+Nhide}为种群中参数标号。
步骤B2:计算ELM网络隐层节点输出值。利用前N个周期的极限学习机网络输入节点归一化数据YN作为训练样本数据,结合极限学习机网络参数解向量Xj,计算每个极限学习机网络参数解向量Xj所对应的网络隐层节点输出值,其计算表达式为:
步骤B3:计算ELM网络隐层节点权值。根据前N个周期的风速实际测量值Ok(k={1,2,…N}),计算极限学习机网络参数解向量Xj所对应的网络隐层节点到输出节点的权值其中βji表示第j个种群中、第i个网络隐层节点到输出节点的权值。网络隐层节点到输出节点的权值计算表达式为:
式中:H+为极限学习机网络隐层节点输出矩阵H的Moore-Penrose增广逆矩阵。
步骤B4:计算样本风速预测值。根据计算得到的极限学习机网络参数解向量Xj所对应的网络隐层节点到输出节点的权值和前N个周期的极限学习机网络输入节点归一化数据YN,计算极限学习机网络参数解向量Xj和第k个极限学习机网络输入节点归一化数据Yk,其中k={1,2,…,N}所对应的极限学习机网络输出风速预测值其计算表达式为:
步骤B5:计算ELM网络参数解适应度值。利用适应度公式计算每个极限学习机网络参数解向量Xj的适应度值,并存储其中最大适应度值所对应的极限学习机网络参数解向量到Xbest变量中,Xbest为极限学习机网络参数最优解向量。适应度值计算公式表示为:
式中:fitj为第j个极限学习机网络参数解向量Xj所对应的适应度值。
步骤B6:生成ELM网络参数候选解。根据已有的极限学习机网络参数解向量Xj,应用候选解公式生成新的候选极限学习机网络参数解向量Vj,并应用适应度值计算公式(5)计算新解Vj的适应度值,其中并在Vj和Xj之间择优选择,如果Vj的适应度值大于Xj,则用Vj替换Xj。候选解公式表示为:
vji=xji+φji(xji-xjm) (6)
式中:vji为新的极限学习机网络参数解;φji为[-1,1]之间的随机数;m为不等于i的正整数。
步骤B7:计算ELM网络参数候选解概率。应用选择概率公式(7)计算与Xj相关的选择概率Pj,并以选择概率Pj选择已有的极限学习机网络参数解向量Xj,根据候选解公式(6)进行相邻域搜索产生新解,并应用适应度值计算公式(5)计算新解适应度值,对Vj和Xj进行择优选择。其中:
式中:Pj为极限学习机网络参数解向量Xj所对应的选择概率;
进一步的,步骤B7中以选择概率Pj选择已有的极限学习机网络参数解向量Xj的方法描述为:对应每个Pj随机生成一个取值区间为[0,1]的随机数rand(j);如果Pj>randj,则Pj所对应的Xj为选中的极限学习机网络参数解向量。
步骤B8:判断迭代次数是否达到门限,输出最优极限学习机网络参数解向量。根据迭代控制数ILimit,判断智能蜂群算法在ILimit次迭代后,是否发现新的适应度值更优的极限学习机网络参数解向量Xj;如果没有发现新的适应度值更优的极限学习机网络参数解向量Xj,则采用随机搜索公式(8)随机搜索产生一个新解替换旧解。记录已搜索到的极限学习机网络参数解向量Xj中的最优解向量到Xbest变量中,智能蜂群算法迭代次数I=I+1,如果I<IMax,则转至步骤B2;否则,输出最优极限学习机网络参数解向量Xbest。随机搜索公式为:
wji=xjmin+rand(0,1)(xjmax-xjmin) (8)
式中:wji为新产生的第j个种群中、第i个极限学习机网络参数解;xjmin表示第j个种群中值最小的极限学习机网络参数解;xjmax表示第j个种群中值最大的极限学习机网络参数解。
如图3所示,横坐标为智能蜂群算法迭代次数,纵坐标为极限学习机网络参数适应度值。利用智能蜂群算法对极限学习机网络参数进行迭代优化,随着迭代次数的增加极限学习机网络参数的适应度逐渐变强。
步骤C:实现利用已优化的极限学习机网络对风速进行实时预测。具体实施过程如下:
步骤C1:利用优化训练后的极限学习网络对风速进行预测。根据智能蜂群算法迭代输出的最优极限学习机网络参数解向量Xbest和当前时刻采集的极限学习机网络输入节点归一化数据Y,并利用公式(4)对下一周期的风速进行在线预测。
如图4所示,横坐标为时间,纵坐标为风速预测值。采用优化训练后的极限学习机网络对风速进行预测的曲线图,与实际的风速变化趋势相一致。其预测误差如图5所示,图5中横坐标同样为时间,纵坐标为风速预测值与实际测量风速间的差值,其误差相对较小,均值为0.4587m/s。
步骤C2:判断当前输入节点参数采集次数Z是否满足极限学习机网络参数在线优化条件,既当前输入节点参数采集次数Z是否为Num的整数倍。若当前输入节点参数采集次数Z是Num的整数倍,则转入步骤B利用智能蜂群算法对极限学习机网络参数进行在线优化;若Z不是Num的整数倍,则转入步骤C1继续对风速进行预测。
此外,由于硬件条件的限制和算法复杂度的实际,当Num取值过小时,极限学习机网络会因为训练过于频繁,而导致风速预测实时性降低;而当Num取值过大时,又会因为极限学习机网络参数适应度降低,而导致风速预测精度降低。因此在本实施例中具体的极限学习机网络参数优化条件阈值Num100,100为极限学习机网络参数优化条件阈值Num的相对最优取值。
通过采用上述技术方案,本发明具有以下优点:本发明采用智能蜂群算法和极限学习机对UCAV着舰时所处空间风速进行预测,可以根据风速的变化对着舰时飞行姿态进行实时调整,从而减小大气扰动对自主着舰的影响、提高自主着舰的成功率;所采用的极限学习机网络对于风速的非线性模型具有较强的适应性,同时也避免了传统神经网络对大样本数据进行训练耗时的缺点,从而提高了风速预测的实时性要求。图6为传统神经网络对风速进行预测前进行50次网络训练的耗时曲线图,横坐标为训练次数,纵坐标为训练耗时,50次训练的平均耗时为0.6577s。而图7为本发明所用方法对风速进行预测前进行50次网络训练的耗时曲线图,其50次训练的平均耗时为0.2243s。
所采用的智能蜂群算法对极限学习机网络隐含层节点参数进行在线优化,避免了隐含层节点参数随机选取所造成的网络预测结果平均精度不高的缺点,从而提高了网络预测的精度。所采用的步进在线式网络参数训练方法,可以有效提高预测网络对大气扰动时变系统的适应性,从而进一步提高了网络预测精度。图8为未采用智能蜂群算法对极限学习机网络隐含层节点参数进行在线优化时,风速预测误差曲线图。图8中横坐标为时间,纵坐标为风速预测误差值,误差均值为0.8725m/s。其误差均值要高于图5中本发明所采用方法时的误差均值0.4576m/s。
Claims (6)
1.一种中小型无人机自主着舰风速预测方法,包括:
步骤A:参数初始化及数据采样存储:对智能蜂群算法和极限学习机网络参数进行初始化设置,参数包括:智能蜂群算法的种群数量Csize、迭代次数I、最大迭代次数IMax、误差阈值EGoal、优化参数尺度Dparameter、迭代控制数ILimit、优化参数下限Rlow、优化参数上限Rup;极限学习机的网络输入节点数Ninput、网络隐层节点数Nhide、输出节点数Nout、训练样本数N、输入节点参数采集次数Z、输入节点参数采集周期Tr、网络参数优化条件阈值Num;
对风速预测所需的极限学习机网络输入节点数据进行采集,并按照采集时间的先后顺序编号将其进行存储,极限学习机网络输入节点数据G=[Vwind,Dwind,Vuav,T,P,H],其中:Vwind为UCAV所处空间点实时风速、Dwind为风向、Vuav为UCAV航速、T为大气温度、P为大气压强、H为大气湿度;
对采集到的极限学习机网络输入节点数据进行归一化处理;归一化处理公式为:
Y=2(G-Gmin)/(Gmax-Gmin)-1 (1)
式中:G为采集的极限学习机网络输入节点数据,Y为采集的极限学习机网络输入节点归一化数据,Gmax为G取值最大值,Gmin为G取值最小值;
步骤B:极限学习机网络参数优化:根据当前时刻的前N个采样周期的极限学习机网络输入节点数据作为极限学习机网络训练样本数据,采用智能蜂群算法对极限学习机网络隐层节点权值和阈值进行优化计算,以优化计算得到权值和阈值作为极限学习机最终网络参数;
其中包括以下步骤:
步骤B1:随机生成ELM网络参数初始解;
步骤B2:计算ELM网络隐层节点输出值;
步骤B3:计算ELM网络隐层节点权值;
步骤B4:计算样本风速预测值;
步骤B5:计算ELM网络参数解适应度值;
步骤B6:生成ELM网络参数候选解;
步骤B7:计算ELM网络参数候选解概率;
步骤B8:判断迭代次数是否达到门限,输出最优极限学习机网络参数解向量;
步骤C:风速预测:采用智能蜂群算法优化后的极限学习机网络,结合当前时刻采样的极限学习机网络输入节点数据,对下一采样周期的风速进行在线预测;同时判断当前时刻是否满足极限学习机网络参数优化条件,当满足条件时进行极限学习机网络参数优化计算;
其中包括以下步骤:
步骤C1:利用优化训练后的极限学习网络对风速进行预测;
步骤C2:判断当前输入节点参数采集次数Z是否为极限学习机网络参数优化条件阈值Num的整数倍,若当前输入节点参数采集次数Z是Num的整数倍,则转入步骤B利用智能蜂群算法对极限学习机网络参数进行在线优化;若Z不是Num的整数倍,则转入步骤C1继续对风速进行预测。
2.如权利要求1所述的一种中小型无人机自主着舰风速预测方法,其特征在于:
步骤B实现极限学习机网络参数优化,包括:
步骤B1:随机生成ELM网络参数初始解:利用智能蜂群算法随机生成Csize个种群的网络隐层节点权值和阈值组合向量Xj作为极限学习机网络参数解向量,其中Xj=[Aj,Bj],为极限学习机网络输入节点到网络隐层节点的权值,为极限学习机网络隐层节点阈值,xji为第j个种群中、第i个极限学习机网络参数解,其中j={1,2,…Csize}为种群标号,i={1,2,…Ninput×Nhide+Nhide}为种群中参数标号;
步骤B2:计算ELM网络隐层节点输出值:利用前N个周期的极限学习机网络输入节点归一化数据YN作为训练样本数据,结合极限学习机网络参数解向量Xj,计算每个极限学习机网络参数解向量Xj所对应的网络隐层节点输出值,其计算表达式为:
步骤B3:计算ELM网络隐层节点权值:根据前N个周期的风速实际测量值Ok(k={1,2,…N}),计算极限学习机网络参数解向量Xj所对应的网络隐层节点到输出节点的权值其中βji表示第j个种群中、第i个网络隐层节点到输出节点的权值,网络隐层节点到输出节点的权值计算表达式为:
式中:H+为极限学习机网络隐层节点输出矩阵H的Moore-Penrose增广逆矩阵;
步骤B4:计算样本风速预测值:根据计算得到的极限学习机网络参数解向量Xj所对应的网络隐层节点到输出节点的权值和前N个周期的极限学习机网络输入节点归一化数据YN,计算极限学习机网络参数解向量Xj和第k个极限学习机网络输入节点归一化数据Yk所对应的极限学习机网络输出风速预测值其中k={1,2,…,N},其计算表达式为:
步骤B5:计算ELM网络参数解适应度值:利用适应度公式计算每个极限学习机网络参数解向量Xj的适应度值,并存储其中最大适应度值所对应的极限学习机网络参数解向量到Xbest变量中,Xbest为极限学习机网络参数最优解向量,适应度值计算公式表示为:
式中:fitj为第j个极限学习机网络参数解向量Xj所对应的适应度值;
步骤B6:生成ELM网络参数候选解:根据已有的极限学习机网络参数解向量Xj,应用候选解公式生成新的候选极限学习机网络参数解向量Vj,并应用适应度值计算公式(5)计算新解Vj的适应度值,其中并在Vj和Xj之间择优选择,如果Vj的适应度值大于Xj,则用Vj替换Xj,候选解公式表示为:
vji=xji+φji(xji-xjm) (6)
式中:vji为新的极限学习机网络参数解,φji为[-1,1]之间的随机数,m为不等于i的正整数;
步骤B7:计算ELM网络参数候选解概率:应用选择概率公式计算与Xj相关的选择概率Pj,并以选择概率Pj选择已有的极限学习机网络参数解向量Xj,根据候选解公式(6)进行相邻域搜索产生新解,并应用适应度值计算公式(5)计算新解适应度值,对Vj和Xj进行择优选择,其中选择概率公式为:
式中:Pj为极限学习机网络参数解向量Xj所对应的选择概率;
步骤B8:判断迭代次数是否达到门限,输出最优极限学习机网络参数解向量:根据迭代控制数ILimit,判断智能蜂群算法在ILimit次迭代后,是否发现新的适应度值更优的极限学习机网络参数解向量Xj,如果没有发现新的适应度值更优的极限学习机网络参数解向量Xj,则采用随机搜索公式随机搜索产生一个新解替换旧解;记录已搜索到的极限学习机网络参数解向量Xj中的最优解向量到Xbest变量中,智能蜂群算法迭代次数I=I+1,如果I<IMax,则转至步骤B2;否则,输出最优极限学习机网络参数解向量Xbest;随机搜索公式为:
wji=xj min+rand(0,1)(xj max-xj min) (8)
式中:wji为新产生的第j个种群中、第i个极限学习机网络参数解,xj min表示第j个种群中值最小的极限学习机网络参数解,xj max表示第j个种群中值最大的极限学习机网络参数解;
步骤C:实现利用已优化的极限学习机网络对风速进行实时预测,包括:
步骤C1:根据智能蜂群算法迭代输出的最优极限学习机网络参数解向量Xbest和当前时刻采集的极限学习机网络输入节点归一化数据Y,并利用公式(4)对下一周期的风速进行在线预测。
3.如权利要求2所述的一种中小型无人机自主着舰风速预测方法,其特征在于:
步骤A1中智能蜂群算法的种群数量取值范围为5≤Csize≤20、最大迭代次数取值范围为30≤IMax≤200、误差阈值取值范围为0≤EGoal≤5、迭代控制数10≤ILimit<IMax;极限学习机的网络隐层节点数取值范围为5≤Nhide≤50、训练样本数取值范围为N≥20、输入节点参数采集周期取值范围为Tr≥0、网络参数优化条件阈值取值范围为Num>0。
4.如权利要求3所述的一种中小型无人机自主着舰风速预测方法,其特征在于:
步骤A1中智能蜂群算法的种群数量Csize=10、最大迭代次数IMax=50、误差阈值EGoal=5、优化参数尺度Dparameter=70、迭代控制数ILimit=20、优化参数下限Rlow=0、优化参数上限Rup=1;极限学习机的网络输入节点数Ninput=6、输出节点数Nout=1、网络隐层节点数Nhide=10、训练样本数N=100、输入节点参数采集周期Tr=1s。
5.如权利要求2所述的一种中小型无人机自主着舰风速预测方法,其特征在于:
步骤B7中以选择概率Pj选择已有的极限学习机网络参数解向量Xj的方法为:对应每个Pj随机生成一个取值区间为[0,1]的随机数rand(j);如果Pj>randj,则Pj所对应的Xj为选中的极限学习机网络参数解向量。
6.如权利要求2至5中任意一个权利要求所述的一种中小型无人机自主着舰风速预测方法,其特征在于:极限学习机网络参数优化条件阈值Num=100。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710253166.8A CN107092987B (zh) | 2017-04-18 | 2017-04-18 | 一种中小型无人机自主着舰风速预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710253166.8A CN107092987B (zh) | 2017-04-18 | 2017-04-18 | 一种中小型无人机自主着舰风速预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107092987A CN107092987A (zh) | 2017-08-25 |
CN107092987B true CN107092987B (zh) | 2020-11-17 |
Family
ID=59636853
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710253166.8A Expired - Fee Related CN107092987B (zh) | 2017-04-18 | 2017-04-18 | 一种中小型无人机自主着舰风速预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107092987B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107609309B (zh) * | 2017-10-13 | 2021-12-17 | 西安思缪智能科技有限公司 | 一种基于防脉冲平均滤波法和一阶滞后滤波法的海洋风速仿真方法 |
CN108710288A (zh) * | 2018-04-19 | 2018-10-26 | 东华大学 | 基于风速预测的旋翼飞行器空中悬停防漂移的控制方法 |
CN109034478B (zh) * | 2018-07-27 | 2021-09-24 | 中南大学 | 一种高速铁路沿线大风迭代竞争高精度预测方法 |
CN109389253B (zh) * | 2018-11-09 | 2022-04-15 | 国网四川省电力公司电力科学研究院 | 一种基于可信性集成学习的电力系统扰动后频率预测方法 |
CN113051723B (zh) * | 2021-03-12 | 2022-02-18 | 南京航空航天大学 | 一种基于时序网络的蜂群无人机故障传播分析方法 |
CN113539386A (zh) * | 2021-06-30 | 2021-10-22 | 淮阴工学院 | 基于clmvo-elm的溶解氧浓度预测方法、装置、设备及存储介质 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106022465A (zh) * | 2016-05-19 | 2016-10-12 | 江南大学 | 改进人工蜂群优化的极限学习机方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8301406B2 (en) * | 2008-07-24 | 2012-10-30 | University Of Cincinnati | Methods for prognosing mechanical systems |
CN105913151A (zh) * | 2016-04-12 | 2016-08-31 | 河海大学常州校区 | 基于自适应变异粒子群与bp网络光伏电站发电量预测方法 |
-
2017
- 2017-04-18 CN CN201710253166.8A patent/CN107092987B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106022465A (zh) * | 2016-05-19 | 2016-10-12 | 江南大学 | 改进人工蜂群优化的极限学习机方法 |
Non-Patent Citations (5)
Title |
---|
基于极限学习机的脉动风速快速预测方法;李春祥等;《上海交通大学学报》;20161130;第50卷(第11期);第1719-1721页 * |
大气扰动及其对无人机自动着陆影响仿真研究;黄华等;《系统仿真学报》;20091130;第21卷(第21期);第6821-6824页 * |
自适应集成极限学习机在故障诊断中的应用;尹刚等;《振动、测试与诊断》;20131031;第33卷(第5期);第897-918页 * |
采用正则化极限学习机的短期风速预测;袁翀等;《电网与清洁能源》;20161130;第32卷(第11期);第62-68页 * |
风场扰动对无人机控制/导航影响效应仿真;黄华等;《解放军理工大学学报(自然科学版)》;20121030;第13卷(第5期);第566-570页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107092987A (zh) | 2017-08-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107092987B (zh) | 一种中小型无人机自主着舰风速预测方法 | |
CN112256056B (zh) | 基于多智能体深度强化学习的无人机控制方法及系统 | |
CN110580496A (zh) | 一种基于熵最小化的深度迁移学习系统及方法 | |
CN110806759B (zh) | 一种基于深度强化学习的飞行器航线跟踪方法 | |
CN107944648B (zh) | 一种大型船舶航速油耗率预测方法 | |
CN111047085B (zh) | 一种基于元学习的混合动力车辆工况预测方法 | |
CN106951695A (zh) | 多工况下的机械设备剩余使用寿命计算方法及系统 | |
CN110244689A (zh) | 一种基于判别性特征学习方法的auv自适应故障诊断方法 | |
CN106056127A (zh) | 一种带模型更新的gpr在线软测量方法 | |
CN111310965A (zh) | 一种基于lstm网络的飞行器航迹预测方法 | |
CN109255726A (zh) | 一种混合智能技术的超短期风功率预测方法 | |
CN112784920A (zh) | 云边端协同的旋转部件对抗域自适应故障诊断方法 | |
CN112215412B (zh) | 溶解氧预测方法及装置 | |
CN114815882A (zh) | 一种基于强化学习的无人飞行器自主编队智能控制方法 | |
Tan et al. | Support vector machine algorithm for artificial intelligence optimization | |
CN109599866A (zh) | 一种预测辅助的电力系统状态估计方法 | |
CN114528638A (zh) | 基于强化学习的船舶运动大多步实时预测混合方法及系统 | |
CN117113870A (zh) | 一种考虑风干扰的固定翼飞行器气动力建模方法 | |
CN107168066A (zh) | 一种温室环境自适应控制方法 | |
CN114384931B (zh) | 一种基于策略梯度的无人机多目标最优控制方法和设备 | |
CN116565876A (zh) | 一种鲁棒强化学习的配网潮流优化方法及计算机可读介质 | |
Priandana et al. | Development of self-organizing maps neural networks based control system for a boat model | |
Yu et al. | Online ship rolling prediction using an improved OS-ELM | |
CN108960406B (zh) | 一种基于bfo小波神经网络的mems陀螺随机误差预测方法 | |
Na et al. | Modified particle swarm optimization based algorithm for BP neural network for measuring aircraft remaining fuel volume |
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: 20201117 |