泰勒展开

1.泰勒定理

若函数f(x)f(x)f(x)在x0x_0x0​处n阶可微,则
f(x)=∑k=0nf(k)(x)k!(x−x0)k+Rn(x)f(x) = \sum_{k=0}^{n}\frac{f^{(k)}(x)}{k!}(x-x_0)^k+R_n(x) f(x)=k=0∑n​k!f(k)(x)​(x−x0​)k+Rn​(x)
其中,Rn(x)R_n(x)Rn​(x)称为f(x)f(x)f(x)的余项,常用的余项公式如下所示:

佩亚诺型余项:Rn(x)=o((x−x0)n)R_n(x) = o((x-x_0)^n)Rn​(x)=o((x−x0​)n)

拉格朗日型余项:
Rn(x)=f(n+1)(ξ)(n+1)!(x−x0)(n+1)R_n(x)=\frac{f^{(n+1)}(\xi)}{(n+1)!}(x-x_0)^{(n+1)} Rn​(x)=(n+1)!f(n+1)(ξ)​(x−x0​)(n+1)
其中,ξ\xiξ是介于x与x0x_0x0​之间的一个数,特别的,当x0=0x_0 = 0x0​=0时的带拉格朗日余项的泰勒公式如下:
f(x)=f(0)+f′(0)x+f′′(0)2!x2+⋯+f(n)(0)n!xn+f(n+1)(ξ)(n+1)!xn+1,(0<ξ<x)f(x)=f(0)+f^{\prime}(0) x+\frac{f^{\prime \prime}(0)}{2 !} x^{2}+\cdots+\frac{f^{(n)}(0)}{n !} x^{n}+\frac{f^{(n+1)}(\xi)}{(n+1) !} x^{n+1},(0<\xi<x) f(x)=f(0)+f′(0)x+2!f′′(0)​x2+⋯+n!f(n)(0)​xn+(n+1)!f(n+1)(ξ)​xn+1,(0<ξ<x)
该方程称为麦克劳林公式。

下面给出几种常用的带拉格朗日余项的泰勒公式展开:

1)
ex=1+x+x22!+...+xnn!+e(θx)(n+1)!x(n+1)e^x=1+x+\frac{x^2}{2!}+...+\frac{x^n}{n!}+\frac{e^{(\theta x)}}{(n+1)!}x^{(n+1)} ex=1+x+2!x2​+...+n!xn​+(n+1)!e(θx)​x(n+1)
2)
sin(x)=x−x33!+...+(−1)(n−1)x(2n−1)(2n−1)!+(−1)nx(2n+1)(2n+1)!sin(x)=x-\frac{x^3}{3!}+...+(-1)^{(n-1)}\frac{x^{(2n-1)}}{(2n-1)!}+(-1)^n\frac{x^{(2n+1)}}{(2n+1)!} sin(x)=x−3!x3​+...+(−1)(n−1)(2n−1)!x(2n−1)​+(−1)n(2n+1)!x(2n+1)​
3)
cos(x)=1−x22!+...+(−1)nx(2n)(2n)!+(−1)(n+1)cosθx(2n+2)!x(2n+2)cos(x)=1-\frac{x^2}{2!}+...+(-1)^n\frac{x^{(2n)}}{(2n)!}+(-1)^{(n+1)}\frac{cos\theta x}{(2n+2)!}x^{(2n+2)} cos(x)=1−2!x2​+...+(−1)n(2n)!x(2n)​+(−1)(n+1)(2n+2)!cosθx​x(2n+2)
4)
ln(1+x)=x−x22+...+(−1)(n−1)xnn+(−1)nx(n+1)(n+1)(1+θx)(n+1)ln(1+x)=x-\frac{x^2}{2}+...+(-1)^{(n-1)}\frac{x^n}{n}+(-1)^n\frac{x^{(n+1)}}{(n+1)(1+\theta x)^{(n+1)}} ln(1+x)=x−2x2​+...+(−1)(n−1)nxn​+(−1)n(n+1)(1+θx)(n+1)x(n+1)​
5)
(1+x)α=1+αx+α(α−1)2!x2+...+α(α−1)...(α−n+1)n!xn+α(α−1)...(α−n+1)(α−n)(n+1)!(1+θx)(α−n−1)x(n+1)(1+x)^\alpha = 1+\alpha x+\frac{\alpha(\alpha-1)}{2!}x^2+...+\frac{\alpha(\alpha-1)...(\alpha-n+1)}{n!}x^n+\frac{\alpha(\alpha-1)...(\alpha-n+1)(\alpha-n)}{(n+1)!}(1+\theta x)^{(\alpha-n-1)}x^{(n+1)} (1+x)α=1+αx+2!α(α−1)​x2+...+n!α(α−1)...(α−n+1)​xn+(n+1)!α(α−1)...(α−n+1)(α−n)​(1+θx)(α−n−1)x(n+1)
以上就是总结的几个常用的泰勒展开。

