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

核科学

高斯气体分布
数值模拟
5 l( A! J$ Z: h7 b1 y0 |作出动态分布图

923050019

function y=shuzhimoni
%syms x y z t;
( c, X& ]- ~6 G+ ?3 b%Q=10;%初始浓度定为10mg/m^3;
6 k+ _2 [1 {5 z& i) D  r! a& Y%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);
%统一并申明一下单位
' f# z" N, z1 Y7 ^2 E%C (mg/m^3)
8 X) N5 d  d, E: r. \%Q (mg/m^3)
%1.32*10^-5 (m^2/s)
%t (s)
%x,y,z (m)
%申明一下扩散的条件如下：
%1.模拟的球形空间的直径D=3.5m 在x方向安装15个，在y方向安装3个
; w5 d* F& X( i; s%2.当经过时间10s时，可视为扩散均匀，其中时间可以通过一次实验或者相关文献查得。
5 O5 i9 ?5 r: D  V5 Jrandnum=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);
$ R' O/ G$ G. F% i    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];
5 Q# S2 [  S) w    zslice=0.25;
! `. U$ }4 ]3 m* _# s7 r% e    subplot(2,2,i);
    slice(x,y,z,C,xslice,yslice,zslice);
    xlabel('x坐标');
: y. }4 w' B# y2 @    ylabel('y坐标');
    zlabel('z坐标');
    title(['当时间为' num2str(randnum(i)) '时的浓度分布情况']);
    colormap hsv
) R! L2 D& p5 c; i4 D$ Z    colorbar
end
4 X5 q( u) z# e3 q/ D1 ^

13ansys

看不懂什么意思  啊