雅克比迭代求解线性方程组MATLAB课设

雅克比迭代求解线性方程组（MATLAB）课设,已经调试成功

function [ynH]=gauseidel（Abx0epsN）
%Gauss-Serdel迭代法的MATLAB函数文件gauseidel.m如下：
%A=[7 2 1 -2;9 15 3 -2;-2 -2 11 5;1 3 2 13];
%b=[4 7 -1 0]‘;x0=[0000]‘;
%[xn]=gauseidel（Ab[000，0]‘1.0e-6）
format long g;
if nargin==5
    digits（abs（log10（eps）））;
elseif nargin==4
    digits（abs（log10（eps）））;
    N=100;
elseif nargin==3
    digits（6）;
    eps=1.0e-6;
    N=100;
elseif nargin<3
    error
    return
end      
D=diag（diag（A））;    %求A的对角矩阵
L=-tril（A-1）;      %求A的下三角阵
U=-triu（A1）;       %求A的上三角阵
G=（D-L）\U;
f=（D-L）\b;
y=G*x0+f;
n=1;                  %迭代次数
while norm（y-x0）>=eps&n    x0=y;
    y=G*x0+f;
    n=n+1;
end
y=vpa（y）;
if n    H=‘迭代收敛‘;
else
    H=‘迭代不收敛‘;
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        727  2012-06-02 12:01  杨勇—设计—雅克比\gauseidel.m

     文件        685  2012-06-02 11:44  杨勇—设计—雅克比\jacobi.m

     文件       3896  2012-06-02 11:46  杨勇—设计—雅克比\untitled.fig

     文件      13812  2012-06-02 11:47  杨勇—设计—雅克比\untitled.m

     目录          0  2012-06-17 20:40  杨勇—设计—雅克比

----------- ---------  ---------- -----  ----

                19120                    5