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

核科学

高斯气体分布
数值模拟
/ J+ N/ _! \/ ]/ i' x, X作出动态分布图

923050019

function y=shuzhimoni
%syms x y z t;
%Q=10;%初始浓度定为10mg/m^3;
  }  q- T& `. I3 j9 |  Q%1.32*10^-5=1.32*10^-5;%扩散系数
- Y7 z) R, J, ?, a* u%C=Q/(4*pi.*1.32*10^-5.*t)^1.5*exp(-(x.^2+y.^2+z.^2)./4.*1.32*10^-5.*t);
%统一并申明一下单位
5 R) d) {1 G3 X5 Z( A%C (mg/m^3)
* Y/ V+ a" c' b9 ^6 |/ K%Q (mg/m^3)
/ _+ z$ s+ h& ?9 y5 G; O) u5 P%1.32*10^-5 (m^2/s)
( u. g, X+ D; ?; \$ ~* ?6 a%t (s)
+ {6 s" ]# s( D. R% T$ O%x,y,z (m)
& j7 b9 ^; x$ o3 N8 E0 z' x%申明一下扩散的条件如下：
9 q2 t) R  z2 ^' x$ ~%1.模拟的球形空间的直径D=3.5m 在x方向安装15个，在y方向安装3个
%2.当经过时间10s时，可视为扩散均匀，其中时间可以通过一次实验或者相关文献查得。
$ p0 G  c- v7 S5 J+ v/ g" |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);
& [: i4 n6 ]% N4 j2 A    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));
+ C" @* v9 e0 u0 P0 e    xslice=[-0.25,1.75];
    yslice=[-0.25,1.25];
& |, b: k" W9 N+ m5 q3 Y4 y- P/ y$ T    zslice=0.25;
' n7 }2 X: M5 z5 _5 a4 e4 {    subplot(2,2,i);
    slice(x,y,z,C,xslice,yslice,zslice);
    xlabel('x坐标');
    ylabel('y坐标');
    zlabel('z坐标');
  W4 P1 ~0 Y7 _+ p; T    title(['当时间为' num2str(randnum(i)) '时的浓度分布情况']);
7 g$ Q( q6 Y+ H( ^7 K" x: t; N0 o    colormap hsv
" @9 |+ v0 O5 h7 q3 q# E    colorbar
$ k, _# X1 t9 x" eend

13ansys

看不懂什么意思  啊