X-Git-Url: https://bilbo.iut-bm.univ-fcomte.fr/and/gitweb/chaos1.git/blobdiff_plain/e041c0c5dd32ff38f533543cc17b6b2e8ce090ce..8f2ca538d4cca28cafc98b78ce0c0e697085f3ee:/main.tex diff --git a/main.tex b/main.tex index aa920cb..545524f 100644 --- a/main.tex +++ b/main.tex @@ -65,14 +65,15 @@ Many research works deal with chaotic neural networks for various fields of application. Unfortunately, up to now these networks are usually claimed to be chaotic without any mathematical proof. The purpose of this paper is to establish, based on a rigorous theoretical -framework, an equivalence between chaotic iterations according to the -Devaney's formulation of chaos and a particular class of neural +framework, an equivalence between chaotic iterations according to +Devaney and a particular class of neural networks. On the one hand we show how to build such a network, on the other hand we provide a method to check if a neural network is a chaotic one. Finally, the ability of classical feedforward multilayer perceptrons to learn sets of data obtained from a chaotic dynamical system is regarded. Various Boolean functions are iterated on finite -states, some of them are proven to be chaotic as it is defined by +states. Iterations of some of them are proven to be chaotic + as it is defined by Devaney. In that context, important differences occur in the training process, establishing with various neural networks that chaotic behaviors are far more difficult to learn. @@ -85,6 +86,7 @@ behaviors are far more difficult to learn. \maketitle \begin{quotation} + Chaotic neural networks have received a lot of attention due to the appealing properties of deterministic chaos (unpredictability, sensitivity, and so on). However, such networks are often claimed @@ -105,70 +107,71 @@ is far more difficult than non chaotic behaviors. \section{Introduction} \label{S1} -Several research works have proposed or used chaotic neural networks -these last years. This interest comes from their complex dynamics and -the various potential application areas. Chaotic neural networks are -particularly considered to build associative memories -\cite{Crook2007267} and digital security tools like hash functions -\cite{Xiao10}, digital watermarking \cite{1309431,Zhang2005759}, or -cipher schemes \cite{Lian20091296}. In the first case, the background -idea is to control chaotic dynamics in order to store patterns, with -the key advantage of offering a large storage capacity. For the -computer security field, the use of chaotic dynamics is motivated by -their unpredictability and random-like behaviors. Indeed, -investigating new concepts is crucial in this field, because new -threats are constantly emerging. As an illustrative example, the -former standard in hash functions, namely the SHA-1 algorithm, has -been recently weakened after flaws were discovered. +Several research works have proposed or run chaotic neural networks +these last years. The complex dynamics of such a networks leads to +various potential application areas: associative +memories~\cite{Crook2007267} and digital security tools like hash +functions~\cite{Xiao10}, digital +watermarking~\cite{1309431,Zhang2005759}, or cipher +schemes~\cite{Lian20091296}. In the former case, the background idea +is to control chaotic dynamics in order to store patterns, with the +key advantage of offering a large storage capacity. For the latter +case, the use of chaotic dynamics is motivated by their +unpredictability and random-like behaviors. Thus, investigating new +concepts is crucial in this field, because new threats are constantly +emerging. As an illustrative example, the former standard in hash +functions, namely the SHA-1 algorithm, has been recently weakened +after flaws were discovered. Chaotic neural networks have been built with different approaches. In the context of associative memory, chaotic neurons like the nonlinear -dynamic state neuron \cite{Crook2007267} are frequently used to build -the network. These neurons have an inherent chaotic behavior, which is -usually assessed through the computation of the Lyapunov exponent. An -alternative approach is to consider a well-known neural network -architecture: the MultiLayer Perceptron. MLP networks are basically -used to model nonlinear relationships between data, due to their -universal approximator capacity. Thus, this kind of networks can be -trained to model a physical phenomenon known to be chaotic such as +dynamic state neuron \cite{Crook2007267} frequently constitute the +nodes of the network. These neurons have an inherent chaotic behavior, +which is usually assessed through the computation of the Lyapunov +exponent. An alternative approach is to consider a well-known neural +network architecture: the MultiLayer Perceptron (MLP). These networks +are suitable to model nonlinear relationships between data, due to +their universal approximator capacity. Thus, this kind of networks can +be trained to model a physical phenomenon known to be chaotic such as Chua's circuit \cite{dalkiran10}. Sometimes, a neural network which -is build using transfer functions and initial conditions that are both +is build by combining transfer functions and initial conditions that are both chaotic, is itself claimed to be chaotic \cite{springerlink:10.1007/s00521-010-0432-2}. What all of these chaotic neural networks have in common is that they are claimed to be chaotic despite a lack of any rigorous mathematical -proof. The purpose of this paper is to fill this gap, using a -theoretical framework based on the Devaney's definition of chaos +proof. The first contribution of this paper is to fill this gap, +using a theoretical framework based on the Devaney's definition of chaos \cite{Devaney}. This mathematical theory of chaos provides both qualitative and quantitative tools to evaluate the complex behavior of a dynamical system: ergodicity, expansivity, and so on. More precisely, in this paper, which is an extension of a previous work -\cite{bgs11:ip}, we establish an equivalence between chaotic -iterations and a class of globally recurrent multilayer perceptrons. -We investigate the converse problem too, that is, the ability for -classical MultiLayer Perceptrons (MLP) to learn a particular family of +\cite{bgs11:ip}, we establish the equivalence between chaotic +iterations and a class of globally recurrent MLP. +The investigation the converse problem is the second contribution: +we indeed study the ability for +classical MultiLayer Perceptrons to learn a particular family of discrete chaotic dynamical systems. This family, called chaotic iterations, is defined by a Boolean vector, an update function, and a -sequence giving the component to update at each iteration. It has -been previously established that such dynamical systems can behave -chaotically, as it is defined by Devaney, when the chosen function has -a strongly connected iterations graph. In this document, we will +sequence giving which component to update at each iteration. It has +been previously established that such dynamical systems is +chaotically iterated (as it is defined by Devaney) when the chosen function has +a strongly connected iterations graph. In this document, we experiment several MLPs and try to learn some iterations of this kind. -We will show that non-chaotic iterations can be learned, whereas it is +We show that non-chaotic iterations can be learned, whereas it is far more difficult for chaotic ones. That is to say, we have discovered at least one family of problems with a reasonable size, -such that artificial neural networks should not be applied in the -presence of chaos, due to their inability to learn chaotic behaviors -in this context. +such that artificial neural networks should not be applied +due to their inability to learn chaotic behaviors in this context. The remainder of this research work is organized as follows. The next -section is devoted to the basics of chaotic iterations and Devaney's +section is devoted to the basics of Devaney's chaos. Section~\ref{S2} formally describes how to build a neural -network that operates chaotically. The following two~sections are -devoted to the dual case: checking whether an existing neural network -is chaotic or not (Section \ref{S3}), and discussion on topological -properties of chaotic neural networks (Section~\ref{S4}). The +network that operates chaotically. Section~\ref{S3} is +devoted to the dual case of checking whether an existing neural network +is chaotic or not. +Topological properties of chaotic neural networks +are discussed in Sect.~\ref{S4}. The Section~\ref{section:translation} shows how to translate such iterations into an Artificial Neural Network (ANN), in order to evaluate the capability for this latter to learn chaotic behaviors. @@ -179,14 +182,14 @@ system. Prediction success rates are given and discussed for the two sets. The paper ends with a conclusion section where our contribution is summed up and intended future work is exposed. -\section{Link between Chaotic Iterations and Devaney's Chaos} +\section{Chaotic Iterations according to Devaney} In this section, the well-established notion of Devaney's mathematical -chaos is firstly presented. Preservation of the unpredictability of +chaos is firstly recalled. Preservation of the unpredictability of such dynamical system when implemented on a computer is obtained by -using some discrete iterations called ``chaotic iterations'', which -are thus introduced. The result establishing the link between chaotic -iterations and Devaney's chaos is finally recalled at the end of this +using some discrete iterations called ``asynchronous iterations'', which +are thus introduced. The result establishing the link between such +iterations and Devaney's chaos is finally presented at the end of this section. In what follows and for any function $f$, $f^n$ means the composition @@ -210,7 +213,7 @@ mathematically this kind of unpredictability is also referred to as deterministic chaos. For example, many weather forecast models exist, but they give only suitable predictions for about a week, because they are initialized with conditions that reflect only a partial knowledge -of the current weather. Even if initially the differences are small, +of the current weather. Even the differences are initially small, they are amplified in the course of time, and thus make difficult a long-term prediction. In fact, in a chaotic system, an approximation of the current state is a quite useless indicator for predicting @@ -219,8 +222,8 @@ future states. From mathematical point of view, deterministic chaos has been thoroughly studied these last decades, with different research works that have provide various definitions of chaos. Among these -definitions, the one given by Devaney~\cite{Devaney} is well -established. This definition consists of three conditions: +definitions, the one given by Devaney~\cite{Devaney} is +well-established. This definition consists of three conditions: topological transitivity, density of periodic points, and sensitive point dependence on initial conditions. @@ -231,13 +234,17 @@ given region. More precisely, \begin{definition} \label{def2} A continuous function $f$ on a topological space $(\mathcal{X},\tau)$ is defined to be {\bf topologically transitive} if for any pair of -open sets $U$, $V \in \mathcal{X}$ there exists $k>0$ such that +open sets $U$, $V \in \mathcal{X}$ there exists +$k \in +\mathds{N}^{\ast}$ + such that $f^k(U) \cap V \neq \emptyset$. \end{definition} This property implies that a dynamical system cannot be broken into -simpler subsystems. It is intrinsically complicated and cannot be -simplified. On the contrary, a dense set of periodic points is an +simpler subsystems. +Intuitively, its complexity does not allow any simplification. +On the contrary, a dense set of periodic points is an element of regularity that a chaotic dynamical system has to exhibit. \begin{definition} \label{def3} @@ -247,7 +254,7 @@ A point $x$ is called a {\bf periodic point} for $f$ of period~$n \in \begin{definition} \label{def4} $f$ is said to be {\bf regular} on $(\mathcal{X},\tau)$ if the set of - periodic points for $f$ is dense in $\mathcal{X}$ ($\forall x \in + periodic points for $f$ is dense in $\mathcal{X}$ ( for any $x \in \mathcal{X}$, we can find at least one periodic point in any of its neighborhood). \end{definition} @@ -261,39 +268,40 @@ whole system. Then, $f$ has {\bf sensitive dependence on initial conditions} if there exists $\delta >0$ such that, for any $x\in \mathcal{X}$ and any neighborhood $V$ of $x$, there exist $y\in V$ and $n > 0$ such that - $d\left(f^{n}(x), f^{n}(y)\right) >\delta $. $\delta$ is called the + $d\left(f^{n}(x), f^{n}(y)\right) >\delta $. The value $\delta$ is called the {\bf constant of sensitivity} of $f$. \end{definition} Finally, \begin{definition} \label{def5} -$f$ is {\bf chaotic according to Devaney} on $(\mathcal{X},\tau)$ if - $f$ is regular, topologically transitive, and has sensitive - dependence to its initial conditions. +The dynamical system that iterates $f$ is {\bf chaotic according to Devaney} +on $(\mathcal{X},\tau)$ if $f$ is regular, topologically transitive, +and has sensitive dependence to its initial conditions. \end{definition} Let us notice that for a metric space the last condition follows from the two first ones~\cite{Banks92}. -\subsection{Chaotic Iterations} +\subsection{Asynchronous Iterations} -This section presents some basics on topological chaotic iterations. +%This section presents some basics on topological chaotic iterations. Let us firstly discuss about the domain of iteration. As far as we -know, no result rules that the chaotic behavior of a function that has -been theoretically proven on $\R$ remains valid on the floating-point +know, no result rules that the chaotic behavior of a dynamical system +that has been theoretically proven on $\R$ remains valid on the +floating-point numbers, which is the implementation domain. Thus, to avoid loss of chaos this work presents an alternative, that is to iterate Boolean maps: results that are theoretically obtained in that domain are preserved in implementations. Let us denote by $\llbracket a ; b \rrbracket$ the following interval -of integers: $\{a, a+1, \hdots, b\}$, where $a~<~b$. A {\it system} +of integers: $\{a, a+1, \hdots, b\}$, where $a~<~b$. +In that section, a system under consideration iteratively modifies a collection of $n$~components. Each component $i \in \llbracket 1; n \rrbracket$ takes its value $x_i$ among the domain $\Bool=\{0,1\}$. A~{\it - configuration} of the system at discrete time $t$ (also said at {\it - iteration} $t$) is the vector + configuration} of the system at discrete time $t$ is the vector %\begin{equation*} $x^{t}=(x_1^{t},\ldots,x_{n}^{t}) \in \Bool^n$. %\end{equation*} @@ -312,11 +320,24 @@ indeed the negation of $x_i$). Intuitively, $x$ and $N(i,x)$ are neighbors. The discrete iterations of $f$ are represented by the oriented {\it graph of iterations} $\Gamma(f)$. In such a graph, vertices are configurations of $\Bool^n$ and there is an arc labeled -$i$ from $x$ to $N(i,x)$ iff $f_i(x)$ is $N(i,x)$. +$i$ from $x$ to $N(i,x)$ if and only if $f_i(x)$ is $N(i,x)$. In the sequel, the {\it strategy} $S=(S^{t})^{t \in \Nats}$ is the -sequence defining the component to update at time $t$ and $S^{t}$ -denotes its $t-$th term. We introduce the function $F_{f}$ that is +sequence defining which component to update at time $t$ and $S^{t}$ +denotes its $t-$th term. +This iteration scheme that only modifies one element at each iteration +is clasically refered as \emph{asynchronous iterations}. +Next section shows the link between asynchronous iterations and +Devaney's Chaos. + +\subsection{On the link between asynchronous iterations and + Devaney's Chaos} + +In this subsection we recall the link we have established between +asynchronous iterations and Devaney's chaos. The theoretical framework is +fully described in \cite{guyeux09}. + +We introduce the function $F_{f}$ that is defined for any given application $f:\Bool^{n} \to \Bool^{n}$ by $F_{f}: \llbracket1;n\rrbracket\times \mathds{B}^{n} \rightarrow \mathds{B}^{n}$, s.t. @@ -357,17 +378,10 @@ X^{k+1}& = & G_{f}(X^{k})\\ where $\sigma$ is the function that removes the first term of the strategy ({\it i.e.},~$S^0$). This definition means that only one component of the system is updated at an iteration: the $S^t$-th -element. But it can be extended by considering subsets for $S^t$. +element. But it can be extended by considering subsets for $S^t$. -Let us finally remark that, despite the use of the term {\it chaotic}, -there is {\it priori} no connection between these iterations and the -mathematical theory of chaos presented previously. -\subsection{Chaotic Iterations and Devaney's Chaos} -In this subsection we recall the link we have established between -chaotic iterations and Devaney's chaos. The theoretical framework is -fully described in \cite{guyeux09}. The {\bf distance} $d$ between two points $(S,x)$ and $(\check{S},\check{x})\in \mathcal{X} = \llbracket1;n\rrbracket^\Nats @@ -389,12 +403,12 @@ d_{s}(S,\check{S})=\frac{9}{2n}\sum_{t=0}^{\infty This distance is defined to reflect the following information. Firstly, the more two systems have different components, -the larger the distance between them. Secondly, two systems with +the larger the distance between them is. Secondly, two systems with similar components and strategies, which have the same starting terms, must induce only a small distance. The proposed distance fulfill these requirements: on the one hand its floor value reflects the difference between the cells, on the other hand its fractional part -measure the difference between the strategies. +measures the difference between the strategies. The relation between $\Gamma(f)$ and $G_f$ is clear: there exists a path from $x$ to $x'$ in $\Gamma(f)$ if and only if there exists a @@ -403,13 +417,13 @@ initial point $(s,x)$ reaches the configuration $x'$. Using this link, Guyeux~\cite{GuyeuxThese10} has proven that, \begin{theorem}%[Characterization of $\mathcal{C}$] \label{Th:Caracterisation des IC chaotiques} -Let $f:\Bool^n\to\Bool^n$. $G_f$ is chaotic (according to Devaney) -if and only if $\Gamma(f)$ is strongly connected. +Let $f:\Bool^n\to\Bool^n$. Iterations of $G_f$ are chaotic according +to Devaney if and only if $\Gamma(f)$ is strongly connected. \end{theorem} Checking if a graph is strongly connected is not difficult. For example, consider the function $f_1\left(x_1,\dots,x_n\right)=\left( -\overline{x_1},x_1,x_2,\dots,x_\mathsf{n}\right)$. As $\Gamma(f_1)$ is +\overline{x_1},x_1,x_2,\dots,x_{n-1}\right)$. As $\Gamma(f_1)$ is obviously strongly connected, then $G_{f_1}$ is a chaotic map. With this material, we are now able to build a first chaotic neural @@ -479,11 +493,11 @@ alternative functions $f$ for $f_0$ through a simple check of their graph of iterations $\Gamma(f)$. For example, we can build another chaotic neural network by using $f_1$ instead of $f_0$. -\section{Checking if a neural network is chaotic or not} +\section{Checking whether a neural network is chaotic or not} \label{S3} We focus now on the case where a neural network is already available, -and for which we want to know if it is chaotic or not. Typically, in +and for which we want to know if it is chaotic. Typically, in many research papers neural network are usually claimed to be chaotic without any convincing mathematical proof. We propose an approach to overcome this drawback for a particular category of multilayer @@ -492,13 +506,13 @@ In spite of this restriction, we think that this approach can be extended to a large variety of neural networks. We plan to study a generalization of this approach in a future work. -We consider a multilayer perceptron of the following form: as inputs -it has $n$ binary digits and one integer value, while it produces $n$ +We consider a multilayer perceptron of the following form: inputs +are $n$ binary digits and one integer value, while outputs are $n$ bits. Moreover, each binary output is connected with a feedback connection to an input one. \begin{itemize} -\item At initialization, the network is feeded with $n$~bits denoted +\item During initialization, the network is seeded with $n$~bits denoted $\left(x^0_1,\dots,x^0_n\right)$ and an integer value $S^0$ that belongs to $\llbracket1;n\rrbracket$. \item At iteration~$t$, the last output vector @@ -524,12 +538,15 @@ $f\left(x_1,x_2,\dots,x_n\right)$ is equal to \left(F\left(1,\left(x_1,x_2,\dots,x_n\right)\right),\dots, F\left(n,\left(x_1,x_2,\dots,x_n\right)\right)\right) \enspace . \end{equation} -Then $F=F_f$ and this recurrent neural network produces exactly the -same output vectors, when feeding it with +Then $F=F_f$. If this recurrent neural network is seeded with $\left(x_1^0,\dots,x_n^0\right)$ and $S \in \llbracket 1;n -\rrbracket^{\mathds{N}}$, than chaotic iterations $F_f$ with initial +\rrbracket^{\mathds{N}}$, it produces exactly the +same output vectors than the +chaotic iterations of $F_f$ with initial condition $\left(S,(x_1^0,\dots, x_n^0)\right) \in \llbracket 1;n -\rrbracket^{\mathds{N}} \times \mathds{B}^n$. In the rest of this +\rrbracket^{\mathds{N}} \times \mathds{B}^n$. +Theoretically speaking, such iterations of $F_f$ are thus a formal model of +these kind of recurrent neural networks. In the rest of this paper, we will call such multilayer perceptrons CI-MLP($f$), which stands for ``Chaotic Iterations based MultiLayer Perceptron''. @@ -558,25 +575,25 @@ that $d\left(f^n(x),f^n(y)\right) \geq \varepsilon$. \begin{definition} \label{def9} A discrete dynamical system is said to be {\bf topologically mixing} if and only if, for any pair of disjoint open sets $U$,$V \neq -\emptyset$, $n_0 \in \mathds{N}$ can be found such that $\forall n -\geq n_0$, $f^n(U) \cap V \neq \emptyset$. +\emptyset$, we can find some $n_0 \in \mathds{N}$ such that for any $n$, +$n\geq n_0$, we have $f^n(U) \cap V \neq \emptyset$. \end{definition} As proven in Ref.~\cite{gfb10:ip}, chaotic iterations are expansive and topologically mixing when $f$ is the vectorial negation $f_0$. Consequently, these properties are inherited by the CI-MLP($f_0$) -recurrent neural network presented previously, which induce a greater +recurrent neural network previously presented, which induce a greater unpredictability. Any difference on the initial value of the input layer is in particular magnified up to be equal to the expansivity constant. -Now, what are the consequences for a neural network to be chaotic -according to Devaney's definition? First of all, the topological +Let us then focus on the consequences for a neural network to be chaotic +according to Devaney's definition. First of all, the topological transitivity property implies indecomposability. \begin{definition} \label{def10} A dynamical system $\left( \mathcal{X}, f\right)$ is {\bf -indecomposable} if it is not the union of two closed sets $A, B +not decomposable} if it is not the union of two closed sets $A, B \subset \mathcal{X}$ such that $f(A) \subset A, f(B) \subset B$. \end{definition} @@ -647,15 +664,32 @@ $\left( \mathcal{X},d\right)$ is compact and the topological entropy of $(\mathcal{X},G_{f_0})$ is infinite. \end{theorem} -We have explained how to construct truly chaotic neural networks, how -to check whether a given MLP is chaotic or not, and how to study its -topological behavior. The last thing to investigate, when comparing -neural networks and Devaney's chaos, is to determine whether -artificial neural networks are able to learn or predict some chaotic -behaviors, as it is defined in the Devaney's formulation (when they +\begin{figure} + \centering + \includegraphics[scale=0.625]{scheme} + \caption{Summary of addressed neural networks and chaos problems} + \label{Fig:scheme} +\end{figure} + +The Figure~\ref{Fig:scheme} is a summary of addressed neural networks and chaos problems. +Section~\ref{S2} has explained how to construct a truly chaotic neural +networks $A$ for instance. +Section~\ref{S3} has shown how to check whether a given MLP +$A$ or $C$ is chaotic or not in the sens of Devaney. +%, and how to study its topological behavior. +The last thing to investigate, when comparing +neural networks and Devaney's chaos, is to determine whether +an artificial neural network $A$ is able to learn or predict some chaotic +behaviors of $B$, as it is defined in the Devaney's formulation (when they are not specifically constructed for this purpose). This statement is studied in the next section. + + + + + + \section{Suitability of Artificial Neural Networks for Predicting Chaotic Behaviors} @@ -694,7 +728,7 @@ Perceptron. We are then left to compute two disjoint function sets that contain either functions with topological chaos properties or not, depending -on the strong connectivity of their iteration graph. This can be +on the strong connectivity of their iterations graph. This can be achieved for instance by removing a set of edges from the iteration graph $\Gamma(f_0)$ of the vectorial negation function~$f_0$. One can deduce whether a function verifies the topological chaos property or @@ -730,7 +764,7 @@ configuration: for instance, 16~(10000) and 15~(01111) are closed in a decimal ordering, but their Hamming distance is 5. This is why Gray codes~\cite{Gray47} have been preferred. -Secondly, how do we deal with strategies. Obviously, it is not +Let us secondly detail how to deal with strategies. Obviously, it is not possible to translate in a finite way an infinite strategy, even if both $(S^t)^{t \in \Nats}$ and $\sigma((S^t)^{t \in \Nats})$ belong to $\{1,\ldots,n\}^{\Nats}$. Input strategies are then reduced to have a @@ -782,10 +816,10 @@ in particular well-known for its universal approximation property. Furthermore, MLPs have been already considered for chaotic time series prediction. For example, in~\cite{dalkiran10} the authors have shown that a feedforward MLP with two hidden layers, and trained -with Bayesian Regulation backpropagation, can learn successfully the +with Bayesian Regulation back-propagation, can learn successfully the dynamics of Chua's circuit. -In these experimentations we consider MLPs having one hidden layer of +In these experiments we consider MLPs having one hidden layer of sigmoidal neurons and output neurons with a linear activation function. They are trained using the Limited-memory Broyden-Fletcher-Goldfarb-Shanno quasi-newton algorithm in combination @@ -884,7 +918,8 @@ obtained results for the non-chaotic case outperform the chaotic ones. Finally, the rates for the strategies show that the different networks are unable to learn them. -For the second coding scheme, Table~\ref{tab2} shows that any network +For the second coding scheme (\textit{i.e.}, with Gray Codes) +Table~\ref{tab2} shows that any network learns about five times more non-chaotic configurations than chaotic ones. As in the previous scheme, the strategies cannot be predicted. @@ -937,7 +972,8 @@ network topologies, the maximum epoch number and the kind of iterations, the configuration success rate is slightly improved. Moreover, the strategies predictions rates reach almost 12\%, whereas in Table~\ref{tab2} they never exceed 1.5\%. Despite of this -improvement, a long term prediction of chaotic iterations appear to be +improvement, a long term prediction of chaotic iterations still +appear to be an open issue. \begin{table}