【matlab】 GMSK的调制与解调【附详尽注释】

• 简介
• code

1 简介

MSK调制是调制指数为0.5的二元数字频率调制，具有很好的特性，如恒包络、相对窄的带宽、并可以相干检测。MSK【最小频移键控】信号在任一码元间隔内，其相位变化为Π/2，而在码元转换时刻保持相位连续。

然而，MSK信号的相位变化是折线，在码元转换时刻会产生尖角，从而使其频谱特性的旁瓣降缓慢，带外辐射相对较大。移动数字通信中采用高速传输速率时，要求邻道带外辐射低于－(60~80)dB，而MSK信号不能满足功率谱在相邻信道的取值低于主瓣峰值60dB以上的要求，所以需寻求进一步压缩带宽的方法。

为了进一步改善MSK【最小频移键控】的频谱特性，有效的办法是对基带信号进行平滑处理，使调制后的信号相位在码元转换时刻不仅连续而且变化平滑，从而达到改善频谱特性的目的。

GMSK作为MSK的改进型，即是以高斯低通滤波器作为预调滤波基带滤波器的MSK方式，所以称为高斯MSK或GMSK。

2 code【附详尽注释】

% 附录：MATLAB程序
%绘制调制波形00101010
%
clear all;
close all;

%******************** Preparation part ************************************
% Ts=1/16000;            %基带信号周期为1/16000s,即为16KHz

Tb=1/32000;            %[码元]输入信号周期为Ts/2=1/32000s,即32KHz【奈奎斯特采样】
BbTb=0.5;              %取BbTb为0.5,调制指数为0.5的二元数字频率调制
Bb=BbTb/Tb;            %3dB带宽              -_-半带宽=码元频率一半-_-!
Fc=32000;              %载波频率为32KHz      -_-载波频率=码元频率-_-!
Fc_sample=64;          %每载波采样64个点
B_num=8;               %基带信号为8个码元
Dt=1/(Fc*Fc_sample);   %采样间隔[载波周期/采样点数][4.88281250000000e-07]
B_sample=Tb/Dt;        %每基带码元采样点数 B_sample=Tb/Dt[输入信号/采样间隔]
t=0:Dt:B_num*Tb-Dt;    %仿真时间离散点[采样间隔，码元数*（时间/码元）]
T=Dt*length(t);        %仿真时间值[采样间隔*512个采样点]
Ak=[0 0 1 0 1 0 1 0];        %产生8个基带信号[8个比特]
Ak=2*Ak-1;    %[多此一举]    %单极性码元—>双极性码元
gt=ones(1,B_sample);         %每码元对应的载波信号[1*64]
Akk=sigexpand(Ak,B_sample);  %码元扩展[64*8->1*512]
temp=conv(Akk,gt);           %码元扩展[卷积向量Akk和gt 512+64-1]
Akk=temp(1:length(Akk));     %码元扩展[取出temp变量中前512个，类似于银行存款，第一年存入的钱会一直享受利息值最后一年，第二年存入的钱会一直享受利息值最后一年]

%************************* Filter initialization **************************

tt=-2.5*Tb:Dt:2.5*Tb-Dt;   %{2.5*码元周期}/{采样间隔&