- 在线时间
- 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
5 N7 ]8 H \ p- h' Y: }0 Z ^%syms x y z t;
3 _) A0 d: f. X' L%Q=10;%初始浓度定为10mg/m^3;, @4 @9 Y, q9 N X" f0 |+ E r6 r
%1.32*10^-5=1.32*10^-5;%扩散系数
6 |, t% m% z# u W( r! c, r% 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);
& w( P- }4 g* p( A%统一并申明一下单位
1 b" E L' l! p; U%C (mg/m^3)/ _' K& r+ N6 D+ x
%Q (mg/m^3)
4 M, e7 j Q; S: u5 S! `4 s%1.32*10^-5 (m^2/s)
, L- {* Z1 P% Z! [' M8 P6 h%t (s)4 N! L7 E2 {5 c% @& f
%x,y,z (m)
( @' U& B( c3 [7 g/ _! o# V6 Z%申明一下扩散的条件如下:
; \% I3 \2 g6 H# ~; c) j8 e%1.模拟的球形空间的直径D=3.5m 在x方向安装15个,在y方向安装3个8 Z: a- A* A6 |& t; L
%2.当经过时间10s时,可视为扩散均匀,其中时间可以通过一次实验或者相关文献查得。5 U9 J9 e2 g" m; a# s' K
randnum=unifrnd(1,10,10,1);%产生一组符合均匀分布的时间% i8 d1 g8 P% V8 U
for i=1:45 i: D; M5 u3 v) {
[x,y,z] = meshgrid(-1.75:.2:1.75,-1.75:.2:1.75,-1.75:.2:1.75);* g# ^9 L7 d0 J% n9 B
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));
$ O1 m7 Z' k1 x3 ^- W xslice=[-0.25,1.75];0 f) R3 n6 f, r8 [! D" [
yslice=[-0.25,1.25];
! o r) o1 d) T% j# u( T" Q zslice=0.25;* K7 y: g E& X7 Q9 A- D
subplot(2,2,i);
, u$ v6 y A) j7 s# @ slice(x,y,z,C,xslice,yslice,zslice);6 x# Z2 o; U0 {: }1 G" X
xlabel('x坐标');
" Y5 I/ r; r. T7 r4 V( a ylabel('y坐标');' [8 i; b$ F; v5 h: B8 g* S
zlabel('z坐标');
8 p4 A9 ?7 V$ S, ] title(['当时间为' num2str(randnum(i)) '时的浓度分布情况']);
0 }4 { _$ |6 r" n7 \$ K! { colormap hsv5 L3 j' j0 f0 w4 f, k' |" b; b% S2 u
colorbar1 K) F0 s( M& w, y& |* x
end7 |/ b! |: p. a
|
|