实验数据的处理方法解析

1、实验数据的处理方法【摘要】物理学是一门实验的科学，物理学中的新概念、新规律的发现都依赖于 反复的实验。而处理实验数据时，需选择适当的实验数据处理方法，才能较准确、 客观的反映实验结果，减小误差。本文介绍了实验数据处理中涉及到的一些基本 概念，重点综述了物理实验中常用的数据处理方法。 并指出了各自适用的条件及 优缺点。【关键词】误差；数据处理；作图法；最小二乘法；逐差法Abstract: Physics is an experimental scienee. New concepts in physics, the discovery of new rules rely on trial and

2、 error, The experimental data processing Need to select the appropriate treatment of the experimental data, To more accurately reflect the objective results， Reduce errors. This article describes the experimental data process ing invo Ived in some of the basic con cepts Summary of experime nts focus

3、ed on the physical data processing methods commonly used. And pointed out the adva ntages and disadva ntages of each applicable con diti on.Keywords: Error; Data Processing Mapping ； Least squares By subtraction【引言】数据处理是指由实验测得的数据，必须经过科学的分析和处理，才能揭示出各物理量之间的关系。我们把从获得原始数据起到得出结论为止的加工过 程称为数据处理。正确的处理实验记录的数

4、据，对我们科学的了解被测量或研究 对象的客观规律，选择恰当的实验数据处理方法，最大限度的减小误差让实验数 据无限接近理想条件下的结果，这是实验数据处理的意义所在。在这方面研究的 文献有很多，例如费业泰的误差理论与数据处理等。要对实验结果进行分析， 根据不同的实验方法，我们可以采用不同的数据处理方法，常用的有作图法、最小二乘法、逐差法等。本文将分别对这些方法进行了介绍。、实验数据处理中涉及到的基本概念对实验结果进行半定量分析， 需要借助许多评价参量， 这里涉及到许多表征 实验数据好坏及数据离散程度的基本概念， 特别是描述处理后数据的可靠性的参 量，尤其具有重要意义。以下将对一些重要概念进行介绍。

5、1 . 真值及约定真值真值有多种定义，如 “被测量本身所具有的真实大小称为真值。 ”1“如果实 验已消除系统误差， 只存在偶然误差。 则无穷多个观测值得平均值， 就是被测物 体的真值。 ”2 “正在研究某量时所处的条件下严格的确定的量值。 ”3由此可见， 真值是客观存在的，但也还是一个理想的概念，通常是不可确切知道的。约定真值被认为是非常接近真真值的， 它们之间的差别可忽略不计， 无系统 误差条件下的算术平均值、标称值、校准值、理论值、公认值等均可作为约定真 值来使用。2. 影响量和干扰量影响量不是测量的对象， 但却影响被测量的量值或仪器示值， 它通常是一种 与待测的量有一定函数关系的另一种性

6、质的量 9。例如在测量电阻时，由于多数 材料的电阻随温度改变， 因此在测量电阻时温度在影响着电阻值的测量结果， 但 它却不是测量对象， 所以温度就是电阻测量中的影响量。 影响量在测量结果中带 来的影响可以在测出影响量的大小后，按其函数关系从测量结果中加以消除。干扰量是一种与待测的量没有必然联系的外界强行渗入量。 例如拍摄全息照 片时外界的振动， 探测器的噪声， 都会影响测量工作造成干扰， 这样的一些量称 为干扰量 9。为了保证测量的准确度，在安排测量条件时，要消除影响量和最大 限度减小干扰量。3. 精度反映测量结果与真实结果接近程度的量， 称为精度，它与误差的大小相对应， 因此可用误差的大小来

7、表示精度的高低，误差小则精度高，误差大则精度高低。精密度：它反映测量结果中随机误差的影响程度 准确度：它反映测量结果中系统误差的影响程度 精确度：它反映测量结果中系统误差和随机误差综合的影响程度， 其定量特 征可用测量的不确定度（或极限误差）来表示。 4 用一种打靶的例子，可以更好 的理解和掌握，如下图：(a)精密度高、准确度低图1打靶的精确度(c)精确度高4. 误差测量值减去真值为测量值的误差，即：测量值-真值二误差上述误差亦称绝对误差，假设测量值为X，真值为a，误差为；，则有误差；=x-a，；与a的比值；= /a称为相对误差。对误差的来源可以概括为四个方面4: 1,测量装置误差，这里面分三

