freqresp returns the phase already returns the phase in units of radians.
Computation the complex response should be:
sjw = (m * np.exp(1j * p)).transpose(2, 0, 1)