具体实施方式

下面结合附图和具体实施例，进一步阐明本发明。

一种基于模糊神经网络的水质COD预测方法，如图1和图2所示，包括以下步骤：

步骤1：采集湖泊水样，确定合适的采样点和采样频率，每天固定上午8点在湖泊固定五处采集湖泊水样，使用透明塑料杯存放并各自贴上标签，水样存储时间不超过24小时；

步骤2：先使用便携式温度pH测定仪测量水样的温度和pH，再用多参数水质测定仪依次测量水样的COD、总磷、总氮、氨氮和高锰酸盐指数含量，得到89组数据集D，

D＝{(D_COD1,D_TP1,D_TN1,D_NH3-N1,D_CODMN1,D_temp1,D_pH1),...,(D_CODn,D_TP1n,D_TNn,D_NH3-Nn,D_CODMNn,D_tempn,D_pHn)}

步骤3：对测量得到的几种指标与COD进行相关性排序，选取五个相关系数较高的指标(氨氮、pH、温度、总氮、高锰酸盐指数)作为预测输入，预测集为D_p：

再对相关性较高的五种水质参数进行预处理即平滑处理和归一化处理，最终得到的数据集为D_p ^*：

步骤4：建立合适的模糊神经网络模型，随机将数据集D_p ^*拆分为73组训练集和16组预测集，五种水质参数作为网络的输入数据，COD含量作为网络的输出数据建立模型训练网络，网络反复迭代10⁵次；

步骤5：用训练好的网络模型对随机训练数据进行预测，使用训练集预测训练集中COD含量，比较COD的真实值y和预测值使用均方误差值相对误差值和相关系数比较实际值和预测值之间的误差大小。用模型对训练集进行预测，相对误差控制在6％以内，训练得到的网络对训练数据本身有着非常好的拟合度，并且没有出现过拟合的现象。训练集的网络实际输出与真实值之间的均方误差值为0.1050，相关系数为0.9669。；

步骤6：用训练好的模糊神经网络模型预测其余的预测数据，使用预测集预测预测集中COD含量，一共预测8次，比较COD真实值y和预测值使用均方误差值相对误差值和相关系数比较实际值和预测值之间的误差大小。我们使用模型对预测数据进行8次实验，由8次的实验数据可知，该模型的相关系数基本都在0.8-0.95之间，相对误差几乎都控制在15％以内且没有出现过拟合的现象。

步骤1中确定合适的采样点和采样频率，力求以最低的采样频次，取得最有时间代表性的样品，既要满足能反映水质状况的要求，又要切实可行。

步骤2中的多参数测量需要用到的器材有多参数水质测定仪、智能多参数消解仪、便携式温度pH测定仪，烧杯，30mm比色皿若干，10mm比色皿若干，150mm密封试管若干，100mm密封试管若干，敞口试管若干，大小移液枪两个，配套试剂。实验包括以下步骤：

步骤2.1：测量水质中COD的含量，使用快速消解分光光度法，将水样与配套试剂在消解器中进行快速氧化还原反应，反应后产生三价铬离子，通过分光光度法测定其浓度。首先打开消解仪和水质测定仪，选择COD消解和COD低量程皿模式，准备数支敞口反应管置于冷却架的空冷却槽上。准确量取2.5mL蒸馏水加入到0号反应管中，分别量取各水样2.5mL依次加入到其他反应管中，依次向各个反应管中加入0.7mL重铬酸钾试剂，再依次向反应管中加入4.8mL的硫酸银试剂，混匀后放入消解孔中165摄氏度消解10分钟。消解完成后拿出放在冷却槽上空气冷却2分钟，依次向各反应管中加入2.5mL蒸馏水并混匀，再水冷却2分钟。将水冷却完成后的溶液依次倒入对应编号的30mm比色皿中，0号比色皿为空白溶液，放入比色池中并盖上盖子，按空白键，再依次将其余溶液放入比色池中测量COD含量；

步骤2.2：测量水质中总磷的含量，使用钼酸铵分光光度法，水中的含磷化合物在过硫酸钾的作用下转变为正磷酸盐。正磷酸盐在酸性介质中,可同钼酸铵和酒石酸氧锑钾反应,生成磷钼杂多酸。磷钼酸能被抗坏血酸还原,生成深色的磷钼蓝，在700nm波长下,测定样品的吸光度，再得到样品浓度。首先打开消解仪开关选择总磷消解模式，打开水质测定仪选择总磷低量程皿模式，准备数支150mm密封反应管置于冷却架的空冷却槽上。准确量取8mL蒸馏水和其他水样依次加入反应管中，依次向各个反应管中加入1mL过硫酸钾试剂，将瓶盖拧紧后摇匀水样，将各反应管放入消解仪中120摄氏度消解30分钟。消解完成后拿出放在冷却槽上空气冷却2分钟，再水冷却2分钟。之后依次加入1mL抗坏血酸溶液和1mL钼酸盐溶液试剂，混匀后静置10分钟。将0号溶液倒入30mm比色皿之后放入测定仪中测量，并置空白样，之后再依次测量其余溶液；

