数值计算实验报告new.docVIP

数值计算方法实验 实验名称： 解线性方程组的迭代法 班 级： 宿 舍： 成 员: 实验三 解线性方程组的迭代法 实验目的 用雅克比和高斯-赛尔迭代法解线性方程组Aｘ＝b,式中A为非奇异矩阵。再给定迭代初值的情况下，进行迭代，知道满足精度要求。 实验方法 （1）雅克比迭代法 实验程序 雅克比迭代法的MATLAB函数文件agui_jacobi.m如下。 function x=agui_jacobi(a,b) n=length(b); N=100; e=1e-4; x0=zeros(n,1); x=x0; x0=x+2*e; k=0; d=diag(diag(a)); l=-tril(a,-1); u=-triu(a,1); while norm(x0-x,inf)ekN k=k+1; x0=x; x=inv(d)*(1+u)*x+inv(d)*b; k disp(x) end if k==N warning(已达最大迭代次数); end 2、雅克比迭代法源程序 在MATLAB命令窗口输入如下。 a=[4 -1 0 -1 0 0;-1 4 -1 0 -1 0;0 -1 4 -1 0 -1;-1 0 -1 4 -1 0;0 -1 0 -1