数学建模社区-数学中国
标题:
急求Matlab程序分析解释
[打印本页]
作者:
范婷
时间:
2010-11-25 15:45
标题:
急求Matlab程序分析解释
有哪位高手可以帮忙分析一下这个程序！！！！！！！！！
% L* Q3 M, l# g: o6 R
r=1;%点光源的能量值，可以修改。%
/ H7 @6 T# m- O* L) s
A=zeros(81,81); %方程的系数矩阵。%
4 A5 n3 x5 e1 q( l
for i=1:81
( G$ T+ f4 y, K# `
if rem(i-1,9)~=0
9 C& `3 A3 d: p7 P. L" O! E& e
A(i,i-1)=1/4;
9 y- d+ M% N; I- z8 T S$ k) s
end
( U1 G# ^3 m8 e) o$ V
A(i,i)=-1;
8 n* l) N) i2 |5 k; }
if rem(i,9)~=0
& L( F6 {- F% X t3 C6 d f
A(i,i+1)=1/4;
$ c( t$ g5 i2 G5 ^" c( i9 r/ \; n" Y
end
0 i: X2 ~: p2 @8 y% ?6 c- d
if (i-9)>0
7 i0 A. l/ H; I" @: D( V1 j
A(i,i-9)=1/4;
$ p) X% w7 _6 P" ?2 o5 p& Z; K3 G! v
end
Y/ {6 N5 \5 h R/ e+ z7 v# D5 g
if (i+9)<82
, p! O- C8 m2 R9 @/ G
A(i,i+9)=1/4;
+ r4 V ^1 C3 Q& S
end
9 Y$ e6 e4 d0 [+ c0 b
end
! W9 O c! J9 ?; b
B=zeros(81,1);
4 {; t: j, J3 ~9 k$ r9 ^! b
B(6)=-r/4;
0 ]! y5 y6 J4 I: x/ x
X=inv(A)*B;
) h, e$ E# @% b+ N1 d
D=reshape(X,9,9);
) _: y% r1 H) g. ~
D=D';
% B4 {4 L- b5 A5 F ~
E=D;
) e; g# x3 t; z$ B
x=zeros(1,11);
- V# _9 R W3 M; w
y=zeros(9,1);
) E7 P' x$ K) \4 N& x3 Q; Y. D
D=[y,D,y];
2 b" m* U1 o9 E
D=[x;D;x];
4 H+ |7 t& r- ]. l( ^
D(1,6)=1;
8 u# F+ X% y9 I& U3 v$ y: P! p
xt=0:10;
2 p& y+ @ E/ v3 D4 y5 F3 \
[xt,yt]=meshgrid(xt);
- ]& Y4 Z( }8 \3 D! K9 M" E
zt=D;
- P* K3 s$ O( i' F6 M3 ~
surf(xt,yt,zt);
& F; J/ U; j7 |0 i* s- V% s7 q* ^4 ~
D %矩阵D上的值对应每个节点上的能量值。%
/ M5 W( A# M2 w! H' O$ z1 V
作者:
hy282694
时间:
2010-11-26 08:20
你要分析啥？带着你看懂源代码？还是。。。？
作者:
fhfhfh555
时间:
2010-11-26 09:41
1111111111111111111
作者:
陈阳康
时间:
2010-11-26 13:19
只有报纸坐了。
作者:
陈阳康
时间:
2010-11-26 13:20
只有报纸了。
9 G$ ]( }0 D# H* b7 z' k- Y7 C, w; F
作者:
陈阳康
时间:
2010-11-26 13:20
地板地板
5 E% f( J/ \# P+ w6 I5 ?0 w7 f/ @! h
作者:
陈阳康
时间:
2010-11-26 13:20
接下来是什么？
作者:
陈阳康
时间:
2010-11-26 13:21
进入正题。MATLAB我不太会额~~~
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5