矩阵指数精细积分方法中参数的自适应选择  

ADAPTIVE SELECTION OF PARAMETERS FOR PRECISE COMPUTATION OF MATRIX EXPONENTIAL BASED ON PADE APPROXIMATION

作　　者：谭述君[1] 吴志刚[1] 钟万勰[1] 

机构地区：[1]大连理工大学工业装备结构分析国家重点实验室,大连116023

出　　处：《力学学报》2009年第6期961-966,共6页Chinese Journal of Theoretical and Applied Mechanics

基　　金：国家自然科学基金(10632030);国家重点基础研究专项经费(2005CB321704);高等学校博士学科点专项科研基金(20070141067)资助项目~~

摘　　要：讨论了基于Pade逼近的矩阵指数精细积分方法中加权系数N和展开项数q的自适应选择问题。参数(N,q)的选择直接影响到矩阵指数计算的精度和效率。采用矩阵函数逼近理论,研究了参数N和q的增加对精度的影响程度,据此,提出了参数(N,q)优化组合的递推自适应选择方法。该方法可以根据矩阵本身的性态选择合适的参数(N,q),而参数选择的计算量与矩阵指数的计算量相比几乎可以忽略,这对于增强矩阵指数精细积分方法的适应性和提高计算效率是很有益处的。算例验证了该方法的正确性和有效性。Adaptive selection is discussed for scaling parameters N and expanded series q in precise integration method （PIM） of matrix exponential based on Pad6 approximation. In general, scaling parameters N and expanded series q play important roles in the numerical accuracy and computational efficiency of matrix exponential. Using the approximation theory of matrix functions, influences of parameters N and q on the computational accuracy and efficiency are firstly studied, and then the iterative adaptive selection algorithm for optimal combination of parameters （N, q） is presented. Appropriate parameters （N, q） can be selected automat- ically depending on the characteristics of the matrix, and the computation amount of adaptive selection can be neglected compared with that of matrix exponential. So it is very important for enhancing the adaptations and increasing the computation efficiencies of matrix exponential. In addition, computational examples are carried out to testify the correctness and validity of the present method.

关 键 词：矩阵指数 Padé级数逼近 精细积分方法 加权平方 误差分析 

分 类 号：O321[理学—一般力学与力学基础]