PPL锁相环matlab

clear all
clc
fo=32;
fi=27;
% fi=length（ui）*1/（n*fo）;
g1=0.005*2;
g2=0.0975*2;
n=64;
t=0:1/（n*fo）:0.5-1/（n*fo）;
ui=sin（2*pi*fi*t）;%输入信号
uo=zeros（1length（ui））;
uo（1）=dds_sin（10n*fofo01）;%Nco输出初始信号
uo1（1）=dds_sin（10n*fofo01）;
% erro=ph（ui（1）uo（1））;%鉴相器输出
sig（1）=ui（1）*uo（1）;
l=30;
hn=fir1（l0.2*fo/（n*fo））;
temp=conv（hnsig（1））;
temp1=temp（2:end）;
erro（1）=temp（16）;
lpu1（1）=erro（1）*（g1+g2）;
lpu2（1）=lpu1（1）