Study on the Business Cycle Model with Fractional-Order Time Delay under Random Excitation

Time delay of economic policy and memory property in a real economy system is omnipresent and inevitable. In this paper, a business cycle model with fractional-order time delay which describes the delay and memory property of economic control is investigated. Stochastic averaging method is applied to obtain the approximate analytical solution. Numerical simulations are done to verify the method. The effects of the fractional order, time delay, economic control and random excitation on the amplitude of the economy system are investigated. The results show that time delay, fractional order and intensity of random excitation can all magnify the amplitude and increase the volatility of the economy system.


Introduction
The business cycle model is one of the most attractive models in applying the concept of nonlinear dynamic to investigate economic phenomenon [1].Many scholars have been working on this model.Li et al. [2] studied chaos prediction and chaos control of the model under random excitation.Li et al. [3] investigated the first-passage failure of a business cycle model under wide-brand random excitation.Time delay in a real economy system is ubiquitous.When a policy decision is made, requires some time to implement the economic policy.Yoshida and Asada [4] investigated the effect of policy lag on macroeconomic stability and revealed that policy lag contributed to chaotic motion.Ma et al. [5] studied the stability of the equilibrium and the existence of Hopf bifurcations of business cycle model with discrete delay.Wu et al. [6] investigated the multi-parameter bifurcations of the Kaldor-Kalecki model of business cycles with delay.Khalid Hattaf et al. [7,8] investigated a delayed business cycle model with general investment function.Liu et al. [9] investigated the stability and Hopf bifurcation for a business cycle model with expectation and time delay.Li et al. [10] investigated the stochastic response of business cycle model with time-delay feedback.
When analyzing the time series, the series often present autocorrelation.How to describe this property in a continuous model is undefined.In recent years, a lot of work has been done on the issue of fractional calculus [11][12][13][14][15][16][17][18][19][20][21].In the research of fractional-order delay, Shen et al. [22,23] investigated the dynamical response of Mathieu-Duffing oscillator with fractional-order delayed feedback and analyzed the Duffing oscillator with time-delayed fractional-order Proportion, Integration, Differentiation (PID) controller.In the field of economic and finance modeling, many papers on fractional modeling are published.Nick Laskin [24] investigated the equation of financial assets with fractional derivative and

The Model
In the literature [1], Goodwin considered the lag between decisions to invest and the corresponding outlays that tend to lag behind decisions.Therefore, we may say where Y is the national income, α the marginal propensity to consume, η the investment delay, δ the construction time of new equipment, φ( . x(t)) the function of decisions to invest, O * (t) the autonomous outlays.It is actually the equation that income equals consumption and investment.
By Taylor expansion and dropping all but the first two terms in each, we can obtain δη ..
To eliminate O * (t), we set Equation ( 2) can be written as δη ..
With some transformations, Equation (3) can be written as .. where Entropy 2017, 19, 354 3 of 14 the model can be written as .. where Considering economic control to control the economic fluctuation, Equation ( 7) can be written as ..
where c .
x(t) denotes economic control.Also, considering time delay and memory property of implementing economic policy and random excitation, Equation ( 8) can be written as ..

Stationary PDF of the Model with Random Excitation
In this section, we will apply the stochastic averaging method to obtain the approximate solution of the amplitude of the economy system (9).If the intensity of the noise D = 0, the response of system ( 9) is a stochastic process which is non-Markovian due to the fractional-order time delay.In this section, approximate method has been applied to obtain the stationary PDF for the amplitude of the response.

Stochastic Averaging Method
Equation ( 9) can be regarded as random spread of periodic motion of the linear differential equation. ..
By differentiating Equation (11b) to t, we can obtain ..
So, we can obtain the following equation Entropy 2017, 19, 354 Applying stochastic averaging method to Equation (15) in time interval [0 T], we can obtain the approximation of the amplitude a and the phase θ as ψ(−aω sin ϕ) is periodic function and we set T = 2π.D α (a cos(ϕ − ωτ)) is aperiodic and we set T = ∞, and we can obtain the first part of Equation (16) . .
So, the second part of Equation ( 16) can be written as . .
In Equation ( 17), we applied Caputo's definition and the expression is as follows where α is the fractional-order derivative .
With the transformation s = t − u and ds = −du, Equation (20) can be written .
With two basic formulas, we can obtain .
Similarly, we can also obtain .
Thus, combining with Equation ( 17), we can obtain Entropy 2017, 19, 354 where ξ 1 (t), ξ 2 (t) are independent normalized sources of Gaussian white noise.Then the averaged Itô equation for a(t) is the following form: where So, we can obtain the Fokker-Planck-Kolmogrov (FPK) equation of amplitude a ∂p(a, t) ∂t and the stationary solution of system ( 9) is where C 0 is the normalization constant and C 0 = σ(a) da))da .