8、个方面来说，a,标准量具误差：以固定形式复现标准量具的器具，如标准量块、标准先问吃、标准电池、标准电阻、标准砝码等， 它们本身体现的量值，不可避免的都还有误差。b,仪器误差：凡用来直接或间接将被测量和已知量进行比较的仪器设备，称为仪器或仪表，如天平、压力表、 温度计等，它们本身都具有误差。c,附件误差：仪器的附件及附属工具，如测 长仪的标准环规，千分尺的调整量棒等的误差，也会引起测量误差。2,环境误差：由于各种环境因素与规定的标准状态不一致而引起的测量装置和被测量本身 的变化所造成的误差，如温度、适度、气压、振动、照明、重力加速度等所引起 的误差。通常仪器仪表在规定的正常工作条件所具有用的误差

9、，而超出此条件时所增加的误差称为附加误差。3,方法误差：由于测量方法不完善所引起的误差， 如测量一个轴的直径d，因近似数二取值的不同，将会引起误差。4,人员误差：由于测量者受分辨能力的限制或是崛起变的生理变化，固有习惯引起的读数误差，以及实验室的疏忽等所引起的误差。总之，在计算测量结果的精度时， 对上述四个方面的误差来源，必须进行全面的分析，力求不遗漏，不重复，特别 是对误差影响较大的哪些因素。4.1系统误差先看两个例子51）用一个2.5级0 1A的安培计测一回路的电流强的I为 0.73A，而用另一个0.5级0 1A的安培计测同一回路电流为 0.716A; （2）用一 天平称一物体的质量，物体

10、在左盘，砝码在右盘，平衡时测量值为 74.2519g,物 体与砝码交换后则为74.2501g；（ 1）是由于仪器自身误差的问题，（2）是由于天 平左右臂长不完全相等引入的系统误差，可将物体放在天平左、右盘上各称一次 取平均值去消除。上述各项测量值的差异在重复测量时依然不变，这表示其误差的符号和大小 是恒定的，此类误差称为系统误差。对系统误差的研究主要是：（1）探索系统误差的来源，设计实验方案消除或消减该项误差。（2） 估计残存系统误差的可能的范围。54.2偶然误差在同一条件下，对同一物理量进行重复测量，各次测量值一般不完全相同， 这是由于测量时存在的偶然误差。一个测得值的偶然误差是多项偶然因素

11、综合作 用的结果，在测量前不能得知测得值将偏大或偏小。图（2）.单摆周期测量100次的统计直方图这里以测单摆周期的实验为例：用手控制数字毫秒计，测量一摆的周期共100次，测量值的大小变化不定，似乎没有规律，其实这种偶然现象服从统计规 律。现将测得值分布的区域分为 9个区间，统计各区间内测量的个数 Ni，以测 量值为横坐标，NN为纵坐标（N为总数）作统计直方图，图2是一次实验的 结果。从图上可以看出，比较多的测量值集中在分布区域的中部， 而区域的左右 两半的测量值个数都接近一半，由此可以设想被测量的真值就在数据比较集中的 部分。在上述测量之后，用光电门控制一台数字毫秒计去测同一个摆的周期，测10

