如何在FPGA中实施4G无线球形检测器？

宽带互联网接入遇到的最大的挑战之一就是移动性，而这正是最新的 WiMAX 标准所要解决的。IEEE 802.16e-2005 介绍了传输和接收过程中多根天线的用法，即 MIMO 概念，又称为多输入多输出，是移动 WiMAX 的一个关键特性。

空分复用 (SDM) MIMO 处理可显著提高频谱效率，进而大幅增加无线通信系统的容量。空分复用 MIMO 通信系统作为一种能够大幅提升无线系统容量和连接可靠性的手段，近来吸引了人们的广泛关注。

MIMO 无线系统最佳硬判决检测方式是最大似然 (ML) 检测器。ML 检测因为比特误码率 (BER)性能出众，非常受欢迎。不过，直接实施的复杂性会随着天线和调制方案的增加呈指数级增强，使 ASIC 或 FPGA 仅能用于使用少数天线的低密度调制方案。

在 MIMO 检测中，既能保持与最佳 ML 检测相媲美的 BER 性能，又能大幅降低计算复杂性的出色方法非球形检测法莫属。这种方法不仅能够降低 SDM 和空分多接入系统的检测复杂性，同时又能保持与最佳 ML 检测相媲美的 BER 性能。实现球形检测器有多种方法，每种方法又有多种不同算法，因此设计人员可以在诸如无线信道的吞吐量、BER 以及实施复杂性等多项性能指标之间寻求最佳平衡。

虽然算法（比如 K-best 或者深度优先搜索）和硬件架构对 MIMO 检测器的最终 BER 性显而易见有极大的影响，不过一般在球形检测之前进行的信道矩阵预处理也会对 MIMO 检测器的最终 BER 性能产生巨大影响。信道矩阵预处理可繁可简，比如根据对信道矩阵进行的方差计算结果 (variance computation)，计算出处理空分复用数据流的优先次序，也可以使用非常复杂的矩阵因子分解方法来确定更为理想（以 BER 衡量）的数据流处理优先次序。

Signum Concepts 是一家总部位于圣地亚哥的通信系统开发公司，一直与赛灵思和莱斯大学（Rice University）开展通力合作，运用 FPGA 设计出了用于 802.16e 宽带无线系统的空分复用MIMO 的MIMO 检测器。该处理器采用信道矩阵预处理器，实现了类似贝尔实验室分层空时 (BLAST)结构上采用的连续干扰抵消处理技术，最终达到了接近最大似然性能。

系统考虑因素
理想情况下，检测过程要求对所有可能的符号向量组合进行 ML 解决方案计算。球形检测器旨在通过使用简单的算术运算降低计算复杂性，同时还能够保持最终结果的数值完整性。我们的方法，第一步是把复杂的数值信道矩阵分解为只有实数的表达式。这个运算增加了矩阵维数，但简化了处理矩阵元的计算。降低计算复杂性的第二个方面体现在，减少检测方案分析和处理的可选符号。其中，对信道矩阵进行 QR 分解是至关重要的一步。

图 1 显示的是如何进行数学转换，得出计算部分欧几里德距离度量法的最终表达式。欧几里德距离度量法是球形检测过程的基础。R代表三角形矩阵，用于处理以矩阵元 rM,M 开始的可选符号的迭代法。其中，M代表信道矩阵以实数表达的维数。该解决方案通过 M 次迭代定义出遍历树结构，树的每层i对应第i根天线的处理符号。

图 1. 用于球形检测器 MIMO 检测的部分欧几里德距离度量方程


球形检测器处理天线的次序对 BER 性能有着极大的影响。因此，在进行球形检测前，我们的设计采用了类似于 V-BLAST 技术的信道重新排序技术。

实现树的遍历有几种可选方法。在我们的实施方案中，则使用了广度优先搜索法，这是因为该方法采用备受欢迎的前馈结构，因此具有硬件友好特征。在每一层，该实施方案只选择K 个距离最小的幸存节点来计算扩展情况。

球形检测器处理天线的次序对 BER 性能有着极大的影响。因此，在进行球形检测前，我们的设计采用了类似于 V-BLAST 技术的信道重新排序技术。

