最小二乘法应用举例.docxVIP

PAGE PAGE 1 补充材料 1 实验数据的处理（上接教材第二章，p.19） 注意：（1）用最小二乘法计算斜率 k 和截距 b 时，不宜用有效数字的运算法则计算中间过程，否则会有较大的计算误差引入。提倡用计算器计算，将所显示的数值均记录下来为佳。（2）如果 y 和 x 的相关性好，可以粗略考虑 b 的有效位数的最后一位与 y 的有效数字最后一位对齐，k 的有效数 字与 y -y 和 x -x 中有效位数较少的相同。（3）确定有效位数的可靠方法是计算 k 和 b 的不确定度。 n 1 n 1 直线拟合的不确定度估算：（以 y ? kx ? b 为例） 斜率 k 和截距 b 是间接测量物理量，分别令测量数据的 A 类和B 类不确定度分量中的一个分量为零，而求得另一个分量比较简单，最后将两个分量按直接测量的合成方法求出合成不确定度，这种方法被称为等效法。 可以证明，在假设只有 y 存在明显随机误差的条件下（且 y 的仪器不确定度远小于其A 类不确 i 定度），k 和 b 的不确定度分别为： S ? x ? x 2 ? ( i y ? x )2 i n k 式中，S y S ? S ? ? x2 i n ?i? ? i ? 2 i n ? 2 ? S n? n? x 2 ? (? ? x2 i i x )2 i ?( ?( y ? kx ? b)2 i i n ? 2 S ? ? y 根据上述公式即可算出各个系数（斜率k 和截距 b）的不确定度值，初看上去计算似乎很麻烦， 但是利用所列的数据表格，由表中求出的那些累加值?即可很容易算得。 最小二乘法应用举例 应用最小二乘法处理物理量的测量数据是相当繁琐的工作，容易出现差错。因此，工作时要十分细心和谨惯。为便于核对，常将各数据及计算结果首先表格化。 例：已知某铜棒的电阻与温度关系为：R t ? R ? ? ? t 。实验测得7 组数据（见表1）如下：试用最 0 小二乘法求出参量 R 、? 以及确定它们的误差。 0 ?表 1t ? 表 1 t / ℃ 19.1 25.1 30.1 36.0 40.0 45.1 50.1 R / ? 76.30 77.80 79.75 80.80 82.35 83.90 85.10 ? ? x y 0此例中只有两个待定的参量 R 和 ，为得到它们的最佳系数，所需要的数据有 n、 、 、 0 ? ? ? i i x2 、 i y 2 和 i x y 六个累加数，为此在没有常用的科学型计算器时，通过列表计算的方式来 i i 进行，这对提高计算速度将会有极大的帮助（参见表2），并使工作有条理与不易出错。其中表内双 线右边的计算是为了确定 R 和?的误差项用的。 0 表 2 i i t / ℃ R / ? t t×t ( x ) ( y ) ( x2 ) R ? R t t ( y2 ) ( x y ) t×R t R 计算 / ? ? / ? i ? 2×10-4 i i i i i i i 1 19.1 76.30 364.8 5821.7 1457.3 76.26 +0.04 16 2 25.1 77.80 630.0 6052.8 1952.8 77.99 -0.19 361 3 30.1 79.75 906.0 6360.1 2400.5 79.43 +0.32 1024 4 36.0 80.80 1296.0 6528.6 2908.8 81.13 -0.33 1089 5 40.0 82.35 1600.0 6781.5 3294.0 82.28 +0.07 49 6 45.1 83.90 2034.0 7039.2 3783.9 83.75 +0.15 225 7 50.1 85.10 2510.0 7242.0 4263.5 85.19 -0.09 81 n ? ? x ??y ?i7 245.5i566.00 n ? ? x ?? y ? i 7 245.5 i 566.00 ? x2 ? i 9340.8 ? y 2 ? ? x y ? i 45826 i i 20060.8 ?? 2 ? i 2845×10-4 7 ? 20060.8 ? 245.5 ? 566.00 1472.6 ? ? k ? ? ? 0.28788?/ 0C 7 ? 9340.8 ? (245.5)2 5115.35 R ? b ? 566.00 ? 0.28788 ? 245.5 0 7 7 ? 70.76078? 把测量数据代入式（13）和（15）（参见教材第二章，p.19）中可求出相关系数 ? x y [? x 2 ? [? x 2 ? 1 (? x ) 2 ] ? [? y 2 ? 1 (? y ) 2 ] i n