数值分析第一章引论.ppt

四、数值计算中的算法设计技术 4. 加权平均的松弛法 采用加权平均的思想提高精度，减少计 算量 刘徽‘割圆术’求得π=3.1416，需割到3072边形，但事实上他利用了松弛技术 四、数值计算中的算法设计技术 4. 加权平均的松弛法 采用加权平均的思想提高精度，减少计 算量 辛普森求积公式 数值方法 数学系 教师信息 授课教师信息 数学系 刘娟 办公室：南海楼317 办公电话邮箱：jnuljstudent@163.com 教材信息 教材 《数值分析》（第五版），李庆扬等 著，清华大学出版社，2012. 参考书 适合自己的才是最好的！需要时去图书馆找《数值分析》、《计算方法》、《数值方法》等包含相关名称的书籍作参考。 学分和成绩评定 学分及授课方式： 3+1学分，课堂讲授48学时，上机实验32学时 成绩评定 期末成绩：60% 实验成绩：20% 平时及期中成绩：20% 第一章 引论 一、数值分析简介 数值分析 研究各种数学问题的数值计算方法。 这些数值计算方法需要：面向计算机、有可靠的理论分析、有好的计算复杂性。 包括“数值代数、数值逼近和微分方程数值解法”三部分内容。 二、数值计算的误差 误差来源 模型误差 截断误差（数值方法本身产生的误差） 如用 近似 所造成的误差 舍入误差（原始数据误差和计算机舍入误差） 如用3.14近似pi，或用计算机近似0.1时所产生的误差 二、数值计算的误差 2. 误差与有效数字 绝对误差：e*=x*-x 相对误差：er*=(x*-x ) ∕ x 有效数字 二、数值计算的误差 2. 误差与有效数字 绝对误差：e*=x*-x 误差限：由于x通常无法确定，只能估计其绝对误差值不超过某整数ε*，则称ε*为绝对误差限。 相对误差：er*=(x*-x ) ∕ x 有效数字 二、数值计算的误差 2. 误差与有效数字 绝对误差：e*=x*-x 相对误差：er*=(x*-x ) ∕ x 有效数字 二、数值计算的误差 2. 误差与有效数字 绝对误差：e*=x*-x 相对误差：er*=(x*-x ) ∕ x 相对误差限：由于x通常无法确定，只能估计其相对误差近似值|x*-x | ∕ |x*|不超过某整数εr*，则称εr* 为相对误差限。 有效数字 二、数值计算的误差 2. 误差与有效数字 绝对误差：e*=x*-x 相对误差：er*=(x*-x ) ∕ x 有效数字 二、数值计算的误差 2. 误差与有效数字 绝对误差：e*=x*-x 相对误差：er*=(x*-x ) ∕ x 有效数字 若近似值x*的误差限不超过某一位的半个单位，该位到x*的第一位非零数字共有n位，就说x*有n位有效数字。x*可表示为 其中， a1, a2, …, an 都是 0 ～ 9 中的任一整数，但 a1≠0。 其绝对误差限满足： 二、数值计算的误差 2. 误差与有效数字 绝对误差：e*=x*-x 相对误差：er*=(x*-x ) ∕ x 有效数字 如取 作为 的近似值， 取 ， 就有3位有效数字， 就有5位有效数字. 取 ， 有几位有效数字？ 1 二、数值计算的误差 2. 误差与有效数字 绝对误差：e*=x*-x 相对误差：er*=(x*-x ) ∕ x 有效数字 反之， 则 至少具有 位有效数字. 定理 若 具有 位有效数字， 则其相对误差限为 若 的相对误差限 ， 二、数值计算的误差 3. 数值运算的误差估计 三、误差分析与避免误差危害 误差分析：通俗来说，由于初始 数据误差以及算法中每一次运算 都受舍入误差影响，需要分析它们对结果的总的影响 三、误差分析与避免误差危害 误差定量分析（了解） 向前误差分析 向后误差分析 概率分析法 三、误差分析与避免误差危害 误差定量分析（了解） 向前误差分析 如fl(AB)=AB+ε 向后误差分析 如fl(AB)=(A+ε)(B+δ) 概率分析法 三、误差分析与避免误差危害 2. 误差定性分析 算法的数值稳定性 病态问题与条件数 三、误差分析与避免误差危害 2. 误差定性分析 算法的数值稳定性 例如关于In=e-1∫01xnexdx (n=0,1,...)的两种算法 一个算法如果输入数据有误差，而在计算过程中舍入误差不增长，则称算法数值稳定，否则称此算法不稳定 病态问题与条件数