数值计算方法Chp1 预备知识2.pptVIP

1.3 误差 误差的运算与传播： 1.3 误差 初始误差通过一系列运算将进行传播。如果初始条件下的小误差对数值计算的最终结果影响很小，则所采用的算法称为稳定算法，否则，称为不稳定算法。 数学定义：设ε表示初始误差，ε(n)表示第n步计算后的误差增长。如果|ε(n)|≈nε，则称误差按线性增长。如果 |ε(n)|≈Knε，则称误差按指数增长。如果K1，则当n→∞时，指数误差的增长无界；如果0K1，则当n→∞时， 指数误差的增长趋于0。 1.3 误差 （1） 避免除数的绝对值远小于被除数的绝对值 若|x||y|，且x≠0，则 因此，当|x|相对太小时，商的绝对误差可能很大。 1.3 误差 （2）避免两个相近的数相减 若y ? x，设z=y-x，则ε(z)= ε(y)+ ε(x)， 因此，y ? x时，z ? 0，相对误差限会很大，导致结果有效数字位数减少。即，精度损失或差额抵消。 1.3 误差 例： 用6位有效数字和舍入法比较f(500)和g(500)的计算结果。其中 1.3 误差 真实值： 11.174755300747… 1.3 误差 （3）防止大数“吃掉”小数现象 计算机表示的数位数有限，很大的数和很小的数相加减时，很小的数会被“吃掉”（舍去）。 1.3 误差 （4）选用数值稳定性好的算法 有时候，模型本身就是病态（系数引入小变化，解产生大变化） 1.3 误差 O(hn)逼近： 设有函数f(h)和g(h),如果存在常数C和c，使 得对于任意h≤c，有 定义1.9： 则称函数f(h)为函数g(h)的big Oh函数，表示为f(h)=O[g(h)]。 big Oh提供了一个描述函数增长率的有效方法 1.3 误差 定义1.10： 则称序列{xn}以序列{yn}为上界，表示为xn=O(yn)。 设有两个序列 ,如果存在常数C和N，使得对于任意n≥N，有 1.3 误差 设函数f(h)的近似函数为p(h),且存在常数M0和正整数n，对于足够小的h,满足 定义1.11： 则称p(h)以近似阶O(hn)来近似f(h)，表示为 函数p(h)近似f(h)的误差限为： M|hn| 1.3 误差 定理： 设 1.3 误差 设f∈Cn+1[a,b]，如果x0和x=x0+h都在区间[a,b]内，则 泰勒定理： Some properties about O(hn): (i) O(hp) + O(hp) = O(hp); (ii) O(hp) + O(hq) = O(hr), where r = min{p, q}; (iii) O(hp)O(hq) = O(hp+q). 1.3 误差 例： 求如下泰勒多项式的展开的和与积的近似阶： 和 O(h4) 1.3 误差 积 O(h4) 1.3 误差 设limn→∞xn=x，有序列 ，且limn→∞rn=0。如果存在常数K0，对于足够大的n满足 序列收敛阶定义： 或表示为： xn→x，收敛阶为O(rn)。 则称 以收敛阶O(rn)收敛于x。可将其表示为 作业 P29: 2，3a，4a，5b、5d，9，13 1.3.10（1，2） 数值计算方法 E-mail： jdwang@uestc.edu.cn 办公室： 计算机大楼东206 电子科技大学物理电子学院 讲授： 王建东 教材、参考书 教材： 《数值方法》（Matlab版）（第五版） John H. Mathews, Kurtis D. Fink著（周璐等译），电子工业出版社 参考书： 《数值分析》（第四版） 李庆扬、王能超、易大义，清华大学出版社（2008） 教学内容、目的 主要内容： 第2—7章、第9章 授课： 40学时（+习题） 上机： 8学时 目的： 掌握一定的算法理论基础 会选用、修改算法用于具体问题 编制程序实现算法 成绩评定和作业要求 平时成绩： 30% 期末考试： 70% 作业要求： 作业本和以电子邮件（主题：姓名+学号）形式递交的程序（能直接运行并输出结果）。 注意： 学院统一出题考试？ 第一章 预备知识 主要内容 计算方法概述 微积分若干定理 误差（重点） 1.1 计算方法概述 实际问题（具体的科学、工程问题） 计算机 数值 计算方法 物理模型 近似解 数学模型 1.1 计算方法概述 计算方法： 是一种研究并