clc;clear?all;close?all;
fs?=?50e6;?%采樣頻率
ts?=?1/fs;?
num?=?2e6;??%數據長度
SNR?=?-15;
real_fc?=?10000500;?%實信號頻率
data?=?sin（2*pi*real_fc*（0:num-1）*ts+pi/4）+sqrt（10^（SNR/10））*randn（1num）;?%科斯塔斯環的輸入信號
fc?=?10000000;?%本地頻率
K=1e4;
n?=?fs/K;?%累積時間為0.1ms
nn?=?[0:n-1];
nf?=?floor（length（data）/n）;%?將輸入數據分成0.1ms的多個數據塊
wfc?=?2*pi*fc;??%本地信號
phi_prv?=?0;
temp?=?0;
frame?=?0;
carrier_phase?=?0;
phase?=?0;

%環路濾波器的參數
K=1e4;
k1=32e3;
bl=0.0001*K;
w=bl/0.53;
T=1/K;
k0=2*pi*T*fs/2^32;
%?k0=1;
%?c1=2*0.707*w*T/k0;
%?c2=（w*T）^2/k0;
c1=1203.9;
c2=16.22;?
%?c1=18237;
%?c2=2433.48;??
%??c1=10;
%??c2=5;
%數據率為4kbps1bit數據采樣點數為12.5K，每5K個點累加一次，相當于4個數據累加清零一次，
%w=0.01*4KT=1/4Kk0=2*pi*T*fs/2^32c1=2*0.707*w*T/k