An Inversion-Free Method for Finding Positive Definite Solution of a Rational Matrix Equation
A new iterative scheme has been constructed for finding minimal solution of a rational matrix equation of the form . The new method is inversion-free per computing step. The convergence of the method has been studied and tested via numerical experiments.
In this paper, we will discuss the following nonlinear matrix equation: where is an nonsingular complex matrix, is the unit matrix of the appropriate size, and is an unknown Hermitian positive definite (HPD) matrix that should be found. It was proved in  that if (1) has an HPD solution, then all its Hermitian solutions are positive definite and, moreover, it has the maximal solution and the minimal solution in the sense that for any HPD solution .
A lot of papers have been published regarding the iterative HPD solutions of such nonlinear rational matrix equations in the literature due to their importance in some practical problems arising in control theory, dynamical problems, and so forth (see [2, 3]).
The most common iterative method for finding the maximal solution of (1) is the following fixed-point iteration :
The maximal solution of (1) can be obtained through , where is the minimal solution of the dual equation .
In 2010, Monsalve and Raydan in  proposed the following iteration method (also known as Newton’s method) for finding the minimal solution: which is an inversion-free scheme. Since, should be computed only once in contrast to the matrix iteration (2). Note that , and similar notations are used throughout.
Remark 1. We remark that there are several other well-known iterative methods for solving (1) rather than Newton’s method (3). To the best of our knowledge, the procedure of extending higher-order iterative methods for finding the solution of (1) has not been exploited up to now. Hence, we hope that this interlink among the fields of root-finding and solving (1) may lead to discovering novel and innovative techniques.
The rest of this paper is organized as follows. In Section 2, we develop and analyze a new inversion-free method for finding roots of a special map . In Section 3, we provide some numerical comparisons by employing some experiments in machine precision. Some concluding remarks will be drawn in Section 4.
2. A New Iterative Method
An equivalent formulation of (1) is to find an HPD matrix such that . Toward this goal, we write . Furthermore, we have Now, using a change of variable as , we could simplify (4) as follows:
In order to obtain an iterative method for finding the minimal HPD solution of (1), it is now enough to solve the well-known matrix equation , at which . One of such ways to challenge this matrix inversion problem is via applying the Schulz-type iteration methods (see, e.g., [6–10]). Applying Chebyshev’s method  yields
We remark that there is a tight relationship between iterative methods for nonlinear systems and the construction of higher-order methods for matrix equations ([12, 13]).
The matrix iteration (6) requires to be computed only once at the beginning of the iteration and this makes the iterative method fall in the category of inversion-free algorithms for solving (1).
In the meantime, it is easy to show that the zeros of the map are equal to the zeros of the map , wherein . To be more precise, we are finding the inverse of the matrix which matches the minimal HPD solution of (1).
Remark 2. Following Remark 1, we applied Chebyshev’s method for (1) in this work and will study its theoretical behavior. The extension of the other well-known root-finding schemes for finding the minimal solution of (1) will remain for future studies.
Lemma 3. The proposed method (6) produces a sequence of Hermitian matrices using the Hermitian initial matrix .
Proof. The initial matrix is Hermitian, and . Thus, is also Hermitian; that is, . Now using inductive argument, we have By considering , () we now show that Note that has been used in (8). Now the conclusion holds for any . Thus, the proof is complete.
Theorem 4. By considering that and are nonsingular matrices, the sequence generated by (6) is convergent to the minimal solution using the initial matrix .
Proof. Let us consider . We therefore have
Taking a generic matrix operator norm from both sides of (9), we obtain
On the other hand, Chebyshev’s method for matrix inversion problem is convergent if the initial approximation reads . That is to say, . This together with the initial matrix gives which is true when is the minimal HPD solution of (1).
Note that since we obtain that ; thus . And subsequently using mathematical induction, it would be observed that tends to .
The only problem that happens in this process is the fact that the convergence order is -linear. In fact, although Chebyshev’s method for matrix inversion has third local order of convergence, this rate will not be preserved for finding the minimal HPD solution of (1).
The reason is that the matrix , which we must compute its inverse by Chebyshev’s method, is dependent on the itself. That is to say, the unknown is located in the essence of the matrix .
Theorem 5. The sequence of matrices produced by (6) satisfies the following error inequality: where and .
Proof. First since and by using (6) we have
Note that we have used the fact that . Relation (14) yields
wherein and + . We remark that .
Consequently, one has the error inequality (13). This shows the -linear order of convergence for finding the minimal HPD solution of (1). We thus have which is guaranteed since
3. Numerical Comparisons
In this section, we mainly investigate the performance of the new method (6) for matrix equation (1). All experiments were run on a Pentium IV computer, using Mathematica 8 . We report the number of required iterations (Iter) for converging. In our implementations, we stop all considered methods when the infinity norm of two successive iterates is less than given tolerance.
Note that recently Zhang in  studied a way to accelerate the beginning of such iterative methods for finding the minimal solution of (1) via applying multiple Newton’s method for matrix inversion. This technique could be given by for any . Subsequently, we could improve the behavior of the new method (6) using (18) as provided in Algorithm 1.
We compare Algorithm 1, denoted by PM, with (2) denoted by M1, (3) denoted by M2, and the method proposed by El-Sayed and Al-Dbiban  denoted by M3, which is a modification of the method presented by Zhan in , as follows:
Example 1 (see ). In this experiment, we compare the results of different methods for finding the minimal solution of (1) when the matrix is defined by
and the solution is
The results are given in Figure 1(a) in terms of the number of iterations when the stopping criterion is .
Note that PM and M2 converge to , whereas all the other schemes converge to ; thus we use other schemes to find the maximal solution of the dual equation in our written codes so as to have fair comparisons. Here has been chosen for PM (with ). This for the number of iterations in the inner finite loop of Algorithm 1 has been considered in the numerical report.
Furthermore, we have chosen this number empirically. In fact, varying shows us that we even can obtain better or worse results than the reported ones in different examples.
In Example 1, we have used . In fact we have chosen this value for since we are solving an operator equation in essence. To be more precise, we wish to consider the solution of the operator equation to be of multiplicity 2. This consideration makes the algorithm converge faster at the initial phase of the process and when we are enough close to the solution, then we flash back to the ordinary methods, that is, treat the solution as a simple zero (solution) of the operator equation.
Example 2 (see ). Applying the stopping criterion , we compare the behavior of various methods for the following test matrix: with the solution as The results are illustrated in Figure 1(b), wherein has been chosen for PM (with ).
We have studied the fact that the minimal HPD solution of (1) is equivalent to the roots of a nonlinear map. This special map has been solved by the well-known Chebyshev method as a matrix inversion problem.
The developed method requires the computation of one matrix inverse at the beginning of the process and it is hence an inversion-free method. The convergence and the rate of convergence have been studied for this scheme. Furthermore, using a proper acceleration technique from the literature, we have further speeded up the process of finding the HPD solution of (1).
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
The authors express their sincere thanks to the referees for the careful and details reading of the manuscript and very helpful suggestions that improved the manuscript substantially. The authors also gratefully acknowledge that this research was partially supported by the University Putra Malaysia under the GP-IBT Grant Scheme having project number GP-IBT/2013/9420100.
J. C. Engwerda, C. M. A. Ran, and A. L. Rijkeboer, “Necessary and sufficient conditions for the existence of a positive definite solution of the matrix equation ,” Linear Algebra and its Applications, vol. 186, pp. 255–275, 1993.View at: Google Scholar
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 Site | Google Scholar | MathSciNet
J. Li, “Solutions and improved perturbation analysis for the matrix equation ,” Abstract and Applied Analysis, vol. 2013, Article ID 575964, 12 pages, 2013.View at: Publisher Site | Google Scholar
J. C. Engwerda, “On the existence of a positive definite solution of the matrix equation ,” Linear Algebra and its Applications, vol. 194, pp. 91–108, 1993.View at: Publisher Site | Google Scholar | MathSciNet
M. Monsalve and M. Raydan, “A new inversion-free method for a rational matrix equation,” Linear Algebra and Its Applications, vol. 433, no. 1, pp. 64–71, 2010.View at: Publisher Site | Google Scholar | MathSciNet
M. Z. Ullah, F. Soleymani, and A. S. Al-Fhaid, “An efficient matrix iteration for computing weighted Moore-Penrose inverse,” Applied Mathematics and Computation, vol. 226, pp. 441–454, 2014.View at: Publisher Site | Google Scholar | MathSciNet
F. Soleymani, E. Tohidi, S. Shateyi, and F. K. Haghani, “Some matrix iterations for computing matrix sign function,” Journal of Applied Mathematics, vol. 2014, Article ID 425654, 9 pages, 2014.View at: Publisher Site | Google Scholar | MathSciNet
F. Soleymani, M. Sharifi, S. Shateyi, and F. Haghani, “A class of Steffensen-type iterative methods for nonlinear systems,” Journal of Applied Mathematics, vol. 2014, Article ID 705375, 9 pages, 2014.View at: Publisher Site | Google Scholar | MathSciNet
F. Soleymani, P. S. Stanimirović, S. Shateyi, and F. K. Haghani, “pproximating the matrix sign function using a novel iterative method,” Abstract and Applied Analysis, vol. 2014, Article ID 105301, 9 pages, 2014.View at: Publisher Site | Google Scholar | MathSciNet
F. Soleymani, M. Sharifi, S. Shateyi, and F. Khaksar Haghani, “An algorithm for computing geometric mean of two Hermitian positive definite matrices via matrix sign,” Abstract and Applied Analysis, vol. 2014, Article ID 978629, 6 pages, 2014.View at: Publisher Site | Google Scholar
J. F. Traub, Iterative Methods for the Solution of Equations, Prentice Hall, Englewood Cliffs, NJ, USA, 1964.View at: MathSciNet
A. Iliev and N. Kyurkchiev, Nontrivial Methods in Numerical Analysis: Selected Topics in Numerical Analysis, LAP LAMBERT Academic Publishing, 2010.
F. Soleymani, T. Lotfi, and P. Bakhtiari, “A multi-step class of iterative methods for nonlinear systems,” Optimization Letters, vol. 8, no. 3, pp. 1001–1015, 2014.View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
M. Trott, The Mathematica Guide-Book for Numerics, Springer, New York, NY, USA, 2006.
L. Zhang, “An improved inversion-free method for solving the matrix equation ,” Journal of Computational and Applied Mathematics, vol. 253, pp. 200–203, 2013.View at: Google Scholar
S. M. El-Sayed and A. M. Al-Dbiban, “A new inversion free iteration for solving the equation ,” Journal of Computational and Applied Mathematics, vol. 181, no. 1, pp. 148–156, 2005.View at: Publisher Site | Google Scholar | MathSciNet
X. Zhan, “Computing the extremal positive definite solutions of a matrix equation,” SIAM Journal on Scientific Computing, vol. 17, no. 5, pp. 1167–1174, 1996.View at: Publisher Site | Google Scholar | Zentralblatt MATH | 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 Site | Google Scholar | MathSciNet
S. M. El-Sayed, “An algorithm for computing positive definite solutions of the nonlinear matrix equation ,” International Journal of Computer Mathematics, vol. 80, no. 12, pp. 1527–1534, 2003.View at: Publisher Site | Google Scholar | MathSciNet