(西安石油大学现代数值计算方法第章.pptVIP

课程名称： 数值分析 或计算方法 (考试时要带上计算器(有函数的)) 第一章 数值计算中的误差分析（绪论） 内容提要 §1.1 计算方法的任务与特点 §1.2 误差与误差估计(误差知识) §1.3 选用算法时应遵循的原则 §1.1计算方法的任务与特点 科学与工程计算过程： 实际问题→数学模型→数值问题→算法→程序→调试→结果 计算方法的特点 严密的科学性、实验的技术性、高度的抽象性、应用的广泛性、提供算法、算法分析、兼顾计算机的特点：有效数字（精度）、运算量、存储量等。 第一章 绪论 第二章 线性方程组求解 第三章 非线性方程求解 第四章 矩阵特征值问题（不讲） 第五章 函数的插值 第六章 曲线拟合 第七章 数值积分和数值微分 第八章 常微分方程数值解法 §1.2 误差知识（误差与数值计算中的误差估计） 内容提要： 误差的来源及其分类 误差的度量（误差与有效数字） 数值计算的误差估计 一、误差来源及其分类 1)模型误差（描述误差） 反映实际问题有关量之间的计算公式（数学模型）通常是近似的。 3）截断误差（方法误差） 数值方法精确解与待求解模型的理论分析解之间的差异。 这是由于我们需要将无穷过程截断为有限过程，而使得算法必须在有限步内执行结束而导致的。 4）舍入误差 以四舍五入为例（也可以五舍六入等） 最多舍去或添加最后一位的半个单位。 注意：与截断误差不同！ 二、误差的度量 绝对误差 相对误差 有效数字 度量间的关系 1. 绝对误差 绝对误差定义：准确值 x 减近似值 x* 2.相对误差 绝对误差限虽然能够刻划对同一真值不同近似的好坏，但它不能刻划对不同真值近似程度的好坏 。 有效数：当x*为四舍五入得到的近似数，则称x*为有效数。有效数的绝对误差限、相对误差限，有效数字位数举例： 例①若有效数x*=1.02,则绝对误差限： η(x*)=0.005相对误差限： δ(x*)= η(x*)/| x*|= 0.0049， x*具有3位有效数字. ③若有效数x*= 25×102，则 η(x*)=50 ，δ(x*)= η(x*)/| x*|= 0.02，x*具有2位有效数字. 3.有效数字（教材第6页中间）： 若近似数x*的绝对误差限是（不超过）某一位的半个单位，则称其精确到这一位，且从该位到x*的第一位非零数字共有n位，则称近似数x*具有n位有效数字. （黑板举例） 举例：x=π=3.1415926…, 近似数 x1*=3.14102, 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有两位有效数字，因为 例题 三、数值运算的误差估计 由于自变量的误差，引起函数值的误差。也称为误差传播。这种函数值的误差往往很难精确描述，只能大体估计。 例如一元函数，二元函数，等。 复习泰勒公式 泰勒公式分析初值误差传播 §1.3 选用算法应遵循的原则 方法2计算结果 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 简化计算步骤以减少运算次数。 例1 例2 秦九韶算法 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 例 尽量避免相近的数相减 例 x=52.127 x*=52.129 四位有效数字 y=52.123 y*=52.121 四位有效数字 A=x-y=0.004 A*=x*-y*=0.008 零