给我一个matlab代码，生成OFDM信号

以下是一个简单的 MATLAB 代码来生成 OFDM 信号：

%% OFDM 信道参数
N = 64; % 子载波数量
cp_len = 16; % 循环前缀长度
fs = 16000; % 采样率
Fd = 2000; % 多径衰落的最大多普勒频率

%% 生成原始数据
data_len = N * 100; % 数据长度
data = randi([0 1], data_len, 1); % 随机生成二进制数据

%% 将原始数据映射到符号上
M = 4; % 星座点数
symbols = qammod(data, M); % 使用 QAM 调制将二进制数据映射到符号上

%% 将符号插入子载波上
tx_data = zeros(N, length(symbols) / N);
for i = 1:size(tx_data, 2)
    tx_data(:,i) = symbols((i-1)*N+1:i*N);
end

%% 添加循环前缀
tx_cp_data = [tx_data(end-cp_len+1:end,:); tx_data];

%% IFFT 变换
tx_ifft_data = ifft(tx_cp_data, N, 1);

%% 生成高斯白噪声信道
t = (0:data_len/fs:((data_len-1)/fs))';
h = (randn(length(t), 1) + 1i * randn(length(t), 1)) .* exp(-1j*2*pi*Fd*t);

%% 信道传输
rx_ifft_data = tx_ifft_data .* h;

%% 接收端去除循环前缀
rx_cp_data = rx_ifft_data(cp_len+1:end,:);

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

%% 将符号解调回原始数据
rx_symbols = reshape(rx_data, [length(symbols), 1]);
rx_data = qamdemod(rx_symbols, M);

%% 误码率计算
ber = sum(data ~= rx_data) / data_len;
fprintf('误码率为: %f\n', ber);

这个代码生成一个长度为 $64$ 的 OFDM 符号，采用 QAM 调制将二进制数据映射到符号上，插入循环前缀，进行 IFFT 变换，通过高斯白噪声信道传输，然后进行 FFT 变换，解调回原始数据，并计算误码率。