Constrained Feedback Stabilization for Bilinear Parabolic Systems

Azzeddine Tsouli^*, Ali Boutoulout

TSI Team, MACS Laboratory, Faculty of Sciences, Moulay Ismail University, Meknes, Morocco.

DOI: 10.4236/ica.2015.62011                

Abstract

In this paper, we shall study the stabilization and the robustness of a constrained feedback control for bilinear parabolic systems defined on a Hilbert state space. Then, we shall show that stabilizing such a system reduces stabilization only in its projection on a suitable subspace. For this purpose, a new constrained stabilizing feedback control that allows a polynomial decay estimate of the stabilized state is given. Also, the robustness of the considered control is discussed. An illustrating example and simulations are presented.

Keywords

Bilinear Parabolic Systems, Constrained Feedback Stabilisation, Decay Estimate, Robustness

1. Introduction

Bilinear systems represent a small, but important subset of nonlinear systems within which linear systems coexist as a special subclass. Adopting a bilinear model retains a well structured framework, which contains the well- known notional concepts such as time constants and steady-state behaviour. When adopting a bilinear approach, these concepts become operation-dependent quantities which can be appropriately modelled. Bilinear system models represent an important class of nonlinear models that are defined to be linear in both state and control when considered independently, with the nonlinearity (or bilinearity) arising from coupled terms involving products of system state and control (see [1] [2] ). By formulating the model appropriately, the bilinear term could also be represented by products of system output and control input, i.e. the output is defined as a system state. There are numerous combinations of product terms that could be considered, thus potentially increasing the model complexity. However, it has been found in practice that a minimal number of product terms can provide an adequate model for the purpose of control. Bilinear model structures are able to represent nonlinear phenomena more accurately than linear models, and thereby extend the range of satisfactory performance. In this paper, we are concerned with the question of the stabilization by a constrained feedback control for bilinear parabolic systems that can be described in the following form:

$\frac{\text{d}y\left(t\right)}{\text{d}t}=Ay\left(t\right)+p\left(t\right)By\left(t\right),y\left(0\right)=y_0$ (1)

on a real Hilbert space $H$ with inner product $\langle .,.\rangle$ ; and corresponding norm $\Vert .\Vert$ , where the linear operator $A$ generates a contraction semigroup ${\left(S\left(t\right)\right)}_{t\ge 0}$ on $H$ and $B\in \mathcal{L}\left(H\right)$ . While the real valued function $p(.)\in L^2\left(0,+\infty ;ℝ\right)$ represents a control. A function $y(.)\in \mathcal{C}\left(0,t_0;H\right)$ , $t_0>0$ , is a mild solution of the system (1) if and only if the solution $y\left(t\right)$ of the system (1) satisfies the variation of parameters formula:

$y\left(t\right)=S\left(t\right)y_0+{\displaystyle {\int }_0^{t}S}\left(t-\tau \right)p\left(\tau \right)By\left(\tau \right)\text{d}\tau ,t\ge 0$ (2)

(see [3] ). By choosing an adequate feedback control $p\left(t\right)$ in such a way, the corresponding solution $y\left(t\right)$ of the system (1) converges to zero when $t\to +\infty$ , for all $y_0$ in $H$ . For finite-dimensional bilinear systems associated to a skew-adjoint matrix $A$ , the question of stabilization has been treated in [4] , under the condition:

$\text{span}\left\{Az,a{d}^0\left(A,B\right)z,a{d}^1\left(A,B\right)z,\cdots ,a{d}^k\left(A,B\right)z,\cdots \right\}={ℝ}^n,\forall z\in {ℝ}^n-\left\{0\right\}$ (3)

where $a{d}^k\left(A,B\right)$ is defined recursively as $a{d}^0\left(A,B\right)=B$ , $a{d}^1\left(A,B\right)=AB-BA$ and $a{d}^{k+1}\left(A,B\right)=a{d}^1\left(A,a{d}^k\left(A,B\right)\right)$ , $\forall k\in \mathbb{N}$ .

Using the following assumptions:

$\langle B\text{exp}\left(tA\right)y,\text{exp}\left(tA\right)y\rangle =0,,\forall t\ge 0\Rightarrow y=0$ (4)

the problem of stabilization has been studied in [5] . In [3] , when the linear operator $B$ is compact and ${\left(S\left(t\right)\right)}_{t\ge 0}$ is a contraction semigroup, then using the quadratic feedback control

$p_0\left(t\right)=-\langle y\left(t\right),By\left(t\right)\rangle$ (5)

a weak stabilization result is obtained under the weak observability condition:

$\langle BS\left(t\right)y,S\left(t\right)y\rangle =0,\forall t\ge 0\Rightarrow y=0$ (6)

In the case where $B$ is sequentially continuous from $H_w$ ( $H$ endowed with the weak topology) to $H$ , the quadratic feedback control (2) weakly stabilizes the system (1), provided that the following weak observability assumption (4) holds (see [3] ). Under the exact observability assumption

${\displaystyle {\int }_0^T\left|\langle BS\left(t\right)y,S\left(t\right)y\rangle \right|\text{d}t}\ge \delta {\Vert y\Vert }^2,\forall y\in H,\left(T,\delta >0\right)$ (7)

The strong stabilization result with the following decay estimate

$\Vert y\left(t\right)\Vert =\mathcal{O}\left(\frac{1}{\sqrt{t}}\right),\text{as}t\to +\infty$ (8)

i.e. $\Vert y\left(t\right)\Vert \le \frac{M}{\sqrt{t}}$ , $M>0$ for $t$ large enough, has been obtained using the quadratic feedback control (5) (see [6] ). However, in this way the convergence of the resulting closed loop state is not better than (8). In [7] the rational decay rates are established i.e. using the following feedback control:

$p_r\left(t\right)=-\frac{\langle y\left(t\right),By\left(t\right)\rangle }{{\Vert y\left(t\right)\Vert }^r},r\in \left(-\infty ,2\right)$

