1.一种高斯人工蜂群优化的稀土矿区地下水氨氮浓度预测方法,其特征在于:包括以下步骤:步骤1,在需要预测的稀土矿区内连续MD天采集地下水样,并检测采集到地下水样的水质指标:水温、pH值、亚硝酸盐氮、硝酸盐氮、总氮、溶解氧、五日生化需氧量、氨氮,将采集到的地下水质指标数据作为样本数据集;然后归一化处理所采集到的稀土矿区地下水质指标样本数据集,并将前70%设置为支持向量机的训练数据集,后30%设置为测试数据集;
步骤2,用户初始化参数,所述初始化参数包括预测跨度天数KD,种群大小Popsize,最大未更新次数Limit,邻域半径NK,最大评价次数MAX_FEs;
步骤3,当前演化代数t=0,当前评价次数FEs=0,并令支持向量机的优化设计参数个数D=3;
步骤4,随机初始化种群 其中:个体下标i=1,2,...,Popsize,并且 表示种群Pt中的第i个个体,其随机产生公式为:其中维数下标j=1,2,3; 表示了支持向量机的3个优化设计参数的值,即 为支持向量机的惩罚因子C, 为支持向量机的径向基核参数g, 为支持向量机的不敏感损失函数中参数ε;rand(0,1)表示在[0,1]之间产生随机实数的函数,LBj和UBj分别表示支持向量机的第j个优化设计参数的搜索下界和上界;
步骤5,计算种群Pt中每个个体 的适应值 其中个体下标i=1,2,...,Popsize,个体 的适应值 的计算方法为:以个体 作为支持向量机的训练参数,并利用训练数据集对支持向量机进行训练,其中支持向量机的输入变量为归一化的一天的稀土矿区地下水质指标数据:水温、pH值、亚硝酸盐氮、硝酸盐氮、总氮、溶解氧、五日生化需氧量、氨氮;
支持向量机的输出为归一化的KD天后的稀土矿区地下水氨氮浓度值;然后计算训练好的支持向量机在测试数据集上的均方误差MEi,并令个体 的适应值步骤6,令种群Pt中所有个体的未更新次数 其中个体下标i=1,2,...,Popsize;
步骤7,令当前评价次数FEs=FEs+Popsize;
步骤8,保存种群Pt中的最优个体Bestt;
步骤9,雇佣蜂执行其搜索操作;
步骤10,根据种群Pt中个体的适应值计算所有个体的选择概率;
步骤11,观察蜂根据种群Pt中每个个体的选择概率选择出个体执行高斯变异操作生成新个体,然后选择优秀个体进入下一代种群并计算个体的未更新次数,具体步骤如下:步骤11.1,令计数器i=1;
步骤11.2,根据种群Pt中每个个体的选择概率采用轮盘赌策略选择出个体 并令新个体步骤11.3,令邻域下标RNI=(SI-NK+Popsize)%Popsize,其中SI表示轮盘赌策略选择出个体的下标,%表示取余运算符;
步骤11.4,令邻域最优个体 邻域均值个体 并令计数器mt=1;
步骤11.5,令邻域下标RNI=(RNI+1)%Popsize;
步骤11.6,如果个体 比RBestt更优,则令 否则保持RBestt不变;
步骤11.7,令邻域均值个体 然后令计数器mt=mt+1;
步骤11.8,如果mt小于或等于NK×2,则转到步骤11.5,否则转到步骤11.9;
步骤11.9,令邻域均值个体RMeant=RMeant/(NK×2+1),然后在[1,D]之间随机产生一个正整数RDI;
步骤11.10,令随机权值RW=rand(0,1);
步骤11.11,令均值
步骤11.12,令标准差
步骤11.13,以GMean为均值,GSD为标准差产生一个高斯随机实数Val,如果Val的值超出了[LBRDI,UBRDI]之间的范围,则采用同样的方法重新产生高斯随机实数Val,直到Val的值不超出[LBRDI,UBRDI]之间的范围,然后令步骤11.14,计算新个体Ut的适应值Fit(Ut),然后在个体 与新个体Ut之间选择优秀个体进入下一代种群,并计算个体 的未更新次数步骤11.15,令计数器i=i+1;
步骤11.16,如果计数器i小于或等于Popsize,则转到步骤11.2,否则转到步骤12;
步骤12,令当前评价次数FEs=FEs+Popsize×2;
步骤13,侦察蜂找出种群Pt中未更新次数最大的个体,并标记该个体为 如果个体的未更新次数小于Limit,则转到步骤14,否则对个体 执行以最优个体为参考点的反向学习策略,其中以最优个体为参考点的反向学习策略的具体步骤如下:步骤13.1,令计数器j=1;
步骤13.2,令反向个体
步骤13.3,令随机重建个体
步骤13.4,令 其中Bestt为种群Pt中的最优个体;
步骤13.5,令RNBj=LBj+rand(0,1)×(UBj-LBj);
步骤13.6,令j=j+1;
步骤13.7,如果j大于D则转到步骤13.8,否则转到步骤13.4;
步骤13.8,计算个体BRB和个体RNB的适应值,并令当前评价次数FEs=FEs+2;
步骤13.9,令个体EXB为个体BRB和个体RNB两者之间的更优秀者;
步骤13.10,用个体EXB替换种群Pt中的个体 然后转到步骤14;
步骤14,保存种群Pt中最优个体Bestt;
步骤15,令当前演化代数t=t+1;
步骤16,重复步骤9至步骤15直至当前评价次数FEs达到MAX_FEs后结束,将执行过程中得到的最优个体Bestt作为支持向量机的训练参数,并用训练数据集来训练支持向量机,将归一化的一天的稀土矿区溪水指标数据:水温、pH、亚硝酸盐氮、硝酸盐氮、总氮、溶解氧、五日生化需氧量、氨氮输入到训练好的支持向量机,计算支持向量机的输出即可预测稀土矿区地下水KD天后的氨氮浓度值。