Discrete Dynamics in Nature and Society

Volume 2014 (2014), Article ID 304297, 8 pages

http://dx.doi.org/10.1155/2014/304297

## Chaotic Behavior of One-Dimensional Cellular Automata Rule 24

^{1}Internet Data Center, Chongqing University of Science and Technology, Chongqing 401331, China^{2}School of Electrical and Information Engineering, Chongqing University of Science and Technology, Chongqing 401331, China^{3}Department of Mathematics and Information Engineering, Chongqing University of Education College, Chongqing 400065, China

Received 21 January 2014; Accepted 11 April 2014; Published 15 May 2014

Academic Editor: Zhen Jin

Copyright © 2014 Zujie Bie et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Wolfram divided the 256 elementary cellular automata rules informally into four classes using dynamical concepts like periodicity, stability, and chaos. Rule 24, which is Bernoulli -shift rule and is member of Wolfram’s class II, is said to be simple as periodic before. Therefore, it is worthwhile studying dynamical behaviors of four rules, whether they possess chaotic attractors or not. In this paper, the complex dynamical behaviors of rule 24 of one-dimensional cellular automata are investigated from the viewpoint of symbolic dynamics. We find that rule 24 is chaotic in the sense of both Li-Yorke and Devaney on its attractor. Furthermore, we prove that four rules of global equivalence of cellular automata are topologically conjugate. Then, we use diagrams to explain the attractor of rule 24, where characteristic function is used to describe the fact that all points fall into Bernoulli-shift map after two iterations under rule 24.

#### 1. Introduction

Cellular automata (CA) was first introduced by von Neumann in 1951 [1]. CA is a mathematical model consisting of large numbers of simple identical components with local interactions [2]. The simple components act together to produce complex global behavior. CA performs complex computation with high degree of efficiency and robustness. Three major factors have resulted in the revival of interest in the behavior of cellular systems [3]. First, the development of powerful computers and microprocessors has made the rapid simulation of CA possible. Second, the use of CA to simulate physical systems has attracted much interest in the scientific community. Third, the advent of VLSI as an implementation medium has focused attention on the communication requirements of successful hardware algorithms. Therefore, many researches about CA have been reported [4–6], especially in cryptography [7–9] and image processing [10, 11]. In addition, the analysis of dynamical behavior about dynamical system has aroused wide public concern [12–16], such as chaotic behavior [17, 18]. Thus, it is also important to analyze chaotic behavior of CA.

Here, we will only consider Boolean automata for which the cellular state . A cellular automaton consists of a number of cells which evolve by a simple local rule (identical rule). The value of each cell in the next stage is determined by the values of the cell and its neighbor cells in the current stage under the local rule [19]. The identical rule contained in each cell is essentially a finite state machine (FSM) which is specified in the form of a rule table [20]. A rule table contains an entry for possible neighborhood which consists of a cell and the adjacent cells. The cellular array is -dimensional, where is used in practice. In this paper, we will concentrate on . For a one-dimensional CA, a cell is connected to local neighbors (cells) on either side, where is referred to as the radius. A one-dimensional CA has cells linked in a line or in a circle. Denote the value in the th cell at the th stage by . For 2-state 3-neighborhood CA (), the evolution of th cell can be represented as a function of the present states of th, th, and th. The local function is a deterministic function to determine the next-stage value of the th cell, . For example, the rule 24 is a one-dimensional CA and its rule table is shown in Table 1. Thus, we have , , , , , , , and . Then, we find that the groups of three bits between parentheses represent all the possible neighborhood states and the single bit after equal sign is the resulting output bit at the next time step. The rule number is obtained by multiplying each output bit by the corresponding power of two and adding the results. Therefore, there exist possible rules with . In this paper, we mainly discuss rule 24. The rule 24 can be used to extract binary image edges. We know that bit strings 011 and 100 must be detected at the edges of binary images. If every row of an image is evolved according to the rule 24, every column of the image is evolved according to the rule 24 and their results do bitwise XOR (exclusive OR); then the image edges can be obtained. Therefore, it is worth researching the dynamical behaviors of rule 24.

