Mathematical Problems in Engineering

Volume 2015 (2015), Article ID 614135, 6 pages

http://dx.doi.org/10.1155/2015/614135

## Controlling Hopf Bifurcation of a New Modified Hyperchaotic Lü System

^{1}College of Mechanical and Vehicle Engineering, Hunan University, Changsha 410082, China^{2}School of Mathematics and Statistics, Minnan Normal University, Zhangzhou 363000, China

Received 23 October 2014; Revised 3 March 2015; Accepted 11 March 2015

Academic Editor: Michael Vynnycky

Copyright © 2015 Ping Cai 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

Controlling Hopf bifurcation of a new modified hyperchaotic Lü system is investigated in this paper. A hybrid control strategy using both state feedback and parameter control is proposed. The control strategy realizes the delay of Hopf bifurcation. Furthermore, by applying the normal form theory, the stability of the bifurcation is determined. Numerical simulation results are given to support the theoretical analysis.

#### 1. Introduction

Chaos and bifurcation are of great importance in many physical, chemical, and biological nonlinear systems. Chaos theory is one of the most significant achievements of nonlinear science. Nowadays complex nonlinear systems are being used in many fields of science and engineering. Most of the recent works have focused on solving complex chaos control, synchronization, and so on [1–4]. Bifurcation control has been a rapidly growing interest by many research works in recent years [5–7]. Hopf bifurcation is a kind of important dynamic bifurcation. Contributions in Hopf bifurcation control mainly focused on amplitude control of limit cycle [8, 9], changing the critical points of an existing bifurcation [10], delaying the onset of an inherent bifurcation, and stabilizing an existing bifurcation [11–13], creating a desired bifurcation at a desired location, which is called anticontrol of bifurcation [14–17]. Among these researches, 3D chaotic systems play a leading role, such as Lorenz system [18], Chen system [19], and Lü system [20]. The study of 4D hyperchaotic systems has recently become a hot topic [21–24]. 4D hyperchaotic system has more complicated dynamical behaviors compared to 3D chaotic system. So the analysis and calculation are also more difficult.

The idea of the work is to design a control law to control Hopf bifurcation in nonlinear system. Consider the following general nonlinear system:where the dot denotes differentiation with respect to time and is an -dimensional state vector, while is bifurcation parameter. Let be equilibria of system (1); that is, for any value of . A hybrid control strategy is added to model (1), and then we obtain the following controlled system: where is a control parameter and is state feedback controller. In order for the controlled system (2) to keep all the equilibria unchanged under the control, the following conditions should be satisfied:A general formula satisfying condition (3) can be constructed as follows:Obviously, we use nonlinear feedback with polynomial function in the controller . Generally speaking, linear part of a control strategy is used to shift the bifurcation value, in order to eliminate or delay an existing bifurcation. The nonlinear part, on the other hand, can be designed to change the stability of bifurcation solutions. Controller involves higher-order terms, which may not be necessary for stability control. It is preferable to have the simplest possible design for engineering applications. In most cases, using fewer components or just one component may be enough to satisfy the predesigned control objectives. So, it is not necessary to take all the components in the controller for practical system. This greatly simplifies the control formula. For example, if system (1) has two equilibria , , then the general controller can be taken as the following simple form:If system (1) has only one equilibrium , the controller can be taken as We also omit the linear term in (6), because parameter has the same control effect. In fact, it even can be more simple asFor Hopf bifurcation control, as a result of the calculation formula of stability index [25], terms up to the third-order term are enough and the second-order term might not be necessary due to the presence of the third-order term for the simplicity of the calculation. The following are the conditions of system (2) undergoing Hopf bifurcation at the equilibrium.

Let be the Jacobian matrix of system (2) evaluated at . By the Hopf theory [25], contains a complex conjugate pair of eigenvalues satisfyingand the remaining eigenvalues of have negative real parts at the critical point . That is to say, when is varied, the pair of the complex conjugates moves to cross the imaginary axis at . The second condition of (8) is usually called the transversality condition, implying that the crossing of the complex conjugate pair at the imaginary axis is not tangent to the imaginary axis. Without loss of generality, assume that when is varied from to , moves from the left-half of complex plane to the right. Thus, a family of limit cycles will bifurcate from the equilibrium at the critical point .

Next, it will be shown that the parameter can change the bifurcation critical value, and the nonlinear state feedback can ensure the stability of bifurcation solutions.

#### 2. Hybrid Control of Hopf Bifurcation

In this paper, the hybrid control is applied to Hopf bifurcation control of a new modified 4D hyperchaotic Lü system of the form [24]where are state variables and are real parameters. System (9) has a hyperchaotic attractor when , , , and [24]. Obviously, system (9) has only one isolated equilibrium when . The controlled system is is also the equilibrium of system (10). The Jacobian matrix of system (10) at is The characteristic equation of is Taking as the Hopf bifurcation parameter and supposing that (12) has a pair of pure imaginary roots , which leads to The other two roots are Thus, the necessary conditions for system (10) to exhibit Hopf bifurcation at are , , , and . Under these conditions, the transversality condition is also satisfied. Therefore, system (10) undergoes Hopf bifurcation at the equilibrium based on Hopf bifurcation theory [25].

*Remark 1. *In particular, if , system (10) is reverted to the original system (9). The Hopf bifurcation value of the original system is . So, parameter can change the Hopf bifurcation value.

*Remark 2. *By formula (13), we notice that does not affect the bifurcation critical value, so we can set in system (10). And for simplicity, we also set in the following discussion; that is to say, we may only choose the first two equations of system (9) under control. In this case, , are obtained.

#### 3. Analysis of Stability of Hopf Bifurcation

In this section, we apply the normal form theory [25] to study the stability of the Hopf bifurcation for system (10).

By the linear transform , wherethen system (10) has the following normal form:whereThen the stability condition of the bifurcated limit circle can be derived [25]:If , the bifurcated periodic solution is orbitally asymptotically stable, and if , it is unstable. The following three special cases are considered.

*Case 1. *If , then . For , , we choose ; then the bifurcation critical value of system (10) satisfies . Moreover, we have .

*Case 2. *If , that is to say, only the first equation of system (9) is under control, then . For , , we choose ; then the bifurcation critical value of system (10) satisfies . Moreover, we have .

*Case 3. *If , that is to say, only the second equation of system (9) is under control, then . For , , we choose ; then the bifurcation critical value of system (10) satisfies . Moreover, we have .

So, theoretical analyses show that the control strategy not only delays Hopf bifurcation but also achieves the stability control of the bifurcation.

#### 4. Numerical Simulations

In this section, numerical simulations are given to illustrate the above theoretical analyses. We choose , , and , and the original system (9) undergoes Hopf bifurcation at . The bifurcation figure of the original system (9) is shown in Figure 1. If we set , the Hopf bifurcation critical value of the controlled system (10) is , and . The bifurcation figure of the controlled system (10) is shown in Figure 2. Time displacement curves and phase space trajectories are shown in Figures 3 and 4, respectively.