Qualitative Study of a 4 D Chaos Financial System

College of Mathematics and Statistics, Chongqing Technology and Business University, Chongqing 400067, China Mathematical Postdoctoral Station, College of Mathematics and Statistics, Southwest University, Chongqing 400716, China School of Mathematics and Statistics, Ankang University, Shaanxi 725000, China College of Mathematics and Statistics, Henan Polytechnic Institute, Nanyang 473000, China College of Electronic and Information Engineering, Southwest University, Chongqing 400716, China


Introduction
Chaotic systems are characterized by their extreme sensitivity both to initial conditions and to system parameters.The study of various aspects of chaotic systems has received great interest among scientists from various fields due to their numerous potential applications in science and engineering .Among such aspects, the problem of investigating the ultimate boundedness of chaotic systems and hyperchaotic systems is an important subject.Since the discovery of the famous Lorenz chaotic attractor [1], ultimate boundedness of the Lorenz attractor has been investigated by Leonov et al. in a series of articles [7,10].Since then many papers have studied ultimate boundedness of other chaotic systems [25][26][27][28][29][30][31][32].However, the approach taken in each is only suitable for that particular chaotic system.It is very difficult to propose a universal approach to estimate the bounds for an arbitrary chaotic system.Zhang et al. studied the ultimate boundedness of the Lü system and the Chen system by using the Lyapunov stability theory and optimization method [25][26][27].Particularly, Liao et al. studied the ultimate boundedness of the Yang-Chen system by using geometric and algebraic methods [31].Zhang et al. studied the ultimate boundedness of a novel finance chaotic system by using the Lyapunov stability theory and iterative method [32].To this end, it is necessary to study the bounds of the new chaotic systems.
Recently, a financial chaos dynamical system is reported as follows [32,33]: System (1) models a financial dynamical system composed of product, money, bond, and labor force.The variables x, y, and z denote the interest rate, the investment demand, and the price index, respectively.The positive parameters a, b, and c denote the saving amount, the per investment cost, and the demand elasticity of commercials, respectively.The factors that induce the change of the interest rate x mainly come from two aspects: the price index and the surplus between the investment demand and the savings amount.The changing rate of y is determined by the benefit rate of investment (we assume that the rate is constant during a certain period), the feedback of the investment demand, and the interest rate.The change of the price index z is controlled by the real interest rate and price index.
According to the financial dynamical system (1), by adding a variable u, one gets a new 4D financial chaos dynamical system as follows: where the variable u denotes control input and economically state intervention to balance the economic environment and a, b, c, d, k, m are positive real parameters of system (2).The Lyapunov exponents of the financial system (2) are calculated numerically for the parameter values a = 0 9, b = 0 2, c = 1 5, d = 1, k = 0 05, and m = 0 005 with the initial state x 0 , y 0 , z 0 , u 0 = 0, 1, … , 0 5, 0 .System (2) has the Lyapunov exponents as λ LE 1 = 0 1740, λ LE 2 = 0 1314, λ LE 3 = 0 0000, and λ LE 4 = −15 6059 for the above parameters (see [8,9] for a detailed discussion of Lyapunov exponents of strange attractors in chaos dynamical systems).
Remark 1.The chaotic attractor is one of the important concepts of dynamical systems.Although system (2) has a chaotic attractor for the positive parameters a = 0 9, b = 0 2, c = 1 5, d = 1, k = 0 05, and m = 0 005, the type of the chaotic attractor of system (2) is still unknown.It is interesting to discuss whether the chaotic attractor of system (2) is hidden or self-excited in the future (see the excellent papers [17,19] for a detailed discussion of the chaotic attractor of dynamical systems).
This paper is organized as follows.In the next section, we will study the bounds for the chaotic attractors (2) using the Lyapunov stability theory and optimization method.To validate the ultimate bound estimation, numerical simulations are also provided.Finally, we will give some concluding remarks in Section 3.

Main Properties
Let X t = x t , y t , z t , u t be an arbitrary solution of system (2).Then, the following set is the ultimate bound set and positively invariant set of system (2), where Proof 1. Define the Lyapunov-like function where X t = x t , y t , z t , u t .Computing the derivative of V X along the trajectory of (2), we have Obviously, the surface Γ 1 , that is, defined by is an ellipsoid in R 4 ∀a > 0, b > 0, c > 0, k > 0, d ≠ 0, and m > 0.
Outside Γ 1 : dV X /dt < 0, while inside Γ 1 : dV X /dt > 0. Thus, the maximum value of V X can only be reached on Γ 1 .Since the V X is a continuous function and Γ 1 is a bounded closed set, then the function ( 5) can reach its maximum value max (2).By solving the following conditional extremum problem, one can get the maximum value of the function ( 5) as follows: That is to say, Let us take mx = x 1 , my = y 1 , mz = z 1 , du = u 1 as the new variables, then (9) transforms into the following form: According to the Lagrange multiplier method, we can easily get the above conditional extremum problem (10) as Finally, it is easy to show that (3) is the ultimate bound set and positively invariant set of system (2).This completes the proof.
Figure 2 shows the chaotic attractor of system (2) in the xOyz space defined by Ω 1 .
(ii) We can figure out that the volume of the set Though Theorem 1 gives the ultimate bound set and positively invariant set of the financial chaos system (2), it does not give the global exponential attractive set of system (2).The global exponential attractive set of system (2) is described by the following Theorem 2.
Theorem 2. Suppose that a > 0, b > 0, c > 0, k > 0, d ≠ 0, and m > 0, and let Then, the estimation, holds for system (2), and thus which is the globally exponential attractive set of system (2), that is, lim t→+∞ V X t ≤ L 2 .Define the Lyapunov function Differentiating the Lyapunov function V X in (18) with respect to time t along the trajectory of system (2), when V X t > L 2 and V X t 0 > L 2 , we have Thus, we have which clearly shows that Ψ = X | V X ≤ L 2 is the globally exponential attractive set of system (2).
The proof is complete.
(i) We can figure out that the volume of the set Ψ in (15) is v Ψ = π 2 L 4 /2m md according to Theorem 2.

Conclusions
In this paper, we considered some dynamics of a new 4D financial chaos system.We obtained the ultimate boundedness and global exponential attractive sets of this system according to Lyapunov stability theory.MATLAB simulations show that the proposed method is effective.Characteristic time scales of the 4D financial system, the homoclinic orbits, and heteroclinic orbits of the 4D financial system will be taken into consideration in the future.