How can I plot amplitude and phase responses of a CIC filter?

I'm trying to plot amplitude and phase responses of a CIC filter (Cascaded integrated Comb filters), given M,R and N. Here is the z transfer function.

https://en.wikipedia.org/wiki/Cascaded_integrator%E2%80%93comb_filter

Answer 1

[Edson Silva]

CIC1 = dsp.CICDecimator(8,1,3);
h=fvtool(CIC1,'Fs',sampling_frequency);
h.NormalizeMagnitudeto1 = 'on';