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

yipeiwu_com5年前Python基础

这里用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设计】。

相关文章

VSCode中自动为Python文件添加头部注释

VSCode中自动为Python文件添加头部注释

在实际编写Python文件时,往往需要为文件添加相关说明,例如文件名称、文件作用、创建时间、作者信息、版本号等等。这些信息往往是固定模板的,因此希望有一种方式可以自动的为我们添加上这些信...

Python实现将HTML转换成doc格式文件的方法示例

本文实例讲述了Python实现将HTML转换成doc格式文件的方法。分享给大家供大家参考,具体如下: 网页上的一些文章,因为有格式的原因,它们在网页上的源码都是带有html标签的,用cs...

Python使用functools模块中的partial函数生成偏函数

python 中提供一种用于对函数固定属性的函数(与数学上的偏函数不一样) # 通常会返回10进制 int('12345') # print 12345 # 使用参数 返回 8...

django定期执行任务(实例讲解)

要在django项目中定期执行任务,比如每天一定的时间点抓取数据,刷新数据库等,可以参考stackoverflow的方法,先编写一个manage.py命令,然后使用crontab来定时执...

python音频处理用到的操作的示例代码

python音频处理用到的操作的示例代码

前言 本文主要记录python下音频常用的操作,以.wav格式文件为例。其实网上有很多现成的音频工具包,如果仅仅调用,工具包是更方便的。 更多pyton下的操作可以参考: 用python...