Python超越函数积分运算以及绘图实现代码
编译环境：ubuntu17.04Python3.5
所需库：numpy、scipy、matplotlib
下面是理想平面的辐射强度计算（课程大作业～～～）
1、超越函数积分运算
defintegral(x,c1,c2,T): return((c1*0.98)/(x**5))*(1/((np.e**(c2/(x*T)))-1)) resut,err=integrate.quad(integral,3,5,args=(c1,c2,T))
2、绘图实现
plt.figure(1) ax1=plt.subplot(211) plt.sca(ax1) plt.plot(fi,functionI(fi,0.5,5,1,e0),label='n=5,ks=0.5') plt.legend(loc='upperright',bbox_to_anchor=(0.9,0.9)) plt.xlabel(u'ψ/rad') plt.ylabel(u'I/(W/sr)') ax2=plt.subplot(212) plt.sca(ax2) plt.plot(fi,functionI(fi,0.5,5,1,e0),label='n=5,ks=0.5') plt.legend(loc='upperright',bbox_to_anchor=(0.9,0.9)) plt.xlabel(u'ψ/rad') plt.ylabel(u'I/(W/sr)') plt.subplots_adjust(wspace=0.5,hspace=0.5) plt.show()
说一下plt.subplots_adjust这个函数，这个是用来调整子图之间的间距的啦
成果图：
以上这篇Python超越函数积分运算以及绘图实现代码就是小编分享给大家的全部内容了，希望能给大家一个参考，也希望大家多多支持毛票票。
声明：本文内容来源于网络，版权归原作者所有，内容由互联网用户自发贡献自行上传，本网站不拥有所有权，未作人工编辑处理，也不承担相关法律责任。如果您发现有涉嫌版权的内容，欢迎发送邮件至：czq8825#qq.com（发邮件时，请将#更换为@）进行举报，并提供相关证据，一经查实，本站将立刻删除涉嫌侵权内容。