An Analysis on the Fractional Asset Flow Differential Equations

The asset flow differential equation (AFDE) is the mathematical model that plays an essential role for planning to predict the financial behavior in the market. In this paper, we introduce the fractional asset flow differential equations (FAFDEs) based on the Liouville-Caputo derivative. We prove the existence and uniqueness of a solution for the FAFDEs. Furthermore, the stability analysis of the model is investigated and the numerical simulation is accordingly performed to support the proposed model.


Introduction
In recent years, the important thing for financial behavior has been understanding and classifying the dynamics of the financial system such as price crashes, bubbles, momentum and liquidity.Over the last three decades, the mathematical model has played an essential role in preparing planning to predict the financial behavior in the market in the future.Thus, a number of mathematical models have been used to predict financial behavior [1][2][3][4][5][6][7][8][9][10][11][12].Moreover, in order to obtain better simulation of future phenomena, nonlinear continuous models have also developed.Researchers proposed several nonlinear continuous models for describing the complicated dynamics of the financial system.For example, in 2008, Fanti et al. [10] proposed the behavior of capital stock by using a dynamic Investment Saving-Liquidity Preference Money Supply (IS-LM) model with delay, and Chatterjee et al. [11] proposed a dynamic Hecker-Ohlin model and considered the equilibrium of the model.Moreover, Yang et al. [12] proposed a herd mechanism of an open financial market for considering volatility and parameters of bubbles and collapses in the market.In this paper, we study the analysis of a nonlinear system that is the asset flow differential equations (AFDEs).AFDEs were proposed and developed by Gaginalp et al. [13].They investigated the solution of autonomous equations with dimension n ≥ 3.This model may describe many behaviors in real financial markets [13][14][15][16].The model variables are explained by the demand, the supply, the market price, the transition rate, the trend-based and value-based of the investor preferences at sudden times [13,17,18].
Fractional calculus was introduced for over three centuries.At the beginning, fractional calculus was only in a theoretical sense.The traditional system of integer-order of differential equations may fail to explain complicated incidents in real phenomena.However, nowadays, fractional order of differential equations has been applied to model the complicated real situations involving many areas-for example, physics, engineering, epidemic, finance and sciences [19][20][21][22][23][24].The concept of a fractional model helps with expressing the data structure for real life situations more than the integer model.The main dominant point of fractional order equations compared to integer order equations is its memory, which is the history of the system [25][26][27].
Since variables in the financial system have long memories, the fractional order system is suitable for studying the behavior in financial systems.There are many papers proposing fractional-order differential equations for the financial system [25][26][27][28] consisting of the parts of products, money, bond and labor force.Here, we propose a system of time-fractional differential equations modified from the Gaginalp's model [13] to describe the behavior of asset flow in the financial market.
A number of researchers studied about the existence and uniqueness for the solution of fractional differential equations [29][30][31][32][33][34].However, there has been a lack of research focusing on the global existence and uniqueness for the solution of the fractional asset flow differential equations (FAFDEs).In this paper, we propose an alternative way to prove the local existence solution of the system of fractional-order differential equations by the Banach fixed point theorem, which helps to explain the solution of the fractional model.Moreover, the global existence for the solution of the fractional model can be express by continuation theorem.
The remainder of this work is organized as follows.In Section 2, the definitions and lemma involving the fractional-order differential equations are presented.In addition, we present the fractionalorder model in this section.In Section 3, the local existence and uniqueness theorem of the model's solution is presented.In Section 4, we propose the continuation theorems for the model to extend the local solution into the global solution.The stability analysis of the fractional model is presented and is followed by the numerical simulation in Section 6.Finally, section 7 concludes and discusses this work.

Fractional Calculus
Definitions of fractional integral and derivatives [19,35] are as follows : Definition 1.The Riemann-Liouville integral of order 0 < α < 1 for a function f is defined as: where Γ is the gamma function.
Definition 2. The Riemann-Liouville derivative of order 0 < α < 1 for a function f : (0, ∞) → R is defined as: (2) The regularization definition of the Riemann-Liouville derivative was introduced by Liouville-Caputo in 1967 [36].Definition 3. The Liouville-Caputo-type fractional derivative of order 0 < α < 1 for a function f : (0, ∞) → R is defined by: The most regular definition is the Liouville-Caputo-type definition because the initial condition for fractional order differential equations with Liouville-Caputo-type derivatives is similar in form to the integer-order differential equation [19][20][21][22].Lemma 1. (Banach fixed point theorem) Let A be a closed convex bounded subset of a Banach space X and V be a contraction mapping from A into A.Then, there exists a unique z ∈ A such that V(z) = z.

Fractional Model
The behavior of price dynamic in the market based on: x 1 (t) the market price of an asset at time t, x 2 (t) the fraction of total asset invested in the stock at time t, x 3 (t) the trend-based component of the investor preference at time t, x 4 (t) the value-based component of the investor preference at time t, P a the fundamental value, D the demand, S the supply.
The value of the market price of an asset at time t can be only a positive number and cannot be approached to infinity, the value of the fraction of total assets is defined between 0 and 1.The demand (D) and the supply (S) can be determined as D = k(t)(1 − x 2 (t)) and S = (1 − k(t))x 2 (t), respectively.k(t) is the transition rate that is dependent on x 3 (t) and x 4 (t) defined by k(t) = 0.5 + 0.5 tanh(x 3 (t) + x 4 (t)), F(t) is an increasing continuous function with F(1) = 0, and P a , a 1 , a 2 , a 3 and a 4 are positive constants.
In this article, the FAFDEs based on the Liouville-Caputo fractional derivative are introduced as follows: where D α t is the Liouville-Caputo fractional derivative of order 0 < α < 1, b j are positive constants for j = 1, 2, 3, 4, 5 and x i (0) = x i,0 are positive constants for i = 1, 2, 3, 4.
We need the following assumptions; let X be a Banach space with its norm || • || X .

