If you are using the Symbolic Math Toolbox to get the Laplace domain transfer function, the sym2poly function will give you the coefficients. Use the freqs function (with a continuous filter) to get the Bode plot (passband and phase response of the filter). You will have to experiment to get the result you want. See the documentation for the various functions.
With the filtfilt function all discrete filters have a maximally flat phase response, so phase distortion is irrelevant if you use it.
