二章误差.ppt

计算物理 2.5 数值计算中需要注意的几个问题 方法：提高数值计算精度或采用别的算法 2.5 数值计算中需要注意的几个问题 3. 要避免用绝对值小的数去除绝对值大的数 数值计算中，除数的绝对值远小于被除数的绝对值，将会使商的数量级增加，甚至造成‘溢出’错误；而且当除数稍有一点误差，就会对计算结果造成很大的误差。如3.1416/0.001＝3141.6，当分母有了0.0001的误差时，也就是变为0.0011则商变为3.1416/0.0011＝2856，商的误差已经变的非常巨大 2.5 数值计算中需要注意的几个问题 4. 注意计算步骤的简化，减少算术运算的次数。 如果直接运算需作n(n+1)/2次乘法运算和n次加法运算。将表达式改写为下式只需作n次乘法和n次加法。 例如，计算多项式 一般来说，在一个物理问题中，通过计算方式的简化，既可以减少工作量，又可减小累计误差。 2.5 数值计算中需要注意的几个问题 5. 算法或公式要稳定 稳定性是指在数值计算中，误差的传播能否得到控制这样一个性质。 用一个算法进行计算，如果初始数据误差在计算中传播 使计算结果的误差增长很快，这个算法就是数值不稳定的. 计算 并估计误差. 由分部积分可得计算 的递推公式 若计算出 ， 代入，可逐次求出 的值. 例题 2.5 数值计算中需要注意的几个问题 而要算出 就要先计算 . 并取 ， 则得 ， 计算过程中小数点后第5位的数字按四舍五入原则舍入. 若用泰勒多项式展开部分和 用4位小数计算， 截断误差 2.5 数值计算中需要注意的几个问题 当初值取为 时，用（3.2）递推的计 算公式为 计算结果见表1-1的 列. 用 近似 产生的误差 就是初值误差， 它对后面计算结果是有影响的. 2.5 数值计算中需要注意的几个问题 2.5 数值计算中需要注意的几个问题 从表中看到 出现负值， 这与一切 相矛盾. 因此，当 较大时，用 近似 显然是不正确的. 实际上，由积分估值得 2.5 数值计算中需要注意的几个问题 计算公式与每步计算都是正确的，计算结果错误的原因 主要就是初值 有误差 ，由此引起以后各步 计算的误差 满足关系 容易推得 这说明 有误差 ，则 就是 的 倍误差. * * 第 二 章　误差 2.1　误差的来源与分类 2.2 误差与有效数字 2.3 函数的误差估计 2.4 近似数的四则运算及数值计算中需要注意的几个问题 从实际工程问题出发,一直到算出问题结果,这其中的每个过程都会产生误差。计算物理中的误差来自四个方面: 模型误差、观测误差、方法误差（截断误差、舍入误差） 计算物理中的误差 2.1 误差的来源与分类 （a）模型误差。将实际问题归结为数学问题时，总要忽略一些主观上认为是次要的因素，附加若干限制。例如点粒子近似。人们对客观事物的认识是逐步深入的，这样建立的“理想化”的数学模型，虽然具有“精确”而“完美”的外表，实质却只是客观现象的近似而粗糙的描述，这种近似描述就隐含着误差，这就是模型误差。如自由落体运动忽略了空气的阻力。 2.1 误差的来源与分类 （b）观测误差。在数学模型中，往往包含有若干参变量，如物体密度，物态方程与本构方程参数，热量交换系数等等。这些参量一般是通过实验观测确定的，因而不可避免会存在观测误差。如自由落体运动中的时间和重力加速度就是观测值。观测值的精度依赖于仪器和人的操作。 2.1 误差的来源与分类 （c）方法误差。在实际解题过程中，数学模型常常比较复杂，不能获得精确解，要用数值方法得到计算结果，模型的准确解和数值方法的准确解之差称为方法误差。 2.1 误差的来源与分类 例如：指数函数ex可展开成下列幂级数形式 但在实际计算时，不可能计算无穷多项，只能截取有限项： 用Sn(x)作为ex的近似值，其截断误差为 2.1 误差的来源与分类 截断误差属于一种方法误差 由于计算机只能进行有限位运算，所以在运算过程中会不断地按某种规则进行舍入，此即舍入误差。舍入误差实质上也是方法误差。 计算物理课程中关注的就是方法误差（截断误差、舍入误差） 2.1 误差的来源与分类 用有限项近似无限项时,因为截去部分项而产生的误差，称为截断误差 用有限位近似无限位时,由舍入产生的误差，称为舍入误差 绝对误差 2.2 误差与有效数字 设 为准确值， 为