Обратное преобразование Фурье

###############ОБРАТНОЕ БПФ ДЛЯ 50 Герц################
# Определяем функцию нахождения 50 Гц гармоники
def garmonika(In_Mass):   
# находим спектр                                                    
    Spec=fft(In_Mass)    
# длина входного массива для итерации                                                     
    M=len(In_Mass)   
# Берем две первые составляющие массива постоянная +1 гармоника , на месте остальных пишем 0.                                                                        
    Filt_spectrum = [Spec[i] if i >= 0 and i <= 1 else 0 for i in range(M)]   
# делаем ОБПФ отфильтрованного спектра по всей длине массива. 
    Filt_signal = ifft(Filt_spectrum, n=M)      
# в качестве результата функции получаем массив фильтрованного сигнала                              
    return Filt_signal                                                        
 
##############Фильтруем Ia###########################  
  
Filt_signal1 = garmonika(Ia)
        ###строим график для Ia###
fig = plt1.figure() 
plt1.plot(arange(M),Filt_signal1)
plt1.xlabel('Time, ms')
plt1.ylabel('Ia, Ampere')
plt1.title('50 Hz harmonic')
plt1.grid(True)

N=48
Ak=list()
for k in range(0,N):
    for i in range(0,N):
        Ak_cycle=(1/N)*np.sum(Ia[i]*np.cos(2*np.pi*k*i/N))
        Ak.append(Ak_cycle)
  
Bk=list()
for k in range(0,N):
    for i in range(0,N):
        Bk_cycle=(-1/N)*np.sum(Ia[i]*np.sin(2*np.pi*k*i/N))
        Bk.append(Bk_cycle)
 
InvF=list()        
for n in range(0,N):
    for k in range(0,N):
        InvF_cycle=np.sum(Ak[k]*np.cos(2*np.pi*k*n/N)-Bk[k]*np.sin(2*np.pi*k*n/N))
        InvF.append(InvF_cycle)
  
plt1.plot(arange(0,len(InvF)),InvF)
plt1.xlabel('Time, ms')
plt1.ylabel('Ia, Ampere')
plt1.title('50 Hz harmonic')
plt1.grid(True) 
 
 
plt1.show()

N=len(Ia)
Ak=list()
for k in range(0,N):
    for i in range(0,N):
        Ak_cycle=(1/N)*np.sum(Ia[i]*np.cos(2*np.pi*k*i/N))
        Ak.append(Ak_cycle)
  
Bk=list()
for k in range(0,N):
    for i in range(0,N):
        Bk_cycle=(-1/N)*np.sum(Ia[i]*np.sin(2*np.pi*k*i/N))
        Bk.append(Bk_cycle)
        
Ck=list()
for k in range(0,N):
    Ck_cycle=Ak[k]+Bk[k]*1j
    Ck.append(Ck_cycle)
 
  
plt1.plot(arange(0,len(Ck)),Ck)
plt1.xlabel('Time, ms')
plt1.ylabel('Ia, Ampere')
plt1.title('50 Hz harmonic')
plt1.grid(True) 
 
 
plt1.show()