北航数值分析b第1章1.pptVIP

数值分析 孙玉泉 sunyq@buaa.edu.cn * * 例 3216… * 例 设有方程组 * * 将系数舍入成 2 位浮点数, 则方程为 * 数值分析或数值计算方法主要是研究如何运用计算机去获得数学问题的数值解的理论和方法. 对那些在经典数学中,用解析方法在理论上已作出解的存在,但要求出他的解析解又十分困难,甚至是不可能的这类数学问题,数值解法就显得不可缺少,同时有十分有效. 研究对象和特点： * 计算机解决科学计算问题时经历的几个过程 实际问题——〉数学模型——〉数值计算方法——〉程序设计——〉上机运行求出解 数值计算方法重点研究：求解的数值方法及与此有关的理论 包括：方法的收敛性，稳定性，误差分析，计算时间的最小,占用内存空间少. 数值逼近 内容 数值代数 微分方程数值解 * 误差来源与种类 1. 来源与分类 从实际问题中抽象出数学模型 —— 模型误差 通过测量得到模型中参数的值 —— 观测误差 求近似解 —— 方法误差 (截断误差 机器字长有限 —— 舍入误差 * 一、绝对误差与相对误差 1、绝对误差 绝对误差(简称误差) , —— 设某量的准确值为x， 是 x 的近似值， 绝对误差限. 例1 通常 x 是未知的，故 未知 ， 但一般地 ? 已知。 绝对误差不是误差的绝对值，即 可正可负。 * 2、相对误差 —— 设某量的准确值为 x， 是 x 的近似值， 称绝对误差与准确值之比 为 的相对误差。 称δ为 的相对误差限。 例2 解 但 是 的一个好的近似， 不是 的好的近似. * ? 通常将 作为 的相对误差。 ? 相对误差限δ是未知的，但可以确定 是一个无量纲的数 * 二、有效数字 1、定义 设 x 为准确值， 为 x 的近似值且 可表示为 若 的误差满足 其中 为 0, 1, … , 9 中的一个数字。 即 有n位有效数字， 2、通俗说法 则称该位到 的第一位非零数 字为 的有效数字， * 例 1位有效数字，即n =1 5位有效数字，即n=5 但 有效数字的位数不能仅考虑 例 * 定理1 * * 反之： 例：4位有效数字舍入运算： 1234+0.4+0.3+0.2+0.1 0.4+0.3+0.2+0.1+1234=1235 原则1：应避免出现“大数吃小数” 数值计算的基本原则： * 原则2：应避免出现“相近的数相减”， 避免出现“很小的数做分母”。 * 例 求下列方程的根 解: 用 8 位浮点数 (有效数字)计算 用 4 位浮数点(有效数字)计算 两接近数相减 损失了有效数字 仍用4 位浮点数计算 * * 略去高阶项： 条件数： * 相对误差 初始数据 的微小误差可引起结果 A 的很大误差 引起 A 的很大误差 病态问题、 坏条件问题 计算问题的条件数 * 例 设有方程组 将系数舍入成 2 位浮点数, 则方程为 * 例 解 3216… 原因 —— 误差的传播与积累 误差传播： * 解： (a) 作矩阵和向量的乘法： (b) 作向量的内积和加法： 乘除法: kn 乘除法: kn2 例 算法效率： * (b) 利用秦九韶算法： (a) 直接计算每一项再求和： 解： * 例3 计算 ln2 的近似值，要求误差小于10??. 解： (b) 用级数 来计算 用前 9 项（即取 m = 8）计算就能达到精度要求: (a) 用级数 来计算 * 数值计算一般