金属腐蚀的元胞自动机模拟代码

1、For i=3:102For j=3:102For k=3:102lf(C(i,j,k)=1 &V(IJ,K)=1)lfC(i-1,j,k)=1 | (i-2,j,k ) =1 &V (i-2,j,k ) =2 |(C (i-l,j-l,k ) =1 &v(i-l,j-l,k)=4)|(c(i-1,j+1,k)=1 &V(i-1,j+1,k) =3)|(C (i-1,j,k-1 ) =1 &V(i-1,j,k-1)=6)|(C(i-1,j,k+1)=1 &V(i-1,j,k+1) =5)D(i,j,k)=1;Else D(i-1,j,k)=1;E

2、ndElse if(C(i,j,k)=1 &V(i,j,k)=2)ifC(i+1,j,k)=1 |(C(i+2,j,k)=1 &V(i+2,j,k)=1)|(C (i+1,j-1,k ) =1 &V(i+1,j-1,k)=4) |(C(i+1,j+1,k)=1 &V(i+1,j+1,k) =3)|(C (i+1,j,k-1 ) =1 &V(i+1,j,k-1)=6) |(C(i+1,j,k+1)=1 &V(i+1,j,k+1) =5)D(i,j,k)=1;Else D(i+1,j,k)=1;EndElse if(C(i,j,k)=1 &V

3、(i,j,k)=3)ifC(i,j-1,k)=1 |(C(i-1,j- 1,k)=1 &V(i-1,j- 1,k)=2)|(C (i+1,j-1,k ) =1 &v(i+1,j-1,k)=1)|(C(i,j-2,k)=1 &V(i,j-2,k) =4)|(C (i,j-1,k-1 ) =1 &V(i,j-1,k+1)=6)|(C(i+1,j,k+1)=1 &V(i,j-1,k+1) =5)D(i,j,k)=1;Else D(i,j-1,k)=1;EndElse if(C(i,j,k)=1 &V(i,j,k)=4)ifC(i,j+1,k)=1 |(

4、C(i-l,j+l,k)=1 &V(i-l,j+l,k)=2)|(C (i+1,j+1,k ) =1&V(i+1,j+1,k)=1)|(C(i,j+2,k)=1 &V(i,j+2,k) =3)|(C (i,j+1,k-1 ) =1 &V(i,j-1,k+1)=6)|(C(i,j+1,k-1)=1 &V(i,j+1,k+1) =5)D(i,j,k)=1;Else D(i,j+1,k)=1;EndElseif(C(i,j,k)=1 &V(i,j,k)=5)If C(i,j,k-1)=1 |(C(i-1,j,k-1)=1 &V(i-1,j,k-

5、1)=2)|(C (i+1,j,k-1 ) =1&V(i+1,j,k-1)=1)|(C(i,j-1,k-1)=1 &V(i,j-1,k-1) =4)|(C (i,j+1,k-1 ) =1 &V(i,j+1,k-1)=3)|(C(i，j ,k-2)=1 &V(i,j,k-2) =6)D(i,j,k)=1;Else D(i,j,k-1)=1;EndElseif (C(i,j,k)=1 &V(i,j,k)=6)If C(i,j,k+1)=1 |(C(i-1,j,k+1)=1 &V(i-1,j,k+1)=2)|(C (i+1,j,k+1 ) =1 &V(i+1,j,k+1)=1)|(C(i,j-1,k+1)=1 &V(i,j-1,k+1) =4)|(C (i,j+1,k+1 ) =1 &V(i,j+1,k+1)=3) |(C(i,j,k+2)=1 &V(i,j,k+2) =5)D(i,