数学建模社区-数学中国
标题:
谁能用MATLAB模拟高斯气体分布
[打印本页]
作者:
核科学
时间:
2010-5-5 18:11
标题:
谁能用MATLAB模拟高斯气体分布
高斯气体分布
1 Y0 V, h- T: o5 ?) j5 j
数值模拟
. ^/ y, C3 U5 _- W* Z) Q
作出动态分布图
作者:
923050019
时间:
2010-8-6 22:56
function y=shuzhimoni
" O' R/ |+ E# Y0 b9 H
%syms x y z t;
- a, [# f5 m0 X C" r
%Q=10;%初始浓度定为10mg/m^3;
& I4 X% Y4 ~) v' \
%1.32*10^-5=1.32*10^-5;%扩散系数
! f" t2 \( O; F. E+ C
%C=Q/(4*pi.*1.32*10^-5.*t)^1.5*exp(-(x.^2+y.^2+z.^2)./4.*1.32*10^-5.*t);
8 S$ Z" [: M( Q/ w
%统一并申明一下单位
1 ]( Y; v' y" x& t3 I
%C (mg/m^3)
c2 e* _1 L, e* `! L
%Q (mg/m^3)
' V" d) c2 L0 p8 K l
%1.32*10^-5 (m^2/s)
/ h& B, f+ N' k: }
%t (s)
9 }% h+ s' v2 v" Q: x
%x,y,z (m)
- j2 I7 P3 W8 B6 [: a
%申明一下扩散的条件如下：
' P* ^; y6 Z& r3 s$ a
%1.模拟的球形空间的直径D=3.5m 在x方向安装15个，在y方向安装3个
6 x- Y+ r9 E+ W& G* J F& W
%2.当经过时间10s时，可视为扩散均匀，其中时间可以通过一次实验或者相关文献查得。
9 |3 n7 |8 t+ M
randnum=unifrnd(1,10,10,1);%产生一组符合均匀分布的时间
) L. v, W+ z( z+ r
for i=1:4
! w% W5 o: e4 _, ~7 o1 c
[x,y,z] = meshgrid(-1.75:.2:1.75,-1.75:.2:1.75,-1.75:.2:1.75);
9 j. l6 E* W* D4 A P* ?; \; D& W: L
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));
: W" ] a9 W8 p* t5 s
xslice=[-0.25,1.75];
5 P& l/ O k" L2 Y3 j* r2 P. L
yslice=[-0.25,1.25];
6 ?8 ?3 N {( I# g
zslice=0.25;
4 L# X8 b) K; x* ^ ?) j
subplot(2,2,i);
8 ~; W6 O! f# F! d* a( q
slice(x,y,z,C,xslice,yslice,zslice);
3 v4 H4 {* @! C; x. \' E/ i
xlabel('x坐标');
* K; q' I+ E% }2 }+ X0 ^% ^
ylabel('y坐标');
( }5 A* c2 H! c. W/ G/ U
zlabel('z坐标');
4 k+ w$ {, _9 ~3 V; z: h
title(['当时间为' num2str(randnum(i)) '时的浓度分布情况']);
. n- ~1 q9 }' O$ V# t5 |
colormap hsv
& G k2 H1 \) ?6 r! J9 c' A8 `
colorbar
5 y; b6 d+ \$ n2 C( {3 D
end
+ J/ y0 y) x% I, s& K
作者:
13ansys
时间:
2014-8-11 14:24
看不懂什么意思 啊
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5