bias.py (386B)
1 import numpy as np 2 import matplotlib.pyplot as plt 3 4 fs = 48000 * 16 5 6 t = np.linspace (0, 1, fs) 7 8 f = 2000 9 I_in = np.sin (2 * np.pi * f * t) 10 11 f_bias = 50000 12 A_bias = 5 13 I_bias = A_bias * np.sin (2 * np.pi * f_bias * t) 14 15 I_head = I_in + I_bias 16 17 plt.plot (t[0:2000], I_head[0:2000]) 18 plt.title ("2 kHz Sine Wave Biased at 50 kHz") 19 plt.xlabel ("Time [s]") 20 plt.ylabel ("Amplitude") 21 plt.show()