I need to calculate the poles and zeros from minimum 4th order (max ~20) filter coefficients. I have limitations that prohibit the companion matrix/eigen decomposition method used by roots() in Matlab.
I can see many methods for general iterative root solving, but would really appreciate some advice whether there are any methods that work well for digital filters.