GMSK的matlab源程序

%GMSK
%绘制调制波形00101010
clear all;
Ts=1/1000000;                         %基带信号周期为1/1000000s即为1MHz
Tb=1/2000000;                         %输入信号周期为Ts/2=1/2000000s即2MHz
BbTb=0.5;                           %取BbTb为0.3
Bb=BbTb/Tb;                         %3dB带宽
Fc=2000000;                           %载波频率为2MHz
F_sample=64;                        %每载波采样64个点
B_num=8;                            %基带信号为8个码元
B_sample=F_sample*Fc*Tb;            %每基带码元采样点数B_sample=Tb/Dt
Dt=1/Fc/F_sample;                   %采样间隔
t=0:Dt:B_num*Tb-Dt;                 %仿真时间
T=Dt*length（t）;                     %仿真时间值
Ak=[0 0 1 0 1 0 1 0];               %产生8个基带信号
Ak=2*Ak-1;
gt=ones（1B_sample）;         %每码元对应的载波信号
Akk=sigexpand（AkB_sample）;   %码元扩展
temp=conv（Akkgt）;                 %码元扩展
Akk=temp（1:length（Akk））;            %码元扩展

tt=-2.5*Tb:Dt:2.5*Tb-Dt;   
%g（t）=Q[2*pi*Bb*（t-Tb/2）/sqrt（log（2））]-Q[2*pi*Bb*（t+Tb/2）/sqrt（log（2））];
%Q（t）=erfc（t/sqrt（2））/2;
gausst=erfc（2*pi*Bb*（tt-Tb/2）/sqrt（log（2））/sqr