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 = 1e3; % 子载波间隔
Fs = 2.56e6; % 时域离散化采样速率

% 生成QPSK调制所需的符号映射表
qpskMod = comm.QPSKModulator('BitInput',true);
qpskDemod = comm.QPSKDemodulator('BitOutput',true);

% 生成随机数据，并进行QPSK调制
data = randi([0 1], N*log2(M)*L, 1);
modData = qpskMod(data);

% 将调制符号序列映射到不同的子载波上
ofdmData = zeros(N, L);
ofdmData(1:N/2, 1) = modData(1:N/2);
ofdmData(N/2+2:N, 1) = modData(N/2+1:N);

% 进行IFFT变换
ifftData = ifft(ofdmData, N, 1);

% 添加循环前缀
cpLength = N/4;
cpData = [ifftData(N-cpLength+1:N, 1); ifftData(:, 1)];

% 将时域信号串并成一串
txData = reshape(cpData, [], 1);

% 显示调制结果
figure;
subplot(2,1,1);
plot(real(modData), imag(modData), 'o');
title('QPSK Modulated Symbols');
xlabel('Real');
ylabel('Imaginary');

subplot(2,1,2);
plot((0:length(txData)-1)/Fs, real(txData));
title('OFDM Modulated Signal');
xlabel('Time (s)');
ylabel('Amplitude');

方法二：使用FFT/IFFT

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

% 生成QPSK调制所需的符号映射表
qpskMod = comm.QPSKModulator('BitInput',true);
qpskDemod = comm.QPSKDemodulator('BitOutput',true);

% 生成随机数据，并进行QPSK调制
data = randi([0 1], N*log2(M)*L, 1);
modData = qpskMod(data);

% 将调制符号序列映射到不同的子载波上
ofdmData = zeros(N, L);
ofdmData(1:N/2, 1) = modData(1:N/2);
ofdmData(N/2+2:N, 1) = modData(N/2+1:N);

% 进行IFFT变换
ifftData = ifft(ofdmData, N, 1);

% 添加循环前缀
cpLength = N/4;
cpData = [ifftData(N-cpLength+1:N, 1); ifftData(:, 1)];

% 将时域信号串并成一串
txData = reshape(cpData, [], 1);

% 计算OFDM信号的频谱
f = (-N/2:N/2-1)*delta_f;
txData_fft = fft(txData);
txData_fft_shifted = fftshift(txData_fft);

% 显示调制结果
figure;
subplot(2,1,1);
plot(real(modData), imag(modData), 'o');
title('QPSK Modulated Symbols');
xlabel('Real');
ylabel('Imaginary');

subplot(2,1,2);
plot(f, abs(txData_fft_shifted));
title('OFDM Modulated Signal Spectrum');
xlabel('Frequency (Hz)');
ylabel('Magnitude');