步骤2.3：测量水质中总氮的含量，使用变色酸分光光度法，水样在碱性条件下，以过硫酸钾作为氧化剂，经高温高压消解后，所有含氮化合物都转化为硝酸盐氮的形式，剩余的过硫酸钾通过加入过量的焦亚硫酸钠反应消除，硝酸盐-氮在强酸性条件下与铬变酸反应生成黄色化合物，显色液的吸光度值与水样中的总氮含量成正比。首先打开消解仪开关选择总氮消解模式，打开水质测定仪选择总氮模式，准备数支150mm密封反应管置于冷却架的空冷却槽上。准确量取蒸馏水和各水样5mL放入反应管中，依次向各个反应管中加入2mL专用试剂NTA，将瓶盖拧紧后摇匀水样，将各反应管依次放入消解仪中122摄氏度消解40分钟。消解完成后将样品放到冷却架上先空气冷却2分钟，再水冷却2分钟。冷却完成后，依次向各个反应管中加入1mL NTB试剂并摇匀。准备数支100mm的密封比色管置于试管架上，准确量取4mL NTC试剂加入到备好的密封比色管中，分别量取各处理好的水样1mL沿壁管加入到已加好NTC试剂的比色管中，加完水样立即拧紧瓶盖翻转摇匀10次，放入冷却槽中，水冷却10分钟。水冷却结束后，将各反应管从冷却槽中拿出，擦拭干净管壁后放入水质测定仪中测量总氮含量；

步骤2.4：测量水质中氨氮的含量，使用纳氏试剂分光光度法，以游离态的氨或铵离子等形式存在的氨氮与纳氏试剂反应生成淡红棕色络合物，于波长420nm处测量吸光度。首先打开水质测定仪选择氨氮模式。准备数支敞口反应管置于冷却架的空冷却槽上。准确量取10mL无氨水和其他水样加入到各反应管中，向各个反应管中加入1mL N3试剂，1mL N2试剂，摇匀后静置10分钟。然后将溶液依次倒入10mm比色皿中，放入水质测定仪中测量各水样氨氮浓度；

步骤2.5：测量水质中高锰酸盐指数含量，使用分光光度法，在酸性加热条件下用过量高锰酸钾氧化水中还原物质，在波长525nm处测定剩余高锰酸钾的量并得到高锰酸盐指数浓度。首先打开消解仪和水质测定仪，准备数支150mm密封反应管置于冷却架的空冷却槽上。准确量取蒸馏水和各水样5mL放入反应管中，依次向各个反应管加入1mL CM1试剂并混匀，再加入1mL CM2试剂并混匀。将各反应管的瓶口拧紧后依次放入消解仪中105摄氏度消解32分钟，消解完后将各样品放到冷却架上空气冷却2分钟，再水冷却2分钟。依次向各反应管中加入1mL CM3试剂，将其混匀后静置2分钟，再将溶液倒入30mm比色皿放入水质测定仪中测量高锰酸盐指数浓度。

步骤3中参数预处理包括以下方法：

使用皮尔逊相关系数对几种指标和COD之间的相关性进行分析并排序，对预测模型中的输入数据进行相关性排序，可以有效筛选出相关性高、对模型有效的指标，简化预测模型，消除模型中的冗余指标。皮尔逊相关系数是一种可以反应两种变量之间相关联程度大小的统计学方法，是两个变量之间的协方差和标准差之间的商，计算公式为：

最终得到氨氮与COD之间的相关系数的绝对值为0.3734，pH与COD之间的相关系数的绝对值为0.1595，温度与COD之间的相关系数的绝对值为0.1498，总氮与COD之间的相关系数的绝对值为0.1393，高锰酸盐指数与COD之间的相关系数的绝对值为0.1272，总磷与COD之间的相关系数的绝对值为0.0295。说明与COD含量相关性较高的五个指标为氨氮、pH、温度、总氮和高锰酸盐指数。

对数据进行归一化处理，将数据归一化至[0，1]内，数据归一化处理可以大幅度提高神经网络模型的训练速度，保证网络的良好性能，归一化相应方程式和反归一化方程为：

X_i＝(X_max-X_min)·X'_i+X_min

采用lowess平滑方法平滑数据，减小实验误差对神经网络模型的干扰。它是使用加权线性最小二乘法和一级多项式模型进行局部回归的一种平滑方法，其相应的方程式为：

其中平滑值是x_i处的(加权)平均值或(加权)回归预测。

步骤4中建立合适的模糊神经网络模型如图3所示，包括以下步骤：

步骤4.1：根据样本的输入输出维数确定网络的网络结构，根据样本的输入维数为5确定模型的输入节点个数为5，样本的输出维数为1确定模型的输出节点为1，根据输入和输出节点个数人为确定隶属度函数个数为8，即隐含节点为8，网络结构为5-8-1。确定迭代次数为10⁵，对模糊神经网络进行初始化，随机初始化模糊隶属度函数中心c，宽度b和参数p₀～p₆；

步骤4.2：根据模糊规则计算各输入变量x_j的隶属度；

步骤三：将各隶属度进行模糊计算，用来匹配模糊规则；

步骤四：根据模糊计算结果计算模糊模型的输出值；

步骤五：计算误差；

步骤六：修正神经网络系数；

步骤七：修正隶属度函数的中心和宽度。

本发明未作详细阐述的内容属于本专业领域技术人员公知的已有技术。