JAM Journal of Applied Mathematics 1687-0042 1110-757X Hindawi Publishing Corporation 518269 10.1155/2013/518269 518269 Research Article Sensitivity Analysis of the Matrix Equation from Interpolation Problems http://orcid.org/0000-0001-9843-5501 Li Jing 1 Zhang Yuhai 2 Alves Carlos J. S. 1 School of Mathematics and Statistics Shandong University Weihai Weihai 264209 China sdu.edu.cn 2 School of Mathematics Shandong University Jinan 250100 China sdu.edu.cn 2013 18 11 2013 2013 22 07 2013 27 09 2013 2013 Copyright © 2013 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.

This paper studies the sensitivity analysis of a nonlinear matrix equation connected to interpolation problems. The backward error estimates of an approximate solution to the equation are derived. A residual bound of an approximate solution to the equation is obtained. A perturbation bound for the unique solution to the equation is evaluated. This perturbation bound is independent of the exact solution of this equation. The theoretical results are illustrated by numerical examples.

1. Introduction

In this paper we consider the Hermitian positive definite solution of the nonlinear matrix equation: (1)X-i=1mAi*X-1Ai=Q, where A1,A2,,Am are n×n complex matrices, m is a positive integer, and Q is a positive definite matrix. Here, Ai* denotes the conjugate transpose of the matrix Ai.

This type of nonlinear matrix equations arises in many practical problems. The equation X-A*X-1A=Q which is a representative of (1) for m=1 comes from ladder networks, dynamic programming, control theory, stochastic filtering, statistics, and so forth . When m>1, (1) comes from the nonlinear matrix equation: (2)X=Q+A*(X^-C)-1A, where Q is an n×n positive definite matrix, C is an mn×mn positive semidefinite matrix, A is an arbitrary mn×n matrix, and X^ is the m×m block diagonal matrix with, on each diagonal entry, the n×n matrix X. In , (2) is recognized as playing an important role in modelling certain optimal interpolation problems. Let C=0 and A=(A1T,A2T,,AmT)T, where Ai,i=1,2,,m, are n×n matrices. Then X=Q+A*(X^-C)-1A can be rewritten as X-i=1mAi*X-1Ai=Q. When we solve the nonlinear matrix equation X-i=1mAi*X-1Ai=Q, we often do not know Ai and Q exactly but have only approximations A~i and Q~ available. Then we only solve the equation X~-i=1mA~i*X~-1A~i=Q~ exactly which gives a different solution X~. We would like to know how the errors of A~i and Q~ influence the error in X~. Motivated by this, we consider in this paper the sensitivity analysis of (1).

For the equation X-A*X-1A=Q and related equations X-i=1mAi*XδiAi=Q  (0<|δi|<1)  and  X-i=1mAi*XrAi=Q (-1r<0 or 0<r<1) there were many contributions in the literature to the solvability and numerical solutions . However, these papers have not examined the sensitivity analysis about the above equations. Hasanov et al. [12, 15] obtained two perturbation estimates of the solutions to the equations X±A*X-1A=Q. Li and Zhang  derived two perturbation bounds of the unique solution to the equation X-A*X-1A=Q. They also obtained the explicit expression of the condition number for the unique positive definite solution. The perturbation analysis about the related equations X+A*X-1A=P,   X-A*X-pA=Q, and Xs±ATX-tA=In were mentioned in papers . Yin and Fang  obtained an explicit expression of the condition number for the unique positive definite solution of (1). They also gave two perturbation bounds for the unique positive definite solution, whereas, to our best knowledge, there have been no backward error estimates and any computable residual bound for (1) in the known literature. In this paper, we obtain the backward error estimates and a residual bound of the approximate solution to (1) as well as evaluate a new relative perturbation bound for (1). This bound does not need any knowledge of the exact solution of (1), which is important in many practical calculations.

