计算方法0213.pptx

计算方法 2月13日 1 教师信息 主讲：鞠颖 E_mail： yju@ 电话： 办公地点： 海韵科研2号楼312 助教： 2 待 定 教材 数值分析 第5版 作者：李庆扬、王能超、易大义 出版社：清华大学出版社 3 课程服务器 作用：保存课件、学习资料，实验 网址： / 选课密码： CPTM 4 课程名称 英文 Computing Method Computation Method Computational Method …… 中文 数值计算方法 数值分析 应用数学计算方法 …… 5 学习方式 课堂讲授：周学时4（每周两次） 课外自学：作业＋编程练习 没有实验课时? 6 成绩组成 笔试——70％ 作业——10％ 编程实验——20% 布置4个实验 选做两个 每个10分 7 第1章 数值分析与科学计算引论 8 1.1 数值分析的对象、作用与特点 英文 Computing Method Computation Method Computational Method …… 中文 数值计算方法 数值分析 应用数学计算方法 …… 9 目标 数值分析和数学分析 数学分析：解析方式求解 数值分析：数值方式求解 为什么需要学习计算方法 工程问题求解，如果解析方式不行（通常不行），用数值方式 麻烦的问题交给计算机 10 用数值的方法能解决哪些问题？ 插值 逼近与拟合 求积分、微分 求解线性方程组 求解非线性方程 求矩阵特征值 …… 11 用数值方法求解的基本思想 目标：只要结果达到要求的精度范围内 方法： 分解成若干个容易解决的问题 用替代品来拼凑 用迭代法逼近 误差控制： 估计误差数量级 用修正量补偿 12 不过…… 虽然我们是采用数值的方法求解，但是依据的还是解析方法的理论，目标也是解析方法的定义 所以，请各位想办法回忆一下之前的高等数学）、线性代数等课程的内容? 13 1.2 数值计算的误差 误差的来源 模型误差 观测误差 截断误差（方法误差） 计算误差（舍入误差） 这里我们只关心截断误差和舍入误差 14 误差的定义方式 定义1　设x*是精确值x的近似值，称 　 e* = x* – x 为近似值x*的绝对误差． 取x = 0.1， y = 9999.9的近似值分别为x* = 0.2， y* = 10000，则显然绝对误差并不能很好得表示误差的程度 相对误差： 15 绝对误差限 若已知|e*|的一个上界，即存在? 0，使得| e*| ? ?，则称?为近似值x*的绝对误差限．通常?亦记为?(x*)。 16 有效数字 定义2　若近似值x*的绝对误差限是某一位的半个单位，该位到x*的第一位非零数字共有n位，就说x*具有n位有效数字。它可表示为 x* = ±10m×(a1+a2×10-1+…+an×10-(n-1)). 其中ai (i=1，…，n) 是0到9中的一个数字,a1≠0, m为整数,且|x-x*|≤1/2×10m-n+1. 简单地说，x*的表示中，由左往从右第一个非０数字开始的数字的个数n称为x*的有效数字的位数． 0注意：运算所得数字并不都是有效数字！！！ 17 有效数字 18 相对误差限 若已知|er*|的一个上界，即存在? 0，使得 | er*| ? ?，则称?为近似值x*的相对误差限．通常?亦记为?(x*)。 注意： 绝对误差是有量纲单位的；而相对误差是无量纲单位的． 通常精确值是不知道的，所以实际计算中相对误差的计算用公式为 19 有效数字的性质 在m相同的情况下，n越大则10m-n+1有效位数越多，绝对误差限越小。 有效位数越多，相对误差限越小 20 21 22 变量误差引起的误差估计 23 证明：按定义立得． 24 自变量误差对函数误差的影响 25 多元函数 （A=f(x1,…,xn)） 所以…… 26 自变量误差对函数误差的影响 函数误差示例 ? 27 分析 ? 28 第一种计算方法的误差限 ? 29 第二种计算方法的误差限 ? 30 是否还可以更多？ 31 ? 1.3 误差定性分析与避免误差危害 一步误差限可以估计，但…… 实际的数值计算方法需要多步运算 误差限一般要远大于误差 32 稳定性 定义3　对于给定的数值方法，假设初始数据有舍入误差?，并设运算过程中产生的误差仅由?引起．若由?引起的误差积累可以得到控制，则称该方法是稳定的，否则称该数值方法不稳定． 病态问题与稳定性的关系 33 病态问题与条件数 一个数值问题，若输入数据有微小扰动（即误差）引起输出数据（问题解）相对误差很大，该问题是一个病态问题． Cp称为计算函数值问题的条件数． 条件数很大的问题是病态问题． 一般情况下，条件数Cp?10的问题就认为是病态的. 34 避免误差危害