• 2022-05-27 问题

    以下函数用fjacobi.m文件名存盘function X=fjacobi(A,b,X0)D=diag(diag(A));% tril(A)求矩阵的左下三角,tril(A,-1)就不包括主对角线。L=-tril(A,-1);U=-triu(A,1);B=【1】;F=D\b;X=B*X0+F;n=1;m=30;while norm(X-X0)>;=1e-5 &amp;&amp; n<;=mX0=【2】;X=【3】;n=n+1;end在命令窗口输入a,b,x0,调用fjacobi函数得到结果>;>; a=&#91;10 -2 -1;-2 10 -1;-1 -2 5&#93;;>;>; b=&#91;3;15;10&#93;;>;>; x0=&#91;0 0 0&#93;';>;>; x=【4】x =1.00002.00003.0000

    以下函数用fjacobi.m文件名存盘function X=fjacobi(A,b,X0)D=diag(diag(A));% tril(A)求矩阵的左下三角,tril(A,-1)就不包括主对角线。L=-tril(A,-1);U=-triu(A,1);B=【1】;F=D\b;X=B*X0+F;n=1;m=30;while norm(X-X0)>;=1e-5 &amp;&amp; n<;=mX0=【2】;X=【3】;n=n+1;end在命令窗口输入a,b,x0,调用fjacobi函数得到结果>;>; a=&#91;10 -2 -1;-2 10 -1;-1 -2 5&#93;;>;>; b=&#91;3;15;10&#93;;>;>; x0=&#91;0 0 0&#93;';>;>; x=【4】x =1.00002.00003.0000

  • 1