gmsk的matlab实现

以下是GMSK调制的MATLAB示例代码：

% 设置参数
Rb = 1000; %比特率
Tb = 1/Rb; %比特长度
Fc = 2*Rb; %载波频率
Fs = 10*Fc; %采样率
Ts = 1/Fs; %采样时间
M = 2; %调制阶数
BT = 0.3; %带宽时间积(BT积)
L = 10; %滤波器长度
alpha = sqrt(2/log2(M))*BT/(Tb*sqrt(pi)); %GMSK的alpha参数
beta = alpha/Tb; %GMSK的beta参数

% 生成随机数字0和1
data = randi([0 1], 1, 1000);

% 比特串转换为符号串
symbols = 2*data-1;

% GMSK调制
t = (0:length(symbols)-1)*Tb;
phase = cumsum((1-alpha)*symbols.*cos(2*pi*Fc*t) + beta*symbols.*sin(2*pi*Fc*t));
s = cos(phase);

% 低通滤波
h = firrcos(L, Fc/2, 0.5/BT, Fs, 'rolloff', 'sqrt');
s_filtered = filter(h, 1, s);

% 绘制输出信号图形
subplot(2,1,1);
plot(t, symbols);
title('Input Symbols');
ylabel('Amplitude');
ylim([-1.5 1.5]);

subplot(2,1,2);
plot(t, s_filtered);
title('GMSK Modulated Signal');
xlabel('Time (s)');
ylabel('Amplitude');

上述代码将生成一个长度为1000的随机比特流，并将其转换为符号串。然后它对符号串进行GMSK调制，并使用带宽时间积为0.3的根号升余弦滤波器进行低通滤波。最后，绘制输入符号和输出信号的图形。