Magic Road近两年做毫米波通信仿真的人越来越多，但一说到“大规模MIMO的混合波束成形”，很多人第一反应是翻论文、看公式，然后卡在代码实现上。尤其是发射端那一整套流程——从信道建模、模拟预编码到数字预编码的设计，再到最后的频谱效率评估，每一步都有不少坑。我去年接手一个相关课题，用Matlab从零搭了一套完整的发射端混合波束成形仿真框架，跑了上千次蒙特卡洛实验，踩了不少坑，今天把这套东西完整梳理出来。
这篇文章适合正在做5G/6G物理层算法验证的同学，或者刚接触毫米波大规模MIMO、想快速跑通一个混合波束成形仿真链路的研究生。文章会从系统模型、算法原理、代码架构到参数调试，把发射端混合波束成形的Matlab实现完整拆开讲，代码思路可以直接迁移到自己的项目里。
1. 大规模MIMO混合波束成形：问题是什么，为什么要这么做
1.1 从全数字波束成形到混合波束成形的必然选择
传统MIMO系统里，每根天线后面都跟着一条完整的射频链路（包括DAC、混频器、功放等），这就是所谓的全数字波束成形。当天线数量达到64、128甚至256根时，全数字方案的成本和功耗会成倍增长——每一根天线都要配备一套独立的射频硬件，这在基站端几乎不可接受。
混合波束成形的思路很简单：把波束成形拆成两级。第一级是模拟波束成形，用移相器网络在射频域实现，只改变信号相位不改变幅度，硬件成本低；第二级是数字波束成形，在基带完成，可以同时调节幅度和相位，灵活性高。这样射频链路数量可以从“天线数量”降到“数据流数量”，比如64根天线只需要8条射频链路，硬件成本大幅下降。
这就像一个大型会议室里，模拟波束成形相当于先把麦克风对准某个方向区域（粗调），数字波束成形则是在这个区域内更精细地分辨不同人的声音（细调）。两级配合，既保证了波束指向的准确性，又极大降低了硬件复杂度。
1.2 发射端混合波束成形的系统模型
本文采用的是窄带块衰落信道模型，发射端配置 (N_t) 根天线、(N_{RF}) 条射频链路，同时传输 (N_s) 路数据流，满足 (N_s \le N_{RF} \le N_t)。发射信号经过数字预编码矩阵 (F_{BB})（维度 (N_{RF} \times N_s)）、模拟预编码矩阵 (F_{RF})（维度 (N_t \times N_{RF})）后，由天线发射出去。
接收端同样配置 (N_r) 根天线和 (N_{RF}^r) 条射频链路，接收信号先经过模拟合并矩阵 (W_{RF})，再经过数字合并矩阵 (W_{BB})，得到最终的估计信号。整个链路的输入输出关系可以写成：
[
y = W_{BB}^H W_{RF}^H H F_{RF} F_{BB} s + W_{BB}^H W_{RF}^H n
]
其中 (s) 是 (N_s \times 1) 的发射符号向量，(H) 是 (N_r \times N_t) 的信道矩阵，(n) 是高斯白噪声。这里有一个关键约束：模拟预编码矩阵 (F_{RF}) 的每个元素必须满足恒模约束，即 (|F_{RF}(i,j)| = 1/\sqrt{N_t})，这是由移相器只能改变相位的物理特性决定的。
这个约束是混合波束成形设计和传统全数字波束成形最大的区别所在——你不能随意设计 (F_{RF})，它的每个元素都被锁死在单位圆上。整个算法的核心目标就是在这个约束下，找到最优的 (F_{RF}) 和 (F_{BB})，使得系统频谱效率最大化。
需要模型API调用？ 免费领10W Token，多模型网关一键接入 Claude、DeepSeek 等主流模型。
2. 算法设计思路：如何一步步逼近最优解
2.1 目标函数与问题建模
混合波束成形的设计目标很明确：最大化系统频谱效率。在完美信道状态信息已知的条件下，可达频谱效率可以写成：
[
R = \log_2 \det\left( I_{N_s} + \frac{\rho}{N_s} R_n^{-1} W_{BB}^H W_{RF}^H H F_{RF} F_{BB} F_{BB}^H F_{RF}^H H^H W_{RF} W_{BB} \right)
]
其中 (\rho) 是平均接收功率，(R_n) 是噪声协方差矩阵。
直接优化这个目标函数非常困难，因为 (F_{RF}) 的恒模约束让整个问题变成了非凸优化问题。学术界的通用做法是把它拆成两个子问题：先设计模拟预编码器 (F_{RF})，再在固定 (F_{RF}) 的前提下设计数字预编码器 (F_{BB})。这种次优解虽然不保证全局最优，但实际性能可以非常接近全数字方案。
2.2 模拟预编码设计：码本逼近法
模拟预编码的设计主要有两种思路。第一种是基于码本的波束扫描方法，预先设计一组覆盖不同方向的波束矢量，然后选择最匹配信道的一组；第二种是基于信道协方差矩阵的特征分解，提取主特征方向作为模拟波束方向。
我采用的是码本逼近法，这是实际工程中最常用也最容易实现的方法。码本设计基于发射天线阵列的响应向量，对于 (N_t) 根天线的均匀线性阵列，第 (k) 个码字定义为：
[
a_t(\theta_k) = \frac{1}{\sqrt{N_t}} [1, e^{j\pi \sin(\theta_k)}, e^{j2\pi \sin(\theta_k)}, \ldots, e^{j(N_t-1)\pi \sin(\theta_k)}]^T
]
其中 (\theta_k) 是离散化的离开角，取值范围从 (-\pi/2) 到 (\pi/2)，均匀划分为 (K) 个角度。码本大小 (K) 直接决定了波束指向的精度，(K) 越大波束越密，但搜索开销也越大。
具体做法是：对信道矩阵 (H) 做奇异值分解 (H = U\Sigma V^H)，取右奇异矩阵 (V) 的前 (N_s) 列作为最优全数字预编码矩阵 (F_{opt})，然后在码本中逐个搜索，选出与 (F_{opt}) 相关性最