Abstract

The problem of finite-time stability of switched genetic regulatory networks (GRNs) with time-varying delays via Wirtinger’s integral inequality is addressed in this study. A novel Lyapunov–Krasovskii functional is proposed to capture the dynamical characteristic of GRNs. Using Wirtinger’s integral inequality, reciprocally convex combination technique and the average dwell time method conditions in the form of linear matrix inequalities (LMIs) are established for finite-time stability of switched GRNs. The applicability of the developed finite-time stability conditions is validated by numerical results.

1. Introduction

In recent years, GRNs have received much research attention, and many interesting results have been reported [18]. Generally, there are two types of gene network models, the Boolean model [9] and differential equation model [10]. In the Boolean model, the state converges to a terminal state via a series of state transitions that is determined by the Boolean rules. In this model, the activity of each gene is expressed in one of two states, ON or OFF, which is determined by a Boolean function by its own and by other related states. Whereas in the differential equation model, the variables describe the concentrations of gene products, such as mRNA and proteins as continuous values of the gene regulation system, and also, this model talks about the concentrations of gene products such as mRNA and proteins as variables in GRNs [1114]. There are many research results on the stability analysis for GRNs with time delay (e.g., [1518]).

Time delays are ubiquitous in many fields because of finite propagation speeds of signals, finite processing times, finite reaction times, and finite switching speed of amplifiers. Since the biological system especially GRNs is a slow process of transcription, translation, and translocation [1923], the time delay cannot be avoided. From the long term investigations, time delay will bring instability of the system, sustained oscillations such as bifurcation [2426]. So, it is of great importance to deal delayed GRNs. For instance, in [27], authors presented a different equation model for GRNs with constant time delays and proposed stability analysis for GRNs with time delays. In [28], authors developed delay dependent criteria for stability of GRNs with delay and free-weighting matrices.

On the other hand, Markovian switching and switched systems have been studied extensively over the past decades, for its capacity in modeling practical systems and its potential applications. Including a variety of subsystems as constituent parts, switched systems are governed by a switching rule to coordinate the switching. Recently, the stability problem of Markovian jump GRNs and switched GRNs has been investigated in [2934]. As we all know, most of gene networks contain some kinds of switching mechanisms. For instance, by increasing stimulation or by changing some regulatory mechanisms, a bistable system can switch from one steady state to the other. In [33, 3539], authors investigated the stability for switched systems with time delays by utilizing an average dwell time approach.

Recently, many kinds of finite-time issues have attracted particular research interests, and there have been some results on finite-time stabilization and synchronization [4047]. However, to the best of the authors knowledge, there have been very few results on the finite-time stability problem for delayed GRNs with time delays [48, 49], and the purpose of this study is therefore to shorten such a gap.

Motivated by the above discussion, in this study, we are concerned with the finite-time stability of switched GRNs, where the parameter values switch from one mode to another. By utilizing the average dwell time approach and by using a novel Lyapunov–Krasovskii functional, it is shown that the finite-time stability problem is solvable if a set of linear matrix inequalities (LMIs) is feasible. Finally, three examples are provided in the end of the study to show the effectiveness of the proposed criteria.

The rest of this study is organized as follows: in Section 2, preliminaries and problem formulation are given. In Section 3, some conditions are established to ensure the finite-time stability of the considered system. In Section 4, three examples are illustrated to show the effectiveness of the obtained theoretical results. And finally, conclusions are given in Section 5.Notations: throughout this study, , , and denote, respectively, the set of all real numbers, real n-dimensional space, and real -dimensional space. denote the Euclidean norms in . For a vector or matrix , denotes its transpose. For a square matrix , and denote the maximum eigenvalue and minimum eigenvalue of matrix , respectively, and sym is used to represent . For simplicity, in symmetric block matrices, we often use to represent the term that is induced by symmetry.

2. Problem Description and Preliminaries

Consider the following nonlinear GRNs with time-varying delays described bywhere , :, are the concentrations of mRNA and protein, respectively; is the regulatory functions of mRNAs, and are the constant matrices, and they are the rates of degradation; represents the translation rate; is the regulative matrix, and and are the time-varying delays.

For obtaining our conclusions, we make the following assumptions.

Assumption 1. , s = 1, 2, …, n are monotonically increasing functions with saturation and satisfywhere , s = 1, 2, …, n are the nonnegative constants.

Assumption 2. and are the time-varying delays satisfyingwhere are the constants. The initial condition of system (1) is assumed to beUse the following transformation:where and constitute an equilibrium point of system (1) and then shift the intended equilibrium point to the origin. In this way, the system equation turns to bewhere , and .
According to Assumption 1 and the definition of , we know that is bounded, that is, , such that , s = 1, 2, …, n and satisfies the following sector condition:Let .
Sometimes, GRNs were described by the continuous time switched system, as in [33]; so system (6) can be described as the switching system with switching signal:where is the switching signal, which is a piece constant function depending on time t. For each , the matrices are constant matrices of appropriate dimensions.
For the switching signal , we have the following switching sequence: ; in other words, when , subsystem is activated. To assume , .
For proving the theorem, we recall the following definition and lemmas.

