企业官网建设流程全解析

简介：WINNER II模型是无线通信领域评估3D-MIMO系统性能的标准化框架，这份MATLAB源代码包面向无线通信研究人员与研究生，可应用于信道建模、天线阵列设计与系统级仿真。资源共52个文件，其中44个m脚本覆盖信道参数生成、路径损耗计算、天线响应等核心流程，另有p函数、mlx脚本及说明文档，压缩包整体仅202KB，便于快速部署和二次开发。已有487人学习下载。代码内含三个实例，结合LOS/NLOS传播、空间分集与复用、LMMSE信道估计、ZF/MMSE检测等知识点，展示从参数配置到误码率与吞吐量评估的完整链路。通过研读这些源码，读者可掌握3D-MIMO信道构建方法，理解天线阵列与随机几何场景设置，为学术研究或工程项目提供可直接运行的参考实现。

1. WINNER II模型源代码：3D MIMO信道仿真从哪一步开始

手头这份“WINNER II模型源代码.zip_3D”，解压出来是一整套MATLAB信道仿真工程。对做MIMO算法验证的人来说，它的价值不是“又一个能生成随机矩阵的脚本”，而是唯一能把多径时延、离开/到达方位角、仰角、极化和三维坐标同时放进一个信道系数矩阵的成熟开源实现。直接拿高斯随机信道测波束成形，性能曲线永远漂亮，一上外场就翻车；WINNER II给每个多径分量配上物理坐标和传播角度，仿真结果才具备和实测对照的意义。这篇笔记从模型原理、源码调用、参数调整到踩坑排查，把整套代码讲透。

2. WINNER II 3D信道模型：为什么基站天线仰角一改，信道容量就跟着变

2.1 多径不是一条路径，而是一簇带角度的散射体

WINNER II是典型的几何随机信道模型（GSCM）。它不像瑞利模型那样直接给一条复高斯序列，而是在基站和终端之间铺满散射体簇（cluster），每个簇里再分布若干条路径（path）。每一条路径都携带六个维度的信息：时延、功率、出发方位角AoD、到达方位角AoA、出发仰角ZoD、到达仰角ZoA。所谓3D版本，核心就是把后两个仰角维度补上了。

这个补齐的意义在MIMO里非常直接。传统2D信道模型只输出方位角，天线阵列里的垂直维度在仿真里是“看不到信号”的。当你用8×4或16×8的平面阵做波束成形，垂直维度能不能区分用户，取决于信道里有没有ZoD/ZoA。WINNER II的3D实现就是解决这个信息缺口，让基站侧仰角扫描、垂直复用这些操作在仿真里能真实反映到容量曲线和误码率曲线上。

模型内部把参数分成两层。大尺度参数LSP包括路径损耗、阴影衰落、Ricean K因子，决定平均接收功率的走势；小尺度参数SSP负责簇数、每簇功率、时延和四个角度，决定空间相关性的细节。LSP通过协方差矩阵在位置上做相关映射，SSP则按场景指定的概率分布做随机抽样。可以用一句话概括整个流程：位置决定场景，场景决定参数分布，参数分布决定每一条多径的空间特征。

2.2 场景编号不是随便选的：A1到D1背后是整套参数表

WINNER II把无线环境分成十余个场景，从A1室内办公室到D1农村宏蜂窝，每个场景都有自己独立的参数表。参数表里规定了路径损耗指数、时延扩展均值、角度扩展均值、阴影衰落标准差、LoS概率曲线等。出错最多的不是算法代码，而是场景和部署对不上。

B1和C2是3D MIMO仿真里最常用的两个起点。微蜂窝场景路径损耗指数大约是2.1到2.5，NLoS下宏蜂窝可能到3.5甚至4.5，同一组天线在不同场景里容量能差一个数量级。

2.3 从参数表到信道系数：WINNER II生成H矩阵的六个步骤

拿到源码包以后，先理解它内部处理顺序，排错时就不至于摸黑：

1. 配置基站和终端的二维坐标或三维坐标，挂上天线阵列结构。
2. 按场景编号加载大尺度参数表，计算路径损耗、阴影衰落和LoS/NLoS状态。
3. 从小尺度参数分布中抽样出簇数、每簇功率、簇时延、方位角AoA/AoD和仰角ZoA/ZoD。
4. 用极化耦合矩阵和天线方向图对每条路径做加权。
5. 叠加多普勒频移和初相，生成随时间变化的快照。
6. 按分组输出信道系数矩阵H。

这套顺序在绝大多数WINNER II MATLAB实现里都成立。搞清楚第3步和第4步的产物分别是什么，后面调参数、查NaN都有据可依。

2.4 为什么不直接用瑞利、COST207或38.901：选型对比

既然市面上有这么多信道模型，为什么要用WINNER II的源码包？按实际仿真目的一一对比：

