ABSTRACT

Consider the nonlinear matrix equation $X-A^{\ast }{X}^{p}A-B^{\ast }{X}^{-q}B=I\left(0<p,q<1\right)$ . By using the fixed point theorem for mixed monotone operator in a normal cone, we prove that the equation with $0<p,q<1$ always has the unique positive definite solution. Two different iterative methods are given, including the basic fixed point iterative method and the multi-step stationary iterative method. Numerical examples show that the iterative methods are feasible and effective.

Keywords:

Nonlinear Matrix Equation, Positive Definite Solution, Iterative Method, Normal Cone

1. Introduction

We consider the matrix equation

$X-A^{*}{X}^{p}A-B^{*}{X}^{-q}B=I,0<p,q<1$ (1.1)

where $I$ is an $n\times n$ identity matrix, $A\mathrm{,}B$ are arbitrary $mn\times n$ matrices. $A^{\mathrm{<mo>\; *\; </mo>}}$ , $B^{*}$ are used to denote the conjugate transpose of the matrix A and B, separately.

Equation (1.1) often arises in dynamic programming, control theory, stochastic filtering, statistics, and so on, see [1] . Equation (1.1) has been investigated in some special cases [1] - [10] . Reurings [2] showed that the matrix equation $X-A^{*}{X}^{p}A=I\left(0<p<1\right)$ and $X-A^{*}{X}^{-q}A=I\left(0<q<1\right)$ always has a unique positive definite solution. Recently, Gao [3] provides a new proof for the uniqueness of the positive definite solution of the matrix equation $X-A^{*}{X}^{q}A=I\left(0<q<1\right)$ via the fixed point theorem. The method is shown much easier than the way of [1] . Based on this, we study the combination of these two kinds of equations (i.e. Equation (1.1)).

In this paper, we prove that the Equation (1.1) always has the unique positive definite solution by using the fixed point theorem for mixed monotone operator in the normal cone $P\left(n\right)$ . Furthermore, two different iterative methods are given, including the basic fixed point iteration method and multi-step stationary iterative method. In addition, numerical examples show that the iterative methods are feasible and effective.

In this paper, we use $H\left(n\right)$ to denote $n\times n$ Hermitian matrices, $\mathcal{P}\left(n\right)$ to denote $n\times n$ positive definite matrices, and $\overline{\mathcal{P}}\left(n\right)$ to denote $n\times n$ positive semidefinite matrices. For $X\mathrm{,}Y\in \mathcal{P}\left(n\right)$ , we write $X\ge Y$ $\left(X>Y\right)$ if $X-Y$ is positive semidefinite(definite).

2. Uniqueness of Positive Definite Solution

In this section, we will prove that the equation $X-A^{*}{X}^{p}A-B^{*}{X}^{-q}B=$ $I\left(0<p,q<1\right)$ always has the unique positive definite solution. In the process, we need to use the fixed point theorem for mixed monotone operator in a normal cone. So we first introduce the relevant theory through several definitions and lemmas. In this paper, we always use P to denote a solid cone of a real Banach space E. $P^0$ to denote the interior points set of P. A cone is said to be a solid cone if $P^0\ne \varphi$ .

Definition 2.1. [11] . A cone $P\subset E$ is said to be normal if there exists a constant $N>0$ such that $0\le x\le y$ implies $\Vert x\Vert \le N\Vert y\Vert$ . That is the norm $\Vert \cdot \Vert$ is semimonotone.

Definition 2.2. [11] . The operator $\Gamma \mathrm{<mo>\; :\; </mo>}D\times D\to E\mathrm{,}D\subset P$ is said to be a mixed monotone operator if $\Gamma \left(x\mathrm{,}y\right)$ is increasing in $x$ and decreasing in $y$ , that is,

$\Gamma \left(x_1,y_1\right)\ge \Gamma \left(x_2,y_2\right),\forall x_1,x_2,y_1,y_2\in D\text{with}{x}_1\le x_2\text{and}{y}_2\le y_1.$

An element $x^{\mathrm{<mo>\; *\; </mo>}}$ is called a fixed point of $\Gamma$ if $\Gamma \left(x^{*},x^{*}\right)=x^{*}$ .

