请告诉我GMSK的MATLAB调制代码

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

% 设置参数
fc = 2e3;  % 载波频率
fs = 10e3; % 采样率
Tb = 1/fs; % 采样时间间隔
T = 1;     % 信号时间长度
Rb = 1e3;  % 数据速率
Ns = T/Tb; % 采样点数
Nb = Rb*T; % 数据比特数
BT = 0.3;  % 带宽时间乘积
f_dev = BT*Rb/2; % 频偏

% 生成随机数据
data = randi([0 1],Nb,1);

% 数据调制
mod_data = gmsk_mod(data, f_dev, Tb, fs);

% 生成时域波形
t = linspace(0, T, Ns);
y = real(mod_data).*cos(2*pi*fc*t) - imag(mod_data).*sin(2*pi*fc*t);

% 画图
figure;
subplot(2,1,1);
plot(t, data, 'LineWidth', 2);
xlabel('Time (s)');
ylabel('Data');
title('Input Data');
grid on;

subplot(2,1,2);
plot(t, y, 'LineWidth', 2);
xlabel('Time (s)');
ylabel('Amplitude');
title('GMSK Modulated Signal');
grid on;

% GMSK调制函数
function mod_data = gmsk_mod(data, f_dev, Tb, fs)
    % 预处理
    Ns = length(data);
    t = linspace(0, Tb*Ns, Ns);
    phi = 0;
    x = zeros(1, Ns);

    % GMSK调制
    for i = 1:Ns
        if i == 1
            x(i) = data(i) - phi;
        else
            x(i) = data(i) - data(i-1);
        end
        mod_data(i) = exp(1j*2*pi*f_dev*t(i))*exp(1j*phi);
        phi = phi + 2*pi*f_dev*Tb*x(i);
    end
end

其中，`data`为待调制的数据，`f_dev`为频偏，`Tb`为采样时间间隔，`fs`为采样率，`mod_data`为调制后的数据。这段代码实现了GMSK调制，并且画出了原始数据和调制信号的时域波形图。