python实现画出e指数函数的图像


Posted in Python onNovember 21, 2019

这里用Python逼近函数y = exp(x);同样使用泰勒函数去逼近:

exp(x) = 1 + x + (x)^2/(2!) + .. + (x)^n/(n!) + ...

#!/usr/bin/python
# -*- coding:utf-8 -*-
 
import numpy as np
import math
import matplotlib as mpl
import matplotlib.pyplot as plt
 
 
def calc_e_small(x):
 n = 10
 f = np.arange(1, n+1).cumprod()
 b = np.array([x]*n).cumprod()
 return np.sum(b / f) + 1
 
 
def calc_e(x):
 reverse = False
 if x < 0: # 处理负数
  x = -x
  reverse = True
 ln2 = 0.69314718055994530941723212145818
 c = x / ln2
 a = int(c+0.5)
 b = x - a*ln2
 y = (2 ** a) * calc_e_small(b)
 if reverse:
  return 1/y
 return y
 
 
if __name__ == "__main__":
 t1 = np.linspace(-2, 0, 10, endpoint=False)
 t2 = np.linspace(0, 3, 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), '(真实值)')
  # print '误差:', y[i] - math.exp(x)
 plt.figure(facecolor='w')
 mpl.rcParams['font.sans-serif'] = [u'SimHei']
 mpl.rcParams['axes.unicode_minus'] = False
 plt.plot(t, y, 'r-', 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()

python实现画出e指数函数的图像

以上这篇python实现画出e指数函数的图像就是小编分享给大家的全部内容了,希望能给大家一个参考,也希望大家多多支持三水点靠木。

Python 相关文章推荐
用Python写的图片蜘蛛人代码
Aug 27 Python
python中while循环语句用法简单实例
May 07 Python
Python连接SQLServer2000的方法详解
Apr 19 Python
如何用itertools解决无序排列组合的问题
May 18 Python
python利用跳板机ssh远程连接redis的方法
Feb 19 Python
PyQt5响应回车事件的方法
Jun 25 Python
Python实现一个数组除以一个数的例子
Jul 20 Python
python单线程下实现多个socket并发过程详解
Jul 27 Python
详解如何从TensorFlow的mnist数据集导出手写体数字图片
Aug 05 Python
Python传递参数的多种方式(小结)
Sep 18 Python
Django分页器的用法你都了解吗
May 26 Python
Python面向对象之成员相关知识总结
Jun 24 Python
如何获取Python简单for循环索引
Nov 21 #Python
Python Web静态服务器非堵塞模式实现方法示例
Nov 21 #Python
Windows10下Tensorflow2.0 安装及环境配置教程(图文)
Nov 21 #Python
使用python代码进行身份证号校验的实现示例
Nov 21 #Python
Python 面向对象之封装、继承、多态操作实例分析
Nov 21 #Python
用python画一只可爱的皮卡丘实例
Nov 21 #Python
Python 私有化操作实例分析
Nov 21 #Python
You might like
PHP无限分类的类
2007/01/02 PHP
深入PHP curl参数的详解
2013/06/17 PHP
两种php去除二维数组的重复项方法
2015/11/04 PHP
PHP实现将优酷土豆腾讯视频html地址转换成flash swf地址的方法
2017/08/04 PHP
yii 框架实现按天,月,年,自定义时间段统计数据的方法分析
2020/04/04 PHP
JQuery实现自定义对话框的代码
2008/06/15 Javascript
原创javascript小游戏实现代码
2010/08/19 Javascript
js完美实现@提到好友特效(兼容各大浏览器)
2015/03/16 Javascript
jQuery+HTML5加入购物车代码分享
2020/10/29 Javascript
javascript表单事件处理方法详解
2016/05/15 Javascript
喜大普奔!jQuery发布 3.0 最终版
2016/06/12 Javascript
微信公众号开发 实现点击返回按钮就返回到聊天界面
2016/12/15 Javascript
vue2组件之select2调用的示例代码
2017/10/12 Javascript
开发Vue树形组件的示例代码
2017/12/21 Javascript
JS实现的抛物线运动效果示例
2018/01/30 Javascript
解决vue中监听input只能输入数字及英文或者其他情况的问题
2018/08/30 Javascript
详解基于iview-ui的导航栏路径(面包屑)配置
2019/02/22 Javascript
js 实现 list转换成tree的方法示例(数组到树)
2019/08/18 Javascript
Vue源码分析之Vue实例初始化详解
2019/08/25 Javascript
关于JS模块化的知识点分享
2019/10/16 Javascript
vue2和vue3的v-if与v-for优先级对比学习
2020/10/10 Javascript
详解Vue2的diff算法
2021/01/06 Vue.js
详解微信小程序(Taro)手动埋点和自动埋点的实现
2021/03/02 Javascript
Python实现的当前时间多加一天、一小时、一分钟操作示例
2018/05/21 Python
Python实现随机创建电话号码的方法示例
2018/12/07 Python
Python自动抢红包教程详解
2019/06/11 Python
pytorch实现Tensor变量之间的转换
2020/02/17 Python
python继承threading.Thread实现有返回值的子类实例
2020/05/02 Python
Numpy实现卷积神经网络(CNN)的示例
2020/10/09 Python
HTML5页面无缝闪开的问题及解决方案
2020/06/11 HTML / CSS
北京银河万佳Java面试题
2012/03/21 面试题
Java里面Pass by value和Pass by Reference是什么意思
2016/05/02 面试题
党的群众路线教育实践活动心得体会
2014/03/03 职场文书
教师一岗双责责任书
2014/04/16 职场文书
电教室标语
2014/06/20 职场文书
委托开发合同书(标准版)
2019/08/07 职场文书