[发明专利]一种基于曲线拟合和毒理分析算法的水质远程在线监测方法

摘要

本发明提供一种基于曲线拟合和毒理分析算法的水质远程在线监测方法，通过水质综合生物毒性分析仪监测采样液中的有毒物质对生物发光杆菌发光强度的抑制程度，并将发光强度信息发送给远程服务器，在服务器中对发光强度信息进行数据处理后确定污水毒性的污染等级，并可以通过远程客户端查询结果信息。相对于现有技术，本发明通过植入改进的曲线拟合算法，以Fabonacci法进一步优化改进的曲线拟合模型的各项系数，实现毒性物质成分和浓度的预测；同时本发明还考虑磁力搅拌时间、培养温度、培养时间、暴露时间及pH值范围对结果测定值的作用下，分析各条件对发光菌生物毒性实验的影响，从而提高分析仪测量精度、稳定性的同时实现远程在线监测。

主权项

1.一种基于曲线拟合和毒理分析算法的水质远程在线监测方法，其特征在于，通过水质综合生物毒性分析仪监测采样液中的有毒物质对生物发光杆菌发光强度的抑制程度，并将发光强度信息发送给远程服务器，在服务器中对发光强度信息进行数据处理后确定污水毒性的污染等级，并可以通过远程客户端查询结果信息；其中，在服务器中对发光强度信息进行数据处理的步骤包括：步骤1：通过基于数值相近原则的曲线拟合算法对发光强度信息进行数据处理，即利用最小二乘解法得到曲线拟合的参数方程：X＝[ATA]‑1ATB，得到模型y1，并求出精确值，其中X为最小二乘法参数，A、B为拟合系数；步骤2：通过基于形态相似原则的曲线拟合算法对发光强度信息进行数据处理：(1)设时间影响因子为λ，反应发光菌活性受存活时间的量；设形态相似拟合曲线为：y2(x)＝a0+a1x+a2x2+...+anxn，其中ai为各阶系数；(2)将拟合曲线结合形系数方程其中x_12i为拟合点值，S₁₂为偏量值，可得：(3)构造出带参数的增广目标函数，当满足约束条件时，等号两边成立；当不满足约束条件范围时，取一个充分大的数μ>0，构造如下的函数：式中a＝{a1,a2,...,an,λ}；(4)通过powell算法的求解函数极值：步骤a：选定初始点x(0)，n个线性无关的向量组，组成初搜索方向组{p0,p1,....pn‑1}，给定精度ε＞0，置k＝0，k为初始极值；步骤b：令y0＝xk，依次沿{p0,p1,....pn‑1}中的方向进行一堆搜索，对应得到辅助迭代点y1,y2,....yn，即式中βj‑1为沿pj‑1方向的步长；如|y(k)‑x(k)|<ε成立，则停止计算，否则执行步骤c；步骤c：构造加速方向，令pn＝yn‑y0，若||pn||≤ε，则停止迭代，输出xk+1＝yn，否则转步骤d；步骤d：确定调整方向：求出m，其中，m为极小值点，使得f(ym‑1)‑f(ym)＝max{|f(ym‑1)‑f(ym)|1≤j≤n}若下式成立：f(y0)‑2f(yn)+f(2yn‑y0)＜2[f(ym‑1)‑f(ym)]，转到步骤f，否则转步骤e；步骤e：令x^k+1＝yⁿ+β_npⁿ，同时，令{p⁰,p¹,....p^n‑1}_k+1＝{p⁰,...,p^m‑1,p^m+1,....p^n‑1,pⁿ}k＝k+1转步骤b；步骤f：令xk+1＝yn，置k＝k+1转步骤b；(5)通过powell算法，求出a1,a2,......,an，考虑到a0为直流分量，不会影响拟合曲线的形态，代入式y2(x)＝a0+a1x+a2x2+...+anxn中，采用最小二乘法求a0，得到基于形态相似准则的曲线拟合算法y2；利用相关系数公式求出精确度：步骤3：从以上的两个步骤，分别得出曲线拟合模型y1，y2，构造改进的曲线拟合模型为：y(x)＝w1y1(x)+w2y2(x)，其中0＜w1≤1，0＜w2≤1，w1、w2依据求出的精确度值；步骤4：采用Fabonacci法进一步的优化已求出的各项系数，使其更加接近最优值；具体步骤如下：1)设第j个离散点与拟合曲线上对应值的偏差为根据上式可以求出，n个离散点中的最大正偏差点和最大负偏差点则第m次方系数a_m的初始区间[A₁,B₁]为2)先判断ai在初始区间[A1,B1]是单峰函数，求出的最佳aT就是ai在区间[A1,B1]中的近似极小值或极大值，即ai在[A1,aT]区间应严格递减或递增，在[aT,B1]上应严格递增或递减，采用Fabonacci法对其进行优化；3)根据最佳拟合曲线的最大正、负偏差的绝对值近似相等，则偏差取两者之和的一半，即由此可得最佳常系数为步骤5：利用得到的改进拟合曲线y结合本系统数据库中预存的专家系统，进行曲线匹配，分别得到改进拟合曲线的拟合系数B与专家系统中发光菌与不同种类、浓度的毒性物质各反应机理曲线的拟合系数A，求取拟合系数B与各拟合系数A之间的距离，系数最小的将为所预测的毒性物质。