数学建模社区-数学中国
标题:
谁能用MATLAB模拟高斯气体分布
[打印本页]
作者:
核科学
时间:
2010-5-5 18:11
标题:
谁能用MATLAB模拟高斯气体分布
高斯气体分布
' D1 q. b5 \# Z. C: ^
数值模拟
; m% b7 Q+ j- ~
作出动态分布图
作者:
923050019
时间:
2010-8-6 22:56
function y=shuzhimoni
1 j; w: u& h6 v! B. w) h. v. b
%syms x y z t;
, l7 w1 f2 X& p3 ]/ c! s
%Q=10;%初始浓度定为10mg/m^3;
o. u6 @# b9 m) r, N' {
%1.32*10^-5=1.32*10^-5;%扩散系数
( F0 [* k* ^9 g( q2 N' y, }' B, w
%C=Q/(4*pi.*1.32*10^-5.*t)^1.5*exp(-(x.^2+y.^2+z.^2)./4.*1.32*10^-5.*t);
( ~7 x- t' O1 p' e
%统一并申明一下单位
& b. u/ U. w7 O8 Y2 e+ [
%C (mg/m^3)
0 W/ D+ N# [4 D+ z0 ^/ S
%Q (mg/m^3)
* P; A2 P3 w2 W: U
%1.32*10^-5 (m^2/s)
, I8 t% \+ I8 W, N
%t (s)
1 a4 Y ~. }8 I
%x,y,z (m)
- V% d$ s n* O& X$ t( x9 F3 J
%申明一下扩散的条件如下：
: Y' {8 ]4 Q/ N8 u2 m4 f7 ]! R
%1.模拟的球形空间的直径D=3.5m 在x方向安装15个，在y方向安装3个
6 l8 D6 a/ y m6 s
%2.当经过时间10s时，可视为扩散均匀，其中时间可以通过一次实验或者相关文献查得。
2 X9 @8 X$ \2 ]: b
randnum=unifrnd(1,10,10,1);%产生一组符合均匀分布的时间
" r! L K( i/ a- o
for i=1:4
( R; |. H" }$ G5 S+ I
[x,y,z] = meshgrid(-1.75:.2:1.75,-1.75:.2:1.75,-1.75:.2:1.75);
9 x( R$ ]* J- j! m8 L) m4 y9 Q% ]
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));
" T$ A4 T* W- C+ o1 A+ O
xslice=[-0.25,1.75];
" p3 ]. t( g7 C% D
yslice=[-0.25,1.25];
/ _/ c% I+ u5 Q- ]$ d9 o. ?1 A7 P
zslice=0.25;
' d+ B6 b/ N" l# i
subplot(2,2,i);
' x+ e& A& {% i7 M1 u0 Z+ `& G
slice(x,y,z,C,xslice,yslice,zslice);
/ |6 d1 P ?; N
xlabel('x坐标');
2 c& Q/ [6 Y3 x4 g( o" u
ylabel('y坐标');
; s$ P% ~+ O; I0 Z
zlabel('z坐标');
0 b4 k, X' U8 D" T% c8 ~3 H
title(['当时间为' num2str(randnum(i)) '时的浓度分布情况']);
% E$ ]% U3 |( R, { J
colormap hsv
# k; X& i0 Z0 s1 k5 U/ x0 `9 x
colorbar
& i; T- Z: U4 a& @4 Q3 V \5 j; X3 S
end
' v1 q1 |/ v0 D- T( B2 U
作者:
13ansys
时间:
2014-8-11 14:24
看不懂什么意思 啊
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5