As a continuation of the previous results, this paper gives some preliminary knowledge that will be needed to develop this work in Section 2. In Section 3, the backward error estimates of an approximate solution for the unique solution to (1) are discussed. In Section 4, we derive a residual bound of an approximate solution for the unique solution to (1). In Section 5, we give a new perturbation bound for the unique solution to (1), which is independent of the exact solution of (1). Finally, several numerical examples are presented in Section 6.

We denote by 𝒞n×n the set of n×n complex matrices, by n×n the set of n×n Hermitian matrices, by I the identity matrix, by i the imaginary unit, by · the spectral norm, by ·F the Frobenius norm, and by λmax(M) and λmin(M) the maximal and minimal eigenvalues of M, respectively. For A=(a1,,an)=(aij)𝒞n×n and a matrix B, AB=(aijB) is a Kronecker product, and vecA is a vector defined by vecA=(a1T,,anT)T. For X,Yn×n, we write XY  (resp., X>Y) if X-Y is Hermitian positive semidefinite (resp., definite).

2. Preliminaries Lemma 1 (see [<xref ref-type="bibr" rid="B8">8</xref>, Theorem 2.1]).

The matrix equation (3)X-i=1mAi*XrAi=Q,-1r<0 always has a unique positive definite solution.

Lemma 2 (see [<xref ref-type="bibr" rid="B8">8</xref>, Theorem 2.3]).

Let X be the unique Hermitian positive definite solution of (3); then X[βI,αI], where the pair (β,α) is a solution of the system (4)β=λmin(Q)+i=1mλmin(Ai*Ai)αr,    α=λmax(Q)+i=1mλmax(Ai*Ai)βr.

3. Backward Error

In this section, applying the technique developed in , we obtain some estimates for the backward error of the approximate solution of (1).

Let X~n×n be an approximation to the unique solution X to (1), and let ΔAi𝒞n×n(i=1,2,,m) and ΔQn×n be the corresponding perturbations of the coefficient matrices Ai(i=1,2,,m) and Q in (1). A backward error of the approximate solution X~ can be defined by (5)η(X~)=min{(ΔA1α1,ΔA2α2,,ΔAmαm,ΔQρ)F:X~-i=1m(Ai+ΔAi)*X~-1(Ai+ΔAi)=Q+ΔQ(ΔA1α1,ΔA2α2,,ΔAmαm,ΔQρ)F}, where α1,α2,,αm and ρ are positive parameters. Taking αi=AiF,i=1,2,,m, and ρ=QF in (5) gives the relative backward error ηrel(X~), and taking αi=1,i=1,2,,m, and ρ=1 in (5) gives the absolute backward error ηabs(X~).

Let (6)R=Q-X~+i=1mAi*X~-1Ai. Note that (7)Q=X~-i=1m(Ai+ΔAi)*X~-1(Ai+ΔAi)-ΔQ. It follows from (6) that (8)-i=1m(ΔAi*X~-1Ai+Ai*X~-1ΔAi)-ΔQ=R+i=1mΔAi*X~-1ΔAi. Let (9)(I(X~-1Ai)*)=Ui1+iΩi1,  ((X~-1Ai)TI)Π=Ui2+iΩi2,vecΔAi=xi+iyi,vecΔQ=q1+iq2,  vecR=r1+ir2,vec(ΔAi*X~-1ΔAi)=ai+ibi,i=1,2,,m,g=(x1Tα1,y1Tα1,,xmTαm,ymTαm,q1Tρ,q2Tρ)T,Ui=(Ui1+Ui2Ωi2-Ωi1Ωi1+Ωi2Ui1-Ui2),T=[-α1U1,-α2U2,,-αmUm,-ρI2n2], where Π is the vec-permutation. Then (8) can be written as (10)Tg=(r1r2)+i=1m(aibi). It follows from ρ>0 that 2n2×2(m+1)n2 matrix T is full row rank. Hence, TT=I2n2, which implies that every solution to the equation (11)g=T(r1r2)+T(i=1m(aibi)) must be a solution to (10). Consequently, for any solution g to (11) we have (12)η(X~)g. Then we can state the estimates of the backward error as follows.

