计算bpsk的循环累积量

%循环累量，载波估计
clear;clc;
len=1000;
st=rand（1len）;
% st=ones（1len）;
st=round（st）;
st=2*st-1;
gt=repmat（st251）;
gt2=reshape（gt125*len）;
%hamming窗
winh=hamming（25）;
winh=repmat（winhlen1）;
winh=winh‘;
%加窗
gt2=gt2.*winh;

n=0:25*len-1;
fc=1000;
fs=25000/8;
Ts=25/fs;%码元宽度
alpha=pi/8;
e=cos（2*pi*fc*n/fs+alpha）;
cst=gt2.*e;
cst2=awgn（cst10‘measured‘）;
% cst2=cst;


%hilbert变换
r=hilbert（cst2）;

% %1循环均值
% a=fc;
% df=150;
% Ma2=zeros（12*df）;
% for k1=a-df:a+df-1
%     ae=exp（-1i*4*pi*k1*n/fs）;
%     M1=cst2.^2.*ae;
%     Ma2（k1-a+df+1）=mean（M1）;
% end
% figure（1）;
% stem（a-df:a+df-1abs（Ma2））;grid on;
% title（‘bpsk二阶循环均值‘）;
% xlabel（‘循环频率‘）;
% ylabel（‘循环均值幅度‘）;


%2四阶循环均值
% Ma4=zeros（12*df）;
% for k1=a-df:a+df-1
%     ae=exp（-1i*8*pi*k1*n/fs）;
%     M1=cst2.^4.*ae;
%     Ma4（k1-a+df+1）=mean（M1）;
% end
% %四阶循环累量
% c40=zeros（12*df）;
% c40=abs（Ma4-3*Ma2.^2）;
% 
% figure（2）;
%