[理学]计算方法课件.ppt

* 1.3.2 机器数系 前述数的浮点表示方法是当今计算机采用的表示方法。 把计算机中浮点数所组成的集合加上“机器零”，记为F，则F被以下4个参数所描述： 集合F称为机器数系。 机器数系F是一个离散的分布不均匀的有限集。 若计算的结果超出了集合F的范围，则称为溢出。 β 基数 t 字长 [ L, U ] 阶码范围 * 1.4 误差危害的防止1.4.1 使用数值稳定的计算公式 什么是数值稳定的计算公式？先看一个例子 例2 建立积分 的递推关系，并研究它的误差传递。 解：由于 * 递推公式 在计算I0时有舍入误差e0，设I0的近似值为 1.4.1 使用数值稳定的计算公式 * 得到 进而 由此看出，误差e0对第n步的影响是误差扩大 倍，当n较大时，误差将淹没真值。 因此，用 近似 是不正确的，也就是说用这种递推公式求解问题是不正确的。 1.4.1 使用数值稳定的计算公式 * 但是如果我们由 得到 则可以先求出 ，再由上式依次得到其它值 使用广义积分中值定理（p10）可得 取 1.4.1 使用数值稳定的计算公式 * 则问题的第二种递推公式可表示为 研究I20的舍入误差 的传递情况 进而 所以第二种方法的误差传递是逐步缩小的。 1.4.1 使用数值稳定的计算公式 * 结论： 把计算过程中舍入误差对计算结果影响不大的算法称为数值稳定的算法，影响严重的算法称为数值不稳定的算法。 如果第n+1步的误差en+1与第n步的误差en满足 则称此计算公式是绝对稳定的，否则称为不是绝对稳定的。 1.4.1 使用数值稳定的计算公式 * 1.4.2 尽量避免两相近数相减 根据两数差的相对误差估计式： 可以看出，当两个相近的数 和 相减时， 可能比 大得多，从而可能严重损失有效数字。 * 例1：a1 = 0.12345，a2 = 0.12346，各有5位有效数字。而 a2 ? a1 = 0.00001，只剩下1位有效数字。 例2：计算 解：记 ，它们的六位有效数字分别为 方法一： 1.4.2 尽量避免两相近数相减 * 1.4.2 尽量避免两相近数相减 方法二： 分析两种方法的精度： (一) 所以第一种方法只能断定其结果至少具有2位有效数字。 * (二) 1.4.2 尽量避免两相近数相减 * 1.4.2 尽量避免两相近数相减 所以按第二种方法结果至少具有6位有效数字。 * 1.4.3 尽量避免用绝对值很大的数作乘数 因为 因为 当乘数 或 的绝对值很大时， 可能很大。 当 接近0时， 也可能很大。 结论：应尽量避免用绝对值很大的数作乘数，或用接近0的数作除数，否则会严重影响结果的精度，减少有效数字位数。 计算方法 西安电子科技大学 * 课程主要内容及教学计划安排 绪论（第1章，讲授+上机： 1+1学时） 方程求根（第2章，讲授+上机：2+2学时） 线性方程组数值解法（第3章，4+2学时） 插值法（第4章，4+2学时） 曲线拟合（第5章， 1+1学时） 数值积分与数值微分（第6章， 4+2学时） 常微分方程数值解法（第7章， 4+2学时） 矩阵的特征值及特征向量的计算（第8章，不讲） 考试：（2学时，开卷）×70%，带计算器 平时成绩：（作业+上机+考勤+大作业）×30% 上机要求：实验报告（7份），语言：尽量用C 第1章 绪论 * 第1章 绪论 内容提要 1.1 计算方法的对象与特点 1.2 误差的来源及误差的基本概念 1.3 机器数系 1.4 误差危害的防止 * 1.1 计算方法的对象与特点 《计算方法》研究什么内容？有什么特点？ 1 计算方法是研究数学问题的数值解及其理论的一个数学分支。 2计算方法主要研究适合于在计算机上使用的数值计算方法及与此相关的理论 计算方法：适合于计算机运算，计算量小。 相关理论：收敛性、稳定性、误差分析等。 3计算方法除具有数学的抽象性和严格性外，还具有应用的广泛性和实际计算的技术性，是一门与计算机密切结合的实用性很强的课程。 * 1.1 计算方法的对象与特点 计算方法解决实际问题的一般步骤： 提出实际问题 建立数学模型 提出数值问题 设计可靠、高效的算法 程序设计、上机计算 计算结果的可视化 数值计算已经被公认为与理论分析、实验分析并列的科学研究三大基本手段之一。 * 1.1 计算方法的对象与特点 理解每个算法建立的数学背景。 熟悉相关公式，学会推导证明方法，提高理论分析和解决实际问题的能力。 重点掌握：问题、方法的思想、性能分析。 进行习题练习。 进行数值计算