数值分析例题I.pptVIP

*/18 数值分析典型例题 I 前四章内容小结 典型问题分析 割线法收敛速度 课外练习题 ? ? ? ? 具有n 位有效数字,则绝对误差满足 1. 有效数位数概念:设 x* 是准确值, x 是 x*的数值逼近. 如果 2. 二分法收敛定理: 设x*是 f(x)=0在[a, b]内的唯一根,且 f(a)·f(b)0,则二分法迭代数列 满足: | xn – x*|≤ (b – a)/ 2n+1 4. Newton迭代格式: 5. 割线法迭代格式: (n = 0, 1, 2 , ·····) 3.不动点迭代收敛定理: 若 , 满足 (1) ; (2) 则 在 [a, b] 有唯一的不动点 x*. 任取 x0∈[a, b], 迭代格式 收敛, 且 7. 收敛阶定理 设 x* 是 的不动点, 且 而 则 p阶收敛 设 , 若存在 a0 , r0 使得 则称数列{xn} r 阶收敛. 6.迭代法收敛速度 8.高斯消元定理: 约化主元 (k=1,···,n) ? 矩阵A的各阶顺序主子式不为零. 9.向量范数(常用三种): 10.矩阵算子范数(常用三种): 11. 矩阵条件数: Cond(A) = ||A ||·||A-1|| 12. 方程组迭代收敛定理:设x*是 Ax=b 的解 若||B||1,则对迭代格式 x(k+1) = B x(k) + f 有 例题结论: 1. 若 Ax = b 的系数矩阵 A 是严格对角占优矩阵,则Jacobi迭代和Seidel迭代均收敛; 2. 若 A 是实对称正定矩阵,则Gauss-Seidel迭法收敛. 13. 迭代法收敛定理: X(k+1) = B X(k) + f 收敛 ? 谱半径ρ(B) 1 14. SOR迭代收敛定理: 若 A 是对称正定矩阵, 则当 0 2 时, SOR 迭代法解方程组 Ax = b 收敛 15. 初等变分原理 设 为实对称正定矩阵, , 则 ? Ax = b . 应用:若矩阵 G 的列向量线性无关, 则 GTG 为实对称正定矩阵, 对于 G x = f , 有 GTGx = GTf ? , p1, p2∈Rn , A 是 n 阶对称正定矩阵 ? (Ap1, p2)=0 p1, p2 共轭 16.向量共轭 17. 共轭梯度法误差估计 设 x* 是方程组 Ax = b 的解，A的特征值为： ?1≥?2 ≥ …≥?n＞ 0 共轭梯度法迭代向量xk , 则有 其中，范数为 例1.计算球体 V 允许其相对误差限为 1%,问测量球半径R 的相对误差限最大为多少? 解：由球体计算公式分析误差传播规律 当球体V 的相对误差限为 1% 时，测量球半径R的相对误差限最大为0.33%。 ? ? 相对误差传播规律 思考: 计算球冠体积 V 允许其相对误差为1%,问应该对R, h 如何限制? 例2. 利用级数 计算无理数? 的近似值。由于交错级数的部分和数列 Sn 在其极限值上下摆动, 为了得到三位有效数近似值，应取多少项求和。 解: 由部分和 只需 ? n 1000时, Sn有三位有效数 例3 构造3阶下三角矩阵求逆算法 解: 设逆矩阵 对比右端左下角 求逆算法: 例4:设 求A-1的数值逼近 证 由迭代格式得 用牛顿迭代格式 设 计算. 证明: . 取X0 =I, 割线法收敛速度:设割线法收敛, 分析其收敛阶 解: 设 x* 是方程的 f(x)=0 解.割线法迭代格式 利用带余项的牛顿插值公式 其中 ? ? ? 设迭代法 p 阶收敛, 则 ? 割线法收敛速度为: 近似为 1.618 阶