Theorem 3.

Let A1,A2,,Am,Q,X~𝒞n×n be given matrices and let η(X~) be the backward error defined by (5). If (13)r<s4t(i=1mαi2), then one has that (14)U(r)η(X~)B(r), where (15)r=T(r1r2),s=T-1,t=X~-1,(16)B(r)=2rss+s2-4rst(i=1mαi2),U(r)=2rs2-4rst(i=1mαi2)s+s2-4rst(i=1mαi2).

Proof.

Let (17)L(g)=T(r1r2)+T(i=1m(aibi)). Obviously, L:𝒞2(m+1)n2×1𝒞2(m+1)n2×1 is continuous. Condition (13) ensures that the quadratic equation (18)x=r+ts(i=1mαi2)x2 in x has two positive real roots. The smaller one is (19)B(r)=2rss+s2-4rst(i=1mαi2). Define Ω={g𝒞2(m+1)n2×1:gB(r)}. Then for any gΩ, we have (20)L(g)r+1si=1m(aibi)=r+1si=1mΔAi*X~-1ΔAiFr+tsi=1mΔAiF2r+ts(i=1mαi2)(ΔA1α1,ΔA2α2,,ΔAmαm)F2r+ts(i=1mαi2)g2r+ts(i=1mαi2)B2(r)=B(r). The last equality is due to the fact that B(r) is a solution to the quadratic equation (18). Thus we have proved that L(Ω)Ω. By the Schauder fixed-point theorem, there exists a g*Ω such that L(g*)=g*, which means that g* is a solution to (11), and hence it follows from (12) that (21)η(X~)g*B(r). Next we derive a lower bound for η(X~). Suppose that (ΔA1min/α1,,ΔAmmin/αm,ΔQmin/ρ) satisfies (22)η(X~)=(ΔA1minα1,,ΔAmminαm,ΔQminρ)F. Then we have (23)Tgmin=(r1r2)+i=1m(aibi), where (24)vec(ΔAimin*X~-1ΔAimin)=ai+ibi,vec(ΔAimin)=xi+iyi,vec(ΔQmin)=q1+iq2,gmin=(x1Tα1,y1Tα1,,xmTαm,ymTαm,q1Tρ,q2Tρ)T. Let a singular value decomposition of T be T=W(E,0)Z*, where W and Z are unitary matrices and E=diag(e1,e2,,e2n2) with e1e2n2>0. Substituting this decomposition into (23), and letting (25)Z*gmin=(v),v𝒞2n2×1, we get (26)v=E-1W*(r1r2)+E-1W*i=1m(aibi). It follows from (22) that (27)η(X~)=gmin=(v)vE-1W*(r1r2)-E-1W*i=1m(aibi)T(r1r2)-T·i=1m(aibi)r-1si=1mΔAimin*X~-1ΔAiminFr-tsi=1mΔAiminF2r-ts(i=1mαi2)(ΔA1minα1,,ΔAmminαm)F2r-ts(i=1mαi2)B2(r). Here we have used the fact that (28)(ΔA1minα1,,ΔAmminαm)F(ΔA1minα1,,ΔAmminαm,ΔQminρ)F=η(X~)B(r). Let (29)U(r)=r-ts(i=1mαi2)B2(r). Since B(r) is a solution to (18), we have that (30)B(r)=r+ts(i=1mαi2)B2(r), which implies that (31)U(r)=r-ts(i=1mαi2)B2(r)=2r-B(r)=2rs2-4rst(i=1mαi2)s+s2-4rst(i=1mαi2)>0. Then η(X~)U(r).

4. Residual Bound

Residual bound reveals the stability of a numerical method. In this section, in order to derive the residual bound of an approximate solution for the unique solution to (1), we first introduce the following lemma.

Lemma 4.

For every positive definite matrix Xn×n, if X+ΔX(1/ν)I>0, then (32)i=1mAi*((X+ΔX)-1-X-1)Ai(ΔX+νΔX2)i=1mX-1Ai2.

