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

核科学

高斯气体分布
3 E3 Z, S! _) w: K" q2 d数值模拟
作出动态分布图

923050019

function y=shuzhimoni
%syms x y z t;
6 m8 D' _3 l: b) P. }%Q=10;%初始浓度定为10mg/m^3;
" f& h- x$ j1 S' ~' H( _9 D%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);
* ?' X! S7 l" X, C8 P/ g%统一并申明一下单位
%C (mg/m^3)
%Q (mg/m^3)
& ]9 V7 u' s: V, `6 P) h; o- c%1.32*10^-5 (m^2/s)
0 i  Q/ a5 ?, f- J' c%t (s)
%x,y,z (m)
4 U2 J! Y' ~$ L2 @1 }9 }1 t%申明一下扩散的条件如下：
0 V; T+ w: Q* W- ^%1.模拟的球形空间的直径D=3.5m 在x方向安装15个，在y方向安装3个
%2.当经过时间10s时，可视为扩散均匀，其中时间可以通过一次实验或者相关文献查得。
( ]4 w. V6 B# _( B6 rrandnum=unifrnd(1,10,10,1);%产生一组符合均匀分布的时间
for i=1:4
6 U+ H" B" G- E# M* V    [x,y,z] = meshgrid(-1.75:.2:1.75,-1.75:.2:1.75,-1.75:.2:1.75);
' z+ b1 v% _6 ^    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];
    yslice=[-0.25,1.25];
" C! H) ~, @) q6 P    zslice=0.25;
    subplot(2,2,i);
    slice(x,y,z,C,xslice,yslice,zslice);
6 Z0 ~, F& i+ E, z; q) M9 h9 ]7 Y    xlabel('x坐标');
    ylabel('y坐标');
( {% z+ Y, z& R# w    zlabel('z坐标');
    title(['当时间为' num2str(randnum(i)) '时的浓度分布情况']);
    colormap hsv
    colorbar
8 U% J) K+ Y- B4 h! Cend
7 j9 ^% q0 o  p" o# l( g$ V2 u- ~

13ansys

看不懂什么意思  啊