Возврат сложной функции с суммированием - PullRequest
Новая
       8

Возврат сложной функции с суммированием

0 голосов
/ 28 февраля 2019

В настоящее время я пытаюсь применить модифицированный БПФ к некоторым необработанным данным, но я либо получаю колебания в данных, либо совершенно прямую линию, которая не должна происходить.Смотрите уравнение ниже.

Уравнение

Вот то, что я пытаюсь. 

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)

Я думаю, что сумма сумма не рассчитывается правильно.Любая помощь будет оценена.

...