Definition 1 (See [50]). The system (8) is said to be finite-time stable with respect to positive real numbers , ifwhere

Definition 2 (See [51]). For any , let denote the switching number of on an interval . Ifholds for given , , then the constant is called the average dwell time and is the chatter bound. Without loss of generality, we choose throughout this study.

Lemma 1 (See [52]). For any constant matrix , , scalars , and vector function such that the integrations concerned are well defined, and the following inequality holds

Lemma 2 (See [53]). Let have positive values in an open subset D of . Then, the reciprocally convex combination of over D satisfiessubjected to

Lemma 3 (See [54]). For a positive definite matrix , the following inequality holds for all continuously differentiable function in :where , and .

3. Main Results

In this section, we present a finite-time stability theorem for switching genetic regulatory networks with interval time-varying delays (8).

Theorem 1. The switched genetic networks (8) is finite-time stable with respect to positive real numbers and constants , and ; if there exist symmetric positive definite matrices for all , the diagonal matrix , and positive scalars and such that the following LMIs holdand the average dwell time of the switching signal p (t) satisfieswhere

Proof. Choose the Lyapunov functional candidate aswhereTaking the derivatives of along the trajectory of system (8), we have thatFrom Lemmas 1 and 3, we haveand from Lemma 2, we can obtainSimilarly, we havewhereMeanwhile, for any , the following inequality is true from Assumption 1.It can written asWhat is more, the following equations are true for any matrices with appropriate dimensions from system (8).From (24) to (39), we have thatwhereBy condition (17), we haveNote thatFor any , let , and we denote and as the switching times on the interval , whereBy integrating (43) for any , we find thatFrom (18), we can obtainBy substituting (46) into (45), we can obtainOn the other hand, it follows from (22) thatThen, we can easily obtainFrom (49) and (50), we obtainTherefore, by Definition 1, we conclude that . This completes the proof.

Remark 1. We introduce in our Lyapunov–Krasovskii functional. The novel Lyapunov–Krasovskii functional can make the stability criteria applicable to both fast and slow time-varying delays directly. Besides, by using the convex combination technique together with the Jensen inequality lemma, less conservative criteria are obtained.
Case: we consider the following genetic regulatory networks without switched term:Based on Theorem 1, the next rate-independent corollary is derived.

Corollary 1. The genetic regulatory network (52) is asymptotically stable with respect to positive real numbers ; if there exist symmetric positive definite matrices , the diagonal matrix such that the following LMIs holdwhere,

Proof. The corollary follows by a similar argument as that in proof of Theorem 1.

4. Numerical Example

In this section, numerical examples are provided to illustrate the validity and the advantage of the proposed finite-time stability of switched GRNs with time-varying delays.

Example 1. Consider switched GRNs with time-varying delay (8) aswithThe activation function is chosen as , and the values of are given as follows:We show the simulation result of the trajectories of the variables and in Figures 1 and 2. It should be point out that the condition is feasible when employing the LMI toolbox in MATLAB, solve LMIs (16)–(20), and then, we can reach feasible solution. Hence, the switched GRNs (8) is finite-time stable.

Example 2. and .
Furthermore, for the parameters listed above, let and .
In order to compare the results in [1518], using Corollary 1, the comparison results are listed in Table 1 for . Clearly, the results proposed in this study provide a larger admissible upper bound delay to guarantee the asymptotically stable system (52). In addition, the trajectories of the genetic regulatory network (52) are shown in Figures 3 and 4.

Remark 2. The discussion in Example 2 illustrates that the conditions in this study (Corollary 1) is less conservative than those in [1518], which shows the superiority of our method compared with that in [1518].

Example 3. Consider the following switched GRNs with time-varying delay:withThe activation function is chosen as , and the values of are given as follows:By employing the LMI toolbox in MATLAB, solve LMIs (16)–(20), and the feasible solutions are then reached.

5. Conclusion

In this study, a finite-time stability analysis for switched GRNs with time-varying delays has been investigated. We utilized the reciprocally convex combination method, Wirtinger’s integral inequality, and new triple integral with exponential function in Lyapunov–Krasovskii functionals; a less conservative LMI-based finite-time stability criterion is obtained with the switched ADT approach to reduce the conservatism of our results, compared with existing ones. A numerical example has been given to demonstrate the effectiveness and the advantage of our proposed methods. State estimation as well as other research topics such as switched lure systems and complex networks [55, 56] and stabilization of probabilistic Boolean networks [57, 58] of the time delay systems will be further investigated based on the methods proposed in this study.

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Authors’ Contributions

B.H. involved in funding acquisition; G.R. conceptualized the study, wrote the original draft, and validated; G.R., B.H., M.S.A., and S.S. developed software; G.R. and B.H. performed formal analysis, developed methodology, and reviewed and edited the manuscript; S.S., B.P., G.K.T., and M.S.A. supervised the study. All authors have read and agreed to the published version of the manuscript.

Acknowledgments

This work was supported by the Basic Sciences and Mathematics, Faculty of Engineering, Thai-Nichi Institute of Technology, Bangkok, 10250, Thailand.