1 Introduction

Switched control systems have drawn considerable attention for the past decades since they play an essential role in numerous applications such as mechanical systems, air traffic control systems, etc.[13]. Recently some attempts on switched systems have been investigated theoretically[46]. In the existing literatures, the study of switched nonlinear systems mainly focuses on stability and stabilization under arbitrary switchings[79] or under some designed switching law[10,11]. The stability for arbitrary switching sequences is, in fact, a kind of robustness. On the other hand, stochastic nonlinear systems have been widely studied in theory and applied in practice. In the constructive methods for stabilization of stochastic nonlinear systems, the usage of Lya-punov functions can be mainly divided into two types: One is quadratic Lyapunov functions[1214], and the other is the quartic Lyapunov function which can sometimes result in a more simple design algorithm than the former one[15,16]. As an intersecting issue of the above two aspects, switched stochastic nonlinear systems have become an emerging hot research topic for the past decade. These researches are not only interesting and challenging, but also theoretically and practically significant[1720].

As is known to all, the output feedback control is one of the most important problems of nonlinear systems since only the plant output can be measured in many cases. Generally speaking, the separation principle does not work for nonlinear systems. This makes the design of output feedback control complicated and difficult. Thus, it motivated many scholars to work on this interesting topic. Nale et al.[21] investigated the stability of switched systems with output feedback control methodology. Deng and Krstic[22] concentrated on the output feedback stabilization of stochastic nonlinear systems. Wu[23] further considered the backstepping control problem of stochastic nonlinear systems with Markovian switching, where the subsystems were assumed to be in the output feedback form. The nonlinear systems in the output feedback form are a kind of significant nonlinear systems. The feedback controllers of such a kind of nonlinear systems can be constructed systematically by the backstepping approach and its variations[2428].

In this paper, we are interested in stability analysis and control synthesis of switched stochastic nonlinear systems under arbitrary switchings. Compared to full state feedback control problem, the design of the output feedback turns out to be much more difficult. By adopting the backstepping design used in [22], we extended the output feedback design problem of the stochastic systems in the strict-feedback form to a more general class of the switched stochastic nonlinear systems. In the switched systems, multiple Lyapunov functions are often employed to satisfy different performance requirements for individual subsystems. And a set of switching laws are then designed to guarantee stability of the overall switched closed-loop system. Although the existence of multiple Lyapunov functions is less conservative, a common Lyapunov function is necessary for asymptotic stability under arbitrary switching. It is uniform over the set of all switching signals for all the subsystems. In addition, the controller can not be dependent on switching signals because the nonlinear systems we studied here are under arbitrary switchings, and the switching signals are unknown a priori. So a common output feedback controller independent of switching signals is constructed in the paper such that the equilibrium at the origin of the closed-loop system is stable.

2 Problem formulation

Consider the following class of switched stochastic nonlinear systems

(1)

where x i , u and y are the state variable, the input and the outputs of the system, respectively. σ (t) : [0, +∞) → {1, 2, ⋯ , m} is a piecewise constant switching signal, ∀jM,i ∈ {1,2, ⋯ ,n}, f j,i (•) are C 1 functions with f j,i (0) = 0, g j,i (•) are r-vector-valued smooth functions with g j,i (0) = 0, and ω is an independent r-dimensional standard Wiener process.

Remark 1. The switched stochastic nonlinear system (1) has the typical structure called strict-feedback form, which has been widely studied recently, such as [29] for the Markovian switching case, [27] for the non-stochastic and non-switched case, etc. Compared to [22], the class of systems we studied is a more general nonlinear system even in the non-switched case.

The goal of this paper is to design a controller for system (1) by backstepping approach such that the zero solution of the closed-loop system is globally asymptotically stable in probability under arbitrary switchings. Further, it will be proved that the zero equilibrium of the closed-loop system is of pth-moment globally exponential stability, i.e., there exist k > 0, c > 0, p > 0 such that \(\operatorname{E} \left[ {\left\| {x(t)} \right\|_p^p} \right] \leqslant k\left\| {x(0)} \right\|_p^p{\operatorname{e} ^{ - ct}},\forall \left\| {x(0)} \right\| \in {{\mathbf{R}}^n}.\)

3 Main results

3.1 Controller design

Since states x 2 , ⋯ , x n are not measured, the observer is designed as

