帮我写一个基于OFDM的matlab代码，可以绘制误码率随信噪比变化的示意图

好的，以下是一个基于OFDM的MATLAB代码，可以绘制误码率随信噪比变化的示意图：

%% 系统参数设置
N_subcarriers = 64; % 子载波数目
CP_length = 16; % 循环前缀长度
FFT_length = 64; % FFT点数
SNR_dB = 0:2:20; % 信噪比范围
num_bits = 1e6; % 发送比特数目
M = 16; % 调制阶数
bit_per_symbol = log2(M); % 每个符号所含比特数
num_symbols = num_bits / bit_per_symbol; % 发送符号数目

% 生成随机数据
data = randi([0, 1], num_bits, 1);

% 调制
mod_data = qammod(data, M);

%% OFDM调制
% 将调制后的数据进行序列分割
data_matrix = reshape(mod_data, [], N_subcarriers);

% 进行FFT变换
ofdm_signal = ifft(data_matrix, FFT_length, 2);

% 添加循环前缀
ofdm_signal_with_cp = [ofdm_signal(:, end-CP_length+1:end), ofdm_signal];

% 将矩阵转换为向量
tx_signal = ofdm_signal_with_cp(:);

%% 信道传输
for snr_idx = 1:length(SNR_dB)
    % 生成高斯白噪声
    noise = sqrt(0.5 / 10^(SNR_dB(snr_idx) / 10)) * (randn(size(tx_signal)) + 1j * randn(size(tx_signal)));
    
    % 添加噪声
    rx_signal = tx_signal + noise;
    
    % 将向量转换为矩阵
    rx_matrix_with_cp = reshape(rx_signal, [], FFT_length+CP_length);
    
    % 去掉循环前缀
    rx_matrix = rx_matrix_with_cp(:, CP_length+1:end);
    
    % 进行FFT变换
    rx_data_matrix = fft(rx_matrix, FFT_length, 2);
    
    % 将矩阵转换为向量
    rx_data = rx_data_matrix(:);
    
    % 解调
    demod_data = qamdemod(rx_data, M);
    
    % 计算误码率
    num_errors(snr_idx) = sum(demod_data ~= data);
    ber(snr_idx) = num_errors(snr_idx) / num_bits;
end

%% 绘制误码率随信噪比变化的示意图
figure;
semilogy(SNR_dB, ber, 'bo-', 'LineWidth', 2);
grid on;
xlabel('信噪比 (dB)');
ylabel('误码率');
title(['OFDM调制 QAM', num2str(M), '调制']);

运行代码后，会生成一个误码率随信噪比变化的示意图。其中，x轴表示信噪比（dB），y轴表示误码率。可以通过修改代码中的参数来改变信道模型和调制参数，进一步验证OFDM系统的性能。