# Hopf Bifurcation, Positively Invariant Set, and Physical Realization of a New Four-Dimensional Hyperchaotic Financial System

Authors: G. Kai; W. Zhang; Z. C. Wei; J. F. Wang; A. Akgul
Journal: Mathematical Problems in Engineering
 (2017)

---

## Abstract

This paper introduces a new four-dimensional hyperchaotic financial system on the basis of an established three-dimensional nonlinear financial system and a dynamic model by adding a controller term to consider the effect of control on the system. In terms of the proposed financial system, the sufficient conditions for nonexistence of chaotic and hyperchaotic behaviors are derived theoretically. Then, the solutions of equilibria are obtained. For each equilibrium, its stability and existence of Hopf bifurcation are validated. Based on corresponding first Lyapunov coefficient of each equilibrium, the analytical proof of the existence of periodic solutions is given. The ultimate bound and positively invariant set for the financial system are obtained and estimated. There exists a stable periodic solution obtained near the unstable equilibrium point. Finally, the dynamic behaviors of the new system are explored from theoretical analysis by using the bifurcation diagrams and phase portraits. Moreover, the hyperchaotic financial system has been simulated using a specially designed electronic circuit and viewed on an oscilloscope, thereby confirming the results of the numerical integrations and its real contribution to engineering.


---

## Body

## 1. Introduction

In recent years, economic dynamics become very prominent in the mainstream economics, and stochastic analysis is a popular way to interpret financial time series based on the existing information. Random exogenous shocks is usually considered to be the reason of a periodic behavior of the economic system by employing stochastic analyses, caused by factors outside the system. The shortcoming of the stochastic approach is incapable to illustrate the dynamics of financial systems [1].An effective and rapid control method is very important for the government when some chaotic phenomenon appears [2–5]. A classical three-dimensional financial dynamic model describing the time variations of state adjustment has been reported [3], which includes the interest rate, the price index, and the investment demand. Later, another three-dimensional financial chaotic risk dynamic system was constructed for management process of financial markets and proved to be controlled effectively [6]. Then Holyst and Urbanowicz [7] used the method of delayed feedback control (DFC) and proved that a chaotic financial system can be stabilized on various periodic orbits [8]. Reference [9] discussed the complex behaviors of a financial system with time-delayed feedback by numerical simulation. Normal form of a financial system with delaying has been derived, which is associated with Hopf and double Hopf bifurcations and makes the financial system more complicated [10, 11]. In addition, the positively invariant sets of the dynamic system have a basic significance in the state constraints and control constraints, which are widely applied in different kinds of field: general 3-body problem [12], delay differential-equations [13], stability of dynamic system [14], robust attitude control schemes [15], relative motion control of the spacecraft [16], interconnected and time-delay systems [17], and permanent magnet synchronous motor system [18].On the other hand, high dimensional control system is now a research hot topic, which has been applied in many fields. Zhang et al. introduced a three-dimensional ultimate bound and positively invariant set in which parameters are positive [19]. Wang et al. showed a smooth four-dimensional quadratic autonomous hyperchaotic system, which generated two novel double-wing periodic, quasi-periodic, and hyperchaotic attractors [20]. Wei and Zhang obtained the ultimate bound and positively invariant set is in a four-dimensional autonomous system [21]. Du et al. proposed hyperchaotic Rikitake system which was a novel four-dimensional autonomous nonlinear system and assured the existence of Hopf bifurcation [22]. A new four-dimensional quadratic autonomous hyperchaotic attractor was presented and the Hopf bifurcation at the equilibrium point was analyzed by Prakash and Balasubramaniam [23]. Therefore, the dynamical behaviors of hyperchaotic systems are more complex than chaotic systems and will have practical meanings for carrying out this research about financial systems [24]. Therefore, it is very necessary to explain complicated phenomena of financial dynamics by a thorough studying on the internal structural characteristics of hyperchaotic systems. On the viewpoint of mathematics and finance, these research points are attractive and rational. This paper constructs a new type of hyperchaotic financial system by nonlinear feedback method.The remainder of the paper is organized as follows. Firstly, a new four-dimensional financial system is introduced by adding feedback controllers to the classic financial system in Section2. And a sufficient condition for nonexistence of chaotic or hyperchaotic behaviors is obtained theoretically. Then the novel hyperchaotic financial system is confirmed numerically from Lyapunov exponents. Characterizations for the four-dimensional Hopf bifurcations are surveyed in Section 3. Section 4 introduces the ultimate bound and positively invariant set. The dynamic properties of the system are showed via bifurcation diagram in Section 5. In Section 6, a real contribution to engineering will be realized by an electronic circuit and oscilloscope in real time. Section 7 gives some conclusions.


## 2. Financial System from Classic Financial Model

### 2.1. Formulation of System

A financial dynamic system with different factors is reported in [3]. With the proper dimensions and appropriate coordinates, a simplified financial model is proposed:(1a)x˙=z+y-ax,(1b)y˙=1-by-x2,(1c)z˙=-x-cz,where “x” shows the interest rate that is defined as the price or cost of money for borrowing and return to lending. “y” implies the investment demand, which will be directly affected by the interest rate. In addition, the supply and demand of the commercial goods and inflation rate cause changes of “z”; let a>0, b>0, and c>0, respectively, be the saving amount, the cost per investment, and the elasticity of demand in commercial market.This paper is dedicated to greatly survey the dynamical behaviors of a four-dimensional financial system which expanded from a known three-dimensional financial system and also proposes a simplified mathematical system. In this paper, we design the controlled system as follows:(2a)x˙=z+y-ax,(2b)y˙=1-by-x2,(2c)z˙=-x-cz+u,(2d)u˙=-dxy-ku-mz,where u denotes control input and economically state intervention to balance the economic environment. For example, United States interest rates are determined by the Federal Reserve with considering short term economic targets. Regulatory agencies of open markets meet at certain intervals to monitor the economic and financial situation and decide on monetary policies. To decrease inflation and increase the purchasing power of the consumers, government will increase the interest rate. As a result of this, control input is the factor that interacts with all variables. This cross relation between interest rate, inflation (this also represents the price and goods and services which are denoted as “z” in the study), and government regulations can be expressed from (2d), in whichd, k,andm mean corresponding amplitudes.Therefore, it will be expected to study some complex dynamical behaviors about the proposed four-dimensional autonomous system (2a), (2b), (2c), and (2d).


### 2.2. Nonexistence of Chaotic or Hyperchaotic Behaviors

