Threshold behaviour of a triple-delay SIQR stochastic epidemic model with Lévy noise perturbation

Notations:

The superscript "$T$" refers to the matrix transpose. The matrix $\mathcal{P}>0$ $(\mathcal{P}\geq 0)$ means that $\mathcal{P}$ is positive definite (positive semi-definite). "$\ast$" is the symmetric terms in a matrix. $W[\alpha,\beta)$ is the space of vector functions $\varphi$ defined over $[\alpha,\beta]$ that are absolutely continuous with a finite ${\rm lim}_{s\rightarrow \beta^-}\varphi(s)$ and square integrable derivatives. The norm of $W[\alpha,\beta)$ is defined as

1.   Introduction

Over the past two decades, finance systems have been extensively studied due to their sophisticated dynamical behaviors such as the chaos and the bifurcation ^{[1,2,3,4,5,6]}. Note that the chaotic characteristics of finance systems could induce the potential uncertainties of the macroeconomic operation. Therefore, in the past decade or so, the stabilization and synchronization problems for the chaotic finance systems have become the most concerned research topics ^{[7,8,9,10,11]}. In particular, various control schemes have been adopted to achieve the effective stabilization and synchronization. For instance, in ^{[12,13,14,15]}, the delayed control scheme has been utilized to stabilize the chaotic and hyperchaotic finance systems and, in ^[16,17], the adaptive controllers have been proposed to realize desirable control and synchronization performance for the finance systems. In ^[18,19], the adaptive sliding control strategy has been employed to stabilize the fractional-order finance systems. In ^[20], a resilient guaranteed cost controller has been designed to control the chaotic finance system. Also, the impulsive controller and intermittent controller have been designed in ^[21,22], respectively.

In the literature ^[23], the control scheme based on discrete-time observations (DTOs) has been proposed to stabilize the continuous-time stochastic differential equations with Markov chain. Compared with the continuous-time control, such a control scheme costs less as the system state is only required to be observed at some discrete-time instants. The results in ^[23] have been further improved in ^[24] by using some new techniques and, in ^[25,26], the time delay has been taken into account in the DTOs-based control scheme. Note that the control setting proposed in ^{[19,20,21,22]} is essentially the same as the sampled-data control encountered in engineering control systems ^[27,28]. For the sampled-data control, it is worth pointing out that the discontinuous Lyapunov-Krasovskii (L-K) functionals and the Wirtinger's inequality have been developed in ^[28] to establish more effective stabilization conditions for linear control systems under a sampled-data controller with the transmission delay.

So far, most existing literature with respect to the control of finance systems have been based on continuous-time controllers. In addition, discontinuous control schemes have been proposed in ^[21,22] to address the synchronization problem for chaotic and hyperchaotic finance systems. However, to our knowledge, the DTOs-based control strategy has not been adopted to discuss the finance systems, not mention to the time delay is involved. In fact, the DTOs-based control is more realistic for finance systems since the financial policies are generally implemented for a period of time and then updated on the basis of the current economic situation. In addition, the time delay of policy implementation is often unavoidable. Thus, the time delay should be considered in designing controller. Unfortunately, the existing results concerning the DTOs-based time-delayed feedback can be only applicable for nonlinear control systems subject to the rigorous linear growth conditions ^[25,26]. Note that the sampled-data control with transmission delay in ^[28] is similar to the DTOs-based time-delayed control. However, it is worth mentioning that the results in ^[28] are only concerned with linear systems, which are no longer applicable for chaotic finance systems due to the existence of nonlinear characteristics.

Inspired by the aforementioned discussions, the paper is devoted to considering the stabilization problem for a hyperchaotic finance system via a time-delayed feedback controller based on DTOs. By incorporating the quadratic system theory, a piecewise augmented discontinuous L-K functional, and some advanced inequalities, a local stabilization criterion is first established by means of linear matrix inequalities (LMIs). As the by-product, a simplified criterion is also provided in the case of no time delay. Moreover, the optimization problems are given to derive the larger domain of attraction (DOA). Finally, simulations show the availability of the derived results. The novelties of the paper are given as below:

1) The DTOs-based time-delayed control scheme is proposed, for the first time, to stabilize the hyperchaotic finance system.

