Journal of Applied Mathematics

Volume 2014 (2014), Article ID 128249, 8 pages

http://dx.doi.org/10.1155/2014/128249

## Notes on the Hermitian Positive Definite Solutions of a Matrix Equation

^{1}School of Mathematics and Statistics, Shandong University, Weihai 264209, China^{2}School of Mathematics, Shandong University, Jinan 250100, China

Received 15 January 2014; Revised 26 March 2014; Accepted 15 April 2014; Published 6 May 2014

Academic Editor: Alexander Timokha

Copyright © 2014 Jing Li and Yuhai Zhang. 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.

#### Abstract

The nonlinear matrix equation, with is investigated. A fixed point theorem in partially ordered sets is proved. And then, by means of this fixed point theorem, the existence of a unique Hermitian positive definite solution for the matrix equation is derived. Some properties of the unique Hermitian positive definite solution are obtained. A residual bound of an approximate solution to the equation is evaluated. The theoretical results are illustrated by numerical examples.

#### 1. Introduction

In this paper we consider the Hermitian positive definite solution of the following nonlinear matrix equation: where , are complex matrices, is an Hermitian positive definite matrix, and is a positive integer. Here, denotes the conjugate transpose of the matrix .

This type of nonlinear matrix equations arises in many practical applications. Equation (1) in the case comes from ladder networks, dynamic programming, control theory, stochastic filtering, statistics, and so forth [1–7]. When , (1) in the case is recognized as playing an important role in modeling certain optimal interpolation problems (see [8, 9] for more details). The Hermitian positive definite solutions of the general-type equation (1) play an important role in connection with a certain system of linear equations in many physical calculations (see [7, 10] for more details). When solving the nonlinear matrix equation (1), we often do not avoid some round-off errors. Then we only get an approximation . After we computed , we would like to know how good our computation was. Motivated by these, we consider in this paper the Hermitian positive definite solutions and the residual bound of (1).

In the last few years, (1) was investigated in some special cases. For the nonlinear matrix equations, [11–17], [18], [19, 20], [21], [22–25], and [9, 26, 27], there were many contributions in the literature to the solvability, numerical solutions, and perturbation analysis. In addition, the related general equations, [28–34], were studied by some scholars.

For the case and changes with , Duan et al. [35] proved that (1) with has a unique Hermitian positive definite solution by fixed point theorems for monotone and mixed monotone operators in a normal cone. Lim [36] showed that (1) with has a unique Hermitian positive definite solution by using a strict contraction for the Thompson metric on the open convex cone of positive definite matrices. Shi et al. [37] studied the existence and uniqueness of solutions of nonlinear matrix equations, and , with . Li [10] gave perturbation analysis for the positive definite solution of (1) with . Duan et al. [38] gave two perturbation estimates for the positive definite solution of (1) with . However, these papers have not considered (1) in the case for some and the approaches in these papers will become invalid in this case. Meanwhile, in some practical problems (e.g., in certain optimal interpolation problems), the case of for some is required. To our best knowledge, there has been no literature paying attention to the Hermitian positive definite solutions and the residual bound for (1) with . By using the integral representation of matrix function and the fixed point theorem, we prove the existence of a unique Hermitian positive definite solution to (1) and consider the residual bound of this equation. Note that the integral representation of matrix function in the case is different from the case . Furthermore, the monotonicity of in the former case differs from the latter. Based on the above arguments, we will consider (1) with in this paper.

The rest of the paper is organized as follows. In Section 2, we give some preliminary lemmas that will be needed to develop this work. In Section 3, a fixed point theorem in partially ordered sets is proved. And then, by means of this fixed point theorem, the existence of a unique Hermitian positive definite solution for the matrix equation (1) with is derived. We propose an iterative method to compute the Hermitian positive definite solution. We also obtain some properties of the unique Hermitian positive definite solution. Furthermore, in Section 4, a residual bound for the unique Hermitian positive definite solution to (1) with is given. Finally, several numerical examples are presented in Section 5.

