Vc=5; % PLL supply voltage
w0=2*pi*100; % VCO center frequency
kw=w0/Vc; % VCO frequency gain
nF=[1]; % LPF numerator
dF=[1 1]; % LPF denumerator

T=1e-4; % sampling time
S=100; % stop time
N=floor(S/T);
k=0:N;
tt=k*T;
qq=(4*pi+pi/4)/S*tt-(2*pi+pi/4);
p=[tt' qq'];

sim('pll')

% phase comparator characteristics
figure
hold on
grid on
axis([-360 360 -5 5])
plot(q*180/pi,v0,'r')
plot(q*180/pi,-2*Vc/pi*sin(q-pi/2),'r:')
plot(q*180/pi,v1,'b')
plot([-2*pi -pi 0 pi 2*pi]*180/pi,[Vc -Vc Vc -Vc Vc],'b:')
plot(q*180/pi,v2,'g')
plot([-2*pi -pi 0 pi 2*pi]*180/pi,[-Vc Vc -Vc Vc -Vc],'g:')
plot(q*180/pi,v3,'y')
plot([-2*pi 0 0 2*pi]*180/pi,[-Vc Vc -Vc Vc],'y:')
plot(q*180/pi,v4,'c')
plot([-2*pi 2*pi]*180/pi,[-Vc Vc],'c:')