cellular programmWord文档下载推荐.docx
- 文档编号:5768628
- 上传时间:2023-05-05
- 格式:DOCX
- 页数:20
- 大小:17.47KB
cellular programmWord文档下载推荐.docx
《cellular programmWord文档下载推荐.docx》由会员分享,可在线阅读,更多相关《cellular programmWord文档下载推荐.docx(20页珍藏版)》请在冰点文库上搜索。
close;
number=uicontrol('
text'
1'
[20,400,50,20]);
%==============================
%CAsetup
n=128;
%initializethearrays
z=zeros(n,n);
cells=z;
sum=z;
%setafewcellstoone
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;
%howlongforeachcasetostabilityorsimpleoscillators
%buildanimageanddisplayit
imh=image(cat(3,cells,z,z));
set(imh,'
erasemode'
none'
)
axisequal
axistight
%indexdefinitionforcellupdate
x=2:
n-1;
y=2:
%Maineventloop
stop=0;
%waitforaquitbuttonpush
run=0;
%waitforadraw
freeze=0;
%waitforafreeze
while(stop==0)
if(run==1)
%nearestneighborsum
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);
%TheCArule
cells=(sum==3)|(sum==2&
cells);
%drawthenewimage
set(imh,'
cdata'
cat(3,cells,z,z))
%updatethestepnumberdiaplay
stepnumber=1+str2num(get(number,'
));
set(number,'
num2str(stepnumber))
end
if(freeze==1)
run=0;
freeze=0;
drawnow%needthisintheloopforcontrolstowork
end
2、SurfaceTension
clf;
clearall
%=============================================
n,y-1)+cells(x+1,y+1)+...
cells(x,y);
cells=~((sum<
4)|(sum==5));
3、PercolationCluster
%PercolationCluster
threshold=.63;
%
ax=axes('
units'
pixels'
[11500400],'
color'
k'
text('
[50,255,0],...
BioNB'
w'
fontname'
helvetica'
100)
[120,120,0],...
441'
initial=getframe(gca);
[a,b,c]=size(initial.cdata);
z=zeros(a,b);
cells=double(initial.cdata(:
:
1)==255);
visit=z;
imh=image(cat(3,z,cells,z));
%return
fori=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)+...
a-2,1:
b-2)+cells(1:
a-2,3:
a,1:
b-2)+cells(3:
a,3:
b);
pick=rand(a,b);
%edgesonly
%cells=(cells&
(sum<
8))|((sum>
=1)&
(pick>
=threshold)&
(visit==0));
cells=cells|((sum>
visit=(sum>
=1);
%&
(pick<
threshold);
cat(3,z,cells,z))
drawnow
return
figure
(2)
image(cat(3,z,cells,z))
4、Excitablemedia(BZreactionorheart)
%CAdriver
%excitablemedia
z=zeros(n,n);
cells=z;
.1;
%cells(n/2,n*.25:
n*.75)=1;
%cells(n*.25:
n*.75,n/2)=1;
sum=z;
x=[2:
n-1];
y=[2:
t=6;
%centervalue=6;
7makesfastpattern;
5analiatingwaves
t1=3;
%centervalue=3
1200
sum(x,y)=((cells(x,y-1)>
0)&
(cells(x,y-1)<
t))+((cells(x,y+1)>
(cells(x,y+1)<
t))+...
((cells(x-1,y)>
(cells(x-1,y)<
t))+((cells(x+1,y)>
(cells(x+1,y)<
((cells(x-1,y-1)>
(cells(x-1,y-1)<
t))+((cells(x-1,y+1)>
(cells(x-1,y+1)<
((cells(x+1,y-1)>
(cells(x+1,y-1)<
t))+((cells(x+1,y+1)>
(cells(x+1,y+1)<
t));
cells=((cells==0)&
(sum>
=t1))+...
2*(cells==1)+...
3*(cells==2)+...
4*(cells==3)+...
5*(cells==4)+...
6*(cells==5)+...
7*(cells==6)+...
8*(cells==7)+...
9*(cells==8)+...
0*(cells==9);
cat(3,z,cells/10,z))
drawnow
5、ForestFire
%forestfire
n=100;
Plightning=.000005;
Pgrowth=.01;
%.01
o=ones(n,n);
veg=z;
imh=image(cat(3,z,veg*.02,z));
%burning->
empty
%green->
burningifoneneigborburningorwithprob=f(lightning)
%empty->
greenwithprob=p(growth)
%veg={empty=0burning=1green=2}
3000
%nearbyfires?
sum=(veg(1:
n,[n1:
n-1])==1)+(veg(1:
n,[2:
n1])==1)+...
(veg([n1:
n-1],1:
n)==1)+(veg([2:
n1],1:
n)==1);
veg=...
2*(veg==2)-((veg==2)&
0|(rand(n,n)<
Plightning)))+...
2*((veg==0)&
rand(n,n)<
Pgrowth);
cat(3,(veg==1),(veg==2),z))
6、Gasdynamics
%HPP-gas
nx=52;
%mustbedivisibleby4
ny=100;
z=zeros(nx,ny);
o=ones(nx,ny);
sand=z;
sandNew=z;
gnd=z;
diag1=z;
diag2=z;
and12=z;
or12=z;
sums=z;
orsum=z;
gnd(1:
nx,ny-3)=1;
%rightgroundline
nx,3)=1;
%leftgroundline
gnd(nx/4:
nx/2-2,ny/2)=1;
%theholeline
gnd(nx/2+2:
nx,ny/2)=1;
gnd(nx/4,1:
ny)=1;
%topline
gnd(3*nx/4,1:
ny)=1;
%bottomline
%filltheleftside
r=rand(nx,ny);
sand(nx/4+1:
3*nx/4-1,4:
ny/2-1)=r(nx/4+1:
ny/2-1)<
0.3;
%sand(nx/4+1:
3*nx/4-1,ny*.75:
ny-4)=r(nx/4+1:
ny-4)<
0.75;
%sand(nx/2,ny/2)=1;
%sand(nx/2+1,ny/2+1)=1;
imh=image(cat(3,z,sand,gnd));
1000
p=mod(i,2);
%margolisneighborhood
%upperleftcellupdate
xind=[1+p:
2:
nx-2+p];
yind=[1+p:
ny-2+p];
%Seeifexactlyonediagonalisones
%only(atmost)oneofthefollowingcanbetrue!
diag1(xind,yind)=(sand(xind,yind)==1)&
(sand(xind+1,yind+1)==1)&
...
(sand(xind+1,yind)==0)&
(sand(xind,yind+1)==0);
diag2(xind,yind)=(sand(xind+1,yind)==1)&
(sand(xind,yind+1)==1)&
(sand(xind,yind)==0)&
(sand(xind+1,yind+1)==0);
%Thediagonalsbothnotoccupiedbytwoparticles
and12(xind,yind)=(diag1(xind,yind)==0)&
(diag2(xind,yind)==0);
%Onediagonalisoccupiedbytwoparticles
or12(xind,yind)=diag1(xind,yind)|diag2(xind,yind);
%foreverygasparticleseeifitneartheboundary
sums(xind,yind)=gnd(xind,yind)|gnd(xind+1,yind)|...
gnd(xind,yind+1)|gnd(xind+1,yind+1);
%celllayout:
%x,yx+1,y
%x,y+1x+1,y+1
%If(nowalls)and(diagonalsarebothnotoccupied)
%thenthereisnocollision,somoveoppositecelltocurrentcell
%If(nowalls)and(onlyonediagonalisoccupied)
%thenthereisacollisionsomoveccwcelltothecurrentcell
%If(awall)
%thendon'
tchangethecell(causesareflection)
sandNew(xind,yind)=...
(and12(xind,yind)&
~sums(xind,yind)&
sand(xind+1,yind+1))+...
(or12(xind,yind)&
sand(xind,yind+1))+...
(sums(xind,yind)&
sand(xind,yind));
sandNew(xind+1,yind)=...
sand(xind,yind+1))+...
sand(xind,yind))+...
sand(xind+1,yind));
sandNew(xind,yind+1)=...
sand(xind+1,yind))+...
sand(xind+1,yind+1))+...
sand(xind,yind+1));
sandNew(xind+1,yind+1)=...
sand(xind,yind))+...
sand(xind+1,yind))+...
sand(xind+1,yind+1));
sand=sandNew;
cat(3,z,sand,gnd))
8、Diffusionlimitedaggregation
%diffusion+dla
nx=200;
%mustb
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- cellular programm