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

舒米牛牛

%运用Jacobi迭代拟合出来的u关于x和y的矩阵
function J=jac(A,b,u0,eps)
if nargin==3
    eps=1.0e-8
elseif nargin<3
    'error'
    return
' n% J) e- i  w+ Rend

& p2 L" S5 S' Y%定义内部节点矩阵u0
/ h% y5 @! L$ x$ }0 P5 Oh=1;k=1;
/ O7 n& W. u" ~0 b- T: S0 {% _x=0:h:17;y=0:k:10;
  |; e5 L$ A0 \" u3 G+ Me=length(x)-2;f=length(y)-2;
u=zeros(e,f);
' h- q) ^7 B9 V: B/ t% u+ t5 Yu0=u;

%定义外部节点p
! _. t4 e8 \1 ?* n. Y9 ~. l' @* Kp=zeros(e+2,e+2);
p(1,1:f+2)=0;p(e+2,1:f+2)=0;      
- o9 ]4 b. P( G; Jp(1:e+2,1)=100;p(1:e+2,f+2)=100;
( O) T6 K7 X) B3 ], Q6 I" ~
! T1 p* ]: m3 \0 Q%定义系数矩阵A
$ t1 w6 ]) u' u. a/ @A=zeros(e*f,e*f);
( W  q. g( W, L- k3 t; gB=mat2cell(A,ones(e*f/e,1)*e,ones(e*f/e,1)*e);
d1=ones(e,1);d2=ones(e-1,1);
1 S3 v5 B: B5 _' |1 ]  E! kM=4*diag(d1)-diag(d2,1)-diag(d2,-1);
- W" o. D1 R, j- K; nN=-eye(e);
+ C- d' }  C4 g* G0 ]2 bB{1,1}=M;B{e,e}=M
for i=2:e-1
    B{i,i}=M;B{i-1,i}=N;B{i+1,i}=N
end
4 ]% a7 v- c3 k6 t; X+ }A=cell2mat(B);  
0 j5 `1 A% ]: ^( {* Q8 Z: ?$ Y& a( c这里总是显示
+ I6 U% u8 G; R) ]2 v) _??? function J=jac(A,b,u0,eps)
    |
Error: Function definitions are not permitted at the prompt or in scripts.

%定义b
( D7 I  C9 Q- rb=ones(e*f,1);
7 s" I- [* T. }, x1 Rfor i=1:e
2 V' k2 ~$ v: T, r: L    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
- [0 J- V& O  m8 Gend
6 ]' U9 W. y& y- L; W+ d%运用Jacobi迭代法计算
D=diag(diag(A));
8 H* W. K. N: }: GD=inv(D);
5 }6 z# \2 V' f0 G; CL=tril(A,-1);
U=triu(A,1);
B=-D*(L+U);
4 I% W7 m! a, bf=D*b;
J=B*u0+f;
while norm(J-u0)>=eps
5 A. ]  M/ [8 @3 M" U. `x0=J;
J=B*u0+f;
end
return

madio

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

BenCam

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

舒米牛牛

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