数学建模社区-数学中国

标题: 急求Matlab程序分析解释 [打印本页]

作者: 范婷    时间: 2010-11-25 15:45
标题: 急求Matlab程序分析解释
有哪位高手可以帮忙分析一下这个程序!!!!!!!!!
% P: p' A* S! K2 pr=1;%点光源的能量值,可以修改。%
$ E! h8 k; w: [# L  \A=zeros(81,81);  %方程的系数矩阵。%
% `; m1 Q4 w2 b% _" qfor i=1:819 X4 R4 A- p  d5 c, A, `
    if rem(i-1,9)~=0$ K" w! ?! q1 Q+ t. v
        A(i,i-1)=1/4;: t4 a. S: H1 g
    end
' j, ~+ `/ T: [8 H# U3 M    A(i,i)=-1;6 I- K" P4 k# h  o9 Y- O
    if rem(i,9)~=0# d  Y! h! r/ J2 ]
        A(i,i+1)=1/4;
" B* O1 h# E& l* B/ \' r) ^; j1 ^2 z    end8 d, b7 I; d$ e  V2 m# ?% x
    if (i-9)>0
! F2 A% q! z! B        A(i,i-9)=1/4;
0 d! ?. q% C) p7 h: H0 J! o$ p& V    end% w. e9 V+ `& g7 T6 v8 N
    if (i+9)<82" Y4 ~( G4 C% _3 w' R
        A(i,i+9)=1/4;
4 k/ n3 B4 a1 b/ H' L    end2 F6 p: W% [7 V, t
end* L) n: }5 S+ `6 e2 ~4 @! n
B=zeros(81,1);7 p/ e) n! k$ A! X: ?
B(6)=-r/4;   
2 J/ T+ g6 C& J% cX=inv(A)*B;
  _5 g; }5 Y3 u8 [1 K! CD=reshape(X,9,9);
0 Y- Y7 x$ [9 V; x/ _# aD=D';' c# c1 T( r7 s" Z5 ~
E=D;
% \: j% }' O/ b" d: R. `x=zeros(1,11);
6 H# e0 W# E1 a. wy=zeros(9,1);8 T* G( ~/ U! v5 W0 a* Q; C; k/ Y
D=[y,D,y];
* w1 i, \- U% j/ {D=[x;D;x];
0 Y. k! i) ^- v; D8 g, ND(1,6)=1;$ O$ j* A  Z* K& X+ [0 W
xt=0:10;
9 {+ @; [: Z( y, a  z. o: f" L: E2 t% ][xt,yt]=meshgrid(xt);+ u1 T9 h& \& ~( _( z7 F4 |
zt=D;8 u! v/ D& V, Y+ p
surf(xt,yt,zt);6 n; R, |, h8 V* R3 f, L; G/ C
D      %矩阵D上的值对应每个节点上的能量值。%
' `1 S  ~1 ^4 s: L! y
作者: hy282694    时间: 2010-11-26 08:20
你要分析啥?带着你看懂源代码?还是。。。?
作者: fhfhfh555    时间: 2010-11-26 09:41
1111111111111111111
作者: 陈阳康    时间: 2010-11-26 13:19
只有报纸坐了。
作者: 陈阳康    时间: 2010-11-26 13:20
只有报纸了。
' y* x$ @% ?# g, A/ n) h
作者: 陈阳康    时间: 2010-11-26 13:20
地板地板
' L  c3 V' s' f
作者: 陈阳康    时间: 2010-11-26 13:20
接下来是什么?
作者: 陈阳康    时间: 2010-11-26 13:21
进入正题。MATLAB我不太会额~~~




欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5