A new projection method for solving large Sylvester equations

被引:41
作者
Bao, Liang
Lin, Yiqin
Wei, Yimin [1 ]
机构
[1] Fudan Univ, Sch Math Sci, Inst Math, Shanghai 200433, Peoples R China
[2] Hunan Univ Sci & Engn, Dept Math & Computat Sci, Yongzhou 425006, Peoples R China
[3] Fudan Univ, Minist Educ, Key Lab Math Nonlinear Sci, Shanghai 200433, Peoples R China
基金
中国国家自然科学基金;
关键词
block Arnoldi; global Arnoldi; Sylvester equation; generalized sylvester equation; block Krylov subspace; matrix Krylov subspace;
D O I
10.1016/j.apnum.2006.07.005
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
In this paper, we propose a new projection method based on global Arnoldi algorithm for solving large Sylvester matrix equations A X + X B + C D-T = 0 and the large generalized Sylvester matrix equations of the form A X B + X + C D T = 0. We show how to extract low-rank approximate solutions to Sylvester matrix equations and generalized Sylvester matrix equations. Some theoretical results are given. Numerical tests report the effectiveness of these methods. (c) 2006 IMACS. Published by Elsevier B.V. All rights reserved.
引用
收藏
页码:521 / 532
页数:12
相关论文
共 20 条
[1]   Krylov subspace methods for the generalized Sylvester equation [J].
Bao, Liang ;
Lin, Yiqin ;
Wei, Yimin .
APPLIED MATHEMATICS AND COMPUTATION, 2006, 175 (01) :557-573
[2]   ALGORITHM - SOLUTION OF MATRIX EQUATION AX+XB = C [J].
BARTELS, RH ;
STEWART, GW .
COMMUNICATIONS OF THE ACM, 1972, 15 (09) :820-&
[3]   On the solution of large Sylvester-observer equations [J].
Calvetti, D ;
Lewis, B ;
Reichel, L .
NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS, 2001, 8 (6-7) :435-451
[4]  
Datta B. N., 1986, Computational and Combinatorial Methods in Systems Theory, P201
[5]  
Datta BN., 2004, Numerical Methods for Linear Control Systems
[6]  
Demmel JW., 1997, APPL NUMERICAL LINEA
[7]   CONTROLLABILITY, OBSERVABILITY AND THE SOLUTION OF AX-XB=C [J].
DESOUZA, E ;
BHATTACHARYYA, SP .
LINEAR ALGEBRA AND ITS APPLICATIONS, 1981, 39 (AUG) :167-188
[8]   Block Krylov subspace methods for solving large Sylvester equations [J].
El Guennouni, A ;
Jbilou, K ;
Riquet, AJ .
NUMERICAL ALGORITHMS, 2002, 29 (1-3) :75-96
[9]   HESSENBERG-SCHUR METHOD FOR THE PROBLEM AX+XB=C [J].
GOLUB, GH ;
NASH, S ;
VANLOAN, C .
IEEE TRANSACTIONS ON AUTOMATIC CONTROL, 1979, 24 (06) :909-913
[10]  
HEARON JZ, 1997, LINEAR ALGEBRA APPL, V16, P57