2) A piecewise augmented discontinuous L-K functional is constructed under which a novel local stabilization criterion is obtained by means of LMIs.

3) The state evolution over the first time-interval is specifically considered in establishing the local stabilization criterion.

2.   Problem formulation

In ^[1,2], the authors have proposed a finance system containing three variables and nine independent parameters. The finance system is composed of four sub-blocks (namely, labor force, production, stock and money) and can be simplified as follows:

where $z_1(t)$, $z_2(t)$ and $z_3(t)$ represent, respectively, the interest rate, the investment demand and the price index; $a$, $b$ and $c$ denote, respectively, the saving amount, the cost per investment and the demand elasticity of commercial markets.

In the literature ^[9], by introducing an additional variable $z_4(t)$ representing the average profit margin, the system (1) has been modified as below:

where $l$ and $m$ are scalars. In ^[9], it has been shown that the finance system (2) displays the complicated hyperchaotic phenomenon for the case that $a = 0.9$, $b = 0.2$, $c = 1.5$, $l = 0.2$ and $m = 0.17$. Under the assumption that

it is easy to verify that the hyperchaotic system (2) has the following equilibrium points:

Denoting $z(t)\triangleq[z_1(t)\; \; z_2(t)\; \; z_3(t)\; \; z_4(t)]^T$ and adding the feedback control $u(t)$ into (2), we have

where

In order to stabilize the continuous-time stochastic hybrid dynamical systems, in ^[23,24], Mao et al. have proposed the following DTOs-based feedback controller:

where $h > 0$ refers to the duration between two consecutive observations, $[t/h]$ is the maximum integer less than or equal to $t/h$. Considering that the time delay is often inevitable in data transmission, in ^[25,26], the controller (5) has been modified as

where $\eta > 0$ is a time delay.

Remark 1. Under the DTOs-based control scheme, it is seen from (5) that the only the state at the discrete instants $0, h, 2h, \cdots$ are needed in designing the controller. Compared with the feedback control using the continuous-time state, it is clear that the DTOs-based scheme costs less. The DTOs-based control is essentially the same as the sampled-data control in engineering systems ^[27,28].

In the paper, we will design the DTOs-based time-delayed feedback controller. Let $s_k\triangleq kh$ ($k = 0, 1, 2, \cdots$) be the state observation instants and $t_k\triangleq kh+\eta$ ($k = 0, 1, 2, \cdots$) be the updating instants of control signals. Then, the DTOs-based time-delayed controller can be described as

Remark 2. In reality, the financial policies are generally implemented over a period of time and then modified on the basis of the current economic situation. Moreover, the time delay of policy implementation is often unavoidable. Compared with the continuous-time delayed feedback ^{[9,12,13,14,15]}, the DTOs-based delayed feedback might be more realistic in controlling the unstable finance systems. Literature survey shows that, this paper is the first time to address the stabilization problem for the unstable finance systems under a DTOs-based delayed control scheme.

Let $z^\ast = [z_1^\ast\; \; z_2^\ast\; \; z_3^\ast\; \; z_4^\ast]^T$ be the unstable equilibrium point of the system (2). Then, it follows that

Moreover, using (4), (6) and (7), one has the closed-loop system

where $e(t)\triangleq z(t)-z^*$, $\bar{A}\triangleq A+F_0$ and

In particular, it is seen that the nonlinearity $\bar{f}(e)$ can be explicitly formulated as

where $F_3 = 0_{4\times4}$, $F_4 = -dF_1$ and

Using (9), the system (8) can be modified as the following form:

Denoting

it is obvious that $\eta\leq\tau(t)\leq h+\eta$ with $\dot{\tau}(t) = 1$ for $t\neq t_k$. Furthermore, the system (10) can be described as the time-varying delay system

Due to the existence of time delay, the control signals will be updated only when $t\geq t_0 = \eta$. In this case, the system (4) should be specially handled within the time interval $[0, \eta)$ in the framework of local stabilization ^[29]. Here, we set $u(t) = 0$ within $[0, \eta)$. Then, one has the open-loop system

The paper aims to design the DTOs-based time-delayed controller (6) such that the resulting closed-loop system (10) is locally asymptotically stable and has a larger estimate of the DOA.

