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

核科学

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

923050019

function y=shuzhimoni
7 u) l2 A2 l$ ]3 f8 ?. n+ `%syms x y z t;
  ?* B, `" X# z1 R1 P%Q=10;%初始浓度定为10mg/m^3;
  w( t/ Y4 q9 q5 B%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);
%统一并申明一下单位
8 R7 x' ?7 n: ~1 W* v3 Z) J%C (mg/m^3)
%Q (mg/m^3)
- o. c+ W( Z' b! ~$ |7 \. q%1.32*10^-5 (m^2/s)
) O' Y: H8 r+ X3 ^( Q; @5 q%t (s)
2 e! A/ }/ O% o" r( I3 g) c% P%x,y,z (m)
3 t5 G! i! b- Q! I/ ^# T%申明一下扩散的条件如下：
" C7 e6 Z6 {* {%1.模拟的球形空间的直径D=3.5m 在x方向安装15个，在y方向安装3个
# D+ S: y2 R. E& A' R# E- R%2.当经过时间10s时，可视为扩散均匀，其中时间可以通过一次实验或者相关文献查得。
1 G3 c. ~4 N3 {0 _+ F3 m' Orandnum=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);
. r! |: E) o: g- E8 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));
& X1 t- I/ Y* j8 p, ^1 t, k- \6 c    xslice=[-0.25,1.75];
    yslice=[-0.25,1.25];
    zslice=0.25;
    subplot(2,2,i);
+ b- Z! q9 O  R+ m8 r; r    slice(x,y,z,C,xslice,yslice,zslice);
, \( u+ D+ A# B: x5 v    xlabel('x坐标');
6 O" Q* H1 f! v* z6 l9 A, C0 E    ylabel('y坐标');
! L. F+ u+ S  ?    zlabel('z坐标');
3 L( l  `" r8 ~7 s" S    title(['当时间为' num2str(randnum(i)) '时的浓度分布情况']);
    colormap hsv
    colorbar
4 \% Z8 u5 M% M& s: U) |end
/ a1 X9 f9 D% A* G* O* p" y+ w* K+ }- ]

13ansys

看不懂什么意思  啊