3GPP 38.901在很多思路上延续了WINNER II，但整体实现更重，适合标准级验证。WINNER II的MATLAB源码足够轻、场景参数开放可改，做算法预研和论文复现是最平衡的选择。

3. 把MATLAB版WINNER II跑起来：最小脚本、路径配置与输出核对

3.1 解压后先看清这三类文件再去跑demo

拿到源码包后别急着双击main，先花两分钟看目录结构。一个完整的WINNER II MATLAB工程里，通常有三类内容：配置与场景脚本（定义布局、场景编号、天线参数）、核心信道传播函数（完成LSP/SSP抽样和H矩阵生成）、场景参数表（通常是.mat文件，存放各场景的统计参数）。

3D扩展版还会额外带上仰角参数相关的抽样函数和垂直阵列初始化脚本。核心函数名在不同zip包里有差异，可能是propagate_channels也可能是winner2_channel，这不重要，重要的是先找到一个能跑的demo脚本，把demo里的场景参数一行一行改成自己需要的，比从零调接口稳妥得多。

3.2 最小调用脚本：三步生成3D MIMO信道系数

下面给出一个最小可用的脚本框架。函数名以你手里的zip实际为准，但流程是通用的：

% winner2_3d_min.m % 最小WINNER II 3D MIMO信道生成脚本 % 适用：带3D扩展的WINNER II MATLAB源码包 clear; close all; clc; % 1) 把源码包根目录和所有子目录挂进MATLAB路径 % 参数表.mat通常在子目录里，只加根目录会报“找不到文件” addpath(genpath('D:\work\WINNER2_3D')); % 2) 定义布局：坐标单位是米，三列对应x/y/z bs_pos = [0, 0, 25]; % 基站：原点，高25米 ms_pos = [50, 20, 1.5; % 终端1：水平距离约54米 -30, 40, 1.5]; % 终端2：另一侧街道 % 3) 场景与传播参数 cfg.scenario = 'B1'; % 城市微蜂窝 cfg.use_3d = 1; % 打开3D模式，生成ZoD/ZoA cfg.center_freq = 2.6e9; % 中心频率，注意单位是Hz cfg.num_paths = 12; % 每条簇内路径数 cfg.num_snapshots = 100; % 时间快照数 cfg.num_ms = 2; % 终端数量 cfg.polarization = 1; % 1=双极化，0=单极化 % 4) 天线阵列：基站8行4列平面阵，终端2元线阵 bs_ant.shape = 'URA'; % 平面阵列 bs_ant.n_row = 8; % 水平方向阵元数 bs_ant.n_col = 4; % 垂直方向阵元数，3D MIMO关键 bs_ant.row_spacing = 0.5; % 单位：波长 bs_ant.col_spacing = 0.7; % 垂直间距，仰角分辨率相关 bs_ant.tilt_deg = 10; % 下倾角 ms_ant.shape = 'ULA'; ms_ant.n_ant = 2; % 5) 调用传播主函数（函数名以包内实际为准） [H, lsp, ssp] = propagate_channels(bs_pos, ms_pos, cfg, bs_ant, ms_ant); % 6) 核对输出维度 fprintf('H维度: %d x %d x %d x %d x %d\n', size(H));

这段脚本里最容易被忽略的是cfg.center_freq的单位。写成2.6而不是2.6e9，路径损耗公式里的频率项会直接出现数量级错误，但不会报错，属于最隐蔽的一类问题。cfg.use_3d控制的是底层是否执行仰角抽样，建议设成1以后打印一下ssp结构，确认ZoD字段确实有值，才算真正进入3D模式。

3.3 输出数组的维度与前两个核对点

H矩阵的返回维度在不同实现里有两种常见约定，一种直接返回五维数组，另一种把快照维放在最后一个位置。拿到输出后先做两个核对：

% 核对1：平均功率是否合理 P_est = mean(abs(H(:)).^2); fprintf('平均功率: %.3f\n', P_est); % 归一化实现里0.5~1.5都正常 % 核对2：终端维是否正确 n_ms_estimated = size(H, 5); % 或按实际维度位置调整 fprintf('终端数: %d\n', n_ms_estimated);

平均功率如果显示0.001以下，大概率是路径损耗和center_freq单位打架；如果出现Inf，先查坐标是否有重合。这两个核对点能挡住后面一大半的“玄学”问题。

3.4 MATLAB版本与运行时的依赖问题

源码包年代较早，遇到新版本MATLAB报legend、set之类的函数警告是常态，不影响结果。真正影响运行的是两点：一是机器学习或通信工具箱版本不同导致某些内置函数签名变化，这时优先把源码包里自带的那份同名函数路径放在最前面；二是MATLAB递归调用限制，部分实现用递归做树形散射体搜索，报“Maximum recursion limit”时执行set(0,'RecursionLimit', 1000)即可。