12、次，测量值分布在1.866s到1.868s的小区域中，由于此时的偶然误差显著小 于前者，可将光电控制测量的平均值To作为手控测量的近真值，对于测量值的偶 然误差作如下的统计，取T=1.8670s,则T-T0v0 (&_ 0)占 48%Ti-To _0 ( 0占 52%多次测量均有同上相似的结果，因而得出如下几点认识：(1) 每次测量的偶然误差是不确定的。(2) 出现正号或负号偶然误差的机会相近。(3) 出现绝对值小的偶然误差的机会多一些。5 .算术平均值与标准偏差5.1算术平均值设n次测量值X1，X2，,Xn的误差为；1，；2，；n，真值为a，则(X1-a)+(X2-a)+(Xn-a)=+之+

13、计将上式展开整理后，两侧除以n,得, 、 1 ,；n)(X1 + X2+ Xn)-a=(；1+；2+ +n它表示算术平均值的误差，等于各测量值误差的平均值，假如各测量值的误差只 是偶然误差，而偶然误差有正有负，相加时可抵消一些，所以 n越大，算术平均 值越接近真值。因此可以用算术平均值作为被测量真值的最佳估计值。又当测量值的误差中包含有已知的系统误差，则相加时它们不能抵消，这时 应当用算术平均值加上修正值为被测量真值的最佳估计值(修正值与系统误差绝 对值相同，符号相反)。5.2标准偏差具有偶然误差的测量值将是分散的，对分散的情况的定量表示用标准偏差 s，它的定义是为s=n_(X -X)2i =

14、1n为测量值个数。例，比如有如下两组数值:表1.两组数值集中程度不同的实验测量数据A2.12.62.82.93.03.23.7B2.12.42.72.93.13.43.7两组数值都在2.1到3.7之间，平均值都是2.9,计算标准偏差为sa=0.4970.50Sb=0.5570.56可以看出A组数比较像中间集中，B组数则稍差，表现除它们 分散上的差异。6 有效数字在做实验时总要记录很多的数据，并对数据进行计算或处理，但在记录时应 取几位，计算后应保留几位，这是实验数据处理的重要问题，必须有一个明确的 认识。实验时处理的数值，应能反应出被测量的实际大小的数值，即记录与运算后保留的应为能传递出被测量

15、实际大小信息的全部数字，这样的数字称为有效数字。例如用一最小分度1mm的尺，测得一物体的长度为7.62cm,其中7和6是 准确度出来的，最后一位数字2是估读的，并且仪器本身也将在这一位出现误差， 即这一位不一定是2,只是近似的，但是还是一位有效数字。在实际取舍时按照 实验条件以及题目要求为参考。使用有效数字规则时的注意事项：5 （1）物理公式中的有些数值，不是实 验测量值。例如，测量圆柱体的直径 d和长度I求体积V的公式V =丄二d2l中的4-不是测量值，在确定 V的有效数字位数时不必考虑-的位数。（2）对数运算44时，首数不算有效数字。（3）首位数是8或9的m位数值在乘除运算中，计算 有效数

16、字时，可多算一位。（4）有多个数值参加运算时，在运算中途应比按有效 数字运算规则规定的多保留一位，以防止由于多次取舍引入计算误差，但运算最 后仍应舍去。（4）数值的修约规则：开始要舍去的第一位数是 1、2、3、4时 就舍去；是6、7、& 9时在舍去时进1。要舍去的一位是5,而保留的最后一 位是奇数，则舍去5进1,是偶数则舍去5不进位，但是5的下一位不是0是仍 然要进位。5、物理实验中常用的数据处理与分析方法实验数据处理的方法有很多，对不同的实验记录数据要选择不一样的数据处 理方法，以下介绍几种常用的实验数据处理方法， 并最后对各自的优缺点进行总 结。1.列表法把数据按一定规律列成表格，可使物理

17、量之间的一一对应关系简明、醒目， 也有助于发现实验中的规律。列表时应注意：（1）表格设计合理、简单明了，便 于观察。（2）各栏目中均应注明物理量的名称和单位。（3）各量排列顺序尽量与 测量数据顺序一致，用有效数字填写。特殊需要可以采用其他规律。这种方法在实验室常用到，如果将实验要记录的数据在表格中一一罗列出来，这样清晰明了，在处理数据是一目了然。如表2表2某物理量的测量物理量名称测量值1232.分组计算法设变量x、y间存在y=a+bx的直线关系，由于测量存在误差，对n次测量有yi = a + V|y2 = a + bx2 + v2 yn 二 a bxn vn式中Vi表示测量误差。现在将n组测量

18、分为前后两部分，从中取对应的两组:略去误差项，解出a、b的近似值:= a + bx=a bx v煤i-2Jy n ybi =i 2x -xi .2(y+yJ b(x + xi)aiii访I2J这样可的-个ai和b，再求a、s(a)和b、s(b)。如n为奇数，中间数可公用23 .最小二乘法5假设变量x、y间存在直线关系y=a+bx，参量a、b分别为y轴截距和斜率，当将测量值(x，yJ带入此式时，由于存在测量误差 yi =a bxi，引入误差项vi后有yi =a+bXj +比， 或 v = yi(a+bxj，对n次测量，可有W = % -(a+bxj V2 = y2 _(a+bx2),Vn = y

19、- -(a+bXn)由于n个方程中有n+2个未知数，所以不能从解联立方程组求出 a、b值。设y 为等精度测量值，最小二乘法原理指出，满足a v：二极小值条件下求出的参量a、b之值、为最佳拟合值，即从A V：y(a b)2AA、 yi -a b A xi解此联立方程组，得 工 Xi yi = ax + bXi2A _、x YXxyian xi -C xj2A n二 xy 二 Xi 二 Y b 二 n x2-C x)2Sxx三E （Xi-声乏艮臣（n、令 Syy三瓦（yi匚$乏y臣（yn；Sxy三瓦（Xix（%y）乏 x y更咨 / y n JAAA则由及，又可得出b = Sy / Sxx , a

20、 =y - bx 为了反映变量X、y间的线性关系的密切程度，常用关联系数r来描述，其估算式为二 （xi；x）（yi JSxy送（Xi X）送（Yi-y） 12 JsxxSxyA从理论上讲，r0就应该承认x、y之间存在一定的相关关系，但是由于值A是从较少的数据中求出的，根据数理统计理论，对于一定的n值，r要在大于某一临界值r临时，才可以认为存在线性相关关系。下表给出了各n值的r临值。表3各n值的r临值5n3456789r临0.9980.9900.9500.9170.8740.8340.798n10111213141516r临0.7650.7350.7080.6840.6610.6410.623n

21、171819202122r临0.6060.5900.5750.5610.5490.537b的标准偏差Sa、Sb和y的标准偏差Sy之间的关系为A AASyf （八葺浒112 - 匕弓如12结合和二式,n 2n 24 作图法研究两个变量的关系就是用曲线或函数式将二个变量的联系变现出来。 那么 首先在这里先介绍实验图线的描绘。常用直角坐标纸为方格纸，制图线时需注意 的问题：(1) 标的横轴为自变量，纵轴为因变量。一般是以被测量为变量，但有时为了使获得的图线是一条直线，而将被测量做某种变换后的数值做变量。 这种变换 不仅是由于直线便于描绘，更重要的是直线的斜率和截距所包涵的物理内容是我 们所需要的。(

22、2) 标的原点，不一定要和变量的零点一致，若变量x的变化范围是从a到b， 则将坐标原点在取在a的附近即可。因为有时a距x的零点很远，如果将原点取 在x的零点出，贝U坐标纸上将出现很大的空白区域，白白浪费了坐标纸。(3) 标轴的分度要和测量的有效数字位数对应，坐标纸的一小格表示为被测量的最后一位的一个单位、两个单位或五个单位比较好，要避免用一小 格表示三、七或九个单位。因为那样不仅坐标点和读数都不方便， 也容易出现错 误。(4) 和y轴二变量的变化范围(a b) (c d)，表现在坐标纸上的长度应该相差不大，最多也不要长过一倍。注明 x、y轴代表的测定量及单位，按 测量数据标出坐标点。4.1在我

23、们讨论变量间的函数关系，这类问题有两种不同的情况：(1) 已知两个变量函数关系的形式，但是其中有位置参量；如果两个变量x、y之间是直线关系，即y=a+bx 贝U可用n组测量值(x，yj 作图，所得直线的截距即参量 a,斜率是参量b。但是实验中两个变量的关系往 往不是直线，例如，弹簧振子的振动周期T和所负载m的关系为 T =2二m囚 ，式中mo为弹簧自身的重量，c和k是待定参量。测量不 同m对应的T,可以作Tm图线，图3为一例。由于它是曲线，因而无法从图 上得出待定参量值。改换变量，将函数关系转变为直线关系，对此周其公式可以 改成 t2 =(4 严) 毎 m ，他表示T2与m间为直线关系，作T2

24、-m直线（图4）,从图上求出截距a,斜率b,则a=严,b年由此可求 出k和c值。即对于非线性函数，要通过变换变量使之成为线性函数，再用作图法求出截距和斜率，进而确定待定参量图弹簧振子的T m曲线图（4）弹簧振子的T2 m曲线非线性函数如何变换要看函数的形式，例如:y 二 aex In y二 ln a _ bxE -x1y 二Ry 一R R Exxy = ax bx2 = a bxx上列式中的括号为变换后的变量。（2）两个变量函数关系的形式尚未知时首先用测量值作图，如果得一直线，则从图上求出截距和斜率，函数关系就完全 确定了；当得到的是曲线时，就要分析曲线的形式，参照已知的函数曲线，给出 假定的

25、函数式，再用上述（1）中处理非线性函数的方法，使之线性化，但这样 做不一定一次就成功，可能要反复几次才可得出较好的结果4。通过电阻伏安特性研究的实验为例，对作图法和最小二乘法进行比较，其公 式为R = U，式中U为电压，I为电流，R微电阻，测量数据见表4,用所测量的结果画图如图（5）表4电阻伏安特性测量表物理量12345678910U（ V）0.6220.7180.8230.9221.0211.1341.2351.3231.4121.5261(A)0.2340.2640.3050.3310.3690.4210.4610.4790.5210.565图（5）电阻的伏安特性曲在直线上取两点，A （0

26、.600,0.223）、B （1.500,0.553）,利用作图法进行处理， 求得斜率K图曰 0.553 - 0.223 “a?，求K图值的标准偏差，由于U2-U,1.500-0.600K图=0.367,所以作图法拟合直线为I=0.367U,将表4中的U代入得表5表5电阻伏安特性也I值物理量12345678910U（ V）0.6220.7180.8230.9221.0211.1341.2351.3231.4121.526I =0.367U0.2280.2630.3020.3080.3750.4160.4530.4840.5180.560I(A)0.2340.2640.3050.3310.369

27、0.4210.4610.4790.5210.565曲 m _i0.0060.0010.003-0.007-0.0060.0050.008-0.0060.0030.005I的偏差6上也120.008別2 0 0 5怎6 K0.的）标6准偏差Vi# n1 V1-01、k图：:一n 门 n 、i100.0057=0.0062图-ur1012.36一115.26 i =1i =1下面再估算读取数据时产生的误差.由于坐标纸最小分格之间人眼无法辨别具体数据，因此存在视觉误差U、寿,K图2匕U2_S推出：仝图二 巴 3 . g -UJ 亠 .巴 W，假设每次读取的误 K 图I2J1 U2_Uil2_llU2

28、_Ui差相同，且为坐标纸最小分格的一半则.U =0.005，J =0.0025代入数据得：0.0025 0.0025 0.005 0.005、.K图二0.367 () =0.0096、0.01可见作图法在图上图0.553 0.2231.5000.600读取数据时的误差是较大的。若坐标分度值取得不当，求斜率K图时，从直线上取的两点距离较近的话，产生的误差就更大。利用最小二乘法处理数据，求斜率K乘值，由U =1.0736，=1.2356,2U =1.1526，I =0.395，U I =0.4546, U I =0.424,得：U I U I 0.4546 0.424十 / /古砧 士一 、吐 /

29、已辛 rK乘20.367 ：0.37。求K乘值的标准偏差，由U2 _U21.2356 -1.1526K乘 =0.37得到最小二乘法拟合曲线为I=0.37U，将表4中的U代入的表6K乘的偏差nn-e u)2i =1M0 12.3-115.26 O.。059 ： 0.。07小结：1图=0.367 0.01, K乘=0.367 0.007,即最小二乘法拟合的直线较精确。i =1表6数据列表物理量12345678910U( V)0.6220.7180.8230.9221.0211.1341.2351.3231.4121.526I =0.367U0.2300.2660.3050.3410.3780.42

30、00.4570.4890.5220.565I(A)0.2340.2640.3050.3310.3690.4210.4610.4790.5210.565.J=I -i0.004-0.0020.000-0.010-0.0090.0010.004-0.0100.0010.000I 的偏差-丄 n21101 O.00319 ： 0.0059 ：。.0065 .逐差法逐差法是为了减小系统误差，而在实验数据处理当中常用的一种方法。 当实 验中的两个物理量满足正比例关系时， 依次记录改变相同的量时的值、 或者某一 研究对象随实验条件周期性变化时，依次记录研究对象达到某一条件（如峰值、 固定相位等）时的值、逐项逐差再求平均 8。但是实验数据的个数却对逐差法有 着影响，是偶数对数据（n=2m）时，我们会很容易地算出斜率；当是奇数对数 据 （ n=2q-1） 时 需 要 将 中 间 数 据 对 （ n=q） 舍 弃， 例 如 11 对 数 据 （n=11,q=6,m=n/2=5.5）