def sinc(n, radius=3, freq=1.0): taps = np.linspace(-radius,radius, n) return freq * np.sinc(freq*taps)