As in ^[29], it is assumed that the initial conditions of (12) is denoted by $e(t) = e_0$, $t\in[-\eta, 0]$. Note that the finance system (2) is a typical quadratic system ^[30,31,32]. For convenience of the subsequent analysis, as in ^[15,30], the following box is introduced:

where $\varepsilon_j > 0$ $(j = 1, 2, 3, 4)$ are some given scalars. The box $\Omega$ can be rewritten as follows:

where $\gamma_j$ is the row vector whose $i$-th element is 1 and others are zero, "${\rm{Co}}$" is the convex hull and

In addition, we define the ellipsoid $\mathscr{E}(P, \rho)$ described as

Next, we will introduce two inequalities, which are of vital importance in establishing our results.

Lemma 1. ^[28] (Wirtinger inequality) Let $x(t)\in W[a, b)$ and $x(a) = 0$. Then, for any $n\times n$ matrix $Z > 0$, the following integral inequality holds:

Lemma 2. ^[33] (Wirtinger-based inequality) Let the differentiable vector function $x(t)$, and the scalars $a$ and $b$ ($b > a$) be given. Then, for any $n\times n$ matrix $Z > 0$, the following inequalities are true:

where

3.   Main results

Here, we will consider the stabilization problem for the hyperchaotic finance system (2) in a local framework by using a piecewise discontinuous L-K functional and the quadratic system theory.

Theorem 1. Let the scalars $h > 0$, $\eta > 0$, $\varepsilon_j > 0$ ($j = 1, 2, 3, 4$), $\alpha > 0$ and $\delta\neq0$ be given. The hyperchaotic finance system (2) is asymptotically stabilized for all initial conditions $e_0$ satisfying the constraint $V(0)\leq {\bf e}^{-\alpha \eta}$ via the controller (6) with the gain $K = X^{-1}Y$, if there exist $8\times8$ matrix $P = (P_{ij})_{2\times2} > 0$, and $4\times4$ matrices $X$, $Y$, $Q > 0$, $R > 0$, $S > 0$, $Z > 0$, such that the LMIs

are satisfied, where

Proof. Construct a piecewise augmented discontinuous L-K functional

where

with

Using Lemma 1 (Wirtinger inequality), it can be seen that $V_R(t)\geq0$. Moreover, noting $V_R(t_k) = 0$, it follows that ${\rm lim}_{t\rightarrow t_k^-}V(t)\geq V(t_k)$. In addition, one can see that $V(t)$ is continuous at $t = t_0$.

By some calculations, we have

Using Lemma 2 and denoting

it follows that

Using the closed-loop system (10), it is seen that

Then, for any scalar $\delta\neq0$, we have the following zero equation:

Adding the left-hand side of (23) to $\dot{V}(t)$ and using (22), one obtains

where

and

with

Denote $Y\triangleq XK$ and notice that $\Phi(e)$ is affine about the states $e_j$, $j = 1, 2, 3, 4$. Then, it can be seen that, if the LMIs (16) holds, the relation $\Phi(e) < 0$ is ensured on $\Omega$. Moreover, on the box $\Omega$, we get

Using the open-system system (12), we have

Similarly, within the first time-interval $[0, \eta)$, we can obtain

where

and

with

If the LMIs (17) are satisfied, the matrix inequality $\hat{\Phi}(e) < 0$ can be guaranteed on the box $\Omega$. Then, one can obtain from (26) that

Moreover, it follows that

Noting that ${\rm lim}_{t\rightarrow t_k^-}V(t)\geq V(t_k)$, and using (25) and (27), it can be seen that

On the other hand, noting (20) and using Jensen integral inequalities ^[34], we have

where the relation

is utilized, and

If the LMI (18) is true, then we can get from (29) that

In addition, it is seen from the LMIs (19) that

For any $e\in \mathscr{E}(S, 1)$, we have

which implies that the following relation is true:

For any initial condition $e_0$ satisfying the constraint $V(0)\leq {\bf e}^{-\alpha \eta}$, from (27), (28) and (30), we have

which means that the system state $e(t)$ is evolved in the ellipsoid $\mathscr{E}(S, 1)$. Moreover, using (31), it is seen that the system state $e(t)$ will be evolved in the box $\Omega$.

