%transmitancja, 1/((s+5)(s+8))
%A=mod(K(jw))
%w=5
%
w=5;
mon=[1,13,40]
sys=tf(3,mon)
nyquist(sys)
Tw=(3/(65i+15))
Fi=angle(Tw)
A=abs(Tw)
a = 10;
t = 0 : 0.01 : 10;
B = A*cos(pi/2 - w*t +Fi);
u = cos(pi/2 -w*t)
plot(t,B);
hold on;
[Y_s,T_s,X_s] = lsim(sys,u,t)
plot(t,Y_s)
hold off;