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

核科学

高斯气体分布
& D. ]4 z& l' v4 g: r数值模拟
' J; ]3 c6 l& W$ k  L8 J作出动态分布图

923050019

function y=shuzhimoni
%syms x y z t;
8 ^( I& Z3 C  X; h: ]6 _%Q=10;%初始浓度定为10mg/m^3;
: c/ G; N: n7 }2 s* L( m% p  Y%1.32*10^-5=1.32*10^-5;%扩散系数
! C- g& O& A: \/ i' M%C=Q/(4*pi.*1.32*10^-5.*t)^1.5*exp(-(x.^2+y.^2+z.^2)./4.*1.32*10^-5.*t);
$ r! }* _9 k8 u; T) ~8 D2 g$ ~. j%统一并申明一下单位
%C (mg/m^3)
%Q (mg/m^3)
%1.32*10^-5 (m^2/s)
! r' K; ~  \6 ?  ^, q8 F" b2 u%t (s)
/ x1 |8 [. N" T+ \  l+ I; u%x,y,z (m)
3 l2 [# X% z, A* u%申明一下扩散的条件如下：
) ?: r8 b$ Z; O% g0 {%1.模拟的球形空间的直径D=3.5m 在x方向安装15个，在y方向安装3个
) x3 e! Q( n6 C6 n2 b- S, b" \%2.当经过时间10s时，可视为扩散均匀，其中时间可以通过一次实验或者相关文献查得。
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);
7 S4 \" `. m- o4 a9 h, W    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];
+ H2 h  _* c/ G/ [8 v& R" N, |" v- U. f    yslice=[-0.25,1.25];
( }& [4 `6 R2 C    zslice=0.25;
5 H0 X. k& J) o" D/ h    subplot(2,2,i);
    slice(x,y,z,C,xslice,yslice,zslice);
    xlabel('x坐标');
    ylabel('y坐标');
5 P6 Q* u# A) {! z# `    zlabel('z坐标');
  ]. r1 I* w8 H. h) Y    title(['当时间为' num2str(randnum(i)) '时的浓度分布情况']);
    colormap hsv
; R, j3 f$ e# ^1 Q    colorbar
end
3 x' @3 O1 E1 }

13ansys

看不懂什么意思  啊