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.

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 [25]. 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 Section 2. 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:where “” shows the interest rate that is defined as the price or cost of money for borrowing and return to lending. “” 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 “”; let , , and , 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:where 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 which d, k, and m 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 of system (2a), (2b), (2c), and (2d) is considered. If the parameters satisfy , , and , then the following conditions are met asand system (2a), (2b), (2c), and (2d) has no bounded chaotic solutions or hyperchaotic solutions.

Proof. From system (2a), (2b), (2c), and (2d), Under assumption (3), (4) becomesThus, we can get the expression thatSince ( is the amount of saving; is the cost per investment) and at least one of the following conditions is satisfied: , , or 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 ; namely, the system has rotated symmetry around the axis .

We consider the equilibrium point of system (2a), (2b), (2c), and (2d) to analyze the equilibria and let

Combining (7a), (7b), and (7c), we obtainSubstituting (8) into (7d) yieldswhere andSubstituting (11) into (10) yields

Hence, we can obtain the following results:

System (2a), (2b), (2c), and (2d) have only one equilibrium point .

If and and , system (2a), (2b), (2c), and (2d) with has three equilibrium points and , where

In the next step, we research the stability situation of equilibria and . By linearizing system (2a), (2b), (2c), and (2d) at the equilibrium , the Jacobian matrix is

Obviously, the characteristic equation at the equilibrium point isThe following symbols are introduced:

By using the criterion initiated by Routh-Hurwitz, the real parts of all the roots are negative if and only if

Therefore, the equilibrium point 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 , , , and .When , system (2a), (2b), (2c), and (2d) has three equilibria for and : is asymptotically stable when . is unstable for as .When , system (2a), (2b), (2c), and (2d) has three equilibria for and : is unstable for as . is asymptotically stable when .

With the target of finding the effect of control parameters of the new dynamic system of four dimensions, the simulation results are obtained by numerical simulation. According to conditions and , 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 be the first Lyapunov coefficient related to Hopf bifurcation. Firstly, consider the differential equationin which denotes vectors representing phase variables and denotes control parameters. Suppose is a class of in and there is an equilibrium point in (18) at . Let variable be , and we can writeaswhere and, for ,

The following considerations are supposed: has a pair of complex eigenvalues on the imaginary axis , which are the only eigenvalues with , is the generalized eigenspace of with regard to , and are vectors such thatin which is the transposition of the matrix .

Vector can be showed as with . The two-dimensional center manifold related to the eigenvalues is parameterized by and , with an immersion of the form , where is a Taylor expansion of the formwith and . Substituting (23) into (19), the following differential equation could be reached:By solving the system of linear equations defined by the coefficients of (19), we can get the complex vectors . Given the coefficient , system equation (19) can be written on the chart for a central manifold aswhere .

Denote the first Lyapunov coefficient asin which .

A Hopf bifurcation point is an equilibrium point of (18), a pair of purely imaginary eigenvalues and another eigenvalue with nonzero real part only exists in the Jacobian matrix . 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 with , 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 formwhere and , , and are real functions with arbitrarily higher order derivatives, which are continuations of , , and the first Lyapunov coefficient at the Hopf point [23]. When , 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 at

Theorem 3. With system (2a), (2b), (2c), and (2d) and , , , , and , the first Lyapunov coefficient at for critical value is given byTherefore, there exists a transversal Hopf point at   of system (2a), (2b), (2c), and (2d), and thus this point is stable. Moreover, for each , but close to , there is a stable limit cycle close to the unstable equilibrium point .

Proof. As to the parameters and , we haveUnder this condition, it is easy to know the transverse condition . Accordingly, there exists a Hopf bifurcation at . The stability of can be decided by the value of the first Lyapunov coefficient and can show the stability of and the occurred periodic orbits. Making use of the notation of the section above, the multilinear symmetrical functions could be written asMoreover, the following results are obtained:Then the following value is computed:Therefore, the theorem is proved. Near the unstable equilibrium point , we can find a stable periodic solution for .

3.2.2. Hopf Bifurcation at

Theorem 4. With system (2a), (2b), (2c), and (2d) and , , , , and , the first Lyapunov coefficient at for critical value is given byTherefore, system (2a), (2b), (2c), and (2d) has a transversal Hopf point at , and thus this point is stable. Furthermore, for each , but close to , a stable limit cycle exists near the unstable equilibrium point .

