最小二乘法在工业纯碱中铁含量测定中的应用

最小二乘法在碳酸钠铁含量测定中的应用摘要：碳酸钠是钻井液中常用的添加剂之一，主要用来调节PH值和对膨润土进行增效。铁含量是评价碳酸钠质量的一项重要指标，其含量大小直接影响碳酸钠的使用效果。测试碳酸钠铁含量的主要方法是邻菲罗啉分光光度法，样品溶液的吸光度随铁含量的不同呈线性变化。在曲线绘制的过程中，人为因素对实验数据的影响较大。运用最小二乘法，对所测得的数据统计分析、处理后进行线性加归，可以较好地提高测试结果的准确性与可靠性。关键词最小二乘法吸光度铁含量线性回归一概述碳酸钠中铁含量的测定方法1、概述铁含量是工业碳酸钠的一个重要质量指标,其含量高低直接决定着纯碱质量及其使用性能。邻菲罗啉分光光度法是测定碳酸钠中微量铁含量的经典方法,样品溶液的吸光度随铁离子浓度的增大而呈上升趋势。平时我们在实验过程中会遇到两量x、y如果存在y=ax+b的线性关系时，其中a、b为线性函数的参数。当实验数据存在这种线性关系时，通常我们运用作图法对其参数进行处理运算、进而求出实验结果。但是作图法很难得到好的结果，而运用最小二乘法可以得到比较好的线性拟合。对其两种方法比较可以最小二乘法的数据处理方法是比较理想的办法。本文运用作图法和最小二乘法对工业纯碱中铁含量的测定的数据进行处理，得出运用最小二乘法可以得到误差比较小的铁含量。在曲线绘制的过程中，人为因素对实验数据的影响较大。运用最小二乘法，对所测得的数据统计分析、处理后进行线性加归，可以较好地提高测试结果的准确性与可靠性。2试验称取10g纯碱样,精确到0.01g,置于烧杯中,加适量水润湿,滴加35mL盐酸溶液(1+1)和12.5mL抗坏血酸,煮沸3~5min,冷却(必要时过滤),移入250mL容量瓶中,加水至刻度,摇匀,用移液管移取50mL该试液于100mL烧杯中;另取7mL盐酸溶液(1+1)和2.5mL抗坏血酸于另一烧杯中,用氨水(2+3)中和后,与试验溶液一并用氨水(1+8)和盐酸溶液(1+3)调节pH值为2(用精密pH试纸检验)。分别移入100mL容量瓶中,各加入10mL乙乙酸铵缓冲溶液和5mL邻菲口罗啉,用水稀释至刻度,摇匀,静置10min。在波长510nm处,用3cm比色皿,以空白溶液为参比,进行比色。3、数据处理根据实验数据列以下表格序号1234567铁标准溶液体积（mL）0.01.02.04.06.08.010.0含铁质量（mg）0.000.010.020.040.060.080.10吸光度0.0000.0700.1260.2500.3750.5360.630根据列表绘制工作曲线：0.000.020.040.060.080.100.00.10.20.30.40.50.60.7absorbancethecontentofFe(mg)要使绘制的曲线尽可能圆滑并且尽可能经过所有的点，但是根据所测数据据绘出的曲线并不呈线性回归，所以只能主观的绘制一条平滑的曲线，根据直线经过的点，代入公式y=ax+b中求得a和b。根据点（0、0）和点（0.01、0.070）我们求出线性方程为y=7x。y代表吸光度x代表铁的含量，在实验中我们测得试验溶液的吸光度为0.082、0.090分别求得x为0.012、0.013。这里出现了以下问题1、空白颜色过深;2、绘制标准曲线时没有线性;3、同一碱样品平行度很差。究竟是何因素影响其测定的准确性?笔者对此进行了研究。在这里我们引入最小二乘法，对所测得的大量数据进行统计、分析。二、最小二乘法的特点分析和统计思想1、特点分析在实际问题中,常常需要根据两个变量的几组实验数据来寻求两个变量之间的函数关系的近似表达式，通常把这样得到的近似表达式称为经验公式,下面介绍一种寻求两个变量之间近似成线性关系的经验公式的方法。2最小二乘法的统计思想设根据实验测得变量x与y的几组数据：x1x2x┄┄ix┄┄nxy1y2y┄┄iy┄┄ny在xOy坐标面上以iiyx,为点的坐标,相应地可以得到坐标平面上的n个点iiiyxP,ni,,2,1，如果这些点在一直线上,当然我们可以认为这两个变量之间的关系为线性的,即可求得.然而，一般说来，实际操作中这种可能性是很小的。在实际中,如果这些点的连线大致接近于一条直线，我们也就可以认为,变量y和x的关系是线性的,并可设baxy，其中a和b是待定常数。也即我们可以考虑用线性函数baxy近似地表达两个变量这间的函数关系。问题是，如何选择a和b才能使直线baxy更接近于以上各点连线,即希望直线上对应于ixxni,,2,1的点的纵坐标baxi与实验数据iy偏差都尽可能的小。为此,我们常常取),,2,1)((nibaxyii的平方和最小来保证baxi与iy间的差距尽可能的小。这样，问题也就归结为确定函数baxy中的常数a和b,使niniiiiibaxybaxybaW1122,为最小,以使直线baxy尽可能地靠近点iiiyxP,ni,,2,1的连线。.这种根据偏差的平方和为最小的条件来选择常数a、b的方法叫做最小二乘法。下面我们利用多元函数求极值的方法求出上最小二乘法中的ba,。对于函数niiibaxybaW12,因为niiiixbaxyaw12niiibaxybw12令00awbw即0011niiiiniiixbaxybaxy亦即niniiiniiiniiniiyxxbxaynbxa111211解上方程组得niniiiniininiiiixxnyxyxna1212111niniiiniiiniiniiniixxnyxxyxb121211112它还可简化为niiniiixxyyxxaxayb121其中niixnx11，niiyny11于是得经验公式baxy。三最小二乘法在纯碱中测定铁含量的应用1、数据的分析与应用序号1234567含铁质量（mg）0.000.010.020.040.060.080.10吸光度0.0000.0700.1260.2500.3750.5360.630根据刚刚绘制的曲线图,这些点大致在一条直线上，我们可以认为变量y和变量x的关系是线性的，由此可以设y=ax+b其中a和b是待定系数，根据最小二乘法提供的公式我们可以将已知数据列表计算ixiyxxiyyiyyxxii2xxi0.010.07-0.041-0.2610.01070.0020.020.126-0.031-0.2050.00630.0010.040.250-0.011-0.0810.00080.00010.060.3750.0090.0440.00040.000080.080.5360.0290.2050.0060.00080.100.6300.0490.2990.0150.0020.311.987331.0051.0yx将所求得数据代入公式可以求出a和b，得到y=6.6x-5.6×10-3这样给出x为铁的含量，y代表吸光度，由y=0.082、0.090(吸光度)进而求解x为0.013、0.014（铁的含量）。2应用效果分析我们选取了代表数据，初步分析、处理代入最小二乘法的计算公式中，求得经验公式。这样避免了人为主观绘制曲线，进而求得回归方程，减小了误差，使我们的测试数据更加科学合理化。而通过人为绘制的曲线，测得数据四、结论与认识在这里我们利用最小二乘法只是解决了一个项目，对于在以后的监测数据中只要两个变量成线性相关我们都可以用最小二乘法，利用最小二乘法可以把实验误差减小到最小，提高实验数据的准确度与可靠性。以此我们可以利用最小二乘法编程，设计成计算机软件，这样只需将我们测得的数据输入计算机中就可以得出经验公式，进而求得未知量。