<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">ICA</journal-id><journal-title-group><journal-title>Intelligent Control and Automation</journal-title></journal-title-group><issn pub-type="epub">2153-0653</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/ica.2015.62011</article-id><article-id pub-id-type="publisher-id">ICA-55169</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Computer Science&amp;Communications</subject></subj-group></article-categories><title-group><article-title>
 
 
  Constrained Feedback Stabilization for Bilinear Parabolic Systems
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Azzeddine</surname><given-names>Tsouli</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Ali</surname><given-names>Boutoulout</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>TSI Team, MACS Laboratory, Faculty of Sciences, Moulay Ismail University, Meknes, Morocco</addr-line></aff><pub-date pub-type="epub"><day>30</day><month>03</month><year>2015</year></pub-date><volume>06</volume><issue>02</issue><fpage>103</fpage><lpage>115</lpage><history><date date-type="received"><day>January</day>	<month>2</month>	<year>2015</year></date><date date-type="rev-recd"><day>March</day>	<month>27</month>	<year>2015</year>	</date><date date-type="accepted"><day>March</day>	<month>30</month>	<year>2015</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  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.
 
</p></abstract><kwd-group><kwd>Bilinear Parabolic Systems</kwd><kwd> Constrained Feedback Stabilisation</kwd><kwd> Decay Estimate</kwd><kwd> Robustness</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>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 [<xref ref-type="bibr" rid="scirp.55169-ref1">1</xref>] [<xref ref-type="bibr" rid="scirp.55169-ref2">2</xref>] ). 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:</p><p>d y ( t ) d t = A y ( t ) + p ( t ) B y ( t ) , y ( 0 ) = y 0 (1)</p><p>on a real Hilbert space H with inner product 〈 . , . 〉 ; and corresponding norm ‖ . ‖ , where the linear operator A generates a contraction semigroup ( S ( t ) ) t ≥ 0 on H and B ∈ L ( H ) . While the real valued function p ( . ) ∈ L 2 ( 0 , + ∞ ; ℝ ) represents a control. A function y ( . ) ∈ C ( 0 , t 0 ; H ) , t 0 &gt; 0 , is a mild solution of the system (1) if and only if the solution y ( t ) of the system (1) satisfies the variation of parameters formula:</p><p>y ( t ) = S ( t ) y 0 + ∫ 0 t S ( t − τ ) p ( τ ) B y ( τ ) d τ , t ≥ 0 (2)</p><p>(see [<xref ref-type="bibr" rid="scirp.55169-ref3">3</xref>] ). By choosing an adequate feedback control p ( t ) in such a way, the corresponding solution y ( t ) of the system (1) converges to zero when t → + ∞ , 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 [<xref ref-type="bibr" rid="scirp.55169-ref4">4</xref>] , under the condition:</p><p>span { A z , a d 0 ( A , B ) z , a d 1 ( A , B ) z , ⋯ , a d k ( A , B ) z , ⋯ } = ℝ n , ∀ z ∈ ℝ n − { 0 } (3)</p><p>where a d k ( A , B ) is defined recursively as a d 0 ( A , B ) = B , a d 1 ( A , B ) = A B − B A and a d k + 1 ( A , B ) = a d 1 ( A , a d k ( A , B ) ) , ∀ k ∈ ℕ .</p><p>Using the following assumptions:</p><p>〈 B exp ( t A ) y , exp ( t A ) y 〉 = 0 , , ∀ t ≥ 0 ⇒ y = 0 (4)</p><p>the problem of stabilization has been studied in [<xref ref-type="bibr" rid="scirp.55169-ref5">5</xref>] . In [<xref ref-type="bibr" rid="scirp.55169-ref3">3</xref>] , when the linear operator B is compact and ( S ( t ) ) t ≥ 0 is a contraction semigroup, then using the quadratic feedback control</p><p>p 0 ( t ) = − 〈 y ( t ) , B y ( t ) 〉 (5)</p><p>a weak stabilization result is obtained under the weak observability condition:</p><p>〈 B S ( t ) y , S ( t ) y 〉 = 0 , ∀ t ≥ 0 ⇒ y = 0 (6)</p><p>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 [<xref ref-type="bibr" rid="scirp.55169-ref3">3</xref>] ). Under the exact observability assumption</p><p>∫ 0 T | 〈 B S ( t ) y , S ( t ) y 〉 | d t ≥ δ ‖ y ‖ 2 , ∀ y ∈ H , ( T , δ &gt; 0 ) (7)</p><p>The strong stabilization result with the following decay estimate</p><p>‖ y ( t ) ‖ = O ( 1 t ) , as t → + ∞ (8)</p><p>i.e. ‖ y ( t ) ‖ ≤ M t , M &gt; 0 for t large enough, has been obtained using the quadratic feedback control (5) (see [<xref ref-type="bibr" rid="scirp.55169-ref6">6</xref>] ). However, in this way the convergence of the resulting closed loop state is not better than (8). In [<xref ref-type="bibr" rid="scirp.55169-ref7">7</xref>] the rational decay rates are established i.e. using the following feedback control:</p><p>p r ( t ) = − 〈 y ( t ) , B y ( t ) 〉 ‖ y ( t ) ‖ r , r ∈ ( − ∞ , 2 )</p><p>It has been shown in [<xref ref-type="bibr" rid="scirp.55169-ref8">8</xref>] 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><p>p ( t ) = − 〈 B y ( t ) , y ( t ) 〉 1 + | 〈 B y ( t ) , y ( t ) 〉 | (9)</p><p>strongly stabilizes the system (1), provided that the assumption (6) holds. It has been established in [<xref ref-type="bibr" rid="scirp.55169-ref7">7</xref>] that, if</p><p>the linear operator A generates a contraction semigroup ( S ( t ) ) t ≥ 0 in H , then the system (1) is strongly</p><p>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 [<xref ref-type="bibr" rid="scirp.55169-ref3">3</xref>] [<xref ref-type="bibr" rid="scirp.55169-ref8">8</xref>] [<xref ref-type="bibr" rid="scirp.55169-ref9">9</xref>] . 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.</p></sec><sec id="s2"><title>2. Stabilization Results</title><p>Let us now recall the following definition concerning the asymptotic behavior of the system (1).</p><sec id="s2_1"><title>2.1. Definition</title><p>The system (1) is weakly (resp. strongly) stabilizable if there exists a feedback control p ( t ) = f ( y ( t ) ) , t ≥ 0 f : H → K : = ℝ , ℂ such that the corresponding mild solution y ( t ) of the system (1) satisfies the properties:</p><p>1. For each initial state y 0 of the system (1) there exists a unique mild solution defined for all t ∈ ℝ + of the system (1),</p><p>2. { 0 } is an equilibrium state of the system (1),</p><p>3. y ( t ) → 0 , weakly (resp. strongly), as t → + ∞ , for all y 0 ∈ H .</p><p>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 [<xref ref-type="bibr" rid="scirp.55169-ref10">10</xref>] - [<xref ref-type="bibr" rid="scirp.55169-ref12">12</xref>] , it has been shown that if the spectrum σ ( A ) of A can be decomposed</p><p>into σ u ( A ) = { λ : R e ( λ ) ≥ − η , η &gt; 0 } and σ s ( A ) = { λ : R e ( λ ) &lt; − η } , then the state space H can be decomposed according to</p><p>H = H u ⊕ H s (10)</p><p>where H u = P u H = vect { φ j , 1 ≤ j ≤ N } , H s = P s H = vect { φ j , j &gt; N } , P u is given by</p><p>P u = 1 2 π i ∫ C ( λ I − A ) − 1 d λ (11)</p><p>where C is a curve surrounding σ u ( A ) , P s = I − P u and for all j ≥ 1 , φ j is the eigenvector associated to</p><p>the eigenvalue λ 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 ( t ) ∈ H , we set y u = P u y and y s = P s y . For linear systems, it</p><p>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:</p><p>lim t → + ∞ ln ( ‖ S s ( t ) ‖ ) t = supRe ( σ ( A s ) ) (12)</p><p>Which is equivalent to:</p><p>‖ S s ( t ) ‖ ≤ M 1 exp ( − η t ) , ∀ t ≥ 0 ( forsome M 1 &gt; 0 ) (13)</p><p>where ( S s ( t ) ) t ≥ 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 [<xref ref-type="bibr" rid="scirp.55169-ref13">13</xref>] ). In the sequel, we suppose that the operator B satisfies</p><p>B H u ⊂ H u and B H s ⊂ H s (14)</p><p>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:</p><p>d y u ( t ) d t = A u y u ( t ) + p ( t ) B u y u ( t ) , y u ( 0 ) = y 0 u ∈ H u (15)</p><p>d y s ( t ) d t = A s y s ( t ) + p ( t ) B s y s ( t ) , y s ( 0 ) = y 0 s ∈ H s (16)</p><p>in the state spaces H u and H S respectively, and B = B u ⊕ B s . It has been proved that stabilizing a linear system turns out to stabilizing its unstable part (see [<xref ref-type="bibr" rid="scirp.55169-ref13">13</xref>] ).</p></sec><sec id="s2_2"><title>2.2. Remark</title><p>For finite-dimensional systems, the conditions (6) and (7) are equivalent (see [<xref ref-type="bibr" rid="scirp.55169-ref5">5</xref>] [<xref ref-type="bibr" rid="scirp.55169-ref10">10</xref>] ). However, in infinite- dimensional case, and if B is compact, then the condition (7) is impossible. Indeed, if ( φ j ) j ≥ 1 is an orthonormal basis of H , then applying (7) for y = φ j and using the fact that φ j → 0 , weakly as j → + ∞ , we obtain the contradiction: δ = 0.</p><p>The following result concerns the strong stabilization of the system (1).</p></sec><sec id="s2_3"><title>2.3. Theorem</title><p>Let</p><p>1. A generates a linear C 0 -contraction semigroup ( S ( t ) ) t ≥ 0 on H ,</p><p>2. A allows the decomposition (10) of H with dim H u &lt; + ∞ such that (13) holds,</p><p>3. B be compact such that</p><p>〈 B S ( t ) y , S ( t ) y 〉 = 0 , ∀ t ≥ 0 ⇒ y = 0 (17)</p><p>Then, the constrained feedback control law:</p><p>p ρ ( t ) = − ρ 〈 y ( t ) , B y ( t ) 〉 1 + | 〈 y ( t ) , B y ( t ) 〉 | , ρ &gt; 0 (18)</p><p>strongly stabilizes the system (1).</p><p>Proof</p><p>The system (1) controlled by (18) possesses a unique mild solution y ( . ) defined on a maximal interval [ 0 , t max [ and given by the variation of constants formula</p><p>y ( t ) = S ( t ) y 0 − ρ ∫ 0 t S ( t − τ ) 〈 y ( τ ) , B y ( τ ) 〉 1 + | 〈 y ( τ ) , B y ( τ ) 〉 | B y ( τ ) d τ (19)</p><p>corresponds to (18) (see [<xref ref-type="bibr" rid="scirp.55169-ref9">9</xref>] ).</p><p>Since ( S ( t ) ) t ≥ 0 is a contraction semigroup, we get:</p><p>d ‖ y ( t ) ‖ 2 d t ≤ − 2 ρ | 〈 y ( t ) , B y ( t ) 〉 | 2 1 + | 〈 y ( t ) , B y ( t ) 〉 | , ∀ y 0 ∈ D ( A ) (20)</p><p>It follows from (20) that</p><p>‖ y ( t ) ‖ ≤ ‖ y 0 ‖ , ∀ t ≥ 0 (21)</p><p>From (19) and using the fact that ( S ( t ) ) t ≥ 0 is a contraction semigroup and the Gronwall inequality, we deduce that the map y 0 → y ( t ) is continuous from H to H . Then (21) holds for all y 0 ∈ H by density argument, and hence t max = + ∞ (see [<xref ref-type="bibr" rid="scirp.55169-ref14">14</xref>] ). Now, let us show that y ( t ) → 0 , weakly as t → + ∞ . Let t n → + ∞</p><p>such that y ( t n ) weakly converges in H and let z ∈ H such that y ( t n ) → z , weakly as n → ∞ . (The ex- istence of the sequence ( t n ) n ∈ ℕ and z are ensured by (21) and by the fact that space H is reflexive.) Taking y ( t n ) as initial state in (19) and using superposition property of the solution, and via the dominated conver-</p><p>gence theorem, we obtain 〈 B S ( t ) z , S ( t ) z 〉 = 0 , ∀ t ≥ 0 . It follows from (17) that z = 0 . Hence y ( t ) → 0 , weakly as t → + ∞ , and since dim H u &lt; + ∞ , we have y u ( t ) → 0 , as t → + ∞ . For the component y S ( t ) of y ( t ) we have</p><p>y s ( t ) = S s ( t ) y 0 s − ρ ∫ 0 t S s ( t − τ ) 〈 y ( τ ) , B y ( τ ) 〉 1 + | 〈 y ( τ ) , B y ( τ ) 〉 | B s y s ( τ ) d τ (22)</p><p>Then for all 0 ≤ t 0 ≤ t , we have:</p><p>y s ( t ) = S s ( t − t 0 ) y s ( t 0 ) − ρ ∫ t 0 t S s ( t − τ ) 〈 y ( τ ) , B y ( τ ) 〉 1 + | 〈 y ( τ ) , B y ( τ ) 〉 | B s y s ( τ ) d τ (23)</p><p>It follows from (13) that</p><p>‖ y s ( t ) ‖ ≤ M 1 e − η ( t − t 0 ) ‖ y s ( t 0 ) ‖ + ρ M 1 ‖ B ‖ ∫ t 0 t e − η ( t − τ ) ‖ y s ( τ ) ‖ d τ (24)</p><p>From Gronwall inequality, we obtain:</p><p>‖ y s ( t ) ‖ ≤ M 1 e ( ρ M 1 ‖ B ‖ − η ) ( t − t 0 ) ‖ y s ( t 0 ) ‖ , ∀ t ≥ t 0</p><p>Taking ρ &lt; η M 1 ‖ B ‖ , we deduce that y s ( t ) → 0 , as t → + ∞ .</p><p>Hence y ( t ) = y u ( t ) + y s ( t ) → 0 , as t → + ∞ .</p></sec></sec><sec id="s3"><title>3. A Decay Rate Estimate of the Stabilized State</title><p>In what follows, we will study the strong stabilizability of the system (1) with the decay estimate (8).</p><p>Before we state our main result, the following lemmas will be needed (see [<xref ref-type="bibr" rid="scirp.55169-ref15">15</xref>] ).</p><sec id="s3_1"><title>3.1. Lemma</title><p>Let ( s k ) k ≥ 0 be a sequence of positive real numbers satisfying</p><p>s k + 1 + C s k + 1 2 + α ≤ s k , ∀ k ≥ 0 (24)</p><p>where C &gt; 0 and α &gt; − 1 are constants. Then there exists a positive constant M 2 (depending on α and C ) such that</p><p>s k ≤ M 2 ( k + 1 ) α + 1 , k ≥ 0 (25)</p><p>Let us now recall the following existing result (see [<xref ref-type="bibr" rid="scirp.55169-ref9">9</xref>] ).</p></sec><sec id="s3_2"><title>3.2. Lemma</title><p>Let A generate a contraction semigroup S ( t ) 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 ( t ) ∈ H for each y 0 ∈ H which satisfies</p><p>∫ 0 T | 〈 S ( s ) y ( t ) , B S ( s ) y ( t ) 〉 | d s = O ( ∫ t t + T | 〈 y ( s ) , B y ( s ) 〉 | 2 1 + | 〈 B y ( s ) , y ( s ) 〉 | d s ) , as t → + ∞ (26)</p><p>For almost all T &gt; 0 .</p><p>Our main result in this section is stated as follows:</p></sec><sec id="s3_3"><title>3.3. Theorem</title><p>Let</p><p>1. A generates a linear C 0 -semigroup ( S ( t ) ) t ≥ 0 such that ( S u ( t ) ) t ≥ 0 is a semigroup of isometries and (13) holds,</p><p>2. A allows the decomposition (10) of H with dim H u &lt; + ∞ ,</p><p>3. B ∈ L ( H ) such that for all y u ∈ H u , we have</p><p>〈 B u S u ( t ) y u , S u ( t ) y u 〉 = 0 , ∀ t ≥ 0 ⇒ y u = 0 (27)</p><p>Then the constrained feedback control law:</p><p>p ρ , u ( t ) = − ρ 〈 y u ( t ) , B u y u ( t ) 〉 1 + | 〈 y u ( t ) , B u y u ( t ) 〉 | , ρ &gt; 0 , (28)</p><p>strongly stabilizes the system (1) with the explicit decay estimate (8).</p><p>Proof</p><p>Let us consider the system:</p><p>d y u ( t ) d t = A u y u ( t ) − ρ 〈 y u ( t ) , B u y u ( t ) 〉 1 + | 〈 y u ( t ) , B u y u ( t ) 〉 | B u y u ( t ) , y u ( 0 ) = y 0 u (29)</p><p>Multiplying the system (29) by y u ( t ) , integrating over Ω and using the fact that ( S u ( t ) ) t ≥ 0 is a semigroup of isometries, we obtain:</p><p>d ‖ y u ( t ) ‖ 2 d t ≤ − 2 ρ | 〈 y u ( t ) , B u y u ( t ) 〉 | 2 1 + | 〈 y u ( t ) , B u y u ( t ) 〉 | (30)</p><p>which proves that the real function t → ‖ y u ( t ) ‖ is decreasing on ℝ + , and we have</p><p>‖ y u ( t ) ‖ ≤ ‖ y u ( 0 ) ‖ , ∀ t ≥ 0 (31)</p><p>Hence, the system (29) admits a unique mild solution defined for almost all t ≥ 0 (see [<xref ref-type="bibr" rid="scirp.55169-ref9">9</xref>] ).</p><p>Integrating now the inequality (30) over the interval [ k T , ( k + 1 ) T ] , for k ∈ ℕ and T &gt; 0 , we get:</p><p>‖ y u ( ( k + 1 ) T ) ‖ 2 − ‖ y u ( k T ) ‖ 2 ≤ − 2 ρ ∫ k T ( k + 1 ) T | 〈 y u ( τ ) , B u y u ( τ ) 〉 | 2 1 + | 〈 y u ( τ ) , B u y u ( τ ) 〉 | d τ</p><p>using now the estimate (26), we deduce that</p><p>‖ y u ( ( k + 1 ) T ) ‖ 2 − ‖ y u ( k T ) ‖ 2 ≤ − M ( ∫ 0 T | 〈 S u ( τ ) y u , B u S u ( τ ) y u 〉 | d τ ) 2 (32)</p><p>for some M &gt; 0 . Using now the fact that dim ( H u ) &lt; + ∞ , then the assumption (27) is equivalent to</p><p>∫ 0 T | 〈 B u S u ( t ) y u , S u ( t ) y u 〉 | d t ≥ δ ‖ y u ‖ 2 , ∀ y u ∈ H u , ( T , δ &gt; 0 ) (33)</p><p>From (32) and (33) we have</p><p>‖ y u ( ( k + 1 ) T ) ‖ 2 − ‖ y u ( k T ) ‖ 2 ≤ − M δ 2 ‖ y u ( k T ) ‖ 4</p><p>using the fact that the map t → ‖ y u ( t ) ‖ is decreasing on ℝ + , we obtain:</p><p>‖ y u ( ( k + 1 ) T ) ‖ 2 − ‖ y u ( k T ) ‖ 2 ≤ − M δ 2 ‖ y u ( ( k + 1 ) T ) ‖ 4</p><p>which implies that</p><p>‖ y u ( ( k + 1 ) T ) ‖ 2 + C ‖ y u ( ( k + 1 ) T ) ‖ 4 ≤ ‖ y u ( k T ) ‖ 2 , C = M δ 2</p><p>Letting s k = ‖ y u ( k T ) ‖ 2 , the last inequality can be written as</p><p>s k + 1 + C s k + 1 2 ≤ s k , ∀ k ≥ 0</p><p>From Lemma.3.1 we have</p><p>s k ≤ M 2 k + 1</p><p>For k = [ t T ] ( [ t T ] designed the integer part of t T ), then we obtain s k ≤ M 3 t , ( M 3 &gt; 0 ) , which gives</p><p>‖ y u ( t ) ‖ 2 ≤ M 3 t</p><p>Hence</p><p>‖ y u ( t ) ‖ = O ( t − 1 2 ) , as t → + ∞ (34)</p><p>For the component y s ( t ) , we shall show that y s ( t ) is defined for all t ≥ 0 and exponentially converges to 0, as t → + ∞ . The system (1) excited by the constrained feedback control (28) admits a unique mild solution defined for almost all t in a maximal interval [ 0 , t max [ defined by</p><p>y ( t ) = S ( t ) y 0 + ∫ 0 t S ( t − τ ) p ρ , u ( τ ) B y ( τ ) d τ</p><p>Thus</p><p>y s ( t ) = S s ( t ) y 0 s + ∫ 0 t S s ( t − τ ) p ρ , u ( τ ) B s y s ( τ ) d τ , ∀ t ∈ [ 0 , t max [ (35)</p><p>It follows from (13) that</p><p>‖ y s ( t ) ‖ ≤ M 1 e − η t ‖ y 0 s ‖ + ρ M 1 ‖ B ‖ ∫ 0 t e − η ( t − τ ) ‖ y s ( τ ) ‖ d τ</p><p>For almost all t ∈ [ 0 , t max [ . The Gronwall inequality then yields:</p><p>‖ y s ( t ) ‖ ≤ M 1 ‖ y 0 s ‖ e ( ρ M 1 ‖ B ‖ − η ) t , ∀ t ≥ 0 (36)</p><p>Taking ρ &lt; η M 1 ‖ B ‖ , it follows from (36) that y s ( t ) is bounded on [ 0 , t max [ so t max = + ∞ , and therefore</p><p>(36) holds for all t ≥ 0 . Hence, from (34) and (36), the solution of (1) satisfies the estimate (8). This completes the proof of Theorem 3.3.</p></sec><sec id="s3_4"><title>3.4. Remark</title><p>1. Since the function t → ‖ y u ( t ) ‖ decreasing in ℝ + , we have</p><p>∃ t 0 ≥ 0 ; y u ( t 0 ) = 0 ⇔ y u ( t ) = 0 , ∀ t ≥ t 0</p><p>In this case, we have</p><p>p ρ , u ( t ) = 0 , ∀ t ≥ t 0 ⇒ y ( t ) = S s ( t − t 0 ) y s ( t 0 ) , ∀ t ≥ t 0</p><p>Hence, using (13) the system (1) is exponentially stable.</p><p>2. The constrained feedback control (28) depends only on the unstable part y u ( t ) and we have</p><p>| p ρ , u ( t ) | &lt; ρ , ∀ t ≥ 0</p><p>3. The constrained feedback control (28) satisfies</p><p>| p ρ , u ( t ) | = O ( 1 t ) , as t → + ∞</p><p>4. We note that (27) is weaker than (6). The converse is not true as we can see taking an orthonormal basis ( ϕ n ) n ≥ 1 of H , A z = − 〈 z , ϕ 1 〉 ϕ 1 and B z = ∑ n = 1 + ∞ 〈 z , ϕ n 〉 n 2 ϕ n .</p><p>5. In the case dim H u = + ∞ and B is nonlinear and locally Lipschitz, such that B ( 0 ) = 0 , then using the same techniques as in [<xref ref-type="bibr" rid="scirp.55169-ref9">9</xref>] , we can obtain the result of Theorem 3.3, if the estimate (7) is changed to (33).</p></sec></sec><sec id="s4"><title>4. Robustness</title><p>In this section, we study the robustness of the controls (18) and (28), under a class of perturbations of the system (1).</p><sec id="s4_1"><title>4.1. Strong Robustness</title><p>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.</p><p>Let us consider the following perturbed system</p><p>d y ( t ) d t = A y ( t ) + p ( t ) B y ( t ) + ξ ( y ( t ) ) , y ( 0 ) = y 0 (37)</p><p>where the linear bounded operator ξ = ξ u + ξ s is such that the system (37) is decomposed into two following subsystems:</p><p>d y u ( t ) d t = A u y u ( t ) + p ( t ) B u y u ( t ) + ξ u ( y u ( t ) ) , y u ( 0 ) = y 0 u ∈ H u (38)</p><p>d y s ( t ) d t = A s y s ( t ) + p ( t ) B s y s ( t ) + ξ s ( y u ( t ) ) , y s ( 0 ) = y 0 s ∈ H s (39)</p><p>The following main result concerns the strong stability of the system (37).</p></sec><sec id="s4_2"><title>4.2. Proposition</title><p>Let</p><p>1. A generates a linear C 0 -contraction semigroup ( S ( t ) ) t ≥ 0 on H such that (13) holds,</p><p>2. The operator B is compact such that (6) holds,</p><p>3. The linear operator ξ is compact and satisfying ‖ ξ s ‖ &lt; η M 1 and 〈 ξ ( y ) , y 〉 ≤ 0 , ∀ y ∈ H .</p><p>Then the system (37) is strongly stabilizable.</p><p>Proof</p><p>First, let us note that 0 remains an equilibrium state of the perturbed system (37), which can be written in the form:</p><p>d y ( t ) d t = A y ( t ) + g ( y ( t ) ) , y ( 0 ) = y 0 (40)</p><p>where g = f + ξ and f ( y ) = { − ρ 〈 B y , y 〉 1 + | 〈 B y , y 〉 | B y , y ≠ 0 ; 0 , y = 0.</p><p>Since f and ξ are locally Lipschitz, so is g . Also g is dissipative: 〈 ξ ( y ) , y 〉 ≤ 0 , ∀ y ∈ H .</p><p>The assumption 〈 ξ ( y ) , y 〉 ≤ 0 , ∀ y ∈ H , together with (6) guarantees the following implication</p><p>〈 g ( S ( t ) y ) , S ( t ) y 〉 = 0 ⇒ y = 0</p><p>Then the weak stability of the perturbed system (37) follows from Theorem 2.4 of Ball [<xref ref-type="bibr" rid="scirp.55169-ref3">3</xref>] , and since dim H u &lt; + ∞ , we have y u ( t ) → 0 , as t → + ∞ . For the component y s ( t ) of the solution y ( t ) of the system (37), and for all 0 ≤ t 0 ≤ t , we have:</p><p>y s ( t ) = S s ( t − t 0 ) y s ( t 0 ) − ρ ∫ t 0 t S s ( t − τ ) 〈 y ( τ ) , B y ( τ ) 〉 1 + | 〈 y ( τ ) , B y ( τ ) 〉 | B s y s ( τ ) d τ + ∫ t 0 t S s ( t − τ ) ξ s ( y s ( τ ) ) d τ (41)</p><p>It follows from (13) that</p><p>‖ y s ( t ) ‖ ≤ M 1 e − η ( t − t 0 ) ‖ y s ( t 0 ) ‖ + ( ρ M 1 ‖ B ‖ + M 1 ‖ ξ s ‖ ) ∫ t 0 t e − η ( t − τ ) ‖ y s ( τ ) ‖ d τ (42)</p><p>From Gronwall inequality, we obtain:</p><p>‖ y s ( t ) ‖ ≤ M 1 e ( ρ M 1 ‖ B ‖ + M 1 ‖ ξ s ‖ − η ) ( t − t 0 ) ‖ y s ( t 0 ) ‖ , ∀ t ≥ t 0</p><p>Taking ρ &lt; η − M 1 ‖ ξ s ‖ M 1 ‖ B ‖ , we obtain y s ( t ) → 0 , as t → + ∞ . Hence, the solution y ( t ) of the system (37) strongly converges to 0 , as t → + ∞ .</p></sec><sec id="s4_3"><title>4.3. A Polynomial Decay Estimate for the Perturbed System</title><p>Our second main result in this section is stated as follows:</p></sec><sec id="s4_4"><title>4.4. Proposition</title><p>Let</p><p>1. A generate a linear C 0 -semigroup ( S ( t ) ) t ≥ 0 such that ( S u ( t ) ) t ≥ 0 is a semigroup of isometries and (13) holds,</p><p>2. A allows the decomposition (10) of H with dim H u &lt; + ∞ ,</p><p>3. B ∈ L ( H ) satisfies (27),</p><p>4. ‖ ξ s ‖ ≤ η M 1 and 〈 ξ u ( y ) , y 〉 ≤ 0 , for all y ∈ H u .</p><p>Then the constrained feedback control (28) strongly stabilizes the system (37) with the explicit decay estimate (8).</p><p>Proof</p><p>Let us consider the system:</p><p>d y u ( t ) d t = A u y u ( t ) − ρ 〈 y u ( t ) , B u y u ( t ) 〉 1 + | 〈 y u ( t ) , B u y u ( t ) 〉 | B u y u ( t ) + ξ u ( y u ( t ) ) , y u ( 0 ) = y 0 u (43)</p><p>Multiplying the system (43) by y u ( t ) and integrating over Ω and using the fact that ( S u ( t ) ) t ≥ 0 is a semigroup of isometries and the hypothesis: 〈 ξ u ( y u ) , y u 〉 ≤ 0 , ∀ y u ∈ H u , we obtain:</p><p>d ‖ y u ( t ) ‖ 2 d t ≤ − 2 ρ | 〈 y u ( t ) , B u y u ( t ) 〉 | 2 1 + | 〈 y u ( t ) , B u y u ( t ) 〉 | (44)</p><p>which gives ‖ y u ( t ) ‖ ≤ ‖ y u ( 0 ) ‖ , ∀ t ≥ 0 . Then, the system (43) admits a unique global mild solution y ( t ) de-</p><p>fined for almost all t ≥ 0 . By the same argument as in the proof of Theorem 3.3, the solution y ( t ) of the system (37) satisfies:</p><p>‖ y ( t ) ‖ = O ( t − 1 2 ) , as t → + ∞</p><p>which completes the proof of Proposition 4.2.</p></sec></sec><sec id="s5"><title>5. Application and Simulations</title><sec id="s5_1"><title>5.1. An Application</title><p>In this part, we will give an illustrating example of the established results.</p>Example<p>Let us consider the following 1-d bilinear heat equation:</p><p>{ ∂ y ( x , t ) ∂ t = ∂ 2 y ( x , t ) ∂ x 2 + p ( t ) B y ( t ) , x ∈ ( 0 , 1 ) , t &gt; 0 , ∂ y ( 0 , t ) ∂ x = ∂ y ( 1 , t ) ∂ x = 0 , ∀ t &gt; 0 , (45)</p><p>where y ( t ) is the temperature profile at time t . We suppose that the system is controlled via the flow of a liquid p ( t ) in an adequate metallic pipeline. Here we take the state space H = L 2 ( 0 , 1 ) and the operator A is defined by A y = ∂ 2 y ∂ x 2 , with D ( A ) = { y ∈ H 2 ( 0 , 1 ) | ∂ y ( 0 , t ) ∂ x = ∂ y ( 1 , t ) ∂ x = 0 } . The domain of A gives the ho-</p><p>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-</p><p>values λ j = − π 2 ( j − 1 ) 2 , j ∈ ℕ ∗ and eigenfunctions φ 1 ( x ) = 1 and φ j ( x ) = 2 cos ( ( j − 1 ) π x ) for all j ≥ 2 . Then the subspace H u is the one-dimensional space spanned by the eigenfunction φ 1 , and we have S u ( t ) y u = 〈 y u , φ 1 〉 φ 1 so S u ( t ) = I H u (the identity) and hence ( S u ( t ) ) t ≥ 0 , is a semigroup of isometries. The operator of control B , is defined by: B y = ∑ j = 1 + ∞ α j 〈 y , φ j 〉 φ j , α j ≥ 0 , ∀ j ≥ 1 , such that ∑ j = 1 + ∞ α j 2 &lt; ∞ (see [<xref ref-type="bibr" rid="scirp.55169-ref16">16</xref>] ). From the relation: 〈 B u S u ( t ) y u , S u ( t ) y u 〉 = α 1 | 〈 y u , φ 1 〉 | 2 , we can see that (27) holds if α 1 &gt; 0 . To examine the estimate (8), remarking for the scalar functions y j ( t ) = 〈 y ( t ) , φ j 〉 , ∀ j ≥ 1 we have</p><p>d y 1 ( t ) d t = d y u ( t ) d t = − ρ α 1 y u 3 ( t ) 1 + α 1 y u 2 ( t ) , ∀ t ≥ 0 (46)</p><p>which implies that</p><p>d d t y u 2 ( t ) = − 2 ρ α 1 y u 4 ( t ) 1 + α 1 y u 2 ( t ) (47)</p><p>Letting x ( t ) = y u 2 ( t ) , we obtain</p><p>x ′ ( t ) x − 2 ( t ) = − 2 ρ α 1 1 + α 1 x ( t ) , ∀ t &gt; 0</p><p>Integrating now the last equality from 0 to t , we get</p><p>x − 1 ( t ) − x − 1 ( 0 ) = ∫ 0 t 2 ρ α 1 1 + α 1 x ( τ ) d τ</p><p>from (47), we deduce that the nonnegative scalar map t → x ( t ) is decreasing for all t ≥ 0 , and we have</p><p>x − 1 ( t ) ≥ 2 ρ α 1 t 1 + α 1 x ( 0 ) , ∀ t &gt; 0</p><p>which means that</p><p>| y u ( t ) | ≤ 1 + α 1 y 0 u 2 2 ρ α 1 t , ∀ t &gt; 0 , y 0 u ≠ 0</p><p>Then</p><p>| y u ( t ) | = O ( 1 t ) , as t → + ∞</p><p>Furthermore, the control in this case is defined by</p><p>p ρ , u ( t ) = − ρ α 1 y u 2 ( t ) 1 + α 1 y u 2 ( t ) , ∀ t ≥ 0 (48)</p><p>For j ≥ 2 , the functions y j ( t ) are characterized by y j ( 0 ) = 〈 y ( 0 ) , φ j 〉 , ∀ j ≥ 2 and satisfy</p><p>∂ ∂ t y j ( t ) = ( λ j − ρ α 1 α j y 1 2 ( t ) 1 + α 1 y 1 2 ( t ) ) y j ( t ) , t ≥ 0</p><p>which implies that</p><p>| y j ( t ) | ≤ e − π 2 ( j − 1 ) 2 t | y j ( 0 ) | , j ≥ 2</p><p>Then</p><p>‖ y s ( t ) ‖ ≤ e − π 2 t ‖ y s ( 0 ) ‖ , ∀ t ≥ 0</p><p>Hence, the system (45) is strongly stable with the decay rate estimate (8).</p><p>Let us reconsider the above example with the perturbation ξ defined by:</p><p>ξ ( y ) = − β y , β ∈ ( 0 , π 2 2 ) , ∀ y ∈ H</p><p>It is clear that the function ξ 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.</p></sec><sec id="s5_2"><title>5.2. Simulations</title><p>In this part, taking in the system (45), the operator B = I and y 0 ( x ) = x + 0.3 .</p><p>Then we obtain the results shown in Figures 1-5.</p></sec></sec><sec id="s6"><title>6. Conclusion</title><p>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</p><p>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.</p></sec></body><back><ref-list><title>References</title><ref id="scirp.55169-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Mohler, R.R. (1973) Bilinear Control Process. Academic Press, New York.</mixed-citation></ref><ref id="scirp.55169-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">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</mixed-citation></ref><ref id="scirp.55169-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">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</mixed-citation></ref><ref id="scirp.55169-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">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</mixed-citation></ref><ref id="scirp.55169-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">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</mixed-citation></ref><ref id="scirp.55169-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">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</mixed-citation></ref><ref id="scirp.55169-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">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</mixed-citation></ref><ref id="scirp.55169-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">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</mixed-citation></ref><ref id="scirp.55169-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">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</mixed-citation></ref><ref id="scirp.55169-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">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</mixed-citation></ref><ref id="scirp.55169-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Kato, T. (1980) Perturbation Theory for Linear Operators. Springer, New York.</mixed-citation></ref><ref id="scirp.55169-ref12"><label>12</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Ouzahra</surname><given-names> M. </given-names></name>,<etal>et al</etal>. (<year>2011</year>)<article-title>Feedback Stabilization of Parabolic Systems with Bilinear Controls</article-title><source> Electronic Journal of Differential Equations</source><volume> 38</volume>,<fpage> 1</fpage>-<lpage>10</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref><ref id="scirp.55169-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">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</mixed-citation></ref><ref id="scirp.55169-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">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</mixed-citation></ref><ref id="scirp.55169-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">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</mixed-citation></ref><ref id="scirp.55169-ref16"><label>16</label><mixed-citation publication-type="other" xlink:type="simple">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</mixed-citation></ref></ref-list></back></article>