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

核科学

高斯气体分布
数值模拟
( f3 Y! `$ u1 {- c( C作出动态分布图

923050019

function y=shuzhimoni
. x) j+ ~0 z+ U' M8 F%syms x y z t;
- ^% j/ o1 r' m1 f%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);
%统一并申明一下单位
%C (mg/m^3)
%Q (mg/m^3)
/ B3 a) U/ i1 l7 z: N+ L/ d%1.32*10^-5 (m^2/s)
%t (s)
, t% u* A& D. U+ ?3 B%x,y,z (m)
%申明一下扩散的条件如下：
: g. E6 p7 Y+ v- K- M! r5 x%1.模拟的球形空间的直径D=3.5m 在x方向安装15个，在y方向安装3个
4 a' a3 ]) p2 x, a! x: Q8 R" E%2.当经过时间10s时，可视为扩散均匀，其中时间可以通过一次实验或者相关文献查得。
6 x* F, f4 B) Rrandnum=unifrnd(1,10,10,1);%产生一组符合均匀分布的时间
7 [' t/ p0 }) w- Q: Gfor i=1:4
, X9 B( |0 ^7 O/ N7 B& d, R1 J    [x,y,z] = meshgrid(-1.75:.2:1.75,-1.75:.2:1.75,-1.75:.2:1.75);
& p' \9 F  I+ [% ?% y" q' k5 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));
7 F" s8 y$ S% {' z- A' y/ f    xslice=[-0.25,1.75];
  w! H' G3 d, B1 V. ?    yslice=[-0.25,1.25];
& q  C) p: ]0 A$ t" ~: u& s    zslice=0.25;
    subplot(2,2,i);
    slice(x,y,z,C,xslice,yslice,zslice);
* N- _9 v7 O, ?, y( U/ J    xlabel('x坐标');
+ B* b9 l& [# N+ ]+ W/ j& P) }  G5 o: ]    ylabel('y坐标');
    zlabel('z坐标');
    title(['当时间为' num2str(randnum(i)) '时的浓度分布情况']);
    colormap hsv
. d2 z0 @" B: K+ R    colorbar
' f& ?' Z3 z7 k9 \6 x6 Q; p% M6 Qend

13ansys

看不懂什么意思  啊