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

核科学

高斯气体分布
数值模拟
作出动态分布图

13ansys

看不懂什么意思  啊

923050019

function y=shuzhimoni
%syms x y z t;
%Q=10;%初始浓度定为10mg/m^3;
- }, [& @& U8 V! Z+ B. t3 Q' b%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);
%统一并申明一下单位
" M6 e; ^" D7 W6 j1 Z( m; G. J%C (mg/m^3)
%Q (mg/m^3)
8 W: H$ ?. k) ]4 X5 @%1.32*10^-5 (m^2/s)
%t (s)
6 Y7 A/ c3 B1 |$ z2 U; \$ I%x,y,z (m)
! A) p  \6 u0 p0 i$ r. r2 a$ d%申明一下扩散的条件如下：
%1.模拟的球形空间的直径D=3.5m 在x方向安装15个，在y方向安装3个
%2.当经过时间10s时，可视为扩散均匀，其中时间可以通过一次实验或者相关文献查得。
randnum=unifrnd(1,10,10,1);%产生一组符合均匀分布的时间
for i=1:4
0 `& P+ N8 w0 s, _' z6 Y    [x,y,z] = meshgrid(-1.75:.2:1.75,-1.75:.2:1.75,-1.75:.2:1.75);
8 _' @' i2 d" C( Y. _    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];
; b1 W) @. y% I4 z, k  Y7 l    yslice=[-0.25,1.25];
4 y- |+ o& x% p. i7 e  M    zslice=0.25;
0 m3 {, Z, y) Y  G    subplot(2,2,i);
    slice(x,y,z,C,xslice,yslice,zslice);
: G& X  _" R. B! E; R) \    xlabel('x坐标');
    ylabel('y坐标');
    zlabel('z坐标');
    title(['当时间为' num2str(randnum(i)) '时的浓度分布情况']);
    colormap hsv
4 v3 H3 C1 O' \% Z% e& ]    colorbar
3 b* R6 I. ^# K# M7 gend
! u2 F) i6 P$ G# ]0 q) S