假設(shè)能譜已經(jīng)算出來(lái)，那么算格林函數(shù)是非常容易的

可以直接用向量化編程實(shí)現(xiàn)，這樣的計(jì)算效率最高。

nk = 1024 nw = 1024 delta = 0.01 ks = np.linspace(-pi,pi,nk) omega = np.linspace(-2,2,nw) Ks, Omega = np.meshgrid(ks,omega) Ek = -2*np.cos(Ks) Gk = np.zeros((nw,nk),dtype='double') Gk = np.power((Omega+1j*delta-Ek),-1) Ak = -Gk.imag/pi plt.pcolormesh(Ks, Omega, Ak) plt.colorbar()

這里積分只對(duì)k進(jìn)行，我們令\omga = 0，這個(gè)公式就可以寫(xiě)成兩個(gè)格林函數(shù)的卷積

如何用卷積的方式解決這個(gè)計(jì)算。

chi0 = signal.fftconvolve(Gk,Gk,'same')/nkx

?

總結(jié)

以上是生活随笔為你收集整理的格林函数数值计算的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。