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

核科学

高斯气体分布
7 E# {+ \& v  _- v+ T. [数值模拟
* c% l( g+ M1 P4 D  x7 e9 ]5 R9 v作出动态分布图

923050019

function y=shuzhimoni
%syms x y z t;
" Z2 S$ N. o! k; O%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);
' Y5 |4 S6 P% t. S# G4 m%统一并申明一下单位
%C (mg/m^3)
% j0 J- P9 j3 E2 \%Q (mg/m^3)
%1.32*10^-5 (m^2/s)
%t (s)
%x,y,z (m)
7 `- ^9 y6 ?& E1 a; p. P%申明一下扩散的条件如下：
, }, R( P5 Y. s) q%1.模拟的球形空间的直径D=3.5m 在x方向安装15个，在y方向安装3个
%2.当经过时间10s时，可视为扩散均匀，其中时间可以通过一次实验或者相关文献查得。
( f( V* c6 \" t0 H& t2 nrandnum=unifrnd(1,10,10,1);%产生一组符合均匀分布的时间
for i=1:4
4 C8 D; \8 X7 _; [& N4 V/ r    [x,y,z] = meshgrid(-1.75:.2:1.75,-1.75:.2:1.75,-1.75:.2:1.75);
    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];
6 m( \1 O" k; v5 W* q    zslice=0.25;
9 |1 j2 M* h: a) x$ X% B" E: c    subplot(2,2,i);
    slice(x,y,z,C,xslice,yslice,zslice);
    xlabel('x坐标');
4 a9 d  U5 Z) x    ylabel('y坐标');
6 h% n& F, u4 r/ t4 }5 ~+ T+ C$ M1 s    zlabel('z坐标');
    title(['当时间为' num2str(randnum(i)) '时的浓度分布情况']);
5 Z" ^% V4 J3 t7 F: _: u    colormap hsv
/ v* U' C- @4 f2 ^1 V# Q    colorbar
end

13ansys

看不懂什么意思  啊