上海大学学报(自然科学版) ›› 2013, Vol. 19 ›› Issue (3): 303-307.doi: 10.3969/j.issn.1007-2861.2013.03.016

• 数理化科学 • 上一篇    下一篇

二元齐次矩阵Padé-型逼近的计算

潘宝珍, 刘永, 潘鹿鹿   

  1. 上海大学 理学院, 上海 200444
  • 收稿日期:2012-12-26 出版日期:2013-06-30 发布日期:2013-06-30
  • 通讯作者: 潘宝珍(1965—), 女, 副教授, 博士, 研究方向为数值有理逼近. E-mail:bzpan@staff.shu.edu.cn
  • 基金资助:

    上海市重点学科建设资助项目(S30104)

Computation of Bivariate Homogeneous Matrix Padé-Type Approximation

PAN Bao-zhen, LIU Yong, PAN Lu-lu   

  1. College of Sciences, Shanghai University, Shanghai 200444, China
  • Received:2012-12-26 Online:2013-06-30 Published:2013-06-30

摘要: 二元齐次矩阵Padé-型逼近的计算比较复杂, 而通过适当的变量代换, 可以将二元齐次矩阵形式幂级数转化为一元含参数形式的矩阵形式幂级数, 从而给出二元齐次矩阵Padé-型逼近构造性的定义. 为提高二元齐次矩阵Padé-型逼近的逼近解精度, 借助于误差公式推导出基于矩阵EMN 的二元齐次矩阵正交多项式Padé-型逼近的分子和分母行列式表达式; 为避免计算高阶行列式, 建立了一种Sylvester-型递推算法. 最后, 通过数值算例验证了该算法的有效性.

关键词: 二元齐次, Padé-型逼近, 递推算法, 矩阵形式幂级数, 正交多项式

Abstract: With appropriate variable replacement, the bivariate homogeneous matrix formal power series is transformed to univariate matrix formal power series with parameters. The bivariate homogeneous matrix Padé-type approximation was defined. To improve computation accuracy, using an error formula, the numerator and denominator in the determinant expressions of bivariate homogeneous matrix orthogonal polynomial Padé-type approximation are given based on the matrix EMN. A Sylvester-type recursive algorithm is presented to avoid computation of high degree determinants. A numerical example shows effectiveness of the algorithm.

Key words: bivariate homogeneous, iterative algorithm, matrix formal power series, orthogonal polynomial, Padé-type approximation

中图分类号: