QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 4650|回复: 2
打印 上一主题 下一主题

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

[复制链接]
字体大小: 正常 放大
核科学 实名认证       

1

主题

3

听众

19

积分

升级  14.74%

该用户从未签到

自我介绍
我我我

窝窝

Ooo
跳转到指定楼层
1#
发表于 2010-5-5 18:11 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
高斯气体分布
0 x/ B, [( N9 U, b: i数值模拟
  ~( r- J" ~( t9 J5 T) V作出动态分布图
zan
转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
923050019 实名认证       

22

主题

8

听众

769

积分

升级  42.25%

  • TA的每日心情
    开心
    2013-5-6 20:58
  • 签到天数: 3 天

    [LV.2]偶尔看看I

    新人进步奖 最具活力勋章

    群组数学建模

    群组LINGO

    群组代数与数论

    群组C 语言讨论组

    function y=shuzhimoni
    3 k5 v2 L6 P! W8 }! r/ d%syms x y z t;" X1 d1 Q+ l9 c) P' r
    %Q=10;%初始浓度定为10mg/m^3;
    1 C, T& L4 E( X6 w%1.32*10^-5=1.32*10^-5;%扩散系数  I4 B: P8 s  r% |
    %C=Q/(4*pi.*1.32*10^-5.*t)^1.5*exp(-(x.^2+y.^2+z.^2)./4.*1.32*10^-5.*t);' J% K$ k+ C( p
    %统一并申明一下单位
    ( B5 O$ [; k9 B8 Q6 A% C; f%C (mg/m^3)6 N* c" J) v, B( a; O
    %Q (mg/m^3)
    6 |0 o! M; L9 ^/ U  W%1.32*10^-5 (m^2/s). Z1 N  w. i: H3 n/ U. G
    %t (s)$ t: r/ o( d$ ?
    %x,y,z (m)
    * D6 {$ y( E; |/ ]; r& s8 _( M" H; N%申明一下扩散的条件如下:" R4 v1 c  _, S* I( A% z* @9 }  B2 {
    %1.模拟的球形空间的直径D=3.5m 在x方向安装15个,在y方向安装3个% }' `8 M1 ^0 j. Y
    %2.当经过时间10s时,可视为扩散均匀,其中时间可以通过一次实验或者相关文献查得。
    * H; Y3 j8 R1 }+ ?0 y8 G& Jrandnum=unifrnd(1,10,10,1);%产生一组符合均匀分布的时间- c# R5 D6 Z" {
    for i=1:47 U- Y" b: P. ?. P; D
        [x,y,z] = meshgrid(-1.75:.2:1.75,-1.75:.2:1.75,-1.75:.2:1.75);
    $ k) U* |4 R8 m) P8 \6 Z7 h! K    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 g1 \: d4 P3 D+ o$ A0 N7 P    xslice=[-0.25,1.75];! w1 s9 y( A3 l: n2 r) ?
        yslice=[-0.25,1.25];
    * E+ K: ^) p" w    zslice=0.25;
    4 M( }" s% ~1 I0 u) O4 j2 A    subplot(2,2,i);
    ( l5 J8 u' n' T( t- @3 |    slice(x,y,z,C,xslice,yslice,zslice);/ @8 j# p& P2 f+ Y9 x! n% J% Y" K; k! W! U
        xlabel('x坐标');3 O) [0 r. ?. Z5 C" \4 s& N( G
        ylabel('y坐标');
    % \# P  x5 i: ], ]- c    zlabel('z坐标');% K2 F, D" P) k# u$ V) C2 N
        title(['当时间为' num2str(randnum(i)) '时的浓度分布情况']); ! O4 ]) |0 r. n5 p  X2 d. C
        colormap hsv
    8 a3 Z% ~0 u' q    colorbar
    0 y; z, {7 B2 \# k# Xend7 S; O/ ?0 \6 T; n5 Y. m
    回复

    使用道具 举报

    13ansys        

    0

    主题

    11

    听众

    42

    积分

    升级  38.95%

  • TA的每日心情
    奋斗
    2014-8-19 10:44
  • 签到天数: 9 天

    [LV.3]偶尔看看II

    自我介绍
    数学建模菜鸟
    回复

    使用道具 举报

    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2026-6-14 14:33 , Processed in 0.415522 second(s), 62 queries .

    回顶部