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

舒米牛牛

%运用Jacobi迭代拟合出来的u关于x和y的矩阵
function J=jac(A,b,u0,eps)
if nargin==3
    eps=1.0e-8
& m) [  N9 N2 r( S  G  Y" C! g  Oelseif nargin<3
) `& h$ o) u6 m& i    'error'
- d+ w2 V- `' m; U0 D7 p. ]( ]    return
% h; o3 I) b6 Eend

%定义内部节点矩阵u0
! M3 q0 R5 I  Rh=1;k=1;
x=0:h:17;y=0:k:10;
e=length(x)-2;f=length(y)-2;
u=zeros(e,f);
u0=u;
3 D1 T3 R% W" ~$ J$ n
%定义外部节点p
p=zeros(e+2,e+2);
- `: l& u3 I' u: S9 sp(1,1:f+2)=0;p(e+2,1:f+2)=0;      
! I' @$ w5 M& Cp(1:e+2,1)=100;p(1:e+2,f+2)=100;

%定义系数矩阵A
A=zeros(e*f,e*f);
# b" Y4 s4 v2 i: y! ^# o' z' q# t  NB=mat2cell(A,ones(e*f/e,1)*e,ones(e*f/e,1)*e);
# T- F+ }* f) q4 Kd1=ones(e,1);d2=ones(e-1,1);
4 Z; v2 ~  c9 l: O' \$ }5 A( }M=4*diag(d1)-diag(d2,1)-diag(d2,-1);
N=-eye(e);
B{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
$ p$ V" ]  A/ g) L1 ZA=cell2mat(B);  
( u4 c& P1 }2 L( ]1 K* p这里总是显示
* W9 D# k0 K' B# ]4 k??? function J=jac(A,b,u0,eps)
    |
! q, K6 m1 z( V: HError: Function definitions are not permitted at the prompt or in scripts.

' L: f1 \2 A7 h  C3 Z%定义b
/ P5 z  t6 j& Y: L* ^b=ones(e*f,1);
) _5 @# m" w5 i# ~for i=1:e
) |3 Q- ~. @: |' w$ @    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
: H9 M. m; O; F  {5 g8 M- Cend
+ Y  |2 y% P) b0 e& h% ~1 M%运用Jacobi迭代法计算
& I* K+ B: b) B1 _; r7 @D=diag(diag(A));
- H6 q! u) I1 a; ND=inv(D);
4 {1 k1 H  [5 \4 R! z+ KL=tril(A,-1);
4 Y& u2 E+ M% `4 r* a; uU=triu(A,1);
5 f# C" x3 Y# Y8 A3 tB=-D*(L+U);
f=D*b;
J=B*u0+f;
. D7 s& Q5 N3 g/ k$ Awhile norm(J-u0)>=eps
x0=J;
J=B*u0+f;
; {! {' j. B+ H" rend
return

舒米牛牛

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

BenCam

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

madio

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