[理学]数值计算方法第2章.pptVIP

一 误差的基本概念 1 误差的概念 一个量的准确值与其近似值之差。 2 绝对误差/相对误差 ① 绝对误差的定义：设x为正确值，x*为x的一个近似值，定义e(x*)=x-x*为近似值x*的绝对误差。 例1：设近似值5000的绝对误差为1，近似值5的绝对误差为0.1，则哪个近似值准确度高？ ③ 绝对误差界：设x为准确值，x*为x的一个近似值，如果能对x*的绝对误差做出估计|e|≤|x-x*|≤ε,则称ε为x*的绝对误差界。 3 有效数字 定义:设x的一个近似值x*,把x*写成规范化的科学计数形式x*=±0.a1a2…an…×10k(有限或无限),其中k为整数,a1,a2,…,an,…是0~9中的一个数字,且a1≠0,如果有 |x-x*|≤0.5×10k-n , 则称x*为x的具有n位有效数字的近似值。 例2：设x=103×0.7136……,问： ① 若取x*=103×0.714,则x*有多少位有效数字？ ②若取x*=103×0.715,则x*有多少位有效数字？ ③若取x*=103×0.7131,则x*有多少位有效数字？ 4 截断误差/舍入误差 ① 截断误差：也称为方法误差，它是指在构造数值计算方法时，用有限过程代替无限过程，其计算结果所存在的误差。 ② 舍入误差:也称计算误差,计算机所能表示的数字位数是有限的,因此需要对原始数据和计算结果进行”舍入“,这时所产生的误差就是舍入误差. 二 误差分析初步 1 数值计算中的重要公式Taylor公式 设f在含有点x0的某个开区间(a,b)内具有n+1阶导数,则 2 函数计算的误差估计 设 x为准确值, x* 为近似值,则 e(f(x*))=f(x)-f(x*) 3 算术运算的误差估计 例3:①设x0,x的相对误差为δ,求㏑x的误差. ②设x0,x的相对误差为2%,求xn的相对误差. 课堂练习 1 我们所熟悉的圆周率π=3.1415926……,问: ① 若取x*=3.14,则该近似值有多少位有效数字? ②若取x*=3.1416,则该近似值有多少位有效数字? 答案 1 解: π=0… ×101 ① |e|=|π-x*|=|0.0015926……| ≤0.002≤0.005=0.5×10-2=0.5×101-3 ∴ x*有3位有效数字 ② |e|=|π-x*|=|-0.0000073……| ≤0.000008≤0.00005=0.5×101-5 ∴ x*有5位有效数字 作业: 已测量得某矩形面积,长a=110m,宽b=8m,若|a-a*|≤0.1m, |b-b*|≤0.1m,试求其面积 s=ab 的绝对误差界和相对误差界. 算术运算的相对误差 三 计算机算术运算规则 1 加减法运算：首先比较加减两数的阶码，将阶码较小的数的阶码、尾数分别移位至两数阶吗相同后，尾数相加减。 2 乘法运算：阶码相加，尾数相乘，得两倍字长尾数，对其舍入成规格化的形式，并冠以乘积符号。 3 除法运算：将被除数的尾数扩大为双倍字长做除法，最后舍入为规格化形式，并冠以数符。 四 算法的数值稳定性 一个数值问题包括以下几个部分： 1 算法数值稳定性的概念 一个算法，如果执行过程中，初始数据误差和其他舍入误差的积累不增长，或者有增长但在一定条件下可控制，就称这个算法是数值稳定的，否则这个算法是不稳定的。 2 ①设由且仅由ε引起的，n步运算之后的误差为en，如果|en|≈cnε(其中，c是与n无关的常数）,即误差是线性增长的，从而可以通过ε来控制en，故算法是稳定的。 1 设x*的绝对误差很小，x为准确值，则e(f(x*))=f(x)-f(x*) ≈f’(x*)e(x*) 2 如果x≠0，x*充分接近x,则 六 数值算法设计原则 1 选择使用数值稳定性好的公式 使用数值稳定性不好的算法，每进行一次（一步）运算，误差都被放大，若干步后，误差可能已将真值“淹没”了，出现了错误的结果；而数值稳定性好的算式，在计算过程中，误差被逐渐缩小，若干步后，误差对结果的影响已经很小了。 2 避免两个相近数相减 3 防止大数“吃掉”小数 在数值计算中，参加运算的数有时数量级相差很大，又要对阶处理，如不注意运算次序，就可能出现大数“吃掉”小数的现象，影响结果的可靠性。 4 注意简化计算步骤，减少运算次数 减少计算次数不仅可以提高计算速度，而且能够减少误差的积累。 例如： 5 避免用很小的数（接近零）做分母 若用很小的数（接近零）做分母进行运算，会导致计算误差的扩散，给计算结果带来影响，而且，还可能造成“溢出”错误。 课堂练习 作业 问题本身 求解的算法 初始数据 （输入数据） 求解结果 （输出数据） 一个问题的求解结果是否可靠，可以由多种