24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

zyj8119

木虫 (著名写手)

应助: 65 (初中生)
贵宾: 0.003
金币: 915.1
散金: 1440
红花: 35
帖子: 2936
在线: 1329.4小时
虫号: 664177
注册: 2008-11-29
性别: GG
专业: 理论和计算化学

[交流] 【讨论】给个模拟气体扩散的MATLAB代码已有4人参与

CODE:

function gas_diffusive

clc;

rand('state',0);set(gcf,'Doublebuffer','on');

N=300;M=120;h=2;xp=round(N/2);D=round(M/12);

B=ones(M,N);B(1:h,:)=0;B(:,1:h)=0;B(:,N-h+1:N)=0;B(M-h+1:M,:)=0;B(:,xp:...

    +h-1)=0;

B(round(M/2)-D+1:round(M/2)+D,xp:xp+h-1)=1;

S=zeros(M,N);S(h+1:M-h,h+1:xp-1)=rand(M-2*h,xp-h-1)>0.5;

Tc=700;C=min(1-S,B);subplot(211);asa=imshow(C,[]);

ti=title('time=0');

rL=sum(sum(S(:,1:xp)))/((xp-h)*(M-2*h));xR=sum(sum(S(:,xp+h:N)))/((N-xp...

-2*h+1)*(M-2*h));

subplot(212),p1=plot(1,rL,'r-');xlim([1,Tc]);hold on;p2=plot(1,xR,'k--');

xlabel('时间');ylabel('时间');legend('左侧密度','右侧密度');

for k=1:Tc

    p=mod(k,2);x1=h+1+p:2:xp-3+p;y1=h+1+p:2:M-h-2+p;x2=xp-1+p:2:xp+h-1+p;

    y2=round(M/2)-D+2:2:round(M/2)+D-1;x3=xp+h+1+p:2:N-h-2+p;y3=h+1+p:2:...

M-h-2+p;

    St=S;St(y1,x1)=S(y1+1,x1+1);St(y1+1,x1+1)=S(y1,x1);St(y1,x1+1)=S(...

    y1+1,x1);

    St(y1+1,x1)=S(y1,x1+1);St(y2,x2)=S(y2+1,x2+1);St(y2+1,x2+1)=S(y2,x2);

    St(y2+1,x2)=S(y2,x2+1);St(y2,x2+1)=S(y2+1,x2+1);St(y3,x3)=S(y3+1,x3+1);

    St(y3+1,x3+1)=S(y3,x3);St(y3+1,x3)=S(y3,x3+1);St(y3,x3+1)=S(y3+1,x3);

    St=changep(St,S,x1,y1);St=changep(St,S,x2,y2);St=changep(St,S,x3,y3);

    S=St;

    set(ti,'string',['time=',num2str(k)]);set(asa,'CData',min(1-S,B));

    xL(k+1)=sum(sum(S(:,1:xp)))/((xp-h)*(M-2*h));

    xR(k+1)=sum(sum(S(:,xp+h:N)))/((N-xp-2*h+1)*(M-2*h));

    set(p1,'XData',1:k+1,'YData',xL);set(p2,'XData',1:k+1,'YData',xR);

    pause(0.2);

end

function St=changep(St,S,x,y)

Sr1=S(y+1,x)+S(y,x);Sc1=S(y,x+1)+S(y+1,x+1);Sr2=S(y,x+1)+S(y,x);

Sc2=S(y+1,x)+S(y+1,x+1);[p,q]=find(Sr1==1&Sc1==1&Sr2==1&Sc2==1);

for k=1:length(p)

    St(y(p(k)),x(q(k)))=S(y(p(k))+1,x(q(k)));St(y(p(k))+1,x(q(k)))...

        =S(y(p(k))+1,x(q(k)));

    St(y(p(k)),x(q(k))+1)=S(y(p(k))+1,x(q(k))+1);St(y(p(k))+1,x(q(k))+1)...

        =S(y(p(k))+1,x(q(k))+1);

end

[ Last edited by robert2020 on 2010-12-31 at 14:23 ]

回复此楼

» 猜你喜欢

参与限项已经有3人回复
假如你的研究生提出不合理要求已经有7人回复
实验室接单子已经有4人回复
全日制（定向）博士已经有4人回复
对氯苯硼酸纯化已经有3人回复
求助:我三月中下旬出站，青基依托单位怎么办？已经有12人回复
不自信的我已经有12人回复
所感已经有4人回复
要不要辞职读博？已经有7人回复
北核录用已经有3人回复

» 本主题相关价值贴推荐，对您同样有帮助:

求助用MS模拟气体在分子筛中扩散系数时碰到的问题已经有20人回复
三维凝固模拟的时候，怎么用matlab画成三维体啊已经有14人回复
如何通过matlab把文件夹中的部分文件拷贝到另一个文件夹中? 已经有4人回复
问个关于matlab里find函数的问题！已经有4人回复
求回归方程参数估计的最优化算法matlab代码已经有9人回复
如何用MATLAB 实现化学反应方程式（写程序代码）? 已经有20人回复
请教个matlab中的rand函数已经有3人回复
关于matlab 2010b 的一个小问题已经有7人回复
【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题已经有21人回复
【求助】一个matlab窗口占25%CPU，如何提高？已经有11人回复
【求助】问一个关于matlab中legend的问题已经有11人回复
【求助】Weickert的各向异性扩散方程的滤波方法的matlab程序代码已经有8人回复

好好学习，天天向上。

1楼 2010-12-29 22:36:45

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

g1w2h31980

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 87.5
帖子: 3
在线: 6.2小时
虫号: 1076295
注册: 2010-08-17

★
小木虫(金币+0.5):给个红包，谢谢回帖交流

请问和信息扩散有什么关联吗？

本文来自: 小木虫论坛 http://muchong.com/bbs/viewthread.php?tid=2750503&fpage=1

赞一下(1人)

回复此楼

4楼2010-12-31 21:31:40

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 6 个回答

zyj8119

木虫 (著名写手)

应助: 65 (初中生)
贵宾: 0.003
金币: 915.1
散金: 1440
红花: 35
帖子: 2936
在线: 1329.4小时
虫号: 664177
注册: 2008-11-29
性别: GG
专业: 理论和计算化学

引用回帖:

Originally posted by zyj8119 at 2010-12-29 22:36:45:
[code]function gas_diffusive
clc;
rand('state',0);set(gcf,'Doublebuffer','on');
N=300;M=120;h=2;xp=round(N/2);D=round(M/12);
B=ones(M,N);B(1:h,