gauss.m

来自「数值分析常用的一些公式」· M 代码 · 共 19 行

function s=Gauss(a,b,x0)
 %a 为系数矩阵，b为右端向量，x0为初始值
 esp=1.0e-3;
 D=diag(diag(a));
 L=tril(a,-1);%求出严格下三角矩阵
 U=triu(a,1); %求出严格上三角矩阵
 C=inv(D+L);
 B=-C*U;
 f=C*b;
 s=B*x0+f;
 while norm(s-x0)>=esp
     x0=s;
     s=B*x0+f;
 end
 return
 %a=[8 -3 2;4 11 -1;6 3 12];
 %b=[20 33 36]';
 %x0=[0 0 0]';
 %Gauss(a,b,x0)