数学建模社区-数学中国
标题:
急求Matlab程序分析解释
[打印本页]
作者:
范婷
时间:
2010-11-25 15:45
标题:
急求Matlab程序分析解释
有哪位高手可以帮忙分析一下这个程序！！！！！！！！！
% P: p' A* S! K2 p
r=1;%点光源的能量值，可以修改。%
$ E! h8 k; w: [# L \
A=zeros(81,81); %方程的系数矩阵。%
% `; m1 Q4 w2 b% _" q
for i=1:81
9 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
end
8 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
end
2 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% c
X=inv(A)*B;
_5 g; }5 Y3 u8 [1 K! C
D=reshape(X,9,9);
0 Y- Y7 x$ [9 V; x/ _# a
D=D';
' c# c1 T( r7 s" Z5 ~
E=D;
% \: j% }' O/ b" d: R. `
x=zeros(1,11);
6 H# e0 W# E1 a. w
y=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, N
D(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