Nonlinear Investment Function
In this section, we set ψ( . x(t).v is the reciprocal of an adjustment coefficient.This means that the induced investment is a cubic nonlinear function.Substituting this formula into Equation ( 9), Equation ( 9) can be written as follows ..
Based on the analysis in Section 3, we can obtain the averaged Itô equation for the amplitude The stationary solution of the system (30) is where C 0 is the normalization constant.

The Effect of Fractional-Order Time Delay
Firstly, we investigate the effect of time delay on the amplitude of the economy system.We set τ = 0 and D = 0, then system (29) is a deterministic system.From Figure 1, we can obtain that the fractional time delay can change the stability of the fixed point.When τ = 0 in Figure 1a,c, the fixed point (0, 0) is asymptotical stability.When τ = 1 in Figure 1b,d, a limit cycle can be obtained.In other words, when the effect of fractional time delay exists, the economy system will still have periodic fluctuations even it is under economic control.It means that the economic control is ineffective.So, the policy maker must consider this when making economic control.When 0 D ≠ , we consider the stochastic situation.Figure 2 is the stationary probability density of amplitude ( ) p a for different values of time delay.We can obtain the numerical solutions and the analytical solutions are in good agreement.Thus, it can prove the validity of the method proposed in this paper.The peak of stationary probability density function represents the stability of the economy system.The higher peak indicates the stronger stability.In Figure 2, the peak is decreasing with the increase of time delay.This means that time delay lowers the stability of the economy system.With the increase of time delay the amplitude on its peak of stationary probability density function of the economy system is magnified.This means that the time delay can magnify the amplitude of the economy system.By Equation (31), we obtain the stationary probability density function related to time delay τ .Thus, the policy maker can evaluate the effect of time delay of implementing economic policy to put the economy under control.When D = 0 we consider the stochastic situation.Figure 2 is the stationary probability density of amplitude p(a) for different values of time delay.We can obtain the numerical solutions and the analytical solutions are in good agreement.Thus, it can prove the validity of the method proposed in this paper.The peak of stationary probability density function represents the stability of the economy system.The higher peak indicates the stronger stability.In Figure 2, the peak is decreasing with the increase of time delay.This means that time delay lowers the stability of the economy system.With the increase of time delay the amplitude on its peak of stationary probability density function of the economy system is magnified.This means that the time delay can magnify the amplitude of the economy system.By Equation (31), we obtain the stationary probability density function related to time delay τ.Thus, the policy maker can evaluate the effect of time delay of implementing economic policy to put the economy under control.When 0 D ≠ , we consider the stochastic situation.Figure 2 is the stationary probability density of amplitude ( ) p a for different values of time delay.We can obtain the numerical solutions and the analytical solutions are in good agreement.Thus, it can prove the validity of the method proposed in this paper.The peak of stationary probability density function represents the stability of the economy system.The higher peak indicates the stronger stability.In Figure 2, the peak is decreasing with the increase of time delay.This means that time delay lowers the stability of the economy system.With the increase of time delay the amplitude on its peak of stationary probability density function of the economy system is magnified.This means that the time delay can magnify the amplitude of the economy system.By Equation (31), we obtain the stationary probability density function related to time delay τ .Thus, the policy maker can evaluate the effect of time delay of implementing economic policy to put the economy under control.

The Effect of the Fractional-Order
In this section, we study the effect of fractional order on the amplitude of the economy system (28).Figure 3 is the stationary probability density of amplitude for different values of fractional order.From Figure 3, we can obtain that the peak is decreasing with the decrease of the fractional-order derivative.In other words, with the increase of the fractional-order derivative, the stability of the economy system is magnified.The smaller fractional-order derivative represents the stronger memory effect of economic control.This means that the memory property can lower the stability of the economy system.On the other hand, the smaller fractional-order derivative has a bigger amplitude on its peak of stationary probability density function.This means that the memory property of economic control can magnify the amplitude of the economy system.This indicates that the memory property of economic control actually plays a negative role in making the economy system under control, and we can compute the effect of memory property by Equation (31).So the policy maker can make effective economic policy to make the economy under control.

The Effect of the Fractional-Order
In this section, we study the effect of fractional order on the amplitude of the economy system (28).Figure 3 is the stationary probability density of amplitude for different values of fractional order.From Figure 3, we can obtain that the peak is decreasing with the decrease of the fractional-order derivative.In other words, with the increase of the fractional-order derivative, the stability of the economy system is magnified.The smaller fractional-order derivative represents the stronger memory effect of economic control.This means that the memory property can lower the stability of the economy system.On the other hand, the smaller fractional-order derivative has a bigger amplitude on its peak of stationary probability density function.This means that the memory property of economic control can magnify the amplitude of the economy system.This indicates that the memory property of economic control actually plays a negative role in making the economy system under control, and we can compute the effect of memory property by Equation (31).So the policy maker can make effective economic policy to make the economy under control.

The Effect of Economic Control
Figure 4 is the stationary probability density of amplitude for different values of intensity of economic control.The peak is decreasing with the decrease of the intensity of economic control.This means that the economic control improves the of the economy system.With the decrease of economic control, the amplitude on its peak of stationary probability density function of the economy system is magnified.This means that the economic control lowers the amplitude of the economy system, and we can compute the effect of economic control by Equation (31).It is helpful to the policy maker to make economic policy to make the economy under control.

The Effect of Economic Control
Figure 4 is the stationary probability density of amplitude for different values of intensity of economic control.The peak is decreasing with the decrease of the intensity of economic control.This means that the economic control improves the stability of the economy system.With the decrease of economic control, the amplitude on its peak of stationary probability density function of the economy system is magnified.This means that the economic control lowers the amplitude of the economy system, and we can compute the effect of economic control by Equation (31).It is helpful to the policy maker to make economic policy to make the economy under control.

The Effect of Random Excitation
Figure 5 is the stationary probability density of amplitude ( ) p a for different values of intensity of random excitation.The peak is decreasing with the increase of the intensity of random excitation.This means that the external random excitation lowers the stability of the economy system.With the increase of the intensity of random excitation, the amplitude on its peak of stationary probability density function of the economy system is magnified.This means that external random excitation can magnify the amplitude of the economy system, and we can compute the effect of random excitation by Equation (31).It is helpful to the policy maker to make economic policy to make the economy under control. (a)

The Effect of Random Excitation
Figure 5 is the stationary probability density of amplitude p(a) for different values of intensity of random excitation.The peak is decreasing with the increase of the intensity of random excitation.This means that the external random excitation lowers the stability of the economy system.With the increase of the intensity of random excitation, the amplitude on its peak of stationary probability density function of the economy system is magnified.This means that external random excitation can magnify the amplitude of the economy system, and we can compute the effect of random excitation by Equation (31).It is helpful to the policy maker to make economic policy to make the economy under control.

The Effect of Random Excitation
Figure 5 is the stationary probability density of amplitude ( ) p a for different values of intensity of random excitation.The peak is decreasing with the increase of the intensity of random excitation.This means that the external random excitation lowers the stability of the economy system.With the increase of the intensity of random excitation, the amplitude on its peak of stationary probability density function of the economy system is magnified.This means that external random excitation can magnify the amplitude of the economy system, and we can compute the effect of random excitation by Equation (31).It is helpful to the policy maker to make economic policy to make the economy under control. (a)

Conclusions
In this paper, we investigate the stochastic response of a business cycle model with fractional-order time delay under random excitation.We applied stochastic averaging method to obtain the approximate analytical solution.Thus, we can get the accurate solution of the stationary probability density function.The effect of fractional time delay on the stability and amplitude of the economy system is investigated.The results show that the time delay can induce the Hopf bifurcation in a deterministic model and lower the stability of the economy system in the stochastic situation, and the amplitude of the economy system is magnified by time delay.Also, the effect of fractional order on the stationary probability density function is studied.The results show that the memory property of economic control lowers the stability and magnifies the amplitude of the economy system.Then, the effect of intensity of economic control is investigated.The results show that strengthened economic control can improve the stability and lower the amplitude of the economy system.Finally, the influence of external random excitation on the stationary probability density function is researched.The results indicate that the stronger external random excitation

Conclusions
In this paper, we investigate the stochastic response of a business cycle model with fractional-order time delay under random excitation.We applied stochastic averaging method to obtain the approximate analytical solution.Thus, we can get the accurate solution of the stationary probability density function.The effect of fractional time delay on the stability and amplitude of the economy system is investigated.The results show that the time delay can induce the Hopf bifurcation in a deterministic model and lower the stability of the economy system in the stochastic situation, and the amplitude of the economy system is magnified by time delay.Also, the effect of fractional order on the stationary probability density function is studied.The results show that the memory property of economic control lowers the stability and magnifies the amplitude of the economy system.Then, the effect of intensity of economic control is investigated.The results show that strengthened economic control can improve the stability and lower the amplitude of the economy system.Finally, the influence of external random excitation on the stationary probability density function is researched.The results indicate that the stronger external random excitation brings a bigger amplitude and lower stability of the economy system.The analysis above is helpful to the policy maker to make economic policy to make the economy under control.

Figure 3 .
Figure 3.The solid line "-" is analytical solutions, the crosses "x" are numerical solutions;

Figure 4 .
Figure 4.The solid line "-" is analytical solutions, the crosses "x" are numerical solutions;

Figure 4 .
Figure 4.The solid line "-" is analytical solutions, the crosses "x" are numerical solutions;