Proof.

According to (33)(X+ΔX)-1-X-1=-X-1ΔX(X+ΔX)-1=-X-1ΔXX-1+X-1ΔXX-1ΔX(X+ΔX)-1, it follows that (34)i=1mAi*((X+ΔX)-1-X-1)Aii=1m(Ai*X-1ΔXX-1Ai+Ai*X-1ΔXX-1ΔX(X+ΔX)-1Ai)(ΔX+νΔX2)i=1mX-1Ai2.

Theorem 5.

Let X~>0 be an approximation to the solution X of (1). If the residual R(X~)Q+i=1mAi*X~-1Ai-X~ satisfies (35)R(X~)<(1-Σ)21+Σ+2Σ  λmin(X~),where  Σi=1mX~-1Ai2<1, then (36)X~-XθR(X~), where (37)θ=(2λmin(X~))×(((1-Σ)λmin(X~)+R(X~))2)1/2(1-Σ)λmin(X~)+R(X~)+(((1-Σ)λmin(X~)+R(X~))2-4λmin(X~)R(X~)((1-Σ)λmin(X~)+R(X~))2)1/2)-1.

Proof.

Let (38)Ψ={ΔXn×n:ΔXθR(X~)}. Obviously, Ψ is a nonempty bounded convex closed set. Let (39)g(ΔX)=i=1mAi*[(X~+ΔX)-1-X~-1]Ai+R(X~). Evidently g:Ψn×n is continuous.

Note that condition (35) ensures that the quadratical equation (40)x2-(λmin(X~)(1-Σ)+R(X~))x+λmin(X~)R(X~)=0 has two positive real roots, and the smaller one is given by (41)μ*=(2λmin(X~)R(X~))×(((1-Σ)λmin(X~)+R(X~))2)1/2(1-Σ)λmin(X~)+R(X~)+(((1-Σ)λmin(X~)+R(X~))2-4λmin(X~)R(X~)((1-Σ)λmin(X~)+R(X~))2)1/2)-1. Next, we will prove that g(Ψ)Ψ.

For every ΔXΨ, we have (42)ΔX-θR(X~)I. Hence (43)X~+ΔXX~-θR(X~)I(λmin(X~)-θR(X~))I. By (36), one sees that (44)θR(X~)2λmin(X~)R(X~)(1-Σ)λmin(X~)+R(X~)=λmin(X~)×(1+R(X~)-(1-Σ)λmin(X~)(1-Σ)λmin(X~)+R(X~)). According to (35), we obtain (45)R(X~)-(1-Σ)λmin(X~)((1-Σ)21+Σ+2Σ-(1-Σ))λmin(X~)-2(1-Σ)λmin(X~)1+Σ<0, which implies that (46)θR(X~)λmin(X~),(λmin(X~)-θR(X~))I>0.

According to Lemma 4, we obtain (47)g(ΔX)(ΔX+ΔX2λmin(X~)-θR(X~))×i=1mX-1Ai2+R(X~)(θR(X~)+(θR(X~))2λmin(X~)-θR(X~))Σ+R(X~)=θR(X~), for ΔXΨ. That is, g(Ψ)Ψ. By Brouwer fixed point theorem, there exists a ΔXΨ such that g(ΔX)=ΔX. Hence X~+ΔX is a solution of (1). Moreover, by Lemma 1, we know that the solution X of (1) is unique. Then (48)X~-X=ΔXθR(X~).

5. Perturbation Bounds

In this section we develop a relative perturbation bound for the unique solution of (1), which does not need any knowledge of the actual solution X of (1) and is easy to calculate.

Here we consider the perturbed equation: (49)X~-i=1mAi~*X~-1Ai~=Q~, where Ai~ and Q~ are small perturbations of Ai and Q in (1), respectively. It follows from Lemma 1 that the solutions of (1) and (49) exist. Then we assume that X and X~ are the solutions of (1) and (49), respectively. Let ΔX=X~-X, ΔQ=Q~-Q, and ΔAi=Ai~-Ai.

