工程数值分析-第1章-绪论.ppt

本章知识结构图 本章典型例题 题1: 指出如下有效数的有效数字位数并计算绝对误差限和相对误差限。 解: 1）x*有2位有效数字， 绝对误差限为： 相对误差限为： 2）y*有3位有效数字， 绝对误差限为： 相对误差限为： 解: 由题意知，近似数x*的绝对误差限 ，相对误差限 。 题2: 已知 试求 的相对误差限。 注意：此处正好有： 解: 题3: 已知桌子长宽近似值 ，并且已知 ，求近似面积 的绝对误差限和相对误差限。 题4：用最小刻度为毫米的卡尺测量直杆甲和直杆乙，分别读出长度a=312mm和b=24mm，问a,b的绝对误差限、相对误差限各是多少？两直杆实际长度x和y在什么范围内？ 解：由于绝对误差限是测量工具最小单位的一半，故 故：两直杆实际长度x和y的取值范围为 即 同理 注：本题x与y有相同的绝对误差限，但相对误差限不同，相对误差限越小越精确。 题5：设a=-2.18和b=2.1200是分别由准确值x和y经过四舍五入而得到的近似值，问： 各是多少？ 解：由于凡是由准确值经过四舍五入而得到的近似值，其绝对误差限等于该近似值末位的半个单位。故 因此相对误差限为 题6：下列近似值的绝对误差限都是0.005， 问：各个近似值有几位有效数字？ 分析 按有效数字的定义，根据绝对误差限，并将其取为某位的半个单位，从中就可看出有效数字的个数。 解：由题意，近似值的绝对误差限都是0.005，由于 故近似值精确到小数点后第二位，根据有效数字的定义 有三位有效数字1,3,7 有一位有效数字3 没有有效数字 题7： 设x=10±5%，试求函数f(x)=x1/n的相对误差限。 分析 这是标准的一元函数误差的传播问题，只需利用公式直接计算。 解 由x=10±5%知近似值x*=10，绝对误差限 由于 故 因此当x*=10时 即x*=10相对误差限 由 注：从本例可以看出，对于函数x1/n，函数值的相对误差限约是自变量相对误差限的1/n倍。 题8： 设 x=1.30±0.005, y=0.871±0.0005， 如果用 作为f(x,y)的近似值，则 能有几位有效数字？ 解： 而 故 所以 将数据代入 得 这说明 精确到小数点后第二位，故能有两位有效数字4，9。 分析 如果逐个相乘要用255次乘法。要 尽可能减少运算量，一种思路是尽可能 运用已经算出的结果。如2256=(2128)2， 当2128计算出来之后，再需一次乘法便可 得2256；而2128= (264)2，同样的，当264计 算出来后，再需一次乘法就可以计算出 2128；依次类推，这样就可得到一种乘法 次数少得多的算法。 题9： 请给一种算法计算2256要求乘法次数尽可能少。 解 这样共计算 8次乘法就可以计算出结果。 3）四则运算误差估计 两个近似数 a 与 b ，其误差限分别为 及 ， 它们进行加、减、乘、除运算得到的相对误差限分别为 例9：用四位有效数字计算如下y 的值 解：y的精确值是0.0158074374…， 算法1：直接计算，则 结果只有一位有效数字，其相对误差大于26%。 算法2：作适当变形后计算 结果具有四位有效数字，其相对误差不超过0.02%。 2）定义：数值稳定性 由于初始数据的误差或计算中的舍入误差在计算过程中的传播，因算法不同而异。一个算法，如果计算结果受误差的影响小，就称该算法具有较好的数值稳定性。 用一个算法进行计算，如果初始数据误差在计算中传播使计算结果的误差增长很快，这个算法就是数值不稳定的。 1.3.1 算法的数值稳定性概念 计算 并估计误差。 若计算出 ， 代入上公式，可逐次求出 的值。 例10： 由分部积分可得计算 的递推公式 注意此公式精确成立 1.3.1 算法的数值稳定性概念 ? 公式一： （1） 而要算出 就要先计算 . 并取 ， 则得 ， 计算过程中小数点后第5位的数字按四舍五入原则舍入. 若用泰勒多项式展开部分和 用4位小数计算， 截断误差 1.3.1 算法的数值稳