Study on the Stability and Entropy Complexity of an Energy-Saving and Emission-Reduction Model with Two Delays

In this paper, we build a model of energy-savings and emission-reductions with two delays. In this model, it is assumed that the interaction between energy-savings and emission-reduction and that between carbon emissions and economic growth are delayed. We examine the local stability and the existence of a Hopf bifurcation at the equilibrium point of the system. By employing System Complexity Theory, we also analyze the impact of delays and the feedback control on stability and entropy of the system are analyzed from two aspects: single delay and double delays. In numerical simulation section, we test the theoretical analysis by using means bifurcation diagram, the largest Lyapunov exponent diagrams, attractor, time-domain plot, Poincare section plot, power spectrum, entropy diagram, 3-D surface chart and 4-D graph, the simulation results demonstrating that the inappropriate changes of delays and the feedback control will result in instability and fluctuation of carbon emissions. Finally, the bifurcation control is achieved by using the method of variable feedback control. Hence, we conclude that the greater the value of the control parameter, the better the effect of the bifurcation control. The results will provide for the development of energy-saving and emission-reduction policies.


Introduction
Energy resources are the backbone of any national economy.Scholars worldwide have been studying energy prices, the low-carbon economy and energy-savings and emission-reduction.Lv and Zhou [1] studied the dynamic behavior of the energy price model with time delay.They explained the reasons why the energy price model was developed and maintained periodically by the bifurcation theory.Tian et al. [2] developed a nonlinear model for the development of oil, gas and other energy resources against the backdrop of the energy structure in Jiangsu, China, one that predominantly relies on the coal consumption, and evidenced empirically the feasibility of the corresponding energy countermeasures.
Despite the fact that energy promotes economic development and improves people's living standards, excessive energy consumption undoubtedly damages the environment, so it is necessary to adopt a low-carbon economy and reduce the possible environmental damage.In fact, the energy-saving and emission-reduction system is a complex nonlinear system, which involves the interactions between a variety of factors such as energy-saving and emission-reduction, carbon emissions, economic growth, energy efficiency, carbon tax and energy intensity, and so forth [3][4][5].In view of this, a novel three-dimensional energy-saving and emission-reduction chaotic system is proposed [6].
Entropy 2016, 18, 371 2 of 20 The system is established in accordance with the complicated relationship among energy-savings and emission-reduction, carbon emissions and economic growth.This system displays a very complex phenomenon by including a special chaotic attractor named the energy-saving and emission-reduction attractor (the EE attractor for short), which is different from the previous chaotic attractor, such as Lorenz attractor [7], Chen attractor [8], Lü attractor [9], Energy resource attractor [10][11][12] and so on.Moreover, drawing up rational and effective policies and laws will ensure a better progress in energy-saving and emission-reduction [13,14].Wang and Xu [15] reported a new four-dimensional energy-saving and emission-reduction chaotic system.The system is obtained in accordance with the complicated relationship among energy-saving and emission-reduction, carbon emission, economic growth and new energy development.The model is described as follows: where x(t) is the time-dependent variable of energy-saving and emission-reduction; y(t) the time-dependent variable of carbon emissions; z(t) the time-dependent variable of economic growth (GDP) and u(t) the time-dependent variable of new energy development.a i , d i , b j , c j , (i = 1, 2, 3; j = 1, 2, 3, 4) are coefficients and M, N, L, K, C, E positive constants [15].Model (1) stands for the development and utilization of new energy sources.Therefore, the four-dimensional model of energy-saving and emission-reduction is a step closer to the actual conditions.In view of the current mode of economic development and technological conditions, economic development will increase the consumption of coal and oil, and in turn this will cause an increase in carbon emissions.However, there are many other approaches to promote economic growth, such as those in the modern service sector that are low energy-consuming but high value-added, so economic development may not immediately cause a substantial energy consumption load.In other words, it might not lead to an obvious increase in carbon emissions in the short term.In this essay, we have included the delay between economic development and the significant growth of carbon emission in this kind of economic model.
Moreover, energy-savings and emission-reduction reduce carbon emissions by saving or reducing energy consumption.However, due to the complexity of the economic system, the implementation of energy-saving and emission-reduction measures normally will not lead to a significant reduction, showing a certain lag.Meanwhile, energy-saving and emission-reduction will slow down the pace of economic development.This will force enterprises to change their mode of economic development, strengthen energy conservation and efficient use, and actively comply with a recycling economy and low carbon economy, so energy-savings and emission-reduction will promote the sound and fast development of economy, but it still demonstrates an obvious delay.Therefore, the impact of energy-saving and emission-reduction on carbon emission and economic development is estimated as delayed.On the basis of Model (1), we propose a two-delay model as follows: where τ 1 represents the delay time between economic development and carbon emission, and τ 2 represents the delay time among energy-saving and emission-reduction, economic development and carbon emission.The rest of this paper is organized as follows: in Section 2, we focus on the local stability and the existence of Hopf bifurcation at the equilibrium point.In Section 3, we study the effects of delays and Entropy 2016, 18, 371 3 of 20 the feedback control on the stability of the system.In Section 4, the effective Hopf bifurcation control of the system is examined.Finally, conclusions are drawn in Section 5.

