数值分析教学资料 第一章 引论.pptVIP

对于多元函数 设准确值 由多元函数Taylor公式，可得误差估计： 相对误差限为： 2. 算术运算的误差估计： 两个近似数x1, x2，其误差限分别为ε(x1)， ε(x2)，它们进行加、减、乘、除运算得到的误差限分别为： 例4: 设a=2.31,b=1.93,c=2.24都是三位有效数字的近似数，令p=a+bc，求ε(p)和 εr(p)，并判断p有几位有效数字。 解 由题知， ε(a) =ε(b) =ε(c)=0.005 ε(p)= ε(a)+ ε(bc)? ε(a)+|b| ε(c)+|c| ε(b) =0.005+1.93?0.005+2.24 ?0.005=0.02585 p=a+bc=2.31+1.93 ?2.24=6.6332 又 故 εr(p)=ε(p)/|p| ? 0.02585/6.6332 ? 0.0039=0.39% 因为 ε(p) ? 0.025850.05=1/2 ?101-2 所以 p=6.6332中只有两位有效数字。 1. 避免除数的绝对值远小于被除数的绝对值； 2. 避免两个相近的数相减； 3. 防止大数“吃掉”小数现象； 4. 简化计算步骤，减少运算次数； 5. 选用数值稳定性好的算法。 1.4 数值运算中的一些基本原则 1. 避免除数的绝对值远小于被除数的绝对值； 设 若|x||y|， 因 表明当|x|相对太小时，商的绝对误差可能很大。 2. 避免两个相近的数相减； 若y ? x，设z=y-x，则ε(z)= ε(y)+ ε(x)， 当y ? x时，z ? 0，相对误差限会很大，导致结果有效数字位数减少。 避免两个相近的数相减方法：常用一些恒等变形实现。 如当x充分大时， 当x绝对值很小时， 3. 防止大数“吃掉”小数现象； 原因：计算机表示的数位数有限，很大的数和很小的数相加减时，很小的数会被“吃掉”（舍去）。 方法：当绝对值悬殊的一系列数相加时，若有 |x1||x2|…|xn|，应按绝对值由小到大的顺序累加。 例 计算方程 的根。 算法一： 利用求根公式 在计算机内， 109存为0.1? 1010 ，1存为0.1? 101。做加法时，两加数的指数先向大指数对齐，再将浮点部分相加。 即1 的指数部分须变为1010 ，则：1 = 0.0000000001 ? 1010 ，取单精度时就成为： 109 +1=0 1010 +0? 1010 =0? 1010 算法二： 先解出 再利用 例2（思考题）： 按从小到大、以及从大到小的顺序分别计算 1 + 2 + 3 + … + 40 + 109 4. 简化计算步骤，减少运算次数； （即减少计算工作量） 如： 计算 若直接求和，计算第k项时需n-k+1次乘法，故共需n(n+1)/2次“×”和n次“＋”。 若将公式改写为： 则只需n次“×”，n次“＋”，即秦九韶算法。 5. 选用数值稳定性好的算法。 将输入数据有误差，但在运算过程中舍入误差不增长的算法称数值稳定的，否则是数值不稳定的。 只有数值稳定的数值方法才能给出可靠的计算结果。 例： 设有递推关系： (n=1,2,…) 若取 （三位有效数字） 试问计算到y10时误差有多大？ 解： 因 而 于是有 类推，有 即计算到y10，其误差限为1010 ?，亦即若在y0处有误差限为?，则y10的误差限将扩大1010 倍，可见这个计算过程是不稳定的。 P17－18. 习题一：1，2，3，4，5，8。 本章作业 《计算方法》 第一章 引论 数学科学学院 计 算 方 法 电子科技大学数学科学学院 教材和参考书 教材： 数值分析，电子科技大学应用数学学院，钟尔杰，黄廷祝主编，高等教育出版社 参考书： 数值方法（MATLAB版）（第三版），John H. Mathews,Kurtis D. Fink 著，电子工业出版社； 数值分析（第四版），李庆扬，王能超，易大义编，清华大学出版社； 计算方法（第二版），邓建中，刘之行，西安交通大学出版社； 数值分析与实验学习指导，蔡大用编，清华大学出版社。 教学要求 了解计算方法研究的主要内容； 掌握计算方法的基本概念和基本原理，进一步提高抽象思维和逻辑推理的能力； 掌握数值计算的各种方法（或算法）的基本思想，进一步提高数值计算能力 ； 能够与实际问题相结合，利用所学算法解决一些实际的数学模型问题 ； 能够利用数学软件编程实现所学算法（可用MATLAB，MATHEMATICA等）。 课程内容 引论（误差、有效数字等） 非线性方程求根 线性方