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

核科学

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

923050019

function y=shuzhimoni
% @- w7 \/ n/ p. n" C8 k& r%syms x y z t;
%Q=10;%初始浓度定为10mg/m^3;
- x  N6 R0 [+ J3 \3 d) f1 z/ Y%1.32*10^-5=1.32*10^-5;%扩散系数
1 T$ X' z$ p( W' W+ F6 n0 |+ p%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 f) N) o" f0 h2 A9 o7 @2 [0 x%1.32*10^-5 (m^2/s)
%t (s)
%x,y,z (m)
%申明一下扩散的条件如下：
%1.模拟的球形空间的直径D=3.5m 在x方向安装15个，在y方向安装3个
8 Q: N3 a* _3 d) {4 g%2.当经过时间10s时，可视为扩散均匀，其中时间可以通过一次实验或者相关文献查得。
randnum=unifrnd(1,10,10,1);%产生一组符合均匀分布的时间
5 S) o; G5 p) m( g9 G0 E; m* a( [for i=1:4
    [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));
* Q3 g5 g" h2 T2 t2 y& n. {3 c    xslice=[-0.25,1.75];
& Y! E/ m7 Z5 R3 l/ M    yslice=[-0.25,1.25];
7 W8 d, L  ?5 ?% Z5 j& z% f    zslice=0.25;
    subplot(2,2,i);
$ S. K5 b& {& t7 y3 M! j2 K4 ?    slice(x,y,z,C,xslice,yslice,zslice);
    xlabel('x坐标');
7 [+ N  m/ |# F7 _7 b    ylabel('y坐标');
" F1 Q$ _+ h) q7 ~" C2 e: D# U$ b    zlabel('z坐标');
8 o  @2 R4 k6 q4 _1 j* C' g2 w2 u  }    title(['当时间为' num2str(randnum(i)) '时的浓度分布情况']);
- z8 p7 B3 w7 x  z9 m    colormap hsv
  S/ Z' Y' Y+ G    colorbar
+ c- `) K% {) W& ~, Mend

13ansys

看不懂什么意思  啊