元胞自动机程序.txt
《元胞自动机程序.txt》由会员分享,可在线阅读,更多相关《元胞自动机程序.txt(16页珍藏版)》请在冰豆网上搜索。
![元胞自动机程序.txt](https://file1.bdocx.com/fileroot1/2022-10/7/0910e4fb-615f-4c5b-bd55-1a89214fa161/0910e4fb-615f-4c5b-bd55-1a89214fa1611.gif)
��Matlabʵ��Ԫ���Զ����������ռ���ת�أ�����2010-03-1313:
04:
28�Ķ�121����0�ֺţ�����С����.
file:
life.m
%%��ʼ��
m=50;
X=zeros(m,m);
X(25,25)=1;
n=[m1:
m-1];
e=[2:
m1];
s=[2:
m1];
w=[m1:
m-1];
%���Ƴ�ʼͼ��
[i,j]=find(X);
figure(gcf);
plothandle=plot(i,j,'.',...
'Color','blue',...
'MarkerSize',12);
axis([0m+10m+1]);
%%�ݻ�
fork=1:
50
%�ھ���
N=X(n,:
)+X(s,:
)+X(:
e)+X(:
w)+...
X(n,e)+X(n,w)+X(s,e)+X(s,w);
%������
RAND=rand(m);
%����
X=X|(N.*RAND>0.99);
%��ͼ
[i,j]=find(X);
set(plothandle,'xdata',i,'ydata',j)
drawnow
pause(0.2)
k
end
file2:
functionsierpinski(n);
%ʹ��Ԫ���Զ�������sierpinskiֱ�ǵ�Ƭ
%Example:
%sierpinski(256);
%%�㷨��:
�ﲩ��,�������㷨�������ƣ���VisualC++ʵ�֡�
ifnargin==0;
n=256;
end
X=ones(n);
X(1,n-1)=0;
H=imshow(X,[]);
set(gcf,'doublebuffer','on');
k=1;
whilekX(k+1,1:
end-1)=xor(X(k,1:
end-1),X(k,2:
end));
X(k+1,n)=1;
set(H,'CData',X);
pause(0.1);
k=k+1;
end
file3:
functionCA_sim_cloud;
%ʹ��Ԫ���Զ���ģ��������ǵ���ͼ
%
%reference:
%Piazza,E.;Cuccoli,F.;
%CellularAutomataSimulationofCloudsinSatelliteImages,
%GeoscienceandRemoteSensingSymposium,2001.IGARSS'01.
%IEEE2001InternationalVolume4,9-13July2001Page(s):
%1722-1724vol.4DigitalObjectIdentifier10.1109/IGARSS.
%2001.977050
time=888;%����ִ�в���
M=240;
N=320;
S=round(rand(M,N)*15);
p=[1,2,1,6,6,1,2,1];
p=sum(tril(meshgrid(p)),2)/20;
rand('state',0);
SS=S;
R=rand(M,N);
G=R;
B=R;
C=cat(3,R,G,B);
fig=figure;
set(fig,'DoubleBuffer','on');
mov=avifile('example2.avi');
cc=imshow(C,[]);
set(gcf,'Position',[13355157194])
x1=(1:
3)+round(M/2);y1=(1:
3)+round(N/3);
x2=(1:
3)+round(M/3);y2=(1:
3)+round(N/2);
x3=(1:
3)+round(M/1.5);y3=(1:
3)+round(N/2);
q=0;
qq=15/4;
whileq