syms T;
f=1:1:1000;
T=5e-3;%delay time is 5e-3
Gdelay=cos(2*pi*f*T)+i*sin(2*pi*f*T);%explore it into the form of complex
subplot(2,1,1),plot(f,abs(Gdelay)),subplot(2,1,2),plot(f,angle(Gdelay)*180/pi)
