福州大学数值计算课件.pptVIP

数值计算方法 第一章 绪论 课程的研究对象与特点 数值计算的误差 误差的基本理论 误差来源与分类 误差来源与分类 误差来源与分类 机器字长有限—舍入误差 误差来源与分类 小结:模型误差.观测误差不是数值分析讨论的内容,计算方法主要研究截断误差和舍入误差在计算过程中的传播和对计算结果的影响,以提高计算的精度. 数值计算的误差 误差的基本理论 误差与有效数字 绝对误差/* Absolute error */ 有效数字 * * 以计算机为工具,求解各种数学模型,都要经历三个过程: 总体设计——模型的细化 详细设计——主要为算法设计 程序设计 计算机数值方法研究的是将数学模型化为数值问题， 并研究求解数值问题的数值方法进而设计数值算法 数值 分析 输入复杂问题或运算 ? ? ? ? 计算机 近似解 1 用计算机进行实际问题的数值计算时,往往求得的是问题的近似解,都存在误差 在建立数学模型过程中,要将复杂的现象抽象归结为数学模型,往往要忽略一些次要因素的影响,而对问题作一些简化,因此和实际问题有一定的区别.—模型误差 在建模和具体运算过程中所用的数据往往是通过观察和测量得到的,由于精度的限制,这些数据一般是近似的,即有观测误差 如: 若将前若干项的部分和作为函数值的近似公式, 由于以后各项都舍弃了,自然产生了误差。 Taylor展开 由于计算机的字长有限,只能对有限位数进行运算,超过的位数按一定规则舍入,产生“舍入误差”. 2 误差是不可避免的,既要允许误差,又要 控制误差.要重视误差分析,分析误差的来源,误差的传播及对误差作出估计 2. 传播与积累 /* Spread Accumulation */ 例：蝴蝶效应 —— 纽约的一只蝴蝶翅膀一拍，风和日丽的北京就刮起台风来了？！ NY BJ 以上是一个病态问题 /* ill-posed problem*/ 关于本身是病态的问题，我们还是留给数学家去头痛吧！ 蝴蝶效应是气象学家洛伦兹1963年提出来的。其大意为：一只南美洲亚马孙河流域热带雨林中的蝴蝶，偶尔扇动几下翅膀，可能在两周后引起美国德克萨斯引起一场龙卷风。其原因在于：蝴蝶翅膀的运动，导致其身边的空气系统发生变化，并引起微弱气流的产生，而微弱气流的产生又会引起它四周空气或其他系统产生相应的变化，由此引起连锁反映，最终导致其他系统的极大变化。此效应说明，事物发展的结果，对初始条件具有极为敏感的依赖性，初始条件的极小偏差，将会引起结果的极大差异。 例：计算 ? 公式一： 注意此公式精确成立 记为 则初始误差 ? ?? ? ! ! ! What happened?! 考察第n步的误差 我们有责任改变。 造成这种情况的是不稳定的算法 /* unstable algorithm */ 迅速积累，误差呈递增走势。 可见初始的小扰动 ? 公式二： 注意此公式与公式一 在理论上等价。 方法：先估计一个IN ,再反推要求的In ( n N )。 可取 取 We just got lucky? 考察反推一步的误差： 以此类推，对 n N 有： 误差逐步递减, 这样的算法称为稳定的算法 /* stable algorithm */ 在我们今后的讨论中，误差将不可回避， 算法的稳定性会是一个非常重要的话题。 定义1. Hey isn’t it simple? Oh yeah? Then tell me the absolute error of Oops! 例 用一把有毫米的刻度的米尺,来测量桌子 的长度,读出的长度x* =1235mm 这是桌子实际长度x一个近似值,由米尺的精度知, 这个近似值的误差不会超过0.5mm(即绝对误差限为1/2mm),则 1234.5≤x≤1235.5 即x∈[1234.5,1235.5] 或x=1235±0.5mm 绝对误差的某个上界 绝对误差限或误差限, 或 且 误差限的大小还不能完全表示近似值的好坏. 注：e* 理论上讲是唯一确定的，可能取正，也可能取负。 e* 0 不唯一，当然 e* 越小越具有参考价值。 显然 I can tell that this part’s diameter is 20cm?1cm. I can tell that distance between two planets is 1 million light year ±1 light year. Of course mine is more accurate ! T