fsk调制仿真

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‘