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

核科学

高斯气体分布
- C- w# {) X* `/ Y数值模拟
作出动态分布图

923050019

function y=shuzhimoni
0 d: ^" a9 L1 D7 [7 Z%syms x y z t;
+ {! c$ f- a/ y%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)
2 |" S# w# R! n6 q2 ]: X1 ^1 R* }' 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个
; P4 Q3 c- n: s$ l& q9 I+ \%2.当经过时间10s时，可视为扩散均匀，其中时间可以通过一次实验或者相关文献查得。
0 {* v4 J4 ~: A! d! frandnum=unifrnd(1,10,10,1);%产生一组符合均匀分布的时间
' Y8 i* S; W* o; Efor i=1:4
  x7 f% {! t0 E% w2 z, o7 ^# I# x    [x,y,z] = meshgrid(-1.75:.2:1.75,-1.75:.2:1.75,-1.75:.2:1.75);
8 S: c2 e* T, ^9 T    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));
# s% u6 E' v. u    xslice=[-0.25,1.75];
    yslice=[-0.25,1.25];
    zslice=0.25;
    subplot(2,2,i);
" S6 x( ~$ @8 J8 g8 z) K    slice(x,y,z,C,xslice,yslice,zslice);
# ~& Y( R( S+ e2 I: C$ V& z    xlabel('x坐标');
4 U/ ~) W. m% s# v  O2 u! Z$ F    ylabel('y坐标');
: a) H4 w2 Z* c    zlabel('z坐标');
  W5 k3 ~" F6 Z. A* |$ P/ B/ _0 ]    title(['当时间为' num2str(randnum(i)) '时的浓度分布情况']);
    colormap hsv
& P- v6 s8 m6 c. K6 I: W2 ~) Z% M/ |    colorbar
. @, I5 a& s- e4 t9 mend

13ansys

看不懂什么意思  啊