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

舒米牛牛

%运用Jacobi迭代拟合出来的u关于x和y的矩阵
5 t" N1 i2 w9 \; I0 A$ xfunction J=jac(A,b,u0,eps)
if nargin==3
    eps=1.0e-8
elseif nargin<3
' ?5 i$ @9 T- L  s8 L" ~    'error'
    return
/ g. S: j% Q. p  ?) V: K! vend
* |8 V4 t2 a/ `0 q; i+ I+ X, X% L
; S* k( R5 v9 b& A1 V3 n- J2 |%定义内部节点矩阵u0
' m3 A" i1 P$ u. s) w" \3 C0 p# yh=1;k=1;
# Q/ o( Q% O& K) a" sx=0:h:17;y=0:k:10;
6 h; X+ @5 x$ ~2 b1 n. L7 M6 e: D3 Q/ ]e=length(x)-2;f=length(y)-2;
u=zeros(e,f);
u0=u;

4 O( x. x% u! N1 z%定义外部节点p
p=zeros(e+2,e+2);
- d: m( A4 H% L$ ]+ I4 ]p(1,1:f+2)=0;p(e+2,1:f+2)=0;      
p(1:e+2,1)=100;p(1:e+2,f+2)=100;
* t$ y9 Q9 O, A- D5 E, N
%定义系数矩阵A
: g! ?3 W5 d7 w, f0 HA=zeros(e*f,e*f);
% E  s- x- F: A- x* ^5 ?1 yB=mat2cell(A,ones(e*f/e,1)*e,ones(e*f/e,1)*e);
( N: G+ m) X: X0 ~d1=ones(e,1);d2=ones(e-1,1);
M=4*diag(d1)-diag(d2,1)-diag(d2,-1);
6 p+ O. T/ K' h% xN=-eye(e);
2 V: l5 T) N, h+ HB{1,1}=M;B{e,e}=M
for i=2:e-1
* f3 }  G  @) T0 L3 N' {  i    B{i,i}=M;B{i-1,i}=N;B{i+1,i}=N
) l, \( x0 Z3 w0 M' J: ]end
A=cell2mat(B);  
4 J7 z' Z; k* G% }# G% F% j这里总是显示
??? function J=jac(A,b,u0,eps)
' m% F5 f! G5 ?& \    |
1 J5 t; X6 U, g" ?! c( P- IError: Function definitions are not permitted at the prompt or in scripts.
+ o& {3 `/ z' R
* m" R. }+ x, Q+ N; X" W%定义b
b=ones(e*f,1);
for i=1:e
/ `& o) W$ A2 `4 x% T& u8 l6 [    for j=1:f
( a. Y; h3 a4 }3 ^% I        b(i+j)=p(i,j+1)+p(i+1,j)+p(i+2,j+1)+p(i+1,j+2)
% \5 W7 J/ @  Z# u0 j" H' V) o    end
end
%运用Jacobi迭代法计算
D=diag(diag(A));
. ]* b% I, B, r! SD=inv(D);
L=tril(A,-1);
6 }8 `: p: S2 SU=triu(A,1);
B=-D*(L+U);
  O% D, c9 [" c/ l4 tf=D*b;
J=B*u0+f;
# R" f+ Y" L8 K  @' W& J& W& @' }while norm(J-u0)>=eps
5 d+ N8 n! @  U8 D3 M! N0 q. `% fx0=J;
J=B*u0+f;
& J( T4 B0 o1 v' Q. L. R1 w: h4 U2 uend
return

舒米牛牛

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

BenCam

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

madio

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