import numpy as np
import matplotlib.pyplot as plt
import matplotlib as mpl
import math#e^x = 1 + x + x^2/2!+...
def calc_e_small(x):n = 10#累乘  cumsum是求和#1! 2! 3! 4! 5!...10!f = np.arange(1,n+1).cumprod()#x x^2 ... x^10b = np.array([x]*n).cumprod()return 1+np.sum(b / f)'''
e^x = ln2 + (e^ln2)/1!*(x-ln2) + (e^ln2)/2!*(x-ln2)^2+...
x = a*ln2 + b   k<= z  |b| <= 1/2ln2
a = ln( int( x/ln2 + 0.5 ) )
b = x-a*ln2
e^x = 2^a + e^b
'''
def calc_e(x):reverse = Falseif x < 0:#处理负数  exp(-x) = 1/exp(x)x = -xreverse = Trueln2 = 0.69314718055994530941723212145818c = x/ln2a = int(c+0.5)b = x-a*ln2#2的a次方乘以e的b次幂y = (2**a)*calc_e_small(b)if reverse:return 1/yreturn yif __name__ == '__main__':#-2到0 十个数t1 = np.linspace(-2,0,10,endpoint=False)#0到2 二十个数t2 = np.linspace(0,2,20)t = np.concatenate((t1,t2))print(t)#横轴数据y = np.empty_like(t)for i,x in enumerate(t):y[i] = calc_e(x)print('e^',x,'=',y[i],'(近似值)\t',math.exp(x))mpl.rcParams['font.sans-serif'] = [u'SimHei']mpl.rcParams['axes.unicode_minus'] = Falseplt.plot(t, y, 'r-', linewidth=2)plt.plot(t, y, 'go', linewidth=2)plt.title(u'Taylor展开式的应用', fontsize=18)plt.xlabel('X', fontsize=15)plt.ylabel('exp(X)', fontsize=15)plt.grid(True)plt.show()

转载于:https://www.cnblogs.com/xiaochi/p/11236023.html

python写泰勒展开式相关推荐

  1. python 编辑数学公式_Jupyter快速编辑高大上数学公式 泰勒展开式

    欢迎点击「算法与编程之美」↑关注我们! 本文首发于微信公众号:"算法与编程之美",欢迎关注,及时了解更多此系列博客. 人工智能的基础是机器学习,而通过之前的几篇博客了解到,机器学习 ...

  2. 【Python】SymPy库——求泰勒展开式+求极限 代码实现

    记录下泰勒展开式和极限的代码实现,多看官方文档,幸福长久远. 以如下函数为例: import sympy"""自定义函数,求其在x=0附近的泰勒展开(4项)"& ...

  3. python泰勒公式法求正弦函数_正弦函数两种泰勒展开式的比较

    正弦函数两种泰勒展开式的比较 张文华,汲守峰 [摘 要]摘要:讨论了正弦函数在两种不同情况下的泰勒公式展开式,并利用余项比较两种展开式 在近似计算中误差的大小区别,解释了正弦函数展开式中经常展开偶数项 ...

  4. 【学习笔记】超简单的多项式求指(含泰勒展开式、牛顿迭代完成证明)

    整理的算法模板合集: ACM模板 目录 多项式求指 泰勒展开式 牛顿迭代 牛顿迭代应用 P4726 [模板]多项式指数函数(多项式 exp) 数组版 (1.64s) vector版 (960 ms) ...

  5. 正弦函数泰勒展开c语言,C++ 学习笔记_0012_函数(泰勒展开式、三角函数表)

    泰勒展开式 先听故事,再编程序.故事是这样的:话说sin和cos是一对夫妇.一天,sin去听相声了,cos在家.过了一会,有人敲门,cos开门一看,是一个不认识的多项式函数.cos问:你是谁啊?他说: ...

  6. 数学之美—泰勒展开式

    1:背景 首先给大家介绍两位数学界泰斗:    麦克劳林,18世纪英国最具有影响的数学家之一.   他以熟练的几何方法和穷竭法论证了流数学说,还把级数作为求积分的方法,以几何形式给出了无穷级数收敛的积 ...

  7. sinx的泰勒展开式

    sinx的泰勒展开式求解过程 前些天发现了一个巨牛的人工智能学习网站,通俗易懂,风趣幽默,忍不住分享一下给大家,(ノ´▽`)ノ♪->点击这里->一个宝藏级人工智能教程网站. 思路: sin ...

  8. 一元和二元的泰勒展开式

    一元函数的泰勒展开式 由高等数学知识可知,对于一元函数f(x)f(x)f(x) 在kkk点,即x=x(k)x=x^{(k)}x=x(k)的泰勒展开式为: f(x)=f(x(k))+f′(x(k))(x ...

  9. 伯努利数、欧拉数与泰勒展开式

    伯努利数和正切函数泰勒展开式 1.伯努利数 2.欧拉数 3.自然幂指数和 4.tan⁡x\tan{x}tanx泰勒展开式 5. 程序应用与计算 参考 1.伯努利数 伯努利数是十八世纪瑞士数学家雅各布· ...

最新文章

  1. Android面试题详细整理系列(一)
  2. 【数据结构与算法】之深入解析“二叉树展开为链表”的求解思路与算法示例
  3. HDU 2177HDU 2176
  4. nginx和apache的伪静态区别
  5. Java的IO:BIO | NIO | AIO
  6. 大前端时代,从前端小工到架构师的进阶锦囊!
  7. 【教程】合成方法的那些事儿(附加练习源文件)
  8. Atitit table的读取解析from html by jsoup java html paresr 目录 1. atitit.gui控件的定位与选择器 1 2. 读取表格流程 遍历表格ta
  9. 关于无法卸载和安装VISIO2010的问题
  10. 考研高等数学张宇30讲笔记——第十一讲 多元函数微分学
  11. 运筹说 第32期 | 对偶理论与灵敏度分析—灵敏度分析
  12. 基于springboot旅游系统
  13. 如何缩减XFS分区格式的根目录
  14. Ant Design of Vue Icon 图标列表
  15. applicationContext.xml的头文件
  16. JS实现简易画图工具
  17. httpclient动态代理ip
  18. 基于javaweb的校园外卖点餐系统(java+ssm+jsp+mysql)
  19. 八年双11背后:新技术新能源成为新基础
  20. 数据仓库系列3-事实表

热门文章

  1. CSS学习之display属性与浮动
  2. ERP使用技巧:仓库分类管理五大误区
  3. 解答私信@田田WX //2020-12-30 左右n阶螺旋矩阵
  4. 时尚红蓝双色惊爆眼球,华为Mate 9玛瑙红/托帕蓝图赏
  5. Android 仿支付宝支付密码输入框的具体实现
  6. MongoDB副本集搭建、连接php
  7. 关于理想其实我和你一样都有
  8. 教育整改进行时,家长有话要说
  9. 历史大讲堂:iPhone为啥这么火?苹果系统历史回顾(下)
  10. 黑谍旺旺群发V3.11绿色注册版