%GMSK
%繪制調(diào)制波形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;????????????%每基帶碼元采樣點數(shù)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];???????????????%產(chǎn)生8個基帶信號
Ak=2*Ak-1;
gt=ones（1B_sample）;?????????%每碼元對應(yīng)的載波信號
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