多策略人工蜂群优化的水体溶解氧浓度预测方法与流程

技术特征：

1.一种多策略人工蜂群优化的水体溶解氧浓度预测方法，其特征在于：包括以下步骤：

步骤1，在需要预测区域内连续LD天采集水样，然后检测采集水样的水质指标：温度、亚硝酸盐、氨氮、总氮、硝酸盐、溶解氧，并将采集到的水质指标数据作为样本数据集；然后归一化处理所采集到的水质指标样本数据集，并将前70％设置为支持向量机的训练数据集，后30％设置为测试数据集；

步骤2，用户初始化预测跨度天数ND，种群规模Popsize，个体最大未改善次数Limit，最大评价次数MAX_FEs；

步骤3，当前演化代数t＝0，当前评价次数FEs＝0，并令支持向量机的优化设计参数个数D＝3；

步骤4，随机产生初始种群其中：个体下标i＝1,2,...,Popsize，并且表示种群P_t中的第i个个体，其随机产生公式为：

$B_{i,j}^t={\mathrm{LO}}_j+rand(0,1)\×({\mathrm{UP}}_j-{\mathrm{LO}}_j);$

其中维数下标j＝1,2,3；表示了支持向量机的3个优化设计参数的值，即是支持向量机的惩罚因子C，是支持向量机的径向基核参数g，是支持向量机的不敏感损失函数中参数ε；rand(0,1)是在[0,1]之间产生随机实数的函数，LO_j和UP_j分别为支持向量机的第j个优化设计参数的搜索区间的下界和上界；

步骤5，计算种群P_t中每个个体的适应值其中个体下标i＝1,2,...,Popsize，个体的适应值的计算方法为：利用个体作为支持向量机的训练参数，并用训练数据集训练支持向量机，其中支持向量机的输入变量为归一化的一天的水质指标数据：温度、亚硝酸盐、氨氮、总氮、硝酸盐、溶解氧；支持向量机的输出为归一化的ND天后的水体溶解氧浓度值；然后计算训练好的支持向量机在测试数据集上的均方误差NE_i，再令个体的适应值

步骤6，令种群P_t中所有个体的未改善次数其中个体下标i＝1,2,...,Popsize，并令当前评价次数FEs＝FEs+Popsize；

步骤7，保存种群P^t中的最优个体Best^t；

步骤8，令策略因子STU＝0.6；

步骤9，雇佣蜂执行其搜索操作；

步骤10，根据种群P_t中个体的适应值计算每个个体的选择概率；

步骤11，观察蜂根据种群P_t中每个个体的选择概率选择出个体执行适应性多策略搜索操作生成新个体，然后执行选择操作并更新个体的未改善次数，具体步骤如下：

步骤11.1，令计数器i＝1，令策略因子列表UList为空；

步骤11.2，根据种群P_t中每个个体的选择概率采用轮盘赌策略选择出个体并令新个体

步骤11.3，在[1,D]之间随机产生一个正整数RJ；

步骤11.4，产生一个高斯随机数GU＝normrnd(STU,0.2)，如果GU的值超出了[0,1]之间的范围，则采用同样的方法重新产生高斯随机实数GU，直到GU的值不超出[0,1]之间的范围，其中normrnd(STU,0.2)表示为以STU为均值，0.2为标准差的高斯随机数产生函数；

步骤11.5，令策略序号STI＝floor(GU×3)％3，其中floor为向下取整函数，％为取余运算符号；

步骤11.6，如果STI等于0则转到步骤11.7；否则判断STI是否等于1，如果STI等于1，则转到步骤11.12，否则转到步骤11.15；

步骤11.7，在[1,Popsize]之间随机产生一个整数RI1，并令随机权值RK＝rand(0,1)；

步骤11.8，令均值

步骤11.9，令标准差

步骤11.10，以MeanU为均值，STDV为标准差产生一个高斯随机实数Val，如果Val的值超出了[LO_RJ,UP_RJ]之间的范围，则采用同样的方法重新产生高斯随机实数Val，直到Val的值不超出[LO_RJ,UP_RJ]之间的范围，然后令

步骤11.11，转到步骤11.16；

步骤11.12，在[1,Popsize]之间随机产生两个不相等的整数RI2和RI3；

步骤11.13，令

步骤11.14，转到步骤11.16；

步骤11.15，令

步骤11.16，计算新个体U^t的适应值Fit(U^t)；

步骤11.17，如果新个体U^t比个体更优，则将GU添加到策略因子列表UList中；

步骤11.18，在个体与新个体U^t之间执行选择操作并更新个体的未改善次数

步骤11.19，令计数器i＝i+1；

步骤11.20，如果计数器i小于或等于Popsize，则转到步骤11.2，否则转到步骤11.21；

步骤11.21，计算策略因子列表UList中数据的平均值MeanLU；

步骤11.22，令STU＝0.9×STU+0.1×MeanLU，然后转到步骤12；

步骤12，令当前评价次数FEs＝FEs+Popsize×2；

步骤13，侦察蜂执行其搜索操作；

步骤14，保存种群P^t中的最优个体Best^t；

步骤15，当前演化代数t＝t+1；

步骤16，重复步骤9至步骤15直至当前评价次数FEs达到MAX_FEs后结束，将执行过程中得到的最优个体Best^t作为支持向量机的训练参数，并用训练数据集来训练支持向量机，然后将归一化的一天的水质指标数据：温度、亚硝酸盐、氨氮、总氮、硝酸盐、溶解氧输入到训练好的支持向量机，计算支持向量机的输出即可预测ND天后的水体溶解氧浓度值。