4.线性方程组迭代法.doc

三、线性方程组迭代法计算实验题（程序）1. 利用Jacobi法求解. 解: 第一步：function x=equjacobi(A,b,x0)% 用Jacobi法解方程组Ax=b，x0为初始向量，epsilon为精度要求(默认1e-6)D=diag(diag(A);U=triu(A,1);L=tril(A,-1);B=-D(L+U);g=Db;x1=B*x0+g;n=1;epsilon=1e-6;while norm(x1-x0,inf)epsilonx0=x1;x1=B*x0+g;n=n+1;endx=x1;n保存为M文件，即equjacobi.m，保存到D:js目录第二步：将工作目录设为D:js在命令窗口输入 a=10 -1 -2;-1 10 -2;-1 -1 5;b= 7.2;8.3;4.2; x0=0;0;0;x=equjacobi(a,b,x0)回车得到：n= 14x = 1.1000 1.2000 1.30002. 利用Gauss-Seidel法求解，与Jacobi法比较。第一步：function x=equseidel(A,b,x0)% 用Gauss-Seidel法解方程组Ax=b，x0为初始向量，epsilon为精度要求(默认1e-6)D=diag(diag(A);U=triu(A,1);L=tril(A,-1);B=-(D+L)U;g=(D+L)b;x1=B*x0+g;n=1;epsilon=1e-6;while norm(x1-x0,inf)epsilonx0=x1;x1=B*x0+g;n=n+1;endx=x1;n保存为M文件，即equseidel.m，保存到D:js目录第二步：将工作目录设为D:js在命令窗口输入 a= a=10 -1 -2;-1 10 -2;-1 -1 5;b= 7.2;8.3;4.2; x0=0;0;0;x=equseidel(a,b,x0)回车得到：n= 9 x = 1.1000 1.2000 1.30003. 利用SOR法求解第一步：function x=equsor(A,b,x0,w)% 用Gauss-Seidel法解方程组Ax=b，x0为初始向量，epsilon为精度要求(默认1e-6)D=diag(diag(A);U=triu(A,1);L=tril(A,-1);B=-(D+L)U;g=(D+L)b;x1=B*x0+g;n=1;epsilon=1e-6;while norm(x1-x0,inf)epsilonx0=x1;x1=B*x0+g;n=n+1;endx=x1;n保存为M文件，即equsor.m，保存到D:js目录第二步：将工作目录设为D:js在命令窗口输入 a= a=10 -1 -2;-1 10 -2;-1 -1 5;b= 7.2;8.3;4.2; x0=0;0;0;x=equseidel(a,b,x0,1)回车得到：n= 9 x = 1.1000 1.2000 1.30004. 利用SOR法求解并列出当松弛因子w=1.0, 1.1,