Theorem 6.

If (50)i=1mAi2<β2, then (51)ΔXXi=1m(2Ai+ΔAi)ΔAi+ΔQβ2-i=1mAi2ξ1.

Proof.

By Lemma 1, we know that X and X~ are the unique solutions to (1) and (49), respectively. Subtracting (1) from (49) we have (52)ΔX+i=1mAi*X~-1ΔXX-1Ai=i=1m(Ai*X~-1ΔAi+ΔAi*X~-1Ai+ΔAi*X~-1ΔAi)+ΔQ. Then (53)ΔX+i=1mAi*X~-1ΔXX-1AiΔX-i=1mAi*X~-1ΔXX-1Ai(1-i=1mAi2X~-1X-1)ΔX. By Lemma 2, it follows that X-11/β and X~-11/β. Then (54)ΔX+i=1mAi*X~-1ΔXX-1Ai1β2(β2-i=1mAi2)ΔX. Condition (50) ensures (β2-i=1mAi2)ΔX>0.

Combining (52) and (54), we obtain (55)ΔXXi=1m(2Ai+ΔAi)ΔAi+ΔQβ2-i=1mAi2.

Remark 7.

Yin and Fang  obtained two perturbation bounds, which were dependent on the exact solution of (1), whereas, in this paper, the relative perturbation bound in Theorem 6 does not need any knowledge of the actual solution X of (1), which is important in many practical calculations.

Remark 8.

With (56)ξ1=i=1m(2Ai+ΔAi)ΔAi+ΔQβ2-i=1mAi2, we get ξ10 as ΔQ0 and ΔAi0  (i=1,2,,m). Therefore (1) is well-posed.

6. Numerical Examples

To illustrate the theoretical results of the previous sections, in this section several simple examples are given, which were carried out using MATLAB 7.1. For the stopping criterion we take εk+1(X)=Xk-i=1mAi*Xk-1Ai-Q<1.0e-10.

Example 1.

In this example, we consider the backward error of an approximate solution for the unique solution X to (1) in Theorem 3. We consider (57)X-A1*X-1A1-A2*X-1A2=Q, with the coefficient matrices (58)A1=15(101-111-1-11),A2=2345A1,Q=X-A1*X-1A1-A2*X-1A2, where X=diag(1,2,3).

Let (59)X~=X+(0.5-0.10.2-0.10.30.60.20.6-0.4)×10-j be an approximate solution to (1). Take α1=A1F, α2=A2F and ρ=QF in Theorem 3. Some results on lower and upper bounds for the backward error η(X~) are displayed in Table 1.

The results listed in Table 1 show that the backward error of X~ decreases as the error X~-XF decreases.

Backward error for Example 1 with different values of  j.

j X ~ - X F r U(r) B(r)
1 0.2298 0.0633 0.0630 0.0637
3 2.3 × 1 0 - 3 6.3391 × 1 0 - 4 6.3387 × 1 0 - 4 6.3394 × 1 0 - 4
5 2.2978 × 1 0 - 5 6.3391 × 1 0 - 6 6.3391 × 1 0 - 6 6.3391 × 1 0 - 6
7 2.2978 × 1 0 - 7 6.3391 × 1 0 - 8 6.3391 × 1 0 - 8 6.3391 × 1 0 - 8
9 2.2978 × 1 0 - 9 6.3391 × 1 0 - 10 6.3391 × 1 0 - 10 6.3391 × 1 0 - 10
Example 2.

This example considers the residual bound of an approximate solution for the unique solution X to (1) in Theorem 5. We consider (60)X-A1*X-1A1-A2*X-1A2=Q, with (61)A1=(1/3)+2×10-2AA,A2=(1/6)+3×10-2AA,Q=A=(2100012100012100012100012).

Choose X~0=A. Let the approximate solution X~k of X be given with the iterative method Xk=Q+i=1mAi*Xk-1-1Ai,X0>0,k=1,2,, where k is the iteration number.