Lemma 2.1. [12] . Let P be a normal and solid cone of a real Banach space E and $\Gamma \mathrm{<mo>\; :\; </mo>}{P}^0\times P^0\to P^0$ be a mixed monotone operator. Assume that for all $0<t<1$ , there exists $0<\beta <1$ such that

$\Gamma \left(tx\mathrm{,}\frac{1}{t}y\right)\ge t^{\beta }\Gamma \left(x\mathrm{,}y\right)$

holds for all $x\mathrm{,}y\in P^0$ . Then $\Gamma$ has exactly one fixed point $x^{\mathrm{<mo>\; *\; </mo>}}$ in $P^0$ . And for any $x_0\mathrm{,}{y}_0\in P^0$ , we have

$\underset{n\to \infty }{\lim }{x}_n=\underset{n\to \infty }{\lim }{y}_n=x^{*},$

where

$x_n=\Gamma \left(x_{n-1},y_{n-1}\right),y_n=\Gamma \left(y_{n-1},x_{n-1}\right),n=0,1,2,\cdots$

Lemma 2.2. [11] . A cone P is normal if and only if $x_n\le z_n\le y_n$ , $x_n\to x$ and $y_n\to x$ imply $z_n\to x$ .

Lemma 2.3. [13] . If $A\ge B>0$ and $0\le \gamma \le 1$ , then $A^{\gamma }\ge B^{\gamma }$ .

We define the spectral norm $\Vert \cdot \Vert$ in $H\left(n\right)$ , then $H\left(n\right)$ is a real Banach space. We all know that $\overline{P}\left(n\right)$ is a cone in $H\left(n\right)$ and the interior points set is $P\left(n\right)$ . Since the spectral norm is monotone, then we have that the set $\overline{P}\left(n\right)$ is normal cone from Definition 2.1. We define F by $F\left(X,Y\right)=I+A^{*}{X}^{p}A+B^{*}{Y}^{-q}B$ , then we have the following theorem.

Theorem 2.1. Eqution (1.1) always has a unique positive definite solution X. And for any $X_0\mathrm{,}{Y}_0\in P^0$ , we have

$\underset{n\to \infty }{\lim }{X}_n=\underset{n\to \infty }{\lim }{Y}_n=X^{*},$

where

$X_n=F\left(X_{n-1},Y_{n-1}\right),Y_n=F\left(Y_{n-1},X_{n-1}\right),n=0,1,2,\cdots$

Proof. Consider that the solution of Equation (1.1) is a fixed point of F. Now we will prove that the operator F satisfies the conditions of Lemma 2.1.

For $\forall X\mathrm{,}Y\mathrm{,}Z\mathrm{,}W\in P\left(n\right)$ with $X\ge Y\mathrm{,}W\ge Z$ , by Lemma 2.3, we have

1) $F\left(X,Y\right)=I+A^{*}{X}^{p}A+B^{*}{Y}^{-q}B\in P\left(n\right)$ , that is $F\mathrm{<mo>\; :\; </mo>}P\left(n\right)\times P\left(n\right)\to P\left(n\right)$ ;

2) $\begin{array}{c}F\left(X,Z\right)=I+A^{*}{X}^{p}A+B^{*}{Z}^{-q}B\\ \ge I+A^{*}{Y}^{p}A+B^{*}{W}^{-q}B\\ =F\left(Y,W\right),\end{array}$

which means that the operator F is mixed monotone;

3) Let $\beta =\max \left\{p,q\right\}$ , then we have $0<\beta <1$ . For $\forall t\in \left(\mathrm{0,1}\right)$ , we get

$\begin{array}{c}F\left(tX,\frac{1}{t}Y\right)=I+A^{*}{t}^p{X}^{p}A+B^{*}{\left(\frac{1}{t}\right)}^{-q}{Y}^{-q}B\\ =I+A^{*}{t}^p{X}^{p}A+B^{*}{t}^q{Y}^{-q}B\\ \ge t^{\beta }I+t^{\beta }{A}^{*}{X}^{p}A+t^{\beta }{B}^{*}{Y}^{-q}B\\ =t^{\beta }F\left(X,Y\right).\end{array}$

According to Lemma 2.1, we get that the operator $F$ has a unique fixed point $X$ in $P\left(n\right)$ , which is the unique positive definite solution of Equation (1.1).