In addition, chaotic solutions do not exist for certain parameter values in system (2a), (2b), (2c), and (2d). We get three positive Lyapunov exponents, and it will be beneficial for us to find hyperchaos. More precisely, the main results are obtained as follows.Theorem 1.
A six-parameter family(a,b,c,d,k,m) of system (2a), (2b), (2c), and (2d) is considered. If the parameters satisfy a>0, b>0, and c>0, then the following conditions are met as(3)m=d1+ac+2d-a2d+acd+d2-ak-ck+adk-cdk+k21+d2,ad+k1+d<0,and system (2a), (2b), (2c), and (2d) has no bounded chaotic solutions or hyperchaotic solutions.Proof.
From system (2a), (2b), (2c), and (2d), (4)dx˙+dk-a1+dz˙+u˙=-ad-sx+d-m-csz+-k+su.
Under assumption (3), (4) becomes(5)dx˙+dk-a1+dz˙+u˙=-ad+k1+ddx+dk-a1+dz+u.
Thus, we can get the expression that(6)dx+dk-a1+dz+u=dx0+dk-a1+dz0+u0e-ad+k/1+dt.
Since(ad+k)/(1+d)<0 (a>0 is the amount of saving; b>0 is the cost per investment) and at least one of the following conditions is satisfied: x(t), z(t), or u(t) is not bounded, then system (2a), (2b), (2c), and (2d) is not chaotic.
The proof is complete.


## 2.1. Formulation of System

A financial dynamic system with different factors is reported in [3]. With the proper dimensions and appropriate coordinates, a simplified financial model is proposed:(1a)x˙=z+y-ax,(1b)y˙=1-by-x2,(1c)z˙=-x-cz,where “x” shows the interest rate that is defined as the price or cost of money for borrowing and return to lending. “y” implies the investment demand, which will be directly affected by the interest rate. In addition, the supply and demand of the commercial goods and inflation rate cause changes of “z”; let a>0, b>0, and c>0, respectively, be the saving amount, the cost per investment, and the elasticity of demand in commercial market.This paper is dedicated to greatly survey the dynamical behaviors of a four-dimensional financial system which expanded from a known three-dimensional financial system and also proposes a simplified mathematical system. In this paper, we design the controlled system as follows:(2a)x˙=z+y-ax,(2b)y˙=1-by-x2,(2c)z˙=-x-cz+u,(2d)u˙=-dxy-ku-mz,where u denotes control input and economically state intervention to balance the economic environment. For example, United States interest rates are determined by the Federal Reserve with considering short term economic targets. Regulatory agencies of open markets meet at certain intervals to monitor the economic and financial situation and decide on monetary policies. To decrease inflation and increase the purchasing power of the consumers, government will increase the interest rate. As a result of this, control input is the factor that interacts with all variables. This cross relation between interest rate, inflation (this also represents the price and goods and services which are denoted as “z” in the study), and government regulations can be expressed from (2d), in whichd, k,andm mean corresponding amplitudes.Therefore, it will be expected to study some complex dynamical behaviors about the proposed four-dimensional autonomous system (2a), (2b), (2c), and (2d).


## 2.2. Nonexistence of Chaotic or Hyperchaotic Behaviors

In addition, chaotic solutions do not exist for certain parameter values in system (2a), (2b), (2c), and (2d). We get three positive Lyapunov exponents, and it will be beneficial for us to find hyperchaos. More precisely, the main results are obtained as follows.Theorem 1.
A six-parameter family(a,b,c,d,k,m) of system (2a), (2b), (2c), and (2d) is considered. If the parameters satisfy a>0, b>0, and c>0, then the following conditions are met as(3)m=d1+ac+2d-a2d+acd+d2-ak-ck+adk-cdk+k21+d2,ad+k1+d<0,and system (2a), (2b), (2c), and (2d) has no bounded chaotic solutions or hyperchaotic solutions.Proof.
From system (2a), (2b), (2c), and (2d), (4)dx˙+dk-a1+dz˙+u˙=-ad-sx+d-m-csz+-k+su.
Under assumption (3), (4) becomes(5)dx˙+dk-a1+dz˙+u˙=-ad+k1+ddx+dk-a1+dz+u.
Thus, we can get the expression that(6)dx+dk-a1+dz+u=dx0+dk-a1+dz0+u0e-ad+k/1+dt.
Since(ad+k)/(1+d)<0 (a>0 is the amount of saving; b>0 is the cost per investment) and at least one of the following conditions is satisfied: x(t), z(t), or u(t) is not bounded, then system (2a), (2b), (2c), and (2d) is not chaotic.
The proof is complete.


## 3. Some Basic Properties and Bifurcation Analysis of the New System (2a), (2b), (2c), and (2d)

### 3.1. Equilibria and Stability

