Suppose that I have a transfer function model
TransferFunctionModel[{{{1.76073542174078*^7 + 15432.189447175271` s}}, 1.9007939847037695`*^7 + 17800.708617496595` s + 1. s^2}, s]
I want to calculate the group delay. So theoretically I need to differentiate the Bode phase diagrams with frequency. However, in this case, one should first unwrap the phase resulting from the Bode phase plot to get the correct result. I have seen other posts where the unwrapping phase in discrete data was discussed but could not find an appropriate method to apply directly on a transfer function (or maybe I missed something?).
How can one use Mathematica to calculate the "group delay" using an appropriate phase unwrapping method?