Factor of 2pi in tfAAA output

Created by: jeffwack111

The following script creates a simple system and tries to recover it using tfAAA. After plotting the frequency response of both systems, they disagree by a factor of 2pi. I used f and F_Hz in all function calls.

from wield.control import SISO, AAA
import matplotlib.pyplot as plt

P = SISO.zpk([],[-10*2*np.pi],1)
f_Hz = np.geomspace(1,1000,500)

originalresponse = P.fresponse(f = f_Hz)

bary = AAA.tfAAA(F_Hz=f_Hz,xfer=originalresponse.mag)
newresponse = SISO.zpk(bary.zpk).fresponse(f=f_Hz)

plt.loglog(f_Hz,originalresponse.mag,label = "original")
plt.loglog(f_Hz,newresponse.mag,label = "recovered")
plt.legend()