现代科学与工程计算-孟大志-绪论1 2.pptVIP

* * * 0.2 数值计算的误差 误差在数值计算中是不可避免的，因此通过数值计算求出的解通常都是近似解，精确解与近似解之间的偏差就是误差。 0.2.1、误差与有效数字 1、误差的来源： （1）模型误差：数学模型与实际问题之间的差异；（忽略次要因素） （2）观测误差：观测值与真值之间的差异 ；(人或测量工具引起) （3）截断误差：近似解与模型精确解之间的差异 ;（算法本身引起） （4）舍入误差：按四舍五入方法取得的近似值与原始数据之间的差异 。 （计算机物理原因引起） 鲁东大学数 学 与 信 息 学 院 误差在我们日常生活中比比皆是， 在日常生活中你能举出几种误差 称为近似值x*关于准确值x的相对误差，记作 注：e*可正可负，为正时近似值偏大，称为强近似； 为负时近似值偏小，称为弱近似。 定义0.2.1 设x为准确值， x*为其近似值，则 称为近似值x* 关于准确值x的绝对误差，简称误差记作e* 鲁东大学数 学 与 信 息 学 院 过失误差：（人为因素引起） 由于模型错误或方法错误引起的误差.这类误差一般可以避免 除了过失误差可以避免外,其余误差都是难以避免的. 数值计算方法主要考虑和分析的是截断误差和舍入误差。 定义0.2.2 设x为准确值， x*为其近似值，则若 称ε*为x*的绝对误差限，简称误差限。若 由误差限可知 ，应用中常用 来刻画准确值x所在的范围。 称 为x*的相对误差限。 鲁东大学数 学 与 信 息 学 院 误差限的大小能不能完全表达近似值的好坏？ 1.45米的 绝对误差=？ 不知道！ 因为准确值x一般是未知，甚至无法知道！从而e*无法求出。 例:若用以厘米为最小刻度的尺去量桌子的长，大约 为1.45米，求1.45米的绝对误差。 鲁东大学数 学 与 信 息 学 院 绝对误差是有量纲的量，作为衡量x*的精度高低， 比较直观，但无法衡量精度的好坏。 相对误差是无量纲的量，也称百分比误差，衡量精度的好坏更合理。 绝对误差限和相对误差限均无穷多,自然越小越好.误差估计的任务就是提供好的误差限,对于任何一个近似值,如果得到一个好的误差限,那么就可以肯定这些数据是准确可靠的! 定义0.2.3 把x*表示成标准形式， 由准确值按四舍五入原则得到的数字是有效数字。有效数字中小数点后的位数越多，误差越小，计算结果越精确。 则称近似值 有n位有效数字。 每一位的数字 ，p是整数，p0时表示 整数部分的位数，p0时表示 小数点后零的个数，n是 正整数。若 鲁东大学数 学 与 信 息 学 院 例1. 各有几位有效数字？ ， 解： 故 有4位有效数字，而 有5位有效数字。 实际应用中也常用有效数字反映近似值的精确程度。 鲁东大学数 学 与 信 息 学 院 0.2.2、舍入误差 在电子计算机上，通常采用二进制实数系统， 并表示成类似十进制的浮点数形式： 称为机器数。其中2称为x的基数，整数m称为x的阶， 它用来确定该数小数点的位置。 字长，表示尾数的位数。L和U分别是m阶的下界和上界， 若在计算过程中m低于L或超过U，分别称为下溢或上溢， 下溢时机器数视为0，上溢时机器数视为无穷大。 小数 称为尾数，自然数t称为机器数的二进制 注意：要控制舍入误差的传播和积累，防止造成准确值被淹没！ 鲁东大学数 学 与 信 息 学 院 ξ介于x 与 之间，取绝对值得 忽略高阶项可得计算函数的误差限 当 时可得计算函数的相对误差 设函数y = f (x)， 有误差 ，估计y的误差可用Taylor展开： 0.2.3、函数值的误差估计 多元函数的误差估计可由多元函数的Taylor展开得到。 鲁东大学数 学 与 信 息 学 院 例2: 已知某场地长和宽分别为x = 110 m, y = 80 m, ε(x*) = 0.2 m，ε(y*) = 0.1 m，试求该场地面积的 绝对误差限和相对误差限。 解： 因为面积 ，所以误差限为 相对误差限为 鲁东大学数 学 与 信 息 学 院 0.2.4 误差的控制 1. 要避免两个相近数相减 在数值计算中，两个相近的数作减法时有效数字会损失。 当x = 1000，y 的准确值为0.01580 则求得的近似值为 y = 0.01581 ，有3位有效数字的精度。 例: 求 的值。 若直接相减 但若将其变形为函数 只剩下一位有效数字， 鲁东大学数 学 与 信 息 学 院 2. 要防止大数吃小数 　　当两个绝对值相差很大的数进行加法或减法运算时,绝对值小的数有可能被绝对值大的数吃掉从而引起计算结果很不可靠. 例如，在8位十进制计算机上计算 在计算机内利用浮点形式进行计