计算方法绪论.pptVIP

杨东武 ydw_1978@126.com 第一章 绪论 研究计算方法的必要性 误差的基本概念 选用和设计算法应注意的问题 第一章 绪论 研究计算方法的必要性 误差的基本概念 选用和设计算法应注意的问题 计算方法研究的主要内容： 如何把数学模型归结为数值问题,如积分 的求取方法； 如何制定较好的数值计算方法,如提高运算速度、节省存贮空间等； 如何估计一个给定算法的精度,分析误差在计算过程中的积累和传播； 如何构造精度更高的算法，并研究算法的收敛性、稳定性等 第一章 绪论 研究计算方法的必要性 误差的基本概念 选用和设计算法应注意的问题 1.误差来源（分类） 模型误差：从实际问题建立的数学模型往往都忽略了许多次要的因素,因此产生的误差称为模型误差. 观测误差：一般数学问题包含若干参数,它们是通过观测得到的,受观测方式、仪器精度以及外部观测条件等多种因素，不可能获得精确值，由此而来产生的误差称为观测误差。 截断误差：数学模型的精确解与数值计算方法的近似解之间的误差称为截断误差。如 ex =1+x+x2/2+x3/6+x4/24+O[x]5，若用级数展开式中的前3项近似替代ex ,则其截断误差为Rn= x3/6+x4/24+O[x]5。 舍入误差：在计算机上运算时受机器字长的限制，一般必须进行舍入，此时产生的误差称为舍入误差。 例 1.1 试求摆长为L的单摆运动周期. 2.误差的基本概念 1．绝对误差。设x为准确值， 为近似值。则称 　 E = x - x*为x*的绝对误差。如果 ,则称 为x*的绝对误差限。 3.有效数位与误差的关系 第一章 绪论 研究计算方法的必要性 误差的基本概念 选用和设计算法应注意的问题 选用算法应遵循的原则 当|a||b|时，尽量避免a+b 。例如，假设计算机只能存放10位尾数的十进制数，则 考核方式: 开卷 平时成绩30%，考试成绩70% * 机电工程学院 什么是计算方法？ 它能够做什么？ 学习和了解科学计算的桥梁 计算方法／数值分析： 利用计算机求解数学问题近似解的方法。 探讨1：积分 如何求解? 高等数学求解积分的牛顿——莱布尼兹公式： 设 ，F(x)是 在[a,b]上的一个原函数，则有 的原函数无法求出，怎么办？ 通过数值分析方法，可求得其近似解。 探讨2：高阶线性方程组如何求解? 线性代数中求解线性方程组的克莱姆法则： 如果线性方程组 A x = b 的系数行列式不等于零，那末，方程组有唯一解 其中，D为方程组的系数行列式，Dj是把系数行列式中第j列的元素用方程组右端的自由项b代替后所得到的n阶行列式。 求解20阶线性方程组，用克莱姆法则要用 次乘法运算,用每秒1亿次的计算机计算,大约需算30多万年;但用消元法只须3000次乘法运算，只要几秒钟。 探讨3：在计算机上根据数学公式编程是否就能得到正确结果? 研究例子:求解线性方程组 其准确解为x1=x2=x3=1 如把方程组的系数舍入 成两位有效数字 它的解为x1 =-6.222... x2=38.25… x3=-33.65... 计算： 由高等数学知识可知： 1、大致趋势： 越大， 越小 2、递推关系： 计算 可使用以下两个公式进行计算： 公式一： 公式二： (N为一个较大数) 计算结果： 公式二 公式一 ? ?? ? ! ! ! 怎么回事?! 计算结果真好！ 〈已知结论〉大致趋势： 越大， 越小 结果分析： 公式一 记为 则初始误差 考察第n步的误差 可见初始误差 在计算过程中迅速积累，误差呈递增走势。 结论： 公式一： 算法不稳定。 结果分析： 公式二 以此类推，对 n N 有： 误差逐步递减, 这样的算法称为稳定的算法。 在我们今后的讨论中，误差将不可回避， 算法的稳定性会是一个非常重要的话题。 问题： 在以上计算过程中存在的误差都包括哪些呢？ 例如：在10位十进制数限制下： 例如：在10位十进制数限制下： 2．相对误差。称 　　为 的相对误差。 实用中，常用 表示 的相对误差。