计算方法课件-第一章 绪论.pptxVIP

计算方法讲义;第一章 绪 论;§1.1 引 言; 鉴于实际问题的复杂性，通常将其具体地分解为一系列子问题进行研究，本课程主要涉及如下几个方面问题的求解算法： 非线性方程的近似求解方法； 线性代数方程组的求解方法； 函数的插值近似和数据的拟合近似； 积分和微分的近似计算方法； 常微分方程初值问题的数值解法； 矩阵特征值与特征向量的近似计算方法；等等;计算机解决科学计算问题的主要过程 实际问题——〉数学模型——〉数值计算方法——〉程序设计——〉上机运行求出解。其中： 实际问题——〉数学模型：由实际问题应用科学知识和数学理论建立数学模型的过程，是应用数学的任务. 数值计算方法——〉程序设计——〉计算结果：根据数学模型提出求解的数值计算方法，直到编出程序上机算出解，是计算数学的任务。;算法分类 分类方法1：若算法包含有一个进程则称其为串行算法，否则为并行算法。 分类方法2：从算法执行所花费的时间角度来讲，若算术运算占绝大多数时间则称其为数值型算法，否则为非数值型算法。 本课程介绍数值型串行算法。（其它类型算法参阅数据结构、并行算法等课程。）;算法的可靠性 可靠性：可靠性包括原问题的适定性算法的收敛性、稳定性等。 适定性：是指解存在、惟一，且解对原始数据具有连续依赖． 收敛性：是研究当允许计算步越来越多时，是否能够得到越来越可靠的结果，也就是研究截断误差是否能够趋于零． 对于给定的算法，稳定性分析是指随着计算过程的逐步向前推进，研究观测误差、舍入误差对计算结果的影响是否很大。 ;本课程的学习方法;§1.2 误差知识;一、误差来源及其分类;3）截断误差（方法误差） 数值方法精确解与待求解模型的理论分析解之间的差异。 这是由于我们需要将无穷过程截断为有限过程，而使得算法必须在有限步内执行结束而导致的。;4）舍入误差 在实现数值方法的过程中，由于计算机表示浮点数采用的是有限字长，因而仅能够区分有限个信息，准确表示某些数，不能准确表示所有实数，这样在计算机中表示的原始输入???据、中间计算数据、以及最终输出结果必然产生误差，称此类误差为舍入误差。 例如：利用计算机计算e的近似值en时，实际上得不到en的精确值，只能得到en的近似e*；这样e*作为e的近似包含有舍入误差和截断误差两部分： ;二、误差的度量;1. 绝对误差;2.相对误差;3.有效数字;有效数：当x* 准确到末位，即n=p，则称x*为有效数。 举例：x=π, x1*=3.141, x2*=3.142;Remark1: 有效数的误差限是末位数单位的一半，可见有效数本身就体现了误差界。 Remark2: 对真值进行四舍五入得到有效数。 Remark3:准确数字有无穷多位有效数字。 Remark4: 从实验仪器所读的近似数（最后一为是估计位）不是有效数，估计最后一位是为了确保对最后一位进行四舍五入得到有效数。 例 从最小刻度为厘米的标尺读得的数据123.4cm是为了得到有效数123.cm,读得数据156.7cm是为了得到有效数157.cm。;4.误差度量间的联系;定理证明;Remark 1、该定理实质上给出了一种求相对误差限的方法。 2、仅从 并不能保证x*一定具有n位有效数字。如 设其近似值a=0.484，其相对误差为： 我们并不能由此断定a有两位有效数字，因为;例题;三、误差的传播;初值误差传播：假设每一步都是准确计算，即不考虑截断误差和由运算进一步引入的舍入误差，仅介绍初始数据的误差传播规律。 研究方法： 泰勒（Taylor）方法 n元函数;复习泰勒公式;泰勒公式分析初值误差传播;相对误差：; 对于一元函数，有如下初值误差传播近似计算公式：;二元函数算术运算误差传播规律;尽量避免相近的数相减 例 x=52.127 x*=52.129 四位有效数字 y=52.123 y*=52.121 四位有效数字 A=x-y=0.004 A*=x*-y*=0.008 零位有效数字 结论：避免相近数相减;一些避免相近数相减示例 当|x|1时;两种算法的相对误差图比较 ;尽可能避免绝对值很小的数做分母，防止出现溢出。;算例 用不同位数的浮点数系统求解如下线性方程组 ;算例; ;方法1计算结果;方法一结果分析;方法2计算结果;总之, 除了算法的正确性之外, 在算法设计中至少还应:;;本章