import scipy.signal as sg import numpy as np import matplotlib.pyplot as plt sample_rate = 48000.0 nsamples = 320 F_1 = 440.0 F_2 = 10000.0 F_3 = 10000.0 nyq_rate = sample_rate / 2.0 cutoff_hz = 5000.0 numtaps = 29 t = np.arange(nsamples) / sample_rate vin = np.sin(2 * np.pi * F_1 * t) vfm = np.sin(2 * np.pi * F_2 * t + 6.0 * -np.cos(2 * np.pi * F_1 * t)) i1 = vfm * np.cos(2 * np.pi * F_3 * t) q1