4. 必调的六个关键参数：场景编号、天线阵列与极化设置

4.1 场景编号：B1/C2/D1背后的路径损耗差一个数量级

场景编号是第一个要动的参数，也是影响结果最粗暴的参数。做3D MIMO波束验证，通常从B1和C2二选一起步。B1是微蜂窝，基站低于屋顶，Lo S概率高，适合验证垂直波束覆盖街道场景；C2是宏蜂窝，基站高于屋顶，NLoS占比高，适合验证大规模天线阵列的空间复用能力。

选择依据可以总结成三条：基站相对屋顶的高度、终端分布密度、Lo S占比需求。B1的Lo S概率在中近距离明显高于C2，如果算法号称能抗NLoS场景，直接用B1测会高估性能。D1这类农村宏蜂窝场景多径稀疏，阵列增益更容易做出来，论文里常用，但和实际城市部署距离较远。

4.2 中心频率与3D开关：不写对这两个，仰角扩展等于白设

中心频率直接影响路径损耗计算，同时决定阵元间距换算成物理尺寸时的基准。源码包默认频率可能是2.45GHz，做3.5GHz频段项目时只改cfg.center_freq不够，还需要确认阵元间距的单位——如果row_spacing填的是0.5λ，而源码用固定米数计算，频率一改实际电尺寸就变了。

cfg.use_3d这个开关容易被误以为“开了就有仰角”。它只是让底层生成ZoD/ZoA，如果天线阵列本身是水平面单排阵元，仰角信息拿到也用不上。3D MIMO仿真里，use_3d必须和平面阵同时生效，缺一个等于白设。

4.3 天线阵列参数：垂直阵元数和阵元间距决定仰角分辨率

这是3D MIMO仿真里最值得花时间调的参数块。

% 天线阵列参数调整：8x4平面阵，垂直方向间距0.7波长 bs_ant.shape = 'URA'; bs_ant.n_row = 8; % 水平覆盖 bs_ant.n_col = 4; % 垂直覆盖，观测仰角的分辨率 bs_ant.row_spacing = 0.5; bs_ant.col_spacing = 0.7; % 垂直间距大于水平间距，压低旁瓣 bs_ant.tilt_deg = 12; % 下倾角，宏蜂窝场景必调

垂直阵元数n_col从1改成4，信道矩阵的空间相关性结构会明显变化，垂直维度开始出现独立的特征方向。col_spacing建议从0.5λ起步，密度高的小区可以放到1λ，换取更窄的垂直波束。这里有个反直觉的点：垂直阵元间距越大，波束越窄，但栅瓣风险越高，调参时要把方向图旁瓣一起盯住，不能只看容量。

4.4 极化、下倾角和快照数：隐藏的“仿真公平性”三件套

这三个参数最容易在对比算法时被忽略，导致结论不公。极化开关cfg.polarization设为1后，同一场景下容量通常比单极化低一截，因为极化交叉耦合把一部分功率打散了。公平对比时，所有算法必须在同一极化配置下跑。

下倾角tilt_deg对3D MIMO的影响比很多人以为的大。在C2宏蜂窝场景里，下倾角从0°改到15°，小区边缘用户的容量可以提升20%以上，同时中心用户可能受损。仿真时如果只报平均容量，调一个有利的下倾角会让算法排名反转。快照数num_snapshots属于统计稳定性参数，建议至少500个起步，少于100个的容量曲线会有明显抖动，出现“这版结果和上一版差很多”的情况，先看快照数而不是怀疑代码。

5. WINNER II常见问题排查：路径报错、内存溢出与NaN信道系数

5.1 函数找不到：addpath(genpath())不够，还有大小写与子目录

现象：运行脚本，第一行就报Unrecognized function or variable 'xxx'，或者在加载参数表时报Unable to open file。

原因：源码包目录层级深，核心函数和参数表.mat分散在多个子目录里，只addpath根目录覆盖不到；另一种可能是MATLAB对大小写敏感，源码文件名里用了大写。

解决：统一用addpath(genpath('源码包根目录'))挂载全部子目录，并把这段写进startup.m避免每次重跑。报错时用which 函数名确认解释器找到的是哪一个路径，如果是重名的MATLAB内置函数，把源码包目录移到PATH最前面。

5.2 终端坐标和配置不一致：输出信道少了一路

现象：H矩阵里终端维长度是1，但明明配置了两个终端。

原因：源码里H的终端维由ms_pos的行数直接决定，cfg.num_ms只做辅助声明。改动ms_pos时忘了同步cfg.num_ms，或者反过来。

解决：以坐标矩阵为准，每行一个终端。检查size(ms_pos,1)和cfg.num_ms是否一致，再把终端坐标打印出来核对。

5.3 快照数一多就Out of Memory：single精度与分段生成

