我自定义函数[x,k,err]=Jacobi(A,b,n)后怎么才能输出k和err？

shuddkk

function [x,k,err]=Jacobi(A,b,n)
5 O9 N+ g2 i; l  K: z1 N& ~%Input A:系数矩阵
+ k  W% |0 {8 j; G1 j+ V%Input b:右端向量
, t  Y0 G3 n: X; ]5 V. K3 ^%Input n:系数矩阵的阶数
%Input epsilon:计算精度
%Input K:最大迭代次数
" V; I4 k" n- D) t, w' F1 @%Input x0:初始向量
. \! b7 L+ k9 d# t%Output x:迭代近似解
epsilon=1e-6;
) z) c; s2 m9 H% [$ i0 {K=5000;
* F& s. [# ~. o: T2 a1 j; M% |/ Gx0=rand(n,1);
) G- {. ?; ]2 W% Z  _; Z$ D/ dk=1;
1 Y# H, h. ?7 d. z# C, Vx=x0;
while k<K
    for i=1:n
        s=b(i);
/ e/ b/ _* _. A* ^        for j=1:n
+ p- R" N% V% [: j$ \  F! K            s=s-A(i,j)*x(j);
: Q3 N: ^& p6 [$ O        end
        s=(s+A(i,i)*x(i))/A(i,i);
+ C. b! t# v9 G; T  K! }5 e        x(i)=s;
    end
    if norm(x-x0,inf)<epsilon
% s2 o" q" D6 Z" C- a5 I        err=norm(x-x0,inf);  
6 i1 l/ c7 P- v9 R$ z6 v& v        break;
    end
    k=k+1;
  r' j/ j, r7 U$ \; @6 `5 J1 ^- D% gend
( }& l; w) [9 D/ ^5 send
' c, c- H" b6 A5 y
* H: }# A  y, Q: M
- q  r% a* x% A! d

8 `1 _  N% L; H! \. X" u
# @! O8 G6 j/ F" G/ y  f在命令窗口输入 Jacobi([10,5,0,0;5,10,-4,0;0,-4,8,-1;0,0,-1,5],[6;25;-11;-11],4)后
' B" X; ^9 z# c  C# N: M只得到了x的值，怎么输出k和err？
1 G" H: F" i4 B8 P
6 @% R7 _& x/ j$ `2 o

Uifhjvv

需要按方括号的格式调用后面的输出
[x,k,err]=Jacobi([10,5,0,0;5,10,-4,0;0,-4,8,-1;0,0,-1,5],[6;25;-11;-11],4)

Uifhjvv

像k = 5000这样的等式可以写在函数参数里吗

cichishia

Uifhjvv 发表于 2020-7-8 13:43
- ?) s) X( m2 c2 w& H% r像k = 5000这样的等式可以写在函数参数里吗

应该可以吧

pTDbn25

来学习一下