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

舒米牛牛

%运用Jacobi迭代拟合出来的u关于x和y的矩阵
4 f7 Q2 Y% k# h3 U5 k, g+ @function J=jac(A,b,u0,eps)
' ]6 @6 Q4 s2 s  I! bif nargin==3
5 Q: ]2 `* @) u$ B. l8 J    eps=1.0e-8
elseif nargin<3
    'error'
. N  a* a+ A' I1 s% J! I    return
* P3 m5 Y5 m' ]end
0 l( R! ~- H9 J( w- k
6 r+ U+ L* @. b  f' i%定义内部节点矩阵u0
h=1;k=1;
- I; q8 }+ ]8 ix=0:h:17;y=0:k:10;
: J% i$ s' N7 v; y( o3 _- W, w2 Ge=length(x)-2;f=length(y)-2;
u=zeros(e,f);
! y) f, {% b' O' `+ z9 K. \u0=u;
4 _+ Y; ^* M! `& f  ~& p
%定义外部节点p
p=zeros(e+2,e+2);
2 B# O4 _" C+ K! P& _, o* Fp(1,1:f+2)=0;p(e+2,1:f+2)=0;      
( ~* Q8 e' E2 k8 F6 G5 Q7 V: T8 Hp(1:e+2,1)=100;p(1:e+2,f+2)=100;
+ l& ?# Q3 Y3 r* K% C
) j' {+ m& U( j%定义系数矩阵A
* ~5 P$ f  h3 B0 dA=zeros(e*f,e*f);
B=mat2cell(A,ones(e*f/e,1)*e,ones(e*f/e,1)*e);
& H" G5 ]8 S8 J# fd1=ones(e,1);d2=ones(e-1,1);
- S9 U- |9 \3 r0 s9 sM=4*diag(d1)-diag(d2,1)-diag(d2,-1);
) B* x, l5 q% t( N9 e% _& X% |N=-eye(e);
B{1,1}=M;B{e,e}=M
* D' {; b/ f' X9 Y9 C' H  J+ `& y7 @for i=2:e-1
    B{i,i}=M;B{i-1,i}=N;B{i+1,i}=N
end
; s: Y; j( {% m( Q" \A=cell2mat(B);  
这里总是显示
) V7 Q" t" V1 g1 P2 y2 c??? function J=jac(A,b,u0,eps)
    |
Error: Function definitions are not permitted at the prompt or in scripts.

%定义b
1 ?7 O0 }, v+ v2 a3 k6 N! [) }! lb=ones(e*f,1);
: z' P0 K8 Z' f# q$ Jfor i=1:e
5 K$ q4 y  w* [3 w    for j=1:f
. u# {: g; X. A! v7 x        b(i+j)=p(i,j+1)+p(i+1,j)+p(i+2,j+1)+p(i+1,j+2)
/ ]6 h& s0 S3 v( J    end
7 {0 F2 D& e% R. Q" V/ aend
* G+ _, h0 c' T0 X%运用Jacobi迭代法计算
/ K2 D: c  m% oD=diag(diag(A));
! A8 k$ h( l4 T) n8 P0 ]& @/ Q# SD=inv(D);
L=tril(A,-1);
U=triu(A,1);
B=-D*(L+U);
f=D*b;
J=B*u0+f;
- b7 s* d0 j1 D) v# cwhile norm(J-u0)>=eps
x0=J;
J=B*u0+f;
end
return

舒米牛牛

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

BenCam

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

madio

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