pin_pian.m

此为频偏校正程序,通过TD_GG算法得到的频偏值还存在几赫兹的残余频偏

function  [sum,yk,xk,wk]=pin_pian(I1,Q1)


M=8;
%M=16;%%%%%%或16次方
sum=0;
fai=atan2(Q1,I1);
xk=exp(j*fai);
yk=exp(j*M*fai);
N=length(yk);
for k=2:length(yk)
    yk1=yk(k)*conj(yk(k-1));
    wk=6*(k-1)*(N-(k-1))/(M*N*(N*N-1));
    sum=sum+angle(yk1)*wk;
end
return