这里用Python逼近函数y = exp(x);同样使用泰勒函数去逼近:
exp(x) = 1 + x + (x)^2/(2!) + .. + (x)^n/(n!) + ...
#!/usr/bin/python# -*- Coding:utf-8 -*-import numpy as npimport mathimport matplotlib as mplimport matplotlib.pyplot as pltdef calc_e_small(x): n = 10 f = np.arange(1,n+1).cumprod() b = np.array([x]*n).cumprod() return np.sum(b / f) + 1def 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 yif __name__ == "__main__": t1 = np.linspace(-2,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,'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实现画出e指数函数的图像所遇到的程序开发问题。
如果觉得内存溢出网站内容还不错,欢迎将内存溢出网站推荐给程序员好友。
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)