数学建模社区-数学中国»论坛 › 【科学软件论坛】数学建模 › MATLAB论坛 › Matlab simulink建模与仿真 › 谁能用MATLAB模拟高斯气体分布

查看: 4712|回复: 2

谁能用MATLAB模拟高斯气体分布

[复制链接]

字体大小: 正常放大

核科学

1 主题	3 听众	19 积分

升级 14.74%

该用户从未签到

自我介绍: 我我我

窝窝

Ooo

电梯直达

1^#

发表于 2010-5-5 18:11 |只看该作者 |倒序浏览

|招呼Ta 关注Ta

高斯气体分布
数值模拟
作出动态分布图

zan

转播0 淘帖0 分享0 收藏0 支持0 反对0 微信

使用道具举报

923050019

22 主题	8 听众	769 积分

升级 42.25%

TA的每日心情

	开心 2013-5-6 20:58

签到天数: 3 天

[LV.2]偶尔看看I

群组: 数学建模

群组: LINGO

群组: 代数与数论

群组: C 语言讨论组

2^#

发表于 2010-8-6 22:56 |只看该作者 |招呼Ta 关注Ta |邮箱已经成功绑定

function y=shuzhimoni
%syms x y z t;
%Q=10;%初始浓度定为10mg/m^3;
%1.32*10^-5=1.32*10^-5;%扩散系数
%C=Q/(4*pi.*1.32*10^-5.*t)^1.5*exp(-(x.^2+y.^2+z.^2)./4.*1.32*10^-5.*t);
%统一并申明一下单位
%C (mg/m^3)
%Q (mg/m^3)
%1.32*10^-5 (m^2/s)
%t (s)
%x,y,z (m)
%申明一下扩散的条件如下：
%1.模拟的球形空间的直径D=3.5m 在x方向安装15个，在y方向安装3个
%2.当经过时间10s时，可视为扩散均匀，其中时间可以通过一次实验或者相关文献查得。
randnum=unifrnd(1,10,10,1);%产生一组符合均匀分布的时间
for i=1:4
[x,y,z] = meshgrid(-1.75:.2:1.75,-1.75:.2:1.75,-1.75:.2:1.75);
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));
xslice=[-0.25,1.75];
yslice=[-0.25,1.25];
zslice=0.25;
subplot(2,2,i);
slice(x,y,z,C,xslice,yslice,zslice);
xlabel('x坐标');
ylabel('y坐标');
zlabel('z坐标');
title(['当时间为' num2str(randnum(i)) '时的浓度分布情况']);
colormap hsv
colorbar
end