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

核科学

高斯气体分布
8 j+ F4 [: u' T9 L9 r, z数值模拟
) E- I3 _& W- p! p3 |6 L作出动态分布图

923050019

function y=shuzhimoni
%syms x y z t;
. z' o1 A* h5 K" @7 A+ w! p%Q=10;%初始浓度定为10mg/m^3;
* I6 n/ I& v4 e$ Z& E%1.32*10^-5=1.32*10^-5;%扩散系数
8 o7 K5 A# G1 }% m5 W%C=Q/(4*pi.*1.32*10^-5.*t)^1.5*exp(-(x.^2+y.^2+z.^2)./4.*1.32*10^-5.*t);
%统一并申明一下单位
0 B! }( }; r- G1 {: D- H6 E%C (mg/m^3)
! p9 k; |, q" U: ^%Q (mg/m^3)
  U& J) s% \# p% M- W8 T& B%1.32*10^-5 (m^2/s)
+ I% L0 @, \- _%t (s)
5 S1 R" j- C7 U3 k/ z1 ]2 V5 q%x,y,z (m)
%申明一下扩散的条件如下：
%1.模拟的球形空间的直径D=3.5m 在x方向安装15个，在y方向安装3个
1 q( d7 K+ U' L# `8 R  e* g%2.当经过时间10s时，可视为扩散均匀，其中时间可以通过一次实验或者相关文献查得。
randnum=unifrnd(1,10,10,1);%产生一组符合均匀分布的时间
3 j$ m: ~. y$ Y/ d8 T; U& d, Cfor 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));
0 M/ L, I2 c3 x( T# w- ]    xslice=[-0.25,1.75];
  I- L8 d7 h' H  q9 a    yslice=[-0.25,1.25];
    zslice=0.25;
    subplot(2,2,i);
/ @; e& T2 L3 K& u) M% E) A    slice(x,y,z,C,xslice,yslice,zslice);
% z/ X" A7 \4 L; B6 K/ m    xlabel('x坐标');
1 y7 Y! G2 L( O9 t    ylabel('y坐标');
    zlabel('z坐标');
2 y: z6 k1 G1 r* }' L- K    title(['当时间为' num2str(randnum(i)) '时的浓度分布情况']);
    colormap hsv
    colorbar
4 x# I% e: v4 Send
/ E8 J) C# k. A

13ansys

看不懂什么意思  啊