数学建模社区-数学中国
标题:
急求Matlab程序分析解释
[打印本页]
作者:
范婷
时间:
2010-11-25 15:45
标题:
急求Matlab程序分析解释
有哪位高手可以帮忙分析一下这个程序！！！！！！！！！
2 [/ B1 J- V# A" W! U' z
r=1;%点光源的能量值，可以修改。%
8 u/ Z5 K% L* l3 G L
A=zeros(81,81); %方程的系数矩阵。%
& W4 f2 B1 K. @# q0 u6 c0 V. n
for i=1:81
1 j" [ } Q- [7 j% U* c
if rem(i-1,9)~=0
* a+ R$ w: \) K% _
A(i,i-1)=1/4;
5 C$ H* H$ g) F+ m1 U: w
end
" D9 ?! l- X V# Q, x5 X
A(i,i)=-1;
& F) y$ ^; K9 P* F0 B2 g
if rem(i,9)~=0
& A' M" L! w0 t3 H
A(i,i+1)=1/4;
( Y0 K9 W. v; p c. K0 ]
end
. ^1 n8 Z; F. d# p3 A
if (i-9)>0
$ u& {! {* D3 [. L
A(i,i-9)=1/4;
# G: N6 q& s. d9 s( p' r
end
$ _$ b% @! ^! [$ T. b
if (i+9)<82
3 M3 j* l& g+ K
A(i,i+9)=1/4;
& Y7 B8 B! l) `8 G; C
end
& o, Y6 P' |, ^: G! L2 _& ^
end
; Q% ?, H* p8 c
B=zeros(81,1);
7 g, w" w1 }3 h+ q4 ~
B(6)=-r/4;
5 J$ p& p% g3 ]$ m6 V
X=inv(A)*B;
" G e/ u* V9 c" @
D=reshape(X,9,9);
3 i; b- |$ \0 _# w
D=D';
& _6 g! _" t' l9 }6 r" K- C5 r/ J
E=D;
5 n4 g/ e# O2 Y; ]: x5 m' c3 Q
x=zeros(1,11);
! v, b9 K+ ^% p6 Q
y=zeros(9,1);
; `+ j- k6 f3 _1 X' S! t
D=[y,D,y];
. k( }* l8 b: ]# d; [
D=[x;D;x];
9 o* W( U* I$ Q" v0 ]1 g
D(1,6)=1;
/ q" f/ H3 W" l3 H
xt=0:10;
4 O0 D+ L# J( ?9 i
[xt,yt]=meshgrid(xt);
, J# \0 Q; b+ W6 X9 \( o
zt=D;
3 U& S4 B# x3 O6 `/ g3 l
surf(xt,yt,zt);
5 h/ p5 j1 H9 J& _
D %矩阵D上的值对应每个节点上的能量值。%
# e3 x/ ?5 d2 _6 Y0 @6 _. [( K! q# q# b
作者:
hy282694
时间:
2010-11-26 08:20
你要分析啥？带着你看懂源代码？还是。。。？
作者:
fhfhfh555
时间:
2010-11-26 09:41
1111111111111111111
作者:
陈阳康
时间:
2010-11-26 13:19
只有报纸坐了。
作者:
陈阳康
时间:
2010-11-26 13:20
只有报纸了。
( ]- H- J; Q1 G& K9 @% g
作者:
陈阳康
时间:
2010-11-26 13:20
地板地板
% }5 M! I& q# j0 j. N! {. Y
作者:
陈阳康
时间:
2010-11-26 13:20
接下来是什么？
作者:
陈阳康
时间:
2010-11-26 13:21
进入正题。MATLAB我不太会额~~~
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5