It has been shown in [8] that, where the resolvent of the operator $A$ is compact, and $B$ is abounded linear, self-adjoint and monotone, the constrained feedback control law

$p\left(t\right)=-\frac{\langle By\left(t\right),y\left(t\right)\rangle }{1+\left|\langle By\left(t\right),y\left(t\right)\rangle \right|}$ (9)

strongly stabilizes the system (1), provided that the assumption (6) holds. It has been established in [7] that, if

the linear operator $A$ generates a contraction semigroup ${\left(S\left(t\right)\right)}_{t\ge 0}$ in $H$ , then the system (1) is strongly

stable with the explicit decay estimate (8), using the control (9), provided that the estimate (7) holds. Here, we will establish an explicit decay estimate of the stabilized state and the robustness of the control (9) for a large class of bilinear systems as considered in [3] [8] [9] . The method used in this paper is based on decomposing the system (1) into two suitable subsystems: the stable part and the unstable one. Then, we will show that one can concentrate on the determination of a stabilizing control for the so-called unstable part which maintains the exponential stability of the stable part. The rest of this article is as follows: in Section 2, we will give the main hypotheses that allow the decomposition of the system (1) into two subsystems. Then, under the compactness hypothesis of the operator $B$ , we will give a weaker variant of the condition (6) which achieves strong stabilization of the system (1). In Section 3, we will show that under a weaker version of (6), we obtain the stabilization with the decay estimate (8). Section 4 concerns the robustness of the stabilizing controls. The last section is devoted to an illustrating example and simulations.

2. Stabilization Results

Let us now recall the following definition concerning the asymptotic behavior of the system (1).

2.1. Definition

The system (1) is weakly (resp. strongly) stabilizable if there exists a feedback control $p\left(t\right)=f\left(y\left(t\right)\right)$ , $t\ge 0$ $f:H\to K:=ℝ$, $ℂ$ such that the corresponding mild solution $y\left(t\right)$ of the system (1) satisfies the properties:

1. For each initial state $y_0$ of the system (1) there exists a unique mild solution defined for all $t\in {ℝ}^{+}$ of the system (1),

2. $\left\{0\right\}$ is an equilibrium state of the system (1),

3. $y\left(t\right)\to 0$ , weakly (resp. strongly), as $t\to +\infty$ , for all $y_0\in H$ .

In the sequel of this section, we will present an appropriate decomposition of the state space $H$ and the system (1) via the spectral properties of the operator $A$ , and we apply this approach to study the stabilization pro- blem of the system (1). In [10] - [12] , it has been shown that if the spectrum $\sigma \left(A\right)$ of $A$ can be decomposed

into ${\sigma }_u\left(A\right)=\left\{\lambda :\mathcal{R}e\left(\lambda \right)\ge -\eta ,\eta >0\right\}$ and ${\sigma }_s\left(A\right)=\left\{\lambda :\mathcal{R}e\left(\lambda \right)<-\eta \right\}$ , then the state space $H$ can be decomposed according to

$H=H_u\oplus H_s$ (10)

where $H_u=P_{u}H=\text{vect}\left\{{\phi }_j,1\le j\le N\right\}$ , $H_s=P_{s}H=\text{vect}\left\{{\phi }_j,j>N\right\}$ , $P_u$ is given by

$P_u=\frac{1}{2\text{π}i}{\displaystyle {\int }_C{\left(\lambda I-A\right)}^{-1}\text{d}\lambda }$ (11)

where $C$ is a curve surrounding ${\sigma }_u\left(A\right)$ , $P_s=I-P_u$ and for all $j\ge 1,{\phi }_j$ is the eigenvector associated to

the eigenvalue ${\lambda }_j$ . The projection operators $P_u$ and $P_u$ commute with $A$ , and we have $A=A_u+A_s$ with $A_u=P_{u}A{P}_u$ and $A_s=P_{s}A{P}_s$ . Also, for all $y\left(t\right)\in H$ , we set $y_u=P_{u}y$ and $y_s=P_{s}y$ . For linear systems, it

has been shown that the initial system can be decomposed into two subsystems on $H_u$ and $H_s$ . If $A_s$ satisfies the spectrum growth assumption:

$\underset{t\to +\infty }{\lim }\frac{\ln \left(\Vert S_s\left(t\right)\Vert \right)}{t}=\text{supRe}\left(\sigma \left(A_s\right)\right)$ (12)

Which is equivalent to:

$\Vert S_s\left(t\right)\Vert \le M_1\text{exp}\left(-\eta t\right),\forall t\ge 0\left(\text{forsome}{M}_1>0\right)$ (13)

where ${\left(S_s\left(t\right)\right)}_{t\ge 0}$ denotes the semigroup generated by $A_s$ in $H_s$ , then stabilizing the whole system turns out to stabilizing its projection on $H_u$ (see [13] ). In the sequel, we suppose that the operator $B$ satisfies

$B{H}_u\subset H_u\text{and}B{H}_s\subset H_s$ (14)

It is easily verified that the condition (14) is equivalent to the fact that the linear operator $B$ commutes with $P_u$ . We note that the condition (14) also holds in the special case: $H_u=H$ . Let us consider that the system (1) can be decomposed in the following two subsystems:

$\frac{\text{d}{y}_u\left(t\right)}{\text{d}t}=A_u{y}_u\left(t\right)+p\left(t\right)B_u{y}_u\left(t\right),y_u\left(0\right)=y_{0u}\in H_u$ (15)

$\frac{\text{d}{y}_s\left(t\right)}{\text{d}t}=A_s{y}_s\left(t\right)+p\left(t\right)B_s{y}_s\left(t\right),y_s\left(0\right)=y_{0s}\in H_s$ (16)

in the state spaces $H_u$ and $H_S$ respectively, and $B=B_u\oplus B_s$ . It has been proved that stabilizing a linear system turns out to stabilizing its unstable part (see [13] ).

2.2. Remark

