数值分析引论要点.pptVIP

计算方法 郑海标 数学学院 Email：hbzheng13@mail.xjtu.edu.cn 理论分析 科学实验 科学计算 计算的目的不在于数据，而在于洞察事物。 －－理查德·哈明 The purpose of computing is insight，not numbers. －－Richard??Wesley??Hamming 计算方法课程简介 教材 《数值计算方法》 李乃成 邓建中编 参考书目 《数值计算方法》 冯康等编 《数值分析简明教程》 王能超编 《Numerical Analysis》 (Seventh Edition) 《数值分析》 （第七版 影印版） Richard L. Burden J. Douglas Faires （高等教育出版社） 计算方法学什么？ 计算方法=数值分析=科学计算 《计算方法》 《数值逼近》 《数值代数》或者《矩阵计算》 《微分方程数值解法》 《非线性数值解法》 成绩的评定方法 期末考试：70% 闭卷 平时成绩：30%=10%平时 +20%上机 上机： 周四 7：00-10:30 理科楼二楼西侧机房 考试： 暂定十月底十一月初 第一讲 引论 引论 数值计算的研究内容 误差与有效数字 数值算法设计要点 1、数值计算的研究内容 研究过程 （理论上有解，而无求解公式或计算量过大难以用手工实现的数学问题） 地球外部大气流动模型 飞机外形优化设计问题 研究内容 线性方程组的直接法 插值方法 数值积分 常微分方程的差分方法 方程求根的迭代法 线性方程组的迭代法 数值分析课程的期望 掌握各种解决数学问题的数值方法 对近似解进行评估 在计算机上实现求解 仿真模拟 2、误差与有效数字 误差无处不在防不胜防 例 (病态问题)（系数保留小数点后4位,再精确计算） x1=x2=x3=1 x1=1.1788, x2= ? 0.0006, x3= ?1.9658. 病态问题:失之毫厘谬之千里 d. 有效数字 3.数值算法设计要点 作业 P8：1-5 P9：1-9（2）（3） P10：1-11 P10：上机1-14 算法设计的几点注意事项 1. 避免相近二数相减 例：a1 = 0.12345，a2 = 0.12346，各有5位有效数字。 而 a2 ? a1 = 0.00001，只剩下1位有效数字。 ? 几种经验性避免方法： 当 | x | 1 时： 2. 避免小分母 : 分母小会造成舍入误差增大 3. 避免大数吃小数 例：用单精度计算 的根。 精确解为 x 3.81574 y 0.0001 = = 38157.4 x 3.81574 y+?y 0.0001+0.00001 = = 34688.5 ? 算法1：利用求根公式 在计算机内，109存为0.1?1010，1存为0.1?101。做加法时，两加数的指数先向大指数对齐，再将浮点部分相加。即1 的指数部分须变为1010，则：1 = 0.0000000001 ? 1010，取单精度时就成为： 109+1=01010+0?1010=0?1010 大数吃小数 算法2：先解出 再利用 注：求和时从小到大相加，可使和的误差减小。 例：按从小到大、以及从大到小的顺序分别计算 1 + 2 + 3 + … + 40 + 109 4. 先化简再计算，减少步骤，避免误差积累。 一般来说，计算机处理下列运算的速度为 例：多项式求值：给定x, 求n次多项式的值 方法1：直接求和法； 方法2：利用前一次求幂的结果进行计算； 方法3：秦九韶方法； 5. 选用稳定的算法。 例：计算 ?公式一： 此公式成立, 因为 记为 则初始误差 ? ?? ? ! ! ! 考察第 n 步的误差 (unstable algorithm), 我们有责任改变。 这种算法是不稳定的算法 迅速积累, 可见初始的小扰动 误差递增. 可取 当 N ? ?时， ? 公式 方法：先估计一个IN ,再反推要求的 In ( n N )。 避免误差积累的方法 取 * * 实际问题 数学模型 数值分析理论 程序设计 上机计算 重点内容 研究