Research Article | Open Access

# Norm Estimates for Solutions of Polynomial Operator Equations

**Academic Editor:**William E. Fitzgibbon

#### Abstract

We consider the equations and , where , , are given linear bounded operators in a Banach space and is to be found. Representations of solutions are derived. In the cases of Euclidean and Hilbert spaces, norm estimates for the solutions are suggested.

#### 1. Introduction and Statement of the Main Result

The objects of this paper are the equationswhere , , , , are given linear bounded operators in a Banach space with a norm and is to be found.

There is a long tradition of finding different expressions for the solution of an operator equation in the form of operator integrals and series, some prominent examples of which occur in the works of E. Heinz, M. G. Krein, M. Rosenblum and R. Bhatia, C. Davis, and A. McIntosh. A comprehensive summary of these is contained in [1] and [2]. This tradition was continued in [3], which deals with the equationClearly, (3) is a particular case of (2).

The special case of (2) is the much studied Sylvester equation, of great interest in operator theory, numerical analysis, and engineering; compare [4–6] and references therein. Let be the operator adjoint to . The equationis the famous Lyapunov equation playing an essential role in the stability theory of differential equations. The equationis an example of (1); in particular, the discrete Lyapunov equationis an important tool in the theory of difference equations.

About other operator equations see [7–9] and references therein. In particular, [9] deals with necessary and sufficient conditions for the existence of solution to the systems of the general solution to a system of adjointable operator equations over Hilbert modules. In [8] the nonlinear operator equations of the form and are considered. It should be noted that in the finite dimensional case the numerical methods are well developed. Besides, the traditional methods convert matrix equations into their equivalent forms by using the Kronecker product, which involve the inversion of the associated large matrix and result in increasing computation and excessive computer memory. The recently suggested gradient based iterative methods are more powerful; compare [10–13] and references therein.

The aim of the present paper is to derive representations of solutions to (1) and (2) and to estimate the norms of these solutions. Such estimates are important, in particular, for the investigations of linear and perturbations of nonlinear differential and difference equations. Norm estimates for solutions of the Sylvester equation whose coefficients are normal operators can be found in [2]. In the finite dimensional case solution estimates for (4) and (6) have been established in [14].

In this paper, in particular, we considerably generalize the main results from [14].

Denote by and the spectrum and spectral radius of , respectively.

Theorem 1. *Let be the roots taken with the multiplicities of the polynomialand letThen (1) has a unique solution which can be represented by and the series strongly converges.*

The proof of this theorem is presented in the next section.

#### 2. Proof of Theorem 1

Following [1, Section I.3], introduce the operators and by and , respectively. Besides, and commute.

Let be a scalar function regular on . Define the operator valued functionwhere , are closed Jordan contours surrounding and , respectively. So for any bounded linear operator Besides,for functions , regular on . Take . Ifthen due to formula from [1, Section I.3] (1) can be written asSincemaking use of (13), we can writeBut . Thus,and therefore we get the following result.

Lemma 2. *Under condition (14) equation (1) has a unique solution defined by*

Furthermore, for , satisfyingconsider the operatorThe series converges andWe thus get the following result.

Lemma 3. *Let condition (20) hold. Then operator defined by (21) is the unique solution to the equation and .*

*Proof of Theorem 1. *PutLemma 2 implies . Due to the previous lemma are solutions to the equationsprovided (9) holds. SoContinuing this process for , according to Lemma 2, we prove the theorem.

#### 3. Representations of Solutions to (2)

Assume that is invertible; then from (2) we haveLet denote the lower spectral radius of : . Taking into account the fact that and applying now Theorem 1 to (26) with instead of , we get the following result.

Corollary 4. *LetThen (2) has a unique solution which can be represented byand the series strongly converges.*

Condition (27) does not allow us to consider the Lyapunov equation, since . Because of this we are going to derive the representation of solutions to (2) under other conditions. To this end putAccording to (19) for a solution to (26) we haveprovidedHere is defined as . ThusNow the invertibility condition for can be removed. We thus have proved the following result.

Lemma 5. *Let the conditionhold. Then (2) has a unique solution which can be represented by (32).*

Lemma 6. *Let the conditionhold. Then (4) has a unique solution , which can be represented as*

