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

舒米牛牛

%运用Jacobi迭代拟合出来的u关于x和y的矩阵
function J=jac(A,b,u0,eps)
- l, g& `% C# v) c1 U& ~$ Cif nargin==3
* Y. G, t: ]- B  [    eps=1.0e-8
elseif nargin<3
2 U" }2 h" v+ h, r$ v# J$ N    'error'
    return
+ Z6 i6 c. x. a& e1 R( Y7 F9 k- Wend

5 r9 m7 s+ v* E5 |& }8 ?%定义内部节点矩阵u0
! z: i+ h" T$ z5 B6 l; F1 Z0 `h=1;k=1;
x=0:h:17;y=0:k:10;
& S" B4 @  k: B3 i: Ee=length(x)-2;f=length(y)-2;
u=zeros(e,f);
u0=u;
7 b6 y/ V: d" P7 @
. H7 i  I- ^$ w/ W( s. H( Q6 W%定义外部节点p
" K2 }: }+ N% H# n2 H9 tp=zeros(e+2,e+2);
* `" e- [! @' B% J" e! j8 |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;

%定义系数矩阵A
A=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);
M=4*diag(d1)-diag(d2,1)-diag(d2,-1);
% q& U! u( q7 W2 m( V$ \N=-eye(e);
B{1,1}=M;B{e,e}=M
4 [: T/ @0 }3 j. K3 Pfor i=2:e-1
    B{i,i}=M;B{i-1,i}=N;B{i+1,i}=N
end
A=cell2mat(B);  
; |% ]6 }& M* _+ W6 a这里总是显示
9 j4 C; O$ r+ B% b??? function J=jac(A,b,u0,eps)
2 |" Q' d& P0 F1 e8 o9 ?    |
Error: Function definitions are not permitted at the prompt or in scripts.

$ V( x' d( [) j- S5 y%定义b
9 c0 @- x7 ~* j1 n0 Yb=ones(e*f,1);
+ t* v9 S' o4 D" _' {. E, d& ufor i=1:e
' y2 ?* H2 C: [    for j=1:f
  r: D1 J/ Y0 ?* {        b(i+j)=p(i,j+1)+p(i+1,j)+p(i+2,j+1)+p(i+1,j+2)
; s) q) Q/ m0 W; T, Y2 q  I+ r0 e    end
+ l5 f6 b7 ~' s, uend
8 G, j9 h0 c2 l# w9 y* B%运用Jacobi迭代法计算
0 M( M" h, c+ i- [: sD=diag(diag(A));
D=inv(D);
2 A( m+ Z( H3 U% |L=tril(A,-1);
" {: {; W0 T8 K: N6 U1 p, J6 sU=triu(A,1);
, `5 F- w) f3 WB=-D*(L+U);
f=D*b;
$ ?7 i% F5 T$ OJ=B*u0+f;
% J) V! c' q& g1 h: H, |while norm(J-u0)>=eps
x0=J;
J=B*u0+f;
end
5 |+ d1 F8 o0 ?! B+ {5 @8 ureturn

madio

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

BenCam

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

舒米牛牛

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