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

shuddkk

function [x,k,err]=Jacobi(A,b,n)
%Input A:系数矩阵
%Input b:右端向量
* l7 ]1 r' i$ \9 \6 Y$ X%Input n:系数矩阵的阶数
%Input epsilon:计算精度
: ^$ u/ ]8 ?; q8 Q) ]%Input K:最大迭代次数
0 d4 J! i0 Q  K5 r' E5 Z%Input x0:初始向量
$ n/ B* l, V5 g- P) n%Output x:迭代近似解
epsilon=1e-6;
% |( B% o/ W5 Q7 @K=5000;
6 E" a/ p7 @* ]& {0 nx0=rand(n,1);
k=1;
x=x0;
while k<K
    for i=1:n
, i# z/ l; t. F4 a        s=b(i);
; u; W) }1 R& T        for j=1:n
            s=s-A(i,j)*x(j);
: I' O0 k/ A: A; U        end
4 E! M+ F, Q7 R+ `" \6 w! O0 Z        s=(s+A(i,i)*x(i))/A(i,i);
        x(i)=s;
  C) D+ U4 I- T" [- K+ ~1 W    end
' ]' G. S3 i, K' }- N3 a    if norm(x-x0,inf)<epsilon
, w: D" S4 F) F% T1 A        err=norm(x-x0,inf);  
        break;
    end
    k=k+1;
9 p; Q4 v9 X% e" ?" _, u8 dend
8 r1 A) i, e3 Zend
. F3 f& B, @0 o0 ^; B
1 J4 M, Z1 M; [8 |# t* c' H4 F
* z6 D, C( C/ d2 Z

$ _3 @9 B0 z+ G$ \8 f" X: L3 [" {- |" z
在命令窗口输入 Jacobi([10,5,0,0;5,10,-4,0;0,-4,8,-1;0,0,-1,5],[6;25;-11;-11],4)后
, ~; r; ~8 j& Q9 U+ d7 z& ]& x* Y只得到了x的值，怎么输出k和err？
1 |. |  P: G. H' q  k

Uifhjvv

需要按方括号的格式调用后面的输出
. R6 B. u6 @+ R% K5 r. {; z7 d[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
像k = 5000这样的等式可以写在函数参数里吗

应该可以吧

pTDbn25

来学习一下