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

核科学

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

923050019

function y=shuzhimoni
%syms x y z t;
%Q=10;%初始浓度定为10mg/m^3;
%1.32*10^-5=1.32*10^-5;%扩散系数
8 f5 `+ I8 ?, F  t& a# e5 J- `* m%C=Q/(4*pi.*1.32*10^-5.*t)^1.5*exp(-(x.^2+y.^2+z.^2)./4.*1.32*10^-5.*t);
# U! E  U# g8 x# _  e0 m%统一并申明一下单位
- C* V# i/ \- J. k/ P# h0 d%C (mg/m^3)
%Q (mg/m^3)
%1.32*10^-5 (m^2/s)
%t (s)
$ n4 P+ [6 t9 T%x,y,z (m)
%申明一下扩散的条件如下：
%1.模拟的球形空间的直径D=3.5m 在x方向安装15个，在y方向安装3个
%2.当经过时间10s时，可视为扩散均匀，其中时间可以通过一次实验或者相关文献查得。
2 u* c" X9 s5 E! e. irandnum=unifrnd(1,10,10,1);%产生一组符合均匀分布的时间
$ n& i7 x; R, F% d9 q; U4 @: ~for i=1:4
+ r+ O/ Y: z* m- J, r    [x,y,z] = meshgrid(-1.75:.2:1.75,-1.75:.2:1.75,-1.75:.2:1.75);
4 J' }9 Q: J4 p' t( V, S) N    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];
+ c: ]8 P6 C" f7 t' e    zslice=0.25;
    subplot(2,2,i);
5 H1 `. f) S. I3 ~! q6 t$ V* v! _! M    slice(x,y,z,C,xslice,yslice,zslice);
    xlabel('x坐标');
7 F2 l( }$ n+ l( w    ylabel('y坐标');
    zlabel('z坐标');
% ]- k. C4 m" D6 ~! R( \% `    title(['当时间为' num2str(randnum(i)) '时的浓度分布情况']);
    colormap hsv
# K5 ^9 u  f/ ~6 O: D$ A* b    colorbar
end
. a8 O$ t: k' d4 Y& W2 r

13ansys

看不懂什么意思  啊