Local Existence and Uniqueness of a Solution
In this section, we show that system (5) has a unique continuous solution by the Banach fixed point theorem.To do so, we introduce a Banach space: with T being a positive constant.
Theorem 1.There exists a finite time h such that system (5) has a unique continuous solution on [0, h].
We also define the operator B : U → X 4 by Bx = (B 1 x, B 2 x, B 3 x, B 4 x) where B i : U → X with i = 1, 2, 3, 4 and B i are given by: It follows from the definition of h (6) that: Consider By the inequality (6), Therefore, from the inequalities ( 7)- (10) and the definition of h as in the inequality (6): Then, this means that B maps U into itself.Next, we want to show that B is a contraction mapping.Consider where Consider where Consider where Consider where We have that from the inequalities ( 11)-( 14) and the definition of h as in the inequality (6): with 0 < L 6 + L 7 + L 8 + L 9 < 1.This shows that B is a contraction mapping.Therefore, by Banach fixed point theorem, system (5) has a unique continuous solution on [0, h].

Continuation Theorems
In this section, we extend the continuation theorem for the system of FAFDEs.Firstly, we give the following lemma.[37]) Suppose γ > 0, a(t) is a nonnegative function that is locally integrable on 0 ≤ t < T ≤ ∞, and g(t) is a nonnegative, nondecreasing continuous bounded function defined on 0 ≤ t < T. If u(t) is nonnegative and locally integrable on 0 ≤ t < T with:
Proof.Suppose that solution x(t) of system (5) has the maximal existence interval defined by [0, β) with β < ∞.Let us consider: for t ∈ [0, β), By the assumption of theorem, It is well-known that the E α is an analytic function [38] and then we can conclude the x 4 is bounded on [0, β).
Next, we consider x 3 on [0, β).Then, for t ∈ [0, β), By Lemma 2, we obtain: Thus, x 3 is bounded on [0, β).Hence, x(t) is bounded on [0, β).Since x are continuous and bounded on [0, β),lim Similar to Theorem 1, we can construct the positive constant h depending on M and h < 1 such that our problem has a solution on [β, β + h].This part is shown in the below.
We define the operator B : U → X 4 by By = ( B 1 y, B 2 y, B 3 y, B 4 y) for y ∈ U and for any y ∈ U, B i is constructed by: To make it clear about the construction of h, without loss of generality, we only show that B 1 maps U into U.Let us consider that for y ∈ U: In this case, the constant h is determined by: We can construct h in this way and follow method as in Theorem 1.
By considering as the Theorem 1, B has a fixed point x, that is x(t) = B x(t) for t ∈ [β, β + h].We next define the function x * on [0, β + h] by: We can see that x * is a continuous solution of our problem on [0, β + h] with β + h > β.
Remark 1.Because x 1 and x 2 are the market price and the fraction of total asset invested, respectively.By assumptions of the model [13], it is possible to assume that x 1 and x 2 are bounded.

Stability Analysis
In order to determine the stability of system (5), we firstly consider the equilibrium point.To find the equilibrium point , we set all equations of system (5) to zero : Then: where: where: where: x * 3,1 = 0, )) for j = 1, 2, 3.It is easy to see that the Jacobian matrix of system ( 5) is: where: , for i = 1, 2, 3, 4 and j = 1, 2, 3.
The characteristic equation of J(C * j ) for j = 1, 2, 3 is: where:   The results in Figures 1-3 indicate that the behavior for the market price of an asset, the fraction of total asset invested in the stock, and the trend-based component of the investor preference have the same pattern at each time for an integer-order system with α = 1.0 and a fractional-order system with α = 0.9 and α = 0.8.To investigate the effect of fractional derivative order on market behavior, the numerical results demonstrate that a decrease in the derivative order is associated with a decrease in the maximum value of x 1 , x 2 and x 3 .Figure 4 shows that decreasing of derivative order leads to a decrease in the value-based component at each time.In addition, the market price as shown in Figure 1 represents the pattern of head and shoulders or rounding top, in which a major high appears by a local maximum.In order to execute planning in the future, investors can buy or sell when the pattern completes.By Theorem 3, the equilibrium point C * = (x * 1 , x * 2 , x * 3 , x * 4 ) is locally stable, where x * 1 = 6.6141, x * 2 = 0.5060, x * 3 = 0 and x * 4 = 0.0121.In other words, the values of x 1 , x 2 , x 3 and x 4 converge to equilibrium point C * at some long-term behaviors.

Discussion and Conclusions
In this paper, we have proposed a fractional order asset flow differential equations model, which is a generalization of an integer order model proposed by Gaginalp [13].Then, we propose and prove the local existence and uniqueness theorems for the solution of the model by applying the Banach fixed point theorem.After that, we extend the local solution, t ∈ (0, h], to a global solution, t ∈ (0, ∞), by proposing the continuation theorem for the model.The stability of the model is also investigated and we obtain the sufficient condition of the parameters.The numerical solutions of the model are illustrated by using an Adams-Bashforth type predictor-corrector method.