- Research
- Open access
- Published:
On nonlinear matrix equations \(X\pm\sum_{i=1}^{m}A_{i}^{*}X^{-n_{i}}A_{i}=I\)
Journal of Inequalities and Applications volume 2015, Article number: 147 (2015)
Abstract
We study the nonlinear matrix equations \(X+\sum_{i=1}^{m}A_{i}^{*}X^{-n_{i}}A_{i}=I\) and \(X-\sum_{i=1}^{m}A_{i}^{*}X^{-n_{i}}A_{i}=I\), where \(n_{i}\) are positive integers for \(i=1,2,\ldots,m\). The iterative algorithms for obtaining positive definite solutions for these equations are proposed. The necessary and sufficient conditions for the existence of positive definite solutions of these equations are derived. Moreover, the rate of convergence of the sequences generated from the algorithms is studied. The efficiency of proposed algorithms is illustrated by numerical examples.
1 Introduction
Consider the nonlinear matrix equations:
and
where X is an unknown square matrix, I is the identity matrix, \(A_{i}\) are square complex matrices and \(n_{i}\) are positive integers for \(i=1,2,\ldots,m\).
Nonlinear matrix equations of type (1.1) and (1.2) have many applications in engineering, control theory, dynamic programming, stochastic filtering, ladder networks, statistics, etc.; see [1–4] and the references therein. When \(m=1\) and \(n_{1}=1\), (1.2) arises in the analysis of stationary Gaussian reciprocal processes over a finite interval [3]. When \(m>1\) and \(n_{i}=1\) for \(i=1,2,\ldots,m\), (1.1) arises in solving a large-scale system of linear equations in many physical calculations [5] and (1.2) is recognized as playing an important role in modeling certain optimal interpolation problems [6, 7].
In the last few years, many authors have been greatly interested in developing the theory and numerical approaches for positive definite solutions to the nonlinear matrix equations of the form (1.1) and (1.2). Similar types of (1.1) and (1.2) have been investigated [8–12]. The matrix equations \(X\pm A^{*}X^{-1}A=Q\) have been studied by several authors [1–4, 13, 14] and different iterative algorithms for computing the positive definite solutions with linear and quadratic rate of convergence are proposed. Ivanov et al. [15] derived sufficient conditions for the existence of positive definite solutions for the matrix equations \(X\pm A^{*}X^{-2}A=I\) and they proposed iterative algorithms for obtaining positive definite solutions of these equations. El-Sayed [16] presented two iterative methods for calculating the positive definite solutions of the matrix equation \(X-A^{*}X^{-n}A=Q\), for the integer \(n\geq1\), the first method is derived for a normal matrix A and for the second method a sufficient condition for convergence is given for \(n = 2^{k}\). El-Sayed and Ran [17] studied the general matrix equation \(X+ A^{*}F(X)A=Q\) where F maps positive definite matrices either into positive definite matrices or into negative definite matrices and satisfies some monotonicity property. Hasanov and Ivanov [18] considered the matrix equations \(X\pm A^{*}X^{-n}A=Q\), they studied the solutions and perturbation analysis of these solutions. They also derived a sufficient condition for the existence of a unique positive definite solution of the equation \(X- A^{*}X^{-n}A=Q\). Hasanov [19] established and proved theorems for the necessary and sufficient conditions of the existence of positive definite solutions for the matrix equations \(X\pm A^{*}X^{-q}A=Q\) with \(0< q\leq1\), he showed that the equation \(X- A^{*}X^{-q}A=Q\) has a unique positive definite solution by using the properties of matrix sequence in Banach space. Also, in [5] some conditions for the existence of positive definite solution of the equation \(X+\sum_{i=1}^{m}A_{i}^{*}X^{-1}A_{i}=I\) have been obtained and two iterative algorithms to find the maximal positive definite solution of this equation have been presented. Duan et al. [20] gave two perturbation estimates for the positive definite solution of the equation \(X-\sum_{i=1}^{m}A_{i}^{*}X^{\delta_{i}}A_{i}=Q\) with \(0 < |\delta_{i}|< 1\). Duan et al. [6] studied the equation \(X-\sum_{i=1}^{m}N_{i}^{*}X^{-1}N_{i}=I\), they used the Thompson metric to prove that the matrix equation always has a unique positive definite solution and they derived a precise perturbation bound for the unique positive definite solution. In addition, other nonlinear matrix equations such as \(X^{s}\pm A^{T}X^{-t}A=I_{n}\) [21], \(AX^{2} + BX+C = 0\) [22], and \(X=Q +A^{H}(I\otimes X-C)^{-\delta}A^{*}\) [23] have been investigated.
In this paper, we study the positive definite solutions of (1.1) and (1.2). We derive the necessary and sufficient conditions for the existence of positive definite solutions. We suggest iterative algorithms for obtaining positive definite solutions of these equations. Moreover, under some conditions we obtain the rates of convergence of the iterative sequences of approximate solutions and the stopping criterions. Finally, we give some numerical examples to ensure the performance and the effectiveness of the suggested iterative algorithms.
The following notations will be used in this paper. \(A^{\star}\) denotes the complex conjugate transpose of A. We write \(A > 0\) (\(A\geq0\)), if matrix A is positive definite (positive semidefinite). If \(A-B\) is positive definite (positive semidefinite), then we write \(A > B\) (\(A\geq B\)). Moreover, we denote \(\rho(A)\) by the spectral radius of A. We use \(\lambda_{\mathrm{max}}(A)\) and \(\lambda _{\mathrm{min}}(A)\) to denote the maximal and minimal eigenvalues of A. \(\| \cdot \|\) and \(\| \cdot \|_{\infty}\) denote the spectral and infinity norm, respectively.
Lemma 1.1
[24]
If \(A\geq B> 0\), then \(A^{-1}\leq B^{-1}\).
Lemma 1.2
[24]
If A and B are positive definite matrices for which \(A-B>0\) and \(AB=BA\) are satisfied, then \(A^{n}-B^{n}>0\).
Lemma 1.3
[25]
If \(A>B>0\) (or \(A\geq B>0\)), then \(A^{\alpha}>B^{\alpha}>0\) (or \(A^{\alpha}\geq B^{\alpha}>0\)), for all \(\alpha\in(0,1]\), and \(0< A^{\alpha}<B^{\alpha}\) (or \(0< A^{\alpha }\leq B^{\alpha}\)), for all \(\alpha\in[-1,0)\).
2 The matrix equation \(X+\sum_{i=1}^{m}A_{i}^{*}X^{-n_{i}}A_{i}=I\)
In this section, we give some necessary and sufficient conditions for the existence of positive definite solutions of (1.1). We present the following iterative algorithm to compute the positive definite solution of (1.1).
Algorithm 2.1
Remark 2.1
Letting \(m=1\) in Algorithm 2.1, we get Algorithm (2.2) in [10] which is proposed for obtaining the positive definite solutions of the matrix equation \(X+A^{*}X^{-n}A=I\). Also, letting \(n_{i}=1\), \(\forall i=1,2,\ldots,m\), in Algorithm 2.1, we get Algorithm 2.7 in [5], which is proposed for obtaining the positive definite solutions of the matrix equation \(X+\sum_{i=1}^{m}A_{i}^{*}X^{-1}A_{i}=I\).
The following theorem provides the necessary condition for the existence of positive definite solutions of (1.1).
Theorem 2.1
If (1.1) has a positive definite solution X, then
Proof
Since X is a positive definite solution of (1.1), then \(X\leq I\) and \(\sum_{i=1}^{m}A_{i}^{*}X^{-n_{i}}A_{i}< I\). Using the inequality \(X\leq I\) and Lemmas 1.1, 1.2, we have
Also, from the inequality \(\sum_{i=1}^{m}A_{i}^{*}X^{-n_{i}}A_{i}< I\), we have \(A_{i}^{*}X^{-n_{i}}A_{i}< I\). Then
which implies that
Using Lemma 1.3, we obtain
This completes the proof. □
Remark 2.2
Letting \(n_{i}=1\), \(\forall i=1,2,\ldots,m\), in (2.1) we get the condition \(A_{i}A^{*}_{i}< X \leq I-\sum_{i=1}^{m}A_{i}^{*}A_{i}\), which is necessary for the existence of positive definite solutions of the matrix equation \(X+\sum_{i=1}^{m}A_{i}^{*}X^{-1}A_{i}=I\) ([5], Theorem 2.1).
Lemma 2.1
If \(A_{i}\), \(i=1,2,\ldots,m\), are hermitian matrices and \(A_{i}A_{j}=A_{j}A_{i}\), for all \(i,j=1,2,\ldots,m\), then
where the sequence \(\{X_{s}\}\), \(s=0,1,2,\ldots\) , is determined by Algorithm 2.1.
Proof
Since \(X_{0}=I\), \(A_{j}X_{0}=X_{0}A_{j}\). By using the condition \(A_{i}A_{j}=A_{j}A_{i}\), we have
We suppose that \(A_{j}X_{s}=X_{s}A_{j}\). Then for \(X_{s+1}\), we have
Hence, the equalities (2.2) are true, for all \(s=0,1,2,\ldots \) . □
Lemma 2.2
If \(A_{i}\), \(i=1,2,\ldots,m\), are hermitian matrices and \(A_{i}A_{j}=A_{j}A_{i}\), for all \(i,j=1,2,\ldots,m\), then
Here the sequences \(\{X_{s}\}\), \(\{X_{r}\}\), \(s,r=0,1,2,\ldots\) , are determined by Algorithm 2.1.
Proof
Since \(X_{0}=I\), \(X_{0}X_{r}=X_{r}X_{0}\), \(\forall r=0,1,2,\ldots\) . According to Lemma 2.1, we have
That is, \(X_{1}X_{r}=X_{r}X_{1}\), \(\forall r=0,1,2,\ldots\) . We suppose that \(X_{s}X_{r}=X_{r}X_{s}\), \(\forall r=0,1,2,\ldots\) . Then for \(X_{s+1}\), we have
Therefore, the equality (2.3) is true, for all \(s,r=0,1,2,\ldots \) . □
Remark 2.3
When we compare Lemmas 2.1 and 2.2 by Lemmas 4 and 5 in [10], we note that the sequence \(\{X_{s}\}\), \(s=0,1,2,\ldots\) (which is defined by Algorithm 2.1) satisfies the same properties of the sequence \(\{X_{s}\}\), \(s=0,1,2,\ldots \) (which is defined by Algorithm (2.2) in [10]).
The following theorem provides the sufficient condition for the existence of positive definite solutions of (1.1).
Theorem 2.2
Let \(A_{i}\), \(i=1,2,\ldots,m\), be hermitian matrices and \(A_{i}A_{j}=A_{j}A_{i}\), for all \(i,j=1,2,\ldots,m\). If \(A^{2}_{i}\leq\frac{(\alpha-1)}{nm\alpha^{(2n_{i}+1)}}I\), where \(\alpha>1\) and \(n=\max_{1 \leq i \leq m}\{n_{i}\}\), then (1.1) has a positive definite solution.
Proof
We consider the sequence \(\{X_{s}\}\) generated from Algorithm 2.1. For \(X_{0}\), we have \(X_{0}=I>\frac{1}{\alpha}I\). For \(X_{1}\), we have
That is,
We suppose that
Using the inequalities (2.4) and Lemmas 1.1, 1.2, and 2.2, we obtain
Also
Therefore the inequalities (2.4) are true, for all \(s=1,2,\ldots \) . That is, the sequence \(\{X_{s}\}\) is monotonically decreasing and bounded below by \(\frac{1}{\alpha}I\). Hence, the sequence \(\{X_{s}\}\) converges to a positive definite solution X of (1.1). □
Theorem 2.3
Let \(A_{i}\), \(i=1,2,\ldots,m\), be hermitian matrices and \(A_{i}A_{j}=A_{j}A_{i}\), for all \(i,j=1,2,\ldots,m\). If \(A^{2}_{i}\leq\frac{(\alpha-1)}{nm\alpha^{(2n_{i}+1)}}I\), then (1.1) has a positive definite solution X which satisfies
where \(\alpha>1\), \(n=\max_{1 \leq i \leq m}\{n_{i}\}\), and \(\{ X_{s}\}\), \(s=0,1,2,\ldots\) , is the sequence determined by Algorithm 2.1.
Proof
By Theorem 2.2, we know that the sequence \(\{X_{s}\}\), \(s=0,1,2,\ldots \) , is convergent to a positive definite solution X of (1.1). We consider the spectral norm of the matrix \(X_{s+1}-X\). We have
From the proof of Theorem 2.2, we obtain \(X^{-n_{i}}_{s}<\alpha ^{n_{i}}I\), \(X^{-n_{i}} \leq\alpha^{n_{i}}I\), and \(X\leq X_{s}\leq I\). Then we have
□
Theorem 2.4
Let \(A_{i}\), \(i=1,2,\ldots,m\), be hermitian matrices and \(A_{i}A_{j}=A_{j}A_{i}\), for all \(i,j=1,2,\ldots,m\). If \(A^{2}_{i}\leq\frac{(\alpha-1)}{nm\alpha^{(2n_{i}+1)}}I\), where \(\alpha>1\), \(n=\max_{1 \leq i \leq m}\{n_{i}\}\), and after s iterative steps of Algorithm 2.1, we have \(\| I-X^{n_{i}}_{s}X^{-n_{i}}_{s-1}\|<\varepsilon\), then
Proof
From Algorithm 2.1, we have
By taking the norm on both sides of the above equation, we have
From the proof of Theorem 2.2, we have \(X_{s}^{-n_{i}}<\alpha ^{n_{i}}I\), then
□
3 The matrix equation \(X-\sum_{i=1}^{m}A_{i}^{*}X^{-n_{i}}A_{i}=I\)
In this section, we give some necessary and sufficient conditions for the existence of positive definite solutions of (1.2). We present the following iterative algorithm to compute the positive definite solution of (1.2).
Algorithm 3.1
Remark 3.1
Letting \(m=1\) in Algorithm 3.1, we get Algorithm (2.1) in [16] which is proposed for obtaining the positive definite solutions of the matrix equation \(X-A^{*}X^{-n}A=I\). Also, letting \(n_{i}=1\), \(\forall i=1,2,\ldots,m\), in Algorithm 3.1, we get Algorithm (2.1) in [6] which is proposed for obtaining the positive definite solutions of the matrix equation \(X-\sum_{i=1}^{m}A_{i}^{*}X^{-1}A_{i}=I\).
The following theorem provides the necessary condition for the existence of positive definite solutions of (1.2).
Theorem 3.1
If (1.2) has a positive definite solution X, then
Proof
Since X is a positive definite solution of (1.2), then \(\sum_{i=1}^{m}A_{i}^{*}X^{-n_{i}}A_{i}\geq0\). Thus we get
Also, from the inequality \(X\geq I\) and Lemmas 1.1 and 1.2, we have
This completes the proof. □
Remark 3.2
The condition (3.1) is the same necessary condition for the existence of positive definite solutions of the matrix equation \(X-\sum_{i=1}^{m}A_{i}^{*}X^{-1}A_{i}=I\) ([6], Remark 2.1). Also, if A is an invertible matrix and \(m=1\) in condition (3.1), then we get the condition \(I< X < I+A^{*}A\), which is necessary for the existence of positive definite solutions of the matrix equation \(X-A^{*}X^{-n}A=I\) ([16], Corollary 2.1).
Lemma 3.1
If \(A_{i}\), \(i=1,2,\ldots,m\), are hermitian matrices and \(A_{i}A_{j}=A_{j}A_{i}\), for all \(i,j=1,2,\ldots,m\), then
where the sequence \(\{X_{s}\}\), \(s=0,1,2,\ldots\) , is determined by Algorithm 3.1.
Proof
The proof is similar to the proof of Lemma 2.1. □
Lemma 3.2
If \(A_{i}\), \(i=1,2,\ldots,m\), are hermitian matrices and \(A_{i}A_{j}=A_{j}A_{i}\), for all \(i,j=1,2,\ldots,m\), then
Here the sequences \(\{X_{s}\}\), \(\{X_{r}\}\), \(s,r=0,1,2,\ldots\) , are determined by Algorithm 3.1.
Proof
The proof is similar to the proof of Lemma 2.2. □
Remark 3.3
When we compare Lemmas 3.1 and 3.2 by Lemmas 2.3 and 2.4 in [16], we note that the sequence \(\{X_{s}\} \), \(s=0,1,2,\ldots\) (which is defined by Algorithm 3.1) satisfies the same properties of the sequence \(\{X_{s}\}\), \(s=0,1,2,\ldots\) (which is defined by Algorithm (2.1) in [16]).
Theorem 3.2
Let \(A_{i}\), \(i=1,2,\ldots,m\), be hermitian matrices and \(A_{i}A_{j}=A_{j}A_{i}\), for all \(i,j=1,2,\ldots,m\). If \(q=\sum_{i=1}^{m}n_{i}\|A_{i}\|^{2} (1+\sum_{i=1}^{m}\|A_{i}\| ^{2} )^{n_{i}-1}<1\), then (1.2) has a positive definite solution X which satisfies
and
where the sequence \(\{X_{s}\}\), \(s=0,1,2,\ldots\) , is determined by Algorithm 3.1.
Proof
We consider the matrix sequence \(\{X_{s}\}\) generated from Algorithm 3.1 and using Lemmas 1.1, 1.2 and 3.2. Since \(A_{i}^{2}\geq0\), then
and
Consequently
We find the relation between \(X_{2}\), \(X_{3}\), \(X_{4}\), \(X_{5}\). Using \(X_{0}\leq X_{2}\leq X_{1}\), we obtain
and
Hence \(X_{2}\leq X_{3}\leq X_{1}\). In the same way we can prove that
We suppose that
Using the inequalities (3.6), we have
Similarly
Therefore, the inequalities (3.6) are true, for all \(s=0,1,2,\ldots \) . Consequently the subsequences \(\{X_{2s}\}\) and \(\{ X_{2s+1}\}\) are convergent to positive definite matrices. To prove that these sequences have a common limit, we consider
From the inequalities (3.6), we have \(X_{2s}, X_{2s-1}\geq I\), and \(X_{2s}, X_{2s-1}\leq I+\sum_{i=1}^{m}A_{i}^{2}\), for all \(s=1,2,3,\ldots \) . Then we have
Hence
that is,
Hence, the two subsequences \(\{X_{2s}\}\) and \(\{X_{2s+1}\}\) have the same limit X, which is a positive definite solution of (1.2). □
From Theorem 3.2, we can deduce the following corollary.
Corollary 3.1
From inequality (3.5), we have the following upper bound:
Remark 3.4
Theorem 3.2 provides the sufficient condition \(q=\sum_{i=1}^{m}n_{i}\|A_{i}\|^{2} (1+ \sum_{i=1}^{m}\|A_{i}\| ^{2} )^{n_{i}-1}<1\) for the existence of positive definite solutions of (1.2), we note that when \(m=1\) we have the condition \(\|A\|^{2}(1+\|A\|^{2})^{n-1}<\frac{1}{n}\), which is sufficient for the existence of positive definite solutions of the matrix equation \(X-A^{*}X^{-n}A=I\) ([16], Theorem 2.1).
Theorem 3.3
Let \(A_{i}\), \(i=1,2,\ldots,m\), be hermitian matrices and \(A_{i}A_{j}=A_{j}A_{i}\), for all \(i,j=1,2,\ldots,m\). If \(q=\sum_{i=1}^{m}n_{i}\|A_{i}\|^{2} (1+\sum_{i=1}^{m}\|A_{i}\|^{2} )^{n_{i}-1}<1\), and after s iterative steps of Algorithm 3.1, we have \(\| I-X^{n_{i}}_{s-1}X^{-n_{i}}_{s}\|<\varepsilon\), then
Proof
From Algorithm 3.1, we have
By taking the norm on both sides of the above equation, we have
From the proof of Theorem 3.2, we have \(X_{s-1}^{-n_{i}}\leq I\), then
□
4 Numerical examples
In this section, we use the iterative Algorithms 2.1 and 3.1 to compute the positive definite solutions of (1.1) and (1.2), respectively. The solutions are computed for different matrices \(A_{i}\), \(i=1,2,\ldots,m\), with different orders. We denote X, the solution obtained by Algorithms 2.1 and 3.1 and \(\epsilon(X_{s})=\|X-X_{s}\|_{\infty}\), \(R_{1}(X_{s})=\|X_{s}+\sum_{i=1}^{m}A_{i}^{*}X_{s}^{-n_{i}}A_{i}-I\|_{\infty}\), \(Y_{s}= I-\sum_{i=1}^{m}A_{i}^{*}A_{i}-X_{s}\), \(Z_{i,s}= X^{n_{i}}_{s}-A_{i}A^{*}_{i}\) (\(i=1,2,\ldots,m\)), \(R_{2}(X_{s})=\Vert X_{s}-\sum_{i=1}^{m}A_{i}^{*}X_{s}^{-n_{i}}A_{i}-I \Vert _{\infty}\).
Example 4.1
Consider the matrix equation
where \(A_{1}\), \(A_{2}\), and \(A_{3} \) are given by
We use Algorithm 2.1 to solve (4.1). After 22 iterations, we get the positive definite solution
and \(R_{1}(X_{22})= 1.11022\times10^{-16}\), \(\lambda_{\mathrm{min}}(Y_{22})= 0.000488382\), \(\lambda_{\mathrm{min}}(Z_{1,22})=0.825427\), \(\lambda_{\mathrm{min}}(Z_{2,22})=0.778478\), \(\lambda _{\mathrm{min}}(Z_{3,22})=0.860326\).
The other results are listed in Table 1.
Example 4.2
Consider the matrix equation
where \(A_{1}\), \(A_{2}\), \(A_{3}\), and \(A_{4} \) are given by
We use Algorithm 2.1 to solve (4.2). After 29 iterations, we get the positive definite solution
and \(R_{1}(X_{29})= 1.11022\times10^{-16}\), \(\lambda_{\mathrm{min}}(Y_{29})= 0.00144519\), \(\lambda_{\mathrm{min}}(Z_{1,29})=0.649354\), \(\lambda_{\mathrm{min}}(Z_{2,29})=0.878163\), \(\lambda_{\mathrm{min}}(Z_{3,29})= 0.517277\), \(\lambda_{\mathrm{min}}(Z_{4,29})= 0.785219\).
The other results are listed in Table 2.
Example 4.3
Consider the matrix equation
where \(A_{1}\) and \(A_{2}\) are given as in Example 3.1 from [5]:
We use Algorithm 2.1 to solve (4.3). After 78 iterations, we get the positive definite solution
and \(R_{1}(X_{78})= 1.11022\times10^{-16}\), \(\lambda_{\mathrm{min}}(Y_{78})= 0.001965\), \(\lambda_{\mathrm{min}}(Z_{1,78})=0.479668\), \(\lambda_{\mathrm{min}}(Z_{2,78})=0.427158\).
The other results are listed in Table 3.
Example 4.4
Consider the matrix equation
where \(A_{1}\), \(A_{2}\), \(A_{3}\), and \(A_{4}\) are given by
We use Algorithm 3.1 to solve (4.4). After 133 iterations, we get the positive definite solution
and \(R_{2}(X_{133})= 2.22045\times10^{-16}\), \(\lambda _{\mathrm{min}}(I+\sum_{i=1}^{4}A_{i}^{*}A_{i}-X_{133})= 0.490411\), \(\lambda_{\mathrm{min}}(X_{133}-I)=0.100459\).
The other results are listed in Table 4.
Example 4.5
Consider the matrix equation
where \(A_{1}\), \(A_{2}\), and \(A_{3} \) are given by
We use Algorithm 3.1 to solve (4.5). After 78 iterations, we get the positive definite solution
and \(R_{2}(X_{78})= 2.22045\times10^{-16}\), \(\lambda _{\mathrm{min}}(I+\sum_{i=1}^{3}A_{i}^{*}A_{i}-X_{78})= 0.0223096\), \(\lambda_{\mathrm{min}}(X_{78}-I)=0.0237773\).
The other results are listed in Table 5.
Example 4.6
Consider the matrix equation
where \(A_{1}\) and \(A_{2} \) are given as in Example 4.1 from [6]:
We use Algorithm 3.1 to solve (4.6). After 792 iterations, we get the positive definite solution
and \(R_{2}(X_{792})= 7.99361\times10^{-15}\), \(\lambda _{\mathrm{min}}(I+\sum_{i=1}^{2}A_{i}^{*}A_{i}-X_{792})= 0.183803\), \(\lambda_{\mathrm{min}}(X_{792}-I)=0.134276\).
The other results are listed in Table 6.
5 Conclusion
In this paper, we investigate the nonlinear matrix equations \(X\pm\sum_{i=1}^{m}A_{i}^{*}X^{-n_{i}}A_{i}=I\), where \(n_{i}\), \(i=1,2,\ldots,m\), are positive integers. Necessary and sufficient conditions for the existence of positive definite solutions are derived. Iterative algorithms are proposed to compute the positive definite solutions of these equations. Moreover, some numerical examples are given to illustrate the effectiveness and rapidly convergence rate (small run time) of the proposed iterative algorithms (see values of \(\epsilon(X_{s})\), \(R_{1}(X_{s})\), and \(R_{2}(X_{s})\)). Also, the values of \(\lambda _{\mathrm{min}}\) show that the solutions of the matrix equations satisfy the necessary conditions.
References
Anderson, WN, Morley, TD, Trapp, GE: Positive solutions to \(X=A-BX^{-1}B^{*}\). Linear Algebra Appl. 134, 53-62 (1990)
Engwerda, JC: On the existence of a positive definite solution of the matrix equation \(X + A^{T}X^{-1}A = I\). Linear Algebra Appl. 194, 91-108 (1993)
Ferrante, A, Levy, BC: Hermitian solutions of the equation \(X = Q + NX^{-1}N^{*}\). Linear Algebra Appl. 247, 359-373 (1996)
Zhan, X: Computing the extremal positive definite solutions of a matrix equation. SIAM J. Sci. Comput. 17, 1167-1174 (1996)
He, YM, Long, JH: On the Hermitian positive definite solution of the nonlinear matrix equation \(X+\sum_{i=1}^{m}A_{i}^{*}X^{-1}A_{i}=I\). Appl. Math. Comput. 216, 3480-3485 (2010)
Duan, X, Wang, Q, Liao, A: On the matrix equation \(X-\sum_{i=1}^{m}N_{i}^{*}X^{-1}N_{i}=I\) arising in an interpolation problem. Linear Multilinear Algebra 61, 1192-1205 (2013)
Sakhnovich, LA: Interpolation Theory and Its Applications. Mathematics and Its Applications. Kluwer Academic, Dordrecht (1997)
Cheng, M, Xu, S: Perturbation analysis of the Hermitian positive definite solution of the matrix equation \(X-A^{*}X^{-2}A=I\). Linear Algebra Appl. 394, 39-51 (2005)
Duan, X, Li, C, Liao, A: Solutions and perturbation analysis for the nonlinear matrix equation \(X+\sum_{i=1}^{m}A_{i}^{*}X^{-1}A_{i}=I\). Appl. Math. Comput. 218, 4458-4466 (2011)
El-Sayed, SM, Al-Dubiban, AM: On positive definite solutions of the nonlinear matrix equation \(X+A^{*}X^{-n}A=I\). Appl. Math. Comput. 151, 533-541 (2004)
Hasanov, VI: Notes on two perturbation estimates of the extreme solutions to the equations \(X\pm A^{*}X^{-1}A=Q\). Appl. Math. Comput. 216, 1355-1362 (2010)
Ivanov, IG: On positive definite solutions of the family of matrix equations \(X+A^{*}X^{-n}A=Q\). J. Comput. Appl. Math. 193, 277-301 (2006)
Lin, WW, Xu, SF: Convergence analysis of structure-preserving doubling algorithms for Riccati-type matrix equations. SIAM J. Matrix Anal. Appl. 28, 26-39 (2006)
Chiang, CY, Chu, EKW, Guo, CH, Huang, TM, Lin, WW, Xu, SF: Convergence analysis of the doubling algorithm for several nonlinear matrix equations in the critical case. SIAM J. Matrix Anal. Appl. 31, 227-247 (2009)
Ivanov, IG, Hasanov, VI, Minchev, BV: On matrix equations \(X\pm A^{*}X^{-2}A=I\). Linear Algebra Appl. 326, 27-44 (2001)
El-Sayed, SM: Two iteration processes for computing positive definite solutions of the equation \(X-A^{*}X^{-n}A=Q\). Comput. Math. Appl. 41, 579-588 (2001)
El-Sayed, SM, Ran, ACM: On an iterative method for solving a class of nonlinear matrix equations. SIAM J. Matrix Anal. Appl. 23, 632-645 (2001)
Hasanov, VI, Ivanov, IG: Solutions and perturbation estimates for the matrix equations \(X\pm A^{*}X^{-n}A =Q\). Appl. Math. Comput. 156, 513-525 (2004)
Hasanov, VI: Positive definite solutions of the matrix equations \(X\pm A^{*}X^{-q}A=Q\). Linear Algebra Appl. 404, 166-182 (2005)
Duan, XF, Wang, QW, Li, CM: 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\). J. Appl. Math. Inform. 30, 655-663 (2012)
Liu, XG, Gao, H: On the positive definite solutions of the matrix equations \(X^{s}\pm A^{T}X^{-t}A=I_{n}\). Linear Algebra Appl. 368, 83-97 (2003)
Long, JH, Hu, XY, Zhang, L: Improved Newton’s method with exact linear searches to solve quadratic matrix equation. J. Comput. Appl. Math. 222, 645-654 (2008)
Yao, G, Liao, A, Duan, X: Positive definite solution of the matrix equation \(X=Q +A^{H}(I\otimes X-C)^{-\delta}A^{*}\). Electron. J. Linear Algebra 21, 76-84 (2010)
Ivanov, IG, El-Sayed, SM: Properties of positive definite solutions of the equation \(X + A^{*}X^{-2}A = I\). Linear Algebra Appl. 279, 303-316 (1998)
Parodi, M: La localisation des valeurs caracterisiques des matrices et ses applications. Gauthier-Villars, Paris (1959)
Acknowledgements
The author is grateful to the editor and the reviewer for important comments and suggestions, which improved the quality of the paper.
Author information
Authors and Affiliations
Corresponding author
Additional information
Competing interests
The author declares to have no competing interests.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
About this article
Cite this article
Al-Dubiban, A.M. On nonlinear matrix equations \(X\pm\sum_{i=1}^{m}A_{i}^{*}X^{-n_{i}}A_{i}=I\) . J Inequal Appl 2015, 147 (2015). https://doi.org/10.1186/s13660-015-0667-7
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13660-015-0667-7