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

舒米牛牛

%运用Jacobi迭代拟合出来的u关于x和y的矩阵
function J=jac(A,b,u0,eps)
/ R+ f8 }( L1 H8 Mif nargin==3
    eps=1.0e-8
elseif nargin<3
    'error'
- R: D! j6 N# j1 N3 \' V    return
( b- |, P$ L- K. s% Nend
+ u- O6 F3 i( q2 q. m4 ]
%定义内部节点矩阵u0
. C$ Z; U/ j4 q- Fh=1;k=1;
x=0:h:17;y=0:k:10;
e=length(x)-2;f=length(y)-2;
u=zeros(e,f);
9 G7 b/ Q! b7 T7 x! x, iu0=u;
5 H% @- u8 w' c& {& d! `
%定义外部节点p
* j  d4 S) k  ?" L" ~' Z( |6 f) X# kp=zeros(e+2,e+2);
' I" Q' g5 O/ j+ U: l6 Dp(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;
: C* |3 r0 V& K$ Z. Q2 z5 j8 d
) E! y' K, W5 ^8 h% F%定义系数矩阵A
9 y6 x4 J' t1 ?( F" L3 U5 AA=zeros(e*f,e*f);
B=mat2cell(A,ones(e*f/e,1)*e,ones(e*f/e,1)*e);
d1=ones(e,1);d2=ones(e-1,1);
$ ]2 q3 s/ Z. iM=4*diag(d1)-diag(d2,1)-diag(d2,-1);
N=-eye(e);
* Y4 n. C. {. W! ^/ KB{1,1}=M;B{e,e}=M
7 o1 S. B: `6 \- Efor i=2:e-1
    B{i,i}=M;B{i-1,i}=N;B{i+1,i}=N
end
A=cell2mat(B);  
  ]& J- h3 K) B% G这里总是显示
??? function J=jac(A,b,u0,eps)
    |
. Y& e6 {* [/ _' p% U* y$ @7 jError: Function definitions are not permitted at the prompt or in scripts.

) T- M9 A+ H2 x  X%定义b
b=ones(e*f,1);
for i=1:e
' J+ X" ^9 O6 A" c6 w' r    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
$ x/ \7 H9 h. h2 k8 |$ a! |7 A%运用Jacobi迭代法计算
D=diag(diag(A));
D=inv(D);
L=tril(A,-1);
U=triu(A,1);
0 S1 u/ b9 d" ~* T% OB=-D*(L+U);
f=D*b;
* q* n( y2 O" J2 zJ=B*u0+f;
while norm(J-u0)>=eps
' n/ H! P" O( ?) i. f2 }x0=J;
J=B*u0+f;
end
return

舒米牛牛

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

BenCam

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

madio

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