该方法通过多次迭代，计算出信道矩阵的伪逆矩阵的行范数，然后确定信道矩阵最佳列检测次序。根据迭代次数，该方法可以选择出范数最大或者最小的行。欧几里德范数最小的逆矩阵行表示天线的影响最强，而欧几里德范数最大的行则表示天线的影响最弱。这种新颖的方法首先处理最弱的数据流，随后依次迭代处理功率从高到低的数据流。

FPGA 硬件应用
为实现上述系统，我们采用了赛灵思 Virtex®-5 FPGA 技术。该设计流程采用赛灵思 System Generator 进行设计捕获、仿真和验证。为了支持各种不同数量的天线／用户和调制次序，我们将检测器设计用于要求最高的 4x4、64-QAM 情况下。

我们的模型假定接收方非常清楚信道矩阵，这可以通过传统的信道估算方法来实现。在信道重新排序和 QR 分解之后，我们开始使用球形检测器。为准备使用软输入、软输出信道解码器（比如 turbo 解码器），我们通过计算检测到的比特的对数似然比 (LLR) 来生成软输出。

该系统的主要架构元素包括数据副载波处理和系统子模块管理功能，以便实时处理所需数量的子载波，同时最大程度地降低处理时延。对每个数据副载波都进行了信道矩阵估算，限定了每个信道矩阵可用的处理时间。对选中的 FPGA 而言，其目标时钟频率为 225MHz，通信带宽为 5MHz（相当于 WiMAX 系统中的 360 个数据子载波），每个信道矩阵间隔可用的处理时钟周期数为 64。

我们采用硬件功能单元精湛的流水线和时分复用 (TDM) 功能，以达到 WiMAX OFDM 符号的实时要求。

除了高数据率外，在架构设计指导过程中控制子模块时延也是一个重要的问题。我们通过引入连续信道矩阵的 TDM 解决了时延问题。这种方法可以延长同一信道矩阵元之间的处理时间，同时还能保持较高的数据吞吐量。构成 TDM 组的信道数会随着子模块的不同而变化。在 TDM 方案中，信道矩阵求逆过程用了 5 个信道，而有 15 个信道在实数 QR 分解模块中进行了时分复用。图 2 是该系统的高级流程图。


图 2. MIMO 802.16e 宽带无线接收器的高级流程图

信道矩阵预处理

信道矩阵预处理器确定了空分复用复合信号每一层的最佳检测次序。该预处理器负责计算信道矩阵的伪逆矩阵范数，并根据这些范数，选择待处理的下一个传输流。伪逆矩阵中范数最小的行对应着最强传输流（检波后噪声放大最小），而范数最大的行对应着质量最差的层（检波后噪声放大最大）。我们的实施方案首先检测最弱的层，然后按最低噪声放大到最高噪声放大的次序逐层检测。对排序过程中的每一步，信道矩阵中相应的列随后会被清空，然后简化后的矩阵进入下一级的天线排序处理流水线。

在预处理算法中，伪逆矩阵的计算要求最高。这个过程的核心是矩阵求逆，通常通过吉文斯(Givens) 旋转进行 QR 分解 (QRD) 来实现。常用的角度估算和平面旋转算法（如 CORDIC）会造成严重的系统时延，对我们的系统来说是不可接受的。因此，我们的目标是运用 FPGA 的嵌入式 DSP 资源（比如 Virtex-5 器件中的 DSP48E），找出矢量旋转和相位估算的替代性解决方案。

QRD 的脉动阵列结构由两种类型的处理单元构成——对角线单元或边界单元和非对角线单元或内部单元。边界单元执行矢量函数，可以生成阵列内部单元使用的旋转角度。要想得到想要的旋转角度，可以把非对角线单元中的值与对角线单元中的共轭复数相乘，然后除以复数的倒数即可。相除实际是用乘法的方式完成的，即在观察到函数接近线性的时候，乘以根据定义的间隔的多项式近似值计算出的倒数。图 3 显示了采用这种近似值在对角线脉动单元中完成这种复杂旋转的信号流程图。


图 3. 对角线脉动单元结构图

