24小时热门版块排行榜    

查看: 1963  |  回复: 5
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

zyj8119

木虫 (著名写手)

[交流] 【讨论】给个模拟气体扩散的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 ]
回复此楼

» 猜你喜欢

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

好好学习,天天向上。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

g1w2h31980

新虫 (初入文坛)


小木虫(金币+0.5):给个红包,谢谢回帖交流
请问和信息扩散有什么关联吗?

本文来自: 小木虫论坛 http://muchong.com/bbs/viewthread.php?tid=2750503&fpage=1
4楼2010-12-31 21:31:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

zyj8119

木虫 (著名写手)

引用回帖:
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,=0;B(:,1:h)=0;B(:,N-h+1:N)=0;B(M-h+1:M,=0;B(:, ...

请大家分析一下这个程序的运作规则,看的有点不太明白。
好好学习,天天向上。
2楼2010-12-29 22:37:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dzl1988

银虫 (初入文坛)


小木虫(金币+0.5):给个红包,谢谢回帖交流
请问 你有这个的 文章或者 资料么 对于changeP 这个函数 起什么作用额
3楼2010-12-30 23:58:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zyj8119

木虫 (著名写手)

引用回帖:
Originally posted by g1w2h31980 at 2010-12-31 21:31:40:
请问和信息扩散有什么关联吗?

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

信息扩散,大体上应该类似吧。
好好学习,天天向上。
5楼2010-12-31 23:53:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 299求调剂 +6 △小透明* 2026-03-17 6/300 2026-03-21 02:42 by JourneyLucky
[考研] 304求调剂 +6 曼殊2266 2026-03-18 6/300 2026-03-21 00:32 by JourneyLucky
[考研] 299求调剂 +4 shxchem 2026-03-20 5/250 2026-03-21 00:24 by 星空星月
[考研] 274求调剂 +10 S.H1 2026-03-18 10/500 2026-03-20 23:51 by JourneyLucky
[考研] 一志愿中海洋材料工程专硕330分求调剂 +8 小材化本科 2026-03-18 8/400 2026-03-20 23:16 by JourneyLucky
[考研] 330求调剂 +4 小材化本科 2026-03-18 4/200 2026-03-20 23:13 by JourneyLucky
[考研] 265求调剂 +12 梁梁校校 2026-03-19 13/650 2026-03-20 21:01 by 无际的草原
[考研] 295复试调剂 +8 简木ChuFront 2026-03-19 8/400 2026-03-20 20:44 by zhukairuo
[考研] 353求调剂 +3 拉钩不许变 2026-03-20 3/150 2026-03-20 19:56 by JourneyLucky
[考研] 一志愿吉林大学材料学硕321求调剂 +11 Ymlll 2026-03-18 15/750 2026-03-20 19:40 by 丁丁*
[考研] 281求调剂(0805) +14 烟汐忆海 2026-03-16 25/1250 2026-03-20 15:47 by yuncha
[考研] 招收调剂硕士 +4 lidianxing 2026-03-19 12/600 2026-03-20 12:25 by lidianxing
[考研] 能源材料化学课题组招收硕士研究生8-10名 +5 脱颖而出 2026-03-16 14/700 2026-03-20 09:30 by kkcoco25
[考研] 求调剂 +3 暗涌afhb 2026-03-16 3/150 2026-03-20 00:28 by 河南大学校友
[考研] 一志愿西安交通大学材料工程专业 282分求调剂 +5 枫桥ZL 2026-03-18 7/350 2026-03-19 14:52 by 功夫疯狂
[考研] 0703化学调剂 +5 pupcoco 2026-03-17 8/400 2026-03-19 13:58 by houyaoxu
[考研] 302求调剂 +4 小贾同学123 2026-03-15 8/400 2026-03-17 10:33 by 小贾同学123
[考研] 东南大学364求调剂 +5 JasonYuiui 2026-03-15 5/250 2026-03-16 21:28 by 木瓜膏
[考研] 304求调剂 +4 ahbd 2026-03-14 4/200 2026-03-16 16:48 by 我的船我的海
[考研] 289求调剂 +4 这么名字咋样 2026-03-14 6/300 2026-03-14 18:58 by userper
信息提示
请填处理意见