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

核科学

高斯气体分布
- @; V( H8 u+ S! ^1 O' `7 \# @数值模拟
+ B2 O% @7 U: {  c作出动态分布图

923050019

function y=shuzhimoni
6 g% h8 _: p  H9 K%syms x y z t;
) K1 Y( b1 i4 i! s%Q=10;%初始浓度定为10mg/m^3;
%1.32*10^-5=1.32*10^-5;%扩散系数
, \  F# j5 i* b7 }%C=Q/(4*pi.*1.32*10^-5.*t)^1.5*exp(-(x.^2+y.^2+z.^2)./4.*1.32*10^-5.*t);
2 w+ C2 A6 U! E2 a%统一并申明一下单位
' A, T/ j  r+ b$ }( g& \5 C%C (mg/m^3)
%Q (mg/m^3)
7 N, o9 S" ~, G& `  K%1.32*10^-5 (m^2/s)
%t (s)
& e% {$ N3 z* g7 U1 a5 D" {%x,y,z (m)
%申明一下扩散的条件如下：
# R) w  f4 w" F6 g%1.模拟的球形空间的直径D=3.5m 在x方向安装15个，在y方向安装3个
, I7 o- ^# p" R9 ?# F/ p& F: k2 G% m; X%2.当经过时间10s时，可视为扩散均匀，其中时间可以通过一次实验或者相关文献查得。
! c/ O1 m, b  k* w& A6 V1 rrandnum=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));
    xslice=[-0.25,1.75];
    yslice=[-0.25,1.25];
    zslice=0.25;
    subplot(2,2,i);
$ ?- R; _; ~0 `( R% }% v7 c$ N    slice(x,y,z,C,xslice,yslice,zslice);
4 i# J" j3 q# E* ?0 [1 k1 `    xlabel('x坐标');
5 R5 y( k7 A# Y+ {5 q6 E/ c    ylabel('y坐标');
1 g! Y7 T7 ], J5 N1 t0 k; `# s; v    zlabel('z坐标');
    title(['当时间为' num2str(randnum(i)) '时的浓度分布情况']);
7 O+ B0 O5 S9 p# u" @# w/ n& T    colormap hsv
    colorbar
end

13ansys

看不懂什么意思  啊