本发明涉及污水处理领域，且特别涉及一种基于灰色神经网络组合模型的污水处理出水参数预测方法。
背景技术：
污水处理是控制水污染的关键一环。由于相关污水监控设备价格昂贵，因此部分中小污水处理厂出于成本问题或出水监测设备维护困难等问题，存在相关参数监控不理想情况。然而一些参数的不精确不仅对出水水质有重要影响，同时也增加了污水处理成本。
目前，对污水处理厂水质预测的方法有很多，但是大都是基于神经网络的预测。污水处理厂水质的变化是有多种影响因素构成的，单一的预测方法仅能对局部的变化规律进行描述，因而预测精度不高。因此，改进单一模型的预测方法便成了提高预测精度的重要手段。
技术实现要素：
本发明为了克服现有单一性预测方法预测准确性低的问题，提供一种预测精度高的基于灰色神经网络组合模型的污水处理出水参数预测方法。
为了实现上述目的，本发明提供一种基于灰色神经网络组合模型的污水处理出水参数预测方法，该方法包括：
获取历史的进水参数数列x(t)并进行累加处理后形成累加数列y(t)；
根据累加数列y(t)建立灰色神经网络并利用灰色神经网络对历史出水参数进行预测；在预测的过程中采用改进的粒子群算法进行灰色神经网络的优化，优化的步骤如下：
步骤一，初始化改进的粒子群算法内的参数；
步骤二，计算每个离子的粒子的适应度f：
其中，yij为测试样本总群x的预测值；tij为测试样本总群x的真实值；n为测试样本的个数；m为输出节点的个数；
步骤三，根据计算的适应度f，找出最小适应度fmin对应的粒子个体xmin，令x*＝xmin，f*＝fmin，则x*为最优粒子，f*为最优适应度；
步骤四，按照以下公式优化每个粒子的位置xi和速度vi；
样本种群x＝(x1，x2，…，xn)；其中粒子i的位置信息表示为xi＝(xi1，xi2，…，xid)t；速度信息表示为vi＝(vi1，vi2，…，vid)t；粒子i的个体极值点为pi＝(pi1，pi2，…，pid)t；种群的全局极值点为g＝(g1，g2，…，gd)t；
步骤五，以变异遗传概率p对粒子xi的位置进行更新，产生新的粒子xi′并线性递减惯性权重，计算新的粒子xi′的适应度f′；比较f和f′的大小，若f＞f′，则令xi＝xi′；
步骤六，对于更新后的种群，找出最小的fmin，及其对应的xmin，若fmin＜f*，则令x*＝xmin，否则不进行变换；
步骤七，重复步骤四至步骤六，当迭代次数达到最大值时结束循环，输出f*，x*。把x*的值赋给参数a，bi，实现灰色神经网络的优化；
根据优化后的神经网络进行当前出水参数的预测；
采用马尔科夫链对计算所得的当前的出水参数进行修正。
于本发明一实施例中，根据累加数列建立灰色神经网络并利用灰色神经网络对历史出水参数进行预测的步骤如下：
根据累加数列y(t)构建灰色神经网络的微分方程：
求解微分方程后得到其时间响应函数：
其中：
将时间响应函数内的参数映射到bp神经网络内相应层的权重上，得到灰色神经网络；
在累加数列y(t)获取测试样本对灰色神经网络进行训练，当网络收敛时提取a，b1，b2，…，bn-1这些时间响应函数内的参数，将时间响应函数进行白化；
利用白化后的时间响应函数对历史出水参数进行预测。
于本发明一实施例中，bp神经网络内相应层的权重的赋值如下：
输入变量的权重w11＝a，输入层与隐含层的权重
y1＝w31c1+w32c2+…+w3ncn-θy1(5)
其中：c1＝hw21，ci＝yi(t)hw2i，i＝2,…,n，
于本发明一实施例中，采用以下公式以变异遗传概率p对粒子xi的位置进行更新：
其中，xij表示第i个粒子的第j个方向的位置信息；f(k)＝r(1-k/kmax)，k表示当前迭代次数，kmax表示最大迭代次数；amax，amin分别为xij的上界和下界，r为[0,1]之间的随机数；
采用以下公式线性递减惯性权重：
ω(k)＝ωs-(ωs-ωe)(kmax-k)/kmax(7)
其中，ωs为初始惯性权重；ωe为最终迭代惯性权重，ω(k)为第k次迭代时的惯性权重。
于本发明一实施例中，采用马尔科夫链对计算所得的当前的出水参数进行修正的步骤如下：
采用预测结果相对误差的绝对值，将预测值序列分为若干个状态区间，求得一步转移概率矩阵
pt+1＝p0[p(1)]t+1(8)
pt+1为t+1时刻预测值的概率分布；p0为初始时刻预测值的概率分布。
综上所述，本发明提供的基于灰色神经网络组合模型的污水处理出水参数预测方法利用改进的粒子群算法与灰色神经网络相融合，用改进粒子群算法优化灰色神经网络权值，利用灰色神经网络的输出的均方差作为粒子适应度，实现了模型参数的优化。再用马尔科夫链对误差进行修正，使得修正值更加接近实际值。多种算法的融合，综合利用各自算法的优点，提高了污水处理水质参数的预测精度，具有良好的扩展性和实际应用价值。
为让本发明的上述和其它目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合附图，作详细说明如下。
附图说明
图1所示为本发明一实施例提供的基于灰色神经网络组合模型的污水处理出水参数预测方法的流程图。
具体实施方式
如图1所示，本实施例提供的基于灰色神经网络组合模型的污水处理出水参数预测方法包括：
步骤s1、获取历史的进水参数数列x(t)并进行累加处理后形成累加数列y(t)。所述进水参数包括生化需氧量(bod)、化学需氧量(cod)以及ph值等参数。
步骤s2、根据累加数列y(t)建立灰色神经网络并利用灰色神经网络对历史出水参数进行预测；在预测的过程中采用改进的粒子群算法进行灰色神经网络的优化。
在该步骤中首先将灰色系统和神经网络进行融合，利用神经网络模型的大样本、非线性处理能力和学习能力来弥补灰色模型贫信息、建模简单及非线性处理能力弱等特征。建立灰色神经网络并利用灰色神经网络对历史出水参数进行预测的步骤如下：
首先，根据累加数列y(t)构建灰色神经网络的微分方程：
接着，求解微分方程后得到其时间响应函数：
其中：
再接着，将时间响应函数内的参数a，bi映射到bp神经网络内相应层的权重上，得到灰色神经网络。在该步骤中，权重的赋予如下：
输入变量的权重w11＝a，输入层与隐含层的权重
y1＝w31c1+w32c2+…+w3ncn-θy1(5)
其中：c1＝hw21，ci＝yi(t)hw2i，i＝2,…,n，
再接着，在累加数列y(t)获取测试样本对灰色神经网络进行训练，当网络收敛时提取a，bi这些时间响应函数内的参数，将时间响应函数进行白化；
最后，利用白化后的时间响应函数对历史出水参数进行预测。
在对历史出水参数进行预测的过程中，采用改进的粒子群算法进行灰色神经网络的优化，优化灰色神经网络的权重，提高预测的准确性。具体的优化的步骤如下：
步骤一，初始化改进的粒子群算法内的参数。改进的粒子群算法的参数包括：加速系数c1和c2；最大迭代次数kmax，种群规模size，最优位置xmax，最大速度vmax，变异遗传概率p等参数进。再由参数a，bi的个数确定种群长度l。
步骤二，采用公式(1)计算每个离子的粒子的适应度f：
其中，yij为测试样本总群x的预测值；tij为测试样本总群x的真实值；n为测试样本的个数；m为输出节点的个数；
步骤三，根据计算的适应度f，找出最小适应度fmin对应的粒子个体xmin，令x*＝xmin，f*＝fmin，则x*为最优粒子，f*为最优适应度；
步骤四，按照公式(2)和(3)优化每个粒子的位置xi和速度vi；
假设在一个d维空间中，种群x＝(x1，x2，…，xn)，其中粒子i的位置信息表示为xi＝(xi1，xi2，…，xid)t，速度信息表示为vi＝(vi1，vi2，…，vid)t，粒子i的个体极值点为pi＝(pi1，pi2，…，pid)t，种群的全局极值点为g＝(g1，g2，…，gd)t。
步骤五，按照公式(6)以变异遗传概率p对粒子xi的位置进行更新，产生新的粒子xi′并按照公式(7)进行线性递减惯性权重，计算其粒子xi′的适应度f′；比较f和f′的大小，若f＞f′，则令xi＝xi′。
其中，xij表示第i个粒子的第j个方向的位置信息；f(k)＝r(1-k/kmax)，k表示当前迭代次数，kmax表示最大迭代次数；amax，amin分别为xij的上界和下界，r为[0,1]之间的随机数。
ω(k)＝ωs-(ωs-ωe)(kmax-k)/kmax(7)
其中，ωs为初始惯性权重；ωe为最终迭代惯性权重。
步骤六，对于更新后的种群，找出最小适应度fmin，及其对应的粒子个体xmin，若fmin＜f*，则令x*＝xmin，否则不进行变换；
步骤七，重复步骤四至步骤六，当迭代次数达到最大值时结束循环，输出f*，x*。把x*的值赋给参数a，bi，参数a，bi与灰色神经网络内输入层和隐含层的权重相关，因此，实现灰色神经网络的优化。
步骤s3、根据优化后的神经网络进行当前出水参数的预测；
步骤s4、采用马尔科夫链对计算所得的当前的出水参数进行修正。修正的方法为：采用预测结果相对误差的绝对值，将预测值序列分为若干个状态区间，求得一步转移概率矩阵
pt+1＝p0[p(1)]t+1(8)
pt+1为t+1时刻预测值的概率分布；p0为初始时刻预测值的概率分布。
综上所述，本发明提供的基于灰色神经网络组合模型的污水处理出水参数预测方法利用改进的粒子群算法与灰色神经网络相融合，用改进粒子群算法优化灰色神经网络权值，利用灰色神经网络的输出的均方差作为粒子适应度，实现了模型参数的优化。再用马尔科夫链对误差进行修正，使得修正值更加接近实际值。多种算法的融合，综合利用各自算法的优点，提高了污水处理水质参数的预测精度，具有良好的扩展性和实际应用价值。
虽然本发明已由较佳实施例揭露如上，然而并非用以限定本发明，任何熟知此技艺者，在不脱离本发明的精神和范围内，可作些许的更动与润饰，因此本发明的保护范围当视权利要求书所要求保护的范围为准。