A LOOK-AHEAD LEVINSON ALGORITHM FOR INDEFINITE TOEPLITZ-SYSTEMS

被引:28
作者
CHAN, TF [1 ]
HANSEN, PC [1 ]
机构
[1] TECH UNIV DENMARK,DANISH COMP CTR RES & EDUC,DK-2800 LYNGBY,DENMARK
关键词
SYMMETRICAL INDEFINITE TOEPLITZ MATRICES; LEVINSON ALGORITHM; DURBIN ALGORITHM; NUMERICAL STABILITY; CONDITION ESTIMATION;
D O I
10.1137/0613032
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
An extension of Levinson's algorithm for solving linear systems with symmetric indefinite Toeplitz matrices is presented. This new algorithm is able to "look ahead" and, if necessary, use block Gaussian elimination to skip all ill-conditioned leading principal submatrices encountered during the recursive processing. This makes the new algorithm numerically stable for a broader class of symmetric Toeplitz matrices than the standard Levinson algorithm. In addition, a reliable condition number estimate is produced. The overhead is typically small.
引用
收藏
页码:490 / 506
页数:17
相关论文
共 22 条
[1]  
ARUSHANIAN OB, 1983, ANL8316 REP
[3]   STABILITY OF METHODS FOR SOLVING TOEPLITZ-SYSTEMS OF EQUATIONS [J].
BUNCH, JR .
SIAM JOURNAL ON SCIENTIFIC AND STATISTICAL COMPUTING, 1985, 6 (02) :349-364
[4]  
CHAN TF, IN PRESS IEEE T ACOU
[5]  
CHAN TF, 1990, CAM9011 U CAL DEP MA
[6]   COMPUTING THE MINIMUM EIGENVALUE OF A SYMMETRICAL POSITIVE DEFINITE TOEPLITZ MATRIX [J].
CYBENKO, G ;
VANLOAN, C .
SIAM JOURNAL ON SCIENTIFIC AND STATISTICAL COMPUTING, 1986, 7 (01) :123-131
[7]   THE NUMERICAL STABILITY OF THE LEVINSON-DURBIN ALGORITHM FOR TOEPLITZ-SYSTEMS OF EQUATIONS [J].
CYBENKO, G .
SIAM JOURNAL ON SCIENTIFIC AND STATISTICAL COMPUTING, 1980, 1 (03) :303-319
[8]   A GENERALIZATION OF THE LEVINSON ALGORITHM FOR HERMITIAN TOEPLITZ MATRICES WITH ANY RANK PROFILE [J].
DELSARTE, P ;
GENIN, YV ;
KAMP, YG .
IEEE TRANSACTIONS ON ACOUSTICS SPEECH AND SIGNAL PROCESSING, 1985, 33 (04) :964-971
[9]  
Golub G.H., 1996, MATH GAZ, VThird
[10]   INVERSION OF TOEPLITZ MATRICES WHICH ARE NOT STRONGLY NON-SINGULAR [J].
GOVER, MJC ;
BARNETT, S .
IMA JOURNAL OF NUMERICAL ANALYSIS, 1985, 5 (01) :101-110