高等微积分笔记 Lesson 21

本文最后更新于 2024年12月10日 晚上

高等微积分 Lesson 21

Taylor 公式的应用

ffII处处有nn阶导,则

f(b)=i=0n1f(i)(a)i!(ba)i+Rnf(b)=\sum_{i=0}^{n-1}\frac{f^{(i)}(a)}{i!}(b-a)^i+R_n

其中RnR_n可以表示为

  • Lagrange 余项,ξ\exist\xi介于a,ba,b之间使得

    Rn=f(n)(ξ)n!(ba)nR_n=\frac{f^{(n)}(\xi)}{n!}(b-a)^n

  • Cauchy 余项,ξ\exist\xi介于a,ba,b之间使得

    Rn=f(n)(ξ)(n1)!(bξ)n1(ba)R_n=\frac{f^{(n)}(\xi)}{(n-1)!}(b-\xi)^{n-1}(b-a)

记忆方法:对于多项式PPdegPn1\deg P\leq n-1,可以推出

P(b)=i=0n1P(i)(a)i!(ba)iP(b)=\sum_{i=0}^{n-1}\frac{P^{(i)}(a)}{i!}(b-a)^i

此亦为 Taylor 公式,对一般的f(x)f(x),可以定义

F(x)=i=0n1f(i)(x)i!(bx)iF(x)=\sum_{i=0}^{n-1}\frac{f^{(i)}(x)}{i!}(b-x)^i

FF和任意的一个g(x)g(x),可以用 Cauchy 中值定理,产生

