Convergence Analysis of Structure-Preserving Doubling Algorithms for Riccati-Type Matrix Equations
|
|
- Prosper Chapman
- 5 years ago
- Views:
Transcription
1 Convergence Analysis of Structure-Preserving Doubling Algorithms for Riccati-Type Matrix Equations Wen-Wei Lin Shu-Fang u Abstract In this paper a structure-preserving transformation of a symplectic pencil is introduced referred to as the doubling transformation and its some basic properties are presented. Based on the nice properties of this kind of transformations a unified convergence theory for the structure-preserving doubling algorithms for solving a class of Riccati-type matrix equations is established by using only the knowledge from elementary matrix theory. Key words. matrix equation structure-preserving doubling algorithm convergence rate AMS subject classfications. 15A24 65H10 93B50 93D15 1 Introduction In this paper we investigate the convergence of the structure-preserving doubling algorithms (SDA) for the computation of the symmetric positive definite or semi-definite solutions to the following Riccati-type matrix equations: Continuous-time algebraic Riccati equation (CARE) G + A T + A + H = 0 (1.1) where AHG R n n with G and H symmetric positive semi-definite. Discrete-time algebraic Riccati equation (DARE) = A T (I + G) 1 A + H (1.2) where AHG R n n with G and H symmetric positive semi-definite. Nonlinear matrix equation with the plus sign (NME-P) 3 where AQ R n n with Q symmetric positive definite. Nonlinear matrix equation with the minus sign (NME-M) 11 where AQ R n n with Q symmetric positive definite. + A T 1 A = Q (1.3) A T 1 A = Q (1.4) Department of Mathematics National Tsing Hua University Hsinchu 300 Taiwan (wwlin@am.nthu.edu.tw). School of Mathematical Sciences Peking University Beijing China (xsf@pku.edu.cn). This research was supported in part by the National Center for Theoretical Sciences in Taiwan. 1
2 The Riccati-type matrix equations occur in many important applications (see and references therein). The nonlinear matrix equations CARE and DARE have been studied extensively (see and 31); and the nonlinear matrix equations NME-P and NME-M has been studied recently by several authors (see and 32). A class of methods referred to as the doubling algorithms has attracted much interests in 70s and 80s (see 2 and references therein). These methods originate from the fixed-point iteration derived from the DARE: k+1 = A T k (I + G k ) 1 A + H. Instead of generating the sequence { k } doubling algorithms generate { 2 k}. Doubling algorithms were largely forgotten in the past decade. Recently they have been revived for DAREs and CAREs because their nice numerical behavior: quadratical convergence rate low cost computational cost per step and good numerical stability (see 6 7 8). Concerning the matrix equations NME-Ps and NME-Ms in 2001 B. Meini proposed an iterative method with the same numerical behavior as the SDA algorithms for solving the DAREs and CAREs (see 26 16). In this paper by employing the same technique as in 8 we derive two SDA algorithms for solving the NME-Ps and NME-Ms as a result we find they are essentially the same as Meini s in 26. In summary we find all the SDA algorithms for solving the four Riccati-type matrix equations can be viewed as applying some special structure-preserving transformations to the associated symplectic pencils repeatedly. Therefore we first introduce a structure-preserving transformation of a symplectic pencil which is referred to as the doubling transformation and then prove that an important feature of this kind of transformations is that it is structure preserving eigenspace preserving and eigenvalue doubling. Finally based on the nice properties of the doubling transformations we develop a unified convergence theory for these SDA algorithms by using only the knowledge from elementary matrix theory. Throughout this pater the symbols 2 denote the matrix spectral norm. For a given n n matrix A we use ρ(a) to denote the spectral radius of A. For real symmetric matrices and Y we write > Y ( Y ) if Y is symmetric positive definite (semi-definite). The paper is organized as follows. In Section 2 we introduce a structure-preserving transformation of a symplectic pencil and develop its basic properties. In Section 3 we do the convergence analysis of the SDA algorithm for solving the DAREs based on the theory established in Section 2. Section 4 derives a SDA algorithm for solving the NME-P (1.3) by using the doubling transformations and establishes the convergence theory of the alorithm. Concluding remarks are given in Section 5. 2 Doubling Transformation In this section we introduce a structure-preserving transformation of a symplectic pencil and develop its some basic properties. We begin with the definition of the transformation. Let M λl R 2n 2n be a symplectic pencil i.e. MJM T = LJL T J = 0 I I 0. (2.1) Define { N(ML) = M L : M L R 2n 2n L } rankm L = 2n M L = 0. (2.2) M L Since rank 2n it follows that N(ML). For any given M M L N(ML) define M = M M L = L L. (2.3) The transformation M λl M λ L 2
3 is called a doubling transformation. An important feature of this kind of transformations is that it is structure preserving eigenspace preserving and eigenvalue doubling as shown in the following theorem. Theorem 2.1. Assume that the pencil M λ L is a doubling transformation of a symplectic pencil M λl. Then we have: (a) The pencil M λ L is still symplectic. U U (b) If M = L S where UV R V V n m and S R m m then M (c) If the pencil M λl has the Kronecker canonical form Jr 0 Ir 0 WMZ = WLZ = 0 I 2n r 0 N 2n r U = V L U S V 2. (2.4) where WZ are nonsingular J r a Jordan matrix corresponding to the finite eigenvalues of M λl and N 2n r a nilpotent Jordan matrix corresponding to the infinite eigenvalues of M λl then there exists a nonsingular matrix Ŵ such that Ŵ MZ J 2 = r 0 Ŵ 0 I LZ Ir 0 = 2n r 0 N2n r 2. (2.5) Proof. (a) Since M L N(ML) implies that M L = L M it follows from (2.1) that MJ M T = M MJM T M T = M LJL T M T = (L M)J(M T L T ) = L LJL T L T = LJ L T which means that the pencil M λ L is symplectic. (b) Using the equality M L = L M again we have U U M L S = L V M S V and hence (c) Let M U = M V M U = M V L U S = L V M U S = L V L Jr 0 Ir 0 ˇM = W Ľ 0 I = W. 2n r 0 N 2n r U S 2 = V L U S 2. V Then we have Jr 0 Jr 0 Ir 0 Jr 0 ˇM LZ = WLZ = = 0 I 2n r 0 I 2n r 0 N 2n r 0 N 2n r Ir 0 Ir 0 Jr 0 Jr 0 Ľ MZ = WMZ = =. 0 N 2n r 0 N 2n r 0 I 2n r 0 N 2n r This together with Z nonsingular implies that ˇM L Ľ = 0. It is clear that rank M ˇM Ľ = M 2n. Notice that (2.4) implies that the pencil M λl is regular it must have rank = 2n. L Thus the rows of M L and ˇM M Ľ form two different bases of the left null space of. L Consequently there exists a nonsingular matrix Ŵ such that ŴM L = ˇM Ľ. It can be easily verified that the matrix Ŵ satisfies (2.5). 3
4 Remark 2.1. A subspace W of R 2n is called a generalized eigenspace of a pencil M λl R 2n 2n U if W is spanned by the columns of W = where UV R V n m W has full column rank and satisfies MW = LWS with S R m m. Therefore (b) of Theorem 2.1 tells us that if W is a generalized eigenspace of a symplectic pencil M λl then it is also a generalized eigenspace of its doubling transformation which is a corner stone for us to establish convergence theory of the SDA algorithms for solving the Riccati-type matrix equations in the next two sections. Remark 2.2. A pencil M λl is called regular if det(m λl) does not vanish identically. It is well-known that a pencil is regular if and only if it has a Kronecker canonical form. Thus (c) of Theorem 2.1 implies that if M λl is regular then its doubling transformation is also regular and that λ is a eigenvalue of M λl if and only if λ 2 is an eigenvalue of its doubling transformation. A symplectic pencil M λl is said to be a first standard symplectic form (SSF-1) if it has the following form A 0 I G M = L = H I 0 A T (2.6) with HG 0; A symplectic pencil M λl is said to be a second standard symplectic form (SSF-2) if it has the following form A 0 P I M = L = Q I A T (2.7) 0 with PQ 0. The following theorem shows that the two standard symplectic forms are preserved by an appropriate choice of the doubling transformations. Theorem 2.2. (a) Let M λl be a SSF-1 form. Then a matrix M L N(ML) can be constructed such that its corresponding doubling transformation M λ L is still a SSF-1. (b) Let M λl be a SSF-2 form. If Q P > 0 and Q A T (Q P) 1 A 0 then a matrix M L N(ML) can be constructed such that its corresponding doubling transformation M λ L is still a SSF-2. L Proof. (a) Applying block Gaussian elimination and row permutation to the matrix we M can get such that M = A(I + GH) 1 0 A T (I + HG) 1 H I I AG(I + HG) 1 L = 0 A T (I + HG) 1 (2.8) M L = L M (2.9) i.e. M L N(ML). Here the Sherman-Morrison-Woodbury formula (see e.g. 12 p. 50) is used. For more details see 8. We then compute L L and M M to produce  0 I Ĝ M = M M = Ĥ I L = L L = 0  T (2.10) where  = A(I + GH) 1 A (2.11) Ĝ = G + AG(I + HG) 1 A T (2.12) Ĥ = H + A T (I + HG) 1 HA. (2.13) It is clear that the resulting pencil is still a SSF-1. (b) Similarly under the condition Q P > 0 we can compute a matrix M L N(ML) with A(Q P) M = 1 0 I A(Q P) 1 A T (Q P) 1 L I = 0 A T (Q P) 1. (2.14) 4
5 Direct calculation gives rise to where  0 M = M M = Q I L = L L = P I  T 0 (2.15)  = A(Q P) 1 A (2.16) Q = Q A T (Q P) 1 A (2.17) P = P + A(Q P) 1 A T. (2.18) The assumption Q A T (Q P) 1 A 0 implies that the resulting pencil is still a SSF-2. Remark 2.3. The proof of Theorem 2.2 provided us with the well defined computation formulae for constructing the special structure preserving doubling transformations which is the base for us to derive the SDA algorithms for solving the Riccati-type matrix equations. 3 SDA Algorithms for Preserving the SSF-1 In this section we shall use the theory established in the last section to develop the convergence theory of the SDA algorithms for for solving the DAREs and CAREs. The SDA algorithms proposed in 8 and 7 will present in the next two subsections by using the doubling transformations defined in the last section. 3.1 SDA Algorithm for Solving the DAREs It is easy to verify that the DARE (1.2) has a symmetric positive semi-definite solution (i.e. 0) if and only if satisfies that I I M = L S (3.1) for some matrix S R n n where M = A 0 H I L = I G 0 A T. (3.2) Notice that the pencil M λl is just a SSF-1 form. Therefore applying the special doubling transformation defined by (2.11) (2.13) repeatedly gives rise to the following structure-preserving doubling algorithm. Algorithm SDA-1. A 0 = A G 0 = G H 0 = H A k+1 = A k (I + G k H k ) 1 A k G k+1 = G k + A k G k (I + H k G k ) 1 A T k H k+1 = H k + A T k (I + H k G k ) 1 H k A k. This is the SDA algorithm described in 8 in which extensive numerical experiments show that this algorithm is efficient and competitive. 5
6 3.2 SDA Algorithm for Solving the CAREs Assume that 0 solves the CARE (1.1). It is well-known that the CARE (1.1) can be rewritten as I I H = R (3.3) where H = A H G A T R = A G. The matrix H is a Hamiltonian matrix i.e. (HJ) T = HJ. Using Cayley transformation with some appropriate γ > 0 we can transform (3.3) into the following form I I M = L S (3.4) where M = H + γi L = H γi S = (R γi) 1 (R + γi). Now assume that we have chosen a γ > 0 such that the matrices A γ = A γi and W γ = A T γ + HA 1 γ G (3.5) are nonsingular. Chu et al. 7 proposed a method for computing γ such that both A γ and W γ are well-conditioned. Let A 1 γ 0 I A 1 T 1 = HA 1 γ GWγ T γ I 2 = 1 0 Wγ 1 (3.6) which are obtained by alternatively applying block Gaussian elimination to the matrices L and M (see 7 for more details). Then direct calculations gives rise to  0 I Ĝ M = T 2 T 1 M = Ĥ I L = T2 T 1 L = 0  T where  = I + 2γW T γ Ĝ = 2γA 1 γ GW 1 γ Ĥ = 2γW 1 γ HA 1 γ. Here the Sherman-Morrison-Woodbury formula is used. Since γ > 0 and HG 0 implies that ĜĤ 0 it follows that the resulting pencil M λ L is a SSF-1. In addition it follows from (3.4) that I M = L I S (3.7) Thus beginning with (3.7) following the same lines as Algorithm SDA-1 for solving the DARE we can construct a matrix sequence to approximate to the unique symmetric positive semi-definite solution to the CARE (1.1). For more details see Convergence Analysis of Algorithm SDA-1 Now we establish the convergence theory of Algorithm SDA-1 based on Theorem 2.1. The main results are listed in the following theorem. Theorem 3.1. Assume that Y 0 satisfy that = A T (I + G) 1 A + H (3.8) Y = AY (I + HY ) 1 A T + G (3.9) 6
7 where HG 0 that is assume that the DARE (3.8) and its dual DARE (3.9) have symmetric positive semi-definite solutions and Y respectively and let S = (I + G) 1 A T = (I + HY ) 1 A T. (3.10) Then the matrix sequences {A k } {G k } and {H k } generated by Algorithm SDA-1 satisfy that (a) A k = (I + G k )S 2k ; (b) H H k H k+1 and H k = (S T ) 2k ( + G k )S 2k (S T ) 2k ( + Y )S 2k ; (3.11) (c) G G k G k+1 Y and Y G k = (T T ) 2k (Y + Y H k Y )T 2k (T T ) 2k (Y + Y Y )T 2k. (3.12) Proof. Notice that UV 0 implies that I+UV is nonsingular and V (I+UV ) 1 (I+UV ) 1 U 0 it follows that Algorithm SDA-1 is well defined and Define H = H 0 H k H k+1 and G = G 0 G k G k+1. (3.13) M k = A k 0 H k I I Gk L k = 0 A T k Then the pencil M k+1 λl k+1 is a doubling transformation of the pencil M k λl k. Since (3.8) implies that I I M 0 = L 0 S (3.14) where S is defined by (3.10) applying (b) of Theorem 2.1 repeatedly we get I I M k = L k S 2k. (3.15) Equating the blocks of (3.15) yields that Combining (3.16) with (3.17) gives rise to A k = (I + G k )S 2k (3.16) H k = A T k S 2k. (3.17) H k = (S T ) 2k ( + G k )S 2k. (3.18) This together with (I + G k ) 0 implies that H k 0 i.e. H k. Similarly (3.9) can be rewritten as M 0 Y I where T is defined by (3.10) and from (3.19) we can derive that. Y T = L 0 (3.19) I Y G k = (T T ) 2k (Y + Y H k Y )T 2k which implies that Y G k. Thus the theorem is proved. 7
8 Let W = L I Y M Z = I I Y I. Noting that M 0 = M L 0 = L and Y 0 it follows from (3.14) and (3.19) that W and Z are nonsingular and satisfy W 1 S 0 MZ = W 1 I 0 LZ =. 0 I 0 T Thus by the spectral feature of a symplectic pencil it follows that if ρ(s) < 1 then it must have ρ(t) = ρ(s) < 1. In addition it is well known that 0 U V implies that U 2 V 2. Consequently from Theorem 3.1 we immediately get the following convergence result of Algorithm SDA-1. Corollary 3.2. Under the hypothesis of Theorem 3.1 if ρ(s) < 1 then we have (a) A k 2 (1 + 2 Y 2 ) S 2k 2 0 as k 0; (b) H k 2 + Y 2 S 2k as k 0; (c) Y G k 2 Y + Y Y 2 T 2k as k 0. Remark 3.1. The similar convergence result were obtained in 8. In contrast to the work of 8 however our analysis is much simpler and our convergence results are much stronger. Remark 3.2. Let G = BR 1 B T 0 with R > 0 and H = C T C 0 in the DARE (3.8) and assume that (AB) is stabilizable and (AC) is detectable. Then it is well known that the DARE (3.8) and its dual DARE (3.9) have symmetric positive semi-definite solutions and Y respectively and ρ(s) < 1 (see for example for details). Thus the conditions in Corollary 3.1 are satisfied. In fact it is easy to verify that if the DARE (3.8) and its dual DARE (3.9) have symmetric positive semi-definite solutions and Y respectively and ρ(s) < 1 with S = (I + G) 1 A then (AB) is stabilizable and (AC) is detectable. Remark 3.3. By Theorem 3.1 the matrix sequences {H k } and {G k } are monotonically increasing and bounded above and hence there exist symmetric positive semi-definite matrices H and Ḡ such that lim H k = H k lim G k = Ḡ. k Corollary 3.2 tells us that if ρ(s) < 1 then = H and Y = Ḡ. Now it is natural to ask if this is still true without the condition ρ(s) < 1. This is a very interesting and worthwhile problem. See the following two simple examples. Example 3.1. Let α1 0 ζ1 0 A = G = H = 0 α 2 0 ζ 2 θ where α 1 α 2 ζ 1 ζ 2 θ 1 R with ζ 2 0 ζ 1 θ 1 > 0 and α2 1 1+ζ 1 + θ 1 = 1. Direct calculation verifies 1 0 that with these given data the DARE (1.2) has = as its symmetric positive semi-definite 0 0 solution and in this case α1 S = (I G) 1 A = 1+ζ α 2 Since α 2 may be any real number ρ(s) can be any positive real number. Applying Algorithm SDA- 1 to these given data we have α (k) A k = 1 0 ζ (k) (k) 0 α (k) G k = 1 0 θ 2 0 ζ (k) H k =
9 where α (0) 1 = α 1 α (k+1) 1 = (α(k) 1 )2 1 + θ (k) α (0) 2 = α 2 α (k+1) 2 = (α (k) 2 )2 1 ζ(k) 1 ζ (0) 1 = ζ 1 ζ (k+1) 1 = ζ (k) 1 (1 + α(k+1) 1 ) ζ (0) 2 = ζ 2 ζ (k+1) 2 = ζ (k) 2 (1 + α(k+1) 2 ) θ (0) 1 = θ 1 θ (k+1) 1 = θ (k) 1 (1 + α(k+1) 1 ). From the proof of Theorem 3.1 we can see that in this case the equality H k = (S T ) 2k ( + G k )S 2k still holds and hence letting λ 1 = α1 1+ζ 1 we have 1 θ (k) 1 = (1 + ζ (k) 1 )λ2k 1 0 (k ) since λ 1 = α1 (k) 1+ζ 1 < 1 and 0 ζ 1 ζ1 θ 1 which can be proved by taking ζ 2 = 0 and using the same method as the proof of (c) of Theorem 3.1. This shows that in this case we still have lim k H k =. However when α 2 > 1 lim k α (k) 2 = +. Example 3.2. Let A = α1 0 G = 0 1 ζ1 0 H = 0 0 θ where α 1 ζ 1 θ 1 R with ζ 1 θ 1 > 0 and α2 1 1+ζ 1 + θ 1 = 1. It easy to verify that with these given 1 0 data for any ξ 0 the matrix ξ = 0 is a symmetric positive semi-definite solution of 0 ξ the DARE (1.2) which means that in this case the DARE (1.2) have infinitely many symmetric positive semi-definite solutions. Simple calculation gives rise to α1 S = (I G ξ ) 1 A = 1+ζ and so we have that ρ(s) = 1 for any ξ 0. We can easily prove that in this case we still have 1 0 ν 0 lim H k = 0 = lim k 0 0 G k = Y = k 0 0 where H k and G k are generated by applying Algorithm SDA-1 to these given data Y is a symmetric positive semi-definite solution of the dual DARE with ν = ζ1 θ 1. 4 SDA Algorithms for Preserving the SSF-2 In this section we shall first use the doubling transformations defined in the last section to derive two SDA algorithms for solving the NME-Ps and NME-Ms. Then we shall use the theory established in the last section to develop the convergence theory of these SDA algorithms. 4.1 SDA Algorithm for Solving the NME-Ps It is easy to verify that the NME-P (1.3) has a symmetric positive definite (i.e. > 0) if and only if satisfies I I M = L S (4.1) 9
10 for some matrix S R n n where A 0 M = Q I L = 0 I A T 0. (4.2) Notice that the pencil M λl is just a SSF-2. Therefore applying the special doubling transformation defined by (2.16) (2.18) repeatedly gives rise to the following structure-preserving doubling algorithm. Algorithm SDA-2. A 0 = A Q 0 = Q P 0 = 0 A k+1 = A k (Q k P k ) 1 A k Q k+1 = Q k A T k (Q k P k ) 1 A k P k+1 = P k + A k (Q k P k ) 1 A T k. Remark 4.1. Of course to ensure that this iteration is well defined the matrix Q k P k must be symmetric positive definite for all k. Below we shall prove that this condition can be guaranteed if the NME-P (1.3) has a symmetric positive solution (see Theorem 4.1). Remark 4.2. It is interesting to note that this algorithm is essentially the same as proposed in 26. In other words Meini s algorithm in 26 is a SDA algorithm. It has pointed out that this algorithm has very nice numerical behavior: having quadratical convergence rate low computational cost per step and good numerical stability. For more details see SDA Algorithm for Solving the NEM-Ms It is proved in 11 that there always exists a unique positive definite solution to the NME-M A T 1 A = Q (4.3) and moreover the spectral radius of 1 A is strictly less than 1. The solution is closely related the generalized eigenspace of the following pencil A 0 0 I M λl = λ Q I A T. (4.4) 0 In fact it is easy to verify that a symmetric positive definite matrix is a solution to the NME-M (4.3) if and only if satisfies that I I M = L S (4.5) for some matrix S R n n. Although the pencil (4.4) is not symplectic we can use the same technique as described in Section 2 to transform it into a symplectic pencil. Take then we have Direct calculation leads to M = AQ 1 0 A T Q 1 I  0 M 0 = M M = Q I L = I AQ 1 0 A T Q 1 M L = L M. (4.6) L0 = L L = P I  T 0 where  = AQ 1 A Q = Q + A T Q 1 A P = AQ 1 A T. (4.7) 10
11 Direct verification shows that M 0 λ L 0 is a symplectic pencil but neither a SSF-1 nor a SSF-2. Assume that > 0 is the unique symmetric positive solution to the NME-M (4.3). Then it satisfies the equality (4.5) with S = 1 A. Similar to the proof of (b) of Theorem 2.1 we can show that I M 0 = L I 0 S 2. (4.8) Now Let I I 0 = P I I Then it follows from (4.8) that M = I M Â 0 Q + P I L = 0 I Â T 0. I = L S 2. (4.9) Clearly now the pencil M λ L is a SSF-2. Thus beginning with (4.9) following the same lines as Algorithm SDA-2 for solving the NME-P (1.3) we can construct a matrix sequence to approximate to. Then the unique symmetric positive definite solution to the NME-M (4.3) can be obtained by computing = P. 4.3 Convergence Analysis of Algorithm SDA-2 Now we establish the convergence theory of Algorithm SDA-2 based on Theorem 2.1. The main results are listed in the following theorem. Theorem 4.1. Assume that > 0 satisfy that + A T 1 A = Q (4.10) where Q > 0 and let S = 1 A. Then the matrix sequences {A k } {Q k } and {P k } generated by Algorithm SDA-2 satisfy that (a) A k = ( P k )S 2k ; (b) 0 P k P k+1 < and (c) Q k+1 Q k Q and Proof. (Apply Mathematical Induction.) Denote Ak 0 M k = I Q k P k = ( P k ) + A T k ( P k ) 1 A k > 0; (4.11) Q k = (S T ) 2k ( P k )S 2k (S T ) 2k S 2k. (4.12) Q k L k = Pk I A T k 0 where P 0 = 0. For k = 1 since Q 0 P 0 = Q > 0 it follows that A 1 Q 1 P 1 are all well defined. Using the equality (4.10) we have A I 0 0 I A T = 1 A Q A T 1 > 0. (4.13) I 0 0 I Simple computation yields that I AQ 1 A I 0 0 I A T Q Q 1 A T = I Combining (4.14) and (4.13) we get AQ 1 A T 0. (4.14) 0 Q P 1 = AQ 1 A T > 0. (4.15) 11
12 On the other hand From (4.10) it is easy to verify that satisfies I I M 0 = L 0 S with S = 1 A. Since M 1 λl 1 is a doubling transformation of M 0 λl 0 Applying (b) of Theorem 2.1 we get I I M 1 = L 1 S 2. (4.16) Equating the blocks of (4.16) gives rise to This together with (4.15) implies that A 1 = ( P 1 )S 2 Q 1 = A T 1 S 2. Q 1 P 1 = ( P 1 ) + A T 1 ( P 1 ) 1 A 1 > 0 (4.17) Q 1 = (S T ) 2 ( P 1 )S 2 0. (4.18) Obviously the inequalities Q = Q 0 Q 1 and 0 = P 0 P 1 hold. Thus we have proved that the theorem is true for k = 1 Next assume that the theorem is true for all positive integers less or equal to k. Consider the case of k + 1. Since Q k P k > 0 it follows that A k+1 Q k+1 P k+1 are all well defined. Similar to the proof of (4.15) using the equality (4.11) we can prove that P k+1 = ( P k ) A k (Q k P k ) 1 A T k > 0. On the other hand since M j+1 λl j+1 is a doubling transformation of M j λl j for j = 01...k By using (b) of Theorem 2.1 k + 1 times we get I I M k+1 = L k+1 S 2k+1. (4.19) From (4.19) following the same lines as the proof of (4.17) and (4.18) it can be proved that Q k+1 P k+1 = ( P k+1 ) + A T k+1( P k+1 ) 1 A k+1 > 0 Q k+1 = (S T ) 2k+1 ( P k+1 )S 2k+1 0. Clearly P k P k+1 and Q k Q k+1 are true. This shows that the theorem is also true for integer k + 1. By induction principle the theorem is true for all positive integers. Remark 4.3. The similar results were obtained in 26 by using some properties of cyclic reduction and the spectral properties of block Toeplitz matrices having nonnegative definite matrix-valued generating functions. In contrast to the work of 26 however our analysis is much simpler and the results are much stronger. It was proved in 10 that if the matrix equation (1.3) has a symmetric positive definite solution then all symmetric solutions are positive definite and it has a maximal and minimal solution + and respectively. Since Theorem 4.1 is true for any symmetric positive definite solution to the NME-P (1.3) the following result follows immediately. Corollary 4.2. Under the hypothesis of Theorem 4.1 we have Q k P k > + 0 for all k where + and are the maximal and minimal solution of (1.3) respectively. In addition from Theorem 3.1 we immediately get the following result. 12
13 Corollary 4.3. Under the hypothesis of Theorem 4.1 if ρ(s) < 1 then we have (a) A k 2 2 S 2k 2 0 as k 0; (b) Q k 2 2 S 2k as k 0. Remark 4.4. By Theorem 4.1 the matrix sequence {Q k } is monotonically decreasing and bounded below by > 0 and hence there exists Q > 0 such that lim k Q k = Q. Corollary 4.3 tells us that if ρ(s) < 1 then = Q. In fact ρ(s) < 1 implies that is the maximal solution of (1.3) and moreover it has been proved that is the maximal solution of (1.3) if and only if ρ(s) 1 (see 16). Now assume that = + is the maximal solution of (1.3) it is natural to ask whether Q = + still holds if ρ(s) = 1. In 16 Guo proved that if ρ(s) = 1 and all eigenvalues of S on the unit circle are semisimple then Q = + is sill true and in this case the convergence is at least linear with rate 1/2. When S has nonsemisimple unimodular eigenvalues it is still open whether Q = +. Remark 4.5. It is proved that the NME-P (1.3) has a symmetric positive definite solution If and only if the nonlinear matrix equation Y + AY 1 A T = Q (4.20) has a symmetric positive solution Y (see for example 26). Assume that the maximal solution of (4.20) is Y +. Then it follows from (4.20) that A 0 I 0 I I T = Q I Q Y + A T (4.21) 0 Q Y + where T = Y 1 + A T. Similar to the proof of (4.12) from (4.21) we can show that 0 Q Y + P k = (T T ) 2k (Q k Q + Y + )T 2k (T T ) 2k Y + T 2k where P k and Q k are generated by Algorithm SDA-2. Since ρ(t) = ρ(y 1 + A T ) = ρ( 1 + A) (see for example 26) where + is the maximal solution of the NME-P (1.3) it follows that under the conditions of Corollary 4.3 we have lim k P k = Q Y +. If A is nonsingular then = Q Y + (see 26) where is the minimal solution of the NME-P (1.3) and so in this case we have lim k P k =. Remark 4.6. Since lim k (Q k P k ) = + if A is nonsingular and ρ(s) < 1 the under bound + in Corollary 4.2 is the best one. However + may be singular indeed it can be zero matrix. For example the NME-P (1.3) with Q = I and A = 1 2 I has + = = 1 2 I. 5 Conclusions In this paper we have introduced a structure-preserving transformation of a symplectic pencil referred to as the doubling transformation and developed its basic properties. Based on these nice properties of the transformation a unified convergence theory for the structure-preserving doubling algorithms for solving a class of Riccati-type matrix equations has been established by using only the knowledge from elementary matrix theory. References 1 G. Ammar and V. Mehrmann On Hamiltonian and symplectic Hessenberg forms Lin. Alg. Appl. 149 (1991) pp B. Anderson Second-order convergent algorithms for the steady-state Riccati equation Int. J. Control 28 (1978) pp W.N. Anderson T.D. Morley and G.E. Trapp Positive solutions to = A B 1 B Lin. Alg. Appl. 134 (1990) pp
14 4 Z. Bai J. Demmel and M. Gu An inverse free parallel spectral divide and conquer algorithm for nonsymmetric eigenproblems Numerische Mathematik 76 (1997) pp P. Benner and R. Byers Evaluating products of matrix pencils and collapsing matrix products Num. Lin. Alg. Appl. 8 (2001) pp E.K.-W. Chu H.-Y. Fan and W.-W. Lin A generalized structure-preserving doubling algorithm for generalized discrete-time algebraic Riccati equations preprint NCTS National Tsing Hua University Hsinchu 300 Taiwan E.K.-W. Chu H.-Y. Fan and W.-W. Lin A structure-preserving doubling algorithm for cotinuous-time algebraic Riccati equations to appear in Lin. Alg. Appl.. 8 E.K.-W. Chu H.-Y. Fan W.-W. Lin and C.-S. Wang A structure-preserving doubling algorithm for periodic descrete-time algebraic Riccati equations preprint NCTS National Tsing Hua University Hsinchu 300 Taiwan J.C. Engwerda On the existence of a positive definite solution of the matrix equation + A T 1 A = I Lin. Alg. Appl. 194 (1993) pp J.C. Engwerda A.C.M. Ran and A.L. Rijkeboer Necessary and sufficient conditions for the existence of a positive definite solution of the matrix equation + A 1 A = Q Lin. Alg. Appl. 186 (1993) pp A. Ferrante and B.C. Levy Hermitian solutions of the Equation = Q+N 1 N Lin. Alg. Appl. 247 (1996) pp G.H. Golub and C.F. Van Loan Matrix Computations 3rd ed. The Johns Hopkins University Press T. Gudmundsson C. Kenney and A. J. Laub Scaling of the discrete-time algebraic Riccati equation to enhance stability of the Schur solution method IEEE Trans. Auto. Control 37 (1992) pp C.-H. Guo Newton s method for discrete algebraic Riccati equations when the closed-loop matrix has eigenvalues on the unit circle SIAM J. Matrix Anal. Appl. 20 (1998) pp C.-H. Guo and P. Lancaster Iterative solution of two matrix equations Math. Comp. 68(1999) pp C.-H. Guo Comvergence rate of an iterative method for a nonlinear matrix equation SIAM J. Matrix Anal. Appl. 23 (2001) pp J.J. Hench and A.J. Laub Numerical solution of the discrete-time periodic Riccati equation IEEE Trans. Auto. Control 39 (1994) pp M. Kimura Convergence of the doubling algorithm for the discrete-time algebraic Riccati equation Int. J. Systems Science 19 (1988) pp D. Lainiotis N. Assimakis and S. Katsikas New doubling algorithm for the discrete periodic Riccati equation Applied Mathematics and Computation 60 (1994) pp P. Lancaster and L. Rodman Algebraic Riccati Equations Oxford University Press Oxford A.J. Laub A Schur method for solving algebraic Riccati equations IEEE Trans. Auto. Control 24 (1979) pp A.J. Laub Invariant subspace methods for numerical solution of Riccati equation in The Raccati Equation S. Bittanti A.J. Laub and J.C. Willems eds Springer-Verlag Berl 14
15 23 L.-Z. Lu and W.-W. Lin An iterative algorithm for the solution of the discrete time algebraic Riccati equations Lin. Alg. Appl. 189 (1993) pp L.-Z. Lu W.-W. Lin and C. E. M. Pearce An efficient algorithm for the discrete-time algebraic Riccati equation IEEE Trans. Auto. Control 44 (1999) pp V. Mehrmann The Autonomous Linear Quadratic Control Problem Lecture Notes in Control and Information Sciences 163 Springer-Verlag Berlin B. Meini Efficient computation of the extreme solutions of + A 1 A = Q and A 1 A = Q Math. Comp. 71 (2001) pp C. Paige and C. Van Loan A Schur decomposition for Hamiltonian matrices Lin. Alg. Appl. 41 (1981) pp T. Pappas A. J. Laub and N. R. Sandell On the numerical solution of the discrete-time algebraic Riccati equation IEEE Trans. Auto. Control 25 (1980) pp J.-G. Sun Sensitivity analysis of the discrete-time algebraic Riccati equation Lin. Alg. Appl. 275/276 (1998) pp J.-G. Sun and S.F. u Perturbation analysis of the maximal solution of the matrix equation + A T 1 A = P.II Lin. Alg. Appl. 362 (2003) pp P. Van Dooren A generalized eigenvalue approach for solving Riccati equations SIAM J. Sci. & Statistical Computing 2 (1981) pp S.F. u On the maximal solution of the matrix equation +A T 1 A = I Acta Scientiarum Naturalium Universitatis Pekinensis 36(2000) pp
Recursive Solutions of the Matrix Equations X + A T X 1 A = Q and X A T X 1 A = Q
Applied Mathematical Sciences, Vol. 2, 2008, no. 38, 1855-1872 Recursive Solutions of the Matrix Equations X + A T X 1 A = Q and X A T X 1 A = Q Nicholas Assimakis Department of Electronics, Technological
More informationA generalized structure-preserving doubling algorithm for generalized discrete-time algebraic Riccati equations
A generalized structure-preserving doubling algorithm for generalized discrete-time algebraic Riccati equations T.-M. Hwang E. K.-W. Chu W.-W. Lin March 1, 2005 Abstract In 10, an efficient structure-preserving
More informationON THE SEMIGROUP OF STANDARD SYMPLECTIC MATRICES AND ITS APPLICATIONS
ON THE SEMIGROUP OF STANDARD SYMPLECTIC MATRICES AND ITS APPLICATIONS M. CHU, N. DEL BUONO, F. DIELE, T. POLITI, AND S. RAGNI Abstract. A matrix Z[ R 2n 2n ] is said [ to be in] the standard symplectic
More informationThe Solvability Conditions for the Inverse Eigenvalue Problem of Hermitian and Generalized Skew-Hamiltonian Matrices and Its Approximation
The Solvability Conditions for the Inverse Eigenvalue Problem of Hermitian and Generalized Skew-Hamiltonian Matrices and Its Approximation Zheng-jian Bai Abstract In this paper, we first consider the inverse
More informationIterative Solution of a Matrix Riccati Equation Arising in Stochastic Control
Iterative Solution of a Matrix Riccati Equation Arising in Stochastic Control Chun-Hua Guo Dedicated to Peter Lancaster on the occasion of his 70th birthday We consider iterative methods for finding the
More informationTheorem Let X be a symmetric solution of DR(X) = 0. Let X b be a symmetric approximation to Xand set V := X? X. b If R b = R + B T XB b and I + V B R
Initializing Newton's method for discrete-time algebraic Riccati equations using the buttery SZ algorithm Heike Fabender Zentrum fur Technomathematik Fachbereich - Mathematik und Informatik Universitat
More informationOn the Hermitian solutions of the
Journal of Applied Mathematics & Bioinformatics vol.1 no.2 2011 109-129 ISSN: 1792-7625 (print) 1792-8850 (online) International Scientific Press 2011 On the Hermitian solutions of the matrix equation
More informationPositive Solution to a Generalized Lyapunov Equation via a Coupled Fixed Point Theorem in a Metric Space Endowed With a Partial Order
Filomat 29:8 2015, 1831 1837 DOI 10.2298/FIL1508831B Published by Faculty of Sciences Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat Positive Solution to a Generalized
More informationThe Kalman-Yakubovich-Popov Lemma for Differential-Algebraic Equations with Applications
MAX PLANCK INSTITUTE Elgersburg Workshop Elgersburg February 11-14, 2013 The Kalman-Yakubovich-Popov Lemma for Differential-Algebraic Equations with Applications Timo Reis 1 Matthias Voigt 2 1 Department
More informationAn SVD-Like Matrix Decomposition and Its Applications
An SVD-Like Matrix Decomposition and Its Applications Hongguo Xu Abstract [ A matrix S C m m is symplectic if SJS = J, where J = I m I m. Symplectic matrices play an important role in the analysis and
More informationThe iterative methods for solving nonlinear matrix equation X + A X 1 A + B X 1 B = Q
Vaezzadeh et al. Advances in Difference Equations 2013, 2013:229 R E S E A R C H Open Access The iterative methods for solving nonlinear matrix equation X + A X A + B X B = Q Sarah Vaezzadeh 1, Seyyed
More informationOn matrix equations X ± A X 2 A = I
Linear Algebra and its Applications 326 21 27 44 www.elsevier.com/locate/laa On matrix equations X ± A X 2 A = I I.G. Ivanov,V.I.Hasanov,B.V.Minchev Faculty of Mathematics and Informatics, Shoumen University,
More informationDefinite versus Indefinite Linear Algebra. Christian Mehl Institut für Mathematik TU Berlin Germany. 10th SIAM Conference on Applied Linear Algebra
Definite versus Indefinite Linear Algebra Christian Mehl Institut für Mathematik TU Berlin Germany 10th SIAM Conference on Applied Linear Algebra Monterey Bay Seaside, October 26-29, 2009 Indefinite Linear
More informationBOUNDS OF MODULUS OF EIGENVALUES BASED ON STEIN EQUATION
K Y BERNETIKA VOLUM E 46 ( 2010), NUMBER 4, P AGES 655 664 BOUNDS OF MODULUS OF EIGENVALUES BASED ON STEIN EQUATION Guang-Da Hu and Qiao Zhu This paper is concerned with bounds of eigenvalues of a complex
More informationThe Eigenvalue Shift Technique and Its Eigenstructure Analysis of a Matrix
The Eigenvalue Shift Technique and Its Eigenstructure Analysis of a Matrix Chun-Yueh Chiang Center for General Education, National Formosa University, Huwei 632, Taiwan. Matthew M. Lin 2, Department of
More informationYONGDO LIM. 1. Introduction
THE INVERSE MEAN PROBLEM OF GEOMETRIC AND CONTRAHARMONIC MEANS YONGDO LIM Abstract. In this paper we solve the inverse mean problem of contraharmonic and geometric means of positive definite matrices (proposed
More informationOn a quadratic matrix equation associated with an M-matrix
Article Submitted to IMA Journal of Numerical Analysis On a quadratic matrix equation associated with an M-matrix Chun-Hua Guo Department of Mathematics and Statistics, University of Regina, Regina, SK
More informationISOLATED SEMIDEFINITE SOLUTIONS OF THE CONTINUOUS-TIME ALGEBRAIC RICCATI EQUATION
ISOLATED SEMIDEFINITE SOLUTIONS OF THE CONTINUOUS-TIME ALGEBRAIC RICCATI EQUATION Harald K. Wimmer 1 The set of all negative-semidefinite solutions of the CARE A X + XA + XBB X C C = 0 is homeomorphic
More informationNumerical Methods in Matrix Computations
Ake Bjorck Numerical Methods in Matrix Computations Springer Contents 1 Direct Methods for Linear Systems 1 1.1 Elements of Matrix Theory 1 1.1.1 Matrix Algebra 2 1.1.2 Vector Spaces 6 1.1.3 Submatrices
More informationFinite dimensional indefinite inner product spaces and applications in Numerical Analysis
Finite dimensional indefinite inner product spaces and applications in Numerical Analysis Christian Mehl Technische Universität Berlin, Institut für Mathematik, MA 4-5, 10623 Berlin, Germany, Email: mehl@math.tu-berlin.de
More informationOn the Local Convergence of an Iterative Approach for Inverse Singular Value Problems
On the Local Convergence of an Iterative Approach for Inverse Singular Value Problems Zheng-jian Bai Benedetta Morini Shu-fang Xu Abstract The purpose of this paper is to provide the convergence theory
More informationImproved Newton s method with exact line searches to solve quadratic matrix equation
Journal of Computational and Applied Mathematics 222 (2008) 645 654 wwwelseviercom/locate/cam Improved Newton s method with exact line searches to solve quadratic matrix equation Jian-hui Long, Xi-yan
More informationIndex. for generalized eigenvalue problem, butterfly form, 211
Index ad hoc shifts, 165 aggressive early deflation, 205 207 algebraic multiplicity, 35 algebraic Riccati equation, 100 Arnoldi process, 372 block, 418 Hamiltonian skew symmetric, 420 implicitly restarted,
More informationOn the computation of the Jordan canonical form of regular matrix polynomials
On the computation of the Jordan canonical form of regular matrix polynomials G Kalogeropoulos, P Psarrakos 2 and N Karcanias 3 Dedicated to Professor Peter Lancaster on the occasion of his 75th birthday
More informationNonlinear palindromic eigenvalue problems and their numerical solution
Nonlinear palindromic eigenvalue problems and their numerical solution TU Berlin DFG Research Center Institut für Mathematik MATHEON IN MEMORIAM RALPH BYERS Polynomial eigenvalue problems k P(λ) x = (
More informationA CONVERGENCE RESULT FOR MATRIX RICCATI DIFFERENTIAL EQUATIONS ASSOCIATED WITH M-MATRICES
A CONVERGENCE RESULT FOR MATRX RCCAT DFFERENTAL EQUATONS ASSOCATED WTH M-MATRCES CHUN-HUA GUO AND BO YU Abstract. The initial value problem for a matrix Riccati differential equation associated with an
More informationarxiv: v1 [math.na] 6 Jan 2011
THE SDA METHOD FOR NUMERICAL SOLUTION OF LUR E EQUATIONS FEDERICO POLONI AND TIMO REIS arxiv:1101.1231v1 math.na 6 Jan 2011 Abstract. We introduce a numerical method for the numerical solution of the so-called
More informationA Note on Simple Nonzero Finite Generalized Singular Values
A Note on Simple Nonzero Finite Generalized Singular Values Wei Ma Zheng-Jian Bai December 21 212 Abstract In this paper we study the sensitivity and second order perturbation expansions of simple nonzero
More informationMatrix Mathematics. Theory, Facts, and Formulas with Application to Linear Systems Theory. Dennis S. Bernstein
Matrix Mathematics Theory, Facts, and Formulas with Application to Linear Systems Theory Dennis S. Bernstein PRINCETON UNIVERSITY PRESS PRINCETON AND OXFORD Contents Special Symbols xv Conventions, Notation,
More informationSOLVING A STRUCTURED QUADRATIC EIGENVALUE PROBLEM BY A STRUCTURE-PRESERVING DOUBLING ALGORITHM
SOLVING A STRUCTURED QUADRATIC EIGENVALUE PROBLEM BY A STRUCTURE-PRESERVING DOUBLING ALGORITHM CHUN-HUA GUO AND WEN-WEI LIN Abstract In studying the vibration of fast trains, we encounter a palindromic
More informationDOUBLING ALGORITHMS WITH PERMUTED LAGRANGIAN GRAPH BASES
DOUBLNG ALGORTHMS WTH PERMUTED LAGRANGAN GRAPH BASES VOLKER MEHRMANN AND FEDERCO POLON Abstract. We derive a new representation of Lagrangian subspaces in the form m Π T, X where Π is a symplectic matrix
More informationDampening Controllers via a Riccati Equation. Approach. Pod vodarenskou vez Praha 8. Czech Republic. Fakultat fur Mathematik
Dampening Controllers via a Riccati Equation Approach J.J. Hench y C. He z V. Kucera y V. Mehrmann z y Institute of Information Theory and Automation Academy of Sciences of the Czech Republic Pod vodarenskou
More informationRITZ VALUE BOUNDS THAT EXPLOIT QUASI-SPARSITY
RITZ VALUE BOUNDS THAT EXPLOIT QUASI-SPARSITY ILSE C.F. IPSEN Abstract. Absolute and relative perturbation bounds for Ritz values of complex square matrices are presented. The bounds exploit quasi-sparsity
More informationCyclic reduction and index reduction/shifting for a second-order probabilistic problem
Cyclic reduction and index reduction/shifting for a second-order probabilistic problem Federico Poloni 1 Joint work with Giang T. Nguyen 2 1 U Pisa, Computer Science Dept. 2 U Adelaide, School of Math.
More informationOn the Skeel condition number, growth factor and pivoting strategies for Gaussian elimination
On the Skeel condition number, growth factor and pivoting strategies for Gaussian elimination J.M. Peña 1 Introduction Gaussian elimination (GE) with a given pivoting strategy, for nonsingular matrices
More informationDiscreteness of Transmission Eigenvalues via Upper Triangular Compact Operators
Discreteness of Transmission Eigenvalues via Upper Triangular Compact Operators John Sylvester Department of Mathematics University of Washington Seattle, Washington 98195 U.S.A. June 3, 2011 This research
More informationMATLAB TOOLS FOR SOLVING PERIODIC EIGENVALUE PROBLEMS 1
MATLAB TOOLS FOR SOLVING PERIODIC EIGENVALUE PROBLEMS 1 Robert Granat Bo Kågström Daniel Kressner,2 Department of Computing Science and HPC2N, Umeå University, SE-90187 Umeå, Sweden. {granat,bokg,kressner}@cs.umu.se
More informationS.F. Xu (Department of Mathematics, Peking University, Beijing)
Journal of Computational Mathematics, Vol.14, No.1, 1996, 23 31. A SMALLEST SINGULAR VALUE METHOD FOR SOLVING INVERSE EIGENVALUE PROBLEMS 1) S.F. Xu (Department of Mathematics, Peking University, Beijing)
More informationELA
POSITIVE DEFINITE SOLUTION OF THE MATRIX EQUATION X = Q+A H (I X C) δ A GUOZHU YAO, ANPING LIAO, AND XUEFENG DUAN Abstract. We consider the nonlinear matrix equation X = Q+A H (I X C) δ A (0 < δ 1), where
More informationFoundations of Matrix Analysis
1 Foundations of Matrix Analysis In this chapter we recall the basic elements of linear algebra which will be employed in the remainder of the text For most of the proofs as well as for the details, the
More informationThe semi-convergence of GSI method for singular saddle point problems
Bull. Math. Soc. Sci. Math. Roumanie Tome 57(05 No., 04, 93 00 The semi-convergence of GSI method for singular saddle point problems by Shu-Xin Miao Abstract Recently, Miao Wang considered the GSI method
More informationPositive systems in the behavioral approach: main issues and recent results
Positive systems in the behavioral approach: main issues and recent results Maria Elena Valcher Dipartimento di Elettronica ed Informatica Università dipadova via Gradenigo 6A 35131 Padova Italy Abstract
More informationCHUN-HUA GUO. Key words. matrix equations, minimal nonnegative solution, Markov chains, cyclic reduction, iterative methods, convergence rate
CONVERGENCE ANALYSIS OF THE LATOUCHE-RAMASWAMI ALGORITHM FOR NULL RECURRENT QUASI-BIRTH-DEATH PROCESSES CHUN-HUA GUO Abstract The minimal nonnegative solution G of the matrix equation G = A 0 + A 1 G +
More informationNumerical Linear Algebra And Its Applications
Numerical Linear Algebra And Its Applications Xiao-Qing JIN 1 Yi-Min WEI 2 August 29, 2008 1 Department of Mathematics, University of Macau, Macau, P. R. China. 2 Department of Mathematics, Fudan University,
More informationPerturbation of Palindromic Eigenvalue Problems
Numerische Mathematik manuscript No. (will be inserted by the editor) Eric King-wah Chu Wen-Wei Lin Chern-Shuh Wang Perturbation of Palindromic Eigenvalue Problems Received: date / Revised version: date
More informationStability and Inertia Theorems for Generalized Lyapunov Equations
Published in Linear Algebra and its Applications, 355(1-3, 2002, pp. 297-314. Stability and Inertia Theorems for Generalized Lyapunov Equations Tatjana Stykel Abstract We study generalized Lyapunov equations
More informationKey words. Strongly eventually nonnegative matrix, eventually nonnegative matrix, eventually r-cyclic matrix, Perron-Frobenius.
May 7, DETERMINING WHETHER A MATRIX IS STRONGLY EVENTUALLY NONNEGATIVE LESLIE HOGBEN 3 5 6 7 8 9 Abstract. A matrix A can be tested to determine whether it is eventually positive by examination of its
More informationPERTURBATION ANAYSIS FOR THE MATRIX EQUATION X = I A X 1 A + B X 1 B. Hosoo Lee
Korean J. Math. 22 (214), No. 1, pp. 123 131 http://dx.doi.org/1.11568/jm.214.22.1.123 PERTRBATION ANAYSIS FOR THE MATRIX EQATION X = I A X 1 A + B X 1 B Hosoo Lee Abstract. The purpose of this paper is
More informationInterlacing Inequalities for Totally Nonnegative Matrices
Interlacing Inequalities for Totally Nonnegative Matrices Chi-Kwong Li and Roy Mathias October 26, 2004 Dedicated to Professor T. Ando on the occasion of his 70th birthday. Abstract Suppose λ 1 λ n 0 are
More informationPreliminary/Qualifying Exam in Numerical Analysis (Math 502a) Spring 2012
Instructions Preliminary/Qualifying Exam in Numerical Analysis (Math 502a) Spring 2012 The exam consists of four problems, each having multiple parts. You should attempt to solve all four problems. 1.
More informationTwo Results About The Matrix Exponential
Two Results About The Matrix Exponential Hongguo Xu Abstract Two results about the matrix exponential are given. One is to characterize the matrices A which satisfy e A e AH = e AH e A, another is about
More informationETNA Kent State University
C 8 Electronic Transactions on Numerical Analysis. Volume 17, pp. 76-2, 2004. Copyright 2004,. ISSN 1068-613. etnamcs.kent.edu STRONG RANK REVEALING CHOLESKY FACTORIZATION M. GU AND L. MIRANIAN Abstract.
More informationp q m p q p q m p q p Σ q 0 I p Σ 0 0, n 2p q
A NUMERICAL METHOD FOR COMPUTING AN SVD-LIKE DECOMPOSITION HONGGUO XU Abstract We present a numerical method to compute the SVD-like decomposition B = QDS 1 where Q is orthogonal S is symplectic and D
More informationThe symmetric linear matrix equation
Electronic Journal of Linear Algebra Volume 9 Volume 9 (00) Article 8 00 The symmetric linear matrix equation Andre CM Ran Martine CB Reurings mcreurin@csvunl Follow this and additional works at: http://repositoryuwyoedu/ela
More informationTHE PERTURBATION BOUND FOR THE SPECTRAL RADIUS OF A NON-NEGATIVE TENSOR
THE PERTURBATION BOUND FOR THE SPECTRAL RADIUS OF A NON-NEGATIVE TENSOR WEN LI AND MICHAEL K. NG Abstract. In this paper, we study the perturbation bound for the spectral radius of an m th - order n-dimensional
More informationSingular-value-like decomposition for complex matrix triples
Singular-value-like decomposition for complex matrix triples Christian Mehl Volker Mehrmann Hongguo Xu December 17 27 Dedicated to William B Gragg on the occasion of his 7th birthday Abstract The classical
More informationSensitivity analysis of the differential matrix Riccati equation based on the associated linear differential system
Advances in Computational Mathematics 7 (1997) 295 31 295 Sensitivity analysis of the differential matrix Riccati equation based on the associated linear differential system Mihail Konstantinov a and Vera
More informationSpectral inequalities and equalities involving products of matrices
Spectral inequalities and equalities involving products of matrices Chi-Kwong Li 1 Department of Mathematics, College of William & Mary, Williamsburg, Virginia 23187 (ckli@math.wm.edu) Yiu-Tung Poon Department
More informationOPTIMAL SCALING FOR P -NORMS AND COMPONENTWISE DISTANCE TO SINGULARITY
published in IMA Journal of Numerical Analysis (IMAJNA), Vol. 23, 1-9, 23. OPTIMAL SCALING FOR P -NORMS AND COMPONENTWISE DISTANCE TO SINGULARITY SIEGFRIED M. RUMP Abstract. In this note we give lower
More informationOn Solving Large Algebraic. Riccati Matrix Equations
International Mathematical Forum, 5, 2010, no. 33, 1637-1644 On Solving Large Algebraic Riccati Matrix Equations Amer Kaabi Department of Basic Science Khoramshahr Marine Science and Technology University
More informationWe first repeat some well known facts about condition numbers for normwise and componentwise perturbations. Consider the matrix
BIT 39(1), pp. 143 151, 1999 ILL-CONDITIONEDNESS NEEDS NOT BE COMPONENTWISE NEAR TO ILL-POSEDNESS FOR LEAST SQUARES PROBLEMS SIEGFRIED M. RUMP Abstract. The condition number of a problem measures the sensitivity
More informationUMIACS-TR July CS-TR 2721 Revised March Perturbation Theory for. Rectangular Matrix Pencils. G. W. Stewart.
UMIAS-TR-9-5 July 99 S-TR 272 Revised March 993 Perturbation Theory for Rectangular Matrix Pencils G. W. Stewart abstract The theory of eigenvalues and eigenvectors of rectangular matrix pencils is complicated
More informationSTRUCTURE PRESERVING DEFLATION OF INFINITE EIGENVALUES IN STRUCTURED PENCILS
Electronic Transactions on Numerical Analysis Volume 44 pp 1 24 215 Copyright c 215 ISSN 168 9613 ETNA STRUCTURE PRESERVING DEFLATION OF INFINITE EIGENVALUES IN STRUCTURED PENCILS VOLKER MEHRMANN AND HONGGUO
More informationThe quadratic eigenvalue problem (QEP) is to find scalars λ and nonzero vectors u satisfying
I.2 Quadratic Eigenvalue Problems 1 Introduction The quadratic eigenvalue problem QEP is to find scalars λ and nonzero vectors u satisfying where Qλx = 0, 1.1 Qλ = λ 2 M + λd + K, M, D and K are given
More informationON THE REDUCTION OF A HAMILTONIAN MATRIX TO HAMILTONIAN SCHUR FORM
ON THE REDUCTION OF A HAMILTONIAN MATRIX TO HAMILTONIAN SCHUR FORM DAVID S WATKINS Abstract Recently Chu Liu and Mehrmann developed an O(n 3 ) structure preserving method for computing the Hamiltonian
More informationStructured Condition Numbers of Symmetric Algebraic Riccati Equations
Proceedings of the 2 nd International Conference of Control Dynamic Systems and Robotics Ottawa Ontario Canada May 7-8 2015 Paper No. 183 Structured Condition Numbers of Symmetric Algebraic Riccati Equations
More informationHomework 2 Foundations of Computational Math 2 Spring 2019
Homework 2 Foundations of Computational Math 2 Spring 2019 Problem 2.1 (2.1.a) Suppose (v 1,λ 1 )and(v 2,λ 2 ) are eigenpairs for a matrix A C n n. Show that if λ 1 λ 2 then v 1 and v 2 are linearly independent.
More informationCharacterization of half-radial matrices
Characterization of half-radial matrices Iveta Hnětynková, Petr Tichý Faculty of Mathematics and Physics, Charles University, Sokolovská 83, Prague 8, Czech Republic Abstract Numerical radius r(a) is the
More informationThe matrix sign function
The matrix sign function 1 Re x > 0, sign(x) = 1 Re x < 0, undefined Re x = 0. Suppose the Jordan form of A is reblocked as A = [ V 1 V 2 ] [ J 1 J 2 ] [ V 1 V 2 ] 1, where J 1 contains all eigenvalues
More informationThe Cayley-Hamilton Theorem and the Jordan Decomposition
LECTURE 19 The Cayley-Hamilton Theorem and the Jordan Decomposition Let me begin by summarizing the main results of the last lecture Suppose T is a endomorphism of a vector space V Then T has a minimal
More informationQUALITATIVE CONTROLLABILITY AND UNCONTROLLABILITY BY A SINGLE ENTRY
QUALITATIVE CONTROLLABILITY AND UNCONTROLLABILITY BY A SINGLE ENTRY D.D. Olesky 1 Department of Computer Science University of Victoria Victoria, B.C. V8W 3P6 Michael Tsatsomeros Department of Mathematics
More informationPerturbation theory for eigenvalues of Hermitian pencils. Christian Mehl Institut für Mathematik TU Berlin, Germany. 9th Elgersburg Workshop
Perturbation theory for eigenvalues of Hermitian pencils Christian Mehl Institut für Mathematik TU Berlin, Germany 9th Elgersburg Workshop Elgersburg, March 3, 2014 joint work with Shreemayee Bora, Michael
More informationMath 102, Winter Final Exam Review. Chapter 1. Matrices and Gaussian Elimination
Math 0, Winter 07 Final Exam Review Chapter. Matrices and Gaussian Elimination { x + x =,. Different forms of a system of linear equations. Example: The x + 4x = 4. [ ] [ ] [ ] vector form (or the column
More informationSKEW-HAMILTONIAN AND HAMILTONIAN EIGENVALUE PROBLEMS: THEORY, ALGORITHMS AND APPLICATIONS
SKEW-HAMILTONIAN AND HAMILTONIAN EIGENVALUE PROBLEMS: THEORY, ALGORITHMS AND APPLICATIONS Peter Benner Technische Universität Chemnitz Fakultät für Mathematik benner@mathematik.tu-chemnitz.de Daniel Kressner
More informationMATH 425-Spring 2010 HOMEWORK ASSIGNMENTS
MATH 425-Spring 2010 HOMEWORK ASSIGNMENTS Instructor: Shmuel Friedland Department of Mathematics, Statistics and Computer Science email: friedlan@uic.edu Last update April 18, 2010 1 HOMEWORK ASSIGNMENT
More informationComponentwise perturbation analysis for matrix inversion or the solution of linear systems leads to the Bauer-Skeel condition number ([2], [13])
SIAM Review 4():02 2, 999 ILL-CONDITIONED MATRICES ARE COMPONENTWISE NEAR TO SINGULARITY SIEGFRIED M. RUMP Abstract. For a square matrix normed to, the normwise distance to singularity is well known to
More informationLinear Algebra and its Applications
Linear Algebra and its Applications 430 (2009) 579 586 Contents lists available at ScienceDirect Linear Algebra and its Applications journal homepage: www.elsevier.com/locate/laa Low rank perturbation
More informationTotal least squares. Gérard MEURANT. October, 2008
Total least squares Gérard MEURANT October, 2008 1 Introduction to total least squares 2 Approximation of the TLS secular equation 3 Numerical experiments Introduction to total least squares In least squares
More informationChapter 3 Transformations
Chapter 3 Transformations An Introduction to Optimization Spring, 2014 Wei-Ta Chu 1 Linear Transformations A function is called a linear transformation if 1. for every and 2. for every If we fix the bases
More informationA Continuation Approach to a Quadratic Matrix Equation
A Continuation Approach to a Quadratic Matrix Equation Nils Wagner nwagner@mecha.uni-stuttgart.de Institut A für Mechanik, Universität Stuttgart GAMM Workshop Applied and Numerical Linear Algebra September
More informationA q x k+q + A q 1 x k+q A 0 x k = 0 (1.1) where k = 0, 1, 2,..., N q, or equivalently. A(σ)x k = 0, k = 0, 1, 2,..., N q (1.
A SPECTRAL CHARACTERIZATION OF THE BEHAVIOR OF DISCRETE TIME AR-REPRESENTATIONS OVER A FINITE TIME INTERVAL E.N.Antoniou, A.I.G.Vardulakis, N.P.Karampetakis Aristotle University of Thessaloniki Faculty
More informationSome inequalities for sum and product of positive semide nite matrices
Linear Algebra and its Applications 293 (1999) 39±49 www.elsevier.com/locate/laa Some inequalities for sum and product of positive semide nite matrices Bo-Ying Wang a,1,2, Bo-Yan Xi a, Fuzhen Zhang b,
More informationAPPLIED NUMERICAL LINEAR ALGEBRA
APPLIED NUMERICAL LINEAR ALGEBRA James W. Demmel University of California Berkeley, California Society for Industrial and Applied Mathematics Philadelphia Contents Preface 1 Introduction 1 1.1 Basic Notation
More informationFINITE-DIMENSIONAL LINEAR ALGEBRA
DISCRETE MATHEMATICS AND ITS APPLICATIONS Series Editor KENNETH H ROSEN FINITE-DIMENSIONAL LINEAR ALGEBRA Mark S Gockenbach Michigan Technological University Houghton, USA CRC Press Taylor & Francis Croup
More informationAn Even Order Symmetric B Tensor is Positive Definite
An Even Order Symmetric B Tensor is Positive Definite Liqun Qi, Yisheng Song arxiv:1404.0452v4 [math.sp] 14 May 2014 October 17, 2018 Abstract It is easily checkable if a given tensor is a B tensor, or
More informationThe norms can also be characterized in terms of Riccati inequalities.
9 Analysis of stability and H norms Consider the causal, linear, time-invariant system ẋ(t = Ax(t + Bu(t y(t = Cx(t Denote the transfer function G(s := C (si A 1 B. Theorem 85 The following statements
More informationAffine iterations on nonnegative vectors
Affine iterations on nonnegative vectors V. Blondel L. Ninove P. Van Dooren CESAME Université catholique de Louvain Av. G. Lemaître 4 B-348 Louvain-la-Neuve Belgium Introduction In this paper we consider
More informationThe Nearest Doubly Stochastic Matrix to a Real Matrix with the same First Moment
he Nearest Doubly Stochastic Matrix to a Real Matrix with the same First Moment William Glunt 1, homas L. Hayden 2 and Robert Reams 2 1 Department of Mathematics and Computer Science, Austin Peay State
More informationMATH36001 Perron Frobenius Theory 2015
MATH361 Perron Frobenius Theory 215 In addition to saying something useful, the Perron Frobenius theory is elegant. It is a testament to the fact that beautiful mathematics eventually tends to be useful,
More informationBare-bones outline of eigenvalue theory and the Jordan canonical form
Bare-bones outline of eigenvalue theory and the Jordan canonical form April 3, 2007 N.B.: You should also consult the text/class notes for worked examples. Let F be a field, let V be a finite-dimensional
More informationOn rank one perturbations of Hamiltonian system with periodic coefficients
On rank one perturbations of Hamiltonian system with periodic coefficients MOUHAMADOU DOSSO Université FHB de Cocody-Abidjan UFR Maths-Info., BP 58 Abidjan, CÔTE D IVOIRE mouhamadou.dosso@univ-fhb.edu.ci
More informationAMS526: Numerical Analysis I (Numerical Linear Algebra for Computational and Data Sciences)
AMS526: Numerical Analysis I (Numerical Linear Algebra for Computational and Data Sciences) Lecture 1: Course Overview; Matrix Multiplication Xiangmin Jiao Stony Brook University Xiangmin Jiao Numerical
More informationAlgebra C Numerical Linear Algebra Sample Exam Problems
Algebra C Numerical Linear Algebra Sample Exam Problems Notation. Denote by V a finite-dimensional Hilbert space with inner product (, ) and corresponding norm. The abbreviation SPD is used for symmetric
More informationA Backward Stable Hyperbolic QR Factorization Method for Solving Indefinite Least Squares Problem
A Backward Stable Hyperbolic QR Factorization Method for Solving Indefinite Least Suares Problem Hongguo Xu Dedicated to Professor Erxiong Jiang on the occasion of his 7th birthday. Abstract We present
More informationExponentials of Symmetric Matrices through Tridiagonal Reductions
Exponentials of Symmetric Matrices through Tridiagonal Reductions Ya Yan Lu Department of Mathematics City University of Hong Kong Kowloon, Hong Kong Abstract A simple and efficient numerical algorithm
More informationMultiplicative Perturbation Bounds of the Group Inverse and Oblique Projection
Filomat 30: 06, 37 375 DOI 0.98/FIL67M Published by Faculty of Sciences Mathematics, University of Niš, Serbia Available at: http://www.pmf.ni.ac.rs/filomat Multiplicative Perturbation Bounds of the Group
More informationA Regularized Directional Derivative-Based Newton Method for Inverse Singular Value Problems
A Regularized Directional Derivative-Based Newton Method for Inverse Singular Value Problems Wei Ma Zheng-Jian Bai September 18, 2012 Abstract In this paper, we give a regularized directional derivative-based
More informationGeometric Mapping Properties of Semipositive Matrices
Geometric Mapping Properties of Semipositive Matrices M. J. Tsatsomeros Mathematics Department Washington State University Pullman, WA 99164 (tsat@wsu.edu) July 14, 2015 Abstract Semipositive matrices
More informationSolving Regularized Total Least Squares Problems
Solving Regularized Total Least Squares Problems Heinrich Voss voss@tu-harburg.de Hamburg University of Technology Institute of Numerical Simulation Joint work with Jörg Lampe TUHH Heinrich Voss Total
More informationIMPROVED METHODS AND STARTING VALUES TO SOLVE THE MATRIX EQUATIONS X ± A X 1 A = I ITERATIVELY
MATHEMATICS OF COMPUTATION Volume 74, Number 249, Pages 263 278 S 0025-5718(04)01636-9 Article electronically published on January 27, 2004 IMPROVED METHODS AND STARTING VAUES TO SOVE THE MATRIX EQUATIONS
More information