正交频分多址接入（Orthogonal Frequency Division Multiple Access, OFDMA）是一种广泛应用于4G LTE和5G无线通信系统中的多址接入技术。它通过将频域划分为多个正交子载波，并为不同用户分配不同的子载波组，实现了高效的频谱利用和多用户接入。本示例将展示如何使用MATLAB/Simulink搭建一个简化的OFDMA系统模型，并通过AWGN信道验证其性能。

一、背景介绍

OFDMA概述：
- OFDMA是OFDM（正交频分复用）的多址扩展，允许在同一时频资源上同时服务多个用户。
- 每个用户被分配一组正交子载波，避免了用户之间的干扰。
- 相比传统的FDMA或TDMA技术，OFDMA具有更高的频谱效率和灵活性。
关键技术点：
- 子载波分配：根据用户的带宽需求和信道条件动态分配子载波。
- 调制与编码：选择合适的调制方式（如QPSK、16-QAM等）和编码方案以适应不同用户的信道条件。
- IFFT/FFT变换：实现时域信号与频域信号之间的转换。
- 循环前缀（Cyclic Prefix, CP）：用于消除符号间干扰（ISI）。
仿真目标：
- 搭建一个简化的OFDMA系统模型，涵盖信号生成、子载波分配、IFFT/FFT变换、AWGN信道及相应的解调过程。
- 分析不同信噪比条件下的误码率（BER）性能，并探讨子载波分配策略对系统性能的影响。

二、所需工具和环境

为了完成此OFDMA系统的仿真实现，你需要以下工具和环境：

MATLAB/Simulink：用于设计系统模型和运行仿真。
Communications Toolbox：提供必要的模块，如调制解调器、AWGN信道、IFFT/FFT等。
Signal Processing Toolbox（推荐）：用于实现信号处理算法。
5G Toolbox（可选）：对于更高级的功能实现，特别是针对现代无线通信标准的支持。

确保你已经安装了上述工具箱，并且拥有有效的许可证。

三、步骤详解

步骤1：创建Simulink模型

首先，在MATLAB中启动Simulink并创建一个新的空白模型。

matlab

深色版本

% 创建新的Simulink模型
modelName = 'OFDMA_Simulation';
new_system(modelName);
open_system(modelName);

步骤2：设计输入信号

我们需要两组随机比特流作为两个用户的输入信息比特。

在 Sources 库中拖拽两个 Random Integer Generator 模块到模型编辑区，分别表示两个用户的信号。
设置参数如下：
- M-ary number：根据所使用的调制方式设置（例如，对于QPSK，设置为4）。
- Initial seed：任意整数值（用于随机数生成器的种子）。
- Sample time：1/Fs（根据设定的采样频率）。

matlab

深色版本

% 设定参数
Fs = 1000; % 采样频率 (Hz)

% 添加随机整数生成器模块（用户1）
add_block('comm/Random Integer Generator', [modelName '/Random_Integer_User1']);
set_param([modelName '/Random_Integer_User1'], 'M', '4'); % 对于QPSK
set_param([modelName '/Random_Integer_User1'], 'InitialSeed', '73');
set_param([modelName '/Random_Integer_User1'], 'SampleTime', '1/1000');

% 添加随机整数生成器模块（用户2）
add_block('comm/Random Integer Generator', [modelName '/Random_Integer_User2']);
set_param([modelName '/Random_Integer_User2'], 'M', '4'); % 对于QPSK
set_param([modelName '/Random_Integer_User2'], 'InitialSeed', '97');
set_param([modelName '/Random_Integer_User2'], 'SampleTime', '1/1000');

步骤3：添加调制模块

选择适当的数字调制方式，这里以QPSK为例。

在 Communications Toolbox > Modulation > Digital Baseband Modulation 库中找到并拖拽两个 QPSK Modulator Baseband 模块到模型编辑区，分别对应两个用户。

matlab

深色版本

% 添加QPSK调制器模块（用户1）
add_block('comm/QPSK Modulator Baseband', [modelName '/QPSK_Modulator_User1']);

