泰勒公式







指数函数y=ex{displaystyle y=e^{x}}y=e^{{x}}(红色实线)与在原点展开的泰勒多项式前四项(绿色虚线)。在这个函数中,泰勒多项式展开的项数越多,曲线拟合的越好。



在数学中,泰勒公式英语:Taylor's Formula)是一个用函数在某点的信息描述其附近取值的公式。這個公式來自於微積分的泰勒定理Taylor's theorem),泰勒定理描述了一個可微函數,如果函数足够光滑的话,在已知函数在某一点的各阶导数值的情况之下,泰勒公式可以用这些导数值做系数构建一个多项式来近似函数在这一点的邻域中的值,這個多項式稱為泰勒多項式Taylor polynomial)。泰勒公式还给出了餘項即这个多项式和实际的函数值之间的偏差。泰勒公式得名于英国数学家布鲁克·泰勒。他在1712年的一封信里首次叙述了这个公式,尽管1671年詹姆斯·格雷高里已经发现了它的特例[1]。拉格朗日在1797年之前,最先提出了帶有餘項的現在形式的泰勒定理。





目录





  • 1 泰勒公式


  • 2 泰勒定理


  • 3 餘項估計


  • 4 多元泰勒公式


  • 5 參閱


  • 6 参考来源





泰勒公式


泰勒公式的初衷是用多项式来近似表示函数在某点周围的情况。比如说,指数函数exx = 0 的附近可以用以下多项式来近似地表示:



ex≈1+x+x22!+x33!+⋯+xnn!.{displaystyle {textrm {e}}^{x}approx 1+x+{frac {x^{2}}{2!}}+{frac {x^{3}}{3!}}+cdots +{frac {x^{n}}{n!}}.}{textrm  {e}}^{x}approx 1+x+{frac  {x^{2}}{2!}}+{frac  {x^{3}}{3!}}+cdots +{frac  {x^{n}}{n!}}.


称为指数函数在0处的n 泰勒展开公式。这个公式只对0附近的x 有用,x0 越远,这个公式就越不准确。实际函数值和多项式的偏差称为泰勒公式的余项



Rn(x)=ex−(1+x+x22!+x33!+⋯+xnn!).{displaystyle R_{n}(x)={textrm {e}}^{x}-left(1+x+{frac {x^{2}}{2!}}+{frac {x^{3}}{3!}}+cdots +{frac {x^{n}}{n!}}right).}R_{n}(x)={textrm  {e}}^{x}-left(1+x+{frac  {x^{2}}{2!}}+{frac  {x^{3}}{3!}}+cdots +{frac  {x^{n}}{n!}}right).



泰勒定理


对于一般的函数,泰勒公式的系数的选择依赖于函数在一点的各阶导数值。这个想法的原由可以由微分的定义开始。微分是函数在一点附近的最佳线性近似:




f(a+h)=f(a)+f′(a)h+o(h){displaystyle f(a+h)=f(a)+f^{prime }(a)h+o(h)}f(a+h)=f(a)+f^{{prime }}(a)h+o(h),其中o(h){displaystyle o(h)}o(h) 是比h 高阶的无穷小。


也就是说f(a+h)≈f(a)+f′(a)h{displaystyle f(a+h)approx f(a)+f^{prime }(a)h}f(a+h)approx f(a)+f^{{prime }}(a)h,或f(x)≈f(a)+f′(a)(x−a){displaystyle f(x)approx f(a)+f^{prime }(a)(x-a)}f(x)approx f(a)+f^{{prime }}(a)(x-a)


注意到f(x){displaystyle f(x)}f(x)f(a)+f′(a)(x−a){displaystyle f(a)+f^{prime }(a)(x-a)}f(a)+f^{{prime }}(a)(x-a)a 处的零阶导数和一阶导数都相同。对足够光滑的函数,如果一个多项式在a 处的前n 次导数值都与函数在a 处的前n 次导数值重合,那么这个多项式应该能很好地近似描述函数在a 附近的情况。以下定理说明这是正确的:




定理


n 是一个正整数。如果定义在一个包含 a 的区间上的函数 fa 点处 n+1 次可导,那么对于这个区间上的任意 x,都有:



f(x)=f(a)+f′(a)1!(x−a)+f(2)(a)2!(x−a)2+⋯+f(n)(a)n!(x−a)n+Rn(x).{displaystyle f(x)=f(a)+{frac {f'(a)}{1!}}(x-a)+{frac {f^{(2)}(a)}{2!}}(x-a)^{2}+cdots +{frac {f^{(n)}(a)}{n!}}(x-a)^{n}+R_{n}(x).}f(x)=f(a)+{frac  {f'(a)}{1!}}(x-a)+{frac  {f^{{(2)}}(a)}{2!}}(x-a)^{2}+cdots +{frac  {f^{{(n)}}(a)}{n!}}(x-a)^{n}+R_{n}(x).[2]


其中的多项式称为函数在a 处的泰勒展开式,剩余的Rn(x){displaystyle R_{n}(x)}R_{n}(x) 是泰勒公式的余项,是(x−a)n{displaystyle (x-a)^{n}}(x-a)^{n} 的高阶无穷小。





Rn(x){displaystyle R_{n}(x)}R_{n}(x) 的表达形式有若干种,分别以不同的数学家命名。


带有皮亚诺型余项的泰勒公式说明了多项式和函数的接近程度:


f(x)=f(a)+f′(a)1!(x−a)+f(2)(a)2!(x−a)2+⋯+f(n)(a)n!(x−a)n+o[(x−a)n]{displaystyle f(x)=f(a)+{frac {f'(a)}{1!}}(x-a)+{frac {f^{(2)}(a)}{2!}}(x-a)^{2}+cdots +{frac {f^{(n)}(a)}{n!}}(x-a)^{n}+o[(x-a)^{n}]}f(x)=f(a)+{frac  {f'(a)}{1!}}(x-a)+{frac  {f^{{(2)}}(a)}{2!}}(x-a)^{2}+cdots +{frac  {f^{{(n)}}(a)}{n!}}(x-a)^{n}+o[(x-a)^{{n}}]

也就是说,当x 无限趋近a 时,余项Rn(x){displaystyle R_{n}(x)}R_{n}(x) 将会是(x−a)n{displaystyle (x-a)^{n}}(x-a)^{{n}} 的高阶无穷小,或者说多项式和函数的误差将远小于(x−a)n{displaystyle (x-a)^{n}}(x-a)^{{n}}[3]。这个结论可以由下面更强的结论推出。


带有拉格朗日型余项的泰勒公式可以视为拉格朗日微分中值定理的推广:


f(x)=f(a)+f′(a)1!(x−a)+f(2)(a)2!(x−a)2+⋯+f(n)(a)n!(x−a)n+f(n+1)(θ)(n+1)!(x−a)(n+1){displaystyle f(x)=f(a)+{frac {f'(a)}{1!}}(x-a)+{frac {f^{(2)}(a)}{2!}}(x-a)^{2}+cdots +{frac {f^{(n)}(a)}{n!}}(x-a)^{n}+{frac {f^{(n+1)}(theta )}{(n+1)!}}(x-a)^{(n+1)}}f(x)=f(a)+{frac  {f'(a)}{1!}}(x-a)+{frac  {f^{{(2)}}(a)}{2!}}(x-a)^{2}+cdots +{frac  {f^{{(n)}}(a)}{n!}}(x-a)^{n}+{frac  {f^{{(n+1)}}(theta )}{(n+1)!}}(x-a)^{{(n+1)}}

Rn(x)=f(n+1)(θ)(n+1)!(x−a)(n+1){displaystyle R_{n}(x)={frac {f^{(n+1)}(theta )}{(n+1)!}}(x-a)^{(n+1)}}R_{n}(x)={frac  {f^{{(n+1)}}(theta )}{(n+1)!}}(x-a)^{{(n+1)}},其中θ(a,x){displaystyle theta in (a,x)}theta in (a,x)[4]


带有积分型余项的泰勒公式可以看做微积分基本定理的推广[5]



Rn(x)=∫axf(n+1)(t)n!(x−t)ndt,{displaystyle R_{n}(x)=int _{a}^{x}{frac {f^{(n+1)}(t)}{n!}}(x-t)^{n},dt,}R_{n}(x)=int _{a}^{x}{frac  {f^{{(n+1)}}(t)}{n!}}(x-t)^{n},dt,



餘項估計


拉格朗日型余项或积分型余项可以帮助估计泰勒展开式和函数在一定区间之内的误差。设函数在区间[ar, a + r]n 次连续可微并且在区间(ar, a + r)n + 1 次可导。如果存在正实数Mn 使得区间(ar, a + r) 里的任意x 都有 |f(n+1)(x)|≤Mn{displaystyle |f^{(n+1)}(x)|leq M_{n}}|f^{{(n+1)}}(x)|leq M_{n},那么:


f(x)=f(a)+f′(a)1!(x−a)+f(2)(a)2!(x−a)2+⋯+f(n)(a)n!(x−a)n+Rn(x),{displaystyle f(x)=f(a)+{frac {f'(a)}{1!}}(x-a)+{frac {f^{(2)}(a)}{2!}}(x-a)^{2}+cdots +{frac {f^{(n)}(a)}{n!}}(x-a)^{n}+R_{n}(x),}f(x)=f(a)+{frac  {f'(a)}{1!}}(x-a)+{frac  {f^{{(2)}}(a)}{2!}}(x-a)^{2}+cdots +{frac  {f^{{(n)}}(a)}{n!}}(x-a)^{n}+R_{n}(x),

其中|Rn(x)|≤Mnrn+1(n+1)!{displaystyle |R_{n}(x)|leq M_{n}{frac {r^{n+1}}{(n+1)!}}}|R_{n}(x)|leq M_{n}{frac  {r^{{n+1}}}{(n+1)!}}。这个上界估计对区间(ar, a + r) 里的任意x 都成立,是一个一致估计。


如果当n 趋向于无穷大时,还有Mnrn+1(n+1)!→0{displaystyle M_{n}{frac {r^{n+1}}{(n+1)!}}rightarrow 0}M_{n}{frac  {r^{{n+1}}}{(n+1)!}}rightarrow 0,那么可以推出 Rn(x)→0{displaystyle R_{n}(x)rightarrow 0}R_{n}(x)rightarrow 0f 是区间(ar, a + r) 上解析函数。f 在区间(ar, a + r) 上任一点的值都等于在这一点的泰勒展开式的极限。



多元泰勒公式


对于多元函数,也有类似的泰勒公式。设B(a, r ) 是欧几里得空间RN 中的开球,ƒ 是定义在B(a, r ) 的闭包上的实值函数,并在每一点都存在所有的n+1 次偏导数。这时的泰勒公式为:



对所有x∈B(a,r){displaystyle xin mathbf {B} (a,r)}xin {mathbf  {B}}(a,r)

f(x)=∑|=0n1α!∂αf(a)∂(x−a)α+∑|=n+1Rα(x)(x−a)α{displaystyle f(x)=sum _{|alpha |=0}^{n}{frac {1}{alpha !}}{frac {partial ^{alpha }f(a)}{partial x^{alpha }}}(x-a)^{alpha }+sum _{|alpha |=n+1}R_{alpha }(x)(x-a)^{alpha }}f(x)=sum _{{|alpha |=0}}^{n}{frac  {1}{alpha !}}{frac  {partial ^{alpha }f(a)}{partial x^{alpha }}}(x-a)^{alpha }+sum _{{|alpha |=n+1}}R_{{alpha }}(x)(x-a)^{alpha }



其中的 α{displaystyle alpha } alpha 是多重指标,即|=α1+α2+...+αn{displaystyle |alpha |=alpha _{1}+alpha _{2}+...+alpha _{n}}|alpha |=alpha _{1}+alpha _{2}+...+alpha _{n} , α!=α1!⋅α2!⋅...⋅αn!{displaystyle alpha !=alpha _{1}!cdot alpha _{2}!cdot ...cdot alpha _{n}!}{displaystyle alpha !=alpha _{1}!cdot alpha _{2}!cdot ...cdot alpha _{n}!}

x=(x1,x2,...,xn){displaystyle x=(x_{1},x_{2},...,x_{n})}x=(x_{1},x_{2},...,x_{n}) ,则记


=x1α1x2α2...xnαn{displaystyle x^{alpha }=x_{1}^{alpha _{1}}x_{2}^{alpha _{2}}...x_{n}^{alpha _{n}}}x^{{alpha }}=x_{1}^{{alpha _{1}}}x_{2}^{{alpha _{2}}}...x_{n}^{{alpha _{n}}}.


αf(a)∂=∂α1+α2+...+αnf(a)∂x1α1x2α2...xnαn{displaystyle {frac {partial ^{alpha }f(a)}{partial x^{alpha }}}={frac {partial ^{alpha _{1}+alpha _{2}+...+alpha _{n}}f(a)}{partial x_{1}^{alpha _{1}}x_{2}^{alpha _{2}}...x_{n}^{alpha _{n}}}}}{displaystyle {frac {partial ^{alpha }f(a)}{partial x^{alpha }}}={frac {partial ^{alpha _{1}+alpha _{2}+...+alpha _{n}}f(a)}{partial x_{1}^{alpha _{1}}x_{2}^{alpha _{2}}...x_{n}^{alpha _{n}}}}}.


其中的余项也满足不等式:



对所有满足 |α| = n + 1 的 α,|Rα(x)|≤supy∈|1α!∂αf(y)∂|{displaystyle |R_{alpha }(x)|leq sup _{yin {bar {B}}}left|{frac {1}{alpha !}}{frac {partial ^{alpha }f(y)}{partial x^{alpha }}}right|}|R_{{alpha }}(x)|leq sup _{{yin {bar  {B}}}}left|{frac  {1}{alpha !}}{frac  {partial ^{alpha }f(y)}{partial x^{alpha }}}right|


特别地,多元形式的泰勒公式可表示为:



f(a+h)=∑k=0m∑|=k∂αf(a)∂+Rm{displaystyle f(a+h)=sum _{k=0}^{m}sum _{|alpha |=k}{frac {partial ^{alpha }f(a)}{partial x^{alpha }}}h^{alpha }+R_{m}}{displaystyle f(a+h)=sum _{k=0}^{m}sum _{|alpha |=k}{frac {partial ^{alpha }f(a)}{partial x^{alpha }}}h^{alpha }+R_{m}}

其中Rm=∑|=m+1∂αf(a+θh)α!hα{displaystyle R_{m}=sum _{|alpha |=m+1}{frac {partial ^{alpha }f(a+theta h)}{alpha !}}h^{alpha }}{displaystyle R_{m}=sum _{|alpha |=m+1}{frac {partial ^{alpha }f(a+theta h)}{alpha !}}h^{alpha }}.


在应用上述公式时,特别重要的是展开式的前三项,即:




f(a+h)=f(a)+∂f∂x1(a)h1+...+∂f∂xn(a)hn+12∑i,j=1n∂2f∂xi∂xj(a)+...{displaystyle f(a+h)=f(a)+{frac {partial f}{partial x_{1}}}(a)h_{1}+...+{frac {partial f}{partial x_{n}}}(a)h_{n}+{frac {1}{2}}sum _{i,j=1}^{n}{frac {partial ^{2}f}{partial x_{i}partial x_{j}}}(a)+...}{displaystyle f(a+h)=f(a)+{frac {partial f}{partial x_{1}}}(a)h_{1}+...+{frac {partial f}{partial x_{n}}}(a)h_{n}+{frac {1}{2}}sum _{i,j=1}^{n}{frac {partial ^{2}f}{partial x_{i}partial x_{j}}}(a)+...}.

运用雅可比矩阵与海森矩阵,则上式可表示为:

f(a+h)=f(a)+Jf(a)h+12(h1,...hn)Hf(a)(h1...hn)+...{displaystyle f(a+h)=f(a)+Jf(a)h+{frac {1}{2}}(h_{1},...h_{n})Hf(a){begin{pmatrix}h_{1}\...\h_{n}\end{pmatrix}}+...}{displaystyle f(a+h)=f(a)+Jf(a)h+{frac {1}{2}}(h_{1},...h_{n})Hf(a){begin{pmatrix}h_{1}\...\h_{n}\end{pmatrix}}+...}

其中Jf(a){displaystyle Jf(a)}{displaystyle Jf(a)} 为雅可比矩阵,Hf(a){displaystyle Hf(a)}{displaystyle Hf(a)} 为海森矩阵.



參閱



  • 泰勒級數

  • 拉格朗日型餘項

  • 佩亞諾型餘項

  • 麥克勞林公式



参考来源




  1. ^ J J O'Connor and E F Robertson. Brook Taylor's Biography. 


  2. ^ Rudin, 第123至124页.


  3. ^ 《微积分(Ⅱ)》第88-90页.


  4. ^ Klein (1998) 20.3; Apostol (1967) 7.7


  5. ^ Protter, Morrey, 第135-136页




  • Apostol, Tom. 《微积分学》(Calculus). Jon Wiley & Sons, Inc. 1967. ISBN 0-471-00005-1. 


  • Klein, Morris. 《微积分学:直观物理方法》(Calculus: An Intuitive and Physical Approach). Dover. 1998. ISBN 0-486-40453-6. 


  • Walter Rudin. 《数学分析原理》(Principles of Mathematical Analysis). Mcgraw-hill Book Company. 1976. ISBN 978-0-070-54235-8. 


  • 清华大学数学科学系微积分编写组. 《微积分(Ⅱ)》. 清华大学出版社. 2000. ISBN 7-302-06917-4. 


  • Murray H. Protter,Charles Bradfield Morrey. 《中等微积分》(Intermediate calculus). Springer. 1986. ISBN 978-0-387-96058-6. 



這個網誌中的熱門文章

請回答1994

FM

龍潭區 (台灣)