Furthermore, for any $X_0\mathrm{,}{Y}_0\in P^0$ , we have

$\underset{n\to \infty }{\lim }{X}_n=\underset{n\to \infty }{\lim }{Y}_n=X^{*},$

where

$X_n=F\left(X_{n-1},Y_{n-1}\right),Y_n=F\left(Y_{n-1},X_{n-1}\right),n=0,1,2,\cdots$

,

3. Iterative Methods

In this section, we propose two different iterative methods for solving Equation (1.1), including the basic fixed point iterative method and multi-step stationary iterative method.

According to Theorem 2.1, sequences $X_n=F\left(X_{n-1},Y_{n-1}\right),Y_n=F\left(Y_{n-1},X_{n-1}\right),n=0,1,2,\cdots$ converge to the unique positive definite solution of Equation (1.1). Let $X_0=Y_0$ , we have $X_n=Y_n=F\left(X_{n-1},X_{n-1}\right)=I+A^{*}{X}_{n-1}^{p}A+B^{*}{X}_{n-1}^{-q}B$ . Consider the above iteration, which is the basic fixed point iterative method. The following are the algorithm.

Take $X_0=I$

$X_{n+1}=I+A^{*}{X}_n^{p}A+B^{*}{X}_n^{-q}B,n=0,1,2,\cdots$ (3.1)

We then consider a kind of multi-step stationary iterative method.

Theorem 3.1. For arbitrary initial matrix $X_1\mathrm{,}{X}_2\in P\left(n\right)$ , the matrix sequence

$X_{n+1}=I+A^{*}{X}_{n-1}^{p}A+B^{*}{X}_n^{-q}B,n=2,3,\cdots$ (3.2)

converges to the unique positive definite solution X of Equation (1.1).

Proof. Since X is the unique positive definite solution of Equation (1.1). For X, $X_1$ and $X_2$ , there exist a positive number $0<\alpha \le 1$ satisfy:

$\alpha X\le X_1,X_2\le {\alpha }^{-1}X.$ (3.3)

For any positive integer k, there exists a unique nonnegative integer K such that $k=2K+1$ or $k=2K+2$ . We will use mathematical induction to prove the following inequality

${\alpha }^{{\beta }^K}X\le X_k\le {\left({\alpha }^{-1}\right)}^{{\beta }^K}X\mathrm{,}k=\mathrm{1,2,}\cdots$ (3.4)

where $\beta =\max \left\{p,q\right\}$ . When $K=0$ , that is $k=1,2$ , the inequality (3.4) holds following (3.3). When $K=n-1$ , that is $k=2n-1,2n$ , Assume (3.4) is true. That is

${\alpha }^{{\beta }^{n-1}}X\le X_{2n-1},X_{2n}\le {\left({\alpha }^{-1}\right)}^{{\beta }^{n-1}}X.$ (3.5)

Now we need to prove (3.4) is true for $K=n$ , at this time $k=2n+1,2n+2$ . From (3.5) and $\beta =\max \left\{p,q\right\}$ , we get that

${\alpha }^{{\beta }^n}{X}^p\le {\alpha }^{{\beta }^{n-1}p}{X}^p\le X_{2n-1}^p\le {\left({\alpha }^{-1}\right)}^{{\beta }^{n-1}p}{X}^p\le {\left({\alpha }^{-1}\right)}^{{\beta }^n}{X}^p,$

$\begin{array}{c}{\alpha }^{{\beta }^n}{X}^{-q}\le {\alpha }^{{\beta }^{n-1}p}{X}^{-q}={\left({\alpha }^{-1}\right)}^{{\beta }^{n-1}\left(-q\right)}{X}^{-q}\le X_{2n}^{-q}\\ \le {\alpha }^{{\beta }^{n-1}\left(-q\right)}{X}^{\left(-q\right)}={\left({\alpha }^{-1}\right)}^{{\beta }^{n-1}q}{X}^{-q}\le {\left({\alpha }^{-1}\right)}^{{\beta }^n}{X}^{-q},\end{array}$

From $0<\alpha \le 1$ , it follows that $0<{\alpha }^{{\beta }^n}\le 1$ , ${\left({\alpha }^{-1}\right)}^{{\beta }^n}\ge 1$ . Then