Equilibrium Points and Local Stability
The calculation formulas of equilibrium points of Equation ( 2) are fairly complicated, so the equilibrium points can be calculated directly by using specific parameters value in numerical simulation section.We assume that Equation (2) has an equilibrium point, called hereafter by E(x * , y * , z * , u * ).It means that x(t), y(t), z(t) and u(t) can be in an equilibrium state through dynamic game.In this paper, we focus on the influence of τ 1 , τ 2 and the feedback control on the stability of system (2) at the equilibrium point.
Proposition 2. Suppose that ∆ = 0, then Re dλ(τ 20 ) and ∆ have the same sign, where: Proof.Substituting λ(τ 2 ) (the root of the characteristic Equation ( 4), which is defined above) into Equation ( 4) and taking the derivative with respect to τ 2 , we obtain: From Equation ( 18), we have: where: and ∆ have the same sign.
Entropy 2016, 18, 371 8 of 20 Thus, according to propositions 1 and 2 and the Hopf bifurcation theorem in [17], we have the following theorem.

Numerical Simulation and Analysis
In this section, the numerical simulation is carried out in order to support the theoretical analysis.Let x(0) = 0.3; y(0) = 0.5; z(0) = 0.7; u(0) = 0.9; We consider the following system with given parameter values: The following equilibrium points can be obtained by simple calculation: According to actual economic meaning, only E 1 is Nash Equilibrium point.So we examine the system stability of the system at the equilibrium point E 1 (−0.6625,1.1410, −0.0578, 0.2100).

