Descomposición de una secuencia discreta en una combinación lineal de impulsos desplazados ponderados.

import numpy as np
import pylab as pl
from impulse_step_ramp import *

#-----------------------------------------------------
def udelta(n):
    y=(n==0)*1 #Convert bool to int [...*1]
    #i=np.where(y==1)
    return y
#-----------------------------------------------------

n=np.arange(-10,10)
s=udelta(n+6)-3*udelta(n+5)-4*udelta(n+4)+1*udelta(n+3)-5*udelta(n+2)+7*udelta(n+1)-1*udelta(n)+\
3*udelta(n-1)+3*udelta(n-2)-5*udelta(n-3)+2*udelta(n-4)-1*udelta(n-5)-1*udelta(n-6)
pl.stem(n,s)
pl.axis([-10,10,-10,10])
pl.grid(True)
pl.xlabel('Tiempo [n]')
for i, txt in enumerate(n):
    pl.annotate(s[i], (n[i],s[i]+0.5))
pl.show

No hay comentarios.:

Publicar un comentario