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

核科学

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

923050019

function y=shuzhimoni
5 u0 \! e. v/ P8 x( |$ H& z% `%syms x y z t;
%Q=10;%初始浓度定为10mg/m^3;
. i1 q! B7 ]: V%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);
2 h- |) m7 ~* t9 Y( @%统一并申明一下单位
%C (mg/m^3)
%Q (mg/m^3)
%1.32*10^-5 (m^2/s)
%t (s)
%x,y,z (m)
%申明一下扩散的条件如下：
" l3 n& U1 X6 N! e2 e6 Q; |% a%1.模拟的球形空间的直径D=3.5m 在x方向安装15个，在y方向安装3个
2 u: A* r4 M6 A: {! h  M0 b%2.当经过时间10s时，可视为扩散均匀，其中时间可以通过一次实验或者相关文献查得。
( Q3 p8 f3 T  \7 qrandnum=unifrnd(1,10,10,1);%产生一组符合均匀分布的时间
4 e. Z9 u: i$ X/ H% p  \for i=1:4
    [x,y,z] = meshgrid(-1.75:.2:1.75,-1.75:.2:1.75,-1.75:.2:1.75);
2 G5 u" J4 z5 v+ N/ {$ N9 L    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));
5 Y. e8 j' Q$ U' \8 g& G    xslice=[-0.25,1.75];
    yslice=[-0.25,1.25];
" [8 @$ W8 q3 S& M    zslice=0.25;
- I& r& |8 y: v+ K; Y2 |1 v    subplot(2,2,i);
7 }4 r% d* p) j& R9 X9 E    slice(x,y,z,C,xslice,yslice,zslice);
! I/ F! q) m$ o    xlabel('x坐标');
    ylabel('y坐标');
9 w6 G  T7 Q# M9 H6 Q' k7 w  U# W    zlabel('z坐标');
7 h8 u2 l/ ?/ e6 f    title(['当时间为' num2str(randnum(i)) '时的浓度分布情况']);
    colormap hsv
    colorbar
  @, T. A0 ~: e+ l( ?end

13ansys

看不懂什么意思  啊