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

核科学

高斯气体分布
数值模拟
5 w0 N# ?8 B6 a$ ?; D$ B: }2 t( V作出动态分布图

923050019

function y=shuzhimoni
, m( V0 a7 D5 J' C! ^" w%syms x y z t;
/ n7 E5 B$ a$ X# o%Q=10;%初始浓度定为10mg/m^3;
- x6 @" M+ v3 Z) H1 O! S%1.32*10^-5=1.32*10^-5;%扩散系数
1 l  i% C0 T5 |0 q%C=Q/(4*pi.*1.32*10^-5.*t)^1.5*exp(-(x.^2+y.^2+z.^2)./4.*1.32*10^-5.*t);
%统一并申明一下单位
- j" P, E( d  ]%C (mg/m^3)
%Q (mg/m^3)
%1.32*10^-5 (m^2/s)
' p; [6 }" }) L) h* _, U%t (s)
%x,y,z (m)
0 o$ L4 ^8 ^, K* P# ~4 ~0 `%申明一下扩散的条件如下：
2 S8 h* L) j" x* H%1.模拟的球形空间的直径D=3.5m 在x方向安装15个，在y方向安装3个
%2.当经过时间10s时，可视为扩散均匀，其中时间可以通过一次实验或者相关文献查得。
randnum=unifrnd(1,10,10,1);%产生一组符合均匀分布的时间
4 P( s0 f3 B$ N0 xfor i=1:4
7 q& w' T/ E. x9 D6 _3 \% B2 q4 y* |; V    [x,y,z] = meshgrid(-1.75:.2:1.75,-1.75:.2:1.75,-1.75:.2:1.75);
2 T4 R7 T* O, G) l9 Y/ C    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];
6 h; E$ @/ c: n  f- X3 M    yslice=[-0.25,1.25];
  p2 d( y+ \) F% m2 Q3 m    zslice=0.25;
5 W& a2 p+ q$ ?. y( X6 N. u    subplot(2,2,i);
    slice(x,y,z,C,xslice,yslice,zslice);
2 x; e7 y1 y& {- W* J    xlabel('x坐标');
2 R( n' B7 t/ i4 X4 X+ O* C$ a) Y    ylabel('y坐标');
& k) g+ O2 u$ b1 n* d    zlabel('z坐标');
    title(['当时间为' num2str(randnum(i)) '时的浓度分布情况']);
    colormap hsv
    colorbar
end
! @: m  a( f0 A/ E+ S% o: s

13ansys

看不懂什么意思  啊