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

核科学

高斯气体分布
数值模拟
8 D1 m4 O$ u$ [: ~) t* w9 L作出动态分布图

923050019

function y=shuzhimoni
! _2 `# g7 X: |%syms x y z t;
%Q=10;%初始浓度定为10mg/m^3;
( A3 s. l5 P- X5 w) e%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);
%统一并申明一下单位
%C (mg/m^3)
%Q (mg/m^3)
6 D' n1 W/ Z$ Y6 a1 z! ~$ m4 {%1.32*10^-5 (m^2/s)
%t (s)
%x,y,z (m)
%申明一下扩散的条件如下：
$ _2 u7 ~" \, v5 h2 P* c%1.模拟的球形空间的直径D=3.5m 在x方向安装15个，在y方向安装3个
%2.当经过时间10s时，可视为扩散均匀，其中时间可以通过一次实验或者相关文献查得。
randnum=unifrnd(1,10,10,1);%产生一组符合均匀分布的时间
for i=1:4
    [x,y,z] = meshgrid(-1.75:.2:1.75,-1.75:.2:1.75,-1.75:.2:1.75);
1 d3 i" k" r1 A! O: h( D, `5 F    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));
# V) Z7 d* v& I0 |3 n3 F# z    xslice=[-0.25,1.75];
! s: C' m; Y/ o$ N    yslice=[-0.25,1.25];
$ A( a3 I5 X, g9 M  m3 Y! T" O    zslice=0.25;
/ W* u3 N" Q" ]    subplot(2,2,i);
' o0 E0 P' n& l' M2 e* S# V0 I    slice(x,y,z,C,xslice,yslice,zslice);
8 W; [2 B! h1 y8 u# P* Z    xlabel('x坐标');
; V# }2 b  F* U6 N9 r. D  @    ylabel('y坐标');
! ~$ S7 Q; y$ n0 V+ D7 G    zlabel('z坐标');
    title(['当时间为' num2str(randnum(i)) '时的浓度分布情况']);
8 M; ?5 n  j5 A  ?6 B$ ]    colormap hsv
    colorbar
end

13ansys

看不懂什么意思  啊