求助：关于分块矩阵的还原，为什么实行不了

舒米牛牛

%运用Jacobi迭代拟合出来的u关于x和y的矩阵
. S0 z5 I9 A# D- o# M# }2 G5 Lfunction J=jac(A,b,u0,eps)
$ M8 r4 y6 h& A9 c7 G/ A( {if nargin==3
    eps=1.0e-8
0 O, b  _/ p, E( nelseif nargin<3
    'error'
    return
1 }5 j8 b8 _0 O( \$ R& Hend
8 T+ G' t2 f; _( M! I8 V- ~  [
; ^& ^: T6 M* c. ~0 n%定义内部节点矩阵u0
7 q) @# m6 m$ m" d+ z6 v+ C' f; |h=1;k=1;
% @1 l  ?: R) M/ ^6 ], E1 W0 [x=0:h:17;y=0:k:10;
e=length(x)-2;f=length(y)-2;
$ q* Z9 j& C1 G$ }u=zeros(e,f);
u0=u;

%定义外部节点p
! b, Q( a2 o- J! Ap=zeros(e+2,e+2);
p(1,1:f+2)=0;p(e+2,1:f+2)=0;      
" Z5 `. }* y8 d! Lp(1:e+2,1)=100;p(1:e+2,f+2)=100;
9 B' V) ~! \+ m! k+ @' x4 b
%定义系数矩阵A
$ q. z6 |2 F7 l' t$ ~A=zeros(e*f,e*f);
9 M, S/ o4 t2 \% F, P0 c5 YB=mat2cell(A,ones(e*f/e,1)*e,ones(e*f/e,1)*e);
d1=ones(e,1);d2=ones(e-1,1);
M=4*diag(d1)-diag(d2,1)-diag(d2,-1);
N=-eye(e);
: l& f- ^9 P0 S# v7 tB{1,1}=M;B{e,e}=M
for i=2:e-1
- k( R% H5 r6 a9 @$ A; T. ?( k    B{i,i}=M;B{i-1,i}=N;B{i+1,i}=N
4 h" H# G; c, p6 ?0 f" Wend
A=cell2mat(B);  
! X; H: ~4 Y7 P( G+ _6 n这里总是显示
4 \( t: C6 J5 ]5 c! y??? function J=jac(A,b,u0,eps)
# v) ^+ ^. T3 [& c    |
! U+ m8 U7 H6 h( oError: Function definitions are not permitted at the prompt or in scripts.
' i9 b9 r0 B9 X" J: P, ^) J
  t. E, g; Q; M* M%定义b
b=ones(e*f,1);
for i=1:e
' H+ i; i& o5 M# l. |% T    for j=1:f
" g4 f/ u( C  x        b(i+j)=p(i,j+1)+p(i+1,j)+p(i+2,j+1)+p(i+1,j+2)
$ l9 U% X. O/ m: Q1 c: a5 a    end
* N2 ?) S  g% A& W- X  U- M3 Eend
+ g1 s$ A. u7 x' g- N; M" C%运用Jacobi迭代法计算
$ D# `; [5 T" \+ ~( U* xD=diag(diag(A));
* n- l, L* P" b( jD=inv(D);
L=tril(A,-1);
U=triu(A,1);
& Y, E2 v9 F6 r6 U! e8 KB=-D*(L+U);
f=D*b;
$ W7 {  b# |: ~5 d* B+ DJ=B*u0+f;
while norm(J-u0)>=eps
+ X! @4 @8 C0 `5 @x0=J;
9 k% K, e  I) c1 iJ=B*u0+f;
end
1 u) `  k. J, g0 Nreturn

舒米牛牛

自己顶一下，拜托哪位高手指点一下，纠结这个矩阵的还原，想了好多方法还是不行~~实在想不出哪里出错了

BenCam

对不起,我也不知道,帮不了忙!

madio

你是不是函数的定义没有放在M文件中？