仪器分析中的计算机方法

仪器分析中的计算机方法---回归分析的原理及应用在分析化学，特别是仪器分析中，常常需要做工作曲线（也叫标准曲线，或校正曲线，或检量线）。例如,原子吸收法中作吸光度和浓度的工作曲线，极谱法中作波高和浓度的工作曲线等等。在分析化学中所使用的工作曲线，通常都是直线。一般是把实验点描在坐标纸上，横坐标X表示被测物质的浓度，叫自变量。大都是把可以精确测量或严格控制的变量（如标准溶液的浓度）作为自变量；纵坐标y表示某种特征性质（如吸光度、波高等）的量，称因变量，一般设因变量是一组相互独立、其误差服从同一正态分布N（Ο，σ2）的随机变量。然后根据坐标纸上的这些散点（实验点）的走向，用直尺描出一条直线。这就是分析工作者习惯的制作工作曲线的方法。若吸光度----浓度的直线能通过所有实验点，在统计上就说溶液的吸光度和浓度有最密切的线性关系。吸光度完全依赖于浓度的改变而变，完全遵循比尔定律。实验条件中的各种偶然因素对它无任何影响（亦即没有实验误差）。我们称这种关系为确定性关系或函数关系。这时做工作曲线图的任务比较简单，借助于一支直尺和一支铅笔，就能完成。但是由于实验中不可避免的有误差存在，实验点全部密集在回归线上的情况通常是极少见的，尤其当误差较大时，实验点比较分散，并不在一条线上，这时作图就有困难了。因为凭直觉很难判断怎样才能使所联的线对干所有实验点来说是误差最小的，亦即难于确定到底哪条线才是最好的回归线。例如，用火焰原子吸收法测定镁，得到下表数据Mg(ppm)0.000.200.400.600.801.00A0.000.2020.4100.5530.6410.736图中的工作曲线是用Excel的方法回归得到，选取的数据点不同，R就不一样。原子吸收测镁y=0.7343x+0.0565R2=0.9669y=0.9335x+0.0112R2=0.993600.10.20.30.40.50.60.70.80.9100.20.40.60.811.2ppmA一．最小二乘法原理若用（χi，yi）表示n个数据点（i=1，2，3，...,n）,而任意一条直线方程可写成:bxay*(1)在（1）式中,采用y*符号,表示这是一条任意的直线,如果用这条直线来代表x和y的关系,即对每个已知的数据点(xi,yi)来说,其误差为iiibxayyy*(2)令各数据点误差的平方的加和(差方和)为Q,则Q是总的误差：(3)回归直线就是在所有直线中,差方和Q最小的一条直线.换句话说,回归直线的系数b及常数项a,应使Q达到极小值.根据微积分求值的原理,要使Q达到极小值,只需将(3)式分别对a,b求偏微商,令它们等于0.于是a,b满足(4)(5)从(4)式可得到21*)(niiyyQ21)(iniibxayabxaybxayaQiiniii)()(21niiibxay10)(2niiiiibbxaybxaybQ1)()(20)(21iniiixbxaynininiiiiixbnaybxay1110)(niniiixbyna11（6）yx,分别代表xi和yi的平均值。从式（5）可得到将（6）式代入，得所以根据差方和关系式，若令xbyaxnbynaniniii11110)(11112nininiiiiiniiiixbxayxxbxay0)(1211niiniiiniiiixbxnxbnyyx2)(1))((12iiiiiixnxbyxnyx22)(1))((1iiiiiixnxyxnyxb22xnxyxnyxiii)2()(222xxxxxxliiixx222xnxxxii222nxnxnxxiiii222)()(22nxnnxxiii22221xnxxnxiii同理，可推出：由观测值（一组样本）算出a，b的值，称为参数a，b的估算值，用符号baˆ,ˆ,表示，于是回归直线方程式便可确定如下：xbayˆˆˆ式中bayˆ,ˆ,ˆ分别表示由样本求得的y，a，b的估算值。如果xx，则有，yyˆ这种方法就称为最小二乘法，即也就是“最小差方和法”。二．回归方程的类型这里的“线性”，是对a，b而言，对y，x并不一定。只要通过适当变化，a，b仅为一次待确定参数，就可使用这种方法求出。例：1.双曲线xbay11(令xxyy1;1)2.抛物线acxby2)((2)(cxx)3.幂函数bdxy]log,log;log[daxxyy4.指数函数dayydeybxln;ln........[]2)(yyliyy2222)(1ynyynyiii))((yyxxliixyyxnyxiilxxlxyxxyyxxxnxyxnyxbiiiiiiiii222)())(()(1))((15.指数函数xbdey/[daxxyyln;1;ln]6.S型曲线];1.......[1xeexyyxbay7.对数曲线]log......[logxxxbay三．回归方程的显著性检验1．关系数R在求回归方程时，假定y与x存在线性关系。怎样判别这种关系的好坏呢？引入R这个相关系数的概念。首先让我们讨论一些有关概念：回归平方和22)]([)ˆ(xxbyyU剩余平方和2)ˆ(yyQi总离差平方和UQLyy令：LxxLyyLRxyLyyULxxLyyLxyR22R的正负号由Lxy的符号决定，即与b同号。R的绝对值为小于1，大于0的无量纲统计量。当|R|≌1时，表明y与x之间线性关系密切。|R|≌0时，表明y与x之间无线性关系。通常使用R2，具有更实际的意义。2．著性检验F)2/()1/(21nfQfUF)2(nLxybLyyLxybf1-回归差和自由度，f2-残余差方和自由度。F＜Fa（临界F值，见表），y与x无线性关系；F＞Fa，表明回归方程是显著性的，假设是可靠的。3．回归线的精度可以使用回归方程得到y的平均值yˆ。那么实际的y离yˆ值偏差多大呢？即回归的精度如何呢？通常规定，剩余平方和Q除以它的fQ，所得商称为剩余方差：22nQS剩余方差的平方根称为剩余标准偏差：2nQS又可得2)ˆ(2nyySi代入R后，2)1(2nLyyRSS值越小，说明精度越高。四．使用Excel回归计算