跳转到内容

2.5 Taylor定理及其应用

王绵森、马知恩主编高等教育出版社《工科数学分析基础(第三版)》大学本科 · 教材双语助读 (Alt+Y)

用已知点的信息来表达未知点信息,用简单函数逼近(近似表示)复杂函数是数学中的重要思想方法。本节将要介绍的 Taylor 定理就是用高阶多项式来逼近具有足够可微性函数所得到的一个基本定理,它在理论研究和近似计算中有重要的应用。

第三节中曾经指出,如果函数 f 在 x₀ 处可微,那么就可用微分来近似计算函数 f 在 x₀ 附近的值。也就是说,当 |x - x₀| 很小时,可用线性函数(一次多项式)来近似表示 f(x),即

f(x) ≈ f(x₀) + f'(x₀)(x - x₀).

这个近似公式具有形式简单、计算方便的优点,但也存在着精度不高、计算误差仅是 x - x₀ 一次幂的高阶无穷小 o(x - x₀) 的缺陷。之所以有这样的缺陷,从几何上看,是由由于这个近似公式是用曲线 y = f(x) 上点 (x₀, f(x₀)) 处的切线(直线)来代替该曲线得到的。为了改进该公式,我们自然想到用曲线来代替曲线。在曲线中,比较简单的是关于 x - x₀ 的高次多项式所表示的曲线,因为多项式也具有形式简单,计算方便(只用到加法和乘法运算)的优点。

现在的问题是:能否找到一个适当的 n (n > 1) 次多项式

Pₙ(x) = a₀ + a₁(x - x₀) + a₂(x - x₀)² + … + aₙ(x - x₀)ⁿ

来逼近 f(x),并使误差为 (x - x₀)ⁿ 的高阶无穷小?如果能够找到,那么使等式

f(x) = Pₙ(x) + o((x - x₀)ⁿ) (5.1)

成立的 Pₙ(x) 与 f(x) 应满足什么条件?上一节例4.10 表明,对于指数函数 f(x) = eˣ,(5.1) 式成立(其中 x₀ = 0),并且 Pₙ(x) 的系数

aₖ = (1)/(k!) f^((k))(0) (k = 0, 1, 2, …, n).

这启发我们,对更一般的函数 f,如果在 x₀ 处 n 阶可微,并且 Pₙ(x) 的系数

aₖ = (1)/(k!) f^((k))(x₀) (k = 0, 1, 2, …, n),

那么 (5.1) 式是否也能成立呢?下面的定理对这个问题作了肯定的回答。

定理 5.1 (带 Peano 余项的 Taylor 定理)

设函数 f 在 x₀ 处 n 阶可微,则

f(x) = ∑ₖ=0ⁿ fracf^((k))(x₀)k! (x - x₀)^k + o((x - x₀)ⁿ). (5.2)
注意

舍弃余项 o((x - x₀)ⁿ),(5.2) 式给出了用 n 次多项式(简单函数)逼近 n 阶可微函数(复杂函数)f(x) 的一个基本定理,本质上是用 f(x) 及其各阶导数在一点 x₀ 处的值(信息)来估算它在 x₀ 的一个小邻域 U(x₀) 内任一点处值(信息)的一个公式,体现了数学中的一种重要的科学思维方法。

证明

为了证明 (5.2) 式成立,我们令

Pₙ(x) = ∑ₖ=0ⁿ fracf^((k))(x₀)k! (x - x₀)^k, (5.3)Rₙ(x) = f(x) - Pₙ(x) (称为 f(x) 的余项),

那么,只要证明 Rₙ(x) = o((x - x₀)ⁿ) 也就是 limₓ → x₀ (Rₙ(x))/((x - x₀)ⁿ) = 0 就行了。由于 f 在 x₀ 处 n 阶可微,所以 f^((n-1)) 在 x₀ 处可微。根据定理1.1 可知 f^((n-1)) 必在 x₀ 处连续,再由连续的定义可知,f 在 x₀ 的某邻域内 n-1 阶可微。因此,Rₙ(x) 在此邻域内也 n-1 阶可微,并且