F(b)F(a)g(b)g(a)=F(ξ)g(ξ)f(b)i=0n1f(i)(a)i!(ba)i=g(b)g(a)g(ξ)f(n)(ξ)(n1)!(bξ)n1\begin{aligned} \frac{F(b)-F(a)}{g(b)-g(a)}&=\frac{F'(\xi)}{g'(\xi)}\\\\ \Longrightarrow f(b)-\sum_{i=0}^{n-1}\frac{f^{(i)}(a)}{i!}(b-a)^i&=\frac{g(b)-g(a)}{g'(\xi)}\frac{f^{(n)}(\xi)}{(n-1)!}(b-\xi)^{n-1} \end{aligned}

得到任何一种余项. 对于 Cauchy 余项,就是取g(x)=xg(x)=x.

为了更好的应用 Taylor 公式,我们从不同角度解读公式的含义,有以下几点解释:

  1. ff近似为多项式;
  2. ff的高阶导数信息刻画 / 描述ff.(微分学的任务)

下面讲讲应用.

/Claim/

xR\forall x\in\R,有

ex=n=0xnn!=limn(1+x11!++xnn!+)e^x=\sum_{n=0}^\infty\frac{x^n}{n!}=\lim_{n\to\infty}(1+\frac{x^1}{1!}+\cdots+\frac{x^n}{n!}+\cdots)

/Claim/

xR\forall x\in\R,有

sinx=xx33!+=n=0(1)nx2n+1(2n+1)!cosx=1x22!+=n=0(1)nx2n(2n)!\begin{aligned} \sin x&=x-\frac{x^3}{3!}+\cdots=\sum_{n=0}^\infty(-1)^n\frac{x^{2n+1}}{(2n+1)!}\\ \cos x&=1-\frac{x^2}{2!}+\cdots=\sum_{n=0}^\infty(-1)^n\frac{x^{2n}}{(2n)!} \end{aligned}

实际上有的偏代数书籍中会用上面的式子定义三角函数,但是这样比较缺乏几何含义.

/Proof/

对称,只证明sinx\sin x的结论. 导函数是周期变化的,有

f(1)=cosx,f(2)=sinx,f(3)=cosx,f(4)=sinxf^{(1)}=\cos x\,,\quad f^{(2)}=-\sin x\,,\quad f^{(3)}=-\cos x\,,\quad f^{(4)}=\sin x

用 Taylor 公式,但是注意到展开至2n2n阶和2n+12n+1阶,多项式部分并无区别,但是余项的阶数相差一阶,所以展开至2n+12n+1阶比较划算.

由 Lagrange 余项得:

R2n+1(x)=f(2n+1)(ξ)(2n+1)!x2n+1=±cosξ(2n+1)!x2n+1R_{2n+1}(x)=\frac{f^{(2n+1)}(\xi)}{(2n+1)!}x^{2n+1}=\frac{\pm\cos\xi}{(2n+1)!}x^{2n+1}

注意ξ\xi依赖于xxnn. 只需证明上式在nn\to\infty时趋于00,理由是

0R2n+1(x)1(2n+1)!x2n+1=bn0\leq|R_{2n+1}(x)|\leq\frac{1}{(2n+1)!}|x|^{2n+1}=b_n

然后考虑到limnbn+1/bn=0\underset{n\to\infty}{\lim}b_{n+1}/b_n=0,所以limnbn=0\underset{n\to\infty}{\lim}b_n=0,根据夹逼定理证明余项趋于00. 得证.

证毕.

/Example/

sinx\sin x的近似表达式,要求x<1\forall|x|<1,误差<103<10^{-3}.

/Solution/

f(x)=sinxf(x)=\sin x的 Taylor 公式,得到

sinx=xx33!++(1)n1x2n1(2n1)!P(x)+f(2n+1)(ξx)(2n+1)!x2n+1\sin x=\underset{P(x)}{\underline{x-\frac{x^3}{3!}+\cdots+(-1)^{n-1}\frac{x^{2n-1}}{(2n-1)!}}}+\frac{f^{(2n+1)}(\xi_x)}{(2n+1)!}x^{2n+1}

Feynman 在 Manhattan 计划期间,觉得非常无聊,于是挑战计划中的各种数学家,要他们出题测试自己,然后有一位数学家让他在 1 分钟内算出21/32^{1/3},Feynman 发现这可以用(1.728+0.272)1/3(1.728+0.272)^{1/3}来算,因为他熟知1728=1231728=12^3,并估算0.272/1.7280.10.272/1.728\sim0.1.

但是当有人问他 tan101000\tan10^{1000}时,他无法回答了. 这位问倒 Feynman 的数学家曾经还用钢尺反复抽出打自己的手,并装作不痛,诱骗 Feynman 做同样的事(?)

为什么 Feynman 会熟知1728=1231728=12^3呢?因为 Ramanujan 的老师曾经生病,住在编号为17291729的房间里,向他抱怨说这个房间号并没有什么特点. 但是 Ramanujan 在 1 分钟之内向他说明这个数字相当有特点,它是“最小的能被表达为两种不同立方和的数”(因为1729=103+931729=10^3+9^3或者1729=123+131729=12^3+1^3. 这个故事被很多人熟知.

我们刚刚在讲什么来着?
——艾神

/Theorem/

(1+x)α=n=0(αn)xn,x(1,1)(1+x)^\alpha=\sum_{n=0}^\infty\begin{pmatrix} \alpha\\n \end{pmatrix}x^n\,,\quad x\in(-1,1)

/Proof/

我们会发现这个定理无法用 Lagrange 余项证明,因为 Lagrange 余项在这里写成:

Rn+1=α(α1)(αn)(n+1)!(1+ξ)αn1xn+1R_{n+1}=\frac{\alpha(\alpha-1)\cdots(\alpha-n)}{(n+1)!}(1+\xi)^{\alpha-n-1}x^{n+1}

但是

(x1+ξ)n+1(\frac{x}{1+\xi})^{n+1}

可能大于11,所以不能解决问题. 我们选择使用 Cauchy 余项,

Rn+1=f(n+1)(ξ)n!(xξ)n(x0)R_{n+1}=\frac{f^{(n+1)}(\xi)}{n!}(x-\xi)^n(x-0)

在叙述 Taylor 公式时,可以用另外一种方式来写出ξ\xi,即定比分点的写法(这就是一种无量纲化!):

ξa=θ(ba),0<θ<1\xi-a=\theta(b-a)\,,0<\theta<1

θ\theta改写ξ\xi.

在这里就可以把 Cauchy 余项写成

Rn+1=α(α1)(αn)n!(1+θx)αn1xn+1(1θ)nR_{n+1}=\frac{\alpha(\alpha-1)\cdots(\alpha-n)}{n!}(1+\theta x)^{\alpha-n-1}x^{n+1}(1-\theta)^n

改写一下:

Rn+1=α(α1)(αn)n!xn+1A(1θ1+θx)nB(1+θx)α+1CR_{n+1}=\underset{A}{\underline{\frac{\alpha(\alpha-1)\cdots(\alpha-n)}{n!}x^{n+1}}}\cdot\underset{B}{\underline{(\frac{1-\theta}{1+\theta x})^n}}\cdot\underset{C}{\underline{(1+\theta x)^{\alpha+1}}}

分开考虑AABBCC三个部分,

  1. AA部分,limnAn+1/An=0\underset{n\to\infty}{\lim}|A_{n+1}|/|A_n|=0,趋于零.
  2. BB部分,因为0<θ(1+x)0<\theta(1+x)0<θ<10<\theta<1,可得0<1θ<1+θx0<1-\theta<1+\theta x. 这导致BB部分的底数小于11,所以BB部分也趋于零.
  3. CC部分,(1+θx)α1(1+\theta x)^{\alpha-1}有界,它介于(1x)α1(1-|x|)^{\alpha-1}(1+x)α1(1+|x|)^{\alpha-1}之间.

所以结合在一起得到Rn+10R_{n+1}\to0.

证毕.

从上面这个题能够看出,两种不同的余项虽然表达的是同样的含义,但是实际使用时并不相同,在不同情况下应该有选择地应用不同表达式.

下面来看一些更加低阶的例子.

/Example/ (一个流传广泛的例子)

ffR\R上处处有二阶导数,且f(x)A|f(x)|\leq Af(x)B|f''(x)|\leq B,则f(x)2AB|f'(x)|\leq\sqrt{2AB}.

/Proof/

fff(2)f^{(2)}控制f(1)f^{(1)},这显然不平凡. 但是注意到这三者均在 Taylor 公式中间出现,考虑二阶 Taylor 公式,得到

f(x+h)=f(x)+f(x)h+f(ξ)2!h2f(x+h)=f(x)+f'(x)h+\frac{f''(\xi)}{2!}h^2

得到

f(x)=f(x+h)f(x)12f(ξ)h2h2Ah+Bh2h>0f(x)min(2Ah+Bh2)\begin{aligned} |f'(x)|&=|\frac{f(x+h)-f(x)-\frac{1}{2}f''(\xi)h^2}{h}|\\\\ &\leq\frac{2A}{h}+\frac{Bh}{2}\quad\forall h>0\\\\ \Longrightarrow|f'(x)|&\leq\min(\frac{2A}{h}+\frac{Bh}{2}) \end{aligned}

这里能用到均值不等式,上式在h=4A/Bh=\sqrt{4A/B}时取最小值,代回知f(x)2AB|f'(x)|\leq2\sqrt{AB}.

但是这个结论不够强,为了强化结论,我们考虑下面的方法:

f(x+h)=f(x)+f(x)h+f(ξ)2h2f(xh)=f(x)f(x)h+f(η)2h2\begin{aligned} f(x+h)&=f(x)+f'(x)h+\frac{f''(\xi)}{2}h^2\\ f(x-h)&=f(x)-f'(x)h+\frac{f''(\eta)}{2}h^2 \end{aligned}

作差得到

f(x)=f(x+h)f(xh)2hf(ξ)f(η)4h2Ah+Bh2|f'(x)|=|\frac{f(x+h)-f(x-h)}{2h}-\frac{f''(\xi)-f''(\eta)}{4}h^2|\leq\frac{A}{h}+\frac{Bh}{2}

获得证明.

证毕.

下面一个例子来自于去年的期末考试.

上次出了这个题之后我意识到同学们的水平是很高的.
——艾神

/Example/

ffR\R上有连续的二阶导,且f(0)=f(1)=0f(0)=f(1)=0. 求证a(0,1)\forall a\in(0,1)f(a)18maxx[0,1]f(x)|f(a)|\leq\frac{1}{8}\underset{x\in[0,1]}{\max}|f''(x)|.

/Proof/

证法 1:由f|f''|连续和最值定理,知道max[0,1]f(x)\underset{[0,1]}{\max}|f''(x)|存在,记为MM,目标是证明M/8f(a)M/8-M/8\leq f(a)\leq M/8,只看f(a)M/8f(a)\leq M/8.

考虑最值在何处取到,如果是在端点,那么肯定满足条件,所以只考虑不在端点取到最值的情况. 最值点是极值点,f(x0)=0f'(x_0)=0,Taylor 公式得到:

f(1)=f(x0)+f(x0)(1x0)+f(ξ)2(1x0)2f(0)=f(x0)+f(x0)(x0)+f(η)2x02\begin{aligned} f(1)&=f(x_0)+f'(x_0)(1-x_0)+\frac{f''(\xi)}{2}(1-x_0)^2\\ f(0)&=f(x_0)+f'(x_0)(-x_0)+\frac{f''(\eta)}{2}x_0^2 \end{aligned}

做差得到

f(x0)=f(ξ)2(1x0)2=f(η)2x02f(x_0)=-\frac{f''(\xi)}{2}(1-x_0)^2=-\frac{f''(\eta)}{2}x_0^2

于是

f(x0)min(12(1x0)2M,12x02M)f(x_0)\leq\min(\frac{1}{2}(1-x_0)^2M,\frac{1}{2}x_0^2M)

但是(1x0)(1-x_0)x0x_0必有一项是[0,12]\in[0,\frac{1}{2}],所以

f(x0)18Mf(x_0)\leq\frac{1}{8}M

证毕.


下节课讲其他证明方法.


高等微积分笔记 Lesson 21
https://physnya.top/2024/12/04/integral21/
作者
菲兹克斯喵
发布于
2024年12月4日
更新于
2024年12月10日
许可协议