We introduce and investigate a four-dimensional hidden hyperchaotic system without equilibria, which is obtained by augmenting the three-dimensional self-exciting homopolar disc dynamo due to Moffatt with an additional control variable. Synchronization of two such coupled disc dynamo models is investigated by active control and sliding mode control methods. Numerical integrations show that sliding mode control provides a better synchronization in time but causes chattering. The solution is obtained by switching to active control when the synchronization errors become very small. In addition, the electronic circuit of the four-dimensional hyperchaotic system has been realized in ORCAD-PSpice and on the oscilloscope by amplitude values, verifying the results from the numerical experiments.
National Natural Science Foundation of China11561069Guangxi Natural Science Foundation of China2016GXNSFBA380170Scientific Research Foundation of the Higher Education Institutions of Guangxi Province of ChinaKY2016YB364Open Foundation for Guangxi Colleges and Universities Key Lab of Complex System Optimization and Big Data Processing2016CSOBDP00032016CSOBDP0202Sakarya University Scientific Research Projects Unit201609-00-0081. Introduction
Hyperchaos is a feature of a chaotic system having more than one positive Lyapunov exponent [1]. Because of potential theoretical and practical applications in technology, such as secure communications, lasers, nonlinear circuits, neural networks, generation, control, and synchronization, hyperchaos has featured as an important research area in nonlinear science [2–5]. The theory about hidden hyperchaos with either only stable or no equilibrium states is still in its infancy and has only recently been understood by mathematicians [6–9].
In 1979, Moffatt identified inconsistencies in the modeling of a simple self-exciting homopolar disc dynamo because of the neglect of induced azimuthal eddy currents, which can be resolved by introducing a segmented disc dynamo [10]. Here we investigate hidden hyperchaos, synchronization, and electronic circuit realization for a higher-dimensional version of the self-exciting homopolar disc dynamo, which was not yet completely well understood.
Since Pecora and Carroll [11] investigated synchronization in chaotic systems in 1990, such behavior has become an important research area in nonlinear science, not only for understanding the complicated phenomena in various fields but also for its potential applications especially in secure communication and image encryption. Two indistinguishable chaotic systems, starting from nonidentical initial values, would evolve in time to completely different trajectories because of the sensitive dependence of chaotic systems to their initial values. The aim of synchronizing chaos is to ensure that the states track the desired trajectory. Many effective methods exist to deal with synchronization of chaotic and hyperchaotic systems. These include active control [12–20], passive control [21], sliding mode control [22–31], adaptive control [32], and backstepping design [33]. Of these, active control is an important simple method used in the synchronization of nonlinear systems. It maintains asymptotic stability at zero error by eliminating the nonlinear terms and making all the eigenvalues have negative real parts. The other commonly preferred method, sliding mode control, maintains the synchronization by enforcing the error system to stay on a constructed sliding surface.
The first model with a simple electronic application was realized by Chua et al. [34]. In the following years, many electronic circuit applications such as simple RLC, RC circuits [35–37], oscillators [38, 39], power circuits [40, 41], and capacitor circuits which show chaotic features were improved upon. On the one hand, numerous electronic circuit realizations with interesting features, which mimic novel chaotic and hyperchaotic systems, have been proposed in recent years [3, 42–45].
Current interest in hidden hyperchaotic attractors motivates us to study an extension about the self-exciting homopolar disc dynamo [10] to 4D homopolar dynamo without equilibria. The existence of hidden hyperchaotic attractors in this new disc dynamo is confirmed. Synchronization of two such coupled 4D self-exciting homopolar disc dynamo systems is analyzed with active and sliding mode control methods. Moreover, we have designed an electronic circuit and have used an oscilloscope to view the hyperchaotic rescaled dynamo without equilibria, implemented in real time.
2. Model and Hidden Hyperchaos of 4D Self-Exciting Homopolar Disc Dynamo System without Equilibria
Dynamo models have been the object of much interest in order to understand both the generation of magnetic fields and their reversals in astrophysics. Moffatt [10] extended the simplest self-exciting Bullard dynamo to include radial magnetic diffusion, to produce the disc dynamo model, written nondimensionally as(1)x˙=ry-x,y˙=mx-1+my+xz,z˙=g1+mx2-1+mxy,where x, y, and z are state variables and r, m, and g are positive parameters. x and y represent the magnetic fluxes due to radial and azimuthal current distributions, respectively. z denotes the angular velocity of the disc. g is the applied torque; r and m are the constants which depend on the inductances and the electrical resistance of the dynamo.
By modifying the characteristics of the segmented disc dynamo (1), hidden chaotic or hyperchaotic spiral attractors have been observed numerically under special initial conditions with two symmetric stable node-foci. This leads to the interesting and striking observation of multiple attractors for a broad range of parameters. As in many nonlinear dynamical systems, the occurrence of multiple attractors implies the existence of multistability in the self-exciting dynamo, with the long-term behavior being fundamentally different depending on which basin of attraction the initial conditions belong. Now, we introduce a dislocated feedback controller to system (1) as a new state w and translate z to z-m to result in following 4D system:(2)x˙=ry-x+w,y˙=-1+my+xz,z˙=g1+mx2-1+mxy,w˙=-ky,where k is a positive parameter. Although system (2) is similar to the algebraic forms of hyperchaotic Lorenz, Chen, Lü, and unified systems, they are not topologically equivalent [46–52]. The proposed system (2) will be a way of understanding the generation of magnetic fields and their reversals in the Earth, the Sun, and other astrophysical bodies. Figure 1 shows a bifurcation diagram exhibiting a period-doubling route to chaos of the peak of z (zmax) of system (2) versus the parameters m∈[0.5,1], r=8, g=35, and k=3. There are some periodic windows in the chaotic region. Plots of the Lyapunov exponents about m∈[0.5,1.5] are shown in Figure 2. Figure 3 indicates that system (2) is indeed hidden hyperchaotic for initial states (1.13, 0.5, 0.8, and 1.5) and parameter m=0.5. Its Lyapunov exponents are 0.4113, 0.2233, 0.0000, and −10.1345) and Kaplan-Yorke dimension is DKY = 3.0626.
Bifurcation diagram of system (2) without equilibria versus parameter m in [0.5, 1.5].
Lyapunov exponents LE1, LE2, LE3, and LE4 of system (2) without equilibria versus parameter m in [0.5, 1.5].
Hyperchaotic attractor of four-dimensional self-exciting homopolar disc dynamo system (2) without equilibria for r=8, g=35, k=3, and m=0.5.
To find the equilibrium states of system (2), we set x˙=y˙=z˙=w˙=0 and solve(3)ry-x+w=0,-1+my+xz=0,g1+mx2-1+mxy=0,-ky=0.
There seems to be two equilibria: E1(1/-m,0,0,r/-m) and E2(-1/-m,0,0,-r/-m); however, because m is positive, hyperchaotic self-exciting homopolar disc dynamo system (2) has no real equilibria.
3. Synchronization of the 4D Self-Exciting Homopolar Disc Dynamo System
For synchronization, two self-exciting homopolar disc dynamo hyperchaotic systems are coupled together with different initial values. The driver system, x, controls the response system, y. They are given, respectively, by(4)x1˙=rx2-x1+x4,x2˙=-1+mx2+x1x3,x3˙=g1+mx12-1+mx1x2,x4˙=-kx2,(5)y1˙=ry2-y1+y4,y2˙=-1+my2+y1y3+u1,y3˙=g1+my12-1+my1y2+u2,y4˙=-ky2,where u1 and u2 in system (5) are the control functions to be determined. The synchronization errors are obtained by subtracting the driver and response systems. Thus, they are defined as ei=yi-xi (for i=1,2,3,4) and the error dynamics become as follows:(6)e1˙=re2-e1+e4,e2˙=-1+me2+y1y3-x1x3+u1,e3˙=gmy12-x12-1+my1y2-x1x2+u2,e4˙=-ke2.Our objective is to make system (6) asymptotically stable at the zero error state.
3.1. Active Control
The nonlinear terms in system (6) can be eliminated by defining the control functions u1 and u2 as in the following:(7)u1=-y1y3+x1x3+v1,u2=-gmy12-x12-1+my1y2-x1x2+v2.Then, substituting (7) into system (6) gives(8)e1˙=re2-e1+e4,e2˙=-1+me2+v1,e3˙=v2,e4˙=-ke2.This implies that the equations in system (8) are linear. Provided that the proper choices of control inputs v1 and v2 stabilize the error system (8), then e1, e2, e3, and e4 will converge to zero as t→+∞. Then, the synchronization of two identical self-exciting homopolar disc dynamo hyperchaotic systems will be achieved. A number of choices are possible for control functions v1 and v2. They are simply taken as(9)v1=1+m-k1e2-re1+ke4,v2=-k2e3,where k1 and k2 are positive control gains. Substituting (9) to system (8) gives the following synchronization error dynamics:(10)e1˙=re2-e1+e4,e2˙=-re1-k1e2+ke4,e3˙=-k2e3,e4˙=-ke2.The associated characteristic matrix of system (10) is(11)A=-r-r01-r-k10k00-k200-k00. For the particular choice of control functions in (9), the closed loop system (10) has all of its eigenspectrum in the negative half plane since all of the parameters r, k, k1, and k2 are positive. So, this choice leads to a stable system where the error states e1, e2, e3, and e4 tend to zero as time t tends to infinity. Synchronization of two identical hyperchaotic self-exciting homopolar disc dynamos is therefore completed with the active control method.
3.2. Sliding Mode Control
It can be seen from system (6) that when e2 and e3 become zero, e4˙ will be zero and then e1˙=-re1. Therefore, when time goes to infinite, e1 will converge to zero, too. Appropriate sliding surfaces for e2 and e3 can be, respectively, designed as(12)s1=e2-k3e4,s2=e3-k4e4,where k3 and k4 are positive control gains.
The attainability condition for sliding mode is ss˙<0. Provided that this condition is satisfied, the sliding mode control functions are(13)u1=1+me2-y1y3+x1x3-k3ke2-k5s1-q1signs1,u2=-gmy12-x12-1+my1y2-x1x2-k4ke2-k6s2-q2signs2,where k5 and k6 are positive control gains. Large values of k5 and k6 decrease the time t taken to reach the sliding surface but lead to chattering; small values of q1 and q2 reduce chattering but increase the time to reach the sliding surface. Here “sign” means the signum function.
The designed control functions (13) guarantee that system (6) is held on the sliding surface s=0. The time derivations of sliding surfaces are(14)s1˙=-k5s1-q1signs1,s2˙=-k6s2-q2signs2.For the Lyapunov function(15)V=12s12+s22,the time derivative of V becomes(16)V˙=-k5s12-q1s1signs1-k6s22-q2s2signs2≤0,where k5,6≥0. These conditions guarantee that the constructed sliding surfaces s1 and s2 would asymptotically stabilize to the zero synchronization error state, and we obtain synchronization between the two identical hyperchaotic self-exciting homopolar disc dynamos via the sliding mode control method.
3.3. Numerical Simulations
We now perform some numerical experiments to show that synchronization occurs. We use an ode45 integration solver function with a variable step size. We take r=8, m=0.5, g=35, and k=3 to ensure that hyperchaotic behavior occurs. The gains of the active controllers are taken to be k1=1 and k2=1. The gains of the sliding mode controllers are taken to be k3=1, k4=1, k5=1, k6=1, q1=0.5, and q2=0.5. We choose the initial conditions for the driver and response systems to be (1.13,0.5,0.8,1.5) and (1.3,2,11.1,1.4), respectively. The controllers are activated when t=10, and we plot the synchronization simulations in Figure 4, while the synchronization errors are plotted in Figure 5.
Time series of driver and response four-dimensional self-exciting homopolar disc dynamo hyperchaotic systems with the controllers are activated at t=10.
Synchronization errors between driver and response four-dimensional self-exciting homopolar disc dynamo hyperchaotic systems with the controllers are activated at t=10: (a) active controllers and (b) sliding mode controllers.
Figure 4 shows that both active and sliding mode controllers achieve the synchronization of the four-dimensional self-exciting homopolar disc dynamo hyperchaotic system. Figure 5 also shows that, after the activation of controllers at t=10, the synchronization errors approach zero asymptotically, thereby validating the theoretical analyses. Synchronization is complete for t≥16 with active control and for t≥14 with sliding mode control. Activation of control at different times gives similar synchronization performances. The comparisons point out that the sliding mode control scheme has the advantage of a faster synchronization time. However, when the results are viewed more closely, as in Figure 6, the chattering problem of sliding mode controllers is evident. A solution is to switch to active control when the mean squared error is less than 0.00001. The new results from switching controllers are presented in Figure 7. It has no chattering and the synchronization performance is similar to that with the sliding mode control scheme.
Synchronization errors between t≥12 and t≤20: (a) active controllers and (b) sliding mode controllers.
Synchronization errors with the proposed switching controllers are activated: (a) t≤20; (b) t≥12, t≤20.
4. Electronic Circuit Implementation of the 4D Self-Exciting Homopolar Disc Dynamo System
Because the values for x, y, and w fall within the interval of (-15,15), while z falls within the interval of (-20,25), z must be rescaled for observations using the oscilloscope. We let X=x, Y=y, Z=z/2, and W=w and consider the system:(17)X˙=rY-X+W,Y˙=-1+mY+2XZ,Z˙=12g1+mX2-1+mXY,W˙=-kY.
We can now design an electronic circuit for the scaled hyperchaotic model (17) with parameters r=8, m=0.5, g=35, and k=3 and initial conditions (1.13,0.5,0.4,1.5) using ORCAD-PSpice. The schematic of the electronic circuit is shown in Figure 8. We used resistors, capacitors, TL081 opamps, and AD633 multipliers with the values R1=R2=50 Kohm, R3=400 Kohm, R4=266.66 Kohm, R5=20 Kohm, R6=133.33 Kohm, R7=342 Kohm, R8=4.6 Kohm, R9=1.5 Kohm, R10=R11=R12=R13=100 Kohm, C1=C2=C3=C4=1 nF, Vn=-15 V, and Vp=15 V. Since the AD633 multiplier IC is limited to between −10 V and +10 V, the output voltage must be divided by 10 V. Real-time application of system (17) was realized with electronic components and shown on the electronic card in Figure 9.
The electronic circuit schematic of the scaled hyperchaotic system (17).
The experimental circuit of the hyperchaotic circuit.
Outputs from the ORCAD-PSpice simulation and oscilloscope phase portraits for the scaled hyperchaotic system (17) with parameters r=8, m=0.5, g=35, and k=3 are given in Figures 10 and 11, respectively. The outputs verify those of the hyperchaotic system, which was modeled using MATLAB.
The phase portraits of scaled hyperchaotic system (17) in ORCAD-PSpice.
The phase portraits of scaled hyperchaotic system (17) with r=8, m=0.5, g=35, and k=3 on the oscilloscope.
5. Conclusion
In this paper, we propose a novel four-dimensional self-exciting homopolar disc dynamo without equilibria, but exhibiting hidden hyperchaos. Furthermore, active control and sliding mode control methods are applied to synchronize two coupled four-dimensional dynamo systems. The feasibility of active controllers is ensured by requiring that the spectrum of eigenvalues of the synchronized error system falls in the left half plane. A Lyapunov function is proposed to guarantee the asymptotic stability of sliding surfaces and convergence to zero synchronization error. Numerical integrations are presented to compare the performances of the two controllers. Sliding mode control scheme gives better results but has chattering. Solution is provided by a switch to active controllers when chattering starts. An electronic circuit for the rescaled system is implemented via the ORCAD-PSpice program. Numerical simulations validated the theoretical analyses. This physical example will form the basis of more systematic studies of hyperchaos without equilibria in a future study.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported by the National Natural Science Foundation of China (Grant no. 11561069), the Guangxi Natural Science Foundation of China (Grant no. 2016GXNSFBA380170), the Scientific Research Foundation of the Higher Education Institutions of Guangxi Province of China (Grant no. KY2016YB364), the Open Foundation for Guangxi Colleges and Universities Key Lab of Complex System Optimization and Big Data Processing (nos. 2016CSOBDP0003 and 2016CSOBDP0202), and Sakarya University Scientific Research Projects Unit (no. 201609-00-008).
RösslerO. E.An equation for hyperchaos1979712-315515710.1016/0375-9601(79)90150-6MR588951ArenaP.BaglioS.FortunaL.ManganaroG.Hyperchaos from cellular neural networks199531425025110.1049/el:199501892-s2.0-0029254491LiQ.TangS.ZengH.ZhouT.On hyperchaos in a small memristive neural network2014782108710992-s2.0-8491100535110.1007/s11071-014-1498-7Zbl1331.94090LiuL.LiuC.ZhangY.Theoretical analysis and circuit implementation of a novel complicated hyperchaotic system201166470771510.1007/s11071-011-9943-3MR2859596LiC. B.PehlivanI.SprottJ. C.AkgulA.A novel four-wing strange attractor born in bistability201512411210.1587/elex.12.201411162-s2.0-84923350748KiselevaM. A.KuznetsovN. V.LeonovG. A.Hidden attractors in electromechanical systems with and without equilibria20164914515510.1016/j.ifacol.2016.07.975LeonovG. A.KuznetsovN. V.Hidden attractors in dynamical systems: From hidden oscillations in Hilbert-Kolmogorov, Aizerman and Kalman problems to hidden chaotic attractor in Chua circuits2013231133000210.1142/S0218127413300024MR3038624WeiZ.YuP.ZhangW.YaoM.Study of hidden attractors, multiple limit cycles from Hopf bifurcation and boundedness of motion in the generalized hyperchaotic Rabinovich system2015821-213114110.1007/s11071-015-2144-8MR3405634WeiZ. C.ZhangW.Hidden hyperchaotic attractors in a modified Lorenz-Stenflo system with only one stable equilibrium20142410145012710.1142/S0218127414501272MR3273806MoffattH. K.A self-consistent treatment of simple dynamo systems197914114716610.1080/030919279082445362-s2.0-0018725484PecoraL. M.CarrollT. L.Synchronization in chaotic systems199064882182410.1103/PhysRevLett.64.821MR10382632-s2.0-0343689904AgizaH. N.YassenM. T.Synchronization of Rossler and Chen chaotic dynamical systems using active control2001278419119710.1016/S0375-9601(00)00777-5MR1827067BaiE.-W.LonngrenK. E.Synchronization of two Lorenz systems using active control199781515810.1016/s0960-0779(96)00060-4Zbl1079.375152-s2.0-0002959919El-DessokyM. M.Synchronization and anti-synchronization of a hyperchaotic Chen system2009394179017972-s2.0-6344913040910.1016/j.chaos.2007.06.053GoksuA.KocamazU. E.UyarogluY.Synchronization and control of chaos in supply chain management20158610711510.1016/j.cie.2014.09.0252-s2.0-84952794163LiG.-H.An active control synchronization for two modified Chua circuits20051434724752-s2.0-2414446777910.1088/1009-1963/14/3/007UçarA.LonngrenK. E.BaiE.-W.Synchronization of the unified chaotic systems via active control2006275129212972-s2.0-2494450161210.1016/j.chaos.2005.04.104VincentU. E.Synchronization of Rikitake chaotic attractor using active control20053431–313313810.1016/j.physleta.2005.06.003Zbl1194.340912-s2.0-21744449163WeiY. Q.JiangN.Synchronization of hyperchaotic Rossler system and hyperchaotic Lorenz system with different structure6988Int. Conf. on Web Information Systems and Mining2011Taiyuan, China147153ZhouX. B.WuY.LiY.XueH. Q.Synchronization of Hyperchaotic Systems via Active ControlProceedings of 5th International Conference on Communications, Circuits and SystemsJuly 2007Kokura1094109810.1109/ICCCAS.2007.4348237ZhouX.KongB.DingH.Synchronization and anti-synchronization of a new hyperchaotic Lü system with uncertain parameters via the passive control technique20128562-s2.0-8486114139210.1088/0031-8949/85/06/065004065004Zbl1282.34056ChenC.-K.YanJ.-J.LiaoT.-L.Sliding mode control for synchronization of Rössler systems with time delays and its application to secure communication20077654364412-s2.0-3624895374210.1088/0031-8949/76/5/005Zbl1125.37311HouY.-Y.LiauB.-Y.ChenH.-C.Synchronization of unified chaotic systems using sliding mode controller201220122-s2.0-8487215233910.1155/2012/632712632712KocamazU. E.GöksuA.TaşkınH.UyaroğluY.Synchronization of chaos in nonlinear finance system by means of sliding mode and passive control methods: a comparative study201544217218110.5755/j01.itc.44.2.77322-s2.0-84936745157AlmeidaD. I. R.AlvarezJ.BarajasJ. G.Robust synchronization of Sprott circuits using sliding mode control200630111182-s2.0-3364615091010.1016/j.chaos.2005.09.011ToopchiY.WangJ.Chaos control and synchronization of a hyperchaotic Zhou system by integral sliding mode control201416126539655210.3390/e16126539MR3299545VaidyanathanS.Sliding mode controller design for the global chaos synchronization of rucklidge chaotic systems201221233110.5121/ijbb.2012.2103VaidyanathanS.Sliding mode controller design for the global chaos synchronization of hyperchaotic Lü systems201223213110.5121/ijcsa.2012.2303VaidyanathanS.SampathS.Global chaos synchronization of hyperchaotic Lorenz systems by sliding mode control205Proceedings of 1st Int. Conf. on Digital Image Processing and Pattern Recognition2011Tirunelveli, India156164VaidyanathanS.SivaperumalS.Global chaos synchronization of the hyperchaotic Qi systems by sliding mode control20113624302437YanJ.-J.LinJ.-S.LiaoT.-L.Synchronization of a modified Chua's circuit system via adaptive sliding mode control2008361455210.1016/j.chaos.2006.06.008MR2372047VaidyanathanS.Adaptive control and synchronization of hyperchaotic Newton-Leipnik system2011132233HuY. A.SunM. M.Synchronization of a class of hyperchaotic systems via backstepping and its application to secure communicationProceedings of 5th Int. Conf. on Instrumentation & Measurement, Computer, Communication, and Control2015Qinhuangdao, ChinaComputer, Communication, and Control10551060ChuaL. O.WuC. W.HuangA.ZhongG.-Q.Universal circuit for studying and generating chaos—part I: routes to chaos1993401073274410.1109/81.2461492-s2.0-0027681883NakagawaS.SaitoT.An RC OTA hysteresis chaos generator199643121019102110.1109/81.5458462-s2.0-0030382096OgorzalekM. J.Order and chaos in a third-order RC ladder network with nonlinear feedback19893691221123010.1109/31.34668MR1009275TamaševičiusA.NamajunasA.ČenysA.Simple 4D chaotic oscillator199632119579582-s2.0-003014198910.1049/el:19960630KawakamiH.Bifurcation of periodic responses in forced dynamic nonlinear circuits: computation of bifurcation values of the system parameters198431324826010.1109/TCS.1984.1085495MR739206SaitoT.Chaotic phenomena in a coupled oscillatorProceedings of European Conf. on Circuit Theory and Design.1987275280MR1052171HamillD. C.JeffriesD. J.Subharmonics and chaos in a controlled switched-mode power converter19883581059106110.1109/31.1858MR953633PoddarG.ChakrabartyK.BanerjeeS.Control of chaos in the boost converter199531118418422-s2.0-002930666710.1049/el:19950580AkgulA.MorozI.PehlivanI.VaidyanathanS.A new four-scroll chaotic attractor and its engineering applications2016127135491549910.1016/j.ijleo.2016.02.066AkgulA.PehlivanI.A new three-dimensional chaotic system without equilibrium points, its dynamical analyses and electronic circuit application201623120921410.17559/TV-201412121259422-s2.0-84958981179AkgulA.HussainS.PehlivanI.A new three-dimensional chaotic system, its dynamical analysis and electronic circuit applications201612718706270712-s2.0-8497101512610.1016/j.ijleo.2016.05.010KaçarS.Analog circuit and microcontroller based RNG application of a new easy realizable 4D chaotic system201612720955195612-s2.0-8497895912410.1016/j.ijleo.2016.07.044JiaQ.Hyperchaos generated from the Lorenz chaotic system and its control2007366321722210.1016/j.physleta.2007.02.024Zbl1203.930862-s2.0-34248386538SiG.-Q.CaoH.ZhangY.-B.A new four-dimensional hyperchaotic Lorenz system and its adaptive control20112012-s2.0-7995242360510.1088/1674-1056/20/1/010509010509GaoT.ChenZ.YuanZ.ChenG.A hyperchaos generated from Chen's system20061744714782-s2.0-3364654634310.1142/S0129183106008625Zbl1103.37020JiaL.-X.DaiH.HuiM.A new four-dimensional hyperchaotic Chen system and its generalized synchronization201019102-s2.0-7864937694410.1088/1674-1056/19/10/100501100501ChenA.LuJ.LüJ.YuS.Generating hyperchaotic Lü attractor via state feedback control200636410311010.1016/j.physa.2005.09.0392-s2.0-33645149365PangS. Q.LiuY. J.A new hyperchaotic system from the Lü system and its control201123582775278910.1016/j.cam.2010.11.029MR2763185WangX. Y.ZhaoG. B.Hyperchaos generated from the unified chaotic system and its control2010242346194637