数值分析第一章.pptVIP

可求得I9 ? 0.017，按改写后的公式可逐次求得 不妨设I9 ? I10，于是由 将公式 变为 I8 ? 0.019 I7 ? 0.021 I6 ? 0.024 I5 ? 0.028 I4 ? 0.034 I3 ? 0.043 I2 ? 0.058 I1 ? 0.088 I0 ? 0.182 稳定的算法 ！ 在我们今后的讨论中，误差将不可回避， 算法的稳定性会是一个非常重要的话题。 注：递推公式(1)的舍入误差以5的幂次增长进行传播，因此是数值不稳定的，而递推公式(2)的舍入误差在一定范围内以0.2的幂次进行传播，随着n的增大，误差逐步减少，因此该算法是数值稳定的。 因此，可以看出数值不稳定的算法是不能使用的，实际计算中对任何输入数据都是数值稳定的算法，称为无条件稳定。而对某些数据数值稳定，对其它数据数值不稳定的算法，称为条件稳定。 1.要避免两个相近的数相减 在数值计算中，两个相近的数作减法时 有效数字会损失。 例: 求 的值。当x = 1000，y 的准确值为0.01580 §6 数值计算中应该注意的一些原则 类似地 (2) 若将原式改写为 则 y = 0.01581 (1)直接相减 有3位有效数字！ 只有1位有效数字 2.尽量避免绝对值太小的数作分母 例： 如分母变为0.0011，也即分母只有0.0001的变化时 结果相差这么大! 3. 避免大数吃小数 精确解为 ? 算法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. 秦九韶方法(即Hornor算法）； 先计算x2, x3, …, xn, 再作线性组合，需做 2n-1次乘法和n次加法。 解法一：直接求和法 解法二：逐项求和法 按顺序依次计算每一项的值再求和，需做 n(n+1)/2次乘法和n次加法。 解法三：秦九韶算法（即Horner算法） 只需做n次乘法和n次加法。且可以递推实现。 计算机上使用的算法常采用递推化的形式， 递推化的基本思想是把一个复杂的计算过程 归结为简单过程的多次重复。这种重复在程序 上表现为循环。递推化的优点是简化结构和节 省计算量。 算法的递推性 例：用秦九韶方法求多项式p(x)在x=-2处的值 解： K a5-K vK 0 0.00833 0.00833 v0 = a5 1 0.04167 0.04 v1 = v0x+a4 2 0.16667 0.15867 v2 = v1x+a3 3 0.5 0.46827 v3 = v2x+a2 4 1 0.90635 v4 = v3x+a1 5 1 0.81873 v5 = v4x+a0 约翰·冯·诺依曼 （John von Neumann,1903-1957）美藉匈牙利人，1930年接受了普林斯顿大学客座教授的职位，西渡美国。 1931年成为该校终身教授。1933年成为新成立的普林斯顿高等研究院的终身研究员。1951年至1953 年任美国数学会主席。 冯·诺依曼是20世纪少有的数学科学通才，在许多领域都有重要的贡献，被西方人誉为“数学奇才、计算机之父”。冯·诺依曼对人类的最大贡献是对计算机科学、计算机技术和数值分析的开拓性工作。 并行计算 一、 电子计算机的并行计算分类 传统计算机一般采用Von Neumann结构，每一时刻 只有一个进程的算法，即串行算法。 并行计算机每一时刻具有2个以上的进程的算法称为 并行算法。并行机必须包含2台以上处理机，按指令流是 单个还是多个并行算法可分为两类： SIMD算法适用于单指令流多数据流系统； MIMD算法适用于多指令流多数据流系统。 按照进程之间是否同步可将并行算法分为: 同步算法:是指在k个进程的算