For finite-dimensional systems, the conditions (6) and (7) are equivalent (see [5] [10] ). However, in infinite- dimensional case, and if $B$ is compact, then the condition (7) is impossible. Indeed, if ${\left({\phi }_j\right)}_{j\ge 1}$ is an orthonormal basis of $H$ , then applying (7) for $y={\phi }_j$ and using the fact that ${\phi }_j\to 0$ , weakly as $j\to +\infty$ , we obtain the contradiction: $\delta =0.$

The following result concerns the strong stabilization of the system (1).

2.3. Theorem

Let

1. $A$ generates a linear $C_0$ -contraction semigroup ${\left(S\left(t\right)\right)}_{t\ge 0}$ on $H$ ,

2. $A$ allows the decomposition (10) of $H$ with $\text{dim}{H}_u<+\infty$ such that (13) holds,

3. $B$ be compact such that

$\langle BS\left(t\right)y,S\left(t\right)y\rangle =0,\forall t\ge 0\Rightarrow y=0$ (17)

Then, the constrained feedback control law:

$p_{\rho }\left(t\right)=-\rho \frac{\langle y\left(t\right),By\left(t\right)\rangle }{1+\left|\langle y\left(t\right),By\left(t\right)\rangle \right|},\rho >0$ (18)

strongly stabilizes the system (1).

Proof

