On the Hermitian positive definite solutions of nonlinear matrix equation Xs + AXtA = Q

https://doi.org/10.1016/j.amc.2010.05.023Get rights and content

Abstract

Nonlinear matrix equation Xs + AXtA = Q, where A, Q are n × n complex matrices with Q Hermitian positive definite, has widely applied background. In this paper, we consider the Hermitian positive definite solutions of this matrix equation with two cases: s  1, 0 < t  1 and 0 < s  1, t  1. We derive necessary conditions and sufficient conditions for the existence of Hermitian positive definite solutions for the matrix equation and obtain some properties of the solutions. We also propose iterative methods for obtaining the extremal Hermitian positive definite solution of the matrix equation. Finally, we give some numerical examples to show the efficiency of the proposed iterative methods.

Introduction

We consider the nonlinear matrix equationXs+AX-tA=Q,where A is a nonsingular complex matrix and Q is a Hermitian positive definite matrix with two cases: s  1, 0 < t  1 and 0 < s  1, t  1.

Nonlinear matrix equations with the form of (1.1) have many applications in: dynamic programming; control theory; stochastic filtering; statistics etc., see [1], [4], [18]. The case that s = t = 1 has been extensively studied by several authors [1], [2], [3], [4], [5], [7], [11], [18], [19], and some existence conditions and properties of its positive definite solutions are obtained. Moreover, effective iterative methods for computing the Hermitian positive definite solutions are proposed. Ivanov and his coauthors [12], [13] consider the case that s = 1, t = 2 and Q = I. Liu and Gao [14] did some research on the matrix equation Xs + AXtA = I with s, t  N, where N denotes the set of natural numbers. Hasanov and Ivanov [8] discuss the Hermitian positive definite solutions of Eq. (1.1) for s = 1, t  N. Yang [17] considers the case s, t  N. In Hasanov [9], [10], the authors investigate the matrix equation when s = 1, t  (0, 1]. Then Peng et al. [16] propose iterative methods for the extremal positive definite solutions of Eq. (1.1) for s = 1 with two cases: 0 < t  1 and t  1.

In this paper, we discuss the Hermitian positive definite solutions of Eq. (1.1) with more general two cases: s  1, 0 < t  1 and 0 < s  1, t  1. The paper is organized as follows. In Section 2, we derive necessary conditions and sufficient conditions for the existence of the Hermitian positive definite solutions for Eq. (1.1) and present some properties of the solutions. Then in Section 3, we propose iterative methods for obtaining the maximal positive definite solution of Eq. (1.1) with the case s  1, 0 < t  1 and the minimal Hermitian positive definite solution of Eq. (1.1) with the case 0 < s  1, t  1. We give some numerical examples in Section 4 to show the efficiency of the proposed iterative methods. Conclusions will be put in Section 5.

The following notations are used throughout this paper. Let σ1(·), σn(·) denote the biggest and the smallest singular value of a n × n complex matrix and λmax(H), λmin(H) the maximal and the minimal eigenvalue of a Hermitian matrix H, respectively. The notation B  0 (B > 0) means that B is a positive semi-definite (definite) matrix. For two Hermitian matrices B and C, the notation B  C (B > C) indicates that B  C is a positive semi-definite (definite) matrix, and X  [B, C] implies that B  X  C. We call a Hermitian positive definite solution XL of the matrix equation the maximal solution if XL  X holds for any Hermitian positive definite solution X of the matrix equation. The minimal solution XM can be defined similarly. Moreover, A∗ stands for the conjugate transpose of the matrix A and I denotes the identity matrix.

Section snippets

Conditions for the existence and properties of the solutions

In this section, we shall derive necessary conditions and sufficient conditions for the existence of the Hermitian positive definite solutions for Eq. (1.1) and discuss some properties of the solutions.

Lemma 2.1

[15]

If B > C > 0 (or B  C > 0), then Bα > Cα > 0 (or Bα  Cα > 0) for all α  (0, 1] and 0 < Bα < Cα (or 0 < Bα  Cα) for all α  [1, 0).

Lemma 2.2

[6]

Let B and C be positive operators on a Hilbert space such that 0 < m1I  B  M1I, 0 < m2I  C  M2I and 0 < B  C. Then

  • (1)

    BpM1m1p-1Cp,

  • (2)

    BpM2m2p-1Cp

hold for any p  1.

Theorem 2.1

Let X be a Hermitian positive definite solution

Iterative methods for the extremal positive solution

In this section, we consider the iterative methods for obtaining the extremal Hermitian positive definite solution of Eq. (1.1).

To compute the maximal Hermitian positive definite solution XL of Eq. (1.1) for s  1 and 0 < t  1, we propose the following algorithm which avoids calculating matrix inversion in the process of iteration.

Algorithm 3.1

  • step1. Input initial matrices:

    • X0=γQ1s,

    • Y0=1γQ-1s,(γ[αˆ,1]), where αˆ is defined in Theorem 2.2.

  • step2. For k = 0, 1, 2,…, compute:

    • Yk+1 = (I  Xk)Yk + I,

    • Xk+1=(Q-AYk+1tA)1s.

Theorem 3.1

If Eq. (

Numerical examples

In this section, we give some numerical examples to illustrate the efficiency of the proposed algorithms. All the tests are performed by MATLAB 7.4 with machine precision around 10−16. We stop the practical iteration when the residual Xks+AXk-tA-QF<1.0e-010, where ∥·∥F denotes the Frobenius matrix norm.

Example 4.1

LetA=200-100-1200-100030-1010020-110-10300-1001-2,Q=310-12-486013562-12168-3854-11610-48-382892216876605422202-16220135-116168-163745062107622050278.

In this case, A is nonsingular and Q is

Conclusions

In this paper, we investigate the Hermitian positive definite solutions of nonlinear matrix equation Xs + AXtA = Q with more general two cases: s  1, 0 < t  1 and 0 < s  1, t  1. We derive some necessary conditions and sufficient conditions for the matrix equation to have a Hermitian positive definite solution. We also propose iterative methods for obtaining its extremal positive definite solution. Numerical results show that the proposed algorithms are quite efficient. As for the Hermitian positive

Acknowledgement

Research supported by the Natural Science Foundation of Zhejiang Province (No. Y6090164), the Natural Science Foundation of Science and Technology Commission of Shanghai Municipality (No. 09ZR1408700) and the Natural Science Foundation of Huzhou (No. 2010YZ05).

References (19)

There are more references available in the full text version of this article.

Cited by (0)

View full text