Firstly, the invariance of system is easily demonstrated with the transformation of coordinate(x,y,z,u)→(-x,-y,z,-ω); namely, the system has rotated symmetry around the axis z.We consider the equilibrium point of system (2a), (2b), (2c), and (2d) to analyze the equilibria and let(7a)z+y-ax=0,(7b)1-by-x2=0,(7c)-x-cz+u=0,(7d)-dxy-ku-mz=0.Combining (7a), (7b), and (7c), we obtain(8)x=±1-by,z=±a-y1-by,u=±1-by1+ac-cy.Substituting (8) into (7d) yields(9)1-byd-ck-my+k+ack+am=0,where fy′=0 and(10)fy′=-b21-byd-ck-my+k+ack+am+1-byd-ck-m=0,(11)y=23b-k+ack+am3d-ck-m.Substituting (11) into (10) yields(12)1-bk+ack+amd-ck-md-ck-mb+k+ack+am=0.Hence, we can obtain the following results:System (2a), (2b), (2c), and (2d) have only one equilibrium point E0(0,1/b,0,0).Ifd-ck-m≠0 and Γ1=-d+ck+m and Γ2=k+ack+am, system (2a), (2b), (2c), and (2d) with 0<Γ1<bΓ2 has three equilibrium points E0(0,1/b,0,0) and E1,2(±x0,y0,∓z0,∓u0), where(13)x0=Γ1-bΓ2Γ1,y0=Γ2Γ1,z0=ad+kΓ1-bΓ2Γ13/2,u0=d+acd-mΓ1-bΓ2Γ13/2.In the next step, we research the stability situation of equilibriaE0 and E1,2. By linearizing system (2a), (2b), (2c), and (2d) at the equilibrium (x∗,y∗,z∗,u∗), the Jacobian matrix is(14)JE=y∗-ax∗10-2x∗-b00-10-c1-dy∗-dx∗-m-k.Obviously, the characteristic equation at the equilibrium point(x∗,y∗,z∗,u∗) is(15)λ4+a+b+c+k-y∗λ3+1+ab+ac+bc+ak+bk+ck+m+2x∗2-by∗-cy∗λ2+b+abc+k+abk+ack+bck+am+bm+2cx∗2-bcy∗+dy∗-bky∗-cky∗-my∗λ+bk+abck+abm-2dx∗2+2ckx∗2+2mx∗2+bdy∗-bcky∗-bmy∗=0.The following symbols are introduced:(16a)a1=a+b+c+k-y∗,(16b)a2=1+ab+ac+bc+ak+bk+ck+m+2x∗2-by∗-cy∗,(16c)a3=b+abc+k+abk+ack+bck+am+bm+2cx∗2-bcy∗+dy∗-bky∗-cky∗-my∗,(16d)a4=bk+abck+abm-2dx∗2+2ckx∗2+2mx∗2+bdy∗-bcky∗-bmy∗.By using the criterion initiated by Routh-Hurwitz, the real parts of all the roots are negative if and only if(17)ai>0,i=1,2,3,4a1a2-a3>0,a1a2a3-a32-a12a4>0.Therefore, the equilibrium point(x∗,y∗,z∗,u∗) is asymptotically stable when (17) are satisfied.Remark 2.
For the sake of simplification of the study of dynamical behavior of system (2a), (2b), (2c), and (2d) in the following subsections, the values of system parameters are fixed as a=0.9, c=1.5, d=0.2, and k=0.05.(A)
Whenm=0.005, system (2a), (2b), (2c), and (2d) has three equilibria for E0 and E1,2:
(A1)
E0 is asymptotically stable when b>0.7068.(A2)
E1,2 is unstable for b∈R+ as a1a2-a3<0.(B)
Whenb=0.2, system (2a), (2b), (2c), and (2d) has three equilibria for E0 and E1,2:
(B1)
E0 is unstable for m∈R+ as a1<0.(B2)
E1,2 is asymptotically stable when m>0.5139.With the target of finding the effect of control parametersb of the new dynamic system of four dimensions, the simulation results are obtained by numerical simulation. According to conditions (A2) and (B1), some dynamic properties of system (2a), (2b), (2c), and (2d) can be analyzed through bifurcation diagrams.


### 3.2. An Analysis of System Bifurcation (2a), (2b), (2c), and (2d)

Now we calculate the first Lyapunov coefficient in relation to Hopf bifurcation by employing the projection method [25]. Let l1 be the first Lyapunov coefficient related to Hopf bifurcation. Firstly, consider the differential equation(18)X˙=fX,μ,in which x∈R4 denotes vectors representing phase variables and μ∈R6 denotes control parameters. Suppose f is a class of C∞ in R4×R6 and there is an equilibrium point in (18) X=X0 at μ=μ0. Let variable X-X0 be X, and we can write(19)FX=fX,μ0as(20)FX=AX+12BX,X+16CX,X,X+OX4,where A=fx(0,μ0) and, for i=1,2,3,(21)BX,Y=∑j,k=13∂2Fiξ∂ξi∂ξkξ=0XjYk,CX,Y,Z=∑j,k,l=13∂3Fiξ∂ξj∂ξk∂ξlξ=0XjYkZl.The following considerations are supposed:A has a pair of complex eigenvalues on the imaginary axis λ2,3=±iω(ω0>0), which are the only eigenvalues with Re⁡λ=0, Tc is the generalized eigenspace of A with regard to λ2,3, and p,q∈C3 are vectors such that(22)Aq=iω0q,ATp=-iω0p,p,q=1,in which AΤ is the transposition of the matrix A.Vectory∈Tc can be showed as y=ωq+ω-q- with ω=p,y∈C. The two-dimensional center manifold related to the eigenvalues λ2,3 is parameterized by ω and ω-, with an immersion of the form X=H(ω,ω-), where H:C2→R3 is a Taylor expansion of the form(23)Hω,ω-=ωq+ω-q-+∑2≤j+k≤31j!k!hjkωjω-k+oω4,with hjk∈C3 and hjk=h-kj. Substituting (23) into (19), the following differential equation could be reached:(24)Hωω′+Hωω-′=FHω,ω-.By solving the system of linear equations defined by the coefficients of (19), we can get the complex vectors hij. Given the coefficient F, system equation (19) can be written on the chart ω for a central manifold as(25)ω˙=iω0ω+12G21ωω2+Oω4,where G21∈C.Denote the first Lyapunov coefficient as(26)l1=12ReG21,in which G21=p,C(q,q,q-)+B(q-,h20)+2B(q,h11).A Hopf bifurcation point(X0,μ0) is an equilibrium point of (18), a pair of purely imaginary eigenvalues ±iω0(ω>0) and another eigenvalue with nonzero real part only exists in the Jacobian matrix A. A two-dimensional center manifold is well defined at a Hopf point, which is invariant under the flow produced by (18) and can continue with arbitrarily high class of differentiability to nearby parameter values.If the parameter-dependent complex eigenvalues cross the imaginary axis with nonzero derivative, then a Hopf point is called transversal. In the area of a transversal Hopf point withl1≠0, the dynamic behavior of system (18) can decrease to a family of parameter-dependent continuations of the center manifold and is topologically orbital equivalent to the complex normal form(27)ω′=η+iωω+l1ωω2,where ω∈C and η, ω, and l1 are real functions with arbitrarily higher order derivatives, which are continuations of 0, ω0, and the first Lyapunov coefficient at the Hopf point [23]. When l1<0(l1>0), we can find one family of stable (unstable) periodic orbits on this manifolded family, contracting to an equilibrium point at the Hopf point.The remaining part of this section employs the four-dimensional Hopf bifurcation theory and uses symbolic computations to carry out the analysis of parametric variations concerning dynamical bifurcations. Because the system has only one equilibrium, the bifurcation of system (2a), (2b), (2c), and (2d) will be our only concern, and then we get Theorem 3.


#### 3.2.1. Hopf Bifurcation atE0

