数学建模社区-数学中国
标题:
急求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:81
0 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
end
6 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- z
end
7 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& e
D=reshape(X,9,9);
# @9 C' V8 Z# {& r# e, W% H
D=D';
" ^* ]( t/ }& Z( x( v: y
E=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 n
D=[x;D;x];
$ U1 ?; u! l: E
D(1,6)=1;
* B) I* L4 V8 D- F d
xt=0:10;
8 s& o9 D% U; F, M! x
[xt,yt]=meshgrid(xt);
/ T/ i0 F5 x, E
zt=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