数值分析实验报告 jacobi迭代和seidel迭代分析.docVIP

数值分析实验报告 实验目的 了解熟悉jacobi迭代法和seidel迭代法的解法 将原理与matlab语言结合起来，编程解决问题 分析实验结果 二、实验题目 设线性方程组为 考察用jacobi迭代法和seidel迭代法求解该线性方程组的收敛情况。如果收敛，给出误差满足的解。 三、实验原理 将A 作如下分解 这里 Jacobi迭代矩阵为 Seidel迭代矩阵为 它们的迭代格式都可化为 则迭代格式对任何初值都瘦脸的充要条件是迭代矩阵谱半径 Jacobi迭代矩阵的特征方程为 Seidel迭代矩阵的特征方程为 四、实验内容 用matlab编写计算jacobi迭代矩程序，建立m文件如下： function[M]=BJ(A) D=diag(diag(A)); L=tril(-A)+D; U=triu(-A)+D; M=inv(D)*(L+U); 输入： A=[20 1 -3 -1;3 18 0 7;-1 2 40 -2;1 0 -1 5]; [M]=BJ(A) M: M = 0 -0.0500 0.1500 0.0500 -0.1667 0 0 -0.3889 0.0250 -0.0500 0 0.050