G = 1/（z-H(k))

其中H(k)一個(gè)方向作傅里葉變換，另一個(gè)方向不做， 直接對(duì)矩陣進(jìn)行求逆。

Matlab代碼

nw = 1000; delta = 0.002; omega = linspace(-1,1,nw); A = zeros(nw,nk); for i=1:nwfor j=1:nkH0 = Hamiltonian_H0(k(j),t,N);H = H0+U*Nupavg-0.5*U*eye(m,m);G = ((omega(i)+1j*delta)*eye(m,m)-H)^-1;A(i,j) = -imag(trace(G));end end [X,Y] = meshgrid(k,omega); mesh(X,Y,A);

測(cè)試通過

總結(jié)

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

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