In this paper, complex dynamical behaviors of rule 24 are studied in detail. We obtain two conditions according to Bernoulli -shift evolution for attractors under rule 24. The corresponding strongly connected graphs of transition matrices of determinative block systems are given. In terms of strongly connected subgraphs, we can predict the elements of attractor . Finally, we have proven that rule 24 is chaotic in the sense of Li-York and Devaney on sets of , which implies that the number of period orbits of rule 24 is infinite.

The rest of the paper is organized as follows. In Section 2, the intent of the paper and notations of symbolic dynamics are introduced. In Section 3, the dynamical behaviors of rule 24 are studied. In Section 4, we prove that four rules of global equivalence are topologically conjugate. In Section 5, characteristic function is used to describe the fact that all points fall into Bernoulli-shift map after two iterations under rule 24 and Lameray diagram is used to show clearly the iterative process of an attractor. Section 6 presents some conclusions.

#### 2. Preliminaries

In 1980s, Wolfram proposed CA as models for physical systems which exhibit complex or even chaotic behaviors based on empirical observations, and he divided the 256 elementary cellular automata (ECA) rules informally into four classes using dynamical concepts like periodicity, stability, and chaos [21–23]. However, some authors [24–28] found that some rules of Bernoulli -shift rules are chaotic in the sense of both Li-York and Devaney, where these rules were said to be simple as periodic by Wolfram. Rule 24 is belonging to Bernoulli -shift rules. Therefore, we need to research rule 24 and to find its some new dynamical properties.

In this paper, we will use some notations about CA as follows.

Chua et al. [29] mentioned that each rule has three globally equivalent local rules determined by three corresponding global transformations, namely, left-right transformation , global complementation , and left-right complementation . Each equivalence class is identified by , where is complexity index and is index of th class. In [30], the authors presented that 112 rules of 256 local rules were Bernoulli -shift rules. Each of the 112 Bernoulli -shift rules has an ID code , where denotes the number of attractors of rule , denotes the slope of the Bernoulli -shift map, and denotes the relevant forward time-. Hence, the space-time evolution of any one of the 112 rules on their attractors can be uniquely predicted by two parameters: and . For example, the attractors of rule 14 are () and ().

Some notations about symbolic dynamics can be referred to in [24–27]

It follows from [31] that the Boolean function of rule 24 is for all , , where “,” “−,” and “” stand for “AND,” “NOT,” and “XOR” logical operation, respectively. Sometimes, “” is omitted for simplicity. The truth table of Boolean function of rule 24 is shown in Table 1.

The subset, denoted by , is derived from the parameter of rule 24: , , ; that is,

#### 3. Dynamical Behaviors of

In the section, we investigate the complexity and chaotic dynamics of . In order to give our results, some definitions need be introduced.

*Definition 1 (see [32]). *A square matrix is irreducible, if for every pair of indices and there is an such that .

*Definition 2 (see [32]). *A square matrix is aperiodic, if there exists , such that , , for all .

*Definition 3 (see [32]). *Suppose that is a continuous mapping, where is a compact topological space. is said to be topologically mixing if, for any two open sets, , , such that , for all .

*Definition 4 (see [26]). *Let and be compact spaces; we say and are topologically conjugate if there is a homeomorphism , such that .

We give conditions according to Bernoulli -shift evolution for rule 24 as follows.

Theorem 5. *For rule 24, there exists a subset which satisfies if and only if, for all , , , and have the following relations. *(i)*If , then , ; , , ; , , .*(ii)*If , then , , , .*

*Proof. *Necessity: suppose that there exists a subset such that . Then for all , we have , for all .(i)If , then ; according to Table 1, we get , ; , , ; , , .(ii)If , then ; according to Table 1, we get , , , .

Sufficiency: suppose that there exists a subset and, for all , the relations between and satisfy the conditions (i) and (ii) in Theorem 5, for all .(i)If , we have .

Therefore,
(ii)If , we have

Hence, .

*Remark 6. *From the definition of subsystem, we know that are subsystems of .