The Influence of τ 1 on the Stability of Equation (19)
Equation ( 19) moves from stable to unstable with the increase in τ 1 and undergoes Hopf bifurcation at τ 1 = 0.4618 when τ 2 = 0. Figure 1a shows the dynamic evolution process of Equation (19).In this paper, we calculate the largest Lyapunov exponent (LLE) by the method of Wolf reconstruction [19].We judge whether the system is stable according to the exponent value.If it is less than 0, the system is stable.If it is greater than 0, the system is unstable.If it is equal to 0, the system appears bifurcation.In Figure 1b, we can know that the system has a bifurcation at τ 1 = 0.4618.This is consistent with the conclusion of theoretical analysis.Equation (19).In this paper, we calculate the largest Lyapunov exponent (LLE) by the method of Wolf reconstruction [19].We judge whether the system is stable according to the exponent value.If it is less than 0, the system is stable.If it is greater than 0, the system is unstable.If it is equal to 0, the system appears bifurcation.In Figure 1b, we can know that the system has a bifurcation at 1 0.4618 τ = .
This is consistent with the conclusion of theoretical analysis.According to Theorem 1, Equation ( 19) is stable when 1 and 2 0 τ = .In this case, ( )  x t , ( ) y t , ( ) z t and ( ) u t will converge to the equilibrium point 1 E through the game.
Figures 2 and 3 show the above properties.According to Theorem 1, Equation ( 19) is stable when τ 1 = 0.4 < τ 10 = 0.4618 and τ 2 = 0.In this case, x(t), y(t), z(t) and u(t) will converge to the equilibrium point E 1 through the game.Figures 2 and 3 show the above properties.
Equation ( 19) is unstable when τ 1 = 0.5 > τ 10 = 0.4618 and τ 2 = 0 by Theorem 1.At this point, the frequency spectrum plot is discrete, which means that the system has a periodic solution.The x(t), y(t), z(t) and u(t) will lie on the basin of attraction [20] through the game.These analyses can be illustrated by Figures 4 and 5.
We discover that economic growth (slowdown) will lead to an increase (reduction) in carbon emissions when other parameters are fixed.The two variables share a consistent pattern, but with a certain delay in time.Based on the above simulation, we find out that when the delay is greater than the bifurcation value, the carbon emissions pattern will not be consistent with the economic growth trend and volatility surfaces.Therefore, we should guarantee the timing of achieving the goal of reducing emissions and make it less than the bifurcation value.Only in such a stable environment, the implementation of energy-saving and emission-reduction policies will intend a favorable effect.Equation (19).In this paper, we calculate the largest Lyapunov exponent (LLE) by the method of Wolf reconstruction [19].We judge whether the system is stable according to the exponent value.If it is less than 0, the system is stable.If it is greater than 0, the system is unstable.If it is equal to 0, the system appears bifurcation.In Figure 1b, we can know that the system has a bifurcation at 1 0.4618 τ = .
This is consistent with the conclusion of theoretical analysis.According to Theorem 1, Equation ( 19) is stable when 1 1 0 0.4 0.4618 τ τ = < = and 2 0 τ = .In this case, ( )  x t , ( ) y t , ( ) z t and ( ) u t will converge to the equilibrium point 1 E through the game.
Figures 2 and 3 show the above properties.the frequency spectrum plot is discrete, which means that the system has a periodic solution.The ( ) x t , ( ) y t , ( ) z t and ( ) u t will lie on the basin of attraction [20] through the game.These analyses can be illustrated by Figures 4 and 5.
(a) (b)  Equation ( 19) is unstable when 1 1 0 0.5 0.4618 and 2 0 τ = by Theorem 1.At this point, the frequency spectrum plot is discrete, which means that the system has a periodic solution.The ( ) x t , ( ) y t , ( ) z t and ( ) u t will lie on the basin of attraction [20] through the game.These analyses can be illustrated by Figures 4 and 5.
(a) (b)    We discover that economic growth (slowdown) will lead to an increase (reduction) in carbon emissions when other parameters are fixed.The two variables share a consistent pattern, but with a certain delay in time.Based on the above simulation, we find out that when the delay is greater than the bifurcation value, the carbon emissions pattern will not be consistent with the economic growth trend and volatility surfaces.Therefore, we should guarantee the timing of achieving the goal of reducing emissions and make it less than the bifurcation value.Only in such a stable environment, the implementation of energy-saving and emission-reduction policies will intend a favorable effect.

The Influence of 1 τ on the Entropy of Equation (19)
Kolmogorov entropy can be used to measure the degree of complexity of the system.Let k be the value of Kolmogorov entropy.If the system is stable and in regular motion, then 0 k = , otherwise 0 k > .The greater the value of k is, the more complex the system is.We have already found that Equation ( 19) displays bifurcation at 1 0.4618 τ = , as is shown in Figure 1.Based on the above analysis, we can get that if 1 0.4618 τ < , then 0 k = , otherwise 0 k > .The more complex the system is, the longer and the more difficult it will be for the system to return to stability.The entropy properties are displayed in Figure 6.

The Influence of τ 1 on the Entropy of Equation (19)
Kolmogorov entropy can be used to measure the degree of complexity of the system.Let k be the value of Kolmogorov entropy.If the system is stable and in regular motion, then k = 0, otherwise k > 0. The greater the value of k is, the more complex the system is.We have already found that Equation ( 19) displays bifurcation at τ 1 = 0.4618, as is shown in Figure 1.Based on the above analysis, we can get that if τ 1 < 0.4618, then k = 0, otherwise k > 0. The more complex the system is, the longer and the more difficult it will be for the system to return to stability.The entropy properties are displayed in Figure 6.We discover that economic growth (slowdown) will lead to an increase (reduction) in carbon emissions when other parameters are fixed.The two variables share a consistent pattern, but with a certain delay in time.Based on the above simulation, we find out that when the delay is greater than the bifurcation value, the carbon emissions pattern will not be consistent with the economic growth trend and volatility surfaces.Therefore, we should guarantee the timing of achieving the goal of reducing emissions and make it less than the bifurcation value.Only in such a stable environment, the implementation of energy-saving and emission-reduction policies will intend a favorable effect.