Theorem 3.
With system (2a), (2b), (2c), and (2d) and a=0.9, c=1.5, d=0.2, k=0.05, and m=0.005, the first Lyapunov coefficient at E0 for critical value b=b0=0.7068 is given by(28)l1=-0.1568<0.
Therefore, there exists a transversal Hopf point atE0  of system (2a), (2b), (2c), and (2d), and thus this point is stable. Moreover, for each b<b0, but close to b0, there is a stable limit cycle close to the unstable equilibrium point E0.Proof.
As to the parameters(a,c,d,k,m)=(0.9,1.5,0.2,0.05,0.005) and b=b0=0.7068, we have(29)λ1=-1.0351,λ2=-0.7068,λ3,4=±0.5309i.
Under this condition, it is easy to know the transverse conditionλ′(b=b0)<0. Accordingly, there exists a Hopf bifurcation at E0. The stability of E0 can be decided by the value of the first Lyapunov coefficient l1 and l1 can show the stability of E0 and the occurred periodic orbits. Making use of the notation of the section above, the multilinear symmetrical functions could be written as(30a)Bx,y=x1y3+x3y1,-2x1y1,0,-dx1y2-dx2y1,(30b)Cx,y,z=0,0,0,0.
Moreover, the following results are obtained:(31a)p=0.142675+0.581567i,0,-0.0353347+0.375204i,-0.70668,(31b)q=0.522813+1.01926i,0,-0.810362-0.247239i,-0.561461+0.218144i,(31c)h11=-0.370497,-3.71333,1.54211,1.94266,(31d)h20=-2.77049+1.20144i,-0.726017-1.92509i,0.494102-1.6501i,-0.277134-0.749036i.
Then the following value is computed:(32)G21=-0.313594-2.65646i,l1=12ReG21=-0.1568.
Therefore, the theorem is proved. Near the unstable equilibrium pointE0, we can find a stable periodic solution for b<b0.


#### 3.2.2. Hopf Bifurcation atE1,2

Theorem 4.
With system (2a), (2b), (2c), and (2d) and a=0.9, c=1.5, d=0.2, k=0.05, and b=0.2, the first Lyapunov coefficient at E1,2 for critical value m=m0=0.5139 is given by(33)l1=-0.3045<0.Therefore, system (2a), (2b), (2c), and (2d) has a transversal Hopf point at E1,2, and thus this point is stable. Furthermore, for each m<m0, but close to m0, a stable limit cycle exists near the unstable equilibrium point E1,2.Proof.
Because of symmetry,E1 will be considered. Since the parameters satisfy (a,c,d,k,b)=(0.9,1.5,0.2,0.05,0.2) and m=m0=0.5139, we have(34)λ1=-0.8279,λ2=-0.3307,λ3,4=±1.4120i.
Under this condition, it is easy to know the transverse conditionλ′(m=m0)<0. Accordingly, Hopf bifurcation appears. There exists a Hopf bifurcation at E0. The value of the first Lyapunov coefficient l1 can decide the stability of E0. It demonstrates the stability of the periodic orbits and the equilibrium point. Making use of the notation obtained from the section above, the multilinear symmetric functions could be written as(35a)Bx,y=x1y3+x3y1,-2x1y1,0,-dx1y2-dx2y1,(35b)Cx,y,z=0,0,0,0.
Then, the following results are obtained:(36a)p=0.741494,-0.0915377+0.445225i,0.331105-0.22946i,-0.154015,-0.239956i(36b)q=0.6928-0.277114i,-0.208195-0.851544i,-0.192858+0.428759i,-0.201872+0.0937186i,(36c)h11=0.542254,-1.02525,-0.674678,-0.469763,(36d)h20=0.218437+0.0754352i,0.304126+0.177492i,-0.0578508+0.0199395i,0.075354-0.0580204i.Then the following values are computed:(37)G21=-0.608777+1.36819i,l1=12ReG21=-0.3045.
Therefore, the theorem is justified. We can find a stable periodic solution close to the unstable equilibrium pointE1 for m<m0.


## 3.1. Equilibria and Stability

Firstly, the invariance of system is easily demonstrated with the transformation of coordinate(x,y,z,u)→(-x,-y,z,-ω); namely, the system has rotated symmetry around the axis z.We consider the equilibrium point of system (2a), (2b), (2c), and (2d) to analyze the equilibria and let(7a)z+y-ax=0,(7b)1-by-x2=0,(7c)-x-cz+u=0,(7d)-dxy-ku-mz=0.Combining (7a), (7b), and (7c), we obtain(8)x=±1-by,z=±a-y1-by,u=±1-by1+ac-cy.Substituting (8) into (7d) yields(9)1-byd-ck-my+k+ack+am=0,where fy′=0 and(10)fy′=-b21-byd-ck-my+k+ack+am+1-byd-ck-m=0,(11)y=23b-k+ack+am3d-ck-m.Substituting (11) into (10) yields(12)1-bk+ack+amd-ck-md-ck-mb+k+ack+am=0.Hence, we can obtain the following results:System (2a), (2b), (2c), and (2d) have only one equilibrium point E0(0,1/b,0,0).Ifd-ck-m≠0 and Γ1=-d+ck+m and Γ2=k+ack+am, system (2a), (2b), (2c), and (2d) with 0<Γ1<bΓ2 has three equilibrium points E0(0,1/b,0,0) and E1,2(±x0,y0,∓z0,∓u0), where(13)x0=Γ1-bΓ2Γ1,y0=Γ2Γ1,z0=ad+kΓ1-bΓ2Γ13/2,u0=d+acd-mΓ1-bΓ2Γ13/2.In the next step, we research the stability situation of equilibriaE0 and E1,2. By linearizing system (2a), (2b), (2c), and (2d) at the equilibrium (x∗,y∗,z∗,u∗), the Jacobian matrix is(14)JE=y∗-ax∗10-2x∗-b00-10-c1-dy∗-dx∗-m-k.Obviously, the characteristic equation at the equilibrium point(x∗,y∗,z∗,u∗) is(15)λ4+a+b+c+k-y∗λ3+1+ab+ac+bc+ak+bk+ck+m+2x∗2-by∗-cy∗λ2+b+abc+k+abk+ack+bck+am+bm+2cx∗2-bcy∗+dy∗-bky∗-cky∗-my∗λ+bk+abck+abm-2dx∗2+2ckx∗2+2mx∗2+bdy∗-bcky∗-bmy∗=0.The following symbols are introduced:(16a)a1=a+b+c+k-y∗,(16b)a2=1+ab+ac+bc+ak+bk+ck+m+2x∗2-by∗-cy∗,(16c)a3=b+abc+k+abk+ack+bck+am+bm+2cx∗2-bcy∗+dy∗-bky∗-cky∗-my∗,(16d)a4=bk+abck+abm-2dx∗2+2ckx∗2+2mx∗2+bdy∗-bcky∗-bmy∗.By using the criterion initiated by Routh-Hurwitz, the real parts of all the roots are negative if and only if(17)ai>0,i=1,2,3,4a1a2-a3>0,a1a2a3-a32-a12a4>0.Therefore, the equilibrium point(x∗,y∗,z∗,u∗) is asymptotically stable when (17) are satisfied.Remark 2.
For the sake of simplification of the study of dynamical behavior of system (2a), (2b), (2c), and (2d) in the following subsections, the values of system parameters are fixed as a=0.9, c=1.5, d=0.2, and k=0.05.(A)
Whenm=0.005, system (2a), (2b), (2c), and (2d) has three equilibria for E0 and E1,2:
(A1)
E0 is asymptotically stable when b>0.7068.(A2)
E1,2 is unstable for b∈R+ as a1a2-a3<0.(B)
Whenb=0.2, system (2a), (2b), (2c), and (2d) has three equilibria for E0 and E1,2:
(B1)
E0 is unstable for m∈R+ as a1<0.(B2)
E1,2 is asymptotically stable when m>0.5139.With the target of finding the effect of control parametersb of the new dynamic system of four dimensions, the simulation results are obtained by numerical simulation. According to conditions (A2) and (B1), some dynamic properties of system (2a), (2b), (2c), and (2d) can be analyzed through bifurcation diagrams.


