数学建模社区-数学中国
标题:
急求Matlab程序分析解释
[打印本页]
作者:
范婷
时间:
2010-11-25 15:45
标题:
急求Matlab程序分析解释
有哪位高手可以帮忙分析一下这个程序！！！！！！！！！
; y3 ^7 B" M% r
r=1;%点光源的能量值，可以修改。%
+ P$ Y2 T4 E/ W* n
A=zeros(81,81); %方程的系数矩阵。%
& U! F) ? l8 t" F5 O
for i=1:81
% V3 S/ U. t, B2 t1 G! s" |
if rem(i-1,9)~=0
1 e/ @! c) _* ]- h4 D9 N2 b, {5 u$ X
A(i,i-1)=1/4;
" [2 Y6 L0 _6 `0 v. a
end
5 T8 R: o) } R& g: i
A(i,i)=-1;
* j+ P0 t+ k, I* S
if rem(i,9)~=0
, P& H, O+ g& o# Q, e2 q; {) ^
A(i,i+1)=1/4;
& N% A; w) B5 k% z6 v9 S2 o% q
end
; i" p! C% U6 k0 B' P
if (i-9)>0
* B8 C5 J& N5 }/ H& h( E
A(i,i-9)=1/4;
1 G) O0 D7 `2 L4 s( K% G
end
6 t5 h- V' h/ G8 L
if (i+9)<82
6 X! a" R5 Z0 x3 o1 y& Y5 J0 O0 I
A(i,i+9)=1/4;
/ y1 \4 r* |6 X5 |% s" a6 [/ @
end
: m& @* F M% c* n$ M% ~, P5 B
end
1 b5 c9 M, Q. m$ d
B=zeros(81,1);
! o9 R X' ~$ G4 M8 ]% y; R8 U
B(6)=-r/4;
; K& L2 h; C: c, w6 W, D4 Y
X=inv(A)*B;
1 M& ^" p2 d: a5 v2 F5 M) O- r
D=reshape(X,9,9);
( h4 U! ~% v3 j3 d% B8 |$ C
D=D';
1 r3 ]& c- {$ W( p: c
E=D;
) ] z1 |8 ~8 Y
x=zeros(1,11);
% ^: O8 d/ v) U* Y. @2 ^
y=zeros(9,1);
- R6 `9 T S- N |/ M1 M' O. i% ^
D=[y,D,y];
9 g* X; R$ s( e8 A
D=[x;D;x];
: m; p: I+ U/ u+ H4 E0 b
D(1,6)=1;
8 @( K9 @( x) f& K, N
xt=0:10;
; _2 L2 l+ L0 g5 M# o+ B1 B0 p$ g
[xt,yt]=meshgrid(xt);
& [0 {( }8 x. \) i
zt=D;
; s1 k% `3 m/ ^ `, A, m
surf(xt,yt,zt);
" o. O D+ L5 N" `) m
D %矩阵D上的值对应每个节点上的能量值。%
) V9 p6 F3 O/ S+ Z7 |
作者:
hy282694
时间:
2010-11-26 08:20
你要分析啥？带着你看懂源代码？还是。。。？
作者:
fhfhfh555
时间:
2010-11-26 09:41
1111111111111111111
作者:
陈阳康
时间:
2010-11-26 13:19
只有报纸坐了。
作者:
陈阳康
时间:
2010-11-26 13:20
只有报纸了。
5 u) ~: B* ]9 Y, w4 t) b% P
作者:
陈阳康
时间:
2010-11-26 13:20
地板地板
' L! a- Y2 \7 y
作者:
陈阳康
时间:
2010-11-26 13:20
接下来是什么？
作者:
陈阳康
时间:
2010-11-26 13:21
进入正题。MATLAB我不太会额~~~
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5