$\begin{array}{c}{X}_{2n+1}=I+A^{*}{X}_{2n-1}^{p}A+B^{*}{X}_{2n}^{-q}B\\ \ge {\alpha }^{{\beta }^n}I+A^{*}{\alpha }^{{\beta }^n}{X}^{p}A+B^{*}{\alpha }^{{\beta }^n}{X}^{-q}B\\ ={\alpha }^{{\beta }^n}\left[I+A^{*}{X}^{p}A+B^{*}{X}^{-q}B\right]\\ ={\alpha }^{{\beta }^n}X\end{array}$

and

$\begin{array}{c}{X}_{2n+1}=I+A^{*}{X}_{2n-1}^{p}A+B^{*}{X}_{2n}^{-q}B\\ \le {\left({\alpha }^{-1}\right)}^{{\beta }^n}I+A^{*}{\left({\alpha }^{-1}\right)}^{{\beta }^n}{X}^{p}A+B^{*}{\left({\alpha }^{-1}\right)}^{{\beta }^n}{X}^{-q}B\\ ={\left({\alpha }^{-1}\right)}^{{\beta }^n}\left[I+A^{*}{X}^{p}A+B^{*}{X}^{-q}B\right]\\ ={\left({\alpha }^{-1}\right)}^{{\beta }^n}X.\end{array}$

Therefore

${\alpha }^{{\beta }^n}X\le X_{2n+1}\le {\left({\alpha }^{-1}\right)}^{{\beta }^n}X\mathrm{.}$ (3.6)

From $0<q\le 1$ , (3.6) implies

${\alpha }^{{\beta }^n}{X}^{-q}\le {\alpha }^{{\beta }^{n}q}{X}^{-q}\le X_{2n+1}^{-q}\le {\left({\alpha }^{-1}\right)}^{{\beta }^{n}q}{X}^{-q}\le {\left({\alpha }^{-1}\right)}^{{\beta }^n}{X}^{-q}.$

From (3.4), we have

${\alpha }^{{\beta }^n}{X}^p\le {\alpha }^{{\beta }^{n-1}p}{X}^p\le X_{2n}^p\le {\left({\alpha }^{-1}\right)}^{{\beta }^{n-1}p}{X}^p\le {\left({\alpha }^{-1}\right)}^{{\beta }^n}{X}^p\mathrm{.}$

Then

$\begin{array}{c}{X}_{2n+2}=I+A^{*}{X}_{2n}^{p}A+B^{*}{X}_{2n+1}^{-q}B\\ \ge {\alpha }^{{\beta }^n}I+A^{*}{\alpha }^{{\beta }^n}{X}^{p}A+B^{*}{\alpha }^{{\beta }^n}{X}^{-q}B\\ ={\alpha }^{{\beta }^n}\left[I+A^{*}{X}^{p}A+B^{*}{X}^{-q}B\right]\\ ={\alpha }^{{\beta }^n}X\end{array}$

and

$\begin{array}{c}{X}_{2n+2}=I+A^{*}{X}_{2n}^{p}A+B^{*}{X}_{2n+1}^{-q}B\\ \le {\left({\alpha }^{-1}\right)}^{{\beta }^n}I+A^{*}{\left({\alpha }^{-1}\right)}^{{\beta }^n}{X}^{p}A+B^{*}{\left({\alpha }^{-1}\right)}^{{\beta }^n}{X}^{-q}B\\ ={\left({\alpha }^{-1}\right)}^{{\beta }^n}\left[I+A^{*}{X}^{p}A+B^{*}{X}^{-q}B\right]\\ ={\left({\alpha }^{-1}\right)}^{{\beta }^n}X.\end{array}$

Therefore

${\alpha }^{{\beta }^n}X\le X_{2n+2}\le {\left({\alpha }^{-1}\right)}^{{\beta }^n}X\mathrm{.}$

Now we have proved that (3.4) is true for $K=n$ , at this time $k=2n+1,2n+2$ .

Hence we get that the inequality (3.4) holds for any positive integer k. Let $k\to \infty$ , then $K\to \infty$ , we have

${\alpha }^{{\beta }^K}\to \mathrm{1,}{\left({\alpha }^{-1}\right)}^{{\beta }^K}\to 1.$