The residual R(X~k)Q+A1*X~k-1A1+A2*X~k-1A2-X~k satisfies the conditions in Theorem 5. By Theorem 5, we can compute the residual bounds for X~k(62)X~k-XθR(X~k), where (63)θ=(2λmin(X~k))×(((1-Σ)λmin(X~)+R(X~))2)1/2(1-Σ)λmin(X~k)+R(X~k)+(((1-Σ)λmin(X~k)+R(X~k))2-4λmin(X~k)R(X~k)((1-Σ)λmin(X~k)+R(X~k))2)1/2)-1. Some results are listed in Table 2.

The results listed in Table 2 show that the residual bound given by Theorem 5 is fairly sharp.

Residual bounds for Example 2 with different values of  k.

k 1 2 3 4
X ~ k - X 5.0268 × 1 0 - 4 5.7662 × 1 0 - 6 6.6162 × 1 0 - 8 7.5024 × 1 0 - 10
θ R ( X ~ k ) 5.1435 × 1 0 - 4 5.9000 × 1 0 - 6 6.7689 × 1 0 - 8 7.7656 × 1 0 - 10
Example 3.

In this example, we consider the corresponding perturbation bound for the solution X in Theorem 6.

We consider the matrix equation (64)X-A1*X-1A1-A2*X-1A2=I, with (65)A1=(1/3)+2×10-2AA,A2=(1/6)+3×10-2AA,  (66)A=(2100012100012100012100012). Suppose that the coefficient matrices A1 and A2 are perturbed to Ai~=Ai+ΔAi,  i=1,2, where (67)ΔA1=10-jCT+C(CT+C),  ΔA2=3×10-j-1CT+C(CT+C) and C is a random matrix generated by MATLAB function randn.

By Theorem 6, we can compute the relative perturbation bound ξ1. The results averaged as the geometric mean of 20 randomly perturbed runs. Some results are listed in Table 3.

The results listed in Table 3 show that the perturbation bound ξ1 given by Theorem 6 is fairly sharp.

Perturbation bounds for Example 3 with different values of  j.

j 4 5 6 7
X ~ - X / X 2.7482 × 1 0 - 5 2.4983 × 1 0 - 6 2.5705 × 1 0 - 7 2.9406 × 1 0 - 8
ξ 1 1.0845 × 1 0 - 4 9.2695 × 1 0 - 6 9.6710 × 1 0 - 7 1.0595 × 1 0 - 7
7. Concluding Remarks

In this paper, we consider the sensitivity analysis of the nonlinear matrix equation X-i=1mAi*X-1Ai=Q. Compared with existing literature, the contributions of this paper are as follows.

A backward error and a computable residual bound of an approximate solution for the unique solution to (1) are derived, which do not appear in other known literature works.

Some results in this paper can cover the work of Li and Zhang  for the matrix equation X-A*X-1A=Q as a special case.

This paper develops a new relative perturbation bound for the solution to (1), which does not need any knowledge of the actual solution X of (1) and could be computed easily.

Acknowledgments

The authors would like to express their gratitude to the referees for their fruitful comments. 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. The authors declare that there is no conflict of interests regarding the publication of this paper.

