Journal of Applied MathematicsVolume 2012 (2012), Article ID 357382, 13 pagesdoi:10.1155/2012/357382

Research Article

Stability and Bifurcation Analysis of a Three-Dimensional Recurrent Neural Network with Time Delay

Yingguo Li

School of Mathematics and Computer Science, Fujian Normal University, Fuzhou 350007, China

Received 14 August 2011; Revised 4 October 2011; Accepted 5 October 2011

Academic Editor: Shiping Lu

Copyright © 2012 Yingguo Li. 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.

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.

1. Introduction

Starting with the work of Hopfield [1] 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 [5], 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 [5]; that is, Ruiz et al. found and analyzed the Hopf bifurcation behavior in system (1.2). In [4], 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.

Figure 1: Class of recurrent neural networks.

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. [12]. We do some computer observations to validate our theoretical results in Section 4.

2. Stability and Existence of Hopf Bifurcation

For most of the models in the literature, including the ones in [5, 7, 8], the activation function is . However, we only make the following assumption on function :(H).

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).

When , characteristic equation (2.3) yields Let , (2.4) reduces to

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 .

Let , then (2.3) becomes We notice that (2.7) and (2.5) have the same coefficients.

Denote the three roots of (2.7) by . Hence, (2.3) is equivalent to

Clearly, is a root of (2.3) if and only if satisfies which implies that This yields Obviously, we have the following lemma.

Lemma 2.3. Equation (2.11) is meaningless when . If (2.7) has a root, denoted by , satisfying , then (2.11) has a positive root given by

Summarizing the discussion above, we have the following.

Lemma 2.4. If , then (2.11) has at most two positive roots. If , then (2.11) has at most three positive roots.

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 has Proof. 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).

Applying Lemmas 2.1–2.5 and Corollary 2.4 of Ruan and Wei [13], we have the following results.

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.

From Lemmas 2.5–2.7 and the Hopf bifurcation theorem for functional differential equations in [14], we have the theorem.

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. [12].

Let , , , then system (1.3) becomes functional differential equation in as whereSetting , we know that is Hopf bifurcation value of system (3.1).

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.

Suppose is the eigenvectors of corresponding to , then . Then from the definition of and (3.3)–(3.5), we have For , then we obtain Similarly, we can obtain the eigenvector of corresponding to , where

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 [12] and using similar computation process in [7], 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).

(a)

(b)

(c)

(d)

Figure 2: When = 1.7 < = 1.8434, the zero equilibrium is asymptotically stable. Here initial value is (1, 1, 1).

(a)

(b)

(c)

(d)

Figure 3: When = 1.9 > = 1.8434, a periodic orbit bifurcates from the zero equilibrium. Here, initial value is (1, 1, 1).

Acknowledgment

The author is very grateful to the referees for their valuable suggestions, which help to improve the paper significantly.

References

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.

B. Gao and W. Zhang, “Equilibria and their bifurcations in a recurrent neural network involving iterates of a transcendental function,” IEEE Transactions on Neural Networks, vol. 19, no. 5, pp. 782–794, 2008.

R. Haschke and J. J. Steil, “Input space bifurcation manifolds of recurrent neural networks,” Neurocomputing, vol. 64, no. 1–4, pp. 25–38, 2005.

F. Maleki, B. Beheshti, A. Hajihosseini, and G. R.R. Lamooki, “The Bogdanov-Takens bifurcation analysis on a three dimensional recurrent neural network,” Neurocomputing, vol. 73, no. 16–18, pp. 3066–3078, 2010.

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.

X. Liao, K. W. Wong, C. S. Leung, and Z. Wu, “Hopf bifurcation and chaos in a single delayed neuron equation with non-monotonic activation function,” Chaos, Solitons and Fractals, vol. 12, no. 8, pp. 1535–1547, 2001.

J. Wei and S. Ruan, “Stability and bifurcation in a neural network model with two delays,” Physica D. Nonlinear Phenomena, vol. 130, no. 3-4, pp. 255–272, 1999.

K. Gopalsamy and I. Leung, “Delay induced periodicity in a neural netlet of excitation and inhibition,” Physica D. Nonlinear Phenomena, vol. 89, no. 3-4, pp. 395–426, 1996.

C. Huang, Y. He, L. Huang, and Y. Zhaohui, “Hopf bifurcation analysis of two neurons with three delays,” Nonlinear Analysis. Real World Applications, vol. 8, no. 3, pp. 903–921, 2007.

X.-P. Yan, “Hopf bifurcation and stability for a delayed tri-neuron network model,” Journal of Computational and Applied Mathematics, vol. 196, no. 2, pp. 579–595, 2006.

D. Fan and J. Wei, “Hopf bifurcation analysis in a tri-neuron network with time delay,” Nonlinear Analysis. Real World Applications, vol. 9, no. 1, pp. 9–25, 2008.

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.

J. K. Hale, Theory of Functional Differential Equations, Springer, New York, NY, USA, 1977.