*Proof. *For the brevity put and . Then and Moreover, . Equation (4) is equivalent to the following one:Due to Theorem 9.2 from [2] a solution of (36) is defined by the equalityas claimed.

Corollary 7. *Let the conditionhold. Then the equationhas a unique solution , which can be represented as*

A solution of (39) is also given by . So under condition (38) we haveNow assume thatThen by (41)Let . Then . So condition (42) means that .

PutThen . Due to the previous lemma are solutions to the equationsprovided (42) holds. SoContinuing this process for , according to Lemma 5, we obtainWe arrive at the following result.

Theorem 8. *Let be the roots of and let condition (42) hold. Then (2) has a unique solution , which can be represented by (47).*

#### 4. Solution Estimates in the Finite Dimensional Case

##### 4.1. Equation (1)

In this section is a Euclidean space. Let be the eigenvalues of a matrix counted with their multiplicities. The following quantity (the departure from normality of ) plays a key role hereafter:where is the Frobenius (Hilbert-Schmidt norm) of .

The following relations are checked in [15, Section 2.1]:where . If is a normal matrix: , then . If and are commuting matrices, then . By the inequality between geometric and arithmetic mean values we haveSo . Obviously .

By Corollary 2.7.2 from [15], one hasNote that . Furthermore, due to Example 1.10.3 from [15] we have

From (21) and Lemma 3 we haveprovided condition (20) holds. Hence, and from (51), it directly followswhere

Now Lemma 2 and (54) imply the following result.

Corollary 9. *Let condition (9) hold. Then a unique solution to (1) in satisfies the inequality*

If is normal, then andButThus, we haveprovided is normal. If both and are normal, thenFurthermore, obviously,Takingwe getwhereNow Corollary 9 implies the following result.

Corollary 10. *Let condition (9) hold. Then a unique solution of (1) in satisfies the inequality*

Theorem 1 and simple calculations imply the following result.

Corollary 11. *Let . Then (7) has a unique solution , which satisfies the inequalities*

From Corollary 11 and (63) it follows that

##### 4.2. Equation (2)

Due to (41) and (52) under (38) we havewhere . ButThus whereprovided . Due to Lemma 5 we arrive at the following result.

Corollary 12. *Let condition (42) hold. Then a unique solution to (2) in satisfies the inequality*

If is normal, then and thereforeIf both and are normal, then

For the Sylvester equation we have the following result.

Corollary 13. *Let . Then (4) in has a unique solution , which satisfies the inequality *

Consider the Lyapunov equation (5). Taking in Corollary 13 instead of and instead of , since and , we get the following result.

Corollary 14. *Let . Then (5) in has a unique solution , which satisfies the inequality *

#### 5. Solution Estimates in the Infinite-Dimensional Case

##### 5.1. Equation (1)

In this section , a separable Hilbert space. It is assumed thatPutRecall that and are the eigenvalues of a matrix counted with their multiplicities. If is a Hilbert-Schmidt operator, then , wherecompare Lemma 6.5.2 [15]. Due to Example 7.10.2 from [15] we can writeAs it was shown in Section 4provided (9) holds. Hence, and from (80), it directly follows thatwhereNow Lemma 2 implies the following result.

Corollary 15. *Let conditions (9) and (77) hold. Then the unique solution to (1) in satisfies the inequality*

If is normal, then and according to (59)If both and are normal, thenFurthermore, according to (63), we get , where where is defined by (62). Now Corollary 15 implies the following result.

Corollary 16. *Let conditions (9) and (77) hold. Then a unique solution of (1) in satisfies the inequality*

By virtue of Corollary 15 we can assert that a unique solution of (7) satisfies the inequalityprovided satisfies conditions (77) and . From Corollary 16 it follows that

##### 5.2. Equation (2)

Due to Example 7.10.3 from [15] we haveprovided satisfies (77). Now let us consider (2), assuming that condition (42) holds. Due to (41) and (91),where . Thuswhereprovided . Due to Lemma 5 we arrive at the following result.

Corollary 17. *Let conditions (42) and (77) hold. Then the unique solution to (2) in satisfies the inequality*

Let . Then, if is normal, then and thereforeIf both and are normal, thenAccording to Corollary 17 we get the following result.

