[发明专利]基于自适应变异果蝇优化算法优化RBF神经网络的海杂波最优软测量仪表及方法

权利要求书

1.一种基于自适应变异果蝇优化算法优化RBF神经网络的海杂波最优软测量仪表，所述的最优软测量仪表，包括雷达、用于测量易测变量的现场智能仪表、用于测量操作变量的控制站、存放数据的现场数据库以及海杂波软测量值显示仪；所述现场智能仪表、控制站与雷达连接，所述现场智能仪表、控制站与现场数据库连接，其特征在于：所述软测量仪表还包括基于自适应变异果蝇优化算法优化RBF神经网络的最优软测量上位机，所述现场数据库与所述基于自适应变异果蝇优化算法优化RBF神经网络的最优软测量上位机的输入端连接，所述基于自适应变异果蝇优化算法优化RBF神经网络的最优软测量上位机的输出端与海杂波软测量值显示仪连接；所述基于自适应变异果蝇优化算法优化RBF神经网络的最优软测量上位机包括：

数据预处理模块，用于将从现场数据库输入的模型训练样本进行预处理，对训练样本中心化，即减去样本的平均值，然后对其进行标准化：

计算均值：

计算方差：

标准化：

其中，TX为训练样本，N为训练样本数，为训练样本的均值，X为标准化后的训练样本，σ_x为为计算方差。

采用前向网络拓扑结构，隐含层的每个单元输出h为:

hi=Ri(X)=Ri(||X-C||/σi)i=1,2,3,...,L---(4)]]>

其中，X是N维输入向量，C_i是与X同维数的向量，R_i为径向基核函数。

R_i(·)通常采用高斯函数，即：

Ri(X)=exp(-||X-Ci||2/σi2)---(5)]]>

因此网络输出层的输出y为

y=Σi=1lwjh(i)---(6)]]>

其中w是连接权重，h是隐含层的单元输出。

自适应变异果蝇优化算法模块，用于采用自适应变异果蝇优化算法对RBF神经网络的核函数宽度进行优化，具体步骤如下：

①确定自适应变异果蝇优化算法的优化参数为RBF神经网络模块的平移因子和伸缩因子、粒子群个体数目popsize、最大循环寻优次数iter_max、第p个粒子的初始位置区域X_axis,Y_axis。

②设定优化目标函数，将其转换为适应度，通过相应的误差函数计算适应度函数，并认为误差大的粒子适应度小，粒子p的适应度函数f表示为：

f_p＝1/(E_p+1)(7)

式中，E_p是RBF神经网络模型的误差函数，表示为：

Ep=1NΣi=1N(y^i-Oi)2---(8)]]>

式中，是RBF神经网络模型的预测输出，O_i为RBF神经网络模型的目标输出；N为训练样本数；

③按照如下公式，粒子进行搜索，

Xi=X_axis+Random valueYi=Y_axis+Random value---(9)]]>

式中，RandomValue为搜索距离；

④对于粒子p，预先估计与原点的距离Dist，再计算味道浓度判定值S，该值为距离倒数：

Dist_i＝(X_i²+Y_i²)^1/2 (10)

S_i＝1/Dist_i(11)

⑤将味道浓度判定值S_i代人味道浓度判定函数(或称为适应度函数fitness function)，用来求出果蝇个体位置的味道浓度Smell_i:

[bestSmellbestIndex]＝min(Smell) (12)

⑥记录最优个体位置及味道浓度值，此时所有果蝇个体利用视觉将向这个位置飞去：

Smellbest=bestSmellX_axis=X(bestIndex)Y_axis=Y(bestIndex)---(13)]]>

⑦按照公式(14)计算该果蝇群体的平均味道浓Smell_ave，然后按照公式(15)计算果蝇群体适应度方差τ²

Smellave=Σi=1sizepopSmelli/sizepop---(14)]]>

τ2=Σi=1sizepop(Smelli-Smellave)2---(15)]]>

⑧假如τ²≤δ₁并且Smellbest>δ₂或者产生分布于[0，1]之间的随机数r<P，则首先按照式(16)复制个最优果蝇个体(X_axism_j,Y_axism_j)(j＝1,2,…,M)；其次按照式(17)对复制的最优果蝇个体进行变异，将复制的最优果蝇个体位置更新为新位置(X_axism_j,Y_axism_j)(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_axism_j,Y_axism_j)与原点的距离Dist'，再按照公式计算新位置味道浓度判定值S':

Distj′=(X_axismj′2+Y_axismj′2)1/2Sj′=1/Distj′---(18)]]>

⑩将再次计算出的味道浓度判定值S'代入味道浓度判定函数，求出新位置的味道浓度Smell‘。

Smell_i'＝Function(S_i') (19)

假如满足Smell_i'<Smellbest的条件，那么Smellbest＝Smell'_j,X_axis＝X_axism'_j,Y_axis＝Y_axism'_j,(j＝1,2,…,M).

判断是否满足性能要求，若是，结束寻优，得到一组优化的RBF神经网络的参数；否则返回步骤③，继续迭代寻优，直至达到最大迭代次数iter_max。

模型更新模块，用于模型的在线更新，定期将离线化验数据输入到训练集中，更新RBF神经网络模型。