The Influence of 1 τ on the Entropy of Equation (19)
Kolmogorov entropy can be used to measure the degree of complexity of the system.Let k be the value of Kolmogorov entropy.If the system is stable and in regular motion, then 0 k = , otherwise 0 k > .The greater the value of k is, the more complex the system is.We have already found that Equation ( 19) displays bifurcation at 1 0.4618 τ = , as is shown in Figure 1.Based on the above analysis, we can get that if 1 0.4618 τ < , then 0 k = , otherwise 0 k > .The more complex the system is, the longer and the more difficult it will be for the system to return to stability.The entropy properties are displayed in Figure 6.The influence of τ 1 and b 3 on y is shown in Figure 7.The system shows a sign of instability at τ 1 = 0.4618.However, b 3 has no obvious impact on y when τ 1 < 0.55.With the increase of b 3 , y becomes larger and then decreases when τ 1 > 0.55.The minimum value of y is −1.448 for (τ 1 , b 3 ) = (0.5, 0.85) while the maximum value of y is 5.124 for (τ 1 , b 3 ) = (0.6, 0.55).The value of y is −0.66 approximately when τ 1 < 0.4618.Thus, the impact of b 3 on y can be ignored when τ 1 < 0.4618.We focus on the influence of .We conclude that the stability of energy-saving and emission-reduction system can be maintained, only controlling for the value of delay parameter 1 τ .We focus on the influence of .We conclude that the stability of energy-saving and emission-reduction system can be maintained, only controlling for the value of delay parameter 1 τ .In this part, we study the effect of 2 τ on the stability of the system when 1 τ is fixed.With an increase in 2 τ , Equation ( 19) will lose stability.Figure 10 displays that Equation ( 19) undergoes bifurcation at 2 0.0622 . This is consistent with the theoretical analysis.Equation ( 19) is stable when 1    In this part, we study the effect of τ 2 on the stability of the system when τ 1 is fixed.With an increase in τ 2 , Equation ( 19) will lose stability.Figure 10 displays that Equation ( 19) undergoes bifurcation at τ 2 = 0.0622 for τ 1 = 0.4.This is consistent with the theoretical analysis.In this part, we study the effect of 2 τ on the stability of the system when 1 τ is fixed.With an increase in 2 τ , Equation ( 19) will lose stability.Figure 10 displays that Equation ( 19) undergoes bifurcation at 2 0.0622 . This is consistent with the theoretical analysis.Equation ( 19) is stable when 1    Equation ( 19) is stable when τ 1 = 0.4, τ 2 = 0.04 < τ 20 = 0.0622 by Theorem 2. x(t), y(t), z(t) and u(t) will converge to the equilibrium point E 1 through the game.The properties are shown in Figures 11 and 12.In this part, we study the effect of 2 τ on the stability of the system when 1 τ is fixed.With an increase in 2 τ , Equation ( 19) will lose stability.Figure 10 displays that Equation ( 19) undergoes bifurcation at 2 0.0622 . This is consistent with the theoretical analysis.Equation ( 19) is stable when 1     .
Meanwhile, the Poincare plot has five discrete points, which means that the system has a periodic solution.( ) x t , ( ) y t , ( ) z t and ( ) u t will lie on the basin of attraction through the game.Figures   On the basis of theorem 2, Equation ( 19) is unstable when τ 1 = 0.4, τ 2 = 0.08 > τ 20 = 0.0622.Meanwhile, the Poincare plot has five discrete points, which means that the system has a periodic solution.x(t), y(t), z(t) and u(t) will lie on the basin of attraction through the game.Figures 13 and 14 show these characteristics.  .
Meanwhile, the Poincare plot has five discrete points, which means that the system has a periodic solution.( ) x t , ( ) y t , ( ) z t and ( ) u t will lie on the basin of attraction through the game.Figures     From Figure 10, we know that Equation ( 19) has a bifurcation when 1 2 0.4 , 0.0622 Equation ( 19) is unstable for 2 0.0622 τ > , then the entropy value is more than 0. As in Section 3.2, the growth in the value of entropy shares the same pattern with that in 2 τ .The change of entropy is shown in Figure 15.With the increase of entropy value, the system will be more unstable.In this case, the implementation effect of energy-saving and emission-reduction is poorly maintained.From Figure 10, we know that Equation ( 19) has a bifurcation when τ 1 = 0.4, τ 2 = 0.0622.Equation ( 19) is unstable for τ 2 > 0.0622, then the entropy value is more than 0. As in Section 3.2, the growth in the value of entropy shares the same pattern with that in τ 2 .The change of entropy is shown in Figure 15.With the increase of entropy value, the system will be more unstable.In this case, the implementation effect of energy-saving and emission-reduction is poorly maintained.From Figure 10, we know that Equation ( 19) has a bifurcation when 1 2 0.4 , 0.0622 Equation ( 19) is unstable for 2 0.0622 τ > , then the entropy value is more than 0. As in Section 3.2, the growth in the value of entropy shares the same pattern with that in 2 τ .The change of entropy is shown in Figure 15.With the increase of entropy value, the system will be more unstable.In this case, the implementation effect of energy-saving and emission-reduction is poorly maintained.Figure 16 shows that the influence of τ 2 on y is greater.With an increase in τ 2 , the system moves from stable to unstable.The system is becoming unstable at τ 2 = 0.06.The minimum value of y is 1.048 for (τ 2 , b 3 ) = (0.1, 0.1) and the maximum value of y is 1.212 for (τ 2 , b 3 ) = (0.09, 0.1).The value of y is stable in the vicinity of 1.1.But b 3 has little effect on y.Only when τ 2 > 0.06, the value of y will turn from large to small with the increase of b 3 .Therefore, we have to make τ 2 < 0.06, so as to guarantee the effective implementation of energy-savings and emission-reduction., so as to guarantee the effective implementation of energy-savings and emission-reduction.It can be seen from Figure18 that the system loses stability if 1 τ shows an increase.In that case, x demonstrates a large range of ups and downs.On the other hand, 2 τ has no impact on x ., so as to guarantee the effective implementation of energy-savings and emission-reduction.It can be seen from Figure18 that the system loses stability if 1 τ shows an increase.In that case, x demonstrates a large range of ups and downs.On the other hand, 2 τ has no impact on x .It can be seen from Figure 18 that the system loses stability if τ 1 shows an increase.In that case, x demonstrates a large range of ups and downs.On the other hand, τ 2 has no impact on x. , τ τ on the Entropy of Equation ( 19) Figure 19 shows that the entropy becomes greater than 0 with an increase in 1 τ .Meanwhile, the system is unstable.According to Section 3.    3.8.The Influence of τ 1 , τ 2 on the Entropy of Equation ( 19) Figure 19 shows that the entropy becomes greater than 0 with an increase in τ 1 .Meanwhile, the system is unstable.According to Section 3.4, (τ 1 , τ 2 ) = (0.4,0.0622) is a bifurcation point, so we can infer that the point (0.4, 0.0622) is on the boundary of the entropy change.τ 2 still has no effect on entropy.So keeping the other parameters constant, the system is stable when τ 1 < 0.4., τ τ on the Entropy of Equation ( 19) Figure 19 shows that the entropy becomes greater than 0 with an increase in 1 τ .Meanwhile, the system is unstable.According to Section 3.

