SMALL-SAMPLE STATISTICAL CONDITION ESTIMATES FOR GENERAL MATRIX FUNCTIONS

被引:67
作者
KENNEY, CS
LAUB, AJ
机构
关键词
CONDITIONING; MATRIX FUNCTIONS;
D O I
10.1137/0915003
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
A new condition estimation procedure for general matrix functions is presented that accurately gauges sensitivity by measuring the effect of random perturbations at the point of evaluation. In this procedure the number of extra function evaluations used to evaluate the condition estimate determines the order of the estimate. That is, the probability that the estimate is off by a given factor is inversely proportional to the factor raised to the order of the method. The ''transpose-free'' nature of this new method allows it to be applied to a broad range of problems in which the function maps between spaces of different dimensions. This is in sharp contrast to the more common power method condition estimation procedure that is limited, in the usual case where the Frechet derivative is known only implicitly, to maps between spaces of equal dimension. A group of examples illustrates the flexibility of the new estimation procedure in handling a variety of problems and types of sensitivity estimates, such as mixed and componentwise condition estimates.
引用
收藏
页码:36 / 61
页数:26
相关论文
共 49 条
[1]  
Abramowitz M., 1965, HDB MATH FUNCTIONS
[2]  
[Anonymous], 1968, INTRO PROBABILITY TH
[3]  
[Anonymous], 1979, MONTE CARLO METHODS
[4]   SOLVING SPARSE LINEAR-SYSTEMS WITH SPARSE BACKWARD ERROR [J].
ARIOLI, M ;
DEMMEL, JW ;
DUFF, IS .
SIAM JOURNAL ON MATRIX ANALYSIS AND APPLICATIONS, 1989, 10 (02) :165-190
[5]  
ATKINSON K. E., 1978, INTRO NUMERICAL ANAL
[7]  
BYERS R, 1985, LINEAR ALGEBRA ITS R, V47, P35
[8]   ESTIMATE FOR THE CONDITION NUMBER OF A MATRIX [J].
CLINE, AK ;
MOLER, CB ;
STEWART, GW ;
WILKINSON, JH .
SIAM JOURNAL ON NUMERICAL ANALYSIS, 1979, 16 (02) :368-375
[9]  
Copson E. T., 1965, ASYMPTOTIC EXPANSION
[10]   ACCURATE SINGULAR-VALUES OF BIDIAGONAL MATRICES [J].
DEMMEL, J ;
KAHAN, W .
SIAM JOURNAL ON SCIENTIFIC AND STATISTICAL COMPUTING, 1990, 11 (05) :873-912