The dynamical behaviors of on the set are shown as follows.

Let be a new state set, where , , , , and , , such that . Furthermore, subshift of is defined as . The transition matrix of the is

Obviously, is a square matrix. A square matrix corresponds to a directed graph. The vertices of the graph are the indices for the rows and columns of . There is an edge from vertex to vertex if . A square matrix is irreducible if and only if the corresponding graph is strongly connected. If is a two-order subshift of finite type, then it is topologically mixing if and only if is irreducible and aperiodic.

We give corresponding graph of the matrix in Figure 1, where vertices are the elements of set . It is obvious that is a strongly connected graph.

*Remark 7. *By definition, we know that is the determinative block system of and is a subshift of finite type.

*Remark 8. *Carefully observing Figure 1, we find that there are several strongly connected subgraphs: , , . The elements of will be composed by all vertices of subgraphs, respectively. For example, and is composed of vertices of subgraph ; then we have and vertices of subgraph will appear in , if , .

*Remark 9. *Let , , , and . In terms of Table 1, we find that, after one iteration, the possible sequences are generated by , , , and as follows: , or 001, , 001, 100 or 101, or 100.We find that, though , after one or two iterations, they belong to . Therefore, we guess that is the global attractor of .

*Theorem 10. is global attractor of .*

*Proof. *To prove that is the global attractor of , we consider two situations.(i). Since is invariant under , we have .(ii). Suppose that there exist such that . Then, we have , , , or by Theorem 5. Without loss of generality, we consider several situations as follows.(a)Let . Since , for all , we especially have
Therefore, it follows from Table 1 that by (7), which implies that . It is in contradiction to (8). Hence, is no ancestor. * *(b)Let . Since , for all , we especially have

Therefore, it follows from Table 1 that by (9). Furthermore, if , we know that is no ancestor. If , then only satisfies (11). However, we know that is no ancestor. The situations of and are similar to those of or . In conclusion, if , will be no ancestor or its ancestor will be no ancestor. Therefore, is global attractor of .

*Remark 11. *From the above conclusion, we find that all binary sequences will be bound to fall into after two iterations under .

*Based on the above definitions and analysis, we have the following results.*

*Theorem 12. (a) and are topologically conjugate.(b) is topologically mixing.(c) is topologically mixing.(d) The topological entropy .*

*Proof. *(a) We find out a homeomorphism from to .

Define
where , for all . In fact, by the definition of , we have , for all ; thus . Then, it is easy to check that is homeomorphism and . Therefore, and are topologically conjugate.

(b) Since and are topologically conjugate, we only need to check that the transition matrix of is irreducible and aperiodic. Actually, , for all . By Definitions 1 and 2, we know that is irreducible and aperiodic. So, in terms of [32, 33], is topologically mixing.

(c) Since and are equal in the set and is topologically mixing, is topologically mixing.

(d) The topological entropy on equals , where is the spectral radius of the transition matrix of the subshift . So, . Because two topological conjugate systems have the same topological entropy, the topological entropy of is equal to that of ; namely, .

*Theorem 13. is chaotic in the sense of both Li-Yorke and Devaney on .*

*Proof. *It follows from [33] that the positive topological entropy implies chaos in the sense of Li-Yorke and topologically mixing implies chaos in the sense of Li-Yorke and Devaney, since rule possesses very rich and complicated dynamical properties on .

*4. The Relationship between Four Rules of Global Equivalence Class *

*In this section, we will discuss the relationship between four rules of global equivalence class . In [31], rules 24, 66, 231, and 189 are partitioned into global equivalence class . Next, we prove that they are topologically conjugate to each other.*

*Theorem 14. (i) and are topologically conjugate.(ii) and are topologically conjugate.(iii) and are topologically conjugate.*

*Proof. *By [31], we have
Then we have
Therefore, we have , , . So, , , and are topologically conjugate to each other.

