Este es el código
T = 1 #Se normaliza respecto de la frecuencia de muestreo
Wo = 1 / float(2 * T) #Frecuencia de Nyquist
taps = 65 #Numero de coeficientes del filtro
K = taps / 2
step = 1 / float(K * T)
f = np.arange(-1 , 1 , step)
alfa = 0.5 #Factor de despliegue
W = alfa * Wo + Wo
H = sample_sqrtHrc(f, W, Wo)
hc = np.fft.ifft(H)
h = np.real(hc)
L=len(h)
# esta linea de codigo es el causante necesito ayuda, que modifico
b = np.concatenate((h[L-K:L], h[0:K]), axis=None)