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）