## 3.2. An Analysis of System Bifurcation (2a), (2b), (2c), and (2d)

Now we calculate the first Lyapunov coefficient in relation to Hopf bifurcation by employing the projection method [25]. Let l1 be the first Lyapunov coefficient related to Hopf bifurcation. Firstly, consider the differential equation(18)X˙=fX,μ,in which x∈R4 denotes vectors representing phase variables and μ∈R6 denotes control parameters. Suppose f is a class of C∞ in R4×R6 and there is an equilibrium point in (18) X=X0 at μ=μ0. Let variable X-X0 be X, and we can write(19)FX=fX,μ0as(20)FX=AX+12BX,X+16CX,X,X+OX4,where A=fx(0,μ0) and, for i=1,2,3,(21)BX,Y=∑j,k=13∂2Fiξ∂ξi∂ξkξ=0XjYk,CX,Y,Z=∑j,k,l=13∂3Fiξ∂ξj∂ξk∂ξlξ=0XjYkZl.The following considerations are supposed:A has a pair of complex eigenvalues on the imaginary axis λ2,3=±iω(ω0>0), which are the only eigenvalues with Re⁡λ=0, Tc is the generalized eigenspace of A with regard to λ2,3, and p,q∈C3 are vectors such that(22)Aq=iω0q,ATp=-iω0p,p,q=1,in which AΤ is the transposition of the matrix A.Vectory∈Tc can be showed as y=ωq+ω-q- with ω=p,y∈C. The two-dimensional center manifold related to the eigenvalues λ2,3 is parameterized by ω and ω-, with an immersion of the form X=H(ω,ω-), where H:C2→R3 is a Taylor expansion of the form(23)Hω,ω-=ωq+ω-q-+∑2≤j+k≤31j!k!hjkωjω-k+oω4,with hjk∈C3 and hjk=h-kj. Substituting (23) into (19), the following differential equation could be reached:(24)Hωω′+Hωω-′=FHω,ω-.By solving the system of linear equations defined by the coefficients of (19), we can get the complex vectors hij. Given the coefficient F, system equation (19) can be written on the chart ω for a central manifold as(25)ω˙=iω0ω+12G21ωω2+Oω4,where G21∈C.Denote the first Lyapunov coefficient as(26)l1=12ReG21,in which G21=p,C(q,q,q-)+B(q-,h20)+2B(q,h11).A Hopf bifurcation point(X0,μ0) is an equilibrium point of (18), a pair of purely imaginary eigenvalues ±iω0(ω>0) and another eigenvalue with nonzero real part only exists in the Jacobian matrix A. A two-dimensional center manifold is well defined at a Hopf point, which is invariant under the flow produced by (18) and can continue with arbitrarily high class of differentiability to nearby parameter values.If the parameter-dependent complex eigenvalues cross the imaginary axis with nonzero derivative, then a Hopf point is called transversal. In the area of a transversal Hopf point withl1≠0, the dynamic behavior of system (18) can decrease to a family of parameter-dependent continuations of the center manifold and is topologically orbital equivalent to the complex normal form(27)ω′=η+iωω+l1ωω2,where ω∈C and η, ω, and l1 are real functions with arbitrarily higher order derivatives, which are continuations of 0, ω0, and the first Lyapunov coefficient at the Hopf point [23]. When l1<0(l1>0), we can find one family of stable (unstable) periodic orbits on this manifolded family, contracting to an equilibrium point at the Hopf point.The remaining part of this section employs the four-dimensional Hopf bifurcation theory and uses symbolic computations to carry out the analysis of parametric variations concerning dynamical bifurcations. Because the system has only one equilibrium, the bifurcation of system (2a), (2b), (2c), and (2d) will be our only concern, and then we get Theorem 3.


### 3.2.1. Hopf Bifurcation atE0

Theorem 3.
With system (2a), (2b), (2c), and (2d) and a=0.9, c=1.5, d=0.2, k=0.05, and m=0.005, the first Lyapunov coefficient at E0 for critical value b=b0=0.7068 is given by(28)l1=-0.1568<0.
Therefore, there exists a transversal Hopf point atE0  of system (2a), (2b), (2c), and (2d), and thus this point is stable. Moreover, for each b<b0, but close to b0, there is a stable limit cycle close to the unstable equilibrium point E0.Proof.
As to the parameters(a,c,d,k,m)=(0.9,1.5,0.2,0.05,0.005) and b=b0=0.7068, we have(29)λ1=-1.0351,λ2=-0.7068,λ3,4=±0.5309i.
Under this condition, it is easy to know the transverse conditionλ′(b=b0)<0. Accordingly, there exists a Hopf bifurcation at E0. The stability of E0 can be decided by the value of the first Lyapunov coefficient l1 and l1 can show the stability of E0 and the occurred periodic orbits. Making use of the notation of the section above, the multilinear symmetrical functions could be written as(30a)Bx,y=x1y3+x3y1,-2x1y1,0,-dx1y2-dx2y1,(30b)Cx,y,z=0,0,0,0.
Moreover, the following results are obtained:(31a)p=0.142675+0.581567i,0,-0.0353347+0.375204i,-0.70668,(31b)q=0.522813+1.01926i,0,-0.810362-0.247239i,-0.561461+0.218144i,(31c)h11=-0.370497,-3.71333,1.54211,1.94266,(31d)h20=-2.77049+1.20144i,-0.726017-1.92509i,0.494102-1.6501i,-0.277134-0.749036i.
Then the following value is computed:(32)G21=-0.313594-2.65646i,l1=12ReG21=-0.1568.
Therefore, the theorem is proved. Near the unstable equilibrium pointE0, we can find a stable periodic solution for b<b0.


### 3.2.2. Hopf Bifurcation atE1,2