*Remark 15. *If there are two systems topologically conjugate, these two systems have the same dynamical properties. Rules , , , and are topologically conjugate, respectively. Therefore, if we know that one of four rules is chaotic in the sense of both Li-Yorke and Devaney in its attractors, we can deem that the other four rules are chaotic in the sense of both Li-Yorke and Devaney in their attractors, respectively.

*5. Using Diagrams to Explain Attractors of Four Rules*

*First, we give a definition on global characteristic function [30]. Given any local rule and binary configuration, for a ECA, where . Then we can uniquely associate the Boolean string with the binary expansion of a real number on the unit interval :
where is the decimal form of Boolean string . The ECA’ characteristic function of rule is defined as
where denotes rational numbers.*

*We choose . Figure 2 shows characteristic functions of rule 24. Figures 2(a), 2(b), and 2(c) describe the fact that all points fall into Bernoulli-shift map after two iterations under rule 24, which shows that is global attractor of . The phenomenon also shows that the forecast in Remark 11 is correct.*

*If we choose different values of for rule 24, we can get different initial binary configurations for the evolution of rule 24. The different initial binary configurations may lead to different attractor periods. If the value of is fixed, we may find the periods of attractors different. Given , we find three attractors of rule 24 shown in Figures 3(a), 3(b), and 3(c). Periods of these attractors are 1, 3, and 6, respectively.*

*Remark 16. *The sequence consists of arbitrary alternations of 0 and 1, and there are different possibilities of choice for .

*Next, we use Lameray diagram [34] to present the evolution process of attractor. In terms of the attractor of Figure 3(b), we get that the values of characteristic function of the attractor are 0.2813, 0.1406, and 0.5625, respectively. Figure 4 shows the iterative process of the attractor. Then, we can associate the period-3 attractor of rule 24 as a period-3 point of a continuous map which we know to be chaotic because “period-3 implies chaos” [35].*

*6. Conclusion*

*In this work, we characterize the global attractor of rule 24. We derive the conditions according to Bernoulli -shift evolution for attractor of rule 24. Then, in terms of the transition matrix of determinative block system of subsystem of rule 24, we obtain the value of topological entropy of subsystem. By corresponding strongly connected graph of transition matrix of determinative block system of subsystem , we guess that is the global attractor of . Furthermore, we prove that is the global attractor of . We find that rule 24 is topologically mixing on . Then, we prove that four rules of global equivalence of ECA are topologically conjugate. So, these four rules are chaotic in the sense of both Li-Yorke and Devaney on their attractors, respectively. We use diagrams to explain the attractor of rule 24, where characteristic function and Lameray diagram are used to describe the fact that all points fall into Bernoulli-shift map after two iterations and to show clearly the iterative process of an attractor, respectively.*

*Conflict of Interests*

*The authors declare that there is no conflict of interests regarding the publication of this paper.*

*Acknowledgments*

*This work was supported in part by the Research Project of Chongqing University of Science and Technology (CK2013B15), Scientific and Technological Research Program of Chongqing Municipal Education Commission (KJ131401, KJ131416), the Natural Science Foundation Project of CQCSTC (cstc2012jjB0095), the National Natural Science Foundation of China (51275547), and Achievement Transfer Program of Institutions of Higher Education in Chongqing (KJ121413).*

*References*

