GMSK 代碼

2024年3月30日發(fā)(作者：我有個桃源仙境)

%繪制調(diào)制波形00101010

clear all;

Ts=1/16000; %基帶信號周期為1/16000s,即為16KHz

Tb=1/32000;

32KHz

BbTb=0.5;

Bb=BbTb/Tb;

Fc=32000;

F_sample=64;

B_num=8;

B_sample=F_sample*Fc*Tb;

Dt=1/Fc/F_sample;

t=0:Dt:B_num*Tb-Dt;

T=Dt*length(t);

%輸入信號周期為Ts/2=1/32000s,即

%取BbTb為0.3

%3dB帶寬

%載波頻率為32KHz

%每載波采樣64個點

%基帶信號為8個碼元

%每基帶碼元采樣點數(shù)B_sample=Tb/Dt

%采樣間隔

%仿真時間

%仿真時間值

Ak=[0 0 1 0 1 0 1 0]; %產(chǎn)生8個基帶信號

Ak=2*Ak-1;

gt=ones(1,B_sample); %每碼元對應(yīng)的載波信號

Akk=sigexpand(Ak,B_sample); %碼元擴展

temp=conv(Akk,gt); %碼元擴展

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))/sqrt(2))/2-erfc(2*pi*Bb*(tt+Tb/2)/s

qrt(log(2))/sqrt(2))/2;

J_g=zeros(1,length(gausst)); %使J_g 的長度和Gausst的一樣

for i=1:length(gausst)

if i==1