1.基于反向精简差分演化算法的水质模型参数识别方法,其特征在于,包括以下步骤:步骤1,根据待研究的水质规律问题的实际情况选择合适的水质模型,并针对选择的水质模型确定待识别的参数,以及待识别参数的个数D;
步骤2,通过现场实验获得水质模型的样本数据,包括水质模型的实际输入数据以及对应的实际输出值;
步骤3,根据所有现场实验采集到的样本数据,针对每一组输入数据,以水质模型的计算输出值与样本数据中对应的实际输出值之间的误差平方和作为反向精简差分演化算法的适应值函数;
步骤4,用户初始化参数,所述初始化参数包括待识别水质模型参数的个数D,种群大小Popsize,反向学习率Po,最大评价次数MAX_FEs;
步骤5,当前演化代数t=0,并设置综合学习率 杂交率 其中下标i=1,...,Popsize,当前评价次数FEs=0;
步骤6,随机 产生初始种 群 其中:下标 i=
1,...,Popsize,并且 为种群Pt中的第i个个体,其随机初始化公式为:其中下标j=1,...,D,并且D表示待识别水质模型参数的个数; 为在种群Pt中的第i个个体,存储了D个待识别水质模型参数的值;rand(0,1)为在[0,1]之间服从均匀分布的随机实数产生函数,Loj和Upj分别为第j个待识别水质模型参数的取值范围的下界和上界;
步骤7,利用所有现场实验采集到的样本数据,计算种群Pt中每个个体的适应值,其中适应值越小表明个体越优秀,计算适应值的方法为:分别将种群Pt中每个个体所表示的水质模型参数代入到选择的水质模型中,并对每一组输入数据,计算出水质模型的理论输出值,将所有水质模型的理论输出值与样本数据中对应的实际输出值之间的误差平方和作为个体的适应值;
步骤8,当前评价次数FEs=FEs+Popsize,并保存种群Pt中适应值最小的个体为最优个体Bestt,设置适应性反向学习种子uOK=0.5;
步骤9,在[0,1]之间产生一个随机实数ra;
步骤10,如果ra小于反向学习率Po,则转到步骤11,否则转到步骤14;
步骤11,执行适应性反向学习策略,对种群Pt产生一个适应性反向种群OPt,并从种群Pt和适应性反向种群OPt中选择出最优的前Popsize个个体进入下一代种群Pt+1,其步骤如下:步骤11.1,设置成功反向学习因子列表SKList为空,并令适应性反向种群OPt为空;
步骤11.2,按以下公式计算种群Pt在每一维上的下界LAj以及上界UBj:其中个体下标i=1,...,Popsize,维数下标j=1,...,D,min为取最小值函数,max为取最大值函数;
步骤11.3,令计数器i=1;
步骤11.4,如果计数器i大于种群大小Popsize,则转到步骤11.15,否则转到步骤
11.5;
步骤11.5,以适应性反向学习种子uOK为期望,0.1为标准差产生一个服从正态分布的随机实数Vk;当Vk的值超出了[0,1]之间的范围,则采用同样的方法重新产生Vk,直到Vk的值不超出[0,1]之间的范围,然后令适应性反向学习因子OK=Vk;
步骤11.6,令个体 的适应性反向个体步骤11.7,令计数器j=1;
步骤11.8,如果计数器j大于D,则转到步骤11.11,否则转到步骤11.9;
步骤11.9,
步骤11.10,令计数器j=j+1,转到步骤11.8;
步骤11.11,计算适应性反向个体 的适应值,并把适应性反向个体 添加到适t应性反向种群OP 中;
步骤11.12,如果适应性反向个体 的适应值小于个体 的适应值,则将OK加入到成功反向学习因子列表SKList中;
步骤11.13,令计数器i=i+1;
步骤11.14,转到步骤11.4;
t t
步骤11.15,从种群P 和适应性反向种群OP 中选择出最优的前Popsize个个体进入t+1下一代种群P ;
步骤11.16,计算成功反向学习因子列表SKList中所有存储数值的平均值MeanOK,并计算适应性反向学习种子uOK=0.9·uOK+0.1·MeanOK;
步骤12,当前评价次数FEs=FEs+Popsize,保存种群Pt中适应值最小的个体为最优t个体Best,当前演化代数t=t+1;
步骤13,转到步骤9;
步骤14,令计数器i=1;
步骤15,如果计数器i大于种群大小Popsize,则转到步骤23,否则转到步骤16;
t t
步骤16,计算个体 的当前杂交率NCri 和当前综合学习率NPri,计算公式如下:其中r1和r2均为在[0,1]之间随机产生的实数;
步骤17,根据个体 的当前综合学习率NPrit,对个体 产生一个综合学习个体 其步骤如下:步骤17.1,令计数器j=1;
步骤17.2,如果计数器j大于D,则转到步骤18,否则转到步骤17.3;
步骤17.3,在[0,1]之间产生一个随机实数r3;如果r3小于个体 的当前综合学习t率NPri 则转到步骤17.7,否则转到步骤17.4;
步骤17.4,在[1,Popsize]之间随机产生两个不相等的正整数RI1,RI2;
步骤17.5,如果个体 的适应值小于个体 的适应值,则 否则步骤17.6,令计数器j=j+1,转到步骤17.2;
步骤17.7, j=j+1,转到步骤17.2;
步骤18,以 为个体 的当前杂交率,采用高斯随机采样方法产生个体 的试验个体 并计算试验个体 的适应值,具体步骤如下:步骤18.1,令计数器j=1;
步骤18.2,在[1,D]之间随机产生一个正整数jRand;
步骤18.3,如果计数器j大于D,则转到步骤18.8,否则转到步骤18.4;
步骤18.4,在[0,1]之间产生一个随机实数r4,如果r4小于个体 的当前杂交率tNCri 或者jRand等于计数器j,则转到步骤18.5,否则转到步骤18.6;
步骤18.5,以 为期望, 为标准差产生一个服从正态分布的随机实数Vg;如果Vg的值超出了[Loj,Upj]之间的范围,则采用同样的方法重新产生Vg,直到Vg的值不超出[Loj,Upj]之间的范围,然后 转到步骤18.7;
步骤18.6,
步骤18.7,令计数器j=j+1,转到步骤18.3;
步骤18.8,计算试验个体 的适应值,转到步骤19;
步骤19,按以下公式在个体 与试验个体 之间选择出个体进入下一代种群:步骤20,按以下公式更新个体 的杂交率 和综合学习率步骤21,令计数器i=i+1;
步骤22,转到步骤15;
步骤23,当前评价次数FEs=FEs+Popsize;保存种群Pt中适应值最小的个体为最优t个体Best,当前演化代数t=t+1;
步骤24,转到步骤9;
步骤25,重复步骤9至步骤24直至当前评价次数FEs达到MAX_FEs后结束,执行过程t中得到的最优个体Best 即为水质模型参数的识别结果。