% 连接随机整数生成器到QPSK调制器
add_line(modelName, [modelName '/Random_Integer_User1'], [modelName '/QPSK_Modulator_User1'], 'autorouting', 'on');

% 添加QPSK调制器模块（用户2）
add_block('comm/QPSK Modulator Baseband', [modelName '/QPSK_Modulator_User2']);

% 连接随机整数生成器到QPSK调制器
add_line(modelName, [modelName '/Random_Integer_User2'], [modelName '/QPSK_Modulator_User2'], 'autorouting', 'on');

步骤4：实现子载波分配

在OFDMA中，每个用户被分配一组特定的子载波。我们可以通过矩阵操作来模拟这一过程。

使用 Matrix Concatenate 模块将两个用户的调制信号合并到频域。

matlab

深色版本

% 添加矩阵拼接模块
add_block('simulink/Math Operations/Matrix Concatenate', [modelName '/Subcarrier_Allocation']);

% 配置矩阵拼接模块
set_param([modelName '/Subcarrier_Allocation'], 'Mode', 'Horizontal'); % 水平拼接

% 连接QPSK调制器输出到子载波分配模块
add_line(modelName, [modelName '/QPSK_Modulator_User1'], [modelName '/Subcarrier_Allocation'], 'autorouting', 'on');
add_line(modelName, [modelName '/QPSK_Modulator_User2'], [modelName '/Subcarrier_Allocation'], 'autorouting', 'on');

步骤5：添加IFFT模块

OFDMA通过IFFT将频域信号转换为时域信号。

在 Communications Toolbox > Signal Operations 库中拖拽 IFFT 模块到模型编辑区。

matlab

深色版本

% 添加IFFT模块
add_block('comm/IFFT', [modelName '/IFFT']);

% 连接子载波分配模块到IFFT模块
add_line(modelName, [modelName '/Subcarrier_Allocation'], [modelName '/IFFT'], 'autorouting', 'on');

步骤6：添加循环前缀（CP）

为了消除符号间干扰（ISI），需要在OFDM符号前添加循环前缀。

使用 Buffer 和 Delay 模块实现循环前缀的插入。

matlab

深色版本

% 添加缓冲区模块
add_block('dsp/Buffer', [modelName '/Buffer']);

% 添加延迟模块
add_block('simulink/Discrete/Delay', [modelName '/Delay']);

% 配置缓冲区模块
set_param([modelName '/Buffer'], 'OutputBufferSize', '64'); % 假设子载波数量为64
set_param([modelName '/Buffer'], 'BufferOverlapLength', '16'); % 循环前缀长度为16

% 连接IFFT模块到缓冲区模块
add_line(modelName, [modelName '/IFFT'], [modelName '/Buffer'], 'autorouting', 'on');

步骤7：添加AWGN信道

为了模拟现实中的噪声环境，我们向叠加后的OFDMA信号中添加AWGN。

在 Communications Toolbox > Channels 库中拖拽 AWGN Channel 模块到模型编辑区。
设置参数如下：
- SNR (dB)：设置适当的信噪比以模拟不同的信道条件。

matlab

深色版本

% 添加AWGN信道模块
add_block('comm/AWGN Channel', [modelName '/AWGN_Channel']);
set_param([modelName '/AWGN_Channel'], 'SNR', '10'); % 根据实际情况调整SNR

% 连接缓冲区模块到AWGN信道
add_line(modelName, [modelName '/Buffer'], [modelName '/AWGN_Channel'], 'autorouting', 'on');

步骤8：接收端处理

在接收端，执行与发送端相反的操作，恢复原始信号。

移除循环前缀：
- 使用 Selector 模块去除循环前缀。
FFT变换：
- 在 Communications Toolbox > Signal Operations 库中拖拽 FFT 模块到模型编辑区。
分离子载波：
- 使用 Demux 模块将频域信号分离为各用户的信号。

matlab

深色版本