发送到非对角线单元中的数据是旋转矢量的同相部分和正交部分除以相应的近似值得出的结果。我们不仅通过在对角线单元和非对角线单元采用流水线架构实现了高数据吞吐量，同时还通过对跨5个信道的硬件进行时分复用的方式控制了近似值模块和复杂乘法器引起的时延。

对 4x4 矩阵，我们使用了 1 个对角线单元和 7 个非对角线单元。分解单个矩阵所花的处理时间为 4x4=16 个数据周期，而该设计交付数据的速度是每三个时钟周期一个样本，因此分解单个矩阵的所用的总时长为 3x4x4=48 个时钟周期（低于可用的 64 个时钟周期）。我们对分解后的矩阵使用了回代法（back substitution），同时以相同的 TDM 方式进一步进行了重新排序操作。

球形检测器
球形检测器采用PED 单元进行范数计算。根据树的层次，我们采用了三种不同类型的 PED 单元。根节点 PED 模块负责计算所有可能的 PED。二级 PED 模块针对上一级计算得出的 8 个幸存路径计算出 8 个可能的 PED。这样在树的下一级索引中，我们就有 64 个生成的 PED。第三种类型的PED模块用于其它树级，负责计算上一级计算出的所有 PED 的最邻近的节点 PED。

球形检测器 (SD) 的流水线架构可以在每个时钟周期中处理数据。其结果就是树的每级只需要一个 PED 模块。因此，对 4x4 64-QAM 系统而言，PED 单元的总数为 8，与树的级数相等。

SD 可以采用硬解码和软解码两种类型的解码技术。硬解码能够用贯穿树的各级的最小距离矩阵度量次序；软解码，用对数似然比来代表输出的每个比特。对数似然比一般被当作优先输入值提供给信道解码器，比如 turbo 解码器。

FPGA资源占用
实施和仿真包括图 2 所示的检测过程，但不包括软输出生成模块。目标芯片是 Virtex-5 XC5VFX130T-2FF1738 FPGA。设计的时钟频率为 225MHz，可用的数据率为 83.965Mb/s。

表 1 显示了设计中每个主要功能单元的资源占用情况。利用率 (%) 表示 FPGA 面积占XC5VFX130T 器件总面积的百分比。


功能
切片数
LUTs/FFs
DSP48
块RAM
信道预处理
9,999
(48%)
20,339/29,954
(24%)
159
(49%)
105
(17%)
RVD QRD
1,715
(8%)
4,418/5,556
(5%)
30
(9%)
27
(4%)
球形检测器
2,445
(11%)
3,113/6,525
(3%)
48
(15%)
12
(2%)

表 1. 按子系统划分的资源占用情况

图 4. 4x4 64-QAM的浮点 MATLAB 仿真（硬判决）、System Generator设计（硬判决）BER 曲线与最大似然曲线相比

System Generator 和基于模型的设

我们使用针对 DSP 设计流程的赛灵思 System Generator 实现了完整的硬判链。设计验证工作不仅使用了 MATLAB®/Simulink® 环境的仿真语义，还有 System Generator 的协同仿真功能。信道矩阵参数的同相部分和正交部分从正常的分布得出，并由 MATLAB 交付给System Generator 建模环境。我们同样使用这种仿真框架进行了比特误码率计算。图 4 对我们的定点硬判决设计 BER 曲线、浮点硬判决设计BER曲线和最佳 ML 参考曲线进行了比较。我们通过对赛灵思 ML510 开发平台进行基于以太网的硬件协仿真，开发出了该设计的硬件演示。信道矩阵参数采用赛灵思 AWGN IP核发送给球形检测器。我们通过把设计嵌入到自同步 BER 测试器来计算 BER。该仪器能够向检测器发送输入并捕获误码。

本文就采用空分复用 MIMO 的通信系统使用的球形检测器进行了简要介绍。我们详细探讨了球形检测器和信道矩阵预处理器的架构情况。实现预处理的方法有许多种，虽然我们的方法在计算上要复杂一点，但得出的 BER 性能接近最大似然。虽然我们的讨论是围绕 WiMAX 进行的，设计人员可以把其中的许多方法用于 3G LTE（长期演进）无线系统。

我们团队的下一步工作是通过采用 turbo 卷积码和软输出生成模块执行迭代软检测来改善BER 性能。