- 在线时间
- 24 小时
- 最后登录
- 2013-5-6
- 注册时间
- 2009-8-13
- 听众数
- 8
- 收听数
- 0
- 能力
- 0 分
- 体力
- 1988 点
- 威望
- 0 点
- 阅读权限
- 50
- 积分
- 769
- 相册
- 1
- 日志
- 25
- 记录
- 41
- 帖子
- 290
- 主题
- 22
- 精华
- 0
- 分享
- 3
- 好友
- 101
升级   42.25% TA的每日心情 | 开心 2013-5-6 20:58 |
|---|
签到天数: 3 天 [LV.2]偶尔看看I
群组: 数学建模 群组: LINGO 群组: 代数与数论 群组: C 语言讨论组 |
2#
发表于 2010-8-6 22:56
|只看该作者
|
|邮箱已经成功绑定
function y=shuzhimoni
* y. A/ V U9 f( w" A%syms x y z t;
& z& B5 n6 h1 ?( d+ u+ j%Q=10;%初始浓度定为10mg/m^3;/ m8 k+ K2 B6 F+ P: D
%1.32*10^-5=1.32*10^-5;%扩散系数1 @7 y6 z3 j8 \+ q
%C=Q/(4*pi.*1.32*10^-5.*t)^1.5*exp(-(x.^2+y.^2+z.^2)./4.*1.32*10^-5.*t);7 O6 T; e& L. g( G( f; N
%统一并申明一下单位( c' }8 ]8 @% O4 j$ t) p" [! C1 U8 T
%C (mg/m^3)
6 U2 \/ {+ z" f( s! {, h( J%Q (mg/m^3)
% R* ~; }" D, b%1.32*10^-5 (m^2/s)
. k4 L; z9 J0 F6 U6 o0 h2 Y%t (s)8 h9 N9 _7 }( p
%x,y,z (m)
* u+ e" f# o: H%申明一下扩散的条件如下:
# ?5 J. I* X4 @" S1 ?%1.模拟的球形空间的直径D=3.5m 在x方向安装15个,在y方向安装3个
- G) _, P# b4 ^! J" E%2.当经过时间10s时,可视为扩散均匀,其中时间可以通过一次实验或者相关文献查得。" S, P/ R9 w7 H: ?- s
randnum=unifrnd(1,10,10,1);%产生一组符合均匀分布的时间
6 r2 g2 R8 R7 r4 S) y7 K1 r. ufor i=1:4
" o1 E( o! g3 n2 M- x' V5 }) q [x,y,z] = meshgrid(-1.75:.2:1.75,-1.75:.2:1.75,-1.75:.2:1.75);
( \, J* @$ \- b( X C=10/(4*pi.*1.32*10^-5.*randnum(i))^1.5*exp(-(x.^2+y.^2+z.^2)./4.*1.32*10^-5.*randnum(i));- R" h6 `) X% J2 m! O- q
xslice=[-0.25,1.75];
, _- j, r1 h8 ?: o yslice=[-0.25,1.25];
: _3 T( I8 l# P+ z zslice=0.25;4 Y/ {) q% H# A$ f
subplot(2,2,i);1 G1 D. t+ F/ l7 ]" W9 W; h
slice(x,y,z,C,xslice,yslice,zslice);
9 h# j$ p4 |8 _, U* L- ^& Q xlabel('x坐标');
4 L8 M7 F3 g# g2 p6 z. l ylabel('y坐标');
+ m4 j" V$ L! k) K1 ? zlabel('z坐标');
& J0 @4 z- _; j! V# ?7 |* L8 z6 L" ? title(['当时间为' num2str(randnum(i)) '时的浓度分布情况']); $ Q& T8 N+ p6 E, a
colormap hsv7 N$ E1 |4 k7 a+ m
colorbar2 t1 n7 X% }8 k
end
; b3 J9 n& D8 H f) E- {/ b7 W8 E- c |
|