[发明专利]基于自适应果蝇优化算法优化最小二乘支持向量机的海杂波最优软测量仪表及方法

摘要

本发明公开了一种基于自适应变异果蝇优化最小二成支持向量机的海杂波最优软测量仪表及方法，包括雷达、现场智能仪表、控制站、存放数据的现场数据库、基于改进的果蝇优化算法优化小波神经网络的最优软测量上位机以及预报软测量值显示仪，现场智能仪表及控制站与雷达相连，与现场数据库相连；最优软测量上位机与现场数据库及软测量值显示仪相连。所述的基于自适应变异果蝇优化最小二成支持向量机的最优软测量上位机，包括数据预处理模块、小波神经网络模块、模型更新模块。本发明实现了海杂波的在线最优软测量，克服了人为因素造成的随机性影响，提高了模型预报的稳定性，降低了模型预报陷入局部最优的可能性。

主权项

一种基于自适应果蝇优化算法优化最小二乘支持向量机的海杂波最优软测量仪表，所述的最优软测量仪表，包括雷达、用于测量易测变量的现场智能仪表、用于测量操作变量的控制站、存放数据的现场数据库以及海杂波软测量值显示仪；所述现场智能仪表、控制站与雷达连接，所述现场智能仪表、控制站与现场数据库连接，其特征在于：所述软测量仪表还包括基于自适应果蝇优化算法优化最小二乘支持向量机的最优软测量上位机，所述现场数据库与所述基于自适应果蝇优化算法优化最小二乘支持向量机的最优软测量上位机的输入端连接，所述基于自适应果蝇优化算法优化最小二乘支持向量机的最优软测量上位机的输出端与海杂波软测量值显示仪连接；所述基于自适应果蝇优化算法优化最小二乘支持向量机的最优软测量上位机包括：数据预处理模块，用于将从现场数据库输入的模型训练样本进行预处理，对训练样本中心化，即减去样本的平均值，然后对其进行标准化：计算均值：计算方差：标准化：其中，TX为训练样本，N为训练样本数，为训练样本的均值，X为标准化后的训练样本，σx为计算方差。最小二乘支持向量机模块，采用最小二乘支持向量机进行建模。设模型第i个标准化后的训练样本Xi的目标输出为Oi，最小二乘支持向量机通过变换把拟合问题等价于如下二次规划问题：minR(w,ξ)=12wTw+12γΣi=1Nωiξi2---(4)]]>其中，R(w,ξ)是优化问题的目标函数，是非线性映射函数，N是训练样本数，ξ＝{ξ1,…,ξN}是松弛变量，w是最小二乘支持向量机超平面的法向量，b是相应的偏移量，而ωi,i＝1,…,N和γ分别是最小二乘支持向量机的权重和惩罚因子，其中是松弛变量ξi标准差的估计，常量c1,c2通常取为c1＝2.5，c2＝3，由此可得第i个标准化后的训练样本Xi的输出为：y^i=Σm=1Nαm×K<Xm,Xi>+b---(7)]]>其中，K<·>是最小二乘支持向量机的核函数，这里K<·>取线性核函数；αm，m＝1,…,N是对应的拉格朗日乘子。自适应果蝇优化算法模块，用于采用自适应果蝇优化算法对最小二乘支持向量机的惩罚因子和误差容限值进行优化，具体步骤如下：①确定自适应果蝇优化算法的优化参数为最小二乘支持向量机模块的惩罚因子和误差容限值、粒子群个体数目popsize、最大循环寻优次数itermax、群体适应度方差阈值δ1，理论最优值δ2、第p个粒子的初始位置区域X_axis,Y_axis。②设定优化目标函数，将其转换为适应度，通过相应的误差函数计算适应度函数，并认为误差大的粒子适应度小，粒子p的适应度函数f表示为：fp＝1/(Ep+1)   (8)式中，Ep是最小二乘支持向量机模型的误差函数，表示为：Ep=1NΣi=1N(y^i-Oi)2---(9)]]>式中，是最小二乘支持向量机模型的预测输出，Oi为最小二乘支持向量机模型的目标输出；N为训练样本数；③按照如下公式，粒子进行搜索，Xi=X_axis+RandomvalueYi=Y_axis+Randomvalue---(10)]]>式中，RandomValue为搜索距离；④对于粒子p，预先估计与原点的距离Dist，再计算味道浓度判定值S，该值为距离倒数：Disti＝(Xi2+Yi2)1/2   (11)Si＝1/Disti   (12)⑤将味道浓度判定值Si代人味道浓度判定函数(或称为适应度函数fitness function)，用来求出果蝇个体位置的味道浓度Smelli:Smellbest=bestSmellX_axism=X(bestIndex)Y_axism=Y(bestIndex)---(13)]]>⑥按照公式(14)计算该果蝇群体的平均味道浓Smellave，然后按照公式(15)计算果蝇群体适应度方差τ2Smellave=Σi=1sizepopSmelli/sizepop---(14)]]>τ2=Σi=1sizepop(Smelli-Smellave)2---(15)]]>⑦假如τ2≤δ1并且Smellbest>δ2或者产生分布于[0，1]之间的随机数r<P，则首先按照式(10)复制个最优果蝇个体(X_axismj,Y_axismj)(j＝1,2,…,M)；其次按照式(11)对复制的最优果蝇个体进行变异，将复制的最优果蝇个体位置更新为新位置(X_axismj,Y_axismj)(j＝1,2,…,M):X_axismj=X_axismj=1,2,...,MY_axismj=Y_axismj=1,2,...,M---(16)]]>X_axismj′=X_axism+10*rand()j=1,2,...,MY_axismj′=Y_axism+10*rand()j=1,2,...,M---(17)]]>⑧按照下式先再次估计新位置(X_axismj,Y_axismj)与原点的距离Dist'，再按照公式计算新位置味道浓度判定值S':Distj′=(X_axismj′2+Y_axismj′2)1/2Sj′=1/Distj′---(18)]]>⑨将再次计算出的味道浓度判定值S'代入味道浓度判定函数，求出新位置的味道浓度Smell‘。Smelli'＝Function(Si')   (19)⑩假如满足Smelli'<Smellbest的条件，那么Smellbest＝Smell'j,X_axis＝X_axism'j,Y_axis＝Y_axism'j,(j＝1,2,…,M).判断是否满足性能要求，若是，结束寻优，得到一组优化的最小二乘支持向量机的参数；否则返回步骤⑤，继续迭代寻优，直至达到最大迭代次数itermax。模型更新模块，用于模型的在线更新，定期将离线化验数据输入到训练集中，更新最小二乘支持向量机模型。