Then, noting (25) and using the relation ${\rm lim}_{t\rightarrow t_k^-}V(t)\geq V(t_k)$, one can conclude that the closed-loop system (10) is asymptotically stable for all $e_0$ satisfying the constraint $V(0)\leq {\bf e}^{-\alpha \eta}$. The proof is completed. □

Remark 3. The proposed L-K functional (20) is continuous at the instant $t = t_0$ and discontinuous at the instants $t_k$, $k = 1, 2, \cdots$ ^[28]. Moreover, the functional (20) is piecewise. Using the functional (20), the local stability of the closed-loop system (10) can be rigorously analyzed by sufficiently considering the evolution of the open-loop system (12) within the first time-interval $[0, \eta)$.

For the case that $\eta = 0$, one can select the simplified discontinuous L-K functional

Then, a simplified stabilization criterion is readily obtained as follows:

Corollary 1. Let the scalars $h > 0$, $\varepsilon_j > 0$ ($j = 1, 2, 3, 4$) and $\delta\neq0$ be given. The conclusion of Theorem 1 holds for the case $\eta = 0$, if there exist matrices $X$, $Y$, $P > 0$, $R > 0$, such that the LMIs

are satisfied, where

Proof. Along the proof of Theorem 1, it is seen that

where

and

with

From the LMIs (33), it is seen that $\Pi(e) < 0$ is ensured on the box $\Omega$. Moreover, from (35), we have

In addition, it is inferred from the LMIs (34) that

For any $e_0$ satisfying $\check{V}(0) = e_0^TPe_0\leq 1$, using (36) and (37), and noting $\check{V}(t)\geq e^{T}(t)Pe(t)$, it is inferred that the trajectory $e(t)$ is contained in the box $\Omega$. Moreover, using (36), it can be concluded that the closed-loop system (10) is locally asymptotically stable, and this completes the proof. □

In the sequel, we will address the estimate of the DOA. Here, we employ the ellipsoid $\mathscr{E}(\mathcal{P}, {\bf e}^{-\alpha \eta})$ as the the estimate of the DOA ^[35,36]. Note that

Then, we have

where

For any initial condition $e_0$ belongs to $\mathscr{E}(\mathcal{P}, {\bf e}^{-\alpha \eta})$, it is clear that the relation $V(0)\leq {\bf e}^{-\alpha \eta}$ is guaranteed. Let us introduce the following LMI:

Then, the optimization with respect to the ellipsoid $\mathscr{E}(\mathcal{P}, {\bf e}^{-\alpha \eta})$ in Theorem 1 can be described as

For the case that $\eta = 0$, we introduce the following matrix inequality:

The optimization problem about the estimate of the DOA (i.e., the ellipsoid $\mathscr{E}(P, 1)$) is given as

Remark 4. In ^[25,26], the DTOs-based time-delayed feedback has been proposed to stabilize the stochastic hybrid differential equations. However, the nonlinearities in ^[25,26] are assumed to satisfy the rigorous linear growth conditions. It is obvious the results in ^[25,26] cannot be applicable for the finance system (2). In ^[27,28], the similar control scheme has been utilized to stabilize linear systems.

Again, it is seen that the results in ^[28] are no longer applicable for nonlinear system (2). Moreover, different from the existing results, our obtained stabilization criteria are in a local framework. In particular, the state evolution within $[0, \eta)$ is specifically taken into account. It is obvious that the proposed results in this paper are essentially the significant supplements of some existing ones.

Remark 5. Over the past two decades, the fractional-order systems have become an extremely active research field ^[38,39]. Different from the integer-order systems, the fractional-order systems can possess memory. Note that some financial variables often possess very long memory. Therefore, it has been identified that the fractional-order models should be more appropriate to describe the dynamical behaviors in financial systems ^[3,10,18,19]. As the further research topic, we would like to address the local stabilization problem for fractional-order financial systems.

Remark 6. Over the past a decade or so, the event-triggered mechanisms have been extensively employed in network-based control systems ^{[40,41,42,43,44]}. Under the event-triggered mechanisms, the necessary data are released only when certain triggering condition is satisfied, thereby significantly decreasing the usage of communication resources. It is obvious that the event-triggered mechanisms can also be applicable to the finance system (2) under which the number of state observation and control implementation can be greatly reduced, which is our further research topic.