Corollary 18. *Let satisfies conditions (77) and . Then the solution to the Lyapunov equation (5) in satisfies the inequality*

#### Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.

#### References

- Y. L. Daleckii and M. G. Krein,
*Stability of Solutions of Differential Equations in Banach Space*, American Mathematical Society, Providence, RI, USA, 1971. - R. Bhatia and P. Rosenthal, “How and why to solve the operator equation $\text{AX}-\text{XB}=Y$,”
*The Bulletin of the London Mathematical Society*, vol. 29, no. 1, pp. 1–21, 1997. View at: Publisher Site | Google Scholar | MathSciNet - R. Bhatia and M. Uchiyama, “The operator equation ${\sum}_{i=1}^{n}{A}^{n-1}X{B}^{i}=Y$,”
*Expositiones Mathematicae*, vol. 27, no. 3, pp. 251–255, 2009. View at: Publisher Site | Google Scholar - M. Dehghan and M. Hajarian, “The reflexive and anti-reflexive solutions of a linear matrix equation and systems of matrix equations,”
*Rocky Mountain Journal of Mathematics*, vol. 40, no. 3, pp. 825–848, 2010. View at: Publisher Site | Google Scholar | MathSciNet - M. Konstantinov, D.-W. Gu, V. Mehrmann, and P. Petkov,
*Perturbation Theory for Matrix Equations*, vol. 9 of*Studies in Computational Mathematics*, North-Holland, Amsterdam, The Netherlands, 2003. - A. G. Mazko, “Matrix equations, spectral problems and stability of dynamic systems,” in
*Stability, Oscillations and Optimization of Systems*, Cambridge Scientific Publishers, 2008. View at: Google Scholar - D. S. Djordjević, “Explicit solution of the operator equation ${A}^{*}X+{X}^{*}A=B$,”
*Journal of Computational and Applied Mathematics*, vol. 200, no. 2, pp. 701–704, 2007. View at: Publisher Site | Google Scholar - B. P. Duggal, “Operator equations $ABA={A}^{2}$ and $BAB={B}^{2}$,”
*Functional Analysis, Approximation and Computation*, vol. 3, no. 1, pp. 9–18, 2011. View at: Google Scholar | MathSciNet - Q.-W. Wang and C.-Z. Dong, “The general solution to a system of adjointable operator equations over Hilbert
*C*^{∗}modules,”*Operators and Matrices*, vol. 5, no. 2, pp. 333–350, 2011. View at: Publisher Site | Google Scholar | MathSciNet - B. Zhou, G.-R. Duan, and Z.-Y. Li, “Gradient based iterative algorithm for solving coupled matrix equations,”
*Systems & Control Letters*, vol. 58, no. 5, pp. 327–333, 2009. View at: Publisher Site | Google Scholar - F. Ding and T. Chen, “Gradient based iterative algorithms for solving a class of matrix equations,”
*IEEE Transactions on Automatic Control*, vol. 50, no. 8, pp. 1216–1221, 2005. View at: Publisher Site | Google Scholar | MathSciNet - L. Xie, J. Ding, and F. Ding, “Gradient based iterative solutions for general linear matrix equations,”
*Computers & Mathematics with Applications*, vol. 58, no. 7, pp. 1441–1448, 2009. View at: Publisher Site | Google Scholar | MathSciNet - F. P. A. Beik, D. K. Salkuyeh, and M. M. Moghadam, “Gradient-based iterative algorithm for solving the generalized coupled Sylvester-transpose and conjugate matrix equations over reflexive (anti-reflexive) matrices,”
*Transactions of the Institute of Measurement and Control*, vol. 36, no. 1, pp. 99–110, 2014. View at: Publisher Site | Google Scholar - M. I. Gil', “Norm estimates for solutions of matrix equations $AX-XB=C$ and $X-AXB=C$,”
*Discussiones Mathematicae: Differential Inclusions, Control & Optimization*, vol. 34, no. 2, pp. 191–206, 2014. View at: Google Scholar - M. I. Gil',
*Operator Functions and Localization of Spectra*, vol. 1830 of*Lecture Notes in Mathematics*, Springer, Berlin, Germany, 2003.

#### Copyright

Copyright © 2015 Michael Gil’. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.