返回具有求和项的复杂函数

Question

我目前正在尝试对某些原始数据应用经过改进的FFT，但是我要么在数据中出现振荡，要么出现了不应该发生的完美直线。参见下面的等式。

Equation

这是我正在尝试的。

def fft(w, j0 ,eta, time, creep):
def k_th(k):
    return (np.exp(-1j*w*time[k-1])-np.exp(-1j*w*time[k]))*(creep[k]-creep[k-1])/(time[k]-time[k-1])

g_complex = 1j*w/((1j*w*j0+(1-np.exp(-1j*w*time[0]))*((creep[0]-j0)/time[0])+(np.exp(-1j*w*time[-1])/eta)) + np.array([k_th(k) for k in range(1,len(time))]).sum())

return g_complex

w=np.arange(0.1, 50) #omega range
j0 = np.array(jt[1:])[0] #const
t2 = np.array(jtime[1:]) #list to iter
c2 = np.array(jt[1:]) # list to iter

g_comp = fft(w, j0, eta, t2, c2)

我认为总和项的计算不正确。任何帮助，将不胜感激。