2.MATLAB实现方法

从上面的例子可以看出,麦克劳林公式实际上就是将f(x)f(x)f(x)表示为xnx^nxn的和的形式,而在MATLAB中实现泰勒展开的函数为taylor,其具体的语法格式如下所示:

实例:

(1)求ex2e^{x^2}ex2的7阶麦克劳林近似展开

syms x
f = exp(x^2);f7 = taylor(f)f7 =x^4/2 + x^2 + 1

(2)求sinxx\frac{sinx}{x}xsinx​的5阶麦克劳林型展开

syms x
f = sin(x)/x;
f5=taylor(f)f5 =x^4/120 - x^2/6 + 1

(3)求f(x)=asin(x)+bcos(x)+ctan(x)f(x)=asin(x)+bcos(x)+ctan(x)f(x)=asin(x)+bcos(x)+ctan(x)的10阶麦克劳林展开

syms x a b c
f = a*sin(x)+b*cos(x)+c*tan(x);
f10 = taylor(f,x,'order',10)f10 =(a/362880 + (62*c)/2835)*x^9 + (b*x^8)/40320 + ((17*c)/315 - a/5040)*x^7 - (b*x^6)/720 + (a/120 + (2*c)/15)*x^5 + (b*x^4)/24 + (c/3 - a/6)*x^3 - (b*x^2)/2 + (a + c)*x + b

如果想了解更多泰勒定理的应用,可以参考这篇文章:

https://zhuanlan.zhihu.com/p/355572556