{ R'ₙ(x) &= f'(x) - ∑ₖ=1ⁿ fracf^((k))(x₀)(k - 1)! (x - x₀)^(k-1), R''ₙ(x) &= f''(x) - ∑ₖ=2ⁿ fracf^((k))(x₀)(k - 2)! (x - x₀)^(k-2), & vdots Rₙ^((n-1))(x) &= f^((n-1))(x) - [f^((n-1))(x₀) + f^((n))(x₀)(x - x₀)]. . (5.4)

故

limₓ → x₀ Rₙ(x) = limₓ → x₀ R'ₙ(x) = … = limₓ → x₀ Rₙ^((n-1))(x) = 0.

下面证明 limₓ → x₀ (Rₙ(x))/((x - x₀)ⁿ) = 0。对此式左端连续使用 n-1 次 L’Hospital 法则可得

limₓ → x₀ (Rₙ(x))/((x - x₀)ⁿ) = limₓ → x₀ fracR'ₙ(x)n(x - x₀)ⁿ-1 = … = (1)/(n!) limₓ → x₀ fracRₙ^((n-1))(x)x - x₀.
想一想

为什么说:“故上式中右端的不定式不能再用 L’Hospital 法则”呢?

由于定理中仅假设 f 在 x₀ 处 n 阶可微,因而 Rₙ(x) 也仅在 x₀ 处 n 阶可微,故上式中右端的不定式不能再用 L’Hospital 法则,为了证明其值为 0,需要另想办法。由 (5.4) 式中最后一个等式及导数的定义,易得

limₓ → x₀ fracRₙ^((n-1))(x)x - x₀ = limₓ → x₀ [ fracf^((n-1))(x) - f^((n-1))(x₀)x - x₀ - f^((n))(x₀) ] = 0,

从而 limₓ → x₀ (Rₙ(x))/((x - x₀)ⁿ) = 0,即 Rₙ(x) = o((x - x₀)ⁿ)。

称由 (5.3) 式所表示的多项式为函数 f 在 x₀ 处的 n 次 Taylor 多项式,其系数称为 f 在 x₀ 处的 Taylor 系数。称公式 (5.2) 为 f 在 x₀ 处的带 Peano 余项的 Taylor 公式,而 Rₙ(x) = o((x - x₀)ⁿ) 称为 f 的 Peano 余项,它就是用 n 次 Taylor 多项式 Pₙ(x) 来近似计算 f(x) 所产生的误差。然而,还应注意,用这个余项来估计误差,仅能说明误差是 (x - x₀)ⁿ 的高阶无穷小,还不能对误差的大小作具体的数值分析,也不能说明用 Pₙ(x) 来计算 f(x) 的近似值达到怎样的精确度。因而,还需要研究余项 Rₙ(x) 的更具体的表达式。为此,下面给出 Taylor 定理的另一表达形式。

定理 5.2 (带 Lagrange 余项的 Taylor 定理)

设函数 f 在区间 I 上 n+1 阶可导,x₀ ∈ I,则对任何 x ∈ I,在 x 与 x₀ 之间(不含 x₀ 与 x)至少存在一点 ξ,使

f(x) = ∑ₖ=0ⁿ fracf^((k))(x₀)k! (x - x₀)^k + fracf^((n+1))(ξ)(n + 1)! (x - x₀)ⁿ+1. (5.5)
证明

为了证明 (5.5) 式成立,只要证明

Rₙ(x) = f(x) - Pₙ(x) = fracf^((n+1))(ξ)(n + 1)! (x - x₀)ⁿ+1,

或

fracRₙ(x)(x - x₀)ⁿ+1 = fracf^((n+1))(ξ)(n + 1)!.

令 g(x) = (x - x₀)ⁿ+1,则

g(x₀) = g'(x₀) = … = g^((n))(x₀) = 0, g^((n+1))(x₀) = (n + 1)!.

由余项 Rₙ(x) 的表达式不难验证

Rₙ(x₀) = R'ₙ(x₀) = … = Rₙ^((n))(x₀) = 0, Rₙ^((n+1))(x) = f^((n+1))(x).

在以 x 与 x₀ 为端点的区间上对 Rₙ(x) 与 g(x) 应用 Cauchy 定理,得

(Rₙ(x))/(g(x)) = (Rₙ(x) - Rₙ(x₀))/(g(x) - g(x₀)) = (R'ₙ(ξ₁))/(g'(ξ₁)) (ξ₁ 介于 x₀ 与 x 之间);

在以 x₀ 与 ξ₁ 为端点的区间上对 R'ₙ(x) 与 g'(x) 应用 Cauchy 定理,又得

(Rₙ(x))/(g(x)) = (R'ₙ(ξ₁) - R'ₙ(x₀))/(g'(ξ₁) - g'(x₀)) = (R''ₙ(ξ₂))/(g''(ξ₂)) (ξ₂ 介于 x₀ 与 ξ₁ 之间);

上述步骤继续进行 n+1 次,不难得知

(Rₙ(x))/(g(x)) &= (R'ₙ(ξ₁))/(g'(ξ₁)) = (R''ₙ(ξ₂))/(g''(ξ₂)) = … = fracRₙ^((n))(ξₙ)g^((n))(ξₙ) = fracRₙ^((n))(ξₙ) - Rₙ^((n))(x₀)g^((n))(ξₙ) - g^((n))(x₀) &= fracRₙ^((n+1))(ξ)g^((n+1))(ξ) = fracf^((n+1))(ξ)(n + 1)!,

其中 ξ 介于 x₀ 与 ξₙ 之间,ξₙ 介于 x₀ 与 ξₙ-1 之间,…,ξ₁ 介于 x₀ 与 x 之间,因而 ξ 也介于 x₀ 与 x(不含 x₀ 与 x)之间,从而定理得证。

(5.5) 式称为 f 在 x₀ 的邻域内带 Lagrange 余项的 n 阶 Taylor 公式,其中

Rₙ(x) = fracf^((n+1))(ξ)(n + 1)! (x - x₀)ⁿ+1

称为 Lagrange 余项,也可写成如下形式:

Rₙ(x) = fracf^((n+1))[x₀ + θ(x - x₀)](n + 1)! (x - x₀)ⁿ+1, θ ∈ (0, 1).

在 (5.5) 式中取 n = 0,就得到 Lagrange 公式。因此,带 Lagrange 余项的 Taylor 公式是 Lagrange 公式的推广。

Lagrange 余项较之于 Peano 余项更便于对误差进行数值估计。事实上,如果 f 在区间 I 上 n+1 阶可导,并且存在常数 M > 0,使得 ∀ x ∈ I, |f^((n+1))(x)| ≤slant M,那么

|Rₙ(x)| = | fracf^((n+1))(ξ)(n + 1)! | |x - x₀|ⁿ+1 ≤slant (M)/((n + 1)!) |x - x₀|ⁿ+1.

又若 I 的左、右端点分别是有限值 a 与 b,则

|Rₙ(x)| ≤slant (M)/((n + 1)!) (b - a)ⁿ+1 → 0 (n → ∞).

这表明,只要 f^((n+1))(x) 在有限区间 I 上有界,用 n 次 Taylor 多项式 Pₙ(x) 来近似代替 f(x),其绝对误差 |Rₙ(x)| 随着 n 的增大可变得任意小,从而可以选取适当的 n 使计算达到要求的任何精确度。

5.2 几个初等函数的 Maclaurin 公式

Section titled “5.2 几个初等函数的 Maclaurin 公式”

如果 x₀ = 0,那么公式 (5.5) 就变成

f(x) = f(0) + f'(0)x + (f''(0))/(2!)x² + … + fracf^((n))(0)n!xⁿ + fracf^((n+1))(θ x)(n + 1)!xⁿ+1, θ ∈ (0, 1), (5.6)

它常被称为 f 的 Maclaurin 公式,是 Taylor 公式的一种常用的特殊情形。

下面,利用 (5.6) 式求出几个常用初等函数的 Maclaurin 公式。

指数函数 f(x) = eˣ 的 Maclaurin 公式 根据例2.14 中的 n 阶导数公式(1),f^((k))(x) = eˣ (k = 0, 1, 2, …, n),从而有 f^((k))(0) = 1,代入 (5.6) 式便得

eˣ = 1 + x + (x²)/(2!) + … + (xⁿ)/(n!) + fracxⁿ+1(n + 1)! e^(θ x), x ∈ (-∞, +∞), θ ∈ (0, 1). (5.7)

正弦函数 f(x) = sin x 的 Maclaurin 公式 根据例2.14 中的 n 阶导数公式(2),

f^((k))(x) = sin(x + k · (π)/(2)), k = 0, 1, 2, …, n.

从而有:当 k = 2m 时,f^((2m))(0) = 0;当 k = 2m + 1 时,f^((2m+1))(0) = (-1)^m,代入公式 (5.6) 得

sin x = x - (x³)/(3!) + (x⁵)/(5!) - (x⁷)/(7!) + … + (-1)^(m-1) fracx²m-1(2m - 1)! + (-1)^m (cos θ x)/((2m + 1)!) x²m+1, x ∈ (-∞, +∞), θ ∈ (0, 1). (5.8)

类似可得余弦函数 f(x) = cos x 的 Maclaurin 公式

cos x = 1 - (x²)/(2!) + (x⁴)/(4!) - (x⁶)/(6!) + … + (-1)^m fracx²m(2m)! + (-1)^(m+1) (cos θ x)/((2m + 2)!) x²m+2, x ∈ (-∞, +∞), θ ∈ (0, 1). (5.9)

对数函数 f(x) = ln(1 + x) 的 Maclaurin 公式 根据例2.14 中的 n 阶导数公式(5),

f^((k))(x) = (-1)^(k-1) ((k - 1)!)/((1 + x)^k), k = 1, 2, …, n.

从而有 f^((k))(0) = (-1)^(k-1)(k - 1)!,代入公式 (5.6) 便得

ln(1 + x) = x - (x²)/(2) + (x³)/(3) - (x⁴)/(4) + … + (-1)ⁿ-1 (xⁿ)/(n) + (-1)ⁿ fracxⁿ+1(n + 1)(1 + θ x)ⁿ+1, x ∈ (-1, +∞), θ ∈ (0, 1). (5.10)

幂函数 f(x) = (1 + x)^α (α ∈ R) 的 Maclaurin 公式 根据例2.14 中的 n 阶导数公式(4),

f^((k))(x) = α(α - 1)…(α - k + 1)(1 + x)^(α - k), k = 0, 1, 2, …, n.

从而有 f^((k))(0) = α(α - 1)…(α - k + 1),代入公式 (5.6) 便得

(1 + x)^α = 1 + α x + (α(α - 1))/(2!) x² + … + (α(α - 1)…(α - n + 1))/(n!) xⁿ + (α(α - 1)…(α - n))/((n + 1)!) fracxⁿ+1(1 + θ x)ⁿ+1-α, x ∈ (-1, +∞), θ ∈ (0, 1). (5.11)

取 α = -1, 1/2, -1/2 就得到三个常用幂函数的 Maclaurin 公式:

(1)/(1 + x) = 1 - x + x² - x³ + … + (-1)ⁿ xⁿ + (-1)ⁿ+1 fracxⁿ+1(1 + θ x)ⁿ+2, √(1 + x) = 1 + 1/2x - 1/8x² + … + (-1)ⁿ-1 (1 · 3 …… (2n - 3))/(2 · 4 …… 2n) xⁿ + (-1)ⁿ (1 · 3 …… (2n - 1))/(2 · 4 …… (2n + 2)) fracxⁿ+1(1 + θ x)ⁿ + 1/2, (1)/(√(1 + x)) = 1 - 1/2x + 3/8x² - … + (-1)ⁿ (1 · 3 …… (2n - 1))/(2 · 4 …… 2n) xⁿ + (-1)ⁿ+1 (1 · 3 …… (2n + 1))/(2 · 4 …… (2n + 2)) fracxⁿ+1(1 + θ x)ⁿ + 3/2.

仿照上面的方法,读者不难写出上述几个初等函数带 Peano 余项的 Maclaurin 公式。

Taylor 公式有重要的理论意义和应用价值。下面仅举例说明它在近似计算函数值、求极限和证明不等式方面的应用。它在研究函数性态(极值的充分条件)方面的应用在第六节中再作说明。

在近似计算函数值的时候,利用带 Lagrange 余项的 Taylor 公式较之利用微分精确度更高,适用的范围更广,而且可以估计误差。如果当 n → ∞ 时,Rₙ(x) → 0,那么可以把函数值计算到任何精度。

例 5.1

近似计算 e 的值,并估计误差。

解

在 eˣ 的 Taylor 公式 (5.7) 中取 x = 1,得

e = 1 + 1 + (1)/(2!) + … + (1)/(n!) + (e^θ)/((n + 1)!), θ ∈ (0, 1).

由于 e^θ < e < 3,故

Rₙ(1) = (e^θ)/((n + 1)!) < (3)/((n + 1)!) → 0 (n → ∞).

因此,只要 n 取得充分大,用近似公式

e ≈ 2 + (1)/(2!) + … + (1)/(n!)

来计算 e 的近似值可以达到所需要的任何精度。例如,要使误差小于 10⁻5,即要

Rₙ(1) < (3)/((n + 1)!) < 10⁻5,

只要取 n = 8 就行了,因为 R₈(1) < (3)/(9!) < 10⁻5。于是 e 的误差小于 10⁻5 的近似值为

e ≈ 2 + (1)/(2!) + … + (1)/(8!) ≈ 2.71828.
数字资源 2.5.2Taylor 定理的应用
*例 5.2

求方程

x⁵ + ε x - 32 = 0

的一个近似实根,其中 ε 是一个很小的实参数。

解

根据本书习题1.5(A) 第 14 题,此方程至少有一个依赖于 ε 的实根 x = f(ε)。由于五次方程没有精确的求根公式,所以只能求根的近似表达式。大家知道,数学中的方程是实际问题的近似抽象。但在实际问题中,对参数 ε 的测量总是有误差的,因而方程根也会因之产生误差。通常人们总希望参数的微小变化所引起的方程根的变化也很微小,否则,很难利用该方程对问题进行有效的研究。假定由该方程所确定的隐函数 x = f(ε)(即该方程的根)是连续的,而且具有足够的可微性。利用 Taylor 公式 (5.6) 该方程的根就近似表示为(仅取 n = 2)

x = f(ε) ≈ f(0) + f'(0)ε + (f''(0))/(2!) ε². (5.12)

为确定上式右端的诸系数,在已知方程中取 ε = 0,得 x = 2,故 f(0) = 2。为求 f'(0) 与 f''(0),利用隐函数求导法对原方程两端关于 ε 求导得

5x⁴ (dx)/(dε) + ε (dx)/(dε) + x = 0, (5.13)

从而有

f'(0) = . (dx)/(dε) |_(ε = 0) = -. (x)/(5x⁴ + ε) |_(ε = 0) = -(1)/(40).

再对方程 (5.13) 关于 ε 求导得

(5x⁴ + ε) (d² x)/(dε²) + 20x³ ( (dx)/(dε) )² + 2(dx)/(dε) = 0,

从而解得

f''(0) = . (d² x)/(dε²) |_(ε = 0) = -. (20x³ ( (dx)/(dε) )² + 2(dx)/(dε))/(5x⁴ + ε) |_(ε = 0) = -(1)/(1600).

于是得到所求方程实根的二阶近似表达式

x ≈ 2 - (ε)/(40) - (ε²)/(3200).

根据实际问题的精度要求,还可以求出该方程更高阶的近似根的表达式。这种利用 Taylor 公式求方程的关于小参数 ε 的近似根的方法就是所谓摄动法。

带 Peano 余项的 Taylor 公式也是求函数极限的一种常用方法。

例 5.3

求极限 limₓ → 0 (sin x - x)/(x³)。

解
原式 = limₓ → 0 (x - (x³)/(3!) + o(x⁴) - x)/(x³) = limₓ → 0 [ -(1)/(3!) + (o(x⁴))/(x³) ] = -1/6.
想一想

在例5.3 与例5.4 中,对分子利用 Taylor 公式时,是根据什么确定 Taylor 公式的阶数的?

例 5.4

求 limₓ → 0 fraccos x - e^((x²)/(2))x⁴。

解
原式 &= limₓ → 0 ([ 1 - (x²)/(2!) + (x⁴)/(4!) + o₁(x⁵) ] - [ 1 - (x²)/(2) + (1)/(2!) (x⁴)/(4) + o₂(x⁴) ])/(x⁴) &= limₓ → 0 [ -(1)/(12) + (o₁(x⁵))/(x⁴) - (o₂(x⁴))/(x⁴) ] = -(1)/(12).
例 5.5

设 f''(x) > 0,当 x → 0 时,f(x) 与 x 是等价无穷小。证明:当 x ≠ 0 时,f(x) > x。

证

由于题中函数二阶可导,所以利用二阶 Taylor 公式来证明可能更方便(若用 Lagrange 定理,可能要用两次)。取 x₀ = 0,由 Taylor 公式 (5.5) 得

f(x) = f(0) + f'(0)x + (f''(ξ))/(2!) x² (ξ 在 x 与 0 之间).

又已知当 x → 0 时,f(x) 与 x 是等价无穷小,故又有 f(x) = x + o(x)。将 f(x) 两个表达式比较可得 f(0) = 0, f'(0) = 1,所以

f(x) = x + (f''(ξ))/(2) x².

又 f''(x) > 0,因此当 x ≠ 0 时有 f(x) > x。

用同样的方法可证,若此例中 f''(x) < 0,则有不等式 f(x) < x。

由例5.5 及上面的结论容易得到下列不等式:

eˣ > 1 + x (x ≠ 0), sin x < x (0 < x ≤slant π); ln(1 + x) < x (x > -1, x ≠ 0), arcsin x > x (-1 < x < 1, x ≠ 0) 等。
共 0 题