import matplotlib.pyplot as plt
import numpy as np
from sk_dsp_comm.sigsys import peaking
from scipy import signal
b,a = peaking(2.0,500)
f = np.logspace(1,5,400)
w,H = signal.freqz(b,a,2*np.pi*f/44100)
plt.semilogx(f,20*np.log10(abs(H)))
plt.ylabel("Power Spectral Density (dB)")
plt.xlabel("Frequency (Hz)")
plt.show()
#
b,a = peaking(-5.0,500,4)
w,H = signal.freqz(b,a,2*np.pi*f/44100)
plt.semilogx(f,20*np.log10(abs(H)))
plt.ylabel("Power Spectral Density (dB)")
plt.xlabel("Frequency (Hz)")