Advances in Matrices, Finite and Infinite, with ApplicationsView this Special Issue
Research Article | Open Access
Perturbation Analysis for the Matrix Equation
We consider the perturbation analysis of the matrix equation . Based on the matrix differentiation, we first give a precise perturbation bound for the positive definite solution. A numerical example is presented to illustrate the sharpness of the perturbation bound.
In this paper, we consider the matrix equation where are complex matrices, is an identity matrix, are nonnegative integers and the positive definite solution is practical interest. Here, and denote the conjugate transpose of the matrices and , respectively. Equation (1.1) arises in solving some nonlinear matrix equations with Newton method. See, for example, the nonlinear matrix equation which appears in Sakhnovich . Solving these nonlinear matrix equations gives rise to (1.1). On the other hand, (1.1) is the general case of the generalized Lyapunov equation whose positive definite solution is the controllability Gramian of the bilinear control system (see [2, 3] for more details) Set then (1.2) can be equivalently written as (1.1) with and .
Some special cases of (1.1) have been studied. Based on the kronecker product and fixed point theorem in partially ordered sets, Reurings  and Ran and Reurings [5, 6] gave some sufficient conditions for the existence of a unique positive definite solution of the linear matrix equations and . And the expressions for these unique positive definite solutions were also derived under some constraint conditions. For the general linear matrix equation (1.1), Reurings [, Page 61] pointed out that it is hard to find sufficient conditions for the existence of a positive definite solution, because the map is not monotone and does not map the set of positive definite matrices into itself. Recently, Berzig  overcame these difficulties by making use of Bhaskar-Lakshmikantham coupled fixed point theorem and gave a sufficient condition for (1.1) existing a unique positive definite solution. An iterative method was constructed to compute the unique positive definite solution, and the error estimation was given too.
Recently, the matrix equations of the form (1.1) have been studied by many authors (see [8–14]). Some numerical methods for solving the well-known Lyapunov equation , such as Bartels-Stewart method and Hessenberg-Schur method, have been proposed in . Based on the fixed point theorem, the sufficient and necessary conditions for the existence of a positive definite solution of the matrix equation have been given in [8, 9]. The fixed point iterative method and inversion-free iterative method were developed for solving the matrix equations in [13, 14]. By making use of the fixed point theorem of mixed monotone operator, the matrix equation was studied in , and derived a sufficient condition for the existence of a positive definite solution. Assume that maps positive definite matrices either into positive definite matrices or into negative definite matrices, the general nonlinear matrix equation was studied in , and the fixed point iterative method was constructed to compute the positive definite solution under some additional conditions.
Motivated by the works and applications in [2–6], we continue to study the matrix equation (1.1). Based on a new mathematical tool (i.e., the matrix differentiation), we firstly give a differential bound for the unique positive definite solution of (1.1), and then use it to derive a precise perturbation bound for the unique positive definite solution. A numerical example is used to show that the perturbation bound is very sharp.
Throughout this paper, we write if the matrix is positive definite (semidefinite). If is positive definite (semidefinite), then we write . If a positive definite matrix satisfies , we denote that . The symbols and denote the maximal and minimal eigenvalues of an Hermitian matrix , respectively. The symbol stands for the set of Hermitian matrices. The symbol denotes the spectral norm of the matrix .
2. Perturbation Analysis for the Matrix Equation (1.1)
Based on the matrix differentiation, we firstly give a differential bound for the unique positive definite solution of (1.1), and then use it to derive a precise perturbation bound for in this section.
Definition 2.1 (, Definition 3.6). Let , then the matrix differentiation of is . For example, let Then
Lemma 2.2 (, Theorem 3.2). The matrix differentiation has the following properties:(1); (2), where is a complex number;(3); (4); (5); (6), where is a constant matrix.
Theorem 2.4. If then (1.1) has a unique positive definite solution , and it satisfies
Proof. Since then consequently, Combining (2.5)–(2.9) we have Then by Lemma 2.3 we obtain that (1.1) has a unique positive definite solution , which satisfies Noting that is the unique positive definite solution of (1.1), then It is known that the elements of are differentiable functions of the elements of and . Differentiating (2.12), and by Lemma 2.2, we have which implies that By taking spectral norm for both sides of (2.14), we obtain that and noting (2.11) we obtain that Then Due to (2.5) we have Combining (2.17), (2.18) and noting (2.19), we have which implies that
Proof. By (2.22) and Theorem 2.4, we know that (1.1) has a unique positive definite solution . And by (2.23) we have similarly, by (2.24) we have By (2.28), (2.29), and Theorem 2.4 we obtain that the perturbed equation (2.25) has a unique positive definite solution .
From Theorem 2.4 it follows that Noting that and combining Mean Value Theorem of Integration, we have
3. Numerical Experiments
Example 3.1. Consider the symmetric linear matrix equation
and its perturbed equation
It is easy to verify that the conditions (2.22)–(2.24) are satisfied, then (3.1) and its perturbed equation (3.2) have unique positive definite solutions and , respectively. From Berzig  it follows that the sequences and generated by the iterative method both converge to . Choose as the termination scalar, that is, By using the iterative method (3.4) we can get the computed solution of (3.1). Since , then the computed solution has a very high precision. For simplicity, we write the computed solution as the unique positive definite solution . Similarly, we can also get the unique positive definite solution of the perturbed equation (3.2).
Some numerical results on the perturbation bounds for the unique positive definite solution are listed in Table 1.
From Table 1, we see that Theorem 2.5 gives a precise perturbation bound for the unique positive definite solution of (3.1).
In this paper, we study the matrix equation (1.1) which arises in solving some nonlinear matrix equations and the bilinear control system. A new method of perturbation analysis is developed for the matrix equation (1.1). By making use of the matrix differentiation and its elegant properties, we derive a precise perturbation bound for the unique positive definite solution of (1.1). A numerical example is presented to illustrate the sharpness of the perturbation bound.
This research was supported by the National Natural Science Foundation of China (11101100; 11171205; 11261014; 11226323), the Natural Science Foundation of Guangxi Province (2012GXNSFBA053006; 2011GXNSFA018138), the Key Project of Scientific Research Innovation Foundation of Shanghai Municipal Education Commission (13ZZ080), the Natural Science Foundation of Shanghai (11ZR1412500), the Ph.D. Programs Foundation of Ministry of Education of China (20093108110001), the Discipline Project at the corresponding level of Shanghai (A. 13-0101-12-005), and Shanghai Leading Academic Discipline Project (J50101). The authors wish to thank Professor Y. Zhang and the anonymous referees for providing very useful suggestions for improving this paper. The authors also thank Dr. Xiaoshan Chen for discussing the properties of matrix differentiation.
- L. A. Sakhnovich, Interpolation Theory and Its Applications, Mathematics and Its Applications, Kluwer Academic Publishers, Dordrecht, The Netherlands, 1997.
- T. Damm, “Direct methods and ADI-preconditioned Krylov subspace methods for generalized Lyapunov equations,” Numerical Linear Algebra with Applications, vol. 15, no. 9, pp. 853–871, 2008.
- L. Zhang and J. Lam, “On model reduction of bilinear systems,” Automatica, vol. 38, no. 2, pp. 205–216, 2002.
- M. C. B. Reurings, Symmetric matrix equations [Ph.D. thesis], Vrije Universiteit, Amsterdam, The Netherlands, 2003.
- A. C. M. Ran and M. C. B. Reurings, “The symmetric linear matrix equation,” Electronic Journal of Linear Algebra, vol. 9, pp. 93–107, 2002.
- 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.
- M. Berzig, “Solving a class of matrix equation via Bhaskar-Lakshmikantham coupled fixed point theorem,” Applied Mathematics Letters, vol. 25, no. 11, pp. 1638–1643, 2012.
- J. Cai and G. Chen, “Some investigation on Hermitian positive definite solutions of the matrix equation ,” Linear Algebra and Its Applications, vol. 430, no. 8-9, pp. 2448–2456, 2009.
- X. Duan and A. Liao, “On the existence of Hermitian positive definite solutions of the matrix equation ,” Linear Algebra and its Applications, vol. 429, no. 4, pp. 673–687, 2008.
- X. Duan, A. Liao, and B. Tang, “On the nonlinear matrix equation ,” Linear Algebra and its Applications, vol. 429, no. 1, pp. 110–121, 2008.
- 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.
- Z. Gajic and M. T. J. Qureshi, Lyapunov Matrix Equation in System Stability and Control, Academic Press, London, UK, 1995.
- V. I. Hasanov, “Notes on two perturbation estimates of the extreme solutions to the equations ,” Applied Mathematics and Computation, vol. 216, no. 5, pp. 1355–1362, 2010.
- Z.-Y. Peng, S. M. El-Sayed, and X.-l. Zhang, “Iterative methods for the extremal positive definite solution of the matrix equation ,” Journal of Computational and Applied Mathematics, vol. 200, no. 2, pp. 520–527, 2007.
- B. R. Fang, J. D. Zhou, and Y. M. Li, Matrix Theory, Tsinghua University Press, Beijing, China, 2006.
Copyright © 2012 Xue-Feng Duan and Qing-Wen Wang. 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.