上海交大数值分析课件数值分析1-3.ppt

故两直杆实际长度x和y的取值范围为 即 同理 注：本题x与y有相同的绝对误差限，但相对误差限不同，相对误差限越小越精确。 例2 设a=-2.18和b=2.1200是分别由准确值x和y经过四舍五入而得到的近似值，问： 各是多少？ 解：由于凡是由准确值经过四舍五入而得到的近似值，其绝对误差限等于该近似值末位的半个单位。故 因此相对误差限为 例3 下列近似值的绝对误差限都是0.005， 问：各个近似值有几位有效数字？ 分析 按有效数字的定义，根据绝对误差限，并将其取为某位的半个单位，从中就可看出有效数字的个数 解：由题意，近似值的绝对误差限都是0.005，由于 故近似值精确到小数点后第二位，根据有效数字的定义 有三位有效数字1,3,8 有一位有效数字3 没有有效数字 例4 按四舍五入原则写出下列各数具有5位有效数字的近似值： 解：由于有效数字的计算是从第一位非零数字算起，故 具有5位有效数字的近似数是 例5 设x=10±5%，试求函数f(x)=x1/n的相对误差限。 分析 这是标准的一元函数误差的传播问题，只需利用公式直接计算。 解 由x=10±5%知近似值x*=10，绝对误差限 由于 故 因此当x*=10时 即x*=10相对误差限 注：从本例可以看出，对于函数x1/n，函数值的相对误差限约是自变量相对误差限的1/n倍。 由 例6 设 x=1.30±0.005, y=0.871±0.0005， 如果用 作为f(x,y)的近似值，则 能有几位有效数字？ 解 而 故 所以 §3 误差定性分析与避免误差危害 一、算法的数值稳定性概念 二、设计算法的若干原则 1.定义：算法 所谓算法，是指对一些数据按某种规定的顺序进行的运算序列。 一、算法的数值稳定性概念 对同一问题，选用不同的算法，所得结果的精度往往大不相同。 请看下例 例：用四位有效数字计算如下y的值 解：y的精确值是0.0158074374…， 算法1：直接计算，则 结果只有一位有效数字，其相对误差大于26%。 算法2：作适当变形后计算 结果具有四位有效数字，其相对误差不超过0.02%。 2. 定义：数值稳定性 由于初始数据的误差或计算中的舍入误差在计算过程中的传播，因算法不同而异。一个算法，如果计算结果受误差的影响小，就称该算法具有较好的数值稳定性 二、设计算法的若干原则 (一) 要避免相近两数相减 例：a1 = 0.12345，a2 = 0.12346，各有5位有效数字。 而 a2 ? a1 = 0.00001，只剩下1位有效数字。 ε很小， ? 几种经验性避免方法： 很小， (二) 要防止大数“吃掉”小数，注意保护重要数据 例：用单精度计算 的根。 精确解为 ? 算法1：利用求根公式 在计算机内，109存为0.1?1010，1存为0.1?101。做加法时，两加数的指数先向大指数对齐，再将浮点部分相加。即1 的指数部分须变为1010，则：1 = 0.0000000001 ? 1010，取单精度时就成为： 109+1=01010+0?1010=0?1010 大数吃小数 算法2：先解出 再利用 注：求和时从小到大相加，可使和的误差减小。 例：在5位浮点十进制计算机上，计算y = 54 321 + 0.4 + 0.3 + 0.4 解：若按从左到右的顺序进行计算，后三位在对阶过程变为 后三个数都在对阶过程中变为零，得出含有较大误差的结果 y = 54321。 但若按从右到左的顺序进行计算，后三位在对阶过程变为 这种算法避免了大数“吃掉”小数！ 一般地，有如下原则 若干数相加，采用绝对值较小者先加的算法，结果的相对误差限较小 (三) 注意简化计算步骤，减少运算次数，避免误差积累 例：计算多项式的值 解：如果先计算各项然后相加，则 乘法次数=4+3+2+1=10，加法次数=4 但如改用下式计算 则只需做4次乘法和4次加法。计算量大大减少！ 注：第二种方法称为“秦九韶算法” 通常，计算如下n次多项式的值 如果先计算各项然后相加，则 乘法次数 =n+(n-1)+···+2+1 = n(n+1)/2 加法次数= n 若采用“秦九韶算法”，则 乘法次数= n 加法次数= n 两种算法的乘法运算次数随n的变化见下表： 5 4 3 2 方法2 15 10 6 3 方法1 n=5 n=4 n=3 n=2 (四) 要避