11 数值算法概论.pptVIP

“计算方法”是计算数学的一个主要部分。而计算数学是数学科学的一个重要分支，它研究用计算机求解数学问题的数值计算方法及其软件实现. 数值计算已经成为计算机处理实际问题的一种关键手段，它使各科学领域从定性分析阶段走向定量分析阶段，从粗糙走向精密。 范数等价的定义: 例5 求矩阵A的各种常用范数 解(续): 由于 特征方程为 例5 求矩阵A的各种常用范数 解(续): 由于 特征方程为 容易计算 对矩阵元素的 变化比较敏感 较少使用 使用最广泛 性质较好 使用最广泛 计算较复杂 定义2.4 --------(9) 显然 Re Im ? ? ? ? ? ? ? ? 设A为n阶方阵,则对任意算子范数 || · || 有 证明： 由算子范数的相容性，得到 将任意一个特征根 ? 所对应的特征向量 代入 定理3 若A对称，则有 证明： 若? 是 A 的一个特征根，则?2 必是 A2 的特征根。 又对称矩阵的特征根为实数，即 ?2(A) 为非负实数， 故得证。 对某个 A 的特征根? 成立 定理4 定理. 证明:略 返回 现 实 世 界 研究 对象 测量 数据 建立数学模型 构成数值算法 数值运算的执行 测量 误差 模型误差 方法误差 舍入误差 结果 1. 3. 误差的种类及来源 总体设计(含数学模型的建立与模型细化等) 详细设计(主要是算法设计),包括： (1) 连续系统的离散化; (2) 离散型方程的数值求解. 实验验证 以计算机为工具解决实际问题需经历三个过程： (1) 模型误差 在建立数学模型过程中, 要将复杂的现象抽象归结为数学模型, 往往要忽略一些次要因素的影响, 而对问题作一些简化, 因此和实际问题有一定的区别。 (2) 观测误差 在建模和具体运算过程中所用的数据往往是通过观察和测量得到的, 由于精度的限制, 这些数据一般是近似的, 即有误差。 (3) 截断误差 由于计算机只能完成有限次算术运算和逻辑运算,因此要将有些需用极限或无穷过程进行的运算有限化,对无穷过程进行截断,这就带来误差。 (4) 舍入误差 在数值计算过程中还会遇到无穷小数,因计算机受到机器字长的限制,它所能表示的数据只能有一定的有限位数,如按四舍五入规则取有限位数,由此引起的误差。 如: 若将前若干项的部分和作为函数值的近似公式, 由于以后各项都舍弃了,自然产生了误差。 Taylor展开 截断误差 数学模型一旦建立,进入具体计算时所考虑和分析的就是截断误差和舍入误差. 舍入误差 绝对误差 其中 x*为精确值，x为x*的近似值。 例如： 工程上常记为 的上限记为 , 称为绝对误差限。 相对误差 x 的相对误差上限 定义为 2. 误差与有效数字 有效数字 用科学计数法，记 其中 若 则称 有n 位有效数字，精确到 。 ， 的截取按四舍五入规则 例4 解: 已知近似值为: 精确值为: 求 的绝对误差、相对误差 绝对误差 相对误差 或 证明： 例5 问： 有几位有效数字？请证明你的结论。 有4 位有效数字， 精确到小数点后第 3 位。 返回 数值算法是从给定的已知量出发，经过有限次四则运算 及规定的运算顺序，最后求出未知量的数值解，这样构成的完整计算步骤称为算法。 数值算法有四个特点: (1)目的明确 算法必须有明确的目的,其条件 和结论均应有清楚的规定 (2)定义精确 对算法的每一步都必须有精确的定义 (3)算法可执行 算法中的每一步操作都是可执行的 (4)步骤有限 算法必须在有限步内能够完成解题 过程 §1.4 题型分析与小结 - 算法和计算量 计算量：一个算法所需的乘除运算总次数，单位是flop. 计算量是衡量算法好坏的一个重要标准。 例1 求 Ax=b， Det(A)≠0，A=(aij)20 ? 20的计算量。 结论：分析算法的效率，选择算法非常重要。 解：2. 使用Gauss消去法， n=20， N≈3060次= O(n3 /3)次. 解：1. 用Cramar法则求解，总的计算量 N = ( ( n+1)(n-1)n!+n) 次 当n=20, N≈9.7 ? 1020 次. 以一台10亿/秒的计算机需约3万年. 例2 计算n 次多项式 的值。 算法: 采用秦九韶算法(1247) (又称为Horner