帮我写一个基于OFDM的matlab程序，可以绘制原始数据和接收数据的波形图和对比图

OFDM是一种调制技术，用于将高速数据流分成多个较低速率的子流，以便在信道中进行传输。下面是一个基于OFDM的MATLAB程序，可以绘制原始数据和接收数据的波形图和对比图：

% OFDM基本参数
N = 64; % 帧长度
cpLen = 16; % 循环前缀长度
Fs = 100e6; % 采样率
Ts = 1/Fs; % 采样时间
M = 16; % 16QAM调制
k = log2(M); % 每个符号的比特数
numSym = 10; % 发送符号的数量

% 生成随机数据
txData = randi([0 1], k*N*numSym, 1);

% 16QAM调制
txData = qammod(txData, M);

% 将数据转换为矩阵
txDataMat = reshape(txData, k*N, numSym);

% IFFT
txOfdm = ifft(txDataMat);

% 添加循环前缀
txCp = [txOfdm(N-cpLen+1:N,:); txOfdm];

% 串并转换
txSerial = reshape(txCp, [], 1);

% 传输信道
rxSerial = awgn(txSerial, 20, 'measured'); % 加入高斯白噪声

% 并串转换
rxCp = reshape(rxSerial, N+cpLen, []);

% 去除循环前缀
rxOfdm = rxCp(cpLen+1:end,:);

% FFT
rxDataMat = fft(rxOfdm);

% 将数据转换为向量
rxData = reshape(rxDataMat, k*N*numSym, 1);

% 解调
rxDataDemod = qamdemod(rxData, M);

% 绘图
figure;
subplot(3,1,1);
plot(real(txData));
hold on;
plot(imag(txData));
title('发送数据');
legend('实部','虚部');
xlabel('符号');
ylabel('幅度');
grid on;

subplot(3,1,2);
plot(real(rxData));
hold on;
plot(imag(rxData));
title('接收数据');
legend('实部','虚部');
xlabel('符号');
ylabel('幅度');
grid on;

subplot(3,1,3);
plot(real(txData));
hold on;
plot(real(rxDataDemod));
title('发送数据与接收数据比较');
legend('发送数据','接收数据');
xlabel('符号');
ylabel('幅度');
grid on;

这个程序可以生成随机数据并进行16QAM调制，然后将数据转换为OFDM符号并添加循环前缀。发送数据通过AWGN信道传输，并在接收端进行解调和去除循环前缀。最后，程序绘制原始数据和接收数据的波形图和对比图。