第1章 - 数值算的基本概念.ppt

第一章 数值计算的基本概念 引言 误差与有效数字 算法的稳定性与病态问题 计算机计算的几个问题 算法设计的原则 * * * §1 引言 一、数值学科的研究内容和特点 解决现代工程技术问题的基本过程（如左图） 数值分析是研究适合于在计算机上使用的实际可行、理论可靠、计算复杂性好的数值计算方法的学科。 分析实际问题 数值计算方法 数学模型 机器求解 数值分析研究的核心：算法。 构造算法的基本手段：近似 研究算法的核心问题：近似对计算的影响—误差分析 特点 构造性： 递推性： 离散化： 近似替代： 计算离散点上的近似值；有可靠的理论分析；算法理论主要是连续系统的离散化数值求解。 方法的构造，解的存在唯一性的证明 复杂计算过程转化成简单的计算过程的多次重复（适合计算机计算） 在误差允许的范围内，无限次的计算能用有限次计算替代。 模拟仿真： 可通过计算机的仿真实验验证实际的工程计算 §2 误差与有效数字 一个物理量的真实值和我们算出的值往往不相等，其差称为误差。 一、误差定义 定义 分析实际问题 数值计算方法 数学模型 机器求解 二、误差种类与来源 观测误差 模型误差 截断误差 舍入误差 数据观测所引入的误差，如电压、天体运行轨道。 数学模型是实际问题的抽象和简化，只是对客观现象的一种近似。其间出现的误差。 近似计算代替精确求解，如此误差称为截断误差，本身固有，也称方法误差。例如圆周率，自然对数。 计算机数系的有限字长，必须进行四舍五入，称为舍入误差。 计算机在处理数据过程中存在计算误差。原因是机器数系所致。这一数系的特点是有限、离散、支离破碎；这和数学上常用的实数系无限、稠密、连续的特点完全不同。机器数的表示方法通常采用浮点数形式，即： 其中 a1≠0，且 a1, a2, …, an 都是整数0~9中的任一个数。 10m称为尾数，尾数的位数n是有限正整数；m称为阶数，阶数也是有界的数。所以，机器数中有最大的数，也有最小的数。用机器数表示实数时，很多情况下都带有误差 。 附：计算机数系的特点 二、误差基本概念 1. 绝对误差 绝对误差不是误差的绝对值，即 e(x) 可正可负。 通常 x 是未知的，故 e(x) 未知 ，但一般地 ? 已知。 设 x*为准确值 x 的近似值，记 称ε为x *的绝对误差限或误差限。 称 e 为近似值x*的绝对值或误差。 定义 若 例 说明 2. 相对误差 设某量的准确值为 x，x*是 x 的近似值，称绝对误差与准确值之比 称δ为x*的相对误差限。 为x*的相对误差。 定义 若 例：设 估计近似数x1,x2的绝对误差与相对误差。 解： 是x1好的近似， 不是x2好的近似。 近似数的相对误差是近似数精确度的基本度量，一个近似数x*的相对误差越小，则近似数越精确。 结论: er(x)是一个无量纲的数，且 x一般是未知的，所以er(x)难求。 考察量 较小时 通常将 作为x*的相对误差。 相对误差限δ是未知的，但可以确定。 说明 3. 有效数字 若近似值x*的误差限是某一位的半个单位，该位到第一位非零数字 共n位，称近似值有n位有效数字。 设 x 为准确值，x*为 x 的近似值且x*表示为 换言之，如果 称x*具有n位有效数字。 定义 一个十进制数近似值的有效数字，不受单位制的影响。 如g* = 9.81m/s2作为g 的近似值，与g* = 0.00981km/s2，均 为3位有效数字。 在有效数意义下，不同的有效数位数的近似值的近似精度是不同的。如10.4200的精度高于10.42： 前者的绝对误差不超过： 后者的绝对误差不超过： 有效数字与绝对误差有一定的关系。对于某量的近似值，如果有n位有效数字，当m一定时，n越大则相对误差越小。 结论 定理 (有效数字与相对误差的关系)： 设近似值x*表示为 若x*具有n位有效数字，则其相对误差限为 反之，若x*的相对误差限 则x*至少具有n位有效数字。 证明： 由x*的表达式可得： 又由定义可知 有效数位越多，相对误差越小 所以 反之，则有： 因此， x*具有n位有效数字。 为使 近似值的相对误差小于1%，问需要取多少位有效数字？ 解： 的近似值的首位数字a1 = 2, 于是由 可解得 因此，可取 n = 3 即 例： §3 算法的稳定性和病态