4.   Numerical simulations

In this section, we will demonstrate the feasibility of the obtained results via numerical simulations. Here, we choose $a = 0.9$, $b = 0.2$, $c = 1.5$, $l = 0.2$ and $m = 0.17$. In Figure 1, we plot the phase portraits of the finance system (2).

Figure 1 shows the hyperchaotic behaviour. From Figure 1, one can see that the finance system (2) displays the sophisticated hyperchaotic behaviour. Then, it is verified that the system (2) has three unstable equilibrium points

First, we will be concerned with the case without time delay. For the equilibrium points $\mathbb{P}_2^*$ and $\mathbb{P}_3^*$, letting $\varepsilon_1 = 7.3$, $\varepsilon_2 = 13$, $\varepsilon_3 = 10$, $\varepsilon_4 = 12$, $h = 1$ and $\delta = 0.05$, and solving problem (2), one obtains

Using the above parameters, we plot the state evolutions of the error system (8) for the case $\eta = 0$. In the simulation, we choose $e_0 = [6, 5, 3, 1]^T\in\mathscr{E}(P, 1)$. From Figures 2 and 3, it is seen that our proposed DTOs-based time-delayed feedback scheme can stabilize the unstable hyperchaotic finance system. Figure 2 shows that the error state converges to the origin. Figure 3 shows that the error state converges to the origin.

Next, we will consider the case with time delay. Letting $\varepsilon_1 = 7.2$, $\varepsilon_2 = 12$, $\varepsilon_3 = 17$, $\varepsilon_4 = 100$, $\eta = 0.5$, $h = 1$, $\alpha = 0.16$ and $\delta = 0.05$, and solving the optimization problem (1), we have

Using the above obtained parameters, the state evolutions of the error system (8) are plotted in Figures 4 and 5, where the initial condition is selected as $e_0 = [6, 4, 3, 1]^T\in\mathscr{E}(\mathcal{P}, {\bf e}^{-\alpha\eta})$. Figures 4 and 5 show again that our proposed control scheme can effectively stabilize the unstable hyperchaotic finance system. However, compared with the case without time delay, it is seen that from Figures 4 and 5 that the convergence rate of the error system (8) becomes slower due to the existence of time delay. Figure 4 shows that the error state converges to the origin. Figure 5 shows that the error state converges to the origin.

In solving problems (1) and (2), we employ the "mincx" solver involved in LMI toolbox in MATLAB to numerically solve the minimization problem of a linear objective function subject to LMI constraints ^[45]. In the simulation, we utilize the Euler method, where the step size is selected as 0.01.

5.   Conclusions

In the paper, we have investigated the local stabilization design for a hyperchaotic finance system via the time-delayed feedback based on DTOs. By incorporating quadratic system theory, a piecewise augmented discontinuous L-K functional, and two advanced inequalities, a local stabilization criterion has been obtained in the framework of LMIs. In the case of no time delay, the corresponding result is also proposed. Then, the optimization problems have been provided to estimate the DOA as large as possible. The feasibility of proposed results has been illustrated by simulation results. The proposed techniques in this paper can be extended to the synchronization control problem ^[8,10,46].

However, it is worth mentioning that the obtained results in this paper are conservative to a certain extent. As the further improvement direction, we can employ the more effective Bessel-Legendre inequality to deal with the time delay ^[47]. In addition, we can the design the nonlinear feedback controller to reduce the potential conservatism ^[32]. On the other hand, the time delay might be time-varying ^[48,49]. Moreover, the external disturbances might be inevitable in the finance system ^[15,50]. As the further research topic, it is also interesting to address the local stabilization problem for the hyperchaotic finance system subject to external disturbances and time-varying delay.

Use of AI tools declaration

The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

Acknowledgement

This work was supported in part by the National Natural Science Foundation of China under Grant (No. 62273132), and in part by the Natural Science Foundation of Henan Province of China under Grant (No. 202300410159).

Conflict of interest

The authors declare that they have no conflicts of interest.