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

舒米牛牛

%运用Jacobi迭代拟合出来的u关于x和y的矩阵
function J=jac(A,b,u0,eps)
4 `1 \- b1 |) x$ C9 y/ }if nargin==3
) O7 g. M, J* ~: K2 Q; g% O8 q    eps=1.0e-8
* ^- y7 _; ^& E4 @3 nelseif nargin<3
    'error'
    return
4 _- _7 Z8 _) ^3 xend

%定义内部节点矩阵u0
9 ~7 Y( Q! R5 Dh=1;k=1;
x=0:h:17;y=0:k:10;
0 ?6 z) |) m4 ]( Me=length(x)-2;f=length(y)-2;
5 ^# k% B' Q1 ^7 S- ju=zeros(e,f);
+ M- E/ ]1 R+ d# w# }0 x+ Zu0=u;
% |  `, R  f- A- b" G
% {, A0 H  q$ k) m: l$ x%定义外部节点p
8 b2 ?% B: a% n( tp=zeros(e+2,e+2);
/ N  I) G5 u# o, q+ fp(1,1:f+2)=0;p(e+2,1:f+2)=0;      
% p; i9 V9 f8 Z; mp(1:e+2,1)=100;p(1:e+2,f+2)=100;
* E: K: Q1 J4 x
%定义系数矩阵A
$ {% G+ q5 x3 D7 M5 l6 BA=zeros(e*f,e*f);
B=mat2cell(A,ones(e*f/e,1)*e,ones(e*f/e,1)*e);
& S4 t2 C# |% G* |; r/ Dd1=ones(e,1);d2=ones(e-1,1);
M=4*diag(d1)-diag(d2,1)-diag(d2,-1);
& r6 L) H' g& `, H9 d6 x/ i/ }, ~! `N=-eye(e);
B{1,1}=M;B{e,e}=M
2 C& t0 ~3 X$ N  @1 H. @for i=2:e-1
    B{i,i}=M;B{i-1,i}=N;B{i+1,i}=N
4 ~0 _6 u9 |5 h% ?end
A=cell2mat(B);  
( h) R4 }% |4 r: j( t) v. x. I这里总是显示
??? function J=jac(A,b,u0,eps)
    |
: x- m# ?, u0 qError: Function definitions are not permitted at the prompt or in scripts.
" U) c8 S9 U, U; _7 O6 g* i
* `: p" x- s* c( i- F0 x  V. Z%定义b
b=ones(e*f,1);
for i=1:e
  r' S/ m  E- P* y0 h+ T    for j=1:f
        b(i+j)=p(i,j+1)+p(i+1,j)+p(i+2,j+1)+p(i+1,j+2)
    end
end
5 \* A: N$ ?0 |! V%运用Jacobi迭代法计算
0 O5 H2 n4 F" O, d+ `7 z0 L  Z( BD=diag(diag(A));
D=inv(D);
L=tril(A,-1);
U=triu(A,1);
- @" o, ^) `! T- N6 _. AB=-D*(L+U);
f=D*b;
J=B*u0+f;
while norm(J-u0)>=eps
* f  G1 w" {$ G7 J6 F6 b/ w  tx0=J;
2 d: [% _1 o! z/ b! y% vJ=B*u0+f;
end
return

舒米牛牛

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

BenCam

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

madio

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