Theorem 4.
With system (2a), (2b), (2c), and (2d) and a=0.9, c=1.5, d=0.2, k=0.05, and b=0.2, the first Lyapunov coefficient at E1,2 for critical value m=m0=0.5139 is given by(33)l1=-0.3045<0.Therefore, system (2a), (2b), (2c), and (2d) has a transversal Hopf point at E1,2, and thus this point is stable. Furthermore, for each m<m0, but close to m0, a stable limit cycle exists near the unstable equilibrium point E1,2.Proof.
Because of symmetry,E1 will be considered. Since the parameters satisfy (a,c,d,k,b)=(0.9,1.5,0.2,0.05,0.2) and m=m0=0.5139, we have(34)λ1=-0.8279,λ2=-0.3307,λ3,4=±1.4120i.
Under this condition, it is easy to know the transverse conditionλ′(m=m0)<0. Accordingly, Hopf bifurcation appears. There exists a Hopf bifurcation at E0. The value of the first Lyapunov coefficient l1 can decide the stability of E0. It demonstrates the stability of the periodic orbits and the equilibrium point. Making use of the notation obtained from the section above, the multilinear symmetric functions could be written as(35a)Bx,y=x1y3+x3y1,-2x1y1,0,-dx1y2-dx2y1,(35b)Cx,y,z=0,0,0,0.
Then, the following results are obtained:(36a)p=0.741494,-0.0915377+0.445225i,0.331105-0.22946i,-0.154015,-0.239956i(36b)q=0.6928-0.277114i,-0.208195-0.851544i,-0.192858+0.428759i,-0.201872+0.0937186i,(36c)h11=0.542254,-1.02525,-0.674678,-0.469763,(36d)h20=0.218437+0.0754352i,0.304126+0.177492i,-0.0578508+0.0199395i,0.075354-0.0580204i.Then the following values are computed:(37)G21=-0.608777+1.36819i,l1=12ReG21=-0.3045.
Therefore, the theorem is justified. We can find a stable periodic solution close to the unstable equilibrium pointE1 for m<m0.


## 3.2.1. Hopf Bifurcation atE0

Theorem 3.
With system (2a), (2b), (2c), and (2d) and a=0.9, c=1.5, d=0.2, k=0.05, and m=0.005, the first Lyapunov coefficient at E0 for critical value b=b0=0.7068 is given by(28)l1=-0.1568<0.
Therefore, there exists a transversal Hopf point atE0  of system (2a), (2b), (2c), and (2d), and thus this point is stable. Moreover, for each b<b0, but close to b0, there is a stable limit cycle close to the unstable equilibrium point E0.Proof.
As to the parameters(a,c,d,k,m)=(0.9,1.5,0.2,0.05,0.005) and b=b0=0.7068, we have(29)λ1=-1.0351,λ2=-0.7068,λ3,4=±0.5309i.
Under this condition, it is easy to know the transverse conditionλ′(b=b0)<0. Accordingly, there exists a Hopf bifurcation at E0. The stability of E0 can be decided by the value of the first Lyapunov coefficient l1 and l1 can show the stability of E0 and the occurred periodic orbits. Making use of the notation of the section above, the multilinear symmetrical functions could be written as(30a)Bx,y=x1y3+x3y1,-2x1y1,0,-dx1y2-dx2y1,(30b)Cx,y,z=0,0,0,0.
Moreover, the following results are obtained:(31a)p=0.142675+0.581567i,0,-0.0353347+0.375204i,-0.70668,(31b)q=0.522813+1.01926i,0,-0.810362-0.247239i,-0.561461+0.218144i,(31c)h11=-0.370497,-3.71333,1.54211,1.94266,(31d)h20=-2.77049+1.20144i,-0.726017-1.92509i,0.494102-1.6501i,-0.277134-0.749036i.
Then the following value is computed:(32)G21=-0.313594-2.65646i,l1=12ReG21=-0.1568.
Therefore, the theorem is proved. Near the unstable equilibrium pointE0, we can find a stable periodic solution for b<b0.


## 3.2.2. Hopf Bifurcation atE1,2

Theorem 4.
With system (2a), (2b), (2c), and (2d) and a=0.9, c=1.5, d=0.2, k=0.05, and b=0.2, the first Lyapunov coefficient at E1,2 for critical value m=m0=0.5139 is given by(33)l1=-0.3045<0.Therefore, system (2a), (2b), (2c), and (2d) has a transversal Hopf point at E1,2, and thus this point is stable. Furthermore, for each m<m0, but close to m0, a stable limit cycle exists near the unstable equilibrium point E1,2.Proof.
Because of symmetry,E1 will be considered. Since the parameters satisfy (a,c,d,k,b)=(0.9,1.5,0.2,0.05,0.2) and m=m0=0.5139, we have(34)λ1=-0.8279,λ2=-0.3307,λ3,4=±1.4120i.
Under this condition, it is easy to know the transverse conditionλ′(m=m0)<0. Accordingly, Hopf bifurcation appears. There exists a Hopf bifurcation at E0. The value of the first Lyapunov coefficient l1 can decide the stability of E0. It demonstrates the stability of the periodic orbits and the equilibrium point. Making use of the notation obtained from the section above, the multilinear symmetric functions could be written as(35a)Bx,y=x1y3+x3y1,-2x1y1,0,-dx1y2-dx2y1,(35b)Cx,y,z=0,0,0,0.
Then, the following results are obtained:(36a)p=0.741494,-0.0915377+0.445225i,0.331105-0.22946i,-0.154015,-0.239956i(36b)q=0.6928-0.277114i,-0.208195-0.851544i,-0.192858+0.428759i,-0.201872+0.0937186i,(36c)h11=0.542254,-1.02525,-0.674678,-0.469763,(36d)h20=0.218437+0.0754352i,0.304126+0.177492i,-0.0578508+0.0199395i,0.075354-0.0580204i.Then the following values are computed:(37)G21=-0.608777+1.36819i,l1=12ReG21=-0.3045.
Therefore, the theorem is justified. We can find a stable periodic solution close to the unstable equilibrium pointE1 for m<m0.


## 4. Four-Dimensional Estimate of Ultimate Bound and Positively Invariant Set

