数值第一章误差探讨.ppt

In I0 I1 I2 I3 I4 I5 I6 I7 I8 I9 算法1 0.6321 0.3679 0.2642 0.2074 0.1704 0.1480 0.1120 0.2160 －0.7280 7.5520 算法2 0.6321 0.3679 0.2642 0.2073 0.1709 0.1455 0.1268 0.1121 0.1035 0.0684 真值 0.6321 0.3679 0.2642 0.2073 0.1709 0.1455 0.1268 0.1124 0.1009 0.0916 ? 对任何 n都应有In0, 但算法1的计算结果显示I80, 可见, 虽然I0的近似误差不超过0.5×10－4, 但随着计算步数的增加, 误差明显增大. 这说明算法1给出的递推公式是数值不稳定的. ? 而对于算法2, 虽然初始给出的 I9没有一位有效数字, 但算至 I6已有4位有效数字. 这说明算法2中误差随着计算过程的深入是逐步递减的，因而是数值稳定的. * 数 值 分 析 林甲富 linjiafu@ 教材 丁丽娟，程杞元，《数值计算方法》，高等教育出版社，2011年. 最后成绩=实验作业成绩(20%)+考试成绩(80%) 实验作业：下列1和2选择一个,希望选2 1.课堂布置的课本上的数值实验题 2.结合所学专业自选题 (1)叙述实际问题 (2)建立数学模型(解常微分方程组，数据拟合等) (3)设计计算方法 (4)程序(matlab) (5)计算结果及分析 3.交打印文件(截止至考试前一日交) 4.注意完全重复的实验作业没有实验作业成绩 数值分析是做什么用的？ 数值 分析 输入复杂问题或运算 ? ? ? ? 计算机 近似解 ? 研究对象 那些在理论上有解而又无法手工计算的数学问题 ? 例 解300阶的线性方程组 求6阶矩阵的全部特征值 主要内容 ? 数值代数 近似求解线性方程组 (直接解法, 迭代解法) 矩阵特征值的计算 ? 数值逼近： 插值法，函数逼近 ? 数值微分与数值积分 ? 微分方程近似求解：常微分方程数值解法 ? 非线性方程求解 第一章 误差 §2 误差的基本概念 §3 数值计算中误差的传播 §4 数值计算中应注意的问题 ? §2 误差的基本概念 ? 误差按来源可分为： ? 模型误差 ? 观测误差 ? 截断误差 ? 舍入误差 误差：精确解与近似解之间的差 ? 模型误差 数学模型通常是由实际问题抽象得到的，一般带有误差，这种误差称为模型误差. ? 观测误差 数学模型中包含的一些参数通常是通过观测和实验得到的，难免带有误差，这种误差称为观测误差. ? 截断误差 求解数学模型所用的数值方法通常是一种近似方法，这种因方法产生的误差称为截断误差或方法误差. 实际计算时只能截取有限项代数和计算，如取前5项有： 这里产生误差 (记作R5 )截断误差 例如，利用 ln(x+1) 的Taylor公式计算 ln2, ? 舍入误差 由于计算机只能对有限位数进行 原则保留有限位，这时产生的误差称为舍入误差。 等都要按舍入 运算，在运算中像 在数值分析中，均假定数学模型是准确的，因而不考虑模型误差和观测误差，只讨论截断误差和舍入误差对计算结果的影响. ? 设x* 是准确值x 的一个近似值，记 e=x?x* 称 e为近似值 x* 的绝对误差，简称误差. 绝对误差一般很难准确计算, 但可以估计上界. ? 绝对误差 e 0 不唯一，当然 e 越小越具有参考价值. 则称 ?为近似值 x* 的绝对误差限，简称误差限. ? 若? 满足 ◆ 绝对误差、相对误差和有效数字 例 用毫米刻度的米尺测量一长度 x, 如读出的长度是 x*=765 mm, 由于误差限是 0.5 mm, 故准确值 ? 精确值x , 近似值 x* 和误差限 ? 之间满足： 通常记为 ? 绝对误差有时并不能完全地反映近似值的好坏，如测量 100 m 和 10 m 两个长度，若它们的绝对误差都是 1 cm，显然前者的测量结果比后者的准确. ? 因此，决定一个量的近似值的精确度，除了要看绝对误差外，还必须考虑该量本身的大小. 称 er 为近似值 x* 的相对误差. ? 记 由于 x 未知，实际使用时总是将 x* 的相对误差取为 ? 相对误差 称为近似值x*的相对误差限. ? 例 设 x*=1.24是由精确值 x 经过四舍五入得到的近似值, 求x*的绝对误差限和相对误差限. 由已知可得: 所以 ? =0.005, 解 ? 一般地, 凡是由