We denote by the set of complex matrices, by the set of Hermitian matrices, by the identity matrix, by the spectral norm, and by and the maximal and minimal eigenvalues of , respectively. For , we write (resp., if is Hermitian positive semidefinite (resp., definite). Further, the sets and are defined by and , respectively.

#### 2. Preliminaries

Lemma 1 (see [39]). *If , then .*

Lemma 2 (see [39]). *If and , then .*

For the sake of completeness we will provide the proof of the next lemma.

Lemma 3 (see [17]). *For every positive definite matrix , if , with , then
*

*Proof. *Suppose that is a positive definite matrix. If , then

*Lemma 4 (see [23]). For every Hermitian positive definite matrix and , it yields that(i);
(ii), . In addition, if and , with , then
*

*3. The Positive Definite Solutions*

*3. The Positive Definite Solutions*

*In this section, we use a new method, which is different from the approaches applied in [35, 36] to prove that (1) with has a unique Hermitian positive definite solution . Meanwhile, we give an iterative method to compute the unique Hermitian positive definite solution for arbitrary initial positive definite matrix. Moreover, we obtain some properties of the Hermitian positive definite solution to (1).*

*Theorem 5. Let be a partially ordered metric space, with the property that, for any two elements and in , there is a positive number such that .Let be a continuous, order reversing map such that there is a with for all . Put . Then maps into itself.Consider . Assume in addition that for there is a number such that, for all ,
Then has a unique fixed point in , and, for every , the iteration , started with , converges to the unique fixed point.*

*Proof. *To show that maps into itself, we only need to show that for we have . In fact, since is order reversing, something much stronger holds; for we have . In particular, maps into . Moreover, maps into . It follows that if there is a fixed point of , then it is in .

The fact that is order reversing means that is order preserving and one can check that there are two matrices, and , such that these two form a periodic orbit which is the attractor of the iteration of for any starting value. In addition, we have and , so that .

It remains to show that, under the extra condition (5), and are equal. In fact, we will show that , which is enough. By the assumption on , there is a such that . Let . Suppose that . Then, using (5),
which contradicts the definition of . So , and in particular .

*The following results are immediate consequences of Theorem 5.*

*Theorem 6. There exists a unique Hermitian positive definite solution of (1) with , and the iteration,
converges to .*

*Proof. *Let . Then the set is partially ordered and for any two elements, and , in there is a positive number such that . A map associated with (1) is defined by
Obviously, is continuous, and a solution of (1) is a fixed point of . Let and . By Lemmas 1 and 2, we obtain that is order reversing. So . Then
and as also . That is, maps into itself and is also order reversing.

By Theorem 5, it remains to prove that, for , there exists a number such that, for all ,
In fact, choose . Note that
For , a calculation gives
which completes the proof.

*Theorem 7. If is an Hermitian positive definite solution of (1), then , where and are, respectively, the solutions of the following equations:
Moreover,
*

*Proof. *Define the sequences and as
From (16), it follows that
Supposing and , then
Hence, for each , we have and , which imply that the sequences and are monotonic and bounded. Therefore, they are convergent to certain positive numbers. Let
Taking limits in (16) yields
which imply that
Therefore and satisfy (13) and (14), respectively. We will prove that for any positive definite solution . According to Lemmas 1 and 2 and the sequences defined by (16), it follows that for each Hermitian positive definite solution . From , it follows that . Hence
Since , it follows that and . Note that inequality (22) implies . By similar induction, it yields that
Taking limits on both sides of inequality (23), we have .

*Corollary 8. Every Hermitian positive definite solution of (1) is in , where and are defined as in Theorem 7.*

*Proof. *We suppose that is the Hermitian positive definite solution of (1). By Theorem 7, it follows that
Using , we obtain . Applying inequality (24) yields .

*Remark 9. *The above estimate of Hermitian positive definite solution of (1) is more precise than that in Theorem 7.

*4. Residual Bound*

*4. Residual Bound*

*The matrix equation (1) with includes the following cases: first, some in and others ; secondly, all in ; thirdly, all . Without loss of generality, let , , and , . Then (1) can be rewritten as
In this section, a residual bound of an approximate solution for the unique solution to (25) is developed.*

*Theorem 10. Let be an approximation to the solution of (25). If
and the residual satisfies
then
*

*Proof. *Let
where . Obviously, is a nonempty bounded convex closed set. Let
Evidently is continuous. We will prove that . For every , we have
Hence
That is,
Using (27), one sees that
which means that .

*According to Lemmas 3 and 4 and inequality (33), we obtain
By Brouwer's fixed point theorem, there exists a such that . Hence is a solution of (25). Moreover, by Theorem 6, we know that the solution of (25) is unique. Then
*

*5. Numerical Examples*

*5. Numerical Examples**To illustrate the results of the previous sections, in this section, two simple examples are given, which were carried out using MATLAB 7.1. For the stopping criterion we take .*

*Example 1. *In this example, we study the following matrix equation:
with
Algorithm (7) needs 10 iterations to obtain the unique positive definite solution
with the residual .

*Example 2. *In this example, we consider the corresponding perturbation bound for the solution in Theorem 10. We consider the following matrix equation:
with
Choose . Let the approximate solution of be given with the iterative method (7), where is the iterative number.

*The residual satisfies the conditions in Theorem 10. By Theorem 10, we can compute residual bounds for as
where
Some results are listed in Table 1.*

*The results listed in Table 1 show that the residual bound given by Theorem 10 is fairly sharp.*

*Conflict of Interests*

*Conflict of Interests**The authors declare that they have no competing interests.*

*Authors’ Contribution*

*Authors’ Contribution**All the authors carried out the proof. All the authors conceived the study and participated in its design and coordination. All the authors read and approved the final paper.*

*Acknowledgments*

*Acknowledgments**The authors would like to express their gratitude to the referees for their fruitful comments, which have led to the present form of Theorems 5 and 6. The work was supported in part by the National Nature Science Foundation of China (11201263), Natural Science Foundation of Shandong Province (ZR2012AQ004), and Independent Innovation Foundation of Shandong University (IIFSDU), China.*

*References*

*References*

- J. Zabczyk, “Remarks on the control of discrete-time distributed parameter systems,”
*SIAM Journal on Control and Optimization*, vol. 12, pp. 721–735, 1974. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - W. N. Anderson, Jr., G. B. Kleindorfer, P. R. Kleindorfer, and M. B. Woodroofe, “Consistent estimates of the parameters of a linear system,”
*Annals of Mathematical Statistics*, vol. 40, pp. 2064–2075, 1969. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - R. S. Bucy, “A priori bounds for the Riccati equation,” in
*Proceedings of the Sixth Berkeley Symposium on Mathematical Statistics and Probability (Univ. California, Berkeley, Calif., 1970/1971), Vol. III: Probability Theory*, pp. 645–656, Univ. California Press, Berkeley, Calif, USA, 1972. View at MathSciNet - D. V. Ouellette, “Schur complements and statistics,”
*Linear Algebra and Its Applications*, vol. 36, pp. 187–295, 1981. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - W. N. Anderson, Jr., T. D. Morley, and G. E. Trapp, “The cascade limit, the shorted operator and quadratic optimal control,” in
*Linear Circuits, Systems and Signal Processing: Theory and Application (Phoenix, AZ, 1987)*, C. I. Byrnes, C. F. Martin, and R. E. Saeks, Eds., pp. 3–7, North-Holland, Amsterdam, The Netherlands, 1988. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - W. Pusz and S. L. Woronowicz, “Functional calculus for sesquilinear forms and the purification map,”
*Reports on Mathematical Physics*, vol. 8, no. 2, pp. 159–170, 1975. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - B. L. Buzbee, G. H. Golub, and C. W. Nielson, “On direct methods for solving Poisson's equations,”
*SIAM Journal on Numerical Analysis*, vol. 7, pp. 627–656, 1970. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - L. A. Sakhnovich,
*Interpolation Theory and Its Applications*, vol. 428 of*Mathematics and Its Applications*, Kluwer Academic Publishers, Dordrecht, The Netherlands, 1997. View at Publisher · View at Google Scholar · View at MathSciNet - X. Duan, Q. Wang, and A. Liao, “On the matrix equation arising in an interpolation problem,”
*Linear and Multilinear Algebra*, vol. 61, no. 9, pp. 1192–1205, 2013. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J. Li, “Perturbation analysis of the nonlinear matrix equation $X-{\sum}_{i=1}^{m}{A}_{i}^{*}{X}^{{p}_{i}}{A}_{i}=Q$,”
*Abstract and Applied Analysis*, vol. 2013, Article ID 979832, 11 pages, 2013. View at Publisher · View at Google Scholar · View at MathSciNet - J. C. Engwerda, A. C. M. Ran, and A. L. Rijkeboer, “Necessary and sufficient conditions for the existence of a positive definite solution of the matrix equation $X+{A}^{*}{X}^{-1}A=Q$,”
*Linear Algebra and Its Applications*, vol. 186, pp. 255–275, 1993. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J. C. Engwerda, “On the existence of a positive definite solution of the matrix equation $X+{A}^{T}{X}^{-1}A=I$,”
*Linear Algebra and Its Applications*, vol. 194, pp. 91–108, 1993. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - C.-H. Guo and P. Lancaster, “Iterative solution of two matrix equations,”
*Mathematics of Computation*, vol. 68, no. 228, pp. 1589–1603, 1999. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A. Ferrante and B. C. Levy, “Hermitian solutions of the equation $X=Q+{NX}^{-1}{N}^{*}$,”
*Linear Algebra and Its Applications*, vol. 247, pp. 359–373, 1996. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - V. I. Hasanov and I. G. Ivanov, “On two perturbation estimates of the extreme solutions to the equations $X\pm {A}^{*}{X}^{-1}A=Q$,”
*Linear Algebra and Its Applications*, vol. 413, no. 1, pp. 81–92, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - V. I. Hasanov, I. G. Ivanov, and F. Uhlig, “Improved perturbation estimates for the matrix equations $X\pm {A}^{*}{X}^{-1}A=Q$,”
*Linear Algebra and Its Applications*, vol. 379, pp. 113–135, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J. Li and Y. H. Zhang, “The Hermitian positive definite solution and its perturbation analysis for the matrix equation $X-{A}^{*}{X}^{-1}A=Q$,”
*Mathematica Numerica Sinica*, vol. 30, no. 2, pp. 129–142, 2008 (Chinese). View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. Vaezzadeh, S. M. Vaezpour, R. Saadati, and C. Park, “The iterative methods for solving nonlinear matrix equation $X+{A}^{*}{X}^{-1}A+{B}^{*}{X}^{-1}B=Q$,”
*Advances in Difference Equations*, vol. 2013, article 229, 10 pages, 2013. View at Publisher · View at Google Scholar · View at MathSciNet - I. G. Ivanov, V. I. Hasanov, and B. V. Minchev, “On matrix equations $X\pm {A}^{*}{X}^{-2}A=I$,”
*Linear Algebra and Its Applications*, vol. 326, no. 1-3, pp. 27–44, 2001. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Y.-H. Zhang, “On Hermitian positive definite solutions of matrix equation $X-{A}^{*}{X}^{-2}A=I$,”
*Journal of Computational Mathematics*, vol. 23, no. 4, pp. 408–418, 2005. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - X.-G. Liu and H. Gao, “On the positive definite solutions of the matrix equations ${X}^{s}\pm {A}^{T}{X}^{-t}A={I}_{n}$,”
*Linear Algebra and Its Applications*, vol. 368, pp. 83–97, 2003. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - V. I. Hasanov, “Positive definite solutions of the matrix equations $X\pm {A}^{*}{X}^{-q}A=Q$,”
*Linear Algebra and Its Applications*, vol. 404, pp. 166–182, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J. Li and Y. Zhang, “Perturbation analysis of the matrix equation $X-{A}^{*}{X}^{-p}A=Q$,”
*Linear Algebra and Its Applications*, vol. 431, no. 9, pp. 1489–1501, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - H. Xiao and J. T. Wang, “On the matrix equation $X-{A}^{*}{X}^{-p}A=Q(p>1)$,”
*Chinese Journal of Engineering Mathematics*, vol. 26, no. 2, pp. 305–309, 2009. View at Google Scholar · View at MathSciNet - J. Li, “Solutions and improved perturbation analysis for the matrix equation $X-{A}^{*}{X}^{-p}A=Q(p>0)$,”
*Abstract and Applied Analysis*, vol. 2013, Article ID 575964, 12 pages, 2013. View at Publisher · View at Google Scholar · View at MathSciNet - X. Yin and L. Fang, “Perturbation analysis for the positive definite solution of the nonlinear matrix equation $X-{\sum}_{i=1}^{m}{A}_{i}^{*}{X}^{-1}{A}_{i}=Q$,”
*Journal of Applied Mathematics and Computing*, vol. 43, no. 1-2, pp. 199–211, 2013. View at Publisher · View at Google Scholar · View at MathSciNet - J. Li and Y. Zhang, “Sensitivity analysis of the matrix equation from interpolation problems,”
*Journal of Applied Mathematics*, vol. 2013, Article ID 518269, 8 pages, 2013. View at Publisher · View at Google Scholar · View at MathSciNet - S. M. El-Sayed and A. C. M. Ran, “On an iteration method for solving a class of nonlinear matrix equations,”
*SIAM Journal on Matrix Analysis and Applications*, vol. 23, no. 3, pp. 632–645, 2001/02. View at Publisher · View at Google Scholar · View at MathSciNet - A. C. M. Ran and M. C. B. Reurings, “On the nonlinear matrix equation $X+{A}^{*}\mathcal{F}\left(X\right)A=Q$: solutions and perturbation theory,”
*Linear Algebra and Its Applications*, vol. 346, pp. 15–26, 2002. View at Publisher · View at Google Scholar · View at MathSciNet - A. C. M. Ran and M. C. B. Reurings, “A fixed point theorem in partially ordered sets and some applications to matrix equations,”
*Proceedings of the American Mathematical Society*, vol. 132, no. 5, pp. 1435–1443, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A. C. M. Ran, M. C. B. Reurings, and L. Rodman, “A perturbation analysis for nonlinear selfadjoint operator equations,”
*SIAM Journal on Matrix Analysis and Applications*, vol. 28, no. 1, pp. 89–104, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. C. B. Reurings, “Contractive maps on normed linear spaces and their applications to nonlinear matrix equations,”
*Linear Algebra and Its Applications*, vol. 418, no. 1, pp. 292–311, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. C. B. Reurings,
*Symmetric matrix equations [Ph.D. thesis]*, VU University, Amsterdam, The Netherlands, 2003. - D. Zhou, G. Chen, G. Wu, and X. Zhang, “On the nonlinear matrix equation ${X}^{s}+{A}^{*}F\left(X\right)A=Q$ with $s\ge 1$,”
*Journal of Computational Mathematics*, vol. 31, no. 2, pp. 209–220, 2013. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - X. Duan, A. Liao, and B. Tang, “On the nonlinear matrix equation $X-{\sum}_{i=1}^{m}{A}_{i}^{*}{X}^{{\delta}_{i}}{A}_{i}=Q$,”
*Linear Algebra and Its Applications*, vol. 429, no. 1, pp. 110–121, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Y. Lim, “Solving the nonlinear matrix equation $X=Q+{\sum}_{i=1}^{m}{M}_{i}{X}^{{\delta}_{i}}{M}_{i}^{*}$ via a contraction principle,”
*Linear Algebra and Its Applications*, vol. 430, no. 4, pp. 1380–1383, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - X. Shi, F. Liu, H. Umoh, and F. Gibson, “Two kinds of nonlinear matrix equations and their corresponding matrix sequences,”
*Linear and Multilinear Algebra*, vol. 52, no. 1, pp. 1–15, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - X.-F. Duan, Q.-W. Wang, and C.-M. Li, “Perturbation analysis for the positive definite solution of the nonlinear matrix equation $X-{\sum}_{i=1}^{m}{A}_{i}^{*}{X}^{{\delta}_{i}}{A}_{i}=Q$,”
*Journal of Applied Mathematics & Informatics*, vol. 30, no. 3-4, pp. 655–663, 2012. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - R. Bhatia,
*Matrix Analysis*, vol. 169 of*Graduate Texts in Mathematics*, Springer, New York, NY, USA, 1997. View at Publisher · View at Google Scholar · View at MathSciNet

*
*