数学建模社区-数学中国

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

作者: 范婷    时间: 2010-11-25 15:45
标题: 急求Matlab程序分析解释
有哪位高手可以帮忙分析一下这个程序!!!!!!!!!- ^; ?' i( d1 J- i
r=1;%点光源的能量值,可以修改。%# G. N8 T; V1 Q& k2 g8 Q5 o
A=zeros(81,81);  %方程的系数矩阵。%
8 ~' a: B% p8 F# \: F$ P, H# J" t* e% _for i=1:810 L+ R3 O% {, X. [3 K( Z1 d
    if rem(i-1,9)~=0
4 N  g, S% n% _+ D8 E5 H/ ~  N6 z        A(i,i-1)=1/4;, C3 d8 ~; w2 o% I- X+ y" m; g5 L
    end
6 C% @- Y) x: U    A(i,i)=-1;" f7 T! C2 Z/ P5 l0 x2 b% b
    if rem(i,9)~=0
# D$ V( T2 D9 I  Z5 Y$ Z: p        A(i,i+1)=1/4;
% V- l2 }% u- {$ R0 Y2 k    end6 O3 O# k0 t  z2 E
    if (i-9)>0
8 x- W+ G1 ^" ^  f3 i# v        A(i,i-9)=1/4;
6 P( a# T: _, l' `6 D0 b    end
/ t, H' ], K( J, X: v    if (i+9)<82
( |( K( p0 I  d- S9 d- w        A(i,i+9)=1/4;- E' L# N. I& @7 P! R  U- x5 {: `
    end
  K0 d9 h/ A( O- B- zend7 v; S2 X( [/ s; s" U0 ?, Q
B=zeros(81,1);) i1 u8 u  L2 H3 H! p
B(6)=-r/4;   * z* e  W7 T2 ?+ w  Q6 ~5 s9 o
X=inv(A)*B;
& V, i* z0 M3 S& eD=reshape(X,9,9);# @9 C' V8 Z# {& r# e, W% H
D=D';
" ^* ]( t/ }& Z( x( v: yE=D;/ f" l5 b0 @0 y- n
x=zeros(1,11);, A1 w' O# Z( ]: S
y=zeros(9,1);5 U& J) g# T+ X7 Q( X' s& ?8 D
D=[y,D,y];
7 F+ V% k8 N. c: j8 t9 nD=[x;D;x];$ U1 ?; u! l: E
D(1,6)=1;
* B) I* L4 V8 D- F  dxt=0:10;8 s& o9 D% U; F, M! x
[xt,yt]=meshgrid(xt);
/ T/ i0 F5 x, Ezt=D;
$ n0 R/ z, n$ j9 `surf(xt,yt,zt);& ]3 x$ C8 c& W, Y) A! M+ m
D      %矩阵D上的值对应每个节点上的能量值。%
3 l  A$ m) l. @) _) q
作者: hy282694    时间: 2010-11-26 08:20
你要分析啥?带着你看懂源代码?还是。。。?
作者: fhfhfh555    时间: 2010-11-26 09:41
1111111111111111111
作者: 陈阳康    时间: 2010-11-26 13:19
只有报纸坐了。
作者: 陈阳康    时间: 2010-11-26 13:20
只有报纸了。$ c' F- l  I$ g" a1 H

作者: 陈阳康    时间: 2010-11-26 13:20
地板地板
  ~/ A% r" @2 J1 z
作者: 陈阳康    时间: 2010-11-26 13:20
接下来是什么?
作者: 陈阳康    时间: 2010-11-26 13:21
进入正题。MATLAB我不太会额~~~




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