CA元胞自动机优化模型原代码.docx

CA优化模型原代码:M=load(‘d:\ca\jlwm’)N=load(‘d:\ca\jlwn.asc’)lindishy=load(‘d:\ca\ldfj3.asc’)caodishy=load(‘d:\ca\cdfj3.asc’)gengdishy=load(‘d:\ca\htfj3.asc’)[m,n]=size(M);Xr=[1 1 -1111 -1 -111;1111 -1 -1111 -1;-1111 -1 -1 -11 -1 -1;111111 -111I;l -1 -111 -1 -1 -111;1 -1 -11 -11 -11 -1 -1;-11 -1 -1 -1 -11 -1 -1 -1;-1111 -11 -11-1 -1;11-111-1 -1 -111;1 -1 -111-1 -1 -111];caodi=0;lindi=0;gengdi=0;for i=1:mforj=l:nifM(i,j)==4caodi=caodi+1;elseifM(i,j)==3lindi=lindi+1;elseifM(i,j)==2gengdi=gengdi+1;endendendfor i=1:mforj=1:nifM(i,j)==4iflindishy(i,j)gengdishy(i,j)iflindishy(i,j)caodishy(i,j)z=0;forP=max(1,i-1):min(i+1,m)forq=max(j-1,1):min(j+1,n)if(M(p,q)~=0)xr(M(p,q),3)==-1z=1;endendendif z== 0caodi=eaodi-1;M(i,j)=3;lindi=lindi+1;endelseiflindishy(i,j)==caodishy(i,j)caoditemp=0;linditemp=0;gengditemp=0;for p=max(1,i-1):min(i+1,m)for q=max(j-1,1):min(i+1,n)if N(p,q)==4caoditemp=caoditemp+1;elseifN(p,q)==3linditemp=linditemp+1;elseifN(p,q)=2gengditemp=gengditemp+1;endendendiflinditemp=max(caoditemp,gengditemp)z=0;for p=max(1,j-1):min(i+1,m)forq=max(j-1,1):min(j+1,n)if(M(p,q)~=0)xr(M(p,q),3)==-1;z=1;endendendifZ==0caodi=caodi-1;M(i,j)=3;lindi=lindi+1;endendendelseif lindishy(i,j)==gengdishy(i,j)if lindishy(i,j)caodishy(i,j)caoditemp=0:linditemp=0;gengditemp=0:for p=max(1,i-1):min(i+1,m)for q=max(j-1,1):min(j+1,n)ifN(p,q)==4caoditemp=caoditemp+1;elseifN(p,q)==3linditemp=linditemp+1;elseifN(p,q)==2gengditemp=gengdltemp+1;endendendiflinditemp=gengditempz=0;for p=max(1,j-1):min(i+1,m)forq=max(j-1,1):min(j+1,n)if(M(p,q)~=0)xr(M(p,q),3)==-1z=1;endendendif z==0caodi=caodi-1;M(i,j)=3;lindi=lindi+1;endendelseiflindishy(i,j)==caodishy(i,j)caoditemp=0;linditemp=0;gengditemp=0;for p=max(i,i-1):min(i+1,m)for q=max(j-1,1):min(j+1,n)if N(p,q)==4caoditemp=caoditemp+1;elseif N(p,q)==3linditemp=linditemp+1;elseifN(p,q)==2gengditemp=gengditemp+1;endendendif linditemp=max(caoditemp,gengditemp)z=0;for p=max(1,i-1):min(i+1,m)for q=max(j-1,1):min(j+1,n)if(M(p,q)~=0)xr(M(p,q),3)==-1z=1;endendendifz==0c