二元齐次矩阵Padé-型逼近的计算比较复杂, 而通过适当的变量代换, 可以将二元齐次矩阵形式幂级数转化为一元含参数形式的矩阵形式幂级数, 从而给出二元齐次矩阵Padé-型逼近构造性的定义. 为提高二元齐次矩阵Padé-型逼近的逼近解精度, 借助于误差公式推导出基于矩阵EMN 的二元齐次矩阵正交多项式Padé-型逼近的分子和分母行列式表达式; 为避免计算高阶行列式, 建立了一种Sylvester-型递推算法. 最后, 通过数值算例验证了该算法的有效性.
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.
[1] Brezinski C. Padé-type approximation and general orthogonal polynomials [M]. Basel: Birkhauser, 1980.
[2] Draux A. Approximants de type Padé et de Pade [M]. Lille: Universié des Science et Technologies de Lille, 1983: 1-89.
[3] Gu C Q. Matrix Padé-type approximant and directional matrix Padé-type approximant in the inner product space [J]. J Comput Appl Math, 2004, 164/165: 365-385.
[4] Zheng C D. Generalized homogeneous multivariate matrix Padé-type approximants and Padé-type approximants [J]. IEEE Transactions on Automatic Control, 2007, 52(11): 2160-2165.
[5] Tao Y T, Gu C Q. A two-dimensional matrix Padé type approximant in the inner product space [J]. J Comput Appl Math, 2009, 231(2): 680-695.
[6] Benouahmane B, Cuyt A. Multivariate orthogonal polynomials, homogeneous Padé approximants and Gaussian cubature [J]. J Numerical Algorithms, 2000, 24(1/2): 1-15.
[7] 潘宝珍. 二元齐次矩阵Padé-型逼近及误差公式[J]. 应用数学与计算数学学报, 2012, 26(1): 113-120.
[8] Baker G, Graves-Morris P. Padé approximants [M]. 2nd ed. New York: Cambridge University Press, 1997.
