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

核科学

高斯气体分布
" k# A$ G! p+ A, w0 W7 t8 [1 @/ s数值模拟
' ?+ a9 V! P+ ], \- ]  ?作出动态分布图

923050019

function y=shuzhimoni
! l% ~$ v- y# I  B6 f: N%syms x y z t;
%Q=10;%初始浓度定为10mg/m^3;
% F: K$ d  l' E7 i5 f+ _- Q& s%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);
$ M2 t' a6 s3 i# ?9 _%统一并申明一下单位
%C (mg/m^3)
2 b2 k% E9 [. S%Q (mg/m^3)
%1.32*10^-5 (m^2/s)
%t (s)
* q: r" l/ v6 ?% C1 V%x,y,z (m)
4 R1 ?0 ^: q5 E. j%申明一下扩散的条件如下：
%1.模拟的球形空间的直径D=3.5m 在x方向安装15个，在y方向安装3个
%2.当经过时间10s时，可视为扩散均匀，其中时间可以通过一次实验或者相关文献查得。
, Y" `& j  d0 ?1 V# Drandnum=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);
8 K" C# u* L% ]/ J' E8 p7 x/ I9 O% [    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 P' ~6 |: w2 V: h# l; D    xslice=[-0.25,1.75];
    yslice=[-0.25,1.25];
8 F7 F" j; f: `- N    zslice=0.25;
    subplot(2,2,i);
: s5 }& F  U/ l* H6 W5 m; n    slice(x,y,z,C,xslice,yslice,zslice);
    xlabel('x坐标');
' F* d% w3 Q; b. z; \    ylabel('y坐标');
' ?9 b1 u0 V& L( S* G    zlabel('z坐标');
' ~8 A1 ^" C3 ?    title(['当时间为' num2str(randnum(i)) '时的浓度分布情况']);
    colormap hsv
    colorbar
end
- w/ a) W* ?" s) F- ]6 K" N8 Q

13ansys

看不懂什么意思  啊