高等微积分 Lesson 21
Taylor 公式的应用
设f在I处处有n阶导,则
f(b)=i=0∑n−1i!f(i)(a)(b−a)i+Rn
其中Rn可以表示为
Lagrange 余项,∃ξ介于a,b之间使得
Rn=n!f(n)(ξ)(b−a)n
Cauchy 余项,∃ξ介于a,b之间使得
Rn=(n−1)!f(n)(ξ)(b−ξ)n−1(b−a)
记忆方法:对于多项式P,degP≤n−1,可以推出
P(b)=i=0∑n−1i!P(i)(a)(b−a)i
此亦为 Taylor 公式,对一般的f(x),可以定义
F(x)=i=0∑n−1i!f(i)(x)(b−x)i
对F和任意的一个g(x),可以用 Cauchy 中值定理,产生
g(b)−g(a)F(b)−F(a)⟹f(b)−i=0∑n−1i!f(i)(a)(b−a)i=g′(ξ)F′(ξ)=g′(ξ)g(b)−g(a)(n−1)!f(n)(ξ)(b−ξ)n−1
得到任何一种余项. 对于 Cauchy 余项,就是取g(x)=x.
为了更好的应用 Taylor 公式,我们从不同角度解读公式的含义,有以下几点解释:
- 把f近似为多项式;
- 用f的高阶导数信息刻画 / 描述f.(微分学的任务)
下面讲讲应用.
/Claim/
∀x∈R,有
ex=n=0∑∞n!xn=n→∞lim(1+1!x1+⋯+n!xn+⋯)
/Claim/
∀x∈R,有
sinxcosx=x−3!x3+⋯=n=0∑∞(−1)n(2n+1)!x2n+1=1−2!x2+⋯=n=0∑∞(−1)n(2n)!x2n
实际上有的偏代数书籍中会用上面的式子定义三角函数,但是这样比较缺乏几何含义.
/Proof/
对称,只证明sinx的结论. 导函数是周期变化的,有
f(1)=cosx,f(2)=−sinx,f(3)=−cosx,f(4)=sinx
用 Taylor 公式,但是注意到展开至2n阶和2n+1阶,多项式部分并无区别,但是余项的阶数相差一阶,所以展开至2n+1阶比较划算.
由 Lagrange 余项得:
R2n+1(x)=(2n+1)!f(2n+1)(ξ)x2n+1=(2n+1)!±cosξx2n+1
注意ξ依赖于x和n. 只需证明上式在n→∞时趋于0,理由是
0≤∣R2n+1(x)∣≤(2n+1)!1∣x∣2n+1=bn
然后考虑到n→∞limbn+1/bn=0,所以n→∞limbn=0,根据夹逼定理证明余项趋于0. 得证.
证毕.
/Example/
求sinx的近似表达式,要求∀∣x∣<1,误差<10−3.
/Solution/
用f(x)=sinx的 Taylor 公式,得到
sinx=P(x)x−3!x3+⋯+(−1)n−1(2n−1)!x2n−1+(2n+1)!f(2n+1)(ξx)x2n+1
Feynman 在 Manhattan 计划期间,觉得非常无聊,于是挑战计划中的各种数学家,要他们出题测试自己,然后有一位数学家让他在 1 分钟内算出21/3,Feynman 发现这可以用(1.728+0.272)1/3来算,因为他熟知1728=123,并估算0.272/1.728∼0.1.
但是当有人问他 tan101000时,他无法回答了. 这位问倒 Feynman 的数学家曾经还用钢尺反复抽出打自己的手,并装作不痛,诱骗 Feynman 做同样的事(?)
为什么 Feynman 会熟知1728=123呢?因为 Ramanujan 的老师曾经生病,住在编号为1729的房间里,向他抱怨说这个房间号并没有什么特点. 但是 Ramanujan 在 1 分钟之内向他说明这个数字相当有特点,它是“最小的能被表达为两种不同立方和的数”(因为1729=103+93或者1729=123+13. 这个故事被很多人熟知.
/Theorem/
(1+x)α=n=0∑∞(αn)xn,x∈(−1,1)
/Proof/
我们会发现这个定理无法用 Lagrange 余项证明,因为 Lagrange 余项在这里写成:
Rn+1=(n+1)!α(α−1)⋯(α−n)(1+ξ)α−n−1xn+1
但是
(1+ξx)n+1
可能大于1,所以不能解决问题. 我们选择使用 Cauchy 余项,
Rn+1=n!f(n+1)(ξ)(x−ξ)n(x−0)
在叙述 Taylor 公式时,可以用另外一种方式来写出ξ,即定比分点的写法(这就是一种无量纲化!):
ξ−a=θ(b−a),0<θ<1
用θ改写ξ.
在这里就可以把 Cauchy 余项写成
Rn+1=n!α(α−1)⋯(α−n)(1+θx)α−n−1xn+1(1−θ)n
改写一下:
Rn+1=An!α(α−1)⋯(α−n)xn+1⋅B(1+θx1−θ)n⋅C(1+θx)α+1
分开考虑A、B、C三个部分,
- A部分,n→∞lim∣An+1∣/∣An∣=0,趋于零.
- B部分,因为0<θ(1+x),0<θ<1,可得0<1−θ<1+θx. 这导致B部分的底数小于1,所以B部分也趋于零.
- C部分,(1+θx)α−1有界,它介于(1−∣x∣)α−1和(1+∣x∣)α−1之间.
所以结合在一起得到Rn+1→0.
证毕.
从上面这个题能够看出,两种不同的余项虽然表达的是同样的含义,但是实际使用时并不相同,在不同情况下应该有选择地应用不同表达式.
下面来看一些更加低阶的例子.
/Example/ (一个流传广泛的例子)
设f在R上处处有二阶导数,且∣f(x)∣≤A,∣f′′(x)∣≤B,则∣f′(x)∣≤2AB.
/Proof/
用f和f(2)控制f(1),这显然不平凡. 但是注意到这三者均在 Taylor 公式中间出现,考虑二阶 Taylor 公式,得到
f(x+h)=f(x)+f′(x)h+2!f′′(ξ)h2
得到
∣f′(x)∣⟹∣f′(x)∣=∣hf(x+h)−f(x)−21f′′(ξ)h2∣≤h2A+2Bh∀h>0≤min(h2A+2Bh)
这里能用到均值不等式,上式在h=4A/B时取最小值,代回知∣f′(x)∣≤2AB.
但是这个结论不够强,为了强化结论,我们考虑下面的方法:
f(x+h)f(x−h)=f(x)+f′(x)h+2f′′(ξ)h2=f(x)−f′(x)h+2f′′(η)h2
作差得到
∣f′(x)∣=∣2hf(x+h)−f(x−h)−4f′′(ξ)−f′′(η)h2∣≤hA+2Bh
获得证明.
证毕.
下面一个例子来自于去年的期末考试.
上次出了这个题之后我意识到同学们的水平是很高的.
——艾神
/Example/
设f在R上有连续的二阶导,且f(0)=f(1)=0. 求证∀a∈(0,1)有∣f(a)∣≤81x∈[0,1]max∣f′′(x)∣.
/Proof/
证法 1:由∣f′′∣连续和最值定理,知道[0,1]max∣f′′(x)∣存在,记为M,目标是证明−M/8≤f(a)≤M/8,只看f(a)≤M/8.
考虑最值在何处取到,如果是在端点,那么肯定满足条件,所以只考虑不在端点取到最值的情况. 最值点是极值点,f′(x0)=0,Taylor 公式得到:
f(1)f(0)=f(x0)+f′(x0)(1−x0)+2f′′(ξ)(1−x0)2=f(x0)+f′(x0)(−x0)+2f′′(η)x02
做差得到
f(x0)=−2f′′(ξ)(1−x0)2=−2f′′(η)x02
于是
f(x0)≤min(21(1−x0)2M,21x02M)
但是(1−x0)和x0必有一项是∈[0,21],所以
f(x0)≤81M
证毕.
下节课讲其他证明方法.