1.一种水质综合生物毒性远程自动分析仪,其特征在于,包括水质综合生物毒性分析装置、数据发送装置以及数据接收和显示装置,其中,水质综合生物毒性分析装置用于监测采样液中的有毒物质对生物发光杆菌发光强度的抑制程度;所述数据发送装置用于将发光强度信息以无线的方式发送给设置在远端的数据接收和显示装置,所述数据接收和显示装置对发光强度信息进行数据处理后确定污水毒性的污染等级,并将结果信息显示出来;
所述水质综合生物毒性分析装置包括仪器壳体、发光菌的培植单元、发光菌发光强度采集和存储单元、发光菌发光强度传输单元,数据接受单元和数据显示单元;
所述数据接收和显示装置对发光强度信息进行数据处理后确定污水毒性的污染等级包括以下步骤:步骤1:通过基于数值相近原则的曲线拟合算法对发光强度信息进行数据处理,即利用最小二乘解法得到曲线拟合的参数方程:X=[ATA]-1ATB,得到模型y1,并求出精确值,其中X为最小二乘法参数,A、B为拟合系数;
步骤2:通过基于形态相似原则的曲线拟合算法对发光强度信息进行数据处理:(1)设时间影响因子为λ,反应发光菌活性受存活时间的量;设形态相似拟合曲线为:y2(x)=a0+a1x+a2x2+...+anxn,其中ai为各阶系数;
(2)将拟合曲线结合形系数方程 其中x12i为拟合点值,S12为偏量值,可得:(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:令xk+1=yn+βnpn, 同时,令
{p0,p1,....pn-1}k+1={p0,...,pm-1,pm+1,....pn-1,pn}k=k+1转步骤b;
步骤f:令xk+1=yn,置k=k+1转步骤b;
(5)通过powell算法,求出a1,a2,......,an,考虑到a0为直流分量,不会影响拟合曲线的
2 n
形态,代入式y2(x)=a0+a1x+a2x+...+anx中,采用最小二乘法求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次方系数am的初始区间[A1,B1]为
2)先判断ai在初始区间[A1,B1]是单峰函数,求出的最佳aT就是ai在区间[A1,B1]中的近似极小值或极大值,即ai在[A1,aT]区间应严格递减或递增,在[aT,B1]上应严格递增或递减,采用Fabonacci法对其进行优化;
3)根据最佳拟合曲线的最大正、负偏差的绝对值近似相等,则偏差取两者之和的一半,即由此可得最佳常系数为
步骤5:利用得到的改进拟合曲线y结合本系统数据库中预存的专家系统,进行曲线匹配,分别得到改进拟合曲线的拟合系数B与专家系统中发光菌与不同种类、浓度的毒性物质各反应机理曲线的拟合系数A,求取拟合系数B与各拟合系数A之间的距离,系数最小的将为所预测的毒性物质。