不使用 Communications Toolbox 中的 GMSK 调制和解调函数，实现matlab中GMSK调制与解调

以下是一个简单的matlab代码实现GMSK调制和解调的示例：

GMSK调制：

% 设置参数
fs = 100e3;             % 采样率
T = 1/fs;               % 采样时间间隔
fc = 500e3;             % 载波频率
B = 50e3;               % 带宽
BT = 0.5;               % 带限乘积
fc_dev = B/2/BT;        % 频偏
Ts = 1/fs;              % 符号时间
fc_mod = 10*fc_dev/Ts;  % 调制信号频率

% 产生数字信号
n = 0:999;              % 采样点数
bits = randi([0 1], 1, length(n)); % 产生随机比特流

% 将数字信号转换为基带信号
t = n*T;
s = 2*bits-1;
bb = cumsum(s).*T;

% 高通滤波
[b, a] = butter(6, 2*fc/fs, 'high');
bb_filt = filter(b, a, bb);

% 相位调制
ph = cumsum(2*pi*fc_mod*bb_filt);
x = cos(2*pi*fc*t + ph);

% 低通滤波
[b, a] = butter(6, 2*B/fs, 'low');
x_filt = filter(b, a, x);

% 绘图
subplot(2,1,1);
plot(t, bits, 'b');
axis([0 t(end) -1.5 1.5]);
title('数字信号');
subplot(2,1,2);
plot(t, x_filt, 'r');
axis([0 t(end) -1.5 1.5]);
title('GMSK调制信号');

GMSK解调：

% 设置参数
fs = 100e3;             % 采样率
T = 1/fs;               % 采样时间间隔
fc = 500e3;             % 载波频率
B = 50e3;               % 带宽
BT = 0.5;               % 带限乘积
fc_dev = B/2/BT;        % 频偏
Ts = 1/fs;              % 符号时间
fc_mod = 10*fc_dev/Ts;  % 调制信号频率

% 产生数字信号
n = 0:999;              % 采样点数
bits = randi([0 1], 1, length(n)); % 产生随机比特流

% GMSK调制
t = n*T;
s = 2*bits-1;
bb = cumsum(s).*T;
[b, a] = butter(6, 2*fc/fs, 'high');
bb_filt = filter(b, a, bb);
ph = cumsum(2*pi*fc_mod*bb_filt);
x = cos(2*pi*fc*t + ph);
[b, a] = butter(6, 2*B/fs, 'low');
x_filt = filter(b, a, x);

% GMSK解调
[b, a] = butter(6, 2*B/fs, 'low');
y_filt = filter(b, a, x_filt);
ph_diff = diff(ph);
ph_diff = [ph_diff(1); ph_diff];
y = y_filt.*sin(ph_diff);

% 绘图
subplot(3,1,1);
plot(t, bits, 'b');
axis([0 t(end) -1.5 1.5]);
title('数字信号');
subplot(3,1,2);
plot(t, x_filt, 'r');
axis([0 t(end) -1.5 1.5]);
title('GMSK调制信号');
subplot(3,1,3);
plot(t, y, 'g');
axis([0 t(end) -1.5 1.5]);
title('解调后的数字信号');

需要注意的是，这只是一个简单的实现方法，实际应用中还需要考虑很多因素，比如调制指数、滤波器设计等。建议您在实现之前先进行相关的理论学习和实验验证。