The system (1) controlled by (18) possesses a unique mild solution $y(.)$ defined on a maximal interval $\left[0,t_{\text{max}}\right[$ and given by the variation of constants formula

$y\left(t\right)=S\left(t\right)y_0-\rho {\displaystyle {\int }_0^{t}S}\left(t-\tau \right)\frac{\langle y\left(\tau \right),By\left(\tau \right)\rangle }{1+\left|\langle y\left(\tau \right),By\left(\tau \right)\rangle \right|}By\left(\tau \right)\text{d}\tau$ (19)

corresponds to (18) (see [9] ).

Since ${\left(S\left(t\right)\right)}_{t\ge 0}$ is a contraction semigroup, we get:

$\frac{\text{d}{\Vert y\left(t\right)\Vert }^2}{\text{d}t}\le -2\rho \frac{{\left|\langle y\left(t\right),By\left(t\right)\rangle \right|}^2}{1+\left|\langle y\left(t\right),By\left(t\right)\rangle \right|},\forall y_0\in \mathcal{D}\left(A\right)$ (20)

It follows from (20) that

$\Vert y\left(t\right)\Vert \le \Vert y_0\Vert ,\forall t\ge 0$ (21)

From (19) and using the fact that ${\left(S\left(t\right)\right)}_{t\ge 0}$ is a contraction semigroup and the Gronwall inequality, we deduce that the map $y_0\to y\left(t\right)$ is continuous from $H$ to $H$ . Then (21) holds for all $y_0\in H$ by density argument, and hence $t_{\text{max}}=+\infty$ (see [14] ). Now, let us show that $y\left(t\right)\to 0$ , weakly as $t\to +\infty$ . Let $t_n\to +\infty$

such that $y\left(t_n\right)$ weakly converges in $H$ and let $z\in H$ such that $y\left(t_n\right)\to z$ , weakly as $n\to \infty$ . (The ex- istence of the sequence ${\left(t_n\right)}_{n\in \mathbb{N}}$ and $z$ are ensured by (21) and by the fact that space $H$ is reflexive.) Taking $y\left(t_n\right)$ as initial state in (19) and using superposition property of the solution, and via the dominated conver-

gence theorem, we obtain $\langle BS\left(t\right)z,S\left(t\right)z\rangle =0$ , $\forall t\ge 0$ . It follows from (17) that $z=0$ . Hence $y\left(t\right)\to 0$ , weakly as $t\to +\infty$ , and since $\text{dim}{H}_u<+\infty$ , we have $y_u\left(t\right)\to 0$ , as $t\to +\infty$ . For the component $y_S\left(t\right)$ of $y\left(t\right)$ we have

$y_s\left(t\right)=S_s\left(t\right)y_{0s}-\rho {\displaystyle {\int }_0^t{S}_s}\left(t-\tau \right)\frac{\langle y\left(\tau \right),By\left(\tau \right)\rangle }{1+\left|\langle y\left(\tau \right),By\left(\tau \right)\rangle \right|}{B}_s{y}_s\left(\tau \right)\text{d}\tau$ (22)

Then for all $0\le t_0\le t$ , we have:

$y_s\left(t\right)=S_s\left(t-t_0\right)y_s\left(t_0\right)-\rho {\displaystyle \underset{t_0}{\overset{t}{\int }}{S}_s}\left(t-\tau \right)\frac{\langle y\left(\tau \right),By\left(\tau \right)\rangle }{1+\left|\langle y\left(\tau \right),By\left(\tau \right)\rangle \right|}{B}_s{y}_s\left(\tau \right)\text{d}\tau$ (23)

It follows from (13) that

$\Vert y_s\left(t\right)\Vert \le M_1{\text{e}}^{-\eta \left(t-t_0\right)}\Vert y_s\left(t_0\right)\Vert +\rho M_1\Vert B\Vert {\displaystyle {\int }_{t_0}^t{\text{e}}^{-\eta \left(t-\tau \right)}\Vert y_s\left(\tau \right)\Vert \text{d}\tau }$ (24)

From Gronwall inequality, we obtain:

$\Vert y_s\left(t\right)\Vert \le M_1{\text{e}}^{\left(\rho M_1\Vert B\Vert -\eta \right)\left(t-t_0\right)}\Vert y_s\left(t_0\right)\Vert ,\forall t\ge t_0$

Taking $\rho <\frac{\eta }{M_1\Vert B\Vert }$ , we deduce that $y_s\left(t\right)\to 0$ , as $t\to +\infty$ .

Hence $y\left(t\right)=y_u\left(t\right)+y_s\left(t\right)\to 0$ , as $t\to +\infty$ .

3. A Decay Rate Estimate of the Stabilized State

In what follows, we will study the strong stabilizability of the system (1) with the decay estimate (8).

Before we state our main result, the following lemmas will be needed (see [15] ).

3.1. Lemma

Let ${\left(s_k\right)}_{k\ge 0}$ be a sequence of positive real numbers satisfying

$s_{k+1}+C{s}_{k+1}^{2+\alpha }\le s_k,\forall k\ge 0$ (24)

where $C>0$ and $\alpha >-1$ are constants. Then there exists a positive constant $M_2$ (depending on $\alpha$ and $C$ ) such that

$s_k\le \frac{M_2}{{\left(k+1\right)}^{\alpha +1}},k\ge 0$ (25)

Let us now recall the following existing result (see [9] ).

3.2. Lemma

Let $A$ generate a contraction semigroup $S\left(t\right)$ on $H$ and let $B$ be linear operator from $H$ into itself. Then the system (1), controlled by (18) possesses a unique mild solution $y\left(t\right)\in H$ for each $y_0\in H$ which satisfies

${\displaystyle {\int }_0^T|}\langle S\left(s\right)y\left(t\right),BS\left(s\right)y\left(t\right)\rangle |ds=\mathcal{O}\left(\sqrt{{\displaystyle {\int }_t^{t+T}\frac{{\left|\langle y\left(s\right),By\left(s\right)\rangle \right|}^2}{1+\left|\langle By\left(s\right),y\left(s\right)\rangle \right|}\text{d}s}}\right),\text{as}t\to +\infty$ (26)

For almost all $T>0$ .

Our main result in this section is stated as follows:

3.3. Theorem

Let

1. $A$ generates a linear $C_0$ -semigroup ${\left(S\left(t\right)\right)}_{t\ge 0}$ such that ${\left(S_u\left(t\right)\right)}_{t\ge 0}$ is a semigroup of isometries and (13) holds,

2. $A$ allows the decomposition (10) of $H$ with $\text{dim}{H}_u<+\infty$ ,

3. $B\in \mathcal{L}\left(H\right)$ such that for all $y_u\in H_u$ , we have

$\langle B_u{S}_u\left(t\right)y_u,S_u\left(t\right)y_u\rangle =0,\forall t\ge 0\Rightarrow y_u=0$ (27)

Then the constrained feedback control law:

$p_{\rho ,u}\left(t\right)=-\rho \frac{\langle y_u\left(t\right),B_u{y}_u\left(t\right)\rangle }{1+\left|\langle y_u\left(t\right),B_u{y}_u\left(t\right)\rangle \right|},\rho >0,$ (28)

strongly stabilizes the system (1) with the explicit decay estimate (8).

Proof

Let us consider the system:

$\frac{\text{d}{y}_u\left(t\right)}{\text{d}t}=A_u{y}_u\left(t\right)-\rho \frac{\langle y_u\left(t\right),B_u{y}_u\left(t\right)\rangle }{1+\left|\langle y_u\left(t\right),B_u{y}_u\left(t\right)\rangle \right|}{B}_u{y}_u\left(t\right),y_u\left(0\right)=y_{0u}$ (29)

Multiplying the system (29) by $y_u\left(t\right)$ , integrating over $\Omega$ and using the fact that ${\left(S_u\left(t\right)\right)}_{t\ge 0}$ is a semigroup of isometries, we obtain:

$\frac{\text{d}{\Vert y_u\left(t\right)\Vert }^2}{\text{d}t}\le -2\rho \frac{{\left|\langle y_u\left(t\right),B_u{y}_u\left(t\right)\rangle \right|}^2}{1+\left|\langle y_u\left(t\right),B_u{y}_u\left(t\right)\rangle \right|}$ (30)

which proves that the real function $t\to \Vert y_u\left(t\right)\Vert$ is decreasing on ${ℝ}^{+}$ , and we have

$\Vert y_u\left(t\right)\Vert \le \Vert y_u\left(0\right)\Vert ,\forall t\ge 0$ (31)

Hence, the system (29) admits a unique mild solution defined for almost all $t\ge 0$ (see [9] ).

Integrating now the inequality (30) over the interval $\left[kT,\left(k+1\right)T\right]$ , for $k\in \mathbb{N}$ and $T>0$ , we get:

${\Vert y_u\left(\left(k+1\right)T\right)\Vert }^2-{\Vert y_u\left(kT\right)\Vert }^2\le -2\rho {\displaystyle {\int }_{kT}^{\left(k+1\right)T}\frac{{\left|\langle y_u\left(\tau \right),B_u{y}_u\left(\tau \right)\rangle \right|}^2}{1+\left|\langle y_u\left(\tau \right),B_u{y}_u\left(\tau \right)\rangle \right|}\text{d}\tau }$

using now the estimate (26), we deduce that

${\Vert y_u\left(\left(k+1\right)T\right)\Vert }^2-{\Vert y_u\left(kT\right)\Vert }^2\le -M{\left({\displaystyle {\int }_0^T\left|\langle S_u\left(\tau \right)y_u,B_u{S}_u\left(\tau \right)y_u\rangle \right|\text{d}\tau }\right)}^2$ (32)

for some $M>0$ . Using now the fact that $\text{dim}\left(H_u\right)<+\infty$ , then the assumption (27) is equivalent to

${\displaystyle {\int }_0^T\left|\langle B_u{S}_u\left(t\right)y_u,S_u\left(t\right)y_u\rangle \right|\text{d}t}\ge \delta {\Vert y_u\Vert }^2,\forall y_u\in H_u,\left(T,\delta >0\right)$ (33)

From (32) and (33) we have

${\Vert y_u\left(\left(k+1\right)T\right)\Vert }^2-{\Vert y_u\left(kT\right)\Vert }^2\le -M{\delta }^2{\Vert y_u\left(kT\right)\Vert }^4$

using the fact that the map $t\to \Vert y_u\left(t\right)\Vert$ is decreasing on ${ℝ}^{+}$ , we obtain:

${\Vert y_u\left(\left(k+1\right)T\right)\Vert }^2-{\Vert y_u\left(kT\right)\Vert }^2\le -M{\delta }^2{\Vert y_u\left(\left(k+1\right)T\right)\Vert }^4$

which implies that

${\Vert y_u\left(\left(k+1\right)T\right)\Vert }^2+C{\Vert y_u\left(\left(k+1\right)T\right)\Vert }^4\le {\Vert y_u\left(kT\right)\Vert }^2,C=M{\delta }^2$

Letting $s_k={\Vert y_u\left(kT\right)\Vert }^2$ , the last inequality can be written as

$s_{k+1}+C{s}_{k+1}^2\le s_k,\forall k\ge 0$

From Lemma.3.1 we have

$s_k\le \frac{M_2}{k+1}$

For $k=\left[\frac{t}{T}\right]$ ( $\left[\frac{t}{T}\right]$ designed the integer part of $\frac{t}{T}$ ), then we obtain $s_k\le \frac{M_3}{t}$ , $\left(M_3>0\right)$ , which gives

${\Vert y_u\left(t\right)\Vert }^2\le \frac{M_3}{t}$

Hence

$\Vert y_u\left(t\right)\Vert =\mathcal{O}\left(t^{-\frac{1}{2}}\right),\text{as}t\to +\infty$ (34)

For the component $y_s\left(t\right)$ , we shall show that $y_s\left(t\right)$ is defined for all $t\ge 0$ and exponentially converges to 0, as $t\to +\infty$ . The system (1) excited by the constrained feedback control (28) admits a unique mild solution defined for almost all $t$ in a maximal interval $\left[0,t_{\text{max}}\right[$ defined by

$y\left(t\right)=S\left(t\right)y_0+{\displaystyle {\int }_0^{t}S\left(t-\tau \right)p_{\rho ,u}\left(\tau \right)By\left(\tau \right)\text{d}\tau }$

Thus

$y_s\left(t\right)=S_s\left(t\right)y_{0s}+{\displaystyle {\int }_0^t{S}_s\left(t-\tau \right)p_{\rho ,u}\left(\tau \right)B_s{y}_s\left(\tau \right)\text{d}\tau },\forall t\in \left[0,t_{\text{max}}\right[$ (35)

It follows from (13) that

$\Vert y_s\left(t\right)\Vert \le M_1{\text{e}}^{-\eta t}\Vert y_{0s}\Vert +\rho M_1\Vert B\Vert {\displaystyle {\int }_0^t{\text{e}}^{-\eta \left(t-\tau \right)}\Vert y_s\left(\tau \right)\Vert \text{d}\tau }$

For almost all $t\in \left[0,t_{\text{max}}\right[$ . The Gronwall inequality then yields:

$\Vert y_s\left(t\right)\Vert \le M_1\Vert y_{0s}\Vert {\mathrm{e}}^{\left(\rho M_1\Vert B\Vert -\eta \right)t},\forall t\ge 0$ (36)

Taking $\rho <\frac{\eta }{M_1\Vert B\Vert }$ , it follows from (36) that $y_s\left(t\right)$ is bounded on $\left[0,t_{\text{max}}\right[$ so $t_{\text{max}}=+\infty$ , and therefore

(36) holds for all $t\ge 0$ . Hence, from (34) and (36), the solution of (1) satisfies the estimate (8). This completes the proof of Theorem 3.3.

3.4. Remark

1. Since the function $t\to \Vert y_u\left(t\right)\Vert$ decreasing in ${ℝ}^{+}$ , we have

$\exists t_0\ge 0;y_u\left(t_0\right)=0\iff y_u\left(t\right)=0,\forall t\ge t_0$

In this case, we have

$p_{\rho ,u}\left(t\right)=0,\forall t\ge t_0\Rightarrow y\left(t\right)=S_s\left(t-t_0\right)y_s\left(t_0\right),\forall t\ge t_0$

Hence, using (13) the system (1) is exponentially stable.

2. The constrained feedback control (28) depends only on the unstable part $y_u\left(t\right)$ and we have

$\left|p_{\rho ,u}\left(t\right)\right|<\rho ,\forall t\ge 0$

3. The constrained feedback control (28) satisfies

$\left|p_{\rho ,u}\left(t\right)\right|=\mathcal{O}\left(\frac{1}{t}\right),\text{as}t\to +\infty$

4. We note that (27) is weaker than (6). The converse is not true as we can see taking an orthonormal basis ${\left({\varphi }_n\right)}_{n\ge 1}$ of $H$ , $Az=-\langle z,{\varphi }_1\rangle {\varphi }_1$ and $Bz={\displaystyle \underset{n=1}{\overset{+\infty }{\sum }}\frac{\langle z,{\varphi }_n\rangle }{n^2}{\varphi }_n}$ .

5. In the case $\text{dim}{H}_u=+\infty$ and $B$ is nonlinear and locally Lipschitz, such that $B\left(0\right)=0$ , then using the same techniques as in [9] , we can obtain the result of Theorem 3.3, if the estimate (7) is changed to (33).

4. Robustness

In this section, we study the robustness of the controls (18) and (28), under a class of perturbations of the system (1).

4.1. Strong Robustness

In this part, we consider the strong robustness of the feedback (18). Then, we will show that the stability property of the system (1) remains invariant under a certain class of bounded perturbations.

Let us consider the following perturbed system

$\frac{\text{d}y\left(t\right)}{\text{d}t}=Ay\left(t\right)+p\left(t\right)By\left(t\right)+\xi \left(y\left(t\right)\right),y\left(0\right)=y_0$ (37)

where the linear bounded operator $\xi ={\xi }_u+{\xi }_s$ is such that the system (37) is decomposed into two following subsystems:

$\frac{\text{d}{y}_u\left(t\right)}{\text{d}t}=A_u{y}_u\left(t\right)+p\left(t\right)B_u{y}_u\left(t\right)+{\xi }_u\left(y_u\left(t\right)\right),y_u\left(0\right)=y_{0u}\in H_u$ (38)

$\frac{\text{d}{y}_s\left(t\right)}{\text{d}t}=A_s{y}_s\left(t\right)+p\left(t\right)B_s{y}_s\left(t\right)+{\xi }_s\left(y_u\left(t\right)\right),y_s\left(0\right)=y_{0s}\in H_s$ (39)

The following main result concerns the strong stability of the system (37).

4.2. Proposition

Let

1. $A$ generates a linear $C_0$ -contraction semigroup ${\left(S\left(t\right)\right)}_{t\ge 0}$ on $H$ such that (13) holds,

2. The operator $B$ is compact such that (6) holds,

3. The linear operator $\xi$ is compact and satisfying $\Vert {\xi }_s\Vert <\frac{\eta }{M_1}$ and $\langle \xi \left(y\right),y\rangle \le 0$ , $\forall y\in H$ .

Then the system (37) is strongly stabilizable.

Proof

First, let us note that 0 remains an equilibrium state of the perturbed system (37), which can be written in the form:

$\frac{\text{d}y\left(t\right)}{\text{d}t}=Ay\left(t\right)+g\left(y\left(t\right)\right),y\left(0\right)=y_0$ (40)

where $g=f+\xi$ and $f\left(y\right)=\{\begin{array}{ll}-\frac{\rho \langle By,y\rangle }{1+\left|\langle By,y\rangle \right|}By,\hfill & y\ne 0;\hfill \\ 0,\hfill & y=0.\hfill \end{array}$

Since $f$ and $\xi$ are locally Lipschitz, so is $g$ . Also $g$ is dissipative: $\langle \xi \left(y\right),y\rangle \le 0$ , $\forall y\in H$ .

The assumption $\langle \xi \left(y\right),y\rangle \le 0$ , $\forall y\in H$ , together with (6) guarantees the following implication

$\langle g\left(S\left(t\right)y\right),S\left(t\right)y\rangle =0\Rightarrow y=0$

Then the weak stability of the perturbed system (37) follows from Theorem 2.4 of Ball [3] , and since $\text{dim}{H}_u<+\infty$ , we have $y_u\left(t\right)\to 0$ , as $t\to +\infty$ . For the component $y_s\left(t\right)$ of the solution $y\left(t\right)$ of the system (37), and for all $0\le t_0\le t$ , we have:

$y_s\left(t\right)=S_s\left(t-t_0\right)y_s\left(t_0\right)-\rho {\displaystyle \underset{t_0}{\overset{t}{\int }}{S}_s\left(t-\tau \right)\frac{\langle y\left(\tau \right),By\left(\tau \right)\rangle }{1+\left|\langle y\left(\tau \right),By\left(\tau \right)\rangle \right|}{B}_s{y}_s\left(\tau \right)\text{d}\tau }+{\displaystyle \underset{t_0}{\overset{t}{\int }}{S}_s\left(t-\tau \right){\xi }_s\left(y_s\left(\tau \right)\right)\text{d}\tau }$ (41)

It follows from (13) that

$\Vert y_s\left(t\right)\Vert \le M_1{\text{e}}^{-\eta \left(t-t_0\right)}\Vert y_s\left(t_0\right)\Vert +\left(\rho M_1\Vert B\Vert +M_1\Vert {\xi }_s\Vert \right){\displaystyle {\int }_{t_0}^t{\text{e}}^{-\eta \left(t-\tau \right)}\Vert y_s\left(\tau \right)\Vert \text{d}\tau }$ (42)

From Gronwall inequality, we obtain:

$\Vert y_s\left(t\right)\Vert \le M_1{\text{e}}^{\left(\rho M_1\Vert B\Vert +M_1\Vert {\xi }_s\Vert -\eta \right)\left(t-t_0\right)}\Vert y_s\left(t_0\right)\Vert ,\forall t\ge t_0$

Taking $\rho <\frac{\eta -M_1\Vert {\xi }_s\Vert }{M_1\Vert B\Vert }$ , we obtain $y_s\left(t\right)\to 0$ , as $t\to +\infty$ . Hence, the solution $y\left(t\right)$ of the system (37) strongly converges to $0,$ as $t\to +\infty$ .

4.3. A Polynomial Decay Estimate for the Perturbed System

Our second main result in this section is stated as follows:

4.4. Proposition

Let

1. $A$ generate a linear $C_0$ -semigroup ${\left(S\left(t\right)\right)}_{t\ge 0}$ such that ${\left(S_u\left(t\right)\right)}_{t\ge 0}$ is a semigroup of isometries and (13) holds,

2. $A$ allows the decomposition (10) of $H$ with $\text{dim}{H}_u<+\infty$ ,

3. $B\in \mathcal{L}\left(H\right)$ satisfies (27),

4. $\Vert {\xi }_s\Vert \le \frac{\eta }{M_1}$ and $\langle {\xi }_u\left(y\right),y\rangle \le 0$ , for all $y\in H_u.$

Then the constrained feedback control (28) strongly stabilizes the system (37) with the explicit decay estimate (8).

Proof

Let us consider the system:

$\frac{\text{d}{y}_u\left(t\right)}{\text{d}t}=A_u{y}_u\left(t\right)-\rho \frac{\langle y_u\left(t\right),B_u{y}_u\left(t\right)\rangle }{1+\left|\langle y_u\left(t\right),B_u{y}_u\left(t\right)\rangle \right|}{B}_u{y}_u\left(t\right)+{\xi }_u\left(y_u\left(t\right)\right),y_u\left(0\right)=y_{0u}$ (43)

Multiplying the system (43) by $y_u\left(t\right)$ and integrating over $\Omega$ and using the fact that ${\left(S_u\left(t\right)\right)}_{t\ge 0}$ is a semigroup of isometries and the hypothesis: $\langle {\xi }_u\left(y_u\right),y_u\rangle \le 0$ , $\forall y_u\in H_u$ , we obtain:

$\frac{\text{d}{\Vert y_u\left(t\right)\Vert }^2}{\text{d}t}\le -2\rho \frac{{\left|\langle y_u\left(t\right),B_u{y}_u\left(t\right)\rangle \right|}^2}{1+\left|\langle y_u\left(t\right),B_u{y}_u\left(t\right)\rangle \right|}$ (44)

which gives $\Vert y_u\left(t\right)\Vert \le \Vert y_u\left(0\right)\Vert$ , $\forall t\ge 0$ . Then, the system (43) admits a unique global mild solution $y\left(t\right)$ de-

fined for almost all $t\ge 0$ . By the same argument as in the proof of Theorem 3.3, the solution $y\left(t\right)$ of the system (37) satisfies:

$\Vert y\left(t\right)\Vert =\mathcal{O}\left(t^{-\frac{1}{2}}\right),\text{as}t\to +\infty$

which completes the proof of Proposition 4.2.

5. Application and Simulations

5.1. An Application

In this part, we will give an illustrating example of the established results.

Example

Let us consider the following 1-d bilinear heat equation:

$\{\begin{array}{ll}\frac{\partial y\left(x,t\right)}{\partial t}=\frac{{\partial }^{2}y\left(x,t\right)}{\partial x^2}+p\left(t\right)By\left(t\right),\hfill & x\in \left(0,1\right),t>0,\hfill \\ \frac{\partial y\left(0,t\right)}{\partial x}=\frac{\partial y\left(1,t\right)}{\partial x}=0,\hfill & \forall t>0,\hfill \end{array}$ (45)

where $y\left(t\right)$ is the temperature profile at time $t$ . We suppose that the system is controlled via the flow of a liquid $p\left(t\right)$ in an adequate metallic pipeline. Here we take the state space $H=L^2\left(0,1\right)$ and the operator $A$ is defined by $Ay=\frac{{\partial }^{2}y}{\partial x^2}$ , with $\mathcal{D}\left(A\right)=\left\{y\in H^2\left(0,1\right)|\frac{\partial y\left(0,t\right)}{\partial x}=\frac{\partial y\left(1,t\right)}{\partial x}=0\right\}$ . The domain of $A$ gives the ho-

mogeneous Neumann boundary condition imposed at the ends of the bar which require specifying how the heat flows out of the bar and means that both ends are insulated. The spectrum of $A$ is given by the simple eigen-

values ${\lambda }_j=-{\text{π}}^2{\left(j-1\right)}^2$ , $j\in {\mathbb{N}}^{\ast }$ and eigenfunctions ${\phi }_1\left(x\right)=1$ and ${\phi }_j\left(x\right)=\sqrt{2}\text{cos}\left(\left(j-1\right)\text{π}x\right)$ for all $j\ge 2$ . Then the subspace $H_u$ is the one-dimensional space spanned by the eigenfunction ${\phi }_1$ , and we have $S_u\left(t\right)y_u=\langle y_u,{\phi }_1\rangle {\phi }_1$ so $S_u\left(t\right)=I_{H_u}$ (the identity) and hence ${\left(S_u\left(t\right)\right)}_{t\ge 0}$ , is a semigroup of isometries. The operator of control $B$, is defined by: $By={\displaystyle \underset{j=1}{\overset{+\infty }{\sum }}{\alpha }_j}\langle y,{\phi }_j\rangle {\phi }_j$ , ${\alpha }_j\ge 0$ , $\forall j\ge 1$ , such that ${\displaystyle \underset{j=1}{\overset{+\infty }{\sum }}{\alpha }_j^2}<\infty$ (see [16] ). From the relation: $\langle B_u{S}_u\left(t\right)y_u,S_u\left(t\right)y_u\rangle ={\alpha }_1{\left|\langle y_u,{\phi }_1\rangle \right|}^2$ , we can see that (27) holds if ${\alpha }_1>0$ . To examine the estimate (8), remarking for the scalar functions $y_j\left(t\right)=\langle y\left(t\right),{\phi }_j\rangle$ , $\forall j\ge 1$ we have

$\frac{\text{d}{y}_1\left(t\right)}{\text{d}t}=\frac{\text{d}{y}_u\left(t\right)}{\text{d}t}=-\frac{\rho {\alpha }_1{y}_u^3\left(t\right)}{1+{\alpha }_1{y}_u^2\left(t\right)},\forall t\ge 0$ (46)

which implies that

$\frac{\text{d}}{\text{d}t}{y}_u^2\left(t\right)=-\frac{2\rho {\alpha }_1{y}_u^4\left(t\right)}{1+{\alpha }_1{y}_u^2\left(t\right)}$ (47)

Letting $x\left(t\right)=y_u^2\left(t\right)$ , we obtain

$x^{\prime }\left(t\right)x^{-2}\left(t\right)=-\frac{2\rho {\alpha }_1}{1+{\alpha }_{1}x\left(t\right)},\forall t>0$

Integrating now the last equality from 0 to $t$ , we get

$x^{-1}\left(t\right)-x^{-1}\left(0\right)={\displaystyle {\int }_0^t\frac{2\rho {\alpha }_1}{1+{\alpha }_{1}x\left(\tau \right)}\text{d}\tau }$

from (47), we deduce that the nonnegative scalar map $t\to x\left(t\right)$ is decreasing for all $t\ge 0$ , and we have

$x^{-1}\left(t\right)\ge \frac{2\rho {\alpha }_{1}t}{1+{\alpha }_{1}x\left(0\right)},\forall t>0$

which means that

$\left|y_u\left(t\right)\right|\le \sqrt{\frac{1+{\alpha }_1{y}_{0u}^2}{2\rho {\alpha }_{1}t}},\forall t>0,y_{0u}\ne 0$

Then

$\left|y_u\left(t\right)\right|=\mathcal{O}\left(\frac{1}{\sqrt{t}}\right),\text{as}t\to +\infty$

Furthermore, the control in this case is defined by

$p_{\rho ,u}\left(t\right)=-\frac{\rho {\alpha }_1{y}_u^2\left(t\right)}{1+{\alpha }_1{y}_u^2\left(t\right)},\forall t\ge 0$ (48)

For $j\ge 2$ , the functions $y_j\left(t\right)$ are characterized by $y_j\left(0\right)=\langle y\left(0\right),{\phi }_j\rangle$ , $\forall j\ge 2$ and satisfy

$\frac{\partial }{\partial t}{y}_j\left(t\right)=\left({\lambda }_j-\frac{\rho {\alpha }_1{\alpha }_j{y}_1^2\left(t\right)}{1+{\alpha }_1{y}_1^2\left(t\right)}\right)y_j\left(t\right),t\ge 0$

which implies that

$\left|y_j\left(t\right)\right|\le {\text{e}}^{-{\text{π}}^2{\left(j-1\right)}^{2}t}\left|y_j\left(0\right)\right|,j\ge 2$

Then

$\Vert y_s\left(t\right)\Vert \le {\mathrm{e}}^{-{\text{π}}^{2}t}\Vert y_s\left(0\right)\Vert ,\forall t\ge 0$

Hence, the system (45) is strongly stable with the decay rate estimate (8).

Let us reconsider the above example with the perturbation $\xi$ defined by:

$\xi \left(y\right)=-\beta y,\beta \in \left(0,\frac{{\text{π}}^{\text{2}}}{2}\right),\forall y\in H$

It is clear that the function $\xi$ satisfies the conditions of Proposition 4.2.1. Then the perturbed closed-loop system remains stable, i.e., the control (48) still stabilizes the perturbed system i.e., the control (48) is strongly robust.

5.2. Simulations

In this part, taking in the system (45), the operator $B=I$ and $y_0\left(x\right)=x+0.3$ .

Then we obtain the results shown in Figures 1-5.

Figure 1. The norm of the free state.

Figure 2. The free state.

Figure 3. The norm of the stabilized state.

Figure 4. The stabilized state.

Figure 5. The stabilizing control.

6. Conclusion

In this work, we have considered the problem of strong stabilization with polynomial decay rate of the stabilized state for bilinear parabolic systems that can be decomposed in the stable and unstable parts (15) and (16) under a weaker condition (27). We have also considered the problem of using a stabilizing feedback control for the unstable part (15) only that can make the whole system (1) stable. Various questions remain open. This is the case

of stabilization for nonlinear systems. Finally, we have studied the robustness problem of the stabilizing controls with respect to a class of perturbations, but a confrontation to more realistic situations remain done. This leads us to consider the stabilization problem for stochastic bilinear system.

Conflicts of Interest

The authors declare no conflicts of interest.

References

[1]	Mohler, R.R. (1973) Bilinear Control Process. Academic Press, New York.
[2]	Tsouli, A. and Boutoulout, A. (2014) Controllability of the Parabolic System via Bilinear Control. Journal of Dynamical and Control Systems, 20, in press. http://dx.doi.org/10.1007/s10883-014-9247-2
[3]	Ball, J. and Slemrod, M. (1979) Feedback Stabilization of Distributed Semilinear Control Systems. Applied Mathematics and Optimization, 5, 169-179. http://dx.doi.org/10.1007/BF01442552
[4]	Jurdjevic, V. and Quinn, J.P. (1978) Controllability and Stability. Journal of Deferential Equations, 28, 381-389. http://dx.doi.org/10.1016/0022-0396(78)90135-3
[5]	Quinn, J.P. (1980) Stabilization of Bilinear Systems by Quadratic Feedback Control. Journal of Mathematical Analysis and Applications, 75, 66-80. http://dx.doi.org/10.1016/0022-247X(80)90306-6
[6]	Ouzahra, M. (2008) Strong Stabilization with Decay Estimate of Semilinear Systems. Systems and Control Letters, 57, 813-815. http://dx.doi.org/10.1016/j.sysconle.2008.03.009
[7]	Ouzahra, M., Tsouli, A. and Boutoulout, A. (2012) Stabilization and Polynomial Decay Estimate for Distributed Semilinear Systems. International Journal of Control, 85, 451-456. http://dx.doi.org/10.1080/00207179.2012.656144
[8]	Bounit, H. and Hammouri, H. (1999) Feedback Stabilization for a Class of Distributed Semilinear Control Systems. Nonlinear Analysis, 37, 953-969. http://dx.doi.org/10.1016/S0362-546X(97)00577-4
[9]	Tsouli, A., Ouzahra, M. and Boutoulout, A. (2014) A Decay Estimate for Constrained Semilinear Systems. Information Sciences Letters, 3, 77-83. http://dx.doi.org/10.12785/isl/030206
[10]	Chen, M.S. (1998) Exponential Stabilization of a Constrained Bilinear System. Automatica, 34, 989-992. http://dx.doi.org/10.1016/S0005-1098(98)00037-5
[11]	Kato, T. (1980) Perturbation Theory for Linear Operators. Springer, New York.
[12]	Ouzahra, M. (2011) Feedback Stabilization of Parabolic Systems with Bilinear Controls. Electronic Journal of Differential Equations, 38, 1-10.
[13]	Triggiani, R. (1975) On the Stabilizability Problem in Banach Space. Journal of Mathematical Analysis and Applications, 52, 383-403. http://dx.doi.org/10.1016/0022-247X(75)90067-0
[14]	Pazy, A. (1983) Semi-Groups of Linear Operators and Applications to Partial Deferential Equations. Springer Verlag, New York. http://dx.doi.org/10.1007/978-1-4612-5561-1
[15]	Ammari, K. and Tucsnak, M. (2001) Stabilization of Second Order Evolution Equations by a Class of Unbounded Feedbacks. ESAIM: Control, Optimisation and Calculus of Variations, 6, 361-386. http://dx.doi.org/10.1051/cocv:2001114
[16]	Ouzahra, M. (2010) Exponential and Weak Stabilization of Constrained Bilinear Systems. SIAM Journal on Control and Optimization, 48, 3962-3974. http://dx.doi.org/10.1137/080739161