[fase,error_fase_bpsk,Yn]=sinc_fase(In,Qn,fase_anterior,error_fase_bpsk_anterior,Yn_anterior);

function [fase,error_fase_bpsk,Yn]=sinc_fase(In,Qn,fase_anterior,error_fase_bpsk_anterior,Yn_anterior);


%INICIALIZACION DE VARIABLES
alfa=0.1;
beta=0.001;

%ERROR DE FASE BPSK
error_fase_bpsk=Qn./In;

%FILTRADO DE PRIMER ORDEN
Yn=(alfa+beta)*error_fase_bpsk-alfa*error_fase_bpsk_anterior+Yn_anterior;

%VCO
fase=Yn+fase_anterior;