趙轉(zhuǎn)萍

(山西大學商務學院 基礎教學部， 山西 太原 030031)

0 引 言

梁方程是一種雙曲方程， 是偏微分方程的重要研究內(nèi)容之一. 在控制領域， 各式各樣的動態(tài)反饋產(chǎn)生了大量的含有梁方程的耦合系統(tǒng)[1-2]. 因此， 耦合的梁方程的研究逐漸得到了人們的重視[3]. 與上述文獻不同， 本文考慮Euler-Bernoulli梁和抽象系統(tǒng)組成的耦合系統(tǒng).

首先在有界區(qū)間上考慮如下Euler-Bernoulli梁的控制問題：梁的一端是簡支的， 另一端有一個動態(tài)反饋

(1)

式中：G是 Hilbert 空間上的自伴的嚴格正算子；z*∈Z是任意給定的非零元素； (w0,w1,z0)是初始狀態(tài)；w表示振動的梁的位移；z表示抽象的動態(tài)控制變量. 容易看出， 不同的空間Z和不同的算子G將產(chǎn)生不同的耦合系統(tǒng). 當G是n階正定矩陣時， 系統(tǒng)(1)變成聲學結(jié)構(gòu)模型[4-6]， 當G是微分算子G=-?xx時， 系統(tǒng)(1)成為一個梁-熱耦合系統(tǒng)[7]. 系統(tǒng)(1)主要考慮線性動態(tài)反饋， 因此與文獻[8]考慮的問題不同.

在狀態(tài)空間(2)中考慮系統(tǒng)(1)，

(2)

其中，H=L2(0,1)， 且

f′(0)=0}.

(3)

系統(tǒng)(1)的能量函數(shù)為

(4)

對能量函數(shù)求導可得

(5)

由于等式(5)右端沒有w的任何信息， 因此用傳統(tǒng)的李雅普諾夫函數(shù)分析方法來證明系統(tǒng)的穩(wěn)定性將會變得非常困難. 受文獻[9]的啟發(fā)， 本文將從控制的角度來考慮系統(tǒng)解的適定性和漸近穩(wěn)定性， 這為偏微分方程的穩(wěn)定性分析提供了一個新思路.

1 適定性和漸近穩(wěn)定性

近年來， 有很多文獻在研究耦合系統(tǒng)解的適定性和漸近穩(wěn)定性， 例如[10]和 [11]， 然而這些研究大多是基于偏微分方程的方法. 與此不同， 本文將從控制的角度來考慮系統(tǒng)解的適定性和漸近穩(wěn)定性. 首先， 將耦合系統(tǒng)(1)分解成兩部分：控制系統(tǒng)和動態(tài)反饋系統(tǒng). 事實上， 耦合系統(tǒng)(1) 恰好是控制系統(tǒng)

(6)

動態(tài)反饋下的閉環(huán)系統(tǒng)

u(t)=〈z(t),z*〉Z,

(7)

其中，z滿足

(8)

定義算子A∶D(A)?L2(0,1)→L2(0,1)為

(9)

其中，A是H=L2(0,1)上的嚴格正定的自伴算子[12]， 并且有如下Gelfand三嵌入

D(A1/2)→H→D(A1/2)′,

(10)

〈A1/2x,A1/2y〉H,?x,y∈D(A1/2)，

(11)

所以, 系統(tǒng)(6)可以寫成抽象系統(tǒng)

(12)

Bu=uδ(x-1), ?u∈U=C,

(13)

這里δ(·)為Dirac分布. 如果定義B的共軛算子B*∈L(D(A1/2),U)為

〈B*z,u〉U=〈z,Bu〉D(A1/2),D(A1/2)′,

?z∈D(A1/2),u∈U，

(14)

則直接計算可得

B*f=f(1), ?f∈D(A1/2).

(15)

令F∈L(Z,U)滿足

Fz=-〈z,z*〉Z, ?z∈Z，

(16)

則對任意的u∈U,F*∈L(U,Z)滿足F*u=-uz*. 由于z*≠0， 因此

KerF*={0}.

(17)

借助于上面定義的算子， 方程(7)和(8)可以寫成

u(t)=-Fz(t)

(18)

和

(19)

所以， 將系統(tǒng)(1)寫成了如下抽象系統(tǒng)

(20)

其中，C=BF∈L(D(G1/2),D(A1/2)′)，C*=F*B*∈L(D(A1/2),D(G1/2)′) 滿足

anomalies over the western Pacific Ocean warm pool and its possible mechanism

〈C*g,z〉Z=〈B*g,Fz〉U=〈g,Cz〉D(A1/2),D(A1/2)′,

?g∈D(A1/2),z∈Z.

(21)

在狀態(tài)空間H中考慮抽象系統(tǒng)(20). 如果定義A∶D(A)(?H)→H為

(22)

則系統(tǒng)(20)進一步可寫成H中的發(fā)展方程

(23)

定理1 設G是Hilbert空間Z上自伴的、 嚴格正的算子， 0≠z*∈Z. 如果G-1在Z是緊的， 則對任意的(w(·,t),wt(·,0),z(0))∈H， 系統(tǒng)(1)在H上存在唯一的C0-半群解. 此外, 該解是漸近穩(wěn)定的

‖w(·,t),wt(·,t),z(t)‖H→0, 當t→∞.

(24)

證明簡單計算可知， 對任意的(f,g,z)∈D(A)， 有

Re〈A(f,g,z),(f,g,z)〉H=

〈Cz,g〉D(A1/2)′,D(A1/2)-〈Gz,z〉Z-〈C*g,z〉Z]=

(25)

這說明算子A在H中是耗散的. 對任意的(φ,φ,ψ)∈H， 解方程A(f,g,z)=(φ,φ,ψ)∈H可得

(26)

φ,-G-1(F*B*φ+ψ)).

(27)

因此,A-1存在且在H上有界. 由 Lumer-Phillips 定理[13]，A生成H上的C0-壓縮半群.

為了證明系統(tǒng)(1)的漸近穩(wěn)定性， 首先考慮

Ker(B*)∩Ker(λ-A)={0}, ?λ>0,

(28)

式(28)等價于：對任意的λ>0， 特征方程(29)只有零解

(29)

直接計算可得： 式(29)只有零解是顯然的， 從而式(28)成立.

另一方面， 由于G-1在Z中緊， 所以， 式(27)意味著A-1在H中緊. 因此， 由文獻 [14]可知， 只需證明算子A在虛軸上沒有譜點即可. 事實上， 如果設A(f,g,z)=iω(f,g,z), (f,g,z)∈D(A)且0≠ω∈R. 則

(30)

方程(30)等價于

(31)

方程(31)的第一個方程兩邊關于f取內(nèi)積， 同時第二個方程兩邊關于iωz取內(nèi)積可得

(32)

方程(32)的第一個方程乘以-ω2， 再加第二個方程可得

(33)

2 應用舉例

由于系統(tǒng)(1)由算子G和Hilbert空間Z決定， 不同的G和Z意味著不同的耦合系統(tǒng). 下面給出具體的例子來說明這一問題.

例1 如果令G=1,Z=R且z*=1， 則抽象系統(tǒng)(1)變?yōu)?/p>

(34)

簡單計算可知， 上面定義的G和Z滿足定理1的所有條件， 因此系統(tǒng)(34)是漸近穩(wěn)定的. 特別地， 系統(tǒng)(34)對應的半群解為eA1t(w(·,0),wt(·,0),z(0))， 其中系統(tǒng)算子A1定義為

(35)

系統(tǒng)(34)的邊界條件稱為聲學邊界條件. 類似問題可查閱文獻[15].

例2 令Z=L2(0,1)，x∈Z， 定義算子G∶D(G)?Z→Z為Gf=-f″， 其中f∈D(G)={f∈H2(0,1)|f(0)=f′(1)=0}， 則系統(tǒng)(1)變?yōu)?/p>

(36)

與例1類似， 簡單計算可知， 上面定義的G是Hilbert空間Z上自伴的、 嚴格正的算子， 滿足定理1的所有條件， 所以系統(tǒng)(36)是漸進穩(wěn)定的. 此外， 系統(tǒng)(36)對應的半群解為eA2t(w(·,0),wt(·,0),z(0))， 其中系統(tǒng)算子A2定義為

(37)

參考文獻：

[1] Wei Guo, Chen Yunlan, Feng Hongyinping. Output feedback stabilization for a Kirchhoff-type nonlinear beam with general corrupted boundary observation[J]. Int. J. Robust Nonlinear Control, 2017， 27(16)： 3280-3295.

[2] Zhou Huachen， Feng Hongyinping. Disturbance estimator based output feedback exponential stabilization for Euler-Bernoulli beam equation with boundary control[J]. Automatica， 2018， 91(5)： 79-88.

[3] Mercier D， Nicaise S， Sammoury M et al. Optimal energy decay rate for rayleigh beam equation with only one dynamic boundary control[J]. Bol. Soc. Paran. Mat， 2017， 35(3)： 131-171.

[4] Abbas Z， Nicaise S. Polynomial decay rate for a wave equation with general acoustic boundary feedback laws[J]. SeMA Journal， 2013， 61(1)： 19-47.

[5] Abbas Z， Nicaise S. The multidimensional wave equation with generalized acoustic boundary conditions II： strong stability[J]. SIAM J. Control Optim， 2015， 53(4)： 2582-2607.

[6] Abbas Z， Nicaise S. The multidimensional wave equation with generalized acoustic boundary conditions I： strong stability[J]. SIAM J. Control Optim， 2015， 53(4)： 2558-2581.

[7] Zhang Qiong， Wang Junmin， Guo Baozhu. Stabilization of the Euler-Bernoulli equation via boundary connection with heat equation[J]. Math. Control Signals Syst， 2014， 26(1)： 77-118.

[8] Le Gorrec Y， Zwart H， Ramirez H. Asymptotic stability of an Euler-Bernoulli beam coupled to non-linear spring-damper systems[J]. IFAC-papersonline， 2017， 50(1)： 5580-5585.

[9] Feng Hongyinping， Guo Baozhu. On stability equivalence between dynamic output feedback and static output feedback for a class of second order in_nite-dimensional systems[J]. SIAM J. on Control and Optimization， 2015， 53(4)： 1934-1955.

[10] Feng Hongyinping， Li Shengjia. The stability for a one-dimensional wave equation with nonlinear uncertainty on the boundary[J]. Nonlinear Anal， 2013， 89(6)： 202-207.

[11] Vicente A， Frota C. On a mixed problem with a nonlinear acoustic boundary condition for a non-locally reacting boundaries[J]. Math. Anal. Appl， 2013， 407(2)： 328-338.

[12] Tucsnak M， Weiss G. Observation and control for operator semigroups[M]. Basel： Birkhauser， 2009.

[13] Pazy A. Semigroups of linear operators and applications to partial differential equations[M]. New York： Springer-Verlag， 1983.

[14] Luo Zhenghua， Guo Baozhu， Morgul O. Stability and stabilization of infinite dimensional systems with applications[M]. London： Springer-Verlag， 1999.

[15] Wehbe A. Rational energy decay rate for a wave equation with dynamical control[J]. Appl. Math. Lett， 2003， 16(3)： 357-364.