Proof. Because of symmetry, will be considered. Since the parameters satisfy and , we haveUnder this condition, it is easy to know the transverse condition . Accordingly, Hopf bifurcation appears. There exists a Hopf bifurcation at . The value of the first Lyapunov coefficient can decide the stability of . 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 asThen, the following results are obtained:Then the following values are computed:Therefore, the theorem is justified. We can find a stable periodic solution close to the unstable equilibrium point for .

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

Let . Define as the solution to system asand at the initial time t0 and initial state , which is denoted as for simplicity. is assumed as a compact set, and the distance between the solution and the set by is defined and denoted as .

Definition 5 (see [26, 27]). Supposing a compact set . If, for every , , that means, for any , there is , such that, for , , and then the set is named an ultimate bound for system (36a), (36b), (36c), and (36d); if for any and all , , 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 asin which , , and are parameters to determine. Deriving along system (38) and deleting the terms with three or higher degrees, we getHere , , and .

If is positive definite and is negative definite, then must be a bounded sphere in ; and now that chaotic systems are bounded, reaches the maximum or minimum values under the solution set of (18), whereas it is necessary that . So, in order to assess the ultimate bounds of system equation (18), the following optimization problem should only be solved:

If the above discussed conditions are satisfied, then the above optimization problem has a solution supposing , and then set 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 .

Theorem 6. Denote , where

If parameters of the system satisfy , , , , and , 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 where and . Calculating the derivative of along with a track of (2a), (2b), (2c), and (2d), we haveIt is easy to see that, for , , , , and , is positive definite and the quadratic principal part of is negative definite. , which indicates the surfaceis an ellipsoid in the space of four dimensions for certain , , , , and . Outside of , , and, inside of , . Thus, system (2a), (2b), (2c), and (2d) would reach its ultimate bound on Γ.
Next we can obtain the maximum of on by expressions (2a), (2b), (2c), and (2d) according to employing the Lagrange multiplier method. Define

(1) If , , and , we have or and(2) If , , , and , we have and (3)If , , , and , we have and (4)If , , , and , we have and (5)If , , , and , we have , , and(6)If , , , and , we have , , and (7) If , , , and , we have , , and (8)If , , , and , we have , , and

Figure 6 performs the corresponding ultimate bound and positively invariant set with or (, , , , , and ). 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 Figure 1 shows the hyperchaotic attractors for the four-dimensional financial system (2a), (2b), (2c), and (2d) with parameter values , , , , , and and the initial condition . Figure 1(a) is the phase portrait in space and , and Figure 1(b) is the phase portrait in space and . It shows that the corresponding Lyapunov exponents are , , , and , 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 , , , , , and , 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 , , , , and , and the point of Hopf bifurcation is .

The bifurcation curve in Figure 4 shows the bifurcation with , , , , , and altering system parameter , and the point of Hopf bifurcation is . Figure 5(a) shows that there is a stable periodic solution near the unstable equilibrium point for and . There is a stable periodic solution near the unstable equilibrium point for shown in Figure 5(b). Figure 6 shows the corresponding ultimate bound and positively invariant set with or (, , , , , and ). Figure 6(a) is the phase portrait in space and . Figure 6(b) is the phase portrait in space and . The trajectories of system (2a), (2b), (2c), and (2d) are all in the domain of the estimated bound.

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

The electronic circuit of the 4D financial system was designed in OrCAD-PSpice. Figure 7 shows the electronic circuit schematic of the system. Experimental electronic circuit of the 4D financial system was implemented for parameters , , , , and and initial conditions . 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, and Vp = 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 , , , , and , are seen in Figures 9 and 10.

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.

Conflicts of Interest

The authors declare that they have no conflicts of interest.


The authors gratefully acknowledge the support of National Natural Science Foundation of China (NNSFC) through Grant nos. 11290152, 11427801, and 11401543, the Funding Project for Academic Human Resources Development in Institutions of Higher Learning under the Jurisdiction of Beijing Municipality (PHRIHLB), Beijing Postdoctoral Research Foundation (no. 2015ZZ17), the China Postdoctoral Science Foundation funded project (no. 2014M560028 and no. 2015T80029), the Government of Chaoyang District Postdoctoral Research Foundation (no. 2015ZZ-7), and Sakarya University Scientific Research Projects Unit under Grants 2016-09-00-008 and 2016-50-01-026.