Anderson W. N. Kleindorfer G. B. Kleindorfer M. B. Woodroofe M. B. Consistent estimates of the parameters of a linear systems The Annals of Mathematical Statistics 1969 40 6 2064 2075 10.1214/aoms/1177697286 Anderson W. N. Morley T. D. Trapp G. E. Christopher Byrnes I. Martin F. C. Richard Saeks E. The cascade limit, the shorted operator and quadratic optimal control Linear Circuits, Systems and Signal Processsing: Theory and Application 1988 New York, NY, USA North-Holland 3 7 Bucy R. S. A priori bounds for the Riccati equation Proceedings of the 6th Berkeley Symposium on Mathematical Statistics and Probability, Volume 3: Probability Theory 1972 Berkeley, Calif, USA University of California Press 645 656 Ouellette D. V. Schur complements and statistics Linear Algebra and Its Applications C 1981 36 187 295 2-s2.0-0041568115 10.1016/0024-3795(81)90232-9 ZBL0455.15012 Pusz W. Woronowitz S. L. Functional caculus for sequlinear forms and purification map Reports on Mathematical Physics 1975 8 2 159 170 10.1016/0034-4877(75)90061-0 Zabcyk J. Remarks on the control of discrete time distributed parameter systems SIAM Journal on Control and Optimization 1974 12 4 721 735 2-s2.0-0016127198 10.1137/0312056 Ran A. C. M. Reurings M. C. B. A nonlinear matrix equation connected to interpolation theory Linear Algebra and Its Applications 2004 379 1–3 289 302 2-s2.0-0742285659 10.1016/S0024-3795(03)00541-X ZBL1039.15007 Duan X. F. Liao A. P. On Hermitian positive definite solution of the matrix equation X-i=1mAi*XrAi=Q Journal of Computational and Applied Mathematics 2009 229 1 27 36 2-s2.0-65049089159 10.1016/j.cam.2008.10.018 Ferrante A. Levy B. C. Hermitian solutions of the equation X=Q+NX-1N* Linear Algebra and Its Applications 1996 247 359 373 2-s2.0-0030520753 10.1016/0024-3795(95)00121-2 Guo C. Lancaster P. Iterative solution of two matrix equations Mathematics of Computation 1999 68 228 1589 1603 2-s2.0-0033443128 10.1090/S0025-5718-99-01122-9 ZBL0940.65036 Hasanov V. I. Positive definite solutions of the matrix equations X±A*X-qA=Q Linear Algebra and Its Applications 2005 404 1–3 166 182 2-s2.0-20344385000 10.1016/j.laa.2005.02.024 Hasanov V. I. Ivanov I. G. Uhlig F. Improved perturbation estimates for the matrix equations X±A*X-1A=Q Linear Algebra and Its Applications 2004 379 1–3 113 135 2-s2.0-0742320632 10.1016/S0024-3795(03)00424-5 Li J. Zhang Y. H. The Hermitian positive definite solutions and perturbation analysis of the matrix equation X-A*X-1A=Q Bulletin of the Institute of Mathematics Academia Sinica 2008 30 129 142 Lim Y. Solving the nonlinear matrix equation X=Q+i=1mMiXδiMi* via a contraction principle Linear Algebra and Its Applications 2009 430 4 1380 1383 2-s2.0-58049195551 10.1016/j.laa.2008.10.034 Hasanov V. I. Ivanov I. G. On two perturbation estimates of the extreme solutions to the equations X±A*X-1A=Q Linear Algebra and Its Applications 2006 413 1 81 92 2-s2.0-29344460370 10.1016/j.laa.2005.08.013 Li J. Zhang Y. Perturbation analysis of the matrix equation X-A*X-qA=Q Linear Algebra and Its Applications 2009 431 9 1489 1501 2-s2.0-68349157686 10.1016/j.laa.2009.05.013 ZBL1173.15006 Liu X. G. Gao H. On the positive definite solutions of the matrix equations Xs±ATX-tA=In Linear Algebra and Its Applications 2003 368 83 87 2-s2.0-0038730582 10.1016/S0024-3795(02)00661-4 Sun J. G. Xu S. F. Perturbation analysis of the maximal solution of the matrix equation X+A*X-1A=P. II Linear Algebra and Its Applications 2003 362 211 228 2-s2.0-0037444951 10.1016/S0024-3795(02)00490-1 Xu S. F. Perturbation analysis of the maximal solution of the matrix equation X+A*X-1A=P Linear Algebra and Its Applications 2001 336 1–3 61 70 2-s2.0-0035888250 10.1016/S0024-3795(01)00300-7 Yin X. Y. Fang L. Perturbation analysis for the positive definite solution of the nonlinear matrix equation X-i=1mAi*X-1Ai=Q Journal of Applied Mathematics and Computing 2013 43 1 199 211 10.1007/s12190-013-0659-z