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

shuddkk

function [x,k,err]=Jacobi(A,b,n)
  V2 }' p  b' e%Input A:系数矩阵
%Input b:右端向量
2 m- y9 X  W( ?1 ^' B* d3 u%Input n:系数矩阵的阶数
%Input epsilon:计算精度
+ p" C" k9 R1 _7 Q- R% v+ N%Input K:最大迭代次数
%Input x0:初始向量
%Output x:迭代近似解
epsilon=1e-6;
K=5000;
x0=rand(n,1);
k=1;
; Z4 a  z3 v5 _& _2 J2 ~9 v6 p5 Cx=x0;
- V: P8 z  Z9 L9 ^& {2 C' y3 owhile k<K
" _0 ^* b% d/ y. R# w! d6 g    for i=1:n
        s=b(i);
        for j=1:n
: W" ~% I7 Z2 A- w            s=s-A(i,j)*x(j);
        end
        s=(s+A(i,i)*x(i))/A(i,i);
        x(i)=s;
' ]  @6 Z& E, A$ l$ u. _    end
# R! i$ v. p6 A' R8 _: s    if norm(x-x0,inf)<epsilon
        err=norm(x-x0,inf);  
/ ~: T, b" ^5 ~" |$ B% G- p- H" ^        break;
8 L3 R$ ~; T& j; O& h1 L    end
    k=k+1;
( y6 y7 i( z2 Kend
5 {! h  A/ `( A4 zend

  \8 @: u7 b8 X! g. b
. K3 G  P% J, D# T3 f
& N$ u7 z& `+ I( }
2 [, G# E; W* @! O3 I
在命令窗口输入 Jacobi([10,5,0,0;5,10,-4,0;0,-4,8,-1;0,0,-1,5],[6;25;-11;-11],4)后
# _% V7 G9 u- m( r, q只得到了x的值，怎么输出k和err？

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
8 K) n, w) E; I/ h像k = 5000这样的等式可以写在函数参数里吗

应该可以吧

pTDbn25

来学习一下