MATLAB学习之泰勒展开(四)相关推荐

  1. matlab学习笔记(四)

    提示:文章写完后,目录可以自动生成,如何生成可参考右边的帮助文档 matlab学习笔记(四) 一.使用MATLAB实现函数的泰勒展开 傅里叶展开 提示:以下是本篇文章正文内容,下面案例可供参考 一.使 ...

  2. matlab中tunx_f是什么意思,MATLAB学习笔记(四)【更新ing】

    斑码教育web前端课程笔记-第十四天-CSS 斑码教育web前端课程笔记-第十四天-CSS 一.复习 定位:position static 默认值 没有定位在文档流 relative:相对定位 自己原 ...

  3. MATLAB学习笔记(十四)

    MATLAB学习笔记(十四) 一.线性方程组求解 1.1 直接法 1.1.1 利用左除运算符 1.1.2 利用矩阵分解 1.2 迭代法 1.2.1 雅可比(Jacobi)迭代法 1.2.2 高斯-赛德 ...

  4. 控制系统仿真与CAD-薛定宇-第四章matlab学习笔记

    控制系统仿真与CAD-薛定宇-第四章matlab学习笔记 04-02传递函数模型 tfdata() 传递函数属性法 04-07典型系统连接计算 pretty 用法 04-08方框图简化 04-09代数 ...

  5. HFSS学习笔记(四)利用MATLAB脚本进行HFSS快速建模

    HFSS学习笔记(四)利用MATLAB脚本进行HFSS快速建模 在HFSS中,我们经常会遇到复杂且重复的模型,比如你要画30个不同的矩形块,每个矩形块都有自己的参数比如,a1,a2,-,a30,那么此 ...

  6. OpenCV学习笔记(四十一)——再看基础数据结构core OpenCV学习笔记(四十二)——Mat数据操作之普通青年、文艺青年、暴力青年 OpenCV学习笔记(四十三)——存取像素值操作汇总co

    OpenCV学习笔记(四十一)--再看基础数据结构core 记得我在OpenCV学习笔记(四)--新版本的数据结构core里面讲过新版本的数据结构了,可是我再看这部分的时候,我发现我当时实在是看得太马 ...

  7. MATLAB学习笔记(一)

    MATLAB 一.安装的目录结构 二.常用的菜单及工具 布局 设置路径 设置工作路径 三.常用命令 四. 基础知识 数据类型 复数类型 `inf 和 NaN` 字符与字符串 函数句柄 结构体 结构体的 ...

  8. 集成学习(ensemble learning)(四)

    文章目录 一.GBDT概述 二.GDBT原理 1.负梯度拟合 2.GBDT回归算法 3.GBDT分类算法 (1)二元GBDT分类算法 (2)多元GBDT分类算法 4.GBDT常用损失函数 (1)分类任 ...

  9. linspace函数matlab_从零开始的matlab学习笔记——(29)泰勒逼近函数

    matlab应用--求极限,求导,求积分,解方程,概率统计,函数绘图,三维图像,拟合函数,动态图....更多内容尽在个人专栏:matlab学习 上一节我们成功制作了能自己转圈的三维螺旋线,这里我们再来 ...

  10. matlab数组平方的计算自定义函数_从零开始的matlab学习笔记——(38)简单数论计算函数:取整,gcd,lcm,质数,全排列...

    matlab应用--求极限,求导,求积分,解方程,概率统计,函数绘图,三维图像,拟合函数,动态图,傅里叶变换,随机数,优化问题....更多内容尽在个人专栏:matlab学习 翻了翻优化工具箱,发现内容 ...

最新文章

  1. java 基础(匿名内部类)
  2. jQuery 插入元素
  3. 删除归档日志_ORA-00257:archiver error.Connect internal only, until freed的问题
  4. C++学习-环境配置
  5. 游泳后精疲力尽_精疲力尽的编程后如何重回正轨
  6. java编译大项目很慢_优化MyEclipse编译速度慢的问题、build、project clean 慢
  7. 现实世界中正在用Java解决的难题
  8. 计算机基础(五):堆和栈区别
  9. 转:lnmp 搭建手册-黑一路人
  10. AliSQL 配置安装
  11. java -jar 详解,Java中jar命令详解
  12. __wakeup绕过版本_PHP__wakeup()方法
  13. Face Swapping: Automatically Replacing Faces in Photographs论文阅读
  14. python解决猴子偷桃问题_猴子偷桃蟠桃园土地是知道的,不举报是不敢吗?
  15. threw ‘java.lang.NullPointerException‘ exception // toString()
  16. 双十一苏宁“狮晚”PK天猫“猫晚”,玩的是娱乐,赢的是流量
  17. 在m1/m2芯片的mac电脑上运行Stable Diffusion的全步骤
  18. python输出语句print格式_print输出语句的格式是什么?
  19. Fabric区块链开发详解
  20. 执法部门的“新助手” 扫二维码进群反映问题

热门文章

  1. 分数混合运算简便方法_分数混合运算和简便算法
  2. Aerial Cactus Identification(空中仙人掌鉴定)
  3. office 2010 激活 toolkit.exe 应用程序正常初始化失败
  4. 《Kotlin从零到精通Android开发》面世啦
  5. 计算机学院教师老带新总结,班主任以老带新工作总结(共6篇)
  6. 公司部门英文缩写简称收集
  7. 织梦木马 data.php,DedeCMS后门木马专杀工具V2.0
  8. 程序员的日常工作是怎样的?
  9. 小公司个人技术成长路线思考
  10. 基于原生Wine6.0.2稳定版打补丁解决ukylin-wine图标显示异常和deepin-wine闪退的问题,Ubuntu运行QQ接近完美(附补丁)。