GMSK 代码

2024年3月30日发(作者：我有个桃源仙境)

%绘制调制波形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个码元

%每基带码元采样点数B_sample=Tb/Dt

%采样间隔

%仿真时间

%仿真时间值

Ak=[0 0 1 0 1 0 1 0]; %产生8个基带信号

Ak=2*Ak-1;

gt=ones(1,B_sample); %每码元对应的载波信号

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