今天小编就为大家分享一篇大蟒实现画出e指数函数的图像,具有很好的参考价值,希望对大家有所帮助。一起跟随小编过来看看吧
这里用计算机编程语言逼近函数y=exp(x);同样使用泰勒函数去逼近:
exp(x)=1 x (x)^2/(2!).(x)^n/(n!).
#!/usr/bin/python
# -*-编码:utf-8 -*-
将numpy作为铭牌导入
导入数学
将matplotlib作为数学编程语言(Mathematical Programming Language)导入
将matplotlib.pyplot作为plt导入
def calc_e_small(x):
n=10
f=np.arange(1,n ^ 1).累计产量()
b=np.array([x]*n).累计产量()
return np.sum(b/f) 1
定义计算值(十):
反向=错误
如果x 0: #处理负数
x=-x
反向=真
ln2=
c=x/ln2
a=int(c 0.5)
b=x - a*ln2
y=(2 ** a) * calc_e_small(b)
如果相反:
返回1/年
回车y
if __name__=='__main__ ':
t1=np.linspace(-2,0,10,endpoint=False)
t2=np.linspace(0,3,20)
t=np.concatenate((t1,t2))
打印(t) #横轴数据
y=np.empty_like(t)
对于我,x在列举(t):
y[i]=calc_e(x)
print('e^',x,'=',y[i],'(近似值)\t ',math.exp(x),'(真实值)')
#打印'误差:',y[i] - math.exp(x)
plt.figure(facecolor='w ')
mpl。RC params[' font。sans-serif ']=[u ' sim hei ']
mpl。RC params['轴。unicode _ MINUS ']=False
plt.plot(t,y,' r-',t,y,' go ',线宽=2)
工厂名称(泰勒展式的应用-指数函数,fontsize=18)
plt.xlabel('X ',fontsize=15)
plt.ylabel('exp(X)',fontsize=15)
工厂网格(真)
plt.show()
以上这篇大蟒实现画出e指数函数的图像就是小编分享给大家的全部内容了,希望能给大家一个参考,也希望大家多多支持我们。
郑重声明:本文由网友发布,不代表盛行IT的观点,版权归原作者所有,仅为传播更多信息之目的,如有侵权请联系,我们将第一时间修改或删除,多谢。