发明名称 一种水环境重金属检测系统的温度补偿方法
摘要 本发明公开了一种水环境重金属检测系统的温度补偿方法,将采集的极化电流和极化电位通过差分滤波进行平滑处理后,建立拟合回归方程逆模型,通过求偏差、偏差平方以及偏导,获取温度补偿模型,并结合相对误差度量因子REF进行判断,在不符合标准时,自动增加拟合次数重新拟合,以构建出精确度更高的高次多元回归逆模型,该温度补偿模型应用于水环境重金属检测,能提高其检测精度。
申请公布号 CN105445344A 申请公布日期 2016.03.30
申请号 CN201511021710.3 申请日期 2015.12.30
申请人 桂林电子科技大学 发明人 黄廷磊;林科;张纪元;刘杰;覃志松;张姿;卢艳梅;甘国宁
分类号 G01N27/26(2006.01)I 主分类号 G01N27/26(2006.01)I
代理机构 桂林市持衡专利商标事务所有限公司 45107 代理人 林培
主权项 一种水环境重金属检测系统的温度补偿方法,其特征在于,包括以下步骤:1)采集重金属浓度为C<sub>j</sub>的待测水环境的极化电流I(j)、极化电位E(j)和环境温度U<sub>Tj</sub>;j为自然数;C<sub>j</sub>为采集的第j个重金属浓度,I(j)为采集的第j个极化电流,E(j)为采集的第j个极化电位;2)对极化电流I(j)和极化电位E(j)进行差分法滤波,获取平滑的E/I曲线,进而获得电化学传感器输出峰值电流I<sub>Cj</sub>:3)建立温度补偿模型:31)采集n个浓度标定点和m个温度标定点的极化电流和极化电位值,第j个标定点浓度计算值为C(I<sub>Cj</sub>,U<sub>Tj</sub>),j、m、n为正整数,拟合回归方程逆模型为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>C</mi><mrow><mo>(</mo><mrow><msub><mi>I</mi><mrow><mi>C</mi><mi>j</mi></mrow></msub><mo>,</mo><msub><mi>U</mi><mrow><mi>T</mi><mi>j</mi></mrow></msub></mrow><mo>)</mo></mrow><mo>=</mo><mrow><mo>(</mo><mrow><msub><mi>d</mi><mn>0</mn></msub><mo>+</mo><msub><mi>d</mi><mn>1</mn></msub><msub><mi>I</mi><mrow><mi>C</mi><mi>j</mi></mrow></msub><mo>+</mo><msub><mi>d</mi><mn>2</mn></msub><msub><mi>U</mi><mrow><mi>T</mi><mi>j</mi></mrow></msub></mrow><mo>)</mo></mrow><mo>+</mo><mrow><mo>(</mo><mrow><msub><mi>d</mi><mn>3</mn></msub><msubsup><mi>I</mi><mrow><mi>C</mi><mi>j</mi></mrow><mn>2</mn></msubsup><mo>+</mo><msub><mi>d</mi><mn>4</mn></msub><msub><mi>I</mi><mrow><mi>C</mi><mi>j</mi></mrow></msub><msub><mi>U</mi><mrow><mi>T</mi><mi>j</mi></mrow></msub><mo>+</mo><msub><mi>d</mi><mn>5</mn></msub><msubsup><mi>U</mi><mrow><mi>T</mi><mi>j</mi></mrow><mn>2</mn></msubsup></mrow><mo>)</mo></mrow><mo>+</mo><mn>...</mn><mo>+</mo><mi>&xi;</mi><mo>,</mo></mrow>]]></math><img file="FDA0000896853480000011.GIF" wi="1222" he="83" /></maths>式中,各阶表达式的系数d<sub>k</sub>为温度补偿常系数,其中k为正整数,且0≤k≤t‑1,方程的拟合项数t为正整数,且t≥3;ξ为高阶无穷小量;拟合项数t与拟合次数r的关系为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>t</mi><mo>=</mo><mfrac><mrow><mo>(</mo><mi>r</mi><mo>+</mo><mn>3</mn><mo>-</mo><mn>1</mn><mo>)</mo><mo>!</mo></mrow><mrow><mi>r</mi><mo>!</mo><mrow><mo>(</mo><mn>3</mn><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>!</mo></mrow></mfrac><mo>=</mo><mfrac><mrow><mo>(</mo><mi>r</mi><mo>+</mo><mn>2</mn><mo>)</mo><mo>(</mo><mi>r</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mn>2</mn></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000896853480000012.GIF" wi="518" he="143" /></maths>32)重金属浓度标定值C<sub>j</sub>相对C(I<sub>Cj</sub>,U<sub>Tj</sub>)的偏差为δ<sub>j</sub>:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msubsup><mi>&delta;</mi><mi>j</mi><mn>2</mn></msubsup><mo>=</mo><msup><mrow><mo>&lsqb;</mo><mrow><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>0</mn></mrow><mi>t</mi></munderover><mrow><msub><mi>d</mi><mi>k</mi></msub><msub><mi>p</mi><mrow><mi>j</mi><mi>k</mi></mrow></msub></mrow><mo>-</mo><msub><mi>C</mi><mi>j</mi></msub></mrow><mo>&rsqb;</mo></mrow><mn>2</mn></msup></mrow>]]></math><img file="FDA0000896853480000013.GIF" wi="462" he="135" /></maths>式中,p<sub>jk</sub>是回归方程逆模型中温度补偿常系数d<sub>k</sub>的各项系数,标定点偏差平方和L<sub>S</sub>为:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>L</mi><mi>S</mi></msub><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>S</mi><mo>=</mo><mi>m</mi><mo>&times;</mo><mi>n</mi></mrow></munderover><msup><mrow><mo>&lsqb;</mo><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>0</mn></mrow><mi>t</mi></munderover><msub><mi>d</mi><mi>k</mi></msub><msub><mi>p</mi><mrow><mi>j</mi><mi>k</mi></mrow></msub><mo>-</mo><msub><mi>C</mi><mi>j</mi></msub><mo>&rsqb;</mo></mrow><mn>2</mn></msup><mo>=</mo><mi>L</mi><mrow><mo>(</mo><msub><mi>d</mi><mn>0</mn></msub><mo>,</mo><mo>...</mo><mo>,</mo><msub><mi>d</mi><mrow><mi>j</mi><mo>-</mo><mi>k</mi></mrow></msub><mo>,</mo><mo>...</mo><mo>,</mo><msub><mi>d</mi><mi>t</mi></msub><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000896853480000014.GIF" wi="1252" he="150" /></maths>式中,S为标定点的总个数,L<sub>S</sub>是温度补偿常系数d<sub>k</sub>系数的t元函数;33)分别对L<sub>S</sub>求关于d<sub>k</sub>的偏导,令<img file="FDA0000896853480000015.GIF" wi="197" he="135" />k=0,1,...,t,则温度补偿模型为:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>S</mi></munderover><mo>&lsqb;</mo><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>0</mn></mrow><mi>t</mi></munderover><msub><mi>d</mi><mi>k</mi></msub><msub><mi>p</mi><mrow><mi>j</mi><mi>k</mi></mrow></msub><mo>&rsqb;</mo><mo>&times;</mo><msub><mi>p</mi><mrow><mi>j</mi><mi>k</mi></mrow></msub><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>t</mi></munderover><msub><mi>C</mi><mi>j</mi></msub><mo>&times;</mo><msub><mi>p</mi><mrow><mi>j</mi><mi>k</mi></mrow></msub></mrow>]]></math><img file="FDA0000896853480000016.GIF" wi="605" he="147" /></maths>式中,j=0,1,…,S,S=m*n;34)用于衡量性能的相对误差度量因子REF为:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mi>R</mi><mi>E</mi><mi>F</mi><mo>=</mo><mn>100</mn><mo>&times;</mo><mfrac><msup><mrow><mo>&lsqb;</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mrow><mo>(</mo><mfrac><mrow><msub><mi>C</mi><mi>j</mi></msub><mo>-</mo><msubsup><mi>C</mi><mi>j</mi><mi>S</mi></msubsup></mrow><msubsup><mi>C</mi><mi>j</mi><mi>S</mi></msubsup></mfrac><mo>)</mo></mrow><mo>&rsqb;</mo></mrow><mfrac><mn>1</mn><mn>2</mn></mfrac></msup><msqrt><mi>s</mi></msqrt></mfrac></mrow>]]></math><img file="FDA0000896853480000021.GIF" wi="486" he="230" /></maths>式中,C<sub>j</sub>为第j次检测到的重金属离子浓度值,C<sub>j</sub><sup>S</sup>为第j次检测时的标准溶液中重金属离子浓度的实际值,s为检测次数;当REF大于预设的最大允许误差时,拟合次数r自动增加1,执行步骤31);否则,本次建立温度补偿模型结束。
地址 541004 广西壮族自治区桂林市七星区金鸡路1号