%循環(huán)累量，載波估計
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循環(huán)均值
%?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二階循環(huán)均值‘）;
%?xlabel（‘循環(huán)頻率‘）;
%?ylabel（‘循環(huán)均值幅度‘）;


%2四階循環(huán)均值
%?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
%?%四階循環(huán)累量
%?c40=zeros（12*df）;
%?c40=abs（Ma4-3*Ma2.^2）;
%?
%?figure（2）;
%?