Needs["Graphics`Graphics`"]

pole[f_, f0_, q0_] := 
      (1/(1 - (((f/f0))^2) + I f/f0/q0))

zero[f_, f0_, q0_] := 1 - (((f/f0))^2) + I f/f0/q0

H = 
      pole[f, 31.8, 1.41] 
        zero[f, 424, Infinity]*
pole[f, 73.9, 5.83] 
        zero[f, 168, Infinity]*
pole[f, 93.2, 25.1] 
        zero[f, 133.4, Infinity]

LogLogPlot[Abs[H], {f, 1, 1000}, 
      PlotRange -> {0.0001, 3}]

LogLogPlot[Abs[H], {f, 100, 1000}, 
      PlotRange -> {0.0001, 0.001}]

LogLinearPlot[Arg[H]/π*180, {f, 1, 1000}, 
      PlotRange -> {-180, 180}]