请你写一个基于OFDM的MATLAB代码，可以绘制原始数据和接受数据的波形以及数据误码率随信噪比的变化图

好的，以下是一个基于OFDM的MATLAB代码，可以绘制原始数据和接收数据的波形以及数据误码率随信噪比的变化图：

%% OFDM基带信号发射与接收
clc; clear all; close all;
%% 系统参数设置
Nfft = 64; % FFT点数
Ncp = 16; % 循环前缀长度
M = 4; % 调制阶数
L = 8; % 并行传输子载波数
snr_db = 0 : 2 : 20; % 信噪比范围
N = 100000; % 仿真发送数据点数
%% 数据生成
data = randi([0, M-1], N, L); % 生成L个QAM调制的随机数据
mod_data = qammod(data, M); % QAM调制
%% IFFT、加CP、串并转换
tx_signal = zeros(N, Nfft); % 发送信号矩阵
for i = 1 : L
    tx_signal(:, i) = ifft(mod_data(:, i), Nfft); % IFFT
    tx_signal(:, i) = [tx_signal(end-Ncp+1:end, i); tx_signal(:, i)]; % 加循环前缀
end
tx_signal = tx_signal.'; % 串并转换
%% 仿真信道传输
rx_signal = zeros(size(tx_signal)); % 接收信号矩阵
for i = 1 : length(snr_db)
    for j = 1 : L
        snr = 10^(snr_db(i)/10); % 信噪比
        noise_var = 1 / snr; % 噪声方差
        noise = sqrt(noise_var/2) * (randn(size(tx_signal(j,:))) + 1j * randn(size(tx_signal(j,:)))); % 产生高斯白噪声
        rx_signal(j,:) = sqrt(snr) * tx_signal(j,:) + noise; % 加噪声
    end
    %% 并串转换、去CP、FFT
    rx_signal = rx_signal.'; % 并串转换
    rx_signal = rx_signal(Ncp+1:end, :); % 去掉循环前缀
    rx_data = zeros(N, L); % 接收数据矩阵
    for j = 1 : L
        rx_data(:, j) = fft(rx_signal(:, j), Nfft); % FFT
    end
    rx_data = qamdemod(rx_data, M); % QAM解调
    %% 统计误码率
    [~, ber(i)] = biterr(data, rx_data); % 统计误码率
end
%% 绘制结果
figure(1);
subplot(211); plot(real(tx_signal(:,1))); hold on; plot(imag(tx_signal(:,1))); hold off; title('发送信号波形');
subplot(212); plot(real(rx_signal(:,1))); hold on; plot(imag(rx_signal(:,1))); hold off; title('接收信号波形');
figure(2); semilogy(snr_db, ber); xlabel('信噪比(dB)'); ylabel('误码率');

我们可以通过修改系统参数来进行不同情况下的仿真，其中：

- Nfft：FFT点数，也即子载波数；
- Ncp：循环前缀长度；
- M：调制阶数，这里为4QAM；
- L：并行传输子载波数，也即每个OFDM符号中携带的数据量；
- snr_db：信噪比范围，以dB为单位，可以自己设定；
- N：仿真发送数据点数。

运行代码后，会得到两张图：一张是发送信号和接收信号的波形，另一张是误码率随信噪比的变化图。需要注意的是，由于随机数的生成具有一定的随机性，因此误码率曲线可能存在一定的波动。