Python超越函數(shù)積分運(yùn)算以及繪圖實(shí)現(xiàn)代碼
編譯環(huán)境:ubuntu17.04 Python3.5
所需庫:numpy、scipy、matplotlib
下面是理想平面的輻射強(qiáng)度計(jì)算(課程大作業(yè)~~~)
1、超越函數(shù)積分運(yùn)算
def integral(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、繪圖實(shí)現(xiàn)
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='upper right',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='upper right',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這個(gè)函數(shù),這個(gè)是用來調(diào)整子圖之間的間距的啦
成果圖:

以上這篇Python超越函數(shù)積分運(yùn)算以及繪圖實(shí)現(xiàn)代碼就是小編分享給大家的全部?jī)?nèi)容了,希望能給大家一個(gè)參考,也希望大家多多支持腳本之家。
相關(guān)文章
在Python中字典按值排序的實(shí)現(xiàn)方法
這篇文章主要介紹了在Python中字典按值排序的實(shí)現(xiàn)方法,文中通過示例代碼介紹的非常詳細(xì),對(duì)大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值,需要的朋友們下面隨著小編來一起學(xué)習(xí)學(xué)習(xí)吧2020-11-11
ActiveMQ:使用Python訪問ActiveMQ的方法
今天小編就為大家分享一篇ActiveMQ:使用Python訪問ActiveMQ的方法,具有很好的參考價(jià)值,希望對(duì)大家有所幫助。一起跟隨小編過來看看吧2019-01-01
Python中scatter函數(shù)參數(shù)及用法詳解
這篇文章主要介紹了Python中scatter函數(shù)參數(shù)及用法詳解,小編覺得挺不錯(cuò)的,現(xiàn)在分享給大家,也給大家做個(gè)參考。一起跟隨小編過來看看吧2017-11-11
Python使用dict.fromkeys()快速生成一個(gè)字典示例
這篇文章主要介紹了Python使用dict.fromkeys()快速生成一個(gè)字典,結(jié)合實(shí)例形式分析了Python基于dict.fromkeys()生成字典的相關(guān)操作技巧,需要的朋友可以參考下2019-04-04
Python3模擬curl發(fā)送post請(qǐng)求操作示例
這篇文章主要介紹了Python3模擬curl發(fā)送post請(qǐng)求操作,結(jié)合實(shí)例形式分析了Python3使用Request請(qǐng)求模擬curl發(fā)送post相關(guān)操作技巧,需要的朋友可以參考下2019-05-05
matplotlib實(shí)現(xiàn)自定義散點(diǎn)形狀marker的3種方法
本文主要介紹了matplotlib實(shí)現(xiàn)自定義散點(diǎn)形狀marker的3種方法,文中通過示例代碼介紹的非常詳細(xì),具有一定的參考價(jià)值,感興趣的小伙伴們可以參考一下2021-10-10