$$\left\{ {\begin{array}{*{20}{l}} {{{\dot \hat x}_i} = {{\hat x}_{i + 1}} + {k_i}(y - {{\hat x}_1}),\quad i = 1,2, \cdots ,n - 1} \\ {{{\dot \hat x}_n} = u + {k_n}(y - {{\hat x}_1})} \end{array}} \right.$$
(2)

where k i > 0 are the design parameters. Define the estimate error as

$$\tilde x = x - \hat x$$
(3)
$$\operatorname{d} \tilde x = \left[ {A\tilde x + {F_{\sigma (t)}}(y)} \right]\operatorname{d} t + {G_{\sigma (t)}}(y)\operatorname{d} \omega $$
(4)

where

and

Here, the design parameters k i > 0 are chosen such that A is a Hurwitz matrix. Now, the entire system can be described by

(5)

Define the following coordinate transformation

$$\left\{ {\;\,\begin{array}{*{20}{l}} {{z_1} = y} \\ {{z_2} = {{\hat x}_2} - {\alpha _1}(y)} \\ {{z_i} = {{\hat x}_i} - {\alpha _{i - 1}}\left( {{{\bar \hat x}_{i - 1}},y} \right),\quad i = 3, \cdots ,n} \end{array}} \right.$$
(6)

where α i –1 (•) are the virtual stabilizing functions to be determined later and \[{\bar \hat x_i} = {\left[ {{{\hat x}_1},{{\hat x}_2}, \cdots ,{{\hat x}_i}} \right]^\operatorname{T} }\]. Then according to \(\operatorname{It} \hat o'\operatorname{s} \) differentiation rule, system (5) can be transformed into the following equalities

(7)

Naturally, the goal is changed to design a controller for system (7) by backstepping approach such that the equilibrium at the origin of the closed-loop system is globally asymptotically stable in probability under arbitrary switchings. Similar to [22], instead of constructing the stabilizing functions α i and u in a step-by-step fashion, we derive them simultaneously. In this paper, the stabilizing functions α i and u are constructed as

$$\left\{ {\begin{array}{*{20}{l}} {{\alpha _1}(y) = {z_1}{s_1}(y)} \\ {{\alpha _i}({{\bar \hat x}_i},y) = {z_i}{s_i}({{\bar \hat x}_i},y),\quad i = 2, \cdots ,n - 1} \\ {u = {z_n}{s_n}({{\bar \hat x}_n},y)} \end{array}} \right.$$
(8)

where si (•) to be determined later are all smooth functions defined on R i .

3.2 Stability analysis

Theorem 1. For system(8), there exists a control law taking the form of (7) such that the equilibrium at the origin of the closed-loop system is fourth-moment exponentially stable.

Proof. Since ∀jM, f j,i (y) and g j,i (y) are smooth functions with f j,i (0) = 0 and g j,i (0) = 0, by the mean value theorem, f j,i (y) and g j,i (y) can be expressed as

$${f_{j,i}}(y) = y{\varphi _{j,i}}(y)$$
(9)
$${g_{j,i}}(y) = g{\psi _{j,i}}(y)$$
(10)

where φ j,i (y) and ψ j,i (y) are both smooth functions.

Because of the term \(\frac{1}{2}\) tr caused by stochastic case, we employ quartic Lyapunov functions. Define a Lyapunov function of the form

(11)

where P is a positive defined matrix which satisfies A T P + PA = – I. Then according to the \(\operatorname{It} \hat o'\operatorname{s} \) differentiation rule, we have

(12)

where

Next, the terms tr and \(z_i^3{k_i}{\tilde x_1}\) are further handled one by one such that their effects on the negativity of LV can be canceled by s i selected properly. To handle these terms, a special case of the well-known Young’s inequality[25] plays a very important role, which states that the inequality

$$xy \leqslant \frac{{{ \in ^p}}}{p}{\left| x \right|^p} + \frac{1}{{{ \in ^q}q}}{\left| y \right|^q},\quad \forall (x,y) \in {{\mathbf{R}}^2}$$
(13)

holds for any constants p > 1 and q > 1 satisfying (p−1 )(q–1) = 1.

By simple calculation and using Youngs inequality, one has

$$\sum\limits_{i = 2}^{n - 1} {z_i^3} {z_{i + 1}} \leqslant \frac{3}{4}\sum\limits_{i = 2}^{n - 1} {\delta _i^{\tfrac{4}{3}}z_i^4 + \frac{1}{4}\sum\limits_{i = 3}^n {\frac{1}{{\delta _{i - 1}^4}}z_i^4} } $$
(14)
$${y^3}{\tilde x_2} \leqslant \frac{3}{4} \in _1^{\tfrac{4}{3}}{y^4} + \frac{1}{{4 \in _1^4}}\tilde x_2^4 \leqslant \frac{3}{4} \in _1^{\tfrac{4}{3}}{y^4} + \frac{1}{{4 \in _1^4}}{\left\| {\tilde x} \right\|^4}$$
(15)
$${y^3}{z_2} \leqslant \frac{3}{4}\delta _1^{\tfrac{4}{3}}{y^4} + \frac{1}{{4\delta _1^4}}z_2^4$$
(16)
$$\begin{array}{*{20}{l}} { - z_i^3\sum\limits_{l = 2}^{i - 1} {\frac{{\partial {\alpha _{i - 1}}}}{{\partial {{\hat x}_l}}}{k_l}{{\hat x}_1} \leqslant } } \\ {\quad \;\;\frac{3}{4}{{\left( {{\kappa _i}\sum\limits_{l = 2}^{i - 1} {\frac{{\partial {\alpha _i} - 1}}{{\partial {{\hat x}_l}}}{k_l}} } \right)}^{\tfrac{4}{3}}}z_i^4 + \frac{1}{{4\kappa _i^4}}\tilde x_1^4 \leqslant } \\ {\quad \;\;\frac{3}{4}{{\left( {{\kappa _i}\sum\limits_{l = 2}^{i - 1} {\frac{{\partial {\alpha _i} - 1}}{{\partial {{\hat x}_l}}}} {k_l}} \right)}^{\tfrac{4}{3}}}z_i^4 + \frac{1}{{4\kappa _i^4}}{{\left\| {\tilde x} \right\|}^4}} \end{array}$$
(17)
(18)
(19)
$$\begin{array}{*{20}{l}} { - z_i^3\frac{{\partial {\alpha _{i - 1}}}}{{\partial y}}{{\hat x}_2} = } \\ {\quad - z_i^3\frac{{\partial {\alpha _{i - 1}}}}{{\partial y}}{z_2} - z_i^3\frac{{\partial {\alpha _{i - 1}}}}{{\partial y}}y{s_1} \leqslant } \\ {\frac{3}{4}z_i^4{{\left( {\frac{{\partial {\alpha _{i - 1}}}}{{\partial y}}} \right)}^{\tfrac{4}{3}}} + \frac{1}{4}z_2^4 + \frac{1}{4}{y^4} + \frac{3}{4}z_i^4{{\left( {{s_1}\frac{{\partial {\alpha _{i - 1}}}}{{\partial y}}} \right)}^{\tfrac{4}{3}}}} \end{array}$$
(20)
$$\begin{array}{*{20}{l}} { - \sum\limits_{i = 2}^n {z_i^3\frac{{\partial {\alpha _{i - 1}}}}{{\partial y}}{{\tilde x}_2} \leqslant } } \\ {\quad \frac{3}{4}\sum\limits_{i = 2}^n {\eta _i^{\tfrac{4}{3}}} {{\left( {\frac{{\partial {\alpha _{i - 1}}}}{{\partial y}}} \right)}^{\tfrac{4}{3}}}z_i^4 + \frac{1}{4}\sum\limits_{i = 2}^n {\frac{1}{{\eta _i^4}}\tilde x_2^4} \leqslant } \\ {\quad \frac{3}{4}\sum\limits_{i = 2}^n {\eta _i^{\tfrac{4}{3}}} {{\left( {\frac{{\partial {\alpha _{i - 1}}}}{{\partial y}}} \right)}^{\tfrac{4}{3}}}z_i^4 + \frac{1}{4}\sum\limits_{i = 2}^n {\frac{1}{{\eta _i^4}}{{\left\| {\tilde x} \right\|}^4}} } \end{array}$$
(21)
$$\begin{array}{*{20}{l}} { - z_i^3\frac{{\partial {\alpha _{i - 1}}}}{{\partial y}}y{\varphi _{\sigma (t),1}}(y) \leqslant } \\ {\quad \;\; - \frac{3}{4}{{\left( {\frac{{\partial {\alpha _{i - 1}}}}{{\partial y}}{\varphi _{\sigma (t),1}}} \right)}^{\tfrac{4}{3}}}z_i^4 - \frac{1}{4}{y^4} \leqslant } \\ {\quad \;\; - \frac{3}{4}{{\left( {\frac{{\partial {\alpha _{i - 1}}}}{{\partial y}}\sum\limits_{j = 2}^m {{\varphi _{j,1}}} } \right)}^{\tfrac{4}{3}}}z_i^4 - \frac{1}{4}{y^4}} \end{array}$$
(22)
(23)
$$\begin{array}{*{20}{l}} { - z_i^3\sum\limits_{l = 2}^{i - 1} {\frac{{\partial {\alpha _{i - 1}}}}{{\partial {{\hat x}_l}}}} {{\hat x}_{l + 1}} = } \\ {\quad \;\; - z_i^3\sum\limits_{l = 2}^{i - 1} {\frac{{\partial {\alpha _{i - 1}}}}{{\partial {{\hat x}_l}}}} ({z_{l + 1}} + {z_l}{s_l}) = } \\ {\quad \;\; - z_i^4\frac{{\partial {\alpha _{i - 1}}}}{{\partial {{\hat x}_{i - 1}}}} - z_i^3\sum\limits_{l = 2}^{i - 2} {{z_{l + 1}}} \frac{{\partial {\alpha _{i - 1}}}}{{\partial {{\hat x}_l}}} - z_i^3\sum\limits_{l = 2}^{i - 1} {{z_l}} \frac{{\partial {\alpha _{i - 1}}}}{{\partial {{\hat x}_l}}}{s_l} \leqslant } \\ {\quad \;\; - z_i^4\frac{{\partial {\alpha _{i - 1}}}}{{\partial {{\hat x}_{i - 1}}}} + \frac{3}{4}z_i^4\sum\limits_{l = 2}^{i - 2} {{{\left( {\frac{{\partial {\alpha _{i - 1}}}}{{\partial {{\hat x}_l}}}} \right)}^{\tfrac{4}{3}}} + \frac{1}{2}\sum\limits_{k = 1}^{i - 1} {z_4^k} } + } \\ {\quad \;\;\frac{3}{4}z_i^4\sum\limits_{l = 2}^{i - 1} {{{\left( {\frac{{\partial {\alpha _{i - 1}}}}{{\partial {{\hat x}_l}}}{s_l}} \right)}^{\tfrac{4}{3}}}} } \end{array}$$
(24)
$$\frac{3}{2}{\left\| {\tilde x} \right\|^4} + \frac{1}{2}\lambda _{\max }^8(P){y^4}\sum\limits_{i = 1}^n {\sum\limits_{j = 1}^m {\varphi _{j,i}^4(y)} } $$
(25)

and

$$z_i^3{k_i}{\tilde x_1} \leqslant \frac{3}{4}\theta _i^{\tfrac{4}{3}}k_i^{\tfrac{4}{3}}z_i^4 + \frac{1}{{4\theta _i^4}}\tilde x_1^4 \leqslant \frac{3}{4}\theta _i^{\tfrac{4}{3}}k_i^{\tfrac{4}{3}}z_i^4 + \frac{1}{{4\theta _i^4}}{\left\| {\tilde x} \right\|^4}.$$
(26)

Based on the above analysis, we have

$$\begin{array}{*{20}{l}} {\mathcal{L}V \leqslant - {{\left\| {\tilde x} \right\|}^4}p + {y^4}({s_1} + A) + z_2^4({s_2} + B) + } \\ {\quad \quad \;\;\sum\limits_{i = 3}^{n - 1} {z_i^4} ({s_i} + C) + z_n^4({s_n} + D)} \end{array}$$
(27)

where ∀jM,

(28)
(29)
(30)
(31)

and the following equality is used.

$$\begin{array}{*{20}{l}} {\frac{1}{2}\sum\limits_{i = 3}^n {\sum\limits_{k = 1}^{i - 1} {z_k^4} } = } \\ {\quad \;\;\;\;\frac{1}{2}\sum\limits_{i = 3}^{n - 1} {z_i^4} (n - i) + \frac{1}{2}z_1^4(n - 2) + \frac{1}{2}z_2^4(n - 2).} \end{array}$$
(32)

In (27), ε, δ, η, θ, κ and ξ are positive constants to be chosen. λ > 0 is the smallest eigenvalue of P. We will choose ∈1 , ∈2 , ∈3 i , κi and η i to satisfy

$$\begin{array}{*{20}{l}} {p = 2\lambda - \frac{3}{2} \in _3^{\tfrac{4}{3}} - \frac{{3n\sqrt n \in _2^2{{\left\| P \right\|}^4}}}{2} - \frac{1}{{4 \in _1^4}} - \frac{1}{4}\sum\limits_{i = 2}^n {\frac{1}{{\eta _i^4}}} - } \\ {\quad \;\;\frac{1}{4}\sum\limits_{i = 2}^n {\frac{1}{{\theta _i^4}}} - \frac{1}{4}\sum\limits_{i = 3}^n {\frac{1}{{\kappa _i^4}} > 0.} } \end{array}$$
(33)

Obviously, it is a simple task to select s i , ∀ i ∈ {1,2, ⋯ ,n} such that

$${s_1} \leqslant - {c_1} - A$$
(34)
$${s_2} \leqslant - {c_2} - B$$
(35)
$${s_i} \leqslant - {c_i} - C,\quad i = 3, \cdots ,n - 1$$
(36)
$${s_n} \leqslant - {c_n} - D$$
(37)

where c i > 0. Then we get

$$\mathcal{L}V \leqslant - \sum\limits_{i = 1}^n {{c_i}z_i^4} - p{\left\| {\tilde x} \right\|^4} \leqslant - cV$$
(38)

where \(c = \min \left\{ {4{c_i},\frac{{2p}}{{\lambda _{\max }^2}}} \right\},\) i = 1,2, ⋯ ,n, λ max > 0 is the maximal eigenvalue of P.

Observe that

(39)

For a given time t, suppose that switchings occur at the time instants t 1 , ⋯ ,t K , thus the time interval can be divided as [0,t 1 ) ∪ ⋯ ∪ [t K –1 , t K ) ∪ [t K , t]. Integrating both sides of the above inequality on each subinterval and taking expectations yield

$$\operatorname{E} \left[ {{\operatorname{e} ^{ct}}V(x(t))} \right] - \operatorname{E} \left[ {{\operatorname{e} ^{c{t_K}}}V(x({t_K}))} \right] \leqslant 0$$
(40)

and

$$\operatorname{E} \left[ {{\operatorname{e} ^{ct_K^ - }}V(x(t_K^ - ))} \right] - \operatorname{E} \left[ {{\operatorname{e} ^{c{t_{k - 1}}}}V(x({t_{k - 1}}))} \right] \leqslant 0$$
(41)

where k = 1, ⋯ , K − 1 and t 0 = 0. Since the system state does not jump at the switching instants and V is a common Lyapunov function for all subsystems

$$\operatorname{E} \left[ {{\operatorname{e} ^{ct_k^ - }}V(x(t_k^ - ))} \right] = \operatorname{E} \left[ {{\operatorname{e} ^{c{t_k}}}V(x({t_k}))} \right].$$
(42)

Then it follows from (40) and (41) that

$$\operatorname{E} \left[ {V(x(t))} \right] \leqslant \operatorname{E} \left[ {V(0)} \right]{\operatorname{e} ^{ - ct}} = V(0){\operatorname{e} ^{ - ct}}.$$
(43)

As we stated previously, the Lyapunov function we employ satisfies

(44)

where Y = [z 1 ,z 2 , ⋯ , z n , x1 , x 2, ⋯ ,x n ]T. Hence,

$$\operatorname{E} \left[ {V(x(t))} \right] \leqslant \frac{{n\left( {1 + 2\lambda _{\max }^2} \right)}}{4}\left\| {Y(0)} \right\|_4^4{\operatorname{e} ^{ - ct}}.$$
(45)

4 An illustrative example

In the following, an example shown in [30] is presented to illustrate the effectiveness of our result. The continuously stirred tank reactor with two-mode feeding stream is molded as the following system

$$\left\{ {\begin{array}{*{20}{l}} {{{\dot x}_1} = {x_2} + {f_{\sigma (t),1}}(y)} \\ {{{\dot x}_2} = u} \end{array}} \right.$$
(46)

with \({f_{1,1}}({x_1}) = \frac{1}{2}{x_1}\) and f 2 ,1 (x 1 ) = 2x 1 . It is assumed that there exists white noise due to the term f σ (t),1 in the above system. As a result, we have the following stochastic nonlinear system

$$\left\{ {\begin{array}{*{20}{l}} {\operatorname{d} {x_1} = \left[ {{x_2} + {f_{\sigma (t),1}}(y)} \right]\operatorname{d} t + {f_{\sigma (t),1}}(y)\operatorname{d} \omega } \\ {\operatorname{d} {x_2} = u\operatorname{d} t} \\ {y = {x_1}.} \end{array}} \right.$$
(47)

For this system, the estimator is

$$\begin{array}{*{20}{l}} {{{\dot \hat x}_1} = {{\hat x}_2} + {k_1}(y - {{\hat x}_1})} \\ {{{\dot \hat x}_2} = u + {k_2}(y - {{\hat x}_1}).} \end{array}$$

It is easy to obtain that \({\varphi _{1,1}} = \frac{1}{2}\), φ 2,1 = 2, \({\psi _{1,1}} = \frac{1}{2}\), ψ 2,1 = 2 and φ j, 2 = ψ j, 2 = 0 for j = 1, 2. Based on(30), a smooth function s 1 should be chosen such that

(46)

According to the first equation of (8), the virtual control is expressed as

$${\alpha _1} = {s_1}y.$$
(47)

Then, based on (37), a smooth function s 2 should be chosen such that

$$\begin{array}{*{20}{l}} {{s_2} = - {c_2} - \left[ {\frac{1}{{4\delta _1^4}} + \frac{3}{4}\theta _2^{\tfrac{4}{3}}k_2^{\tfrac{4}{3}} + \frac{3}{4}{{\left( {\frac{{\partial {\alpha _1}}}{{\partial y}}} \right)}^{\tfrac{4}{3}}} + } \right.} \\ {\quad \;\;\;\frac{3}{4}{{\left( {{s_1}\frac{{\partial {\alpha _1}}}{{\partial y}}} \right)}^{\tfrac{4}{3}}} + \frac{3}{8}{{\left( {\frac{{17{\tau _2}}}{{64}}\frac{{{\partial ^2}{\alpha _1}}}{{\partial {y^2}}}} \right)}^{\tfrac{4}{3}}} + } \\ {\left. {\quad \;\;\;\frac{3}{4}\eta _2^{\tfrac{4}{3}} + \frac{3}{4}\frac{1}{{\xi _2^2}}{{\left( {\frac{{\partial {\alpha _1}}}{{\partial y}}} \right)}^4} - \frac{3}{4}{{\left( {\frac{{5\partial {\alpha _1}}}{{2\partial y}}} \right)}^{\tfrac{4}{3}}}} \right]} \end{array}$$

which leads to the following control law by the third equation of (8)

$$u = \left[ {{{\hat x}_2} - {\alpha _1}(y)} \right]{s_2}.$$
(48)

In the numerical simulation, we chose k 1 = 3, k 2 = 4.5, C1 = − 0.01, c2 = − 0.1, δ1 = 1, ξ2 = 1, T 2 = 1, ∈1 = 3, ∈2 = 0.125, ∈3 = 1,η 2 = 3, and θ2 = 3 and set the initial condition to be x 1 (0) = 0.1, x 2 (0) = 0, \({\hat x_1}(0) = 0\) and \({\hat x_2}(0) = 0\). The switching signal generated randomly is shown in Fig. 1. The trajectories of the states and the estimates are shown in Figs. 2 and 3, respectively. From Fig. 2, one can see that the outputs of the closed-loop system converges to zero asymptotically. The simulation results show the effectiveness of the proposed control method.

Fig. 1
figure 1

The switching signal generated randomly

Fig. 2
figure 2

The curves of the state one and its estimation

Fig. 3
figure 3

The curves of the state two and its estimation

5 Conclusion

This paper has investigated the problem of global stabilization of switched stochastic output feedback nonlinear systems under arbitrary switchings by backstepping approach. A common output feedback controller independent of switching signals is constructed to guarantee that the equilibrium at the origin of the closed-loop system is fourth-moment exponentially stable. And an example is employed to show the effectiveness of the proposed method.