QQ登录

只需要一步,快速开始

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

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

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

1

主题

3

听众

19

积分

升级  14.74%

该用户从未签到

自我介绍
我我我

窝窝

Ooo
跳转到指定楼层
1#
发表于 2010-5-5 18:11 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
高斯气体分布2 d, Q7 G5 _% P3 M/ B
数值模拟& U: Z- T( U3 T0 K9 Z
作出动态分布图
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
    0 U4 T$ T: x1 z: y. f  H+ B- O%syms x y z t;
    1 i, C: l/ ?8 d2 k/ F6 _( j& X%Q=10;%初始浓度定为10mg/m^3;% i- x& K0 F. d
    %1.32*10^-5=1.32*10^-5;%扩散系数; c0 Z6 U- E( E9 h( Y
    %C=Q/(4*pi.*1.32*10^-5.*t)^1.5*exp(-(x.^2+y.^2+z.^2)./4.*1.32*10^-5.*t);
    / ^) a6 e# _% c0 P2 s7 w' o9 P%统一并申明一下单位) Q9 Z7 u! s6 X5 j' l# G, t4 F$ S
    %C (mg/m^3)7 ?. C* K8 j. l: f
    %Q (mg/m^3)+ I) {2 _1 P7 D# |( R- c
    %1.32*10^-5 (m^2/s)
    / K* O0 j$ ^* k5 i( v%t (s)
    1 A% i, R' O1 ?! W7 s5 o%x,y,z (m)' ]& a  E) ^1 A2 ^
    %申明一下扩散的条件如下:
    9 d  c9 g9 T; k+ b7 d$ H+ @& L. b* h%1.模拟的球形空间的直径D=3.5m 在x方向安装15个,在y方向安装3个
    + Z, D- o- t* a%2.当经过时间10s时,可视为扩散均匀,其中时间可以通过一次实验或者相关文献查得。) A+ t% f9 y+ x9 j8 L
    randnum=unifrnd(1,10,10,1);%产生一组符合均匀分布的时间
    9 h) C" F$ f) u  T8 \8 p# b! ?for i=1:4! {" x" d+ Y- T0 Q: M, d
        [x,y,z] = meshgrid(-1.75:.2:1.75,-1.75:.2:1.75,-1.75:.2:1.75);
    1 }' N" n$ H1 l3 x3 }& n! w) 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));
    , k4 \7 q0 N4 I- }    xslice=[-0.25,1.75];% R0 e8 B9 g  F( i
        yslice=[-0.25,1.25];2 s5 `3 _* T* L) N" b5 l
        zslice=0.25;
    ; [$ `" u  C% B' R- g' f: |    subplot(2,2,i);; s# Z* o( s& r3 G8 b
        slice(x,y,z,C,xslice,yslice,zslice);' }' C7 V$ U5 ]  `
        xlabel('x坐标');, O9 m4 S$ ~4 ?4 L7 g  Y) k
        ylabel('y坐标');  ~; h! y' [- b3 V
        zlabel('z坐标');
    % P' d( \& w& b8 m    title(['当时间为' num2str(randnum(i)) '时的浓度分布情况']);
    . ]7 R4 G/ g: u  t) u( t" M7 U' z& i5 x    colormap hsv
    ! x- S; L  o8 K, g    colorbar8 g$ c& X0 t) X/ M' w) W6 C
    end: {# w0 {" `# t* F7 B! ^
    回复

    使用道具 举报

    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-8-4 16:39 , Processed in 0.405055 second(s), 62 queries .

    回顶部