请你写一个基于OFDM的MATLAB代码，可以绘制QPSK信号接受前后的波形图

以下是一个基于OFDM的MATLAB代码，可以绘制QPSK信号接受前后的波形图：

%% OFDM信号的产生与接收
clear all;
close all;

%% 参数设置
Nq = 2;             % 每个子载波的比特数
N = 64;             % 子载波数
cp_len = 16;        % 循环前缀长度
snr = 20;           % 信噪比

%% 信号产生
% 随机产生N个QPSK符号
data = randi([0 3],1,N*Nq);
% QPSK调制
mod_data = exp(1j*pi/4*(2*data+1));
% 矩阵重塑
mod_data = reshape(mod_data,N,Nq);

%% IFFT变换
ifft_data = ifft(mod_data,N,1);

%% 循环前缀添加
cp_data = [ifft_data(N-cp_len+1:N,:);ifft_data];

%% 并串转换
tx_data = cp_data(:);

%% 信道模型
% AWGN信道
rx_data = awgn(tx_data,snr,'measured');
% 串并转换
rx_data = reshape(rx_data,cp_len+N,Nq);

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

%% FFT变换
fft_data = fft(rx_data,N,1);

%% QPSK解调
rx_mod_data = reshape(fft_data,N*Nq,1);
rx_data = (rx_mod_data<0)*3+(real(rx_mod_data)>0)*1;

%% 画图
figure;
subplot(2,1,1);
plot(real(tx_data),'r');hold on;
plot(real(rx_data),'b');
title('接收前后的波形图(实部)');
legend('Tx','Rx');
subplot(2,1,2);
plot(imag(tx_data),'r');hold on;
plot(imag(rx_data),'b');
title('接收前后的波形图(虚部)');
legend('Tx','Rx');