%GMSK
%繪制調(diào)制波形00101010
clear?all;
Ts=1/1000000;?????????????????????????%基帶信號(hào)周期為1/1000000s即為1MHz
Tb=1/2000000;?????????????????????????%輸入信號(hào)周期為Ts/2=1/2000000s即2MHz
BbTb=0.5;???????????????????????????%取BbTb為0.3
Bb=BbTb/Tb;?????????????????????????%3dB帶寬
Fc=2000000;???????????????????????????%載波頻率為2MHz
F_sample=64;????????????????????????%每載波采樣64個(gè)點(diǎn)
B_num=8;????????????????????????????%基帶信號(hào)為8個(gè)碼元
B_sample=F_sample*Fc*Tb;????????????%每基帶碼元采樣點(diǎn)數(shù)B_sample=Tb/Dt
Dt=1/Fc/F_sample;???????????????????%采樣間隔
t=0:Dt:B_num*Tb-Dt;?????????????????%仿真時(shí)間
T=Dt*length（t）;?????????????????????%仿真時(shí)間值
Ak=[0?0?1?0?1?0?1?0];???????????????%產(chǎn)生8個(gè)基帶信號(hào)
Ak=2*Ak-1;
gt=ones（1B_sample）;?????????%每碼元對(duì)應(yīng)的載波信號(hào)
Akk=sigexpand（AkB_sample）;???%碼元擴(kuò)展
temp=conv（Akkgt）;?????????????????%碼元擴(kuò)展
Akk=temp（1:length（Akk））;????????????%碼元擴(kuò)展

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