OFDM初始调制参数： 星座调制：QPSK 子载波数量：16 OFDM符号个数：1 子载波间隔：1kHz 时域离散化采样速率：2.56MHz （2）请使用模拟调制与FFT/IFFT两种方法完成OFDM调制的MATLAB程序编写

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

方法一：使用模拟调制

%% OFDM初始调制参数
M = 4; % 星座调制：QPSK
N = 16; % 子载波数量
L = 1; % OFDM符号个数
delta_f = 1000; % 子载波间隔
fs = 2.56e6; % 时域离散化采样速率

%% 生成随机QPSK数据
data = randi([0 M-1], N*L, 1);
dataMod = qammod(data, M);

%% 将数据映射到子载波上
dataModMat = reshape(dataMod, N, L);
ofdmSig = ifft(dataModMat, N);

%% 加上循环前缀
cpLen = N/4;
ofdmSig_cp = [ofdmSig((N-cpLen+1):N,:); ofdmSig];

%% 转换为时域信号
t = (0:(N+cpLen-1))/fs;
ofdmSig_t = zeros(length(t), L);
for i = 1:L
    ofdmSig_t(:,i) = real(ofdmSig_cp(:,i)*exp(1j*2*pi*delta_f*(i-1)*t));
end

%% 画图展示OFDM信号
figure;
plot(t, ofdmSig_t);
xlabel('Time (s)');
ylabel('Amplitude');
title('OFDM Signal (Using Analog Modulation)');

方法二：使用FFT/IFFT

%% OFDM初始调制参数
M = 4; % 星座调制：QPSK
N = 16; % 子载波数量
L = 1; % OFDM符号个数
delta_f = 1000; % 子载波间隔
fs = 2.56e6; % 时域离散化采样速率

%% 生成随机QPSK数据
data = randi([0 M-1], N*L, 1);
dataMod = qammod(data, M);

%% 将数据映射到子载波上
dataModMat = reshape(dataMod, N, L);
ofdmSig_fft = ifft(dataModMat, N);

%% 加上循环前缀
cpLen = N/4;
ofdmSig_cp = [ofdmSig_fft((N-cpLen+1):N,:); ofdmSig_fft];

%% 进行FFT变换
ofdmSig_freq = zeros(N+cpLen, L);
for i = 1:L
    ofdmSig_freq(:,i) = fft(ofdmSig_cp(:,i));
end

%% 画图展示OFDM信号
f = (-N/2:(N/2-1))*delta_f;
figure;
plot(f, abs(ofdmSig_freq));
xlabel('Frequency (Hz)');
ylabel('Amplitude');
title('OFDM Signal (Using FFT)');

注意：这两种方法生成的OFDM信号是等效的，只是实现方式不同。