- 在线时间
- 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$ P6 u- b7 A9 L& ?
%syms x y z t;* e8 k. S- P" n u) X, n, v
%Q=10;%初始浓度定为10mg/m^3;% a6 r' e: N* Z3 K7 o# j
%1.32*10^-5=1.32*10^-5;%扩散系数. B8 ^$ J, r+ e$ d" S
%C=Q/(4*pi.*1.32*10^-5.*t)^1.5*exp(-(x.^2+y.^2+z.^2)./4.*1.32*10^-5.*t);2 E. C$ P: I" l$ n; |/ I
%统一并申明一下单位
, Z' _- f" M7 N( e" F7 z( F0 ^%C (mg/m^3): u5 l5 ^& X9 ?# O- F
%Q (mg/m^3)7 k0 d0 M+ u! J; _% Q
%1.32*10^-5 (m^2/s)) r' t1 U* m4 P/ u
%t (s)& Q! e% v, H+ u3 Q/ V5 i7 V# h) x& N/ Z
%x,y,z (m); p8 q1 B! ~8 N
%申明一下扩散的条件如下:2 O4 [( z9 n* ?! U! \5 T. t" x
%1.模拟的球形空间的直径D=3.5m 在x方向安装15个,在y方向安装3个
1 g/ g2 z) t: _, L% g%2.当经过时间10s时,可视为扩散均匀,其中时间可以通过一次实验或者相关文献查得。
3 D- e9 j& C" V) B4 Y: r; _2 vrandnum=unifrnd(1,10,10,1);%产生一组符合均匀分布的时间
! j, Y0 e' q) o& Gfor i=1:4
2 C- d5 S8 I4 D2 \3 l$ } [x,y,z] = meshgrid(-1.75:.2:1.75,-1.75:.2:1.75,-1.75:.2:1.75);+ y7 U$ r& d C2 A& a
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));
' Q8 a4 B# i5 s" A" P. j: p xslice=[-0.25,1.75];
+ c# h$ d2 {9 I( E yslice=[-0.25,1.25];* U3 {9 o7 `. E2 T
zslice=0.25;! b( c; H+ g' J7 H. J4 M( S
subplot(2,2,i);4 U! E$ V; ]7 C2 e% o
slice(x,y,z,C,xslice,yslice,zslice);6 n; t1 ?% V- U4 Y1 t2 p! o: |
xlabel('x坐标');2 q% f; ?2 a* A# J
ylabel('y坐标');1 x1 J; }* n( h) |8 e7 h. w3 }
zlabel('z坐标');# Z( X+ I- U% J5 I
title(['当时间为' num2str(randnum(i)) '时的浓度分布情况']); ! H& u+ T" ^# D! e2 }
colormap hsv; u( M9 R- e* h8 f
colorbar- [+ t2 m4 [0 `3 e3 A
end- i! ?8 f# _+ O' F9 G6 U5 f9 g0 H
|
|