Bifurcation Control
The fact that the system loses stability at the bifurcation will seriously impact the effectiveness of the implementation of the policy.Thus we have to take measures to ensure the stability of the system.The system can be controlled by some approaches, for example, modified straight-line stabilization method [21], pole placement method [22], the OGY method (a control method of chaos proposed by Ott, Grebogi and Yorke) [23], time-delayed feedback method [24], and so forth.Here, we adopt the method of variable feedback control (see, e.g., [25], and references therein) to control for the Hopf bifurcation.The controlled system is given as follows: where k is a feedback control parameter.The state of Equation ( 20) can be changed by adjusting the value of k.If k is greater than the critical value, then Equation ( 20) returns to the stable state, otherwise the system is still in an unstable state.
Figure 20 shows the dynamic evolution process of Equation ( 20) to k for (τ 1 , τ 2 ) = (0.5, 0).We find out that Equation ( 20) undergoes bifurcation at k = 0.1689.In other words, Equation ( 20) is unstable when k < 0.1689 and is stable when k > 0.1689.where k is a feedback control parameter.The state of Equation ( 20) can be changed by adjusting the value of k.If kis greater than the critical value, then Equation ( 20) returns to the stable state, otherwise the system is still in an unstable state.

Bifurcation Value of Equation (20) to k
In case 1, we know that Equation ( 2 Figure 20 shows the dynamic evolution process of Equation ( 20) to k for 1 2 ( , ) (0.5,0) τ τ = . We find out that Equation ( 20      where k is a feedback control parameter.The state of Equation ( 20) can be changed by adjusting the value of k.If kis greater than the critical value, then Equation ( 20) returns to the stable state, otherwise the system is still in an unstable state.

Bifurcation Value of Equation (20) to k
In case 1, we know that Equation ( 2 Figure 20 shows the dynamic evolution process of Equation ( 20) to k for 1 2 ( , ) (0.5,0) τ τ = . We find out that Equation ( 20 20) are displayed in Figure 22.We can see that Equation (20) maintains stable and x, y and z tend to equilibrium point.Therefore, bifurcation control is successful as k > 0.1689.We can conclude that the control effect will be larger if the control parameter k stands higher. .We can conclude that the control effect will be larger if the control parameter k stands higher.Based on the above analysis, we assume that the unstable system can return to stability if effective control is exerted.This means that when problems arise from energy-saving and emission-reduction policies and lead to market volatility, we can take control measures to ensure that the system returns to a stable state.

Conclusions
In this paper, we examine the stability of a four dimensional energy-saving and emission-reduction model with two delays.We focus on the impacts of delays and the feedback control on the stability of the system.The system shifts from stable to unstable when the delays are larger than the bifurcation values.In this case, a large fluctuation shows up in the system and affects the implementation of energy-saving and emission reduction policies.However, we can control the bifurcation of the system by adopting the variable feedback control approach and the system will return to a stable state when the control parameters are set.
The results show that time delays play an important role in the stability of the system.Only when the energy-savings and emission-reduction system is stable, can the purpose of reducing carbon emissions be achieved.Therefore, we must ensure that the delay parameters are in a stable region and then add other parameters to maintain a long-term stable system.
According to the above analysis, if we want to keep the energy-saving and emission-reduction system running in a stable state, we can take the following strategies when other things are equal.First, when 2 0 τ = , we must make sure that 1 0.4618  Based on the above analysis, we assume that the unstable system can return to stability if effective control is exerted.This means that when problems arise from energy-saving and emission-reduction policies and lead to market volatility, we can take control measures to ensure that the system returns to a stable state.

Conclusions
In this paper, we examine the stability of a four dimensional energy-saving and emission-reduction model with two delays.We focus on the impacts of delays and the feedback control on the stability of the system.The system shifts from stable to unstable when the delays are larger than the bifurcation values.In this case, a large fluctuation shows up in the system and affects the implementation of energy-saving and emission reduction policies.However, we can control the bifurcation of the system by adopting the variable feedback control approach and the system will return to a stable state when the control parameters are set.
The results show that time delays play an important role in the stability of the system.Only when the energy-savings and emission-reduction system is stable, can the purpose of reducing carbon emissions be achieved.Therefore, we must ensure that the delay parameters are in a stable region and then add other parameters to maintain a long-term stable system.
According to the above analysis, if we want to keep the energy-saving and emission-reduction system running in a stable state, we can take the following strategies when other things are equal.First, when τ 2 = 0, we must make sure that τ 1 < 0.4618 or b 4 > 0.5.Second, when τ 1 = 0.4, we must keep τ 2 < 0.0622.

Figure 2 .
Figure 2. Time-domain plot when 1 Figure 3.The EE attractor when
Figure 5.The EE attractor when

b
, y becomes larger and then decreases when 1

Figure 8 shows that 4 b 4 (.
Figure 8 shows that 4 b has a greater impact on y .With an increase in 4 b , the system shifts from unstable to stable.When 4 0.2 b =

4 b 4 b
leads only to a larger fluctuation but when 1 0.4618 τ > , the fluctuation trend is basically symmetrical.However, 3 b and have no effect on y when 1 0.4618 τ <

Figure 8 b
Figure8shows that b 4 has a greater impact on y.With an increase in b 4 , the system shifts from unstable to stable.When b 4 = 0.2, the system is beginning to become stable.The minimum value of y is −1.838 for (τ 1 , b 4 ) = (0.5, 0.15) and the maximum value of y is 4.207 for (τ 1 , b 4 ) = (0.6, 0).The value of y is stable at 2.889, but τ 1 has little effect on y.Therefore, we should consider the effect of b 4 on the stability of the system.

Figure 8 shows that 4 b 4 (.
Figure 8 shows that 4 b has a greater impact on y .With an increase in 4 b , the system shifts from unstable to stable.When 4 0.2 b =

4 b
trend is basically symmetrical.However, 3 b and have no effect on y when 1 0.4618 τ <
u t will converge to the equilibrium point 1 E through the game.The properties are shown in Figures11 and 12 .
u t will converge to the equilibrium point 1 E through the game.The properties are shown in Figures11 and 12 .
u t will converge to the equilibrium point 1 E through the game.The properties are shown in Figures11 and 12 .
Figure 12.The EE attractor when
Figure 14.The EE attractor when

3. 5 .
The Influence of τ 2 on the Entropy of Equation(19) Figure 14.The EE attractor when

Figure 16 shows that the influence of 2 τ 2 0
Figure16shows that the influence of 2 τ on y is greater.With an increase in 2 τ , the system

Figure 16 shows that the influence of 2 τ 2 0
Figure16shows that the influence of 2 τ on y is greater.With an increase in 2 τ , the system
can be seen from Figures4 and 5.

4. 2 .
Equation (20) is Unstable When 0.1689 k < Keep the values of other parameters unchanged and let 0.05 0.1689 k = <, and the time-domain plot and the EE attractor of Equation (20) are shown in Figure21.We find out that the system (20) is unstable and has a basin of attraction.It fails to achieve bifurcation control.

4. 3 .
Equation (20) is Stable When 0.1689 k > When the other parameters are consistent with the original, let 0.2 0.1689 k = > and the time-domain plot and the EE attractor of Equation (20) are displayed in Figure22.We can see that Equation (20) maintains stable and x , y and z tend to equilibrium point.Therefore, bifurcation

4. 2 .
Equation (20) is Unstable When k < 0.1689 Keep the values of other parameters unchanged and let k = 0.05 < 0.1689, and the time-domain plot and the EE attractor of Equation (20) are shown in Figure 21.We find out that the system (20) is unstable and has a basin of attraction.It fails to achieve bifurcation control.Entropy 2016, 18, 371 19 of 21 can be seen from Figures4 and 5.

4. 2 .
Equation (20) is Unstable When 0.1689 k < Keep the values of other parameters unchanged and let 0.05 0.1689 k = <, and the time-domain plot and the EE attractor of Equation (20) are shown in Figure21.We find out that the system (20) is unstable and has a basin of attraction.It fails to achieve bifurcation control.

4. 3 .
Equation (20) is Stable When 0.1689 k > When the other parameters are consistent with the original, let 0.2 0.1689 k = > and the time-domain plot and the EE attractor of Equation (20) are displayed in Figure22.We can see that Equation (20) maintains stable and x , y and z tend to equilibrium point.Therefore, bifurcation

4. 3 .
Equation (20) is Stable When k > 0.1689 When the other parameters are consistent with the original, let k = 0.2 > 0.1689 and the time-domain plot and the EE attractor of Equation (