Serie de Fourier de una onda cuadrada.

import pylab
import numpy as np
#Con SINC
k=np.arange(-30,30)
ak=(1/8.0)*np.sinc(k/8.0)*np.exp(-1j*k*np.pi/8.0)
m=np.abs(ak)
f=np.angle(ak)
#pylab.figure(2)
#pylab.stem(k,m,markerfmt=" ")
#pylab.show()
#Reconstruccion
t=np.linspace(-2,2,10000)
s_per=0
for k in range(-30,30):     
    #print s_per
    s_per=s_per+(1/8.0)*np.sinc(k/8.0)*np.exp(-1j*k*np.pi/8.0)*np.exp(1j*k*2*np.pi*t)    
pylab.plot(t,s_per), pylab.show()