- 在线时间
- 326 小时
- 最后登录
- 2019-9-18
- 注册时间
- 2014-8-5
- 听众数
- 36
- 收听数
- 9
- 能力
- 0 分
- 体力
- 4485 点
- 威望
- 0 点
- 阅读权限
- 60
- 积分
- 1854
- 相册
- 1
- 日志
- 0
- 记录
- 0
- 帖子
- 996
- 主题
- 413
- 精华
- 0
- 分享
- 3
- 好友
- 98
升级 85.4% TA的每日心情 | 开心 2019-9-18 21:55 |
---|
签到天数: 258 天 [LV.8]以坛为家I
群组: 2015国赛冲刺 群组: 2016美赛公益课程 群组: 国赛讨论 群组: 第三届数模基础实训 群组: Matlab讨论组 |
- <p>%渗流集群
- %Percolation Cluster
- clc
- clear
- close all
- threshold=0.63;
- %最初的图形是由图形操作决定的,以下程序设定坐
- %标系为一个固定的尺寸,在坐标系里写入文本,然后
- %获得并返回坐标系内容,并用getframe函数把它们
- %写入一个矩阵。
- ax=axes( 'units','pixels','position',[1,1,1050,400],'color','k' );%1050,400画布大小!
- text( 'units','pixels', 'position',[50 270 0],...
- 'string','ZhangWeiDong','color','w','fontname','helvetica',...
- 'fontsize',50 )
- text( 'units','pixels', 'position',[120 150 0],...
- 'string','1302034220','color','w','fontname','helvetica',...
- 'fontsize',50 )
- initial=getframe( gca );</p><p>[a,b,c]=size(initial.cdata);
- z=zeros(a,b);
- cells=double( initial.cdata(:,:,1)==255 );
- visit=z;
- sum=z;
- imh=image( cat(3,z,cells,z) );
- set( imh,'erasemode','none' )</p><p>for i=1:100
- sum(2:a-1,2:b-1)=cells(2:a-1,1:b-2)+cells(2:a-1,3:b)+...
- cells(1:a-2,2:b-1)+cells(3:a,2:b-1)+...
- cells(1:a-2,1:b-2)+cells(1:a-2,3:b)+...
- cells(3:a,1:b-2)+cells(3:a,3:b);
- pick=rand(a,b);
- %edge only
- cells=cells | ( (sum>=1)&(pick>=threshold)&(visit==0) );
- visit=(sum>=1);
-
- set( imh,'cdata',cat(3,z,cells,z) )
- drawnow
- %pause(0.5)
- end</p><p>return
- figure(2)
- image( cat(3,z,cells,z) )</p><p><p><p><p><p><p>
- </p>
复制代码
|
zan
|