数值计算方法1(免费阅读).pptVIP

华长生制作 华长生制作 数值稳定算法 在计算过程中产生的舍入误差能被控制在一定的范围内，且对最后的结果影响不大的算法称为稳定算法。不是数值稳定的算法称为数值不稳定算法。 数值不稳定算法会导致计算结果失真， 对数值不稳定的算法常采用转化成相应的数值稳定的算法来处理 。 例:方程x2-(109+1)x+109=0其精确解为x1=109, x2=1 字长为8位的计算器求解 当被减数和减数相差很大时，大数的相对误差起主要作用 很小 约为1 两个相近数相减,易失有效位 两正数之差 C=x-y的相对误差是 因为x和y的前几位有效数字必然相同，相减之后有效数字位会大大减少，使有效数字严重损失。 例如：cos20=0.9994,1- cos20=0.0006 避免这种情况,可以使用转换公式;或者增加字长,维持一定有效位,保证精度 乘积运算 乘积运算的相对误差为各乘数的相对误差之和，其相对误差限等于各乘数相对误差限之和 求c=12.2*73.56的相对误差限和绝对误差限 商运算 商运算的相对误差限等于除数与被除数的相对误差限之和 例 求c=25.7/3.6的相对误差限和绝对误差限 注意 (1) 当分母很小时，|dc|可能很大 舍入误差 放大了106倍 (2) 当分母为两个相近数相减时,会因有效数字丧失而出现(1)的情况 这里分子的误差被扩大104倍 指数运算 x的p次幂的相对误差是x本身的相对误差的p倍 若令p=1/q，则可得x的q次根的相对误差是x本身相对误差的1/q倍 例 求c=(12.2)2的绝对误差限和相对误差限 例 设正方形面积s=12.34,其绝对误差限|Δs| ? 0.01,问边长a具有多大的相对误差限和多少位有效数字？ 有4位有效数字 有6位有效数字 四、有效数字 定义3. 有8位有效数字 只有5位有效数字 且 因此,可根据上述分析对有效数字有如下结果: 或 五、数值方法的稳定性与算法设计原则 例8. 计算定积分 解: * 第一章 引论 第一章 绪论 1.3 误差 1.1 数值计算的研究对象与特点 1.2 数值问题与数值方法 本章要点: 绝对误差(限)和相对误差(限) 有效数字位数及其与误差的关系 数值问题的性态与误差的关系 数值算法设计原则 1.1 计算机数值方法的研究对象与特点 以计算机为工具,求解各种数学模型,都要经历三个过程: 总体设计——模型的细化 详细设计——主要为算法设计 程序设计 计算机数值方法研究的是将数学模型化为数值问题， 并研究求解数值问题的数值方法进而设计数值算法 数值问题： 输入数据与输出数据之间关系 即： 输入与输出的都是数值的数学问题 如求解线性方程组 求解二次方程 是数值问题 一、数值问题 1.2 数值问题与数值算法 求解微分方程 不是数值问题 将其变成数值问题，即将其“离散化” “离散化”是将非数值问题的数学模型化为数值问题 的主要方法，这也是计算方法的任务之一 二、数值方法 数值方法: 是指解数值问题的在计算机上 可执行的系列计算公式 在计算机上可执行的公式 是指只含有加减乘除的公式 现在的计算机中几乎都含有关于开方的标准函数sqrt() 常见的在计算机上不能直接运行的计算有: 开方、极限、超越函数、微分、积分等等 要在计算机上实行上述运算需将其化为可执行的等价 或近似等价运算 应化为 如求根公式 应化为公式 研究数值方法的主要任务: 1.将计算机上不能执行的运算化为在计算机上可 执行的运算 2.针对所求解的数值问题研究在计算机上可执行 的且有效的计算公式 3.因为可能采用了近似等价运算,故要进行误差分析, 即数值问题的性态及数值方法的稳定性 本课程的重点就是对线性方程组、微积分、微分方程、 矩阵特征值及回归拟合等问题寻找行之有效的数值方法 三、数值算法 数值算法是指有步骤地完成解数值问题的过程. 数值算法有四个特点: 1.目的明确 算法必须有明确的目的,其条件和结论 均应有清楚的规定 2.定义精确 对算法的每一步都必须有精确的定义 3.可执行 算法中的每一步操作都是可执行的 4.步骤有限 算法必须在有限步内能够完成解题过程 例1. 给出等差数列1,2,3,…,10000的求和算法 解: 记数器置零 1.3 误差 一、误差的种类及来源 模型误差 描述误差 在建立数学模型过程中,要将复杂的现 象抽象归结为数学模型,往往要