Therefore,

$X_k\to X\mathrm{.}$

4. Numerical Examples

We now use numerical examples to illustrate our results. All computations were performed using MATLAB, version 7.01. we denote $\epsilon \left(X\right)={\Vert X-I-A^{\mathrm{<mo>\; *\; </mo>}}{X}^{p}A-B^{\mathrm{<mo>\; *\; </mo>}}{X}^{-q}B\Vert }_{\infty }$ , and use the stopping criterion $\epsilon \left(X\right)<1.0\times {10}^{-10}$ .

Example 4.1. Consider Equation (1.1) with $p=\frac{1}{3},q=\frac{1}{2}$ and

$A=\left(\begin{array}{cc}0.1953& 0.2310\\ 0.1835& 0.1796\end{array}\right),B=\left(\begin{array}{cc}0.9501& 0.6068\\ 0.2311& 0.4860\end{array}\right).$

We use the basic fixed point iterative method (3.1) to solve Equation (1.1). Let $X_0=I$ , after 14 iterations we get the following result:

$X\approx X_{14}=\left(\begin{array}{cc}1.7834& 0.5628\\ 0.5628& 1.5167\end{array}\right)$

and $\epsilon \left(X_{14}\right)={\Vert X_{14}-I-A^{*}{X}_{14}^{p}A-B^{*}{X}_{14}^{-q}B\Vert }_{\infty }=4.0992\times {10}^{-11}.$

Considering the multi-step stationary iterative method (3.2) with $X_1=X_2=I$ , after 19 iterations we get the following result:

$X\approx X_{19}=\left(\begin{array}{cc}1.7834& 0.5628\\ 0.5628& 1.5167\end{array}\right)$

and $\epsilon \left(X_{19}\right)={\Vert X_{19}-I-A^{*}{X}_{19}^{p}A-B^{*}{X}_{19}^{-q}B\Vert }_{\infty }=3.1824\times {10}^{-11}.$

Example 4.2. Consider Eq.(1.1) with $p=\frac{1}{2},q=\frac{1}{3}$ and

$A=\left(\begin{array}{ccc}0.2913& 0.5185& 0.6154\\ 0.2621& 0.4214& 0.7919\\ 0.0565& 0.0447& 0.8218\end{array}\right),B=\left(\begin{array}{ccc}0.9382& 0.5355& 0.3936\\ 0.4763& 0.0969& 0.9579\\ 0.9057& 0.0013& 0.8929\end{array}\right),$

Considering the basic fixed point iterative method (3.1) with $X_0=I$ , after 13 iterations we get the following result:

$X\approx X_{13}=\left(\begin{array}{ccc}2.3605& 0.8354& 1.8154\\ 0.8354& 2.0042& 1.6716\\ 1.8154& 1.6716& 6.3423\end{array}\right)$

and $\epsilon \left(X_{13}\right)={\Vert X_{13}-I-A^{\mathrm{<mo>\; *\; </mo>}}{X}_{13}^{p}A-B^{\mathrm{<mo>\; *\; </mo>}}{X}_{13}^{-q}B\Vert }_{\infty }=9.7863\times {10}^{-11}\mathrm{.}$

Considering the multi-step stationary iterative method (3.2) with $X_1=I,X_2=2I$ , after 35 iterations we get the following result:

$X\approx X_{35}=\left(\begin{array}{ccc}2.3605& 0.8354& 1.8154\\ 0.8354& 2.0042& 1.6716\\ 1.8154& 1.6716& 6.3423\end{array}\right)$

and $\epsilon \left(X_{35}\right)={\Vert X_{35}-I-A^{*}{X}_{35}^{p}A-B^{*}{X}_{35}^{-q}B\Vert }_{\infty }=5.8041\times {10}^{-11}.$

Acknowledgements

Sincere thanks to the members of ALAMT for their professional performance. The work was supported by the Natural Science Foundation of Shandong Province of China (ZR2015AL017, ZR2014AM032), the Science Reseach Foundation of Heze University (XYKJ07).

Cite this paper

Gao, D.J. (2017) Iterative Methods for Solving the Nonlinear Matrix Equation Advances in Linear Algebra & Matrix Theory, 7, 72-78. https://doi.org/10.4236/alamt.2017.73007

References