数值分析9.pptVIP

*/18 大型稀疏矩阵的背景 Jacobi迭代与Seidel迭代 迭代法的矩阵表示 迭代法数值实验 《数值分析》 9 ? ? ? ? 边值问题: 令 h = 1/(n+1) , xi= ih, yj = jh ( i , j = 0,1, ···, n+1 ) 记 ui,j= u(xi , yj ), ( i , j = 0,1, ···, n+1 ) ? ( n=4 ) AU = F ? 例4.1 特点:系数矩阵主对角元均不为零 计算格式 X(1)=B X(0) + f 取 X(0) = X* 1.0000 1.0000 1.0000 X(0) 0 0 0 X(1) 0.7778 0.8000 0.8667 X(2) 0.9630 0.9644 0.9778 X(3) 0.9929 0.9935 0.9952 计算格式: X(k+1)=BX(k)+f 准确解 ? X(4) ········ 0.9987 0.9988 0.9991 雅可比迭代法 (i = 1,2,…n; k=1,2,……) 取初始向量X(0)=[x1(0) x2(0) ··· xn(0)]T, 迭代计算 (i = 1,2,…,n) 迭代法适用于解大型稀疏方程组 (万阶以上的方程组,系数矩阵中零元素占很大比例