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

核科学

高斯气体分布
数值模拟
) L% C3 W1 [& j3 M作出动态分布图

923050019

function y=shuzhimoni
* b3 l1 s7 b0 v' V9 x) d6 G%syms x y z t;
%Q=10;%初始浓度定为10mg/m^3;
%1.32*10^-5=1.32*10^-5;%扩散系数
+ z+ U8 ]6 \$ K+ d%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)
% x& I9 [; Q6 p6 A%1.32*10^-5 (m^2/s)
%t (s)
%x,y,z (m)
1 U+ c+ y# y5 r) V$ h$ M) R% ]) p%申明一下扩散的条件如下：
+ x% Y5 a% ]% Y0 r  _" G3 y%1.模拟的球形空间的直径D=3.5m 在x方向安装15个，在y方向安装3个
9 R2 R6 N& v2 n# }  Q  W& e%2.当经过时间10s时，可视为扩散均匀，其中时间可以通过一次实验或者相关文献查得。
' C- D5 t) B5 A& Orandnum=unifrnd(1,10,10,1);%产生一组符合均匀分布的时间
1 w( d( G4 D* @% `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];
) ?1 E% _1 H6 j+ v9 C1 d/ H3 Q    zslice=0.25;
    subplot(2,2,i);
    slice(x,y,z,C,xslice,yslice,zslice);
$ I! s/ S8 J' f. _  g. z    xlabel('x坐标');
    ylabel('y坐标');
    zlabel('z坐标');
) e9 P% V. h% q7 _    title(['当时间为' num2str(randnum(i)) '时的浓度分布情况']);
    colormap hsv
; R+ r+ r/ `' R; x' L    colorbar
end

13ansys

看不懂什么意思  啊