数值分析10迭代法的收敛性分析.pptVIP

宝山壁画 宝山壁画是引人注目的昂贵文物。此壁画发现于阿鲁科尔沁旗东沙布乡境内。1994年列为“全国十大考古新发现”之一。宝山壁画中最引人注目的是《杨贵妃教鹦鹉图》。该画高0.7米、宽2.3米，用于笔重彩绘制，最突出的表现了 晚唐风格。唐代擅长绘贵妇仕女的大师周昉绘制了《杨贵妃教鹦鹉图》，不仅享誉中原，而且还影响全国各地。发现于阿旗宝山古墓里的这幅画，就是契丹人聘请中原画家按照周氏风格绘制的， 技法深得周氏画风的真传。在唐人真迹稀如星风的今天，能够从中完整了解唐代人物画的杰出成就，堪称美术史研究的辛事。这幅壁画现今保存在阿鲁科尔沁旗博物馆，历经千年，恍如新绘，是该馆的镇馆之宝。 欢迎大家观看！ * 迭代法的收敛性分析 迭代矩阵谱半径 误差估计定理 对角占优矩阵 《数值分析》 10 平面点列: (x1 , y1),(x2 ,y2), ···, (xk , yk),······ Xk∈Rn : X1, X2, ···, Xk , ······· 利用向量范数等价性, 对任意范数 || · || 2/18 原方程: A X = b 记 ?(k) = X(k) – X* ( k = 0, 1, 2, 3, ······ ) 则有 ?(k+1) = B ?(k) ?(k) = B ?(k-1) ( k = 1, 2, 3, ······ ) 计算格式: X(k+1) = B X(k) + f X(k+1) – X*= B(X(k) – X*) 设方程组的精确解为 X*,则有 X* = B X* + f ? 3/18 (1) ?(k) = B ?(k-1)=B2 ?(k-2)=···=Bk ?(0) (2) 迭代格式 X(k+1) = B X(k) + f 收敛 ！ 4/18 证: 由?(k) = B ?(k-1),得 || ?(k)|| ≤ || B|| || ?(k-1)|| ( k = 1, 2, 3, ······ ) 所以 命题 若||B||1,则迭代法 X(k+1) =B X(k) +f 收敛 || ?(k)|| ≤ || B||k || ?(0)|| || B|| 1 5/18 矩阵B 的谱 设n阶方阵B 的n个特征值为: 则称集合 为B 的谱. 记为 ch B 矩阵B的谱半径 注1: 当B是对称矩阵时, ||B||2 = ?(B) 注2: 对 Rn×n 中的范数|| · ||,有 ?(B) ≤ || B || 特征值取模最大 6/18 定理4.1 迭代法 X(k+1) = B X(k) + f 收敛 ? 谱半径ρ(B) 1 证: 对任何 n 阶矩阵B都存在非奇矩阵P使 B = P –1 J P 其中, J 为B的 Jordan 标准型 其中, Ji 为Jordan块 7/18 其中,λi 是矩阵B的特征值, 由 B = P –1 J P B k = (P –1 J P) (P –1 J P) ··· (P –1 J P)= P –1 J k P 迭代法 x(k+1) = B x(k) + f 收敛 = (i = 1, 2,···, r) (i = 1, 2,···, r) 谱半径 ?(B) 1 8/18 注1: AX = b ? X = BX + f ? ( I – B )X = f ? X = ( I – B )-1 f 注2: 若 则 ( I - B)-1 = I + B + B2 + ······ + Bk + ······ 事实上 ( I - B)( I + B + B2 + ······ + Bk ) =I – Bk+1 注3: X(k) =B X(k-1) + f = B(B X(k-2) + f) + f =···· = Bk X(0) + ( I + B + ····+ Bk-1)f ≈ ( I – B )-1 f 9/18 Ans= 1.2604e-005 例 线性方程组 A X = b, 分别取系数矩阵为 试分析Jacobi 迭代法和 Seidel 迭代法的敛散性 D=diag(diag(A1)); B1=D\(D-A1); max(abs(eig(B1))) (1) A1=[1,2,-