用MATLAB写一个GMSK解调后的误码率曲线图

以下是一个简单的MATLAB示例代码，用于生成GMSK解调后的误码率曲线图。

% 设置GMSK参数
Fs = 1e6; % 采样率
fc = 100e3; % 载波频率
BT = 0.3; % 带宽-符号时间乘积
fcutoff = BT/2 * Fs; % GMSK滤波器截止频率
span = 10; % 滤波器时间跨度（以符号时间为单位）

% 生成GMSK信号
data = randi([0,1],1,10000); % 随机生成数据
mod_sig = gmskmod(data,Fs,fc,BT); % 进行GMSK调制

% 接收端解调
demod_sig = gmskdemod(mod_sig,Fs,fc,BT); % 进行GMSK解调

% 计算误码率
EbNo = 0:2:20; % 信噪比范围
ber = zeros(size(EbNo)); % 初始化误码率
for n = 1:length(EbNo)
    snr = EbNo(n) + 10*log10(log2(2*BT)); % 计算信噪比
    rx_sig = awgn(demod_sig,snr,'measured'); % 加入高斯白噪声
    rx_data = gmskdemod(rx_sig,Fs,fc,BT)>0.5; % 解调接收信号
    [num_err,ber(n)] = biterr(data,rx_data); % 计算误码率
end

% 绘制误码率曲线图
semilogy(EbNo,ber,'-o');
grid on;
xlabel('Eb/No (dB)');
ylabel('Bit Error Rate');
title('GMSK Demodulation BER Curve');

该代码生成了长度为10000的随机二进制数据，将其GMSK调制，并进行GMSK解调。然后在不同的信噪比下，通过添加高斯白噪声来模拟信道，计算误码率，并绘制误码率曲线图。