Qs = (Q, 1)
for k in range(2):
axs[k].plot(np.arange(T)/T, phi_f[0], "b")
for psi_f in psis_f[k]:
axs[k].plot(np.arange(T)/T, psi_f[0], "b")
axs[k].set_xlim(0, 0.5)
After Change
plt.xlim(0, 0.5)
plt.xlabel(r"$\omega$", fontsize=18)
plt.ylabel(r"$\hat\psi_j(\omega)$", fontsize=18)
plt.title("First-order filters (Q = {})".format(Q), fontsize=18)
//////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
// Do the same plot for the second-order filters. Note that since here `Q = 1`,
