Listing 10.5. costasloop.m simulate the Costas loop with input from pulrecsig.m r= r s c ; fl=500, ff=10 .01 .02 1]; fa=[1 1 0 0]; h-firpm (fl,ff , fa ) mu-.003 f0 1000:; theta-zeros (1 , length(t)); theta (1)=0; zs=zeros ( 1,f1+1); zc=zeros (1, fl+1); for k-1length (t)-1 % rse from pulreesig.m %LPF design % algorithm stepsize % freq at receiver % estimate vector % buffers for LPFs % z contains past inputs zs [zs (2: f 1), 2 r (k)*sin 2* pi Eo*t (k)+theta (k)) zc zc (2: f1), 2*r (k)*cos (2* pi fo*t (k)+theta (k))] lpfszfl pl r ( h ) * zs,; 1pfc efl ipl r (h)*zc"; % output of theta (k+1)= t heta (k)-m u* I pfs *lpfc; filters % algorithm update end Listing 10.5. costasloop.m simulate the Costas loop with input from pulrecsig.m r= r s c ; fl=500, ff=10 .01 .02 1]; fa=[1 1 0 0]; h-firpm (fl,ff , fa ) mu-.003 f0 1000:; theta-zeros (1 , length(t)); theta (1)=0; zs=zeros ( 1,f1+1); zc=zeros (1, fl+1); for k-1length (t)-1 % rse from pulreesig.m %LPF design % algorithm stepsize % freq at receiver % estimate vector % buffers for LPFs % z contains past inputs zs [zs (2: f 1), 2 r (k)*sin 2* pi Eo*t (k)+theta (k)) zc zc (2: f1), 2*r (k)*cos (2* pi fo*t (k)+theta (k))] lpfszfl pl r ( h ) * zs,; 1pfc efl ipl r (h)*zc"; % output of theta (k+1)= t heta (k)-m u* I pfs *lpfc; filters % algorithm update end