b, a = scipy.signal.butter(2, np.array(lfreq)/(sampling_rate/2.), btype = 'low')