- J. von Neumann,
*The General and Logical Theory of Automata*, Pergamon Press, Lendon, UK, 1951. - H. Beigy and M. R. Meybodi, “Cellular learning automata with multiple learning automata in each cell and its applications,”
*IEEE Transactions on Systems, Man, and Cybernetics, Part B*, vol. 40, no. 1, pp. 54–65, 2010. View at Publisher · View at Google Scholar · View at Scopus - W. Pries, A. Thanailakis, and H. C. Card, “Group properties of cellular automata and VLSI applications,”
*IEEE Transactions on Computers*, vol. C-35, no. 12, pp. 1013–1024, 1986. View at Google Scholar · View at Scopus - Q.-X. Liu, Z. Jin, and M.-X. Liu, “Spatial organization and evolution period of the epidemic model using cellular automata,”
*Physical Review E*, vol. 74, no. 3, Article ID 031110, 6 pages, 2006. View at Publisher · View at Google Scholar · View at MathSciNet - G.-Q. Sun, Z. Jin, and L. Li, “Emergent turing pattern in epidemic spreading using cellular automaton,”
*International Journal of Modern Physics B*, vol. 25, no. 32, pp. 4605–4613, 2011. View at Publisher · View at Google Scholar · View at Scopus - G.-Q. Sun, Z. Jin, L.-P. Song, A. Chakraborty, and B.-L. Li, “Phase transition in spatial epidemics using cellular automata with noise,”
*Ecological Research*, vol. 26, no. 2, pp. 333–340, 2011. View at Publisher · View at Google Scholar · View at Scopus - S. Nandi, B. K. Kar, and P. Pal Chaudhuri, “Theory and applications of cellular automata in cryptography,”
*IEEE Transactions on Computers*, vol. 43, no. 12, pp. 1346–1357, 1994. View at Publisher · View at Google Scholar · View at MathSciNet - A. Abdo, S. Lian, I. Ismail, M. Amin, and H. Diab, “A cryptosystem based on elementary cellular automata,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 18, no. 1, pp. 136–147, 2013. View at Google Scholar - J.-C. Jeon and K.-Y. Yoo, “Elliptic curve based hardware architecture using cellular automata,”
*Mathematics and Computers in Simulation*, vol. 79, no. 4, pp. 1197–1203, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Z. Eslami, S. H. Razzaghi, and J. Z. Ahmadabadi, “Secret image sharing based on cellular automata and steganography,”
*Pattern Recognition*, vol. 43, no. 1, pp. 397–404, 2010. View at Publisher · View at Google Scholar · View at Scopus - R.-J. Chen and S.-J. Horng, “Novel SCAN-CA-based image security system using SCAN and 2-D von Neumann cellular automata,”
*Signal Processing: Image Communication*, vol. 25, no. 6, pp. 413–426, 2010. View at Publisher · View at Google Scholar · View at Scopus - L. Feng, X. Liao, Q. Han, and L. Song, “Modeling and analysis of peer-to-peer botnets,”
*Discrete Dynamics in Nature and Society*, vol. 2012, Article ID 865075, 18 pages, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - L. Feng, X. Liao, Q. Han, and H. Li, “Dynamical analysis and control strategies on malware propagation model,”
*Applied Mathematical Modelling*, vol. 37, no. 16-17, pp. 8225–8236, 2013. View at Publisher · View at Google Scholar · View at MathSciNet - L.-P. Song, Z. Jin, G.-Q. Sun, J. Zhang, and X. Han, “Influence of removable devices on computer worms: dynamic analysis and control strategies,”
*Computers & Mathematics with Applications*, vol. 61, no. 7, pp. 1823–1829, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - W. M. Van Ballegooijen and M. C. Boerlijst, “Emergent trade-offs and selection for outbreak frequency in spatial epidemics,”
*Proceedings of the National Academy of Sciences of the United States of America*, vol. 101, no. 52, pp. 18246–18250, 2004. View at Publisher · View at Google Scholar · View at Scopus - G.-Q. Sun, Q.-X. Liu, Z. Jin, A. Chakraborty, and B.-L. Li, “Influence of infection rate and migration on extinction of disease in spatial epidemics,”
*Journal of Theoretical Biology*, vol. 264, no. 1, pp. 95–103, 2010. View at Publisher · View at Google Scholar · View at MathSciNet - Q. Han, C. D. Li, and J. J. Huang, “Anticipating synchronization of chaotic systems with time delay and parameter mismatch,”
*Chaos*, vol. 19, Article ID 013104, 10 pages, 2009. View at Google Scholar - Q. Han, C. Li, and T. Huang, “Anticipating synchronization of a class of chaotic systems,”
*Chaos*, vol. 19, no. 2, Article ID 023105, 10 pages, 2009. View at Publisher · View at Google Scholar · View at MathSciNet - F. Bao, “Cryptanalysis of a partially known cellular automata cryptosystem,”
*IEEE Transactions on Computers*, vol. 53, no. 11, pp. 1493–1497, 2004. View at Publisher · View at Google Scholar · View at Scopus - M. Tomassini and M. Perrenoud, “Cryptography with cellular automata,”
*Applied Soft Computing*, vol. 1, no. 2, pp. 151–160, 2001. View at Google Scholar · View at Scopus - S. Wolfram, “Universality and complexity in cellular automata,”
*Physica D*, vol. 10, no. 1-2, pp. 1–35, 1984. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. Wolfram,
*Theory and Applications of Cellular Automata*, vol. 1 of*Advanced Series on Complex Systems*, World Scientific Publishing, Singapore, 1986. View at MathSciNet - S. Wolfram,
*A New Kind of Science*, Wolfram Media, Champaign, Ill, USA, 2002. View at MathSciNet - F.-Y. Chen, W.-F. Jin, G.-R. Chen, F.-F. Chen, and L. Chen, “Chaos of elementary cellular automata rule 42 of Wolfram's class II,”
*Chaos*, vol. 19, no. 1, Article ID 013140, 6 pages, 2009. View at Publisher · View at Google Scholar · View at MathSciNet - F.-F. Chen and F.-Y. Chen, “Complex dynamics of cellular automata rule 119,”
*Physica A*, vol. 388, no. 6, pp. 984–990, 2009. View at Publisher · View at Google Scholar · View at MathSciNet - L. Chen, F. Y. Chen, F. Chen, and W. Jin, “Complex symbolic dynamics of bernoulli shift cellular automata rule,” in
*Proceedings of the 9th International Conference for Young Computer Scientists (ICYCS '08)*, pp. 2868–2873, November 2008. View at Publisher · View at Google Scholar · View at Scopus - Q. Han, X. Liao, C. Li, and L. Feng, “Complex dynamics behaviors in cellular automata rule 35,”
*Journal of Cellular Automata*, vol. 6, no. 6, pp. 487–504, 2011. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Q. Han, X. Liao, and C. Li, “Complex dynamic behaviors in cellular automata rule 14,”
*Discrete Dynamics in Nature and Society*, Article ID 258309, 12 pages, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - L. O. Chua, V. I. Sbitnev, and S. Yoon, “A nonlinear dynamics perspective of Wolfram's new kind of science. III. Predicting the unpredictable,”
*International Journal of Bifurcation and Chaos*, vol. 14, no. 11, pp. 3689–3820, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - L. O. Chua, V. I. Sbitnev, and S. Yoon, “A nonlinear dynamics perspective of Wolfram's new kind of science. IV. From Bernoulli shift to $1/f$ spectrum,”
*International Journal of Bifurcation and Chaos*, vol. 15, no. 4, pp. 1045–1183, 2005. View at Publisher · View at Google Scholar · View at MathSciNet - J. Guan, S. Shen, C. Tang, and F. Chen, “Extending Chua's global equivalence theorem on Wolfram's new kind of science,”
*International Journal of Bifurcation and Chaos*, vol. 17, no. 12, pp. 4245–4259, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - B. P. Kitchens,
*Symbolic Dynamics: One-Sided, Two-Sided and Countable State Markov Shifts*, Universitext, Springer, Berlin, Germany, 1998. View at Publisher · View at Google Scholar · View at MathSciNet - Z. Zhou,
*Symbolic Dynamics*, Shanghai Scientific and Technological Education Publishing House, Shanghai, China, 1997, (Chinese). - L. P. Shilnikov, A. L. Shilnikov, D. V. Turaev, and L. O. Chua,
*Methods of Qualitative Theory in Nonlinear Dynamics. Part I*, vol. 4 of*World Scientific Series on Nonlinear Science. Series A: Monographs and Treatises*, World Scientific Publishing, River Edge, NJ, USA, 1998. View at Publisher · View at Google Scholar · View at MathSciNet - T. Y. Li and J. A. Yorke, “Period three implies chaos,”
*The American Mathematical Monthly*, vol. 82, no. 10, pp. 985–992, 1975. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet

*
*