clear?all;
clc;

bk0?=?[1?0?1?1?0?0?0?1];
a=bk0*2-1;
omiga1=2*pi*1;
omiga2=2*pi*2;
omigac=（omiga1+omiga2）/2;
omigad=abs（omiga1-omiga2）/2;
fd=omigad/（2*pi）;
Tb=1;
Rb=1/Tb;
h=2*fd/Rb;

phik=[pi?pi/2?pi/3?pi/5?pi/7?pi/9?pi/11?pi/13];
dt?=?0.002;
t=0:dt:（8-dt）;
?Sfsk=cos（omigac*t+（a（1floor（t+1））*omigad）.*t+phik（1floor（t+1）））;
?
subplot（411）
bk=bk0（1floor（t+1））;
plot（tbk‘LineWidth‘