数值分析--颜庆津-北京航空航天大学出版社-2000精选.pptVIP

数值分析 Numerical Analysis 第 1 章 绪 论 研究内容 线性方程组的数值解； 矩阵特征值与特征向量计算； 非线性方程的数值解； 数值逼近； 数值积分； 常微、偏微方程的数值解； 研究方法 理论分析 算法分析 误差分析 收敛性分析 收敛速度讨论 §1.2 误差知识与算法知识 1.2.1 误差的来源与分类 模型误差 处理实际问题时，要建立数学模型，通常模型只是近似的。由此产生的数学模型解与实际问题的解 之间的误差叫模型误差。 例如 是实际问题的解，而若数学模型的解是 由此产生的误差叫作模型误差。 观测误差 数学模型中包含某些变量，如时间、长度、电压等，它们一般是通过观测来获得。由于观测得到的数据与实际数据之间有误差，这种误差叫观测误差。 舍入误差 由于计算机的字长有限，参加运算的数据及其运算结果在计算机中存放会产生误差。这种误差叫舍入误差 或 计算误差。 例如，在 16 位微机上计算，单精度实数存放仅有 7 位有效数字。在其上运算，会有 1 ? 3 ? 0.333 333 3, (1.000 002)2 ?1.000 004 ? 0; 后者的准确结果是 4 ?10?12。 1.2.2 误差传播与积累——算法的稳定性 例如，a ? 3.140 作为圆周率? 的一个近似值，它的绝对误差是 e ? ? ?3.14 又 |e| ? 0.002. 故， a ? 3.14 作为? 的一个近似值，它的一个绝对误差限是 ? ? 0.002。 下面定义的 相对误差 更全面地刻画了近似值与精确值之间差异的情况. 例1 用最小刻度为毫米的卡尺测量直杆甲和直杆乙，分别读出长度 a ?312 mm, b ? 24 mm. 问：? (a), ? (b), ?r(a), ?r(b) 各是多少？两直杆实际长度 x 和 y 在什么范围？ 解 由最小刻度为毫米知, ? (a) ? ? (b) ?0.5 mm, 例2 设 a ? ? 2.18 和 b ? 2.120 0 分别由准确值 x 和 y 经过四舍五入而得到的近似值. 问：? (a), ? (b), ?r(a), ?r(b) 各是多少？ 解 ? (a) ? 0.005，? (b) ?0.000 05， 有效数字 ( significant digits) ?四舍五入带来的绝对误差限 凡是由准确值 x 经四舍五入而得到近似值 a，其绝对误差限等于该近似值 末位 的半个单位。 例如， ?* ? 3.142 作为圆周率? ?3.14159265358 … 的一个近似值，它的绝对误差限是 误差限是 10?2 位的半个单位，从 10?2 位到 ?* 的第一位非零数字共有 3 位，故依定义知, 用 3.142 去近似 ? 时，具有 3 位有效数字。?* 精确 ? 到小数点后第 2 位。 定理1 把准确数 x 的近似值 a 写成 a ? ? 0. a1a2?ak ?10 m, a1 ?0, m?Z. 又绝对误差限 ? (a) ? 0.5 ?10 t, 则 a 的有效数字位数 n ? m ? t 。 证明 ? (a) ? 0.5 ?10t, 仅证情况 m ? 0 ，其它同理。因为, 例3 下列近似值的绝对误差都是 0.005, a ?1.38, b ? ? 0.0312, c ? 0.86 ?10?4, 问，各个近似数有几位有效数字？(提示,使用定理1) 解 因为， a ? 0.138 ? 101, ? (a) ? 0.5 ? 10?2, 知 a 有 3 位有效数字，精确到小数点后第 2 位. b ? 0.312 ? 10?1, ? (b) ? 0.5 ? 10?2, 知 b 有 1 位有效数字：3, 精确到小数点后第 2 位. 对于 c, m ? ?4, t ? ? 2, n ? m ? t ? ? 2, 故 c 无有效数字。 1.2.5 函数的误差估计 一元函数的误差 设 u = f (x) 有足够高阶的导数, a 是自变量 x 的近似值, 则 u* = f (a) 是函数值 u = f (x) 的近似值. 如果 ??(a) ? 0, 且当 a 与 x 很接近时, 有误差估计 e(u*) ? f (x) ? f (a) ? ??(?)(x ? a) ? ??(a)e(a), 因 |e(u*)|? |??(a)|?|e(a)| ? |??(a)|?? (a), 故 ? (u*) ? |??(a)|?? (a). 如果 ??(a) ???(a) ?… ??(