编写matlab代码，实现%基带信号周期为1/16000s,即为16KHz%输入信号周期为Ts/2=1/32000s，即32KHz%取BbTb为0.5%3dB带宽%载波频率为32KHz%每载波采样64个点%基带信号为8个码元%每基带码元采样点数B_sample=Tb/Dt%采样间隔%仿真时间%仿真时间值%产生8个基带信号 。使用GMSK调制完成卫星通信中的调制方法(绘制调制波形00101010)

以下是实现该功能的 MATLAB 代码：

%% 参数设置
fc = 32000;                  % 载波频率
fs = 64*fc;                  % 采样频率
T = 8/fc;                    % 基带信号周期
Ts = 1/(2*fc);               % 输入信号周期
BbTb = 0.5;                  % 3dB带宽
B = BbTb/T;                  % 带宽
B_sample = round(T/(1/fs));  % 每基带码元采样点数
Dt = 1/fs;                   % 采样间隔
t = 0:Dt:T-Dt;               % 基带信号采样时间
tb = 0:Dt:T/B-Dt;            % 基带码元采样时间
N = length(tb);              % 每个基带码元采样点数
data = [0 0 1 0 1 0 1 0];    % 基带信号

%% 生成基带信号
bb = zeros(1, length(t));
for i = 1:length(data)
    if data(i) == 0
        bb((i-1)*N+1:i*N) = -1;
    else
        bb((i-1)*N+1:i*N) = 1;
    end
end

%% GMSK调制
kf = B/T;                            % 调频灵敏度
phi = cumsum(2*pi*fc*Dt + 2*pi*kf*bb.*cos(2*pi*fc*t)); % 累积相位
s = cos(phi);                        % GMSK调制信号

%% 绘制调制波形
figure;
plot(t, s);
xlabel('时间(s)');
ylabel('幅度');
title('GMSK调制波形');

运行该代码，即可得到GMSK调制的波形图。