Stability and Bifurcation Analysis of a Three-Dimensional Recurrent Neural Network with Time Delay
We consider the nonlinear dynamical behavior of a three-dimensional recurrent neural network with time delay. By choosing the time delay as a bifurcation parameter, we prove that Hopf bifurcation occurs when the delay passes through a sequence of critical values. Applying the nor- mal form method and center manifold theory, we obtain some local bifurcation results and derive formulas for determining the bifurcation direction and the stability of the bifurcated periodic solution. Some numerical examples are also presented to verify the theoretical analysis.
Starting with the work of Hopfield  on neural networks, recurrent neural networks including Hopfield neural networks, Cohen-Grossberg neural networks, and cellular neural networks have been used extensively in different areas such as signal processing, pattern recognition, optimization, and associative memories. Many researchers studied the dynamical behavior of Recurrent neural network systems, and most of papers are devoted to the stability of equilibrium, existence and stability of periodic solutions, bifurcation, and chaos [2–5]. In , Ruiz et al. considered a particular configuration of a recurrent neural network, illustrated in Figure 1. In Figure 1, is the input and is the output of the network. This recurrent neural network can be described by the following system: Here, is the state, , are the network parameters or weights, is a smooth input, and is the output. The transfer function of the neurons is taken as . A three-node network of the form (1.1) in feedback configuration, with , has been studied in ; that is, Ruiz et al. found and analyzed the Hopf bifurcation behavior in system (1.2). In , Maleki et al. considered system (1.1) with a transfer function , where for odd and for even. The authors analyzed the Bogdanov-Takens bifurcation in the system.
It is well known that there exist time delays in the information processing of neurons. The delayed axonal signal transmissions in the neural network models make the dynamical behaviors become more complicated and may destabilize the stable equilibria and admit periodic oscillation, bifurcation, and chaos. Therefore, the delay is an important control parameter in living nervous system: different ranges of delays correspond to different patterns of neural activities (see, e.g., [6–11]).
In the present paper, we consider the following three-dimensional recurrent neural network model with time delay By choosing the time delay as a bifurcation parameter, we prove that Hopf bifurcation occurs in the neuron and study the properties of periodic solutions of this model.
The organization of this paper is as follows. In Section 2, by analyzing the characteristic equation of the linearized system of system (1.3) at the equilibrium, we discuss the stability of the equilibrium and the existence of the Hopf bifurcation occurring at the equilibrium. In Section 3, the formulae determining the direction of the Hopf bifurcations and the stability of bifurcating periodic solutions on the center manifold are obtained by using the normal form theory and the center manifold theorem due to Hassard et al. . We do some computer observations to validate our theoretical results in Section 4.
2. Stability and Existence of Hopf Bifurcation
Clearly, is equilibrium of system (1.3). Linearization of (1.3) at the zero equilibrium yields whose characteristic equation is that is, The zero equilibrium is stable if all roots of (2.3) have negative real parts and unstable if at least one root has positive real part. Therefore, in order to study the local stability of the zero equilibrium of system (2.3), we need to investigate the distribution of the roots of (2.3).
Denote From Cardano formula for the third-degree algebra equation, we have the following lemma.
Lemma 2.1. (1) If , then (2.5) has a real root and a pair of conjugate complex roots . Furthermore, the roots of (2.4) are given by and .
(2) If , then (2.5) has a simple root and a multiple root with the multiplicity of 2. Furthermore, the roots of (2.4) are given by and . Meanwhile, if , that is, , then (2.4) has a multiple root −1 with the multiplicity of 3.
(3) If , then (2.5) has three real roots , , and . Furthermore, the roots of (2.4) are given by , , and .
Applying Lemma 2.1, we have the following lemma.
Lemma 2.2. All roots of (2.4) have negative real parts if one of the following holds.
(1) and .
(2) and .
(3) and .
Clearly, is a root of (2.3) if and only if satisfies which implies that This yields Obviously, we have the following lemma.
Summarizing the discussion above, we have the following.
Without loss of generality, one assumes that (2.11) has three positive roots .
Define Then, is a pair of purely imaginary roots of (2.3) with .
Let be the root of (2.3) near satisfying We have the following.
Lemma 2.5. One hasProof. Differentiating both sides of (2.8) with respect to , we have Note that , therefore We have This completes the proof.
For convenience, we let , such that where and is defined in (2.13).
Lemma 2.6. All roots of (2.3) have negative real parts if one of the following holds:, and , and , and .
Lemma 2.7. Suppose that one of the following hypothesis is satisfied:, and .Then, there exists a sequence values of defined by (2.19) such that all roots of (2.3) have negative real parts for all , and (2.3) has at least one root with positive real part when , and (2.3) exactly has a pair of purely imaginary roots when , where and are defined by (2.12) and (2.13), respectively.
Theorem 2.8. (1) If one of the hypothesis , , is satisfied, then the zero solution of system (1.3) is asymptotically stable for all .
(2) If or is satisfied, then the zero solution of system (1.3) is asymptotically stable for and unstable for , and system (1.3) undergoes a Hopf bifurcation at the origin when .
3. Direction of Hopf Bifurcations and Stability of the Bifurcating Periodic Orbits
In this section, we will study the direction of the Hopf bifurcation and stability of bifurcating periodic solutions by using the normal theory and the center manifold theorem due to Hassard et al. .
For , let By the Riesz representation theorem, there exists a function of bounded variation for , such that In fact, we can choose where denotes the Dirac delta function. For , define where Then system (3.1) is equivalent to where , for .
For , define and a bilinear inner product where . Then, and are adjoint operators. By the discussion in Section 2, we know that are eigenvalues of and corresponding to and , respectively.
In order to assure , we need to determine the value of . By (3.10), we haveTherefore, we can choose as
Following the algorithms given in  and using similar computation process in , we can get that the coefficients which will be used to determine the important quantities: where moreover , satisfy the following equations, respectively, Therefore, all in (3.16) can be expressed in terms of parameters. And we can compute the following values: which determine the qualities of bifurcating periodic solution in the center manifold at the critical values , that is, determines the directions of the Hopf bifurcation: if , then the Hopf bifurcation is supercritical (subcritical) and the bifurcating periodic solutions exist for ; determines the stability of the bifurcating periodic solutions: the bifurcating periodic solutions are stable (unstable) if ; determines the period of the bifurcating periodic solutions: the period increases (decreases) if .
4. Computer Simulation
In this section, we will confirm our theoretical analysis by numerical simulation. We give an example of system (3.1) with , and . Then, and .
Equation (1.3) becomes From (2.6), we have . By Lemma 2.1, we know (2.7) has roots and . Clearly, , . From (2.12), it follows that , and, from (2.13), we get . Thus, the zero equilibrium is asymptotically stable when as is illustrated by the computer simulations (see Figures 2(a)–2(d)). When passes through the critical value , zero equilibrium loses its stability and a Hopf bifurcation occurs, that is, a family of periodic solutions bifurcates from the origins , which are depicted in Figures 3(a)–3(d).
The author is very grateful to the referees for their valuable suggestions, which help to improve the paper significantly.
J. J. Hopfield, “Neurons with graded response have collective computational properties like those of two-state neurons,” Proceedings of the National Academy of Sciences of the United States of America, vol. 81, no. 10, pp. 3088–3092, 1984.View at: Google Scholar
A. C. Ruiz, D. H. Owens, and S. Townley, “Existence, learning, and replication of periodic motions in recurrent neural networks,” IEEE Transactions on Neural Networks, vol. 9, no. 4, pp. 651–661, 1998.View at: Google Scholar
B. Hassard, D. Kazarinoff, and Y. Wan, Theory and Applications of Hopf Bifurcation, Cambridge University Press, Cambridge, Mass, USA, 1981.
S. Ruan and J. Wei, “On the zeros of transcendental functions with applications to stability of delay differential equations with two delays,” Dynamics of Continuous, Discrete & Impulsive Systems. Series A. Mathematical Analysis, vol. 10, no. 6, pp. 863–874, 2003.View at: Google Scholar | Zentralblatt MATH
J. K. Hale, Theory of Functional Differential Equations, Springer, New York, NY, USA, 1977.