2.一种如上所述的基于自适应变异果蝇优化算法优化RBF神经网络的海杂波最优软测量仪表实现的软测量方法，其特征在于：所述软测量方法包括以下步骤：

1)、对雷达对象，根据特性分析和气候分析，选择操作变量和易测变量作为模型的输入，操作变量和易测变量由现场数据库获得；

2)、将从现场数据库输入的模型训练样本进行预处理，对训练样本中心化，即减去样本的平均值，然后对其进行标准化，使得其均值为0，方差为1。该处理采用以下算式过程来完成：

2.1)计算均值：

2.2)计算方差：

2.3)标准化：

其中，TX为训练样本，N为训练样本数，为训练样本的均值，X为标准化后的训练样本，σ_x为计算方差。

3)对从数据预处理模块传过来的训练样本，采用RBF神经网络进行建。RBF神经网络用径向基函数作为隐单元的“基”构成隐含层空间，将输人矢量直接映射到隐含层空间.隐含层空间到输出空间的映射是线性的。从总体上看，网络由输入到输出的映射是非线性的，而网络输出对可调参数而言却又是线性的。理论已经证明，RBF网络具有全局和最佳逼近性能，是前馈神经网络中完成映射功能最优的网络。只要有足够多的隐层神经元，RBF能以任意精度近似任何连续函数。采用前向网络拓扑结构，隐含层的每个单元输出h为:

hi=Ri(X)=Ri(||X-C||/σi)i=1,2,3,...,L---(4)]]>

其中，X是N维输入向量，C_i是与X同维数的向量，R_i为径向基核函数。

R_i(·)通常采用高斯函数，即：

Ri(X)=exp(-||X-Ci||2/σi2)---(5)]]>

因此网络输出层的输出y为

y=Σi=1lwjh(i)---(6)]]>

其中w是连接权重，h是隐含层的单元输出。

4)、采用自适应变异果蝇优化算法对RBF神经网络的核函数宽度进行优化，具体步骤如下：

①确定自适应变异果蝇优化算法的优化参数为RBF神经网络模块的平移因子和伸缩因子、粒子群个体数目popsize、最大循环寻优次数iter_max、第p个粒子的初始位置区域X_axis,Y_axis。

②设定优化目标函数，将其转换为适应度，通过相应的误差函数计算适应度函数，并认为误差大的粒子适应度小，粒子p的适应度函数f表示为：

f_p＝1/(E_p+1)(7)

式中，E_p是RBF神经网络模型的误差函数，表示为：

Ep=1NΣi=1N(y^i-Oi)2---(8)]]>

式中，是RBF神经网络模型的预测输出，O_i为RBF神经网络模型的目标输出；N为训练样本数；

③按照如下公式，粒子进行搜索，

Xi=X_axis+Random valueYi=Y_axis+Random value---(9)]]>

式中，RandomValue为搜索距离；

④对于粒子p，预先估计与原点的距离Dist，再计算味道浓度判定值S，该值为距离倒数：

Dist_i＝(X_i²+Y_i²)^1/2 (10)

S_i＝1/Dist_i(11)

⑤将味道浓度判定值S_i代人味道浓度判定函数(或称为适应度函数fitness function)，用来求出果蝇个体位置的味道浓度Smell_i:

[bestSmell bestIndex]＝min(Smell) (12)

⑥记录最优个体位置及味道浓度值，此时所有果蝇个体利用视觉将向这个位置飞去：

Smellbest=bestSmellX_axis=X(bestIndex)Y_axis=Y(bestIndex)---(13)]]>

⑦按照公式(14)计算该果蝇群体的平均味道浓Smell_ave，然后按照公式(15)计算果蝇群体适应度方差τ²

Smellave=Σi=1sizepopSmelli/sizepop---(14)]]>

τ2=Σi=1sizepop(Smelli-Smellave)2---(15)]]>

⑧假如τ²≤δ₁并且Smellbest>δ₂或者产生分布于[0，1]之间的随机数r<P，则首先按照式(16)复制个最优果蝇个体(X_axism_j,Y_axism_j)(j＝1,2,…,M)；其次按照式(17)对复制的最优果蝇个体进行变异，将复制的最优果蝇个体位置更新为新位置(X_axism_j,Y_axism_j)(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_axism_j,Y_axism_j)与原点的距离Dist'，再按照公式计算新位置味道浓度判定值S':

Distj′=(X_axismj′2+Y_axismj′2)1/2Sj′=1/Distj′---(18)]]>

⑩将再次计算出的味道浓度判定值S'代入味道浓度判定函数，求出新位置的味道浓度Smell‘。

Smell_i'＝Function(S_i') (19)

假如满足Smell_i'<Smellbest的条件，那么Smellbest＝Smell'_j,X_axis＝X_axism'_j,Y_axis＝Y_axism'_j,(j＝1,2,…,M).

判断是否满足性能要求，若是，结束寻优，得到一组优化的RBF神经网络的参数；否则返回步骤③，继续迭代寻优，直至达到最大迭代次数iter_max。

定期将离线化验数据输入到训练集中，更新RBF神经网络模型。