LetX=[x1,x2,x3,x4]Τ. Define X(t,t0,X0) as the solution to system as(38)X˙=fX,μ,X∈R4,and X(t,t0,X0)=X0 at the initial timet0 and initial state X0, which is denoted as X(t) for simplicity. Ω∈R4 is assumed as a compact set, and the distance between the solution X(t) and the set Ω by ρ(X(t),Ω)=infY∈Ω⁡X(t)-y is defined and denoted as Ωε=X∣ρ(X,Ω)<ε.Definition 5 (see [26, 27]).
Supposing a compact setΩ∈R4. If, for every X0∈R4/Ω, limt→∞ρ(X(t),Ω)=0, that means, for any ε>0, there is T>t0, such that, for t>T, X(t,t0,X0)∈Ωε, and then the set Ω is named an ultimate bound for system (36a), (36b), (36c), and (36d); if for any X0∈Ω and all t≥t0, X(t,t0,X0)⊂Ω, then Ω is called the positively invariant set for system (38).Since quadratic polynomial type Lyapunov functions are simple both in form and in application, we mainly probe into the bounds of chaotic systems by using conicoids. So we choose Lyapunov function as(39)pX=XΤPX+ζX+a0,in which P=(aij)4×4∈R4×4, ζ=(ζ1,ζ2,ζ3,ζ4)∈R4×4, and a0∈R are parameters to determine. Deriving p(X) along system (38) and deleting the terms with three or higher degrees, we get(40)p˙=XΤQX+ηX+b0.Here Q=(bij)4×4∈R4×4, η=(η1,η2,η3,η4)∈R4×4, and b0∈R.IfP is positive definite and Q is negative definite, then p˙=0 must be a bounded sphere in R4; and now that chaotic systems are bounded, p(X) reaches the maximum or minimum values under the solution set of (18), whereas it is necessary that p˙=0. So, in order to assess the ultimate bounds of system equation (18), the following optimization problem should only be solved:(41)MaximumorminimumpX=XΤPX+ζX+a0XΤQX+ηX+b0=0.If the above discussed conditions are satisfied, then the above optimization problem has a solution supposingpmin≤p(X)≤pmax, and then set Ω=X∈R4∣pmin≤pX≤pmax would be our desired ultimate bounds.Making use of the above method, we talk about the ultimate bound and positively invariant set for system (2a), (2b), (2c), and (2d) with d=0.Theorem 6.
DenoteΩ=(x,y,z,u)∣mx2+my-h2+mz2+u2≤R2, where(42)R2=m4ab-a,b≥2am4kb-k,b≥2km4cb-c,b≥2cmb2,b<2a,b<2k,b<2c.If parameters of the system satisfym>0, a>0, b>0, c>0, and k>0, then system (38) would have an ultimate bound and positively invariant set: hyperellipsoid Ω.Proof.
First the following generalized positive definite and radially unbounded Lyapunov function are defined as(43)Vx,y,z,u=mx2+my2+mz2+u2,where m>0 and h∈R. Calculating the derivative of V(x,y,z,u) along with a track of (2a), (2b), (2c), and (2d), we have(44)12Vx,y,z,udt=-amx2-mby-12b2-cmz-ku2+m4b.
It is easy to see that, form>0, a>0, b>0, c>0, and k>0, V(x,y,z,u) is positive definite and the quadratic principal part of V(x,y,z,u)/dt is negative definite. V(x,y,z,u)=0, which indicates the surface(45)Γ=x,y,z,u∣amx2+mby-12b2+cmz+ku2=m4bis an ellipsoid in the space of four dimensions for certain σ, k, b, d, and k. Outside of Γ, V(x,y,z,u)/dt<0, and, inside of Γ, V(x,y,z,u)/dt>0. Thus, system (2a), (2b), (2c), and (2d) would reach its ultimate bound on Γ.
Next we can obtain the maximum ofV on Γ by expressions (2a), (2b), (2c), and (2d) according to employing the Lagrange multiplier method. Define(46)F=mx2+my-h2+mz2+u2+τamx2+mby-12b2+cmz+ku2-m4b,(47a)12Fx′=m1+aτx=0,(47b)12Fy′=my+τby-12=0,(47c)12Fz′=m1+cτz=0,(47d)12Fu′=1+kτu=0,(47e)12Fτ′=amx2+mby-12b2+cmz+ku2-m4b=0.(1)
Ifτ≠-1/a, τ≠-1/c, and τ≠-1/k, we have (x,y,z,u)=(0,0,0,0) or (0,1/b,0,0) and(48a)Vx,y,z,umax=mb2.
(2)
Ifτ=-1/a, τ=-1/c, τ≠-1/k, and b≥2a, we have (x,y,z,u)=(±1/2b-2a/aa-b2,1/2(b-a),0,0) and (48b)Vx,y,z,umax=m4ab-a.
(3)
Ifτ=-1/c, τ=-1/a, τ≠-1/k, and b≥2c, we have (x,y,z,u)=(0,1/2(b-c),±1/2(c-b)(b-2c)/c,0) and (48c)Vx,y,z,umax=m4ab-a.
(4)
Ifτ=-1/k, τ≠-1/a, τ≠-1/c, and b≥2a, we have (x,y,z,u)=(0,1/2(b-k),0,±1/2(k-b)m(b-2k)/k) and (48d)Vx,y,z,umax=m4kb-k.
(5)
Ifτ=-1/a, a=k, k≠c, and b≥2a, we have (y,u)=(1/2(b-a),0), x2+z2=(b-2a)/4a(a-b)2, and(48e)Vx,y,z,umax=m4ab-a.
(6)
Ifτ=-1/a, a=c, k≠c, and b≥2a, we have (y,u)=(1/2(b-a),0), mx2+u2=m(b-2a)/4a(a-b)2, and(48f)Vx,y,z,umax=m4ab-a.
(7)
Ifτ=-1/k, a=k=q, k≠a, and b≥2k, we have y=1/2(b-k), x2+z2=(b-2a)/4a(a-b)2, and (48g)Vx,y,z,umax=m4kb-k.
(8)
Ifτ=-1/k, a=c, k≠c, and b≥2a, we have (y,u)=(1/2(b-a),0), u2+mx2+mz2=m(b-2k)/4k(b-k)2, and (48h)Vx,y,z,umax=m4kb-k.Figure6 performs the corresponding ultimate bound and positively invariant set Ω with u=0 or x=0 (a=0.9, b=0.2, c=1.5, d=0, k=0.05, and m=0.005). We see that the estimated bounds in Theorem 6 are feasible.


## 5. Numerical Simulation

