load TF.mat

Ts = 1/20000;

TF = tf(num,den,Ts)

freq=[0.01:1:2*pi*(0.5/Ts)];

bode(TF,freq)
