科学计算方法分析.ppt

雅可比迭代法 (i = 1,2,…n; k=1,2,……) 取初始向量X(0)=[x1(0) x2(0) ··· xn(0)]T, 迭代计算 (i = 1,2,…,n) 雅可比 迭代法的矩阵表示 将方程组AX = b 的系数矩阵 A 分解 A = D – U – L AX = b = DX(k+1) = (U+L)X(k) + b X(k+1)=D-1(U+L)X(k)+D-1b 记BJ = D-1(U+L) X(k+1)=BJX(k)+fJ */46 科学计算的背景 非线性方程求根算法 线性方程组求解直接法 线性方程组求解迭代法 《科学计算方法》 ? ? ? ? 科学计算方法与计算机有机结合 构造出强有力的工作平台 数值分析——研究用计算机求解 1969年, Apollo 登月计划实现 1981年,Columbia号航天飞机发射成功 数学问题的方法(算法)和理论 方程组求解、方程求根、数据插值、 数据拟合、数值积分、微分方程求解 von Neumann 1994年, GPS完全投入使用 例1: 圆内接正多边形边长计算Pi方法 评价算法的主要指标:速度和精度 简单迭代算法: n L error 192 3.1414524 1.4e-004 384 3.1415576 3.5e-005 3.1415926 4.6e-010 例2. 通信卫星覆盖地球面积 数学模型 实际问题 获取数据 数值方法、程序 数据结果 将地球考虑成一 个球体, 设R为地 球半径,h为卫星 高度,D为覆盖面 在切痕平面上的 投影(积分区域) 假设某一数据的准确值为 x*,其近似值 为 x，则称 而称 为 x 的相对误差 误差的有关概念 e(x)= x - x* 为 x 的绝对误差 如果存在一个适当小的正数ε ，使得 则称ε为绝对误差限。 称εr为相对误差限。 如果存在一个适当小的正数εr ，使得 十进制浮点数表示 一台微机价格:￥3999.00, 浮点数表示:0.3999×104 地球半径: 6378137m, (6.378137e+006) 浮点数表示: 0.6378137×107 光速: 2+008 浮点数表示: 0.299792458×109 尾数部 阶码部 有效数字概念: 取? 的有限位数如下(? ≈3.1415926) 取 x1 = 3，误差限不超过0.5; 取 x2 = 3.14,误差限不超过0.005 ; 若近似值 x 的绝对误差限是某一位上的半个 单位，该位到 x 的第一位非零数字一共有 n 位，则称近似值 x 有 n 位有效数字. 取 x3 = 3.1416,误差限不超过0.00005 ; r d 例4.水中浮球问题 有一半径r =10 cm的球体,密度? =0.638.球体浸入水中后,浸入水中的深度d 是多少？ 根据阿基米德定律,物体排开水的质量就是水对物体的浮力。 整理得: d 3 – 3 r d 2 + 4 r 3? = 0 由? =0.638, r = 10.代入,得d 3 – 30 d 2 + 2552 = 0 令 f (x) = x 3 – 30 x 2 + 2552 ,函数图形如下所示 求解方程 f(x)=0,即 是求函数 f(x)的零 点. f(x) 的零点所 在区间为:[0, 20] roots([1 -30 0 2552]) ans = 26.3146 11.8615 -8.1761 第一步:对根进行隔离,找出隔根区间,或在隔根区间内确定一个解的近似值x0; 设f(x) = 0的根为 x*,通过迭代计算,产生序列: x0 ? x1 ? x2 ? ··· ? xn········· 用数值方法求非线性方程的根,分两步进行: 第二