In the previous twenty years, different tools were developed to tap the periodic responses of nonlinear dynamical systems. The periodic responses can be spotted by several traditional criteria such as bifurcation diagrams and phase portraits. In this section, we use these techniques to illustrate the existence of the periodic motions for the financial system. Based on the four-dimensional nonlinear financial system (2a), (2b), (2c), and (2d), numerical simulations have been conducted in this paper. Using MATLAB software, a numerical approach is utilized to show the nonlinear dynamic behavior in the financial system.The red in Figure1 shows the hyperchaotic attractors for the four-dimensional financial system (2a), (2b), (2c), and (2d) with parameter values a=0.9, b=0.2, c=1.5, d=0.2, k=0.05, and m=0.005 and the initial condition (0,1,-0.5,0). Figure 1(a) is the phase portrait in (x,y,z) space and u=0, and Figure 1(b) is the phase portrait in (y,z,u) space and x=0. It shows that the corresponding Lyapunov exponents are L1=0.03003, L2=0.01448, L3=-0.0003, and L4=-1.2318, and the system has complex dynamics. Figure 2 implies the bifurcation set of system (2a), (2b), (2c), and (2d) when the values of parameters are set as a=0.9, c=1.5, d=0.2, k=0.05, b∈[-1,1], and m∈[-1,1], and the bifurcation set divides the space into three regions with three solutions above the bifurcation set, one solution in the bifurcation set, and none under the set, which is consistent with the theoretical results. Figure 3 illustrates the bifurcation curve when a=0.9, c=1.5, d=0.2, k=0.05, and m=0.005, and the point of Hopf bifurcation is b=0.68.Figure 1
Hyperchaotic attractor for the four-dimensional financial system (2a), (2b), (2c), and (2d) with parameter values a=0.9, b=0.2, c=1.5, d=0.2, k=0.05, and m=0.005 and the initial condition (0,1,-0.5,0). (a) The phase portrait in (x,y,z) space and u=0; (b) the phase portrait in (y,z,u) space and x=0.
(a)
(b)Figure 2
The bifurcation set of system (2a), (2b), (2c), and (2d) with parameter values a=0.9, c=1.5, d=0.2, and k=0.05, and b∈[-1,1] and m∈[-1,1].Figure 3
The bifurcation curve of system (2a), (2b), (2c), and (2d) with parameter values a=0.9, c=1.5, d=0.2, k=0.05, and m=0.005, and the bifurcation point occurs when b=0.68.The bifurcation curve in Figure4 shows the bifurcation with a=0.9, c=1.5, d=0.2, k=0.05, b=0.2, and altering system parameter m, and the point of Hopf bifurcation is m=0.49. Figure 5(a) shows that there is a stable periodic solution near the unstable equilibrium point E0 for (a,c,d,k)=(0.9,1.5,0.2,0.05) and b=0.68. There is a stable periodic solution near the unstable equilibrium point E1 for m=0.21 shown in Figure 5(b). Figure 6 shows the corresponding ultimate bound and positively invariant set Ω with u=0 or x=0 (a=0.9, b=0.3, c=1.5, d=0, k=0.05, and m=0.005). Figure 6(a) is the phase portrait in (x,y,z) space and u=0. Figure 6(b) is the phase portrait in (y,z,u) space and x=0. The trajectories of system (2a), (2b), (2c), and (2d) are all in the domain of the estimated bound.Figure 4
The bifurcation curves of system (2a), (2b), (2c), and (2d) with parameter values a=0.9, c=1.5, d=0.2, k=0.05, and b=0.2, and the bifurcation point occurs when m=0.21.Figure 5
The stable periodic solution of system (2a), (2b), (2c), and (2d) with (a,c,d,k)=(0.9,1.5,0.2,0.05) and (a) b=0.68; (b) m=0.49.
(a)
(b)Figure 6
The projections of the bound estimated for system (2a), (2b), (2c), and (2d) with a=0.9, b=0.3, c=1.5, d=0, k=0.05, and m=0.005, (a) the phase portrait in (x,y,z) space and u=0, and (b) the phase portrait in (y,z,u) space and x=0.
(a)
(b)


## 6. Physical Realization for the 4D Financial System (2a), (2b), (2c), and (2d)

In this section, the 4D hyperchaotic financial system will be applied for its application using the electronic circuit which has been designed on the oscilloscope. The hyperchaotic system is studied by using an electronic circuit design. The numerical simulation and the oscilloscope outputs were obtained for the similar shaped phase portraits [28–30].The electronic circuit of the 4D financial system was designed in OrCAD-PSpice. Figure7 shows the electronic circuit schematic of the system. Experimental electronic circuit of the 4D financial system was implemented for parameters a=0.9, c=1.5, d=0.2, k=0.05, and m=0.005 and initial conditions (0,1,-0.5,0).R1= 443 Kohm,R2 =R7 =R9 = 400 Kohm,R3 =R6 = 40 Kohm,R4 = 6000 Kohm,R5 = 200 Kohm,R8 = 266 Kohm,R10 =R11 =R15 =R16 = 100 Kohm,R12 = 8000 Kohm,R13 = 800 Kohm,R14 = 20 Kohm,C1 =C2 =C3 = 1 nF,C4 = 10 nF,Vn = −15 V, andVp = 15 V were chosen. Real-time application of the 4D financial system was realized with electronic components on electronic card in Figure 8. The OrCAD-PSpice simulation outputs and oscilloscope outputs of 4D financial system, for parameters a=0.9, c=1.5, d=0.2, k=0.05, and m=0.005, are seen in Figures 9 and 10.Figure 7
The electronic circuit schematic of the 4D financial system (2a), (2b), (2c), and (2d).Figure 8
The experimental circuit of the hyperchaotic circuit.Figure 9
The phase portraits of 4D financial system (2a), (2b), (2c), and (2d) with a=0.9, b=0.2, c=1.5, d=0.2, k=0.05, and m=0.005 in ORCAD-PSpice.Figure 10
The phase portraits of 4D financial system (2a), (2b), (2c), and (2d) with a=0.9, b=0.2, c=1.5, d=0.2, k=0.05, and m=0.005 on the oscilloscope.


## 7. Conclusion

In our paper, the proposed financial system, which is generated from three-dimensional classical financial system, greatly expands the list of hyperchaotic financial attractors. The sufficient conditions for nonexistence of chaotic and hyperchaotic behaviors are derived theoretically, and the solutions of equilibria are obtained. For each equilibrium, its stability and existence of Hopf bifurcation are validated. Besides, based on its corresponding first Lyapunov coefficient, the analytic proof of the existence of periodic solutions is exhibited. Then the research has got the ultimate bound and positively invariant set for the proposed hyperchaotic financial system. The direction and stability of the bifurcating periodic solutions can be determined, and some numerical solutions are obtained to verify the theoretical results. In addition, the hyperchaotic financial system is extended to an electronic circuit implementation for real-time application.We can understand the functions of financial policies and reveal the true geometrical structure of the attractors by the results of this paper. Since the global dynamics and geometrical structure of this system are not presented completely, more detailed theoretical simulation and investigations are expected in the forthcoming study.



---

*Source: 2490580-2017-04-10.xml*