matlab仿真的载波同步算法-Fitz算法

matlab仿真的载波同步算法-Fitz算法，性能很好，接近MCRB，可广泛应用于卫星通信等行业。

clc;clear all;
bit_number=200;
P=32; 
D=68;
B=2;
d=P/B;
p=D/B;
a=1/2;
M=4;
symbol_number=bit_number/log2（M）;  
pi=3.14159;      
id=sqrt（-1）;
T=1/10^（4）;
fs=1/T;
n_t=[-3 3];
v=100;

%-------------------------------------
% add noise and average 
snr=-8:1:8;  
times=2000;     % average times
sigma=zeros（1length（snr））;
ve1=zeros（timeslength（snr））;

avve1=zeros（1length（snr））;

var1=zeros（1length（snr））;
for t1=1:length（snr）        
  sigma（t1）=sqrt（1/（2*10^（snr（t1）/10）））; 
 for k=1:times
     
   Data_input=randint（12*D）; %生成BPSK信号。。
   
   [Data]=qpsk（Data_input）;    %  qpsk调制
   
   Pilot_input=zeros（12*P）; 
   [Pilot]=qpsk（Pilot_input）; 
   
   %ck=[Pilot（1:P/2） Data Pilot（P/2+1:P）];
   ck=[Pilot Data ];

for i=1:length（ck）
    r（i）=ck（i）*exp（id*（2*pi*v*i/fs））;
end
% add noise
y=r+sigma（t1）*（randn（1length（r））+id*randn（1length（r）））;
% choose point from y1
% calculate z
z=zeros（1length（y））;%一般赋值时先给出空间在赋值。。。
for i=1:length（y）
    z（i）=y（i）*conj（ck（i））;
end
%Ip=[1:16 85:100];%前置-后置
Ip=[1:32];%前置
% Fitz  mothed
for m=1:P-1
    sum0=0;
    for k=m:P
       sum0=sum0+z（Ip（k））*conj（z（Ip（k）-m+1））;
    end
    RM（m）=sum0/（P-m）;
end
sum1=0;
for i=1:P-1
    sum1=sum1+angle（RM（i））;
end
f_est=sum1/（pi*（P-1）*P*T）; 
var1（t1）=var1（t1）+（（f_est-v）*T）^2/times;

%  Kay methed
% gama=zeros（1（L0-1））;
% for i=1:L0-1
%      gama（i）=3/2*L0/（L0^2-1）*（1-（（2*i-L0）/L0）^2）;
% end
% 
% for i=2:L0 %不懂。。
%     ve1（（t1-1）*times+k）=ve1（（t1-1）*times+k）+gama（i-1）*angle（z（i）*conj（z（i-1）））;
%  end
% ve1（（t1-1）*times+k）=ve1（（t1-1）*times+k）/（2*pi*T）;
% avve1（t1）=avve1（t1）+ve1（（t1-1）*times+k）/times;
% var1（t1）=var1（t1）+（（ve1（（t1-1）*times+k）-v）*T）^2/times;


 end
end
% MCRB
    MCRB=zeros（1length（snr）） ;
    
    for k=1:length（snr）
        EsN0（k）=10^（snr（k）/10）;%snr和ESN之间的转化。。
        MCRB（k）=3/（2*pi^2*P^3*EsN0（k））  ; %T^2*
    end
    
figure （1）

% semilogy（snrMCRB‘- k‘）;
% hold on
semilogy（snrvar1‘-^ k‘）;
hold on

xlabel（‘       Es/N0 dB‘）;
ylabel（‘  频差归一化均方差‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       2123  2012-11-07 18:22  Fitz\Fitz.m

     文件        523  2012-07-12 08:37  Fitz\qpsk.m

     目录          0  2014-06-14 11:42  Fitz

----------- ---------  ---------- -----  ----

                 2646                    3