pade逼近.docVIP

5．2 pade逼近方法的简介 5.2.1 泰勒级数问题 一个函数的泰勒级数展开式的系数，和这个函数的值之间的关系，既是一个深奥的数学问题，又是一个重要的实际问题。关于它的研究，是基于数学分析和物理、生物科学中的自然数学模型的实际计算的基础之上的。关于这个问题的某些方面已经有过一些研究，但是还有一些问题需要在今后继续研究。规范的解释是：如果泰勒级数展开式绝对收敛，那么它唯一确定了一个任意次可微的函数。相反的，如果一个函数是任意次可微的，那它也只有一个对应的泰勒级数展开式。实际上，我们把函数近似为一个尽可能长的多项式。然而，这种方法在实际运算中有一些很不理想的局限性[9]。考虑以下例子： （5.18） 很容易看出，当x0.5时，泰勒级数表达式均不收敛，尽管在时，是一个保持在的光滑适当的函数。规范方法是要在一个新点（），应用原来的表达式计算，作出一个新的泰勒级数表达式。这个新的表达式可以满足在x很大时的情况，但是不包括处。实际上，用这种方法我们永远不能达到，并且在这个方向上的任何进展都是非常冗长的。对于以上这个例子，我们可以应用一个特殊的技巧，将级数转变为一个较长的多项式。假设我们做一个变量代换： x=w/(1-2w) 或者 w=x/(1+2x) （5.19） 就有： （5.20） 所以，在这个变量代换下，就变为w=0.5。很容易看出，泰勒级数表达式（5.20）收敛于w=0.5，即。于是的最初连续估计值是： 1，1.25，1.34375，1.38281，1.39990…… （5.21） 可以看出它收敛于=1.414…。根据最初的变量x，（5.20）可以表示为： （5.22） 均是关于x的有理分式。 Pade逼近 Pade逼近是一种关于函数值的特殊类型的有理分式逼近法。它的思想是以尽量快的速度与泰勒级数展开式相匹配[10]。举例来说，对于例（5.18），我们选择一个逼近模式： （a+bx）/(c+dx) （5.23） 可见当x趋于无穷时，上式是有界的。如果我们应用（5.18）式的前三项系数做逼近，就可以得到以下的表达式： （5.24） 对于此式，当时，它的值是1.4，这个结果比（5.21）式中的任意一个结果都要理想。同理可算出下一个逼近结果是： （5.25） 这个结果与已经非常接近了。继续以此方法作逼近，随着所选系数数量的增加，结果的收敛性很好，结果如下： 1.414201183, 1.414213198, 1.414213552, … （5.26） 其中的最后一个结果应用了（5.18）式中的前11项系数，且与标准结果只相差在处。 我们可以用同样的方法对（5.20）式做逼近，即将式（5.19）带入以上所求结果，有： 1， ， （5.27） 以上三个结果分别是应用（5.20）式的第一项、前三项、前五项得出的。可以注意到，当我们将w=0.5(即对应于)带入（5.27）式中，就得到结果1，1.4，41/29，……。这些结果与我们计算关于x的级数展开式时所得到的结果对应一致。这种同一性是Pade逼近的一个普遍而且重要的特性，也正因为这种特性，可以在我们的例子中做出关于x的逼近，并且甚至在处都可以得到很理想的结果。 我们还可以观察到，（5.25）和（5.26）式的连续逼近值是单调的。虽然这不是一个普遍的特性，但是在多种情况下都可以证明它的成立。在实际问题中，它是一个非常重要的性质，通过它通常可以证明，Pade逼近在高范围和低范围内都收敛。 现在我们将给Pade逼近和Pade表格做出定义： 定义：我们将关于A(x)的Pade逼近记为： （5.28） 其中，是一个次数最高为L的多项式，是一个次数最高为M的多项式。 由幂级数： （5.29） 通过公式（5.30）可以确定多项式和的系数。 （5.30） 由于我们可以增加分子和分母的项数并且保持不变，于是定义标准化条件为： （5.31） 最后需要注意，和必须是没有公因子的多项式。 如果我们将和的