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

核科学

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

923050019

function y=shuzhimoni
0 U4 T$ T: x1 z: y. f  H+ B- O%syms x y z t;
1 i, C: l/ ?8 d2 k/ F6 _( j& X%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);
/ ^) a6 e# _% c0 P2 s7 w' o9 P%统一并申明一下单位
%C (mg/m^3)
%Q (mg/m^3)
%1.32*10^-5 (m^2/s)
/ K* O0 j$ ^* k5 i( v%t (s)
1 A% i, R' O1 ?! W7 s5 o%x,y,z (m)
%申明一下扩散的条件如下：
9 d  c9 g9 T; k+ b7 d$ H+ @& L. b* h%1.模拟的球形空间的直径D=3.5m 在x方向安装15个，在y方向安装3个
+ Z, D- o- t* a%2.当经过时间10s时，可视为扩散均匀，其中时间可以通过一次实验或者相关文献查得。
randnum=unifrnd(1,10,10,1);%产生一组符合均匀分布的时间
9 h) C" F$ f) u  T8 \8 p# b! ?for i=1:4
    [x,y,z] = meshgrid(-1.75:.2:1.75,-1.75:.2:1.75,-1.75:.2:1.75);
1 }' N" n$ H1 l3 x3 }& n! w) K    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));
, k4 \7 q0 N4 I- }    xslice=[-0.25,1.75];
    yslice=[-0.25,1.25];
    zslice=0.25;
; [$ `" u  C% B' R- g' f: |    subplot(2,2,i);
    slice(x,y,z,C,xslice,yslice,zslice);
    xlabel('x坐标');
    ylabel('y坐标');
    zlabel('z坐标');
% P' d( \& w& b8 m    title(['当时间为' num2str(randnum(i)) '时的浓度分布情况']);
. ]7 R4 G/ g: u  t) u( t" M7 U' z& i5 x    colormap hsv
! x- S; L  o8 K, g    colorbar
end

13ansys

看不懂什么意思  啊