请选择 进入手机版 | 继续访问电脑版

QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 1153|回复: 0

死亡游戏(元胞自动机)代码求助

[复制链接]
字体大小: 正常 放大
690938218        

4

主题

13

听众

199

积分

升级  49.5%

  • TA的每日心情
    奋斗
    2016-9-22 13:29
  • 签到天数: 49 天

    [LV.5]常住居民I

    群组数学建模培训课堂2

    发表于 2014-8-27 20:29 |显示全部楼层
    |招呼Ta 关注Ta
    %Conway's life with GUI
    clf
    clear all
    %=============================================
    %build the GUI
    %define the plot button
    plotbutton=uicontrol('style','pushbutton',...
       'string','Run', ...
       'fontsize',12, ...
       'position',[100,400,50,20], ...
       'callback', 'run=1;');
    %define the stop button
    erasebutton=uicontrol('style','pushbutton',...
       'string','Stop', ...
       'fontsize',12, ...
       'position',[200,400,50,20], ...
       'callback','freeze=1;');
    %define the Quit button
    quitbutton=uicontrol('style','pushbutton',...
       'string','Quit', ...
       'fontsize',12, ...
       'position',[300,400,50,20], ...
       'callback','stop=1;close;');
    number = uicontrol('style','text', ...
        'string','1', ...
       'fontsize',12, ...
       'position',[20,400,50,20]);
    %=============================================
    %CA setup
    n=128;
    %initialize the arrays
    z = zeros(n,n);
    cells = z;
    sum = z;
    %set a few cells to one
    cells(n/2,.25*n:.75*n) = 1;
    cells(.25*n:.75*n,n/2) = 1;
    %cells(.5*n-1,.5*n-1)=1;
    %cells(.5*n-2,.5*n-2)=1;
    %cells(.5*n-3,.5*n-3)=1;
    cells = (rand(n,n))<.5 ;
    %how long for each case to stability or simple oscillators
    %build an image and display it
    imh = image(cat(3,cells,z,z));
    set(imh, 'erasemode', 'none')
    axis equal
    axis tight
    %index definition for cell update
    x = 2:n-1;
    y = 2:n-1;
    %Main event loop
    stop= 0; %wait for a quit button push
    run = 0; %wait for a draw
    freeze = 0; %wait for a freeze
    while (stop==0)
        if (run==1)
            %nearest neighbor sum
            sum(x,y) = cells(x,y-1) + cells(x,y+1) + ...
                cells(x-1, y) + cells(x+1,y) + ...
                cells(x-1,y-1) + cells(x-1,y+1) + ...
                cells(3:n,y-1) + cells(x+1,y+1);
            % The CA rule
            cells = (sum==3) | (sum==2 & cells);      
            %draw the new image
            set(imh, 'cdata', cat(3,cells,z,z) )
            %update the step number diaplay
            stepnumber = 1 + str2num(get(number,'string'));
            set(number,'string',num2str(stepnumber))
        end
        if (freeze==1)
            run = 0;
            freeze = 0;
        end
        drawnow  %need this in the loop for controls to work
    end
    能不能解释下 cells = (sum==3) | (sum==2 & cells);和三个set?
    谢谢
    zan
    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

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

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

    蒙公网安备 15010502000194号

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

    GMT+8, 2024-3-29 22:48 , Processed in 0.428784 second(s), 49 queries .

    回顶部