G=tf(8,[1 3 2]);
H=tf(1,[1 3]);
C=tf(0.987,[.267^2 2*.267 1]);
T=G*C+H;

abs(freqresp(T,4))