Estimating transfer function from frd data

[ymohammadi]

I'm new to python and trying to find a way to fit a transfer function to FRF data. Something similar to sys = tfest(data,np,nz) in matlab where the inputs are the FRF data (in frd form), number of poles (np) and zeros (nz). The output contains the numerator and denominator coefficients of the transfer function.

I checked control library but couldn't find such a function.

Thanks for your help.