帮我写一个基于OFDM的matlab代码，包括绘制随机数据前后的波形图和误码率随信噪比变化图

好的，以下是一个基于OFDM的Matlab代码，包括随机数据的波形图和误码率随信噪比变化图：

% OFDM系统仿真
clear all;
close all;

%% 参数设置
N = 64; % 子载波数
cp_len = 16; % 循环前缀长度
SNR_dB = 0:2:20; % 信噪比范围
M = 16; % 调制阶数
num_bit = N*log2(M); % 每个OFDM符号传输的比特数
num_sym = 1000; % 发送OFDM符号数
BER = zeros(1,length(SNR_dB)); % 误码率

%% OFDM信号产生
for snr_idx = 1:length(SNR_dB)
    % 产生随机数据
    tx_bit = randi([0 1], 1, num_sym*num_bit);
    tx_bit_matrix = reshape(tx_bit, num_bit, num_sym).'; % 分组
    
    % 星座点映射
    tx_sym_index = bi2de(reshape(tx_bit_matrix, num_sym, num_bit/num_sym), 'left-msb');
    tx_sym = qammod(tx_sym_index, M);
    
    % IFFT
    tx_ifft = ifft(tx_sym, N, 2);
    
    % 循环前缀
    tx_cp = [tx_ifft(:,N-cp_len+1:N), tx_ifft];
    
    % 并行发送
    tx_signal = tx_cp.';
    
    %% 信道传输
    % 添加高斯白噪声
    SNR = 10^(SNR_dB(snr_idx)/10);
    noise_var = 1/SNR;
    noise = sqrt(noise_var/2)*(randn(size(tx_signal)) + 1j*randn(size(tx_signal)));
    rx_signal = tx_signal + noise;
    
    %% 接收端处理
    % 并行接收
    rx_cp = rx_signal.';
    
    % 去掉循环前缀
    rx_ifft = rx_cp(:,cp_len+1:cp_len+N);
    
    % FFT
    rx_sym = fft(rx_ifft, N, 2);
    
    % 星座点解映射
    rx_sym_index = qamdemod(rx_sym, M);
    rx_bit_matrix = de2bi(rx_sym_index, num_bit/num_sym, 'left-msb');
    rx_bit = reshape(rx_bit_matrix.', [], 1).'; % 转为行向量
    
    %% BER计算
    bit_err = sum(xor(tx_bit,rx_bit));
    BER(snr_idx) = bit_err/length(tx_bit);
end

%% 绘图
% 随机数据波形图
subplot(2,1,1);
plot(real(tx_signal(:)), '-r'); hold on;
plot(real(rx_signal(:)), '-b'); hold off;
xlabel('采样点'); ylabel('幅度');
title('OFDM随机数据波形图');
legend('发送信号', '接收信号');

% 误码率随信噪比变化图
subplot(2,1,2);
semilogy(SNR_dB, BER, '-o');
xlabel('信噪比(dB)'); ylabel('误码率');
title('误码率随信噪比变化图');
grid on;

希望对你有所帮助！