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

核科学

高斯气体分布
9 @) @1 y7 M5 [+ I, n% `数值模拟
作出动态分布图

13ansys

看不懂什么意思  啊

923050019

function y=shuzhimoni
%syms x y z t;
%Q=10;%初始浓度定为10mg/m^3;
%1.32*10^-5=1.32*10^-5;%扩散系数
+ {/ S( U5 ], o- \! k& t3 A%C=Q/(4*pi.*1.32*10^-5.*t)^1.5*exp(-(x.^2+y.^2+z.^2)./4.*1.32*10^-5.*t);
: T3 x7 h2 ]# A' B%统一并申明一下单位
%C (mg/m^3)
: y1 z1 a' y$ }%Q (mg/m^3)
%1.32*10^-5 (m^2/s)
%t (s)
1 g5 r3 K8 h/ z%x,y,z (m)
%申明一下扩散的条件如下：
: V" p( C( L" ]7 b& i( m7 j%1.模拟的球形空间的直径D=3.5m 在x方向安装15个，在y方向安装3个
%2.当经过时间10s时，可视为扩散均匀，其中时间可以通过一次实验或者相关文献查得。
3 P+ e" V4 x7 D1 [5 Erandnum=unifrnd(1,10,10,1);%产生一组符合均匀分布的时间
- f+ r8 r* M, s! w/ ffor i=1:4
4 `2 s6 K9 f: T8 _! b8 p) d* k; w0 {    [x,y,z] = meshgrid(-1.75:.2:1.75,-1.75:.2:1.75,-1.75:.2:1.75);
8 S/ J; Z; z; n; `4 G5 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));
4 O1 f) [2 n* q" ?, ]4 C" _( Y    xslice=[-0.25,1.75];
& n! H; c! Q* |    yslice=[-0.25,1.25];
5 g7 s, [9 {( K( L1 |5 P/ J; K    zslice=0.25;
    subplot(2,2,i);
) R: R1 Q$ |7 }5 J    slice(x,y,z,C,xslice,yslice,zslice);
    xlabel('x坐标');
$ g0 q/ O" s  S  E  H    ylabel('y坐标');
    zlabel('z坐标');
    title(['当时间为' num2str(randnum(i)) '时的浓度分布情况']);
    colormap hsv
' v+ h" r5 \6 }6 d7 ]7 z: l$ l    colorbar
$ b8 d' [1 `2 k1 M- j# uend
8 W/ v; w2 a% j