数学建模社区-数学中国
标题:
谁能用MATLAB模拟高斯气体分布
[打印本页]
作者:
核科学
时间:
2010-5-5 18:11
标题:
谁能用MATLAB模拟高斯气体分布
高斯气体分布
$ p3 k, u: h1 @, w1 F
数值模拟
/ l1 e( A( [! H1 j% H+ u: m$ b
作出动态分布图
作者:
923050019
时间:
2010-8-6 22:56
function y=shuzhimoni
6 ?$ ?; e/ I' s
%syms x y z t;
: z) o& ?3 H2 R: \' m4 u, I n
%Q=10;%初始浓度定为10mg/m^3;
* q6 o( ?6 I2 }
%1.32*10^-5=1.32*10^-5;%扩散系数
) p n) T/ a$ K) s) \8 B" d
%C=Q/(4*pi.*1.32*10^-5.*t)^1.5*exp(-(x.^2+y.^2+z.^2)./4.*1.32*10^-5.*t);
# i0 r, u2 s3 e+ N4 r4 {
%统一并申明一下单位
% k& Y$ j6 {5 u: _) s
%C (mg/m^3)
/ y' C4 g- {( Y8 d3 L( g& f# T
%Q (mg/m^3)
* C- }3 t! o. d
%1.32*10^-5 (m^2/s)
; K1 H0 m0 X- A9 i6 X6 [& R
%t (s)
) ~8 s/ l; l* @% S
%x,y,z (m)
. S# f9 `* h+ c7 t9 V3 {
%申明一下扩散的条件如下：
: z6 G4 J' W6 D8 F! C% g
%1.模拟的球形空间的直径D=3.5m 在x方向安装15个，在y方向安装3个
: i: b4 e8 {4 R" s8 e; M
%2.当经过时间10s时，可视为扩散均匀，其中时间可以通过一次实验或者相关文献查得。
0 y) H- T3 [6 p1 c8 _% B
randnum=unifrnd(1,10,10,1);%产生一组符合均匀分布的时间
. h) w, |) w8 L
for i=1:4
/ Y7 \2 Z4 f4 H: S
[x,y,z] = meshgrid(-1.75:.2:1.75,-1.75:.2:1.75,-1.75:.2:1.75);
: e, V; @6 C* T8 E1 l. \( T* k
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' E/ d( R$ s3 u+ F
xslice=[-0.25,1.75];
4 U" U) \& z1 B6 F
yslice=[-0.25,1.25];
0 R! `+ f7 m6 l: W! `
zslice=0.25;
: o9 N; R2 \ `8 G/ K) i# A& x
subplot(2,2,i);
2 j, \9 M1 f b% E ?* h. T
slice(x,y,z,C,xslice,yslice,zslice);
* t( C3 i! a3 V* y0 I9 Z
xlabel('x坐标');
' u9 m5 P8 i* K3 ?
ylabel('y坐标');
" O8 B* R" v: T* ]
zlabel('z坐标');
* L! Q& y, Q1 l! `- j8 y. X
title(['当时间为' num2str(randnum(i)) '时的浓度分布情况']);
, W, d; D5 `5 ]% `8 j8 |3 F
colormap hsv
0 C" F- h0 h3 X5 e0 e4 w- E
colorbar
; \/ x; T) I. @' T1 q1 N+ o
end
) A* |( |8 D' v! I. N
作者:
13ansys
时间:
2014-8-11 14:24
看不懂什么意思 啊
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5