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

核科学

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

923050019

function y=shuzhimoni
. A4 q5 h& l& Z8 i) E$ w- P9 t: v%syms x y z t;
6 U/ E$ C, ^3 V( [* t%Q=10;%初始浓度定为10mg/m^3;
%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);
%统一并申明一下单位
7 R# C% [, T  a' f/ m%C (mg/m^3)
%Q (mg/m^3)
%1.32*10^-5 (m^2/s)
%t (s)
) j/ y% m. v$ t. ?) i4 r( |%x,y,z (m)
%申明一下扩散的条件如下：
0 W7 s3 V2 W3 |( [  ]. _1 q%1.模拟的球形空间的直径D=3.5m 在x方向安装15个，在y方向安装3个
5 b2 [" B) n5 T5 s$ W% B; u%2.当经过时间10s时，可视为扩散均匀，其中时间可以通过一次实验或者相关文献查得。
randnum=unifrnd(1,10,10,1);%产生一组符合均匀分布的时间
/ E6 k6 B* v1 p" v: s' w9 e, A4 Dfor i=1:4
! [# x" j* o2 H+ V* r& J- r    [x,y,z] = meshgrid(-1.75:.2:1.75,-1.75:.2:1.75,-1.75:.2:1.75);
$ C3 `2 z, H/ i& H    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));
8 z( w! u9 M6 z7 p1 }5 A0 U6 r" y    xslice=[-0.25,1.75];
    yslice=[-0.25,1.25];
% k: j: A0 z( Y    zslice=0.25;
    subplot(2,2,i);
    slice(x,y,z,C,xslice,yslice,zslice);
    xlabel('x坐标');
. }- c0 P+ Q) u/ m    ylabel('y坐标');
    zlabel('z坐标');
( w, X' p" ~+ N1 q3 p( b6 N8 N    title(['当时间为' num2str(randnum(i)) '时的浓度分布情况']);
    colormap hsv
    colorbar
end

13ansys

看不懂什么意思  啊