请使用模拟调制与FFT/IFFT两种方法完成OFDM调制的MATLAB程序编写（可参考《详解Matlab/Simulink通信系统建模及仿真》中例9.2（P.264）），

以下是使用模拟调制与FFT/IFFT两种方法完成OFDM调制的MATLAB程序：

% OFDM参数
N = 64; % 子载波数
cp_len = 16; % 循环前缀长度
Fs = 100e6; % 采样率
Ts = 1/Fs; % 采样时间间隔
fc = 2e9; % 载波频率
fd = 10e3; % 多径衰落频率
Ts_fd = 1/fd; % 多径时延
fd_amp = 0.8; % 多径衰落幅度

% 生成随机数据
data = randi([0 1], N, 1);

% 模拟调制方法
t = [0:Ts:N*Ts-1e-9]; % 时间向量
ifft_data = ifft(data); % IFFT
ifft_data_cp = [ifft_data(N-cp_len+1:N); ifft_data]; % 添加循环前缀
ofdm_signal = reshape(ifft_data_cp, [], 1); % 转换为列向量
mult_path = fd_amp*exp(-1j*2*pi*fd*t); % 多径衰落信道
rx_signal_ana = ofdm_signal .* mult_path.' .* exp(1j*2*pi*fc*t); % 加载多径衰落并进行载波调制

% FFT/IFFT方法
data_mat = reshape(data, [], 4); % 将数据分组
ifft_data_mat = ifft(data_mat, N, 2); % 对每组数据进行IFFT
ifft_data_mat_cp = [ifft_data_mat(:,N-cp_len+1:N), ifft_data_mat]; % 添加循环前缀
ofdm_signal_fft = reshape(ifft_data_mat_cp.', [], 1); % 转置并转换为列向量
mult_path_fft = repmat(fd_amp*exp(-1j*2*pi*fd*t), 1, 4); % 将多径衰落信道扩展为与数据组数相同的矩阵
rx_signal_fft = ofdm_signal_fft .* mult_path_fft .* exp(1j*2*pi*fc*t); % 加载多径衰落并进行载波调制

% 绘制时域波形
figure;
subplot(2,1,1);
plot(t, real(rx_signal_ana));
hold on;
plot(t, imag(rx_signal_ana));
title('模拟调制方法');
xlabel('时间(s)');
ylabel('幅度');
legend('实部','虚部');

subplot(2,1,2);
plot(t, real(rx_signal_fft));
hold on;
plot(t, imag(rx_signal_fft));
title('FFT/IFFT方法');
xlabel('时间(s)');
ylabel('幅度');
legend('实部','虚部');

该程序使用两种方法生成OFDM信号，并进行多径衰落和载波调制。其中，模拟调制方法直接对整个OFDM符号进行IFFT和循环前缀添加，而FFT/IFFT方法是将数据分组后对每组数据进行IFFT，再添加循环前缀，最后将所有OFDM符号组合成一个列向量。可以通过绘制时域波形来比较这两种方法生成的OFDM信号。