Carrier Recovery Part 3: Phase-Locked Loop and Costas Loop

[@SaintSabel]

good stuff Liang, keep up the good work!

[@EvandsonDantas]

Professor Dong, i'm experiencing a phase error cause truncation. Is my matlab code correct?
t = 0:0.01:10; % Time second
Omega0 = 2*pi*2; % 2 Hz
Omega1 = 2*pi*3; % 3 Hz
theta0 = pi/4; % 45 degree
theta1 = pi/6; % 30 degree
y0 = cos(Omega0*t + theta0);
y1 = cos(Omega1*t + theta1);
plot(t, acos(y0) - acos(y1), t, (Omega0*t + theta0) - (Omega1*t + theta1));