现象：num_snapshots设到10000，跑一会儿MATLAB直接报Out of memory。

原因：H矩阵存的是复数double，8×4×2×12×10000快照，算下来单个数组接近几十GB量级，内存扛不住。

解决：把输出精度降为单个浮点，或者分段生成。分段做法是把快照切成每批500到1000个，跑完一批只保留统计量，不保留原始H。这个处理比调大MATLAB内存上限可靠得多，内存设得再大也经不住几组快照翻倍。

5.4 信道系数出现NaN：坐标重合与单位写错

现象：H里大量NaN或Inf，平均功率没法看。

原因：基站和某个终端坐标重复导致距离为0，路径损耗公式求倒数直接Inf；或者center_freq写成2.6而不是2.6e9，频率项偏差太大；参数表.mat没加载成功也会让抽样函数拿到空数据。

解决：逐一检查bs_pos和每行ms_pos的欧氏距离，必须大于一个波长；把center_freq显式写成2.6e9，运行前打印确认；参数表加载失败时，看MATLAB工作区里变量是否为空。

注意：如果坐标和单位都没问题，H里还是有零星NaN，检查随机种子是否使用了极端值，把播种语句挪到外层，观测结果是否跟随种子分布，属于抽样端的随机性问题。

5.5 3D模式开了但容量曲线没变化：垂直阵元数才是仰角开关

现象：use_3d=1也设了，天线阵列也是面阵，但垂直维扫描容量曲线几乎没有变化，甚至重合。

原因：阵列配置里垂直阵元数n_col被写成了1，底层实际跑的还是一维线阵；或者col_spacing设置的间距和仰角扩展的量级不匹配，垂直方向相关性太高，区分不出角度。

解决：确认bs_ant.n_col>=4，col_spacing从0.5λ开始调整；打印ssp里的ZoD数组，看仰角分布是否随快照变化。ZoD字段如果一直是固定值，说明统计抽样没真正运行，回查参数表加载和3D开关的优先级。

6. 验证WINNER II输出可信的三个技巧：相关性、容量与对照

6.1 平均功率与天线相关性检查

信道矩阵不是黑匣子，拿到手先做两个快速检查。平均功率应该落在模型归一化设计的量级内，通常0.5到1.5之间，偏离三个数量级就要回头查单位。随后做天线相关性检查：取参考阵元的信道系数向量，依次计算它与各阵元的复相关系数，相关系数应随阵元间距增大而下降，垂直方向的相关性周期由col_spacing决定，如果出现水平垂直一条曲线，阵列配置多半写错了。

6.2 瞬时容量与理论对照

容量是最常用的综合指标，固定SNR算瞬时容量，然后和i.i.d.瑞利信道的容量上界做对照。WINNER II信道因为有空间相关性和多径簇结构，容量应当低于独立同分布瑞利上界，如果高出很多，说明信道抽样的空间信息没有进入H矩阵。

% 验证脚本：拿到H以后先别急着喂算法 % 1) 平均功率核对 P_est = mean(abs(H(:)).^2); fprintf('平均功率: %.3f\n', P_est); % 2) 瞬时容量抽查，固定SNR=10dB snr_lin = 10^(10/10); H_flat = reshape(H, [], size(H,3), size(H,4)); H_eff = squeeze(mean(H_flat, 1)); % 对多径维平均 [n_tx, n_rx] = size(H_eff); C = log2(det(eye(min(n_tx,n_rx)) + ... (snr_lin/n_tx) * (H_eff' * H_eff))); fprintf('瞬时容量: %.2f bps/Hz\n', real(C)); % 3) 相关性：相邻阵元系数对比 h0 = squeeze(H(:, :, 1, 1, 1)); for k = 2:size(H,3) hk = squeeze(H(:, :, k, 1, 1)); rho(k) = abs(h0(:)' * hk(:)) / (norm(h0(:)) * norm(hk(:))); end plot(0.5 * (0:length(rho)-1), rho); xlabel('阵元间距(波长)'); ylabel('相关系数');

容量代码里的n_tx和n_rx顺序按你zip包返回的实际维度调整，关键是H_eff的行列要和eye的单位矩阵维度对齐，否则行列式算出来没有物理意义。

6.3 用一组固定种子做交叉对照

最后一个技巧是给整个仿真流程固定随机种子，跑两组场景做增量对照：先跑2D模式记录容量基线，再开3D模式跑同一布局，观察容量和垂直扫描曲线的差异。差异符合预期，说明3D路径确实在起作用；差异为零，回查n_col。这个习惯帮我挡过不少“仿真结果好坏其实是模型设置差异”的争论。

拿到任何信道源码包，我都会先花十分钟做这三个验证再投算法。输出可信，后面波束设计、预编码方案的结论才站得住。希望帮到你。

本文还有配套的精品资源，点击获取