- 在线时间
- 24 小时
- 最后登录
- 2013-5-6
- 注册时间
- 2009-8-13
- 听众数
- 8
- 收听数
- 0
- 能力
- 0 分
- 体力
- 1988 点
- 威望
- 0 点
- 阅读权限
- 50
- 积分
- 769
- 相册
- 1
- 日志
- 25
- 记录
- 41
- 帖子
- 292
- 主题
- 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
- O" }2 C& }: B# p: _3 e%syms x y z t;5 k) A X9 p1 K& ^ i5 U
%Q=10;%初始浓度定为10mg/m^3;
- V1 Z% |/ h3 a, a) Q/ ^%1.32*10^-5=1.32*10^-5;%扩散系数% C. W9 k/ Q* d$ Y; `
%C=Q/(4*pi.*1.32*10^-5.*t)^1.5*exp(-(x.^2+y.^2+z.^2)./4.*1.32*10^-5.*t);
$ j7 G H/ \% f4 L%统一并申明一下单位! d5 \9 v* `6 |6 X$ I4 H3 F- v/ }
%C (mg/m^3)
2 r% R+ C# a. S, w%Q (mg/m^3)$ ?: [. w3 b$ O8 x5 X. h* V
%1.32*10^-5 (m^2/s)
3 O5 n G) p# q2 G, e" x%t (s)
! q n6 `. L) Y0 q3 m5 _%x,y,z (m)# y2 R6 i$ M' V; a+ v+ x9 L2 H% T
%申明一下扩散的条件如下:4 B. P- N4 ^* r( Y$ W& ?+ E
%1.模拟的球形空间的直径D=3.5m 在x方向安装15个,在y方向安装3个
5 o" V( k1 G }" C0 ]' E0 H%2.当经过时间10s时,可视为扩散均匀,其中时间可以通过一次实验或者相关文献查得。' x) R1 f0 V5 r& K
randnum=unifrnd(1,10,10,1);%产生一组符合均匀分布的时间
6 ?: u3 W' V+ {' ]* efor i=1:4 a2 [, N0 A5 P4 \* B* Q2 U/ e
[x,y,z] = meshgrid(-1.75:.2:1.75,-1.75:.2:1.75,-1.75:.2:1.75);
+ r) T# ]8 N4 W% p 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));
& ]' V" o8 q. o) W xslice=[-0.25,1.75];
; O0 o4 r* p3 o" ~4 z0 r0 l yslice=[-0.25,1.25];/ U- T, H" d, Q
zslice=0.25;7 d2 |) {6 ?0 {& Q) b# P p& G @
subplot(2,2,i);
( f0 G! C+ K, U1 ^, L slice(x,y,z,C,xslice,yslice,zslice);
1 f' j) E0 ^4 V1 ^7 }7 Q xlabel('x坐标');
8 S2 V- j& J8 F5 C ylabel('y坐标');
* U+ ]* q0 I" [4 C$ S zlabel('z坐标');
, q* u; l0 ?) k* u. n title(['当时间为' num2str(randnum(i)) '时的浓度分布情况']); 4 ^5 o; W0 N# M. X4 U
colormap hsv
+ I" B1 [- X" j+ Z" K0 p3 R# d colorbar
3 f: a0 w5 n7 \; V) [7 Z& R& Send7 D o2 s6 ^5 x9 C
|
|