% 添加选择器模块
add_block('simulink/Signal Routing/Selector', [modelName '/Remove_CP']);

% 添加FFT模块
add_block('comm/FFT', [modelName '/FFT']);

% 添加Demux模块
add_block('simulink/Signal Routing/Demux', [modelName '/Demux']);

% 连接AWGN信道到接收端模块
add_line(modelName, [modelName '/AWGN_Channel'], [modelName '/Remove_CP'], 'autorouting', 'on');
add_line(modelName, [modelName '/Remove_CP'], [modelName '/FFT'], 'autorouting', 'on');
add_line(modelName, [modelName '/FFT'], [modelName '/Demux'], 'autorouting', 'on');

步骤9：添加解调模块

为了恢复原始数据信号，我们在接收端使用QPSK解调器。

在 Communications Toolbox > Modulation > Digital Baseband Modulation 库中找到并拖拽两个 QPSK Demodulator Baseband 模块到模型编辑区，分别对应两个用户。

matlab

深色版本

% 添加QPSK解调器模块（用户1）
add_block('comm/QPSK Demodulator Baseband', [modelName '/QPSK_Demodulator_User1']);

% 添加QPSK解调器模块（用户2）
add_block('comm/QPSK Demodulator Baseband', [modelName '/QPSK_Demodulator_User2']);

% 连接Demux模块到QPSK解调器
add_line(modelName, [modelName '/Demux/1'], [modelName '/QPSK_Demodulator_User1'], 'autorouting', 'on');
add_line(modelName, [modelName '/Demux/2'], [modelName '/QPSK_Demodulator_User2'], 'autorouting', 'on');

步骤10：计算误码率（BER）

为了评估整个系统的性能，我们需要计算每个用户的误码率（BER）。

在 Communications Toolbox > Measurements and Statistics 库中拖拽两个 Error Rate Calculation 模块到模型编辑区。

matlab

深色版本

% 添加误码率计算模块（用户1）
add_block('comm/Error Rate Calculation', [modelName '/Error_Rate_User1']);

% 连接原始信号和最终解调信号到误码率计算模块
add_line(modelName, [modelName '/Random_Integer_User1'], [modelName '/Error_Rate_User1/Input Port']);
add_line(modelName, [modelName '/QPSK_Demodulator_User1'], [modelName '/Error_Rate_User1/Output Port'], 'autorouting', 'on');

% 添加误码率计算模块（用户2）
add_block('comm/Error Rate Calculation', [modelName '/Error_Rate_User2']);

% 连接原始信号和最终解调信号到误码率计算模块
add_line(modelName, [modelName '/Random_Integer_User2'], [modelName '/Error_Rate_User2/Input Port']);
add_line(modelName, [modelName '/QPSK_Demodulator_User2'], [modelName '/Error_Rate_User2/Output Port'], 'autorouting', 'on');

步骤11：添加显示模块

为了观察每个用户的误码率变化情况，我们需要添加两个输出模块。

在 Sinks 库中拖拽两个 Display 模块到模型编辑区。

matlab

深色版本

% 添加显示模块（用户1）
add_block('simulink/Sinks/Display', [modelName '/Display_User1']);

% 连接误码率计算模块到显示模块
add_line(modelName, [modelName '/Error_Rate_User1'], [modelName '/Display_User1'], 'autorouting', 'on');

% 添加显示模块（用户2）
add_block('simulink/Sinks/Display', [modelName '/Display_User2']);

% 连接误码率计算模块到显示模块
add_line(modelName, [modelName '/Error_Rate_User2'], [modelName '/Display_User2'], 'autorouting', 'on');

步骤12：设置仿真参数

根据需要调整仿真时间、求解器类型和其他相关参数。

matlab

深色版本

% 设置仿真参数
set_param(modelName, 'StopTime', '10'); % 模拟运行时间为10秒
set_param(modelName, 'Solver', 'Fixed-step'); % 设置固定步长求解器
set_param(modelName, 'FixedStep', '1/1000'); % 设置固定步长为采样周期