不能调用函数，实现GMSK

好的，如果您不想使用 Communications Toolbox 中的 GMSK 调制和解调函数，您可以手动实现 GMSK 调制和解调。下面给出一个基本的 MATLAB 代码实现：

% 设置调制参数
Fs = 100e3;            % 采样率
fc = 5e3;              % 载波频率
BT = 0.3;              % 带宽时间积
data = randi([0 1], 1000, 1); % 待调制数据

% 常数设置
Ts = 1/Fs;             % 采样时间间隔
Tb = 1/BT/Fs;          % Bit时间
N = length(data);      % 数据长度
Nsym = N;              % 符号数
phi = 0;               % 初始相位
beta = 0.3;            % GMSK指数

% GMSK 调制
s = zeros(Nsym,1);     % 初始化调制信号
for i=1:N
    % 计算相位变化
    delta_phi = 2*pi*fc*Tb*beta*(2*data(i)-1);
    phi = phi + delta_phi;
    % 计算调制信号
    s((i-1)*Tb/Ts+1:i*Tb/Ts) = sqrt(2/Tb)*cos(phi+(pi/2)*i);
end

% 添加高斯白噪声
SNR = 10; % dB
s_noisy = awgn(s, SNR, 'measured');

% GMSK 解调
data_demod = zeros(N,1); % 初始化解调数据
phi_demod = 0;           % 初始相位
for i=1:N
    % 计算相位变化
    delta_phi_demod = atan2(imag(s_noisy(i)),real(s_noisy(i))) - phi_demod;
    % 限制相位变化在-pi/2到pi/2之间
    if delta_phi_demod > pi/2
        delta_phi_demod = delta_phi_demod - pi;
    elseif delta_phi_demod < -pi/2
        delta_phi_demod = delta_phi_demod + pi;
    end
    phi_demod = phi_demod + delta_phi_demod;
    % 解调数据
    data_demod(i) = (phi_demod > 0);
end

% 验证解调数据是否正确
isequal(data, data_demod)

这个代码实现了 GMSK 调制和解调，您可以根据需要修改参数并进行测试。