上面的代码不是余弦调制滤波器组的，错了

非常抱歉，我再次犯了错误。这是基于余弦调制滤波器组的OFDM系统的Matlab代码：

% OFDM系统参数
N = 64;              % 子载波数
cp_len = 16;         % 循环前缀长度
mod_order = 16;      % 调制阶数（16QAM）
num_bits = 1e6;      % 发送比特数
snr_vec = 0:2:20;    % 信噪比范围
alpha = 0.5;         % 余弦调制滤波器组的滚降系数
L = 8;               % 滤波器个数

% 生成随机比特流
tx_bits = randi([0 1], num_bits, 1);

% 调制
tx_symbols = qammod(tx_bits, mod_order);

% 将符号插入OFDM帧
tx_frames = reshape(tx_symbols, N, []);
tx_ifft = ifft(tx_frames, N, 1);

% 余弦调制滤波器组
filter_bank = zeros(L, N);
for l = 0:L-1
    filter_bank(l+1, :) = sqrt(2/L) * cos(pi/L * (0:N-1) + pi/L * (l+0.5));
end

% 将OFDM帧通过余弦调制滤波器组
tx_filtered = zeros(L*N, size(tx_frames, 2));
for i = 1:size(tx_frames, 2)
    tx_filtered(:,i) = reshape(filter_bank * tx_frames(:,i), [], 1);
end

% 将滤波后的OFDM帧添加循环前缀
tx_cp = [tx_filtered(end-cp_len+1:end,:); tx_filtered];

% OFDM传输信道
rx_cp = awgn(tx_cp, snr_vec(end));

% 去除循环前缀
rx_filtered = rx_cp(cp_len+1:end,:);

% 将滤波后的OFDM帧通过逆余弦调制滤波器组
rx_frames = zeros(N, size(rx_filtered, 2));
for i = 1:size(rx_filtered, 2)
    rx_frames(:,i) = filter_bank' * reshape(rx_filtered(:,i), [], L);
end

% FFT解调
rx_ifft = fft(rx_frames, N, 1);

% 将OFDM符号串转换为比特串
rx_symbols = rx_ifft(:);
rx_bits = qamdemod(rx_symbols, mod_order);

% 计算误比特率
ber_vec = zeros(length(snr_vec), 2);
for i = 1:length(snr_vec)
    % 基于普通OFDM系统的误码率
    rx_bits_ofdm = qamdemod(awgn(qammod(tx_bits, mod_order), snr_vec(i)), mod_order);
    [~, ber_vec(i,1)] = biterr(tx_bits, rx_bits_ofdm);
    
    % 基于余弦调制滤波器组的OFDM系统的误码率
    [~, ber_vec(i,2)] = biterr(tx_bits, rx_bits);
end

% 绘制误码率曲线
semilogy(snr_vec, ber_vec(:,1), 'bo-', snr_vec, ber_vec(:,2), 'rx-');
xlabel('信噪比 (dB)');
ylabel('误比特率');
legend('普通OFDM系统', '余弦调制滤波器组的OFDM系统');

请注意，这只是一个简单的例子，实际上，OFDM系统的误码率性能取决于许多因素，如子载波数、循环前缀长度、调制阶数、滤波器组大小等。因此，您需要根据您的具体应用程序进行调整。