生成GMSK仿真的MATLAB代码

这是一个简单的MATLAB代码，用于生成GMSK信号的仿真：

% GMSK信号仿真
fc = 2e6;           % 载频频率
fs = 10e6;          % 采样频率
BT = 0.5;           % 带限时间
Tb = 1/fs;          % 抽样时间
Ts = 1/(2*BT);      % 时钟间隔
fdev = BT/Ts;       % 频偏
fcutoff = 5*BT/Tb;  % 截止频率
N = 1000;           % 仿真数据点数

% 生成载频信号
t = (0:N-1)*Tb;
carrier = cos(2*pi*fc*t);

% 生成调制信号
data = randi([0 1], 1, N);
modulated = gmskmod(data, fc, fdev, fs, 'InputType', 'bit');

% 生成带噪声信号
SNR = 10;           % 信噪比
noisy = awgn(modulated, SNR, 'measured');

% 显示信号图像
subplot(3,1,1);
plot(t, data);
title('原始数据');

subplot(3,1,2);
plot(t, modulated);
title('GMSK调制信号');

subplot(3,1,3);
plot(t, noisy);
title(['带噪声的GMSK信号 (SNR = ' num2str(SNR) 'dB)']);

这段代码中，我们首先设定了载频频率 `fc` 和采样频率 `fs`，以及带限时间 `BT`，时钟间隔 `Ts`，频偏 `fdev` 和截止频率 `fcutoff`。接着，我们生成一个长度为 `N` 的随机二进制数据序列，并使用 `gmskmod` 函数将其调制成GMSK信号。然后，我们添加噪声并显示信号图像。