泰勒展开及其应用

您所在的位置:网站首页 常用泰勒展开式公式的应用 泰勒展开及其应用

泰勒展开及其应用

2024-07-09 03:43| 来源: 网络整理| 查看: 265

泰勒展开[1]

在实际应用中对于具有复杂形式的函数我们常常希望用较为简单的函数形式表示他,而多项式就是这种简单的形式。比如对于指数函数、三角函数,我们可以使用多项式来逼近。

为了逼近(或者说是仿造)目标函数曲线f(x),首先选择一个切入点(x0,f(f0)),然后让此处的增减性相同,即一阶导数相同。再使其凹凸性相同(二阶导数相同)。然后让更高级的特性相同。。。 因此,整体思路就是让仿造的函数g(x)与f(x)的初始点和从一阶到高阶的导数都相等。然而,在实际操作中并不能无限阶求导,只需要按需(或经验)选择阶数。 假设我们只需要算到n阶,多项式函数为: \(g(x)=a_0+a_1(x-x_0)+a_2(x-x_0)^2+……+a_n(x-x_0)^n\) 由初始点相等: \(g(x_0)=f(x_0)\),得\(a_0=f(x_0)\) 由n阶导数相等: \(g^n(x_0)=a_n n!=f^n(x_0)\),得\(a_n={f^n(x_0)\over n!}\) 因此求得多项式:

\[g(x)\approx f(x_0)+\frac{f^1(x_0)}{1!}(x-x_0)+\frac{f^2(x_0)}{2!}(x-x_0)^2+……+\frac{f^n(x_0)}{n!}(x-x_0)^n \]

求泰勒展开式需要求g(x)的各项系数,计算比较复杂,但是一旦求出展开式之后就可以很方便地使用多项式计算了, 相当于一劳永逸.

误差项,即余项: R(x)=f(x)-g(x)=\(\frac{f^{(n+1)}(x_0)}{(n+1)!}(x-x_0)^{n+1}+\frac{f^{(n+2)}(x_0)}{(n+2)!}(x-x_0)^{n+2}+\cdots +\frac{f^{(\infty)}(x_0)}{\infty!}(x-x_0)^{\infty}\) 利用柯西中值定理,可以求出\(R(x)=\frac{f^{(n+1)}(\xi)}{(n+1)!}(x-x_0)^{n+1}\)

应用举例 优化算法之牛顿法

牛顿法原理:

将 \(f(x)\) 在 \(x_n\) 处二阶 taylor 展开, 有:

\[f(x)\approx f(x_n)+f'(x_n)(x-x_n)+{1\over 2}f''(x_n)(x-x_n)^2 \\ f(x_n+\Delta x)\approx f(x_n)+f'(x_n)\Delta x+{1\over 2}f''(x_n)\Delta x^2 \]

我们的目的是选择合适的 \(\Delta x\) 最小化目标函数 \(f(x_n+\Delta x)\) . 求导令导数为0得到:

\[{d f(x_n+\Delta x) \over d\Delta x}=f'(x_n)+f''(x_n)\Delta x=0 \\ \Longrightarrow \Delta x=-{f'(x_n)\over f''(x_n)} \\ x_{n+1}=x_n+\Delta x=x_n-[f''(x_n)]^{-1}f'(x_n) \]

其中的二阶项 \(f''(x_n)\) 就是 Hessian 矩阵.

基尼指数

用基尼指数代替概率分布的信息熵

概率分布的信息熵等于 \(H(p)=-\sum_{k=1}^K p_k\log p_k\). 将 log(p) 在 p=0.5 处泰勒展开, 得到\(\log(p)\approx \log(0.5)+\log'(0.5)\cdot (p-0.5)+O(p^2)\approx 2(p-1)\) 将 p=0.5 处泰勒展开式代入 H(p) 可得\(H(p)\approx 2\cdot \rm{Gini}(p)\).

平方根 [2]

求正数的平方根, \(f(x)=\sqrt x\), 让f(x)在x=1处泰勒展开得到: \(f(x)=\sqrt x \approx 1+{1\over 2}(x-1)-{1\over 4}{(x-1)^2\over 2!}+{3\over 8}{(x-1)^3\over 3!}-\cdots\) 这样我们可以进行在一定精度内的逼近, 但是这里存在收敛区间的问题。只有 0 < x



【本文地址】


今日新闻


推荐新闻


CopyRight 2018-2019 办公设备维修网 版权所有 豫ICP备15022753号-3