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

核科学

高斯气体分布
* m3 U  `* [# A数值模拟
作出动态分布图

923050019

function y=shuzhimoni
- R, [- e. I4 g" @%syms x y z t;
6 K# i1 E; J$ C- k% d" ]%Q=10;%初始浓度定为10mg/m^3;
% _- h! E: n  n8 Q. Z- D%1.32*10^-5=1.32*10^-5;%扩散系数
& C% G  s$ P0 j+ p4 L/ F8 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);
: z. V7 K* Q) ?& J" ]%统一并申明一下单位
%C (mg/m^3)
0 @$ y. ~; Q  A7 c3 \( p%Q (mg/m^3)
%1.32*10^-5 (m^2/s)
7 L1 _/ f1 d- M0 a1 g# G%t (s)
" K/ z. X3 i# ?, R: t%x,y,z (m)
%申明一下扩散的条件如下：
%1.模拟的球形空间的直径D=3.5m 在x方向安装15个，在y方向安装3个
$ v* _- r+ ]) Q  ]. r%2.当经过时间10s时，可视为扩散均匀，其中时间可以通过一次实验或者相关文献查得。
* L3 R7 |; ?- _& P; trandnum=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));
5 |# {7 L9 [$ s- W    xslice=[-0.25,1.75];
    yslice=[-0.25,1.25];
    zslice=0.25;
5 j1 [* p1 g( d- D/ f% H    subplot(2,2,i);
    slice(x,y,z,C,xslice,yslice,zslice);
5 _) s7 G0 f1 F0 s* {2 s9 ?    xlabel('x坐标');
    ylabel('y坐标');
    zlabel('z坐标');
    title(['当时间为' num2str(randnum(i)) '时的浓度分布情况']);
    colormap hsv
' Y9 L2 o4 s" O4 P    colorbar
4 Q% H* c! q: f' S7 h, Dend
+ r+ ]/ t& d2 J; K% J& n# j

13ansys

看不懂什么意思  啊