部分解线性方程组迭代法.pptVIP

第六部分 解线性方程组的迭代法 求解 迭代法 从一个初始向量出发,按照一定的递推格式,产生逼近方程组的近似解序列。 思路 与解f (x)=0 的不动点迭代相似 ， 将方程组 等价改写成 形式，从而建立迭代格式　 ，从 出发，生成迭代序列 例1：利用Jacobi和Gauss-Seidel迭代法求解方程组 解： Jacobi迭代格式 Jacobi迭代法 一、 Jacobi迭代法 设方程组 将系数矩阵分裂为： 其中 如果 原方程组可化为 其中 相应的迭代格式 上述方法称为Jacobi迭代法 分量形式： A=[20 2 3;1 8 1;2 -3 15]; b=[24;12;30]; B=zeros(3); for i=1:3 B(i,i)=A(i,i); A(i,i)=0; end B A B = 20 0 0 0 8 0 0 0 15 A = 0 2 3 1 0 1 2 -3 0 x(1:3,1)=[1.2 1.5 2]; for i=2:7 x(1:3,i)=inv(B)*(-A)*x(1:3,i-1)+inv(B)*b; end x x = 1.2000 0.75