第三章数值计算方法.pptVIP

《计算机在冶金中的应用》 第3章 数值计算方法 2.绝对误差与相对误差 3.有效数字 2.算术运算的误差传播估计 3-4 非线性方程求解 用计算机求解非线性方程的步骤 4. 牛顿迭代法 5. 消元法 4. 迭代法求解举例 计算机在稀土工程中的应用 计算机在稀土工程中的应用 第三章 数值计算方法 计算方法（也称数值分析、数值求解）：研究求解数学模型的算法及相关理论，它是随着计算机技术的发展而发展起来的应用数学分支。 在工程实际和科学研究中，绝大多数数学模型很难求得其解析解，或根本就没有解析解，只能用计算机求出其数值解。从应用角度来看，这也就足够了。 在数值求解中，常使用如下技术: 1.迭代求解技术 2.连续问题离散化技术 3.离散数据连续化技术 评价算法优劣的标准： 1.速度: 速度涉及计算量、算法收敛速度等 2.精度: 涉及计算误差 3-1 误差与有效数字 数值分析，常常给人一种不严格、不准确的感觉。实际上，误差是不可避免的，近似是正常的，根本不存在绝对的严格和准确。问题是我们如何分析误差，控制误差，保证计算结果误差在允许范围内。 1.误差的来源 在用数模处理实际问题时，误差来源主要有以下四个方面： 模型误差、观测误差、截断误差、舍入误差 1 模型误差：模型的建立往往要忽略一些次要因素，所以模型的解与实际结果必然存在一定的误差，这类误差称为模型误差； 2 观测误差：数学模型常常要用到一些参量数据，如温度、初始浓度等，这些数据又往往是由观测获得，由观测数据带来的误差称为观测误差； 3 截断误差：许多函数常用无限级数表示，而用计算机进行计算时，只能取前有限项，由此引起的误差称为截断误差；比如函数: 4 舍入误差：由于计算机表示的实数有效位数有限，超出部分只能按四舍五入法处理，这部分误差称为舍入误差。尽管一次舍入误差极小，但当运算次数极多时，误差的积累也可能不可忽视。 设 x* 为某数据的准确值，x 为近似值 绝对误差： 若 称为绝对误差限 2 相对误差： 当准确值未知时，可用近似值代替： 称为相对误差限 设近似值 x 的绝对误差限不大于某一位上的半个单位，且从该位到的第一非零数字位共有 n 位，则近似值 x 有n 位有效数字。 设数据近似值为：0.0座机电话号码 若绝对误差限为：0.00006，则有效数字有3位，即0.0276; 若绝对误差限为：0.00004，则有效数字有4位，即0.02765; 0.00006 0.0005 0.5×10-3 0.00004 0.00005 0.5×10-4 相对误差与有效数字位数，在最保守的情况下有如下关系： |相对误差| 3-2 数值计算的误差估计 数值计算中误差的产生与传播的情况非常复杂，对误差估计也比较困难，本节主要介绍采用函数的泰勒 Taylor 级数展开来估计误差，这是一种常用方法。 1.函数运算的误差传播估计 同理，对多元函数 设两个近似数 x1 , x2 的误差限分别为 ，则其两个数的四则运算误差，可用前述多元函数运算误差估计公式计算： 3-3 数值计算中的一些基本原值 在数值计算中，为使误差在传播过程中不增大，应遵循一些基本原则: 1.避免绝对值小的数作除数 2.避免两个相近的数据相减 3.防止大数“吃”小数 4.尽量减少计算工作量 5.选用数值稳定性好的算法 1. 避免绝对值小的数作除数 由除法的误差估计公式可知，除数的绝对值越小，商的绝对误差越大。所以要尽量避免采用绝对值很小的数作除数。 2.避免两个相近的数据相减 当两个相近的数据相减时，差的有效数字位数大大减少，相对误差大大增加。 相对误差 为了避免两个相近的数据直接相减，常用恒等式将其变形后计算，如： 绝对误差 3.防止大数“吃”小数 一个绝对值很大的数与一个绝对值很小的数相加时，很容易发生大数“吃”小数的现象。这是因为计算机表示的实数有效数字位数有限所至。 当绝对值很小的数个数极多时，对计算结果的影不可忽视。在处理这类累加问题时，应按绝对值从小到大的顺序进行累加。 4.尽量减少计算工作量 在考虑算法时应注意尽量减小运算次数，计算机完成算法运算所花费的时间与算术运算量有关，主要取决于乘法、除法的运算次数，减少乘法、除法的运算次数，不仅缩短了解题时间，而且误差的积累相应也要小一些。 比如，多项式求和可用两种算法处理： 1 直接累加 2 秦九韶算法 共进行n次乘法运算 共进行n n+1 /2次乘法运算 5.选用数值稳定性好的算法 对同一问题，采用不同数值求解算法，对计算结果的误差影响也不一样。舍入误差对计算结果影响不大的算法称为数值稳定算法。实际上，分析一个算法在计算过程中是否稳定，就是考察其误差是否增大。 例：计算定积分 经推导，有下式成立： 算法1: 正向递推 由此可见，该递