In this manuscript, we study a thin and narrow plate equation that models the deck of a suspension bridge that is subject to a Balakrishnan-Taylor damping and a strong damping. First, by using the Faedo Galerkin method, we prove the existence of both global weak and regular solutions. Second, we prove the exponential stability of the energy for regular solutions by combining the multiplier method and a well-known result of Komornik.
Citation: Zayd Hajjej. On the exponential decay of a Balakrishnan-Taylor plate with strong damping[J]. AIMS Mathematics, 2024, 9(6): 14026-14042. doi: 10.3934/math.2024682
[1] | Patarawadee Prasertsang, Thongchai Botmart . Improvement of finite-time stability for delayed neural networks via a new Lyapunov-Krasovskii functional. AIMS Mathematics, 2021, 6(1): 998-1023. doi: 10.3934/math.2021060 |
[2] | Sheng Dong, Qingwen Wang, Lei Hou . Determinantal inequalities for block Hadamard product and Khatri-Rao product of positive definite matrices. AIMS Mathematics, 2022, 7(6): 9648-9655. doi: 10.3934/math.2022536 |
[3] | Ali Algefary . Diagonal solutions for a class of linear matrix inequality. AIMS Mathematics, 2024, 9(10): 26435-26445. doi: 10.3934/math.20241286 |
[4] | Wentao Le, Yucai Ding, Wenqing Wu, Hui Liu . New stability criteria for semi-Markov jump linear systems with time-varying delays. AIMS Mathematics, 2021, 6(5): 4447-4462. doi: 10.3934/math.2021263 |
[5] | Yanmei Xue, Jinke Han, Ziqiang Tu, Xiangyong Chen . Stability analysis and design of cooperative control for linear delta operator system. AIMS Mathematics, 2023, 8(6): 12671-12693. doi: 10.3934/math.2023637 |
[6] | Rupak Datta, Ramasamy Saravanakumar, Rajeeb Dey, Baby Bhattacharya . Further results on stability analysis of Takagi–Sugeno fuzzy time-delay systems via improved Lyapunov–Krasovskii functional. AIMS Mathematics, 2022, 7(9): 16464-16481. doi: 10.3934/math.2022901 |
[7] | Huahai Qiu, Li Wan, Zhigang Zhou, Qunjiao Zhang, Qinghua Zhou . Global exponential periodicity of nonlinear neural networks with multiple time-varying delays. AIMS Mathematics, 2023, 8(5): 12472-12485. doi: 10.3934/math.2023626 |
[8] | Erfeng Xu, Wenxing Xiao, Yonggang Chen . Local stabilization for a hyperchaotic finance system via time-delayed feedback based on discrete-time observations. AIMS Mathematics, 2023, 8(9): 20510-20529. doi: 10.3934/math.20231045 |
[9] | Zahra Eidinejad, Reza Saadati, Donal O'Regan, Fehaid Salem Alshammari . Measure of quality and certainty approximation of functional inequalities. AIMS Mathematics, 2024, 9(1): 2022-2031. doi: 10.3934/math.2024100 |
[10] | Ahmad Y. Al-Dweik, Ryad Ghanam, Gerard Thompson, M. T. Mustafa . Algorithms for simultaneous block triangularization and block diagonalization of sets of matrices. AIMS Mathematics, 2023, 8(8): 19757-19772. doi: 10.3934/math.20231007 |
In this manuscript, we study a thin and narrow plate equation that models the deck of a suspension bridge that is subject to a Balakrishnan-Taylor damping and a strong damping. First, by using the Faedo Galerkin method, we prove the existence of both global weak and regular solutions. Second, we prove the exponential stability of the energy for regular solutions by combining the multiplier method and a well-known result of Komornik.
A matrix A∈Rn×n is said to be stable if all its eigenvalues lie in the open left half plane, i.e., all the eigenvalues of A have negative real parts. It is well-known that a matrix A is stable if and only if there is a positive definite matrix P such that ATP+PA is negative definite. This implies that V(x)=xTPx serves as a quadratic Lyapunov function for the asymptotically stable linear system
˙x=Ax. |
In this paper, we consider only real square matrices. Let A be a real n×n matrix. A≻0 (A⪰0, resp.) means A is a symmetric positive definite (semidefinite, resp.) matrix.
For the sake of convenience, we will adopt the concept of positive stability. A matrix A∈Rn×n is defined as positive stable if all its eigenvalues possess positive real parts. Clearly, if A is positive stable, then −A is stable. Therefore, results in positive stability can be translated into terms of stability.
It is well-established that a matrix A∈Rn×n is positive stable if and only if there exists P≻0 in Rn×nsuch that
ATP+PA≻0. | (1.1) |
In this case, P is known as a Lyapunov solution for A or to the Lyapunov inequality (1.1). Several numerical methods have been developed to address the problem of finding such matrices P [1,2,3].
A particular case emerges from (1.1) when a positive diagonal matrix D satisfies the Lyapunov inequality. If so, D is called a Lyapunov diagonal solution for A. Furthermore, we say that A is a Lyapunov diagonally stable matrix. The problem of Lyapunov diagonal stability is well investigated in the literature ([4,5,6,7,8,9] and the references therein). The importance of this problem is due to its applications in, most significantly, population dynamics [10], communication networks [11], and systems theory [12].
Another case of (1.1), known as Lyapunov α-scalar stability, appeard in [13]. For a partition α={α1,…,αs} of the set {1,…,n}, the diagonal solution D has an α-scalar structure, i.e. D[αi]=ciI, ci∈R, i=1,…,s, where D[αi] is the principal submatrix of D on row and column indices αi. A set α={α1,…,αs}, 1≤s≤n, is said to be a partition of {1,…,n} if for all i,j∈{1,…,s}, αi≠∅, αi∩αj=∅, and αi∪⋯∪αs={1,…,n}. We assume, without loss of generality, that these αi's are taken to have contiguous indices because our results are applicable with simultaneous row and column permutations.
For brevity, if A∈Rk×k is a Lyapunov diagonally stable matrix, we will write A∈LDSk. Similarly, we write A∈LDSαk if A is Lyapunov α-scalar stable.
A recent generalization of Lyapunov diagonal stability to a family of real matrices of the same size, A={A(i)}ri=1, has drawn significant interest [14,15,16,17,18]. This extension studies the existence of a diagonal matrix D≻0 satisfying
(A(i))TD+DA(i)≻0, | (1.2) |
i=1,…,r. If such matrix D exists, it is known as a common Lyapunov diagonal solution for A or to (1.2). Consequently, we say A has common Lyapunov diagonal stability. From this definition, it is clear that common Lyapunov diaognal stability can be interpreted as simultaneous Laypunov diagonal stability for the matrices in A. The existence of a common Lyapunov diagonal solution D for A implies that V(x)=xTDx acts as a common Lyapunov diagonal function for the collection of asymptotically stable linear systems
˙x=A(i)x,i=1,…,r. |
An immediate observation here is that when A=A, i.e., A is a singleton, A∈CLDS is equivalent to A∈LDS, and A∈CLDSα is equivalent to A∈LDSα. Additionally, it is worth mentioning that the cardinality of A is not relevant. For convenience, we shall fix it to be r throughout the rest of this note.
Applications of common Lyapunov diagonal stability have been found in the fields of large-scale dynamics [19,20,21,22], as well as in the study of interconnected time-varying and switched systems [18]. Beyond these practical applications, common Lyapunov diagonal stability is also a significant research topic in itself, as evidenced by works such as [14,16,17,23].
Let A∈Rn×n be a nonsingular matrix. In [24], Redheffer proved that A∈LDSn if and only if the (n−1)×(n−1) leading principal submatrices of A and A−1 have a common Lyapunov diagonal solution. This result has been restated in [16,23] using the notion of Schur complements. The new statement is free of the nonsingularity condition. Specifically, it was shown that a matrix A∈LDSn if and only if ann>0 and the (n−1)×(n−1) leading principal submatrix of A and its Schur complement have a common Lyapunov diagonal solution.
For any vectors u,v∈Rn, when we write u≫v, it means ui>vi for all i∈{1,…,n}. For a matrix A∈Rn×n, the vector u∈Rn with ui=aii, for i=1,…,n, is denoted by diag(A). We denote the identity matrix I∈Rk×k by Ik and the matrix of all ones in Rk×k by Jk.
Let A=[aij],B=[bij]∈Rn×n and C∈Rm×m. The Hadamard product of A and B is denoted by A∘B=[aijbij]∈Rn×n. The Kronecker product of A and C is denoted by A⊗C=[aijC]∈Rnm×nm.
Let n,m∈N and for i=1,…,n, let mi∈N such that m=m1+⋯+mn, where n≤m. Then, a matrix S∈Rm×m is an n by n block matrix if it is partitioned into blocks that conform with mi, i=1,…,n. Moreover, we denote each mj by mk block of S as Sjk. Similarly, a vector u∈Rm is called an n-block vector if it is partitioned into n subvectors, i.e., uT=[uTm1…uTmn], where umi∈Rmi, i=1,…,n. Throughout this note, it is assumed that n,m and all mi, i=1,…,n, are natural number with n≤m and m=m1+⋯+mn.
The Khatri-Rao product of a matrix A=[aij]∈Rn×n and an n by n block matrix B=[Bij]∈Rm×m is defined as A⋆B=[aijBij]∈Rm×m. Similarly, let v=[vi]∈Rn and u=[umi]∈Rm be n-block vector, then v⋆u=[viumi]∈Rm.
Suppose that ∅≠α⊆{1,…,k}, |α| is the cardinality of α and αc={1,…,n}∖α. We denote the principal submatrix of A obtained by selecting rows and columns indexed by α as A[α]. Similarly, for a vector u∈Rk, u[α] represents the subvector of u containing only the elements indexed by α.
Lemma 1.1. ([25,Corollary 4.2.13]) Let A∈Rn×n and B∈Rm×m. If A and B are both positive semidefinite matrices then A⊗B∈Rnm×nm is also a positive semidefinite matrix.
Lemma 1.2. ([26,Theorem 3.1]) Let A∈Rn×n be a positive definite matrix. If B∈Rm×m is a positive semidefinite n by n block matrix with Bii≻0, i=1,…,m, then A⋆B≻0.
For the remainder of this note, Q=[qij]∈Rm×m denotes the nonzero n by n block matrix defined as
(Qij)kl={1if k=l0if k≠l. | (1.3) |
Additionally, for an n by n block matrix B∈Rm×m, we define the matrix T(B)∈Rn×n such that (T(B))ij=(Bij)11.
Now, let us recall the definition of a P-matrix. A matrix whose principal minors are all positive is known as a P-matrix. A well-known characterization for P-matrices in the context of real matrices is given next.
Lemma 1.3. ([27,Theorem 3.3]) A matrix A∈Rn×n is a P-matrix if and only if ui(Au)i>0 for all nonzero u∈Rn.
Motivated by Lemma 1.3, a generalization of the concept of P-matrices to Pα-matrices has been developed in [13].
Definition 1.1. Let A∈Rn×n and α={α1,…,αs} be a partition of {1,…,n}. Then, A is a Pα-matrix if there is some k∈{1,…,s} such that u[αk]T(Au)[αk]>0 for all nonzero u∈Rn.
For A∈Rk×k, A∈Pk indicates that A is a P-matrix, while A∈Pαk means that A is Pα-matrix.
Using the characterization in Lemma 1.3 and Definition 1.1, the P-matrix and Pα-matrix properties were extended in [17,28], respectively, to consider a family of real matrices.
Definition 1.2. Let A={A(i)}ri=1 be a family of real n×n matrices. Then, A is called a P-set and write A∈Pn if for any family of vectors {u(i)}ri=1 in Rn, not all being zero, there is some k∈{1,…,n} such that
r∑i=1u(i)k(A(i)u(i))k>0. |
Definition 1.3. Let A={A(i)}ri=1 be a family of real n×n matrices and α={α1,…,αs} be a partition of {1,…,n}. Then, A is called a Pα-set and write A∈Pαn if for any family of vectors {u(i)}ri=1 in Rn, not all being zero, there is some k∈{1,…,s} such that
r∑i=1u(i)[αk]T(A(i)u(i))[αk]>0. |
Theorem 1.4. ([14,Theorem 2]) Let A={A(i)}ri=1 be a family of matrices such that A(i)∈Rn×n, i=1,…,r. Then, A∈CLDSn if and only if the matrix
r∑i=1A(i)H(i) |
has a positive diagonal entry for any H(i)⪰0 in Rn×n, i=1,…,r, not all being zero.
Theorem 1.5. ([17,Theorem 2.5]) Let A={A(i)}ri=1 be a family of matrices such that A(i)∈Rn×n, i=1,…,r. Then, the following are equivalent:
(i) A∈CLDSn.
(ii) {A(i)∘S(i)}ri=1∈CLDSn for all S(i)⪰0, with diag(S(i))≫0 for i=1,…,r.
(iii) {A(i)∘S(i)}ri=1∈CLDSn for all S(i)⪰0, with diag(S(i))=e for i=1,…,r.
(iv) {A(i)∘S(i)}ri=1∈Pn for all S(i)⪰0, with diag(S(i))≫0 for i=1,…,r.
(v) {A(i)∘S(i)}ri=1∈Pn for all S(i)⪰0, with diag(S(i))=e for i=1,…,r.
The above two theorems provide characterizations for common Lyapunov diagonal stability. Theorem 1.4 extends Theorem 1 from [4], while Theorem 1.5 is inspired by the work of Kraaijevanger [8]. The primary objective of our work is to offer additional characterizations that enhance and unify the existing results in the literature.
We begin this section with a lemma that gives a necessary condition for the common Lyapunov diagonal stability.
Lemma 2.1. ([17,Theorem 2.3]) Let A={A(i)}ri=1 be a family of real n×n matrices. If A∈CLDSn, then A∈Pn.
Next, we demonstrate that if a family of matrices A of the same size forms a P-set, then any family of principal submatrices of A obtained by deleting the same rows and columns also forms a P-set.
Lemma 2.2. Let A={A(i)}ri=1 be a family of real n×n matrices and ∅≠α⊆{1,…,n}. If A∈Pn, then B={A[α]}ri=1∈P|α|.
Proof. For i=1,…,r, let v(i)∈R|α|, not all being zero. Then, for each i, construct u(i)∈Rn to be such that u(i)[α]=v(i) and u(i)[αc]=0. Clearly, not all these u(i)'s are zero vectors since not all v(i)'s are zero. Hence, since A∈Pn, there is some k∈{1,…,n} such that
r∑i=1u(i)k(A(i)u(i))k>0. |
Observe that for each i, u(i)k=v(i)l and (A(i)u(i))k=(A[α](i)v(i))l for some l∈α. Otherwise, the above summation equals zero. From this observation, we obtain that
r∑i=1v(i)l(A(i)[α]v(i))l>0. |
Therefore, by Definition 1.2, B∈P|α|.
We are now ready to present our main theorem.
Theorem 2.3. Let A={A(i)}ri=1 be a family of matrices such that A(i)∈Rn×n for i=1,…,r. Then, the following are equivalent:
(i) A∈CLDSn.
(ii) {A(i)⋆S(i)}ri=1∈CLDSm for all n by n block matrices S(i)⪰0 in Rm×m with S(i)jj≻0 for all i=1,…,r and j=1,…,n.
(iii) {A(i)⋆S(i)}ri=1∈CLDSm for all n by n block matrices S(i)⪰0 in Rm×m with S(i)jj=Imj for all i=1,…,r and j=1,…,n.
(iv) {A(i)⋆S(i)}ri=1∈Pm for all n by n block matrices S(i)⪰0 in Rm×m with S(i)jj≻0 for all i=1,…,r and j=1,…,n.
(v) {A(i)⋆S(i)}ri=1∈Pm for all n by n block matrices S(i)⪰0 in Rm×m with S(i)jj=Imj for all i=1,…,r and j=1,…,n.
Proof. It is trivial to see that (ii) implies (iii) and (iv) implies (v). Moreover, from Lemma 2.1, it is clear that (ii) implies (iv) and (iii) implies (v). Hence, to finish the proof, we show that (i) implies (ii) and (v) implies (i).
(i)⇒(ii): Suppose that D≻0 in Rn×n is a common Lyapunov diagonal solution for A. Then, for i=1,…,r, we have (A(i))TD+DA(i)≻0. Let {S(i)}ri=1 be any family of positive semidefinite n by n block matrices in Rm×m with S(i)jj≻0 for j=1,…,n and i=1,…,r. Hence, according to Lemma 1.2, we have
((A(i))TD+DA(i))⋆S(i)≻0, | (2.1) |
for i=1,…,r. Since we have
((A(i))TD+DA(i))⋆S(i)=((A(i))TD)⋆S(i)+(DA(i))⋆S(i), |
it follows from (2.1) that
((A(i))TD)⋆S(i)+(DA(i))⋆S(i)≻0. | (2.2) |
Now, observe that
(DA(i))⋆S(i)=(D⋆Im)(A(i)⋆S(i)), |
and
((A(i))TD)⋆S(i)=(A(i)⋆S(i))T(D⋆Im) |
for each i, where Im∈Rm×m is the identity matrix partitioned into n by n blocks. Using these observations, it follows from (2.2) that
(A(i)⋆S(i))T(D⋆Im)+(D⋆Im)(A(i)⋆S(i))≻0, |
for i=1,…,r. Clearly, the diagonal matrix D⋆Im∈Rm×m is positive definite. Hence, (ii) follows.
(v)⇒(i): For i=1,…,r, let X(i)=[x(i)kl]⪰0 in Rn×n, not all being zero. Now, set D(i) to be the diagonal matrices whose diagonal elements d(i)kk=√x(i)kk for all i=1,…,r and k=1,…,n. Thus, for each i, we can write X(i)=D(i)S(i)D(i) for some S(i)=[s(i)kl]⪰0 in Rn×n with skk=1, k=1,…,n. Next, let us fix p=max{m1,…,mn}. Then, by Lemma 1.1, S(i)⊗Ip⪰0 in Rnp×np, i=1,…r. Observe that for each i, S(i)⋆Q∈Rm×m is a principal submatrix of S(i)⊗Ip, where Q is a matrix defined as in (1.3). Therefore, we conclude that S(i)⋆Q⪰0, with (S(i)⋆Q)jj=Imj, i=1,…,r, j=1,…,n. By (v), {A(i)⋆(S(i))⋆Q)}ri=1∈Pm. So, we obtain from Lemma 2.2 that {T(A(i)⋆(S(i)⋆Q))}ri=1∈Pn. Now, let u(i)∈Rn×n, i=1,…,r, be such that u(i)k=d(i)kk, k=1,…,n. It is clear that not all u(i) are zero vectors. Thus, from the definition of P-sets, we must have
r∑i=1u(i)q[(T(A(i)⋆(S(i)⋆Q)))u(i)]q>0 |
for some q∈{1,…,n}. Hence, it follows that
r∑i=1u(i)q[(T(A(i)⋆S(i)⋆Q))u(i)]q=r∑i=1d(i)qqn∑k=1(T(A(i)⋆S(i)⋆Q))qkd(i)kk=r∑i=1d(i)qqn∑k=1(T((A(i)∘S(i))⋆Q))qkd(i)kk=r∑i=1d(i)qqn∑k=1(a(i)qks(i)qkQqk)11d(i)kk=r∑i=1d(i)qqn∑k=1a(i)qks(i)qkd(i)kk=r∑i=1n∑k=1a(i)qkd(i)qqs(i)qkd(i)kk=r∑i=1n∑k=1a(i)qkx(i)qk=r∑i=1n∑k=1a(i)qkx(i)kq=(r∑i=1A(i)X(i))qq>0. |
From this last inequality and by Theorem 1.4, (i) holds.
The proof is complete now.
To demonstrate the validity of Theorem 2.3, consider the following example.
Example 2.1. Let n=2, m=3, m1=2, and m2=1. Then, consider the family A={A(i)}2i=1, where
A(1)=[2−103]andA(2)=[1−104]. |
According to Theorem 2.3, to show that A∈CLDS2 it suffices to show that {A(i)⋆S(i)}2i=1∈P3 for any 2 by 2 block matrices S(i)⪰0 in R3×3 with S(i)jj=Imj for all i=1,2 and j=1,2. Now, consider the matrices
S(1)=[10s(1)1301s(1)23s(1)13s(1)231]andS(2)=[10s(2)1301s(2)23s(2)13s(2)231]. |
Hence, we have
A(1)⋆S(1)=[20−s(1)1302−s(1)23003]andA(2)⋆S(2)=[10−s(2)1301−s(2)23004]. |
Next, for i=1,2, let u(i) be any vectors in R3. Thus, a simple calculation shows that
(A(1)⋆S(1))u(1)=[2u(1)1−s(1)13u(1)32u(1)2−s(1)23u(1)33u(1)3]and(A(2)⋆S(2))u(2)=[u(2)1−s(2)13u(2)3u(2)2−s(2)23u(2)34u(2)3]. |
If at least one of u(1)3 and u(2)3 is nonzero, then ∑2i=1u(i)3(A(i)u(i))3>0. Otherwise, we must have
(A(1)⋆S(1))u(1)=[2u(1)12u(1)20]and(A(2)⋆S(2))u(2)=[u(2)1u(2)20]. |
Since u(1) and u(2) are not both zero vectors, then we must have k∈{1,2} such that ∑2i=1u(i)k(A(i)u(i))k>0. That means {A(i)⋆S(i)}2i=1∈P3. Therefore, from Theorem 2.3, this implies that A∈CLDS2. In fact, we found that
D=[21] |
is a common Lyapunov diagonal solution for A.
We emphasize here that Theorem 2.6 is equivalent to Theorem 1.5 when m=n. Before we proceed with the presentation of further results, we cite the following two lemmas from [28].
Lemma 2.4. ([28,Lemma 3.2]) Let A={A(i)}ri=1 be a family of real n×n matrices and α be any partition of {1,…,n}. If A∈CLDSαn, then A∈Pαn.
Lemma 2.5. ([28,Proposition 4.1]) Let A={A(i)}ri=1 be a family of real n×n matrices and α be any partition of {1,…,n}. If A∈Pαn, then A∈Pn.
For the remainder of this paper, let α={α1,…,αn} is a partition of {1,…,m} such that α1={1,…,m1},α2={m1+1,…,m1+m2},…,αn={m−mn+1,…,m}. With this notation established, we now provide another characterization of common Lyapunov diagonal stability.
Theorem 2.6. Let A={A(i)}ri=1 be a family of matrices such that A(i)∈Rn×n for i=1,…,r. Then, the following are equivalent:
(i) A∈CLDSn.
(ii) {A(i)⋆S(i)}ri=1∈CLDSαm for all n by n block matrices S(i)⪰0 in Rm×m with S(i)jj≻0 for all i=1,…,r and j=1,…,n.
(iii) {A(i)⋆S(i)}ri=1∈CLDSαm for all n by n block matrices S(i)⪰0 in Rm×m with S(i)jj=Imj for all i=1,…,r and j=1,…,n.
(iv) {A(i)⋆S(i)}ri=1∈Pαm for all n by n block matrices S(i)⪰0 in Rm×m with S(i)jj≻0 for all i=1,…,r and j=1,…,n.
(v) {A(i)⋆S(i)}ri=1∈Pαm for all n by n block matrices S(i)⪰0 in Rm×m with S(i)jj=Imj for all i=1,…,r and j=1,…,n.
Proof. It is clear that (ii)⇒(iii) and (iv)⇒(v). In addition, according to Lemma 2.4, (ii)⇒(iv) and (iii)⇒(v).
(i)⇒(ii): Let {S(i)}ri=1 be a family of positive semidefinite matrices given as in (ii) and D be a common Lyapunov diagonal solution for A. Then, D⋆Im is a positive α-scalar matrix, where Im∈Rm×m is the identity matrix partitioned into n by n blocks. Thus, as we have seen in the proof of Theorem 2.3, we have
(A(i)⋆S(i))T(D⋆Im)+(D⋆Im)(A(i)⋆S(i))≻0, |
for i=1,…,r.
(v)⇒(i): Using Lemma 2.5, we can see that (v) here implies (v) in Theorem 2.3. Therefore, (i) holds.
Now, we extend Theorem 2.6 by considering different partitions of {1,…,m}. Before presenting our next result, let us set the stage first.
Let us define a bijective function τ:αi→βi that maps each element j∈αi to some βi for i∈{1,…,n}. Hence, τ is a permutation of {1,…,m}, and β is a partition of {1,…,m}. Clearly, for every i, the cardinality of αi is the same as the cardinality of βi. For the remainder of this section, β denotes such partitions. In addition, construct the permutation matrix P such that Pjτ(j)=1 for all j=1,…,m and zero everywhere else. For any permutation matrix P, we write CP=PCPT, where C∈Rm×m. Then, the following observation can be easily verified.
Observation 2.1. Let P be a permutation matrix associated with some partition β. Then, we have
(1) S⪰0 (S≻0, resp.) if and only if SP⪰0 (SP≻0, resp.), S∈Rm×m.
(2) D is β-scalar matrix if and only if DP is α-scalar matrix, D∈Rm×m.
(3) {A(i)}ri=1∈Pβm if and only if {A(i)P}ri=1∈Pαm, where A(i)∈Rm×m, i=1,…,r.
Lemma 2.7. Let A={A(i)}ri=1 be a family of matrices such that A(i)∈Rn×n for i=1,…,r and P be a permutation matrix associated with some partition β. Then, {A(i)}ri=1∈CLDSβn if and only if {A(i)P}ri=1∈CLDSαn.
Proof. The conclusion follows directly from observation 1 and noting that
((A(i))TD+DA(i))P=(A(i))TPDP+DPA(i)P, |
for all i.
Theorem 2.8. Let A={A(i)}ri=1 be a family of matrices such that A(i)∈Rn×n for i=1,…,r and P be a permutation matrix associated with some partition β. Then, the following are equivalent:
(i) A∈CLDSn.
(ii) {(PT(A(i)⋆Jm)P)∘S(i)}ri=1∈CLDSβm for all matrices S(i)⪰0 in Rm×m with S(i)[βj]≻0 for all i=1,…,r and j=1,…,n.
(iii) {(PT(A(i)⋆Jm)P)∘S(i)}ri=1∈CLDSβm for all matrices S(i)⪰0 in Rm×m with S(i)[βj]=Imj for all i=1,…,r and j=1,…,n.
(iv) {(PT(A(i)⋆Jm)P)∘S(i)}ri=1∈Pβm for all matrices S(i)⪰0 in Rm×m with S(i)[βj]≻0 for all i=1,…,r and j=1,…,n.
(v) {(PT(A(i)⋆Jm)P)∘S(i)}ri=1∈Pβm for all matrices S(i)⪰0 in Rm×m with S(i)[βj]=Imj for all i=1,…,r and j=1,…,n.
Proof. Clearly, the condition (ii) gives (iii) and (iv) gives (v). Moreover, (ii) leads to (iv) and (iii) to (v) by Lemma 2.4.
(i)⇒(ii): For i=1,…,r, let S(i)⪰0 in Rm×m be given as in (ii). Then, for each i, S(i)P⪰0 n by n block matrix with (S(i)P)jj≻0, j=1,…,n. Hence, it follows from Theorem 2.6 that {A(i)⋆S(i)P}ri=1∈CLDSαm. Now, by observing that
A(i)⋆S(i)P=(A(i)⋆Jm)∘S(i)P=((PT(A(i)⋆Jm)P)∘S(i))P | (2.3) |
for each i, we conclude that {((PT(A(i)⋆Jm)P)∘S(i))P}ri=1∈CLDSαm. Hence, by Lemma 2.7, (ii) follows.
(v)⇒(i): From observation 1, {((PT(A(i)⋆Jm)P)∘S(i))P}ri=1∈Pαm. This, by (2.3), means {A(i)⋆S(i)P}ri=1∈Pαm. Finally, using Theorem 2.6, we obtain (i). This completes the proof.
A final remark before moving to the next section is that Theorems 2.3, 2.6, and 2.8 are equivalent to Theorems 4, 9, and 10 in [29] when A is a singleton.
In this section, we generalize the main results of Section 2 to provide more characterizations for common Lyapunov α-scalar stability. In this section, let γ={γ1,…,γs} be any partition of {1,…,n}. Then, δ={δ1,…,δs}, where
δ1=|γ1|⋃i=1αi,δ2=|γ1|+|γ2|⋃i=|γ1|+1αi,…,δs=n⋃i=n−|γs|+1αi |
is a partition of {1,…,m}, where α={α1,…,αn} is the partition of {1,…,m} defined in Section 2.
Lemma 3.1. Let A={A(i)}ri=1 be a family of matrices such that A(i)∈Rn×n for i=1,…,r and γ={γ1,…,γs} be any partition of {1,…,n}. If {A(i)⋆Q}ri=1∈Pδm, then {T(A(i)⋆Q)}ri=1∈Pγn.
Proof. Let u(i)∈Rn, i=1,…,r, not all being zero. Then, let v=[vmj]∈Rm be the nonzero n-block vector defined as follows
(vmj)k={1if k=10if k≠1. |
Then, for each i, choose z(i)=u(i)⋆v. Clearly, z(i)∈Rm, i=1,…,r, and not all being zero vectors since not all u(i) are zero. Furthermore, we have T(z(i))=u(i) and T(z(i)[δl])=u(i)[γl] for l∈{1,…,s}. Now, because {A(i)⋆Q}ri=1∈Pδm, there is some l∈{1,…s} such that
r∑i=1z(i)[δl]T((A(i)⋆Q)z(i))[δl]>0 | (3.1) |
Now, observe that
r∑i=1z(i)[δl]T((A(i)⋆Q)z(i))[δl]=r∑i=1T(z(i)[δl]T)T(((A(i)⋆Q)z(i))[δl]). |
Consequently, it follows from (3.1) that
0<r∑i=1z(i)[δl]T((A(i)⋆Q)z(i))[δl]=r∑i=1T(z(i)[δl]T)T(((A(i)⋆Q)z(i))[δl])=r∑i=1u(i)[γl]T(T(A(i)⋆Q)T(z(i)))[γl]=r∑i=1u(i)[γl]T(T(A(i)⋆Q)u(i))[γl]. |
Therefore, by Definition 1.3, the result follows.
Lemma 3.2. ([28,Corollary 2.1]) Let A={A(i)}ri=1 be a family of matrices such that A(i)∈Rn×n for i=1,…,r and γ={γ1,…,γs} be a partition of {1,…,n}. Then, A∈CLDSγn if and only if there is l∈{1,…,s} such that
trr∑i=1(A(i)X(i))[γl]>0 |
for any X(i)⪰0, X(i)∈Rn×n, i=1,…,r, not all being zero matrices.
Theorem 3.3. Let A={A(i)}ri=1 be a family of matrices such that A(i)∈Rn×n for i=1,…,r and γ={γ1,…,γs} be any partition of {1,…,n}. Then, the following are equivalent:
(i) A∈CLDSγn.
(ii) {A(i)⋆S(i)}ri=1∈CLDSδm for all n by n block matrices S(i)⪰0 in Rm×m with S(i)jj≻0 for all i=1,…,r and j=1,…,n.
(iii) {A(i)⋆S(i)}ri=1∈CLDSδm for all n by n block matrices S(i)⪰0 in Rm×m with S(i)jj=Imj for all i=1,…,r and j=1,…,n.
(iv) {A(i)⋆S(i)}ri=1∈Pδm for all n by n block matrices S(i)⪰0 in Rm×m with S(i)jj≻0 for all i=1,…,r and j=1,…,n.
(v) {A(i)⋆S(i)}ri=1∈Pδm for all n by n block matrices S(i)⪰0 in Rm×m with S(i)jj=Imj for all i=1,…,r and j=1,…,n.
Proof. It is clear that (ii) implies (iii) and (iv) implies (v). Besides, from Lemma 2.4, it is clear that (ii) implies (iv) and (iii) implies (v).
(i)⇒(ii): Clearly, if D is a positive γ-scalar matrix, then D⋆Im is a positive δ-scalar matrix. Moreover, if D is a common Lyapunov γ-scalar solution for A, then, by Lemma 1.2, for any S(i)'s given as in (ii), we have
((A(i))TD+DA(i))⋆S(i)=(A(i)⋆S(i))T(D⋆Im)+(D⋆Im)(A(i)⋆S(i))≻0, |
for i=1,…,r. This last inequality means that {A(i)⋆S(i)}ri=1 has (D⋆Im) as a common Lyapunov δ-scalar solution.
(v)⇒(i): Let X(i)=[xkl]⪰0, X(i)∈Rn×n i=1,…,r, not all being zero. Next, let D(i)∈Rn×n, i=1,…,r, be a diagonal matrix such that d(i)kk=√x(i)kk for j=1,…,n. Thus, for each i, there is S(i)=[skl]⪰0, S(i)∈Rn×n with s(i)kk=1, k=1,…,n, such that X(i)=D(i)S(i)D(i). By setting p=max{m1,…,mn}, we conclude, using to Lemma 1.1, that S(i)⊗Ip⪰0 in Rnp×np, i=1,…,r. Since, for each i, S(i)⋆Q is a principal submatrix of S(i)⊗Ip, then S(i)⋆Q⪰0, Q here is the matrix in (1.3). Furthermore, each diagonal block (S(i)⋆Q)jj=Imj, i=1,…,r. So, {A(i)⋆(S(i)⋆Q)}ri=1={(A(i)∘S(i))⋆Q}ri=1∈Pδm, by (v). Consequently, according to Lemma 3.1, {T((A(i)∘S(i))⋆Q)}ri=1∈Pγn. Set u(i)=D(i)e, i=1,…,r, where e is the vector of all ones in Rn. By the construction of these u(i)'s, it is easy to see that not all of them are zero vectors. Therefore, there is some index q∈{1,…,s} such that
r∑i=1u(i)[γq]T((T((A(i)∘S(i))⋆Q))u(i))[γq]=r∑i=1(D(i)e)[γq]T((T((A(i)∘S(i))⋆Q))(D(i)e))[γq]=r∑i=1e[γq]TD(i)[γq]((T((A(i)∘(S(i)D(i)))⋆Q))e)[γq]=r∑i=1e[γq]T((T((A(i)∘(D(i)S(i)D(i)))⋆Q))e)[γq]=r∑i=1e[γq]T((T((A(i)∘X(i))⋆Q))e)[γq]=r∑i=1e[γq]T((A(i)∘X(i))e)[γq]=trr∑i=1(A(i)X(i))[γq]>0. |
Therefore, (i) follows by Lemma 3.2. This finishes the proof.
This last Theorem can be generalized to provide more characterizations for common Lyapunov α-scalar stability. Recall that in Section 2, we defined β to be a partition of {1,…,m} obtained from α through a permutation function τ. Using this notation and the definition of δ above, for any partition γ={γ1,…,γs} of {1,…,n}, we define another partition of {1,…,m} called ϵ={ϵ1.…,ϵs}, where
ϵ1=|γ1|⋃i=1βi,ϵ2=|γ1|+|γ2|⋃i=|γ1|+1βi,…,ϵs=n⋃i=n−|γs|+1βi. |
Clearly, if we replace α with δ and β with ϵ, Observation 1 will hold true for a permutation matrix P associated with β. Now, we have the following theorem, whose proof follows the lines of the proof of Theorem 2.8 and is therefore omitted.
Theorem 3.4. Let A={A(i)}ri=1 be a family of matrices such that A(i)∈Rn×n for i=1,…,r and γ={γ1,…,γs} be any partition of {1,…,n}. In addition, let P be a permutation matrix associated with some partition β. Then, the following are equivalent:
(i) A∈CLDSγn.
(ii) {(PT(A(i)⋆Jm)P)∘S(i)}ri=1∈CLDSϵm for all matrices S(i)⪰0 in Rm×m with S(i)[ϵj]≻0 for all i=1,…,r and j=1,…,n.
(iii) {(PT(A(i)⋆Jm)P)∘S(i)}ri=1∈CLDSϵm for all matrices S(i)⪰0 in Rm×m with S(i)[ϵj]=Imj for all i=1,…,r and j=1,…,n.
(iv) {(PT(A(i)⋆Jm)P)∘S(i)}ri=1∈Pϵm for all matrices S(i)⪰0 in Rm×m with S(i)[ϵj]≻0 for all i=1,…,r and j=1,…,n.
(v) {(PT(A(i)⋆Jm)P)∘S(i)}ri=1∈Pϵm for all matrices S(i)⪰0 in Rm×m with S(i)[ϵj]=Imj for all i=1,…,r and j=1,…,n.
We remark here that Theorems 3.3 and 3.4 are the same as Theorems 2.6 and 2.8, respectively, when γ={{1},{2},…,{n}}. Additionally, when r=1, i.e., A is a singleton, these last two theorems reduce to the following corollaries, whose proofs shall be omitted.
Corollary 3.5. Let A∈Rn×n and γ={γ1,…,γs} be any partition of {1,…,n}. Then, the following are equivalent:
(i) A∈LDSγn.
(ii) A⋆S∈LDSδm for all n by n block matrices S⪰0 in Rm×m with Sjj≻0, j=1,…,n.
(iii) A⋆S∈LDSδm for all n by n block matrices S⪰0 in Rm×m with Sjj=Imj, j=1,…,n.
(iv) A⋆S∈Pδm for all n by n block matrices S⪰0 in Rm×m with Sjj≻0, j=1,…,n.
(v) A⋆S∈Pδm for all n by n block matrices S⪰0 in Rm×m with Sjj=Imj, j=1,…,n.
Corollary 3.6. Let A∈Rn×n and γ={γ1,…,γs} be any partition of {1,…,n}. In addition, let P be a permutation matrix associated with some partition β. Then, the following are equivalent:
(i) A∈LDSγn.
(ii) (PT(A⋆Jm)P)∘S∈LDSϵm for all matrices S⪰0 in Rm×m with S[ϵj]≻0, j=1,…,n.
(iii) (PT(A⋆Jm)P)∘S∈LDSϵm for all matrices S⪰0 in Rm×m with S[ϵj]=Imj, j=1,…,n.
(iv) (PT(A⋆Jm)P)∘S∈Pϵm for all matrices S⪰0 in Rm×m with S[ϵj]≻0, j=1,…,n.
(v) (PT(A⋆Jm)P)∘S∈Pϵm for all matrices S⪰0 in Rm×m with S[ϵj]=Imj, j=1,…,n.
Motivated by the work in [29], we have presented new characterizations for common Lyapunov diagonal stability using the Khatri-Rao product. The notions of P-sets and Pα-sets have been used to formulate these results. Moreover, these characterizations have been extended to the notion of common Lyapunov α-scalar stability. Our work here extends and broadens the scope of results in [17,28].
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
The Researchers would like to thank the Deanship of Graduate Studies and Scientific Research at Qassim University for financial support (QU-APC-2024-9/1). We would like to extend my sincere gratitude to the anonymous reviewers for their valuable comments and suggestions, which have significantly contributed to improving the quality of this paper.
The author does not have any conflict of interest.
[1] |
E. Berchio, A. Falocchi, A positivity preserving property result for the biharmonic operator under partially hinged boundary conditions, Ann. Mat. Pur. Appl., 200 (2021), 1651–1681. https://doi.org/10.1007/s10231-020-01054-6 doi: 10.1007/s10231-020-01054-6
![]() |
[2] |
J. M. Ball, Stability theory for an extensible beam, J. Differ. Equations, 14 (1973), 399–418. https://doi.org/10.1016/0022-0396(73)90056-9 doi: 10.1016/0022-0396(73)90056-9
![]() |
[3] |
E. Emmrich, M. Thalhammer, A class of integro-differential equations incorporing nonlinear and nonlocal damping with applications in nonlinear elastodynamics: Existence via time discretization, Nonlinearity, 24 (2011), 2523–2546. https://doi.org/10.1088/0951-7715/24/9/008 doi: 10.1088/0951-7715/24/9/008
![]() |
[4] |
H. R. Clark, Elastic membrane equation in bounded and unbounded domains, Electron. J. Qual. Theo., 11 (2002), 1–21. https://doi.org/10.14232/ejqtde.2002.1.11 doi: 10.14232/ejqtde.2002.1.11
![]() |
[5] |
Y. You, Inertial manifolds and stabilization of nonlinear beam equations with Balakrishnan-Taylor damping, Abstr. Appl. Anal., 1 (1996), 83–102. https://doi.org/10.1155/S1085337596000048 doi: 10.1155/S1085337596000048
![]() |
[6] |
E. H. G. Tavares, M. A. J. Silva, V. Narciso, Long-time dynamics of Balakrishnan-Taylor extensible beams, J. Dyn. Differ. Equ., 32 (2020), 1157–1175. https://doi.org/10.1007/s10884-019-09766-x doi: 10.1007/s10884-019-09766-x
![]() |
[7] |
J. Glover, A. C. Lazer, P. J. Mckenna, Existence and stability of of large scale nonlinear oscillation in suspension bridges, Z. Angew. Math. Phys., 40 (1989), 172–200. https://doi.org/10.1007/BF00944997 doi: 10.1007/BF00944997
![]() |
[8] |
P. J. McKenna, W. Walter, Nonlinear oscillations in a suspension bridge, Arch. Ration. Mech. An., 98 (1987), 167–177. https://doi.org/10.1007/BF00251232 doi: 10.1007/BF00251232
![]() |
[9] |
A. Ferrero, F. Gazzola, A partially hinged rectangular plate as a model for suspension bridges, Discrete Cont. Dyn.-A, 35 (2015), 5879–5908. https://doi.org/10.3934/dcds.2015.35.587 doi: 10.3934/dcds.2015.35.587
![]() |
[10] | F. Gazzola, Mathematical models for suspension bridges: Nonlinear structural instability, modeling, simulation and applications, 1 Eds., New York: Springer-Verlag, 2015. https://doi.org/10.1007/978-3-319-15434-3 |
[11] |
V. S. Guliyev, M. N. Omarova, M. A. Ragusa, Characterizations for the genuine Calderon-Zygmund operators and commutators on generalized Orlicz-Morrey spaces, Adv. Nonlinear Anal., 12 (2023), 20220307. https://doi.org/10.1515/anona-2022-0307 doi: 10.1515/anona-2022-0307
![]() |
[12] | H. Y. Li, B. W. Feng, Exponential and polynomial decay rates of a porous elastic system with thermal damping, J. Funct. Space., 2023 (2023). https://doi.org/10.1155/2023/3116936 |
[13] | N. Taouaf, Global existence and exponential decay for thermoelastic system with nonlinear distributed delay, Filomat, 37 (2023), 8897–8908. |
[14] |
M. Al-Gwaiz, V. Benci, F. Gazzola, Bending and stretching energies in a rectangular plate modeling suspension bridges, Nonlinear Anal.-Theor., 106 (2014), 181–734. https://doi.org/10.1016/j.na.2014.04.011 doi: 10.1016/j.na.2014.04.011
![]() |
[15] |
G. Crasta, A. Falocchi, F. Gazzola, A new model for suspension bridges involving the convexification of the cables, Z. Angew. Math. Phys., 71 (2020), 93. https://doi.org/10.1007/s00033-020-01316-6 doi: 10.1007/s00033-020-01316-6
![]() |
[16] | S. A. Messaoudi, S. E. Mukiawa, A suspension bridge problem: Existence and stability, In: International Conference on Mathematics and Statistics, Cham: Springer International Publishing, 2017,151–165. https://doi.org/10.1007/978-3-319-46310-0_9 |
[17] |
M. M. Cavalcanti, W. J. Corrêa, R. Fukuoka, Z. Hajjej, Stabilization of a suspension bridge with locally distributed damping, Math. Control Signal., 30 (2018), 39. https://doi.org/10.1007/s00498-018-0226-0 doi: 10.1007/s00498-018-0226-0
![]() |
[18] |
A. D. D. Cavalcanti, M. Cavalcanti, W. J. Corrêa, Z. Hajjej, M. S. Cortés, R. V. Asem, Uniform decay rates for a suspension bridge with locally distributed nonlinear damping, J. Franklin I., 357 (2020), 2388–2419. https://doi.org/10.1016/j.jfranklin.2020.01.004 doi: 10.1016/j.jfranklin.2020.01.004
![]() |
[19] |
W. Liu, H. Zhuang, Global existence, asymptotic behavior and blow-up of solutions for a suspension bridge equation with nonlinear damping and source terms, Nonlinear Differ. Equ. Appl., 24 (2017), 67. https://doi.org/10.1007/s00030-017-0491-5 doi: 10.1007/s00030-017-0491-5
![]() |
[20] |
Y. Wang, Finite time blow-up and global solutions for fourth-order damped wave equations, J. Math. Anal. Appl., 418 (2014), 713–733. https://doi.org/10.1016/j.jmaa.2014.04.015 doi: 10.1016/j.jmaa.2014.04.015
![]() |
[21] |
S. A. Messaoudi, S. E. Mukiawa, Existence and stability of fourth-order nonlinear plate problem, Nonauton. Dyn. Syst., 6 (2019), 81–98. https://doi.org/10.1515/msds-2019-0006 doi: 10.1515/msds-2019-0006
![]() |
[22] |
Z. Hajjej, General decay of solutions for a viscoelastic suspension bridge with nonlinear damping and a source term, Z. Angew. Math. Phys., 72 (2021), 90. https://doi.org/10.1007/s00033-021-01526-6 doi: 10.1007/s00033-021-01526-6
![]() |
[23] |
E. Berchio, A. Falocchi, About symmetry in partially hinged composite plates, Appl. Math. Opt., 84 (2021), 2645–2669. https://doi.org/10.1007/s00245-020-09722-y doi: 10.1007/s00245-020-09722-y
![]() |
[24] |
E. Berchio, A. Falocchi, Maximizing the ratio of eigenvalues of nonhomogeneous partially hinged plates, J. Spectr. Theor., 11 (2021), 743–780. https://doi.org/10.4171/JST/355 doi: 10.4171/JST/355
![]() |
[25] |
D. Bonheure, F. Gazzola, I. Lasiecka, J. Webster, Long-time dynamics of a hinged-free plate driven by a nonconservative force, Ann. I. H. Poincaré-An., 39 (2022), 457–500. https://doi.org/10.4171/aihpc/13 doi: 10.4171/aihpc/13
![]() |
[26] | V. Komornik, Exact controllability and stabilization: The multiplier method, Paris: Masson-John Wiley, 1994. |
[27] |
V. Ferreira, F. Gazzola, E. M. dos Santos, Instability of modes in a partially hinged rectangular plate, J. Differ. Equations, 261 (2016), 6302–6340. https://doi.org/10.1016/j.jde.2016.08.037 doi: 10.1016/j.jde.2016.08.037
![]() |
[28] |
J. M. Ball, Initial-boundary value problems for an extensible beam, J. Math. Anal. Appl., 42 (1973), 61–90. https://doi.org/10.1016/0022-247X(73)90121-2 doi: 10.1016/0022-247X(73)90121-2
![]() |
[29] | J. L. Lions, Quelques methodes de resolution des problemes aux limites non lineaires, Paris: Dunod, 2002. |
[30] | M. T. L. Sonrier, Distrubutions espace de Sobolev application, Ellipses/Edition Marketing S.A, 1998. |
[31] |
S. Yayla, C. L. Cardozo, M. A. J. Silva, V. Narciso, Dynamics of a Cauchy problem related to extensible beams under nonlocal and localized damping effects, J. Math. Anal. Appl., 494 (2021), 124620. https://doi.org/10.1016/j.jmaa.2020.124620 doi: 10.1016/j.jmaa.2020.124620
![]() |
1. | Ali Algefary, Diagonal solutions for a class of linear matrix inequality, 2024, 9, 2473-6988, 26435, 10.3934/math.20241286 |