clc,clear all;
%%
y=[0 0 3.045767 4.0824 -0.91515 -7.9588 -10.4576 -12.0412...
    -14.4249 -17.7211 -18.4164 -20];
x=[3.141593 6.283185 9.424778 12.56637 15.70796 18.84956...
    21.99115 25.13274 28.27433 31.41593 34.55752 37.69911];

figure(1)
subplot(2,1,1);
semilogx(x,y);
xlim([4*1e0 4*1e1]);
ylim([-20 5]);
grid on;
hold on;
%%
tv=1/6.238;
t=1/12.57;

s = tf('s');
H = ((s*tv+1)^2)/((s*t + 1)^5)

P = bodeoptions; 
P.PhaseVisible = 'off';
P.XLim=[4*1e0 4*1e1];
P.YLim=[-20 5];
subplot(2,1,2);
h = bodeplot(H,P); 
grid on;