简介WINNER II 信道模型的 MATLAB 源代码包面向无线通信、5G/6G 及 MIMO 方向的研究者与工程师用于 3D-MIMO 场景下的信道建模与系统级仿真。压缩包共 52 个文件约 202KB以 44 个 .m 源文件为主涵盖信道建模、路径损耗、天线响应、信道参数与阵列配置等核心模块覆盖从传播路径计算、参数生成到天线响应的完整链路另有 4 个 .p 加密函数、1 个 .mlx 实时脚本及说明文档便于对照学习与二次开发。包内提供三个实例可系统演示 WINNER II 模型在视线与非视线传播、空间分集与波束赋形、信道估计、信号检测及性能评估等环节的实现流程涉及路径损耗、角度与延迟生成、LMMSE 估计、ZF/MMSE 检测以及 BER、SER、吞吐量等指标计算。适合需要搭建 3D-MIMO 仿真平台、深入理解信道特性或开展课题预研的研究生与工程师参考。已有 487 人次浏览学习可直接作为学术研究或工程项目的起点。1. 第一眼看到这份 WINNER II 源码包先明白 3D-MIMO 仿真缺的是什么如果你要做 5G/6G 信道仿真或者导师扔给你一句“用 WINNER II 搭个 3D-MIMO 链路”这份源码包基本是绕不开的起点。它不是点开就能出图的 demo而是一整套信道模型工程天线阵列配置、场景映射、LOS/NLOS 判定、H 矩阵生成全散落在几十个 .m 与 .p 文件里。适合两类人写论文需要真实信道、不想再拿复高斯噪声凑数的研究生做链路仿真、想给系统级平台换 3D 信道模型的工程师。先说结论这套代码能跑通但想改对必须从 wimparset、layoutparset、wim 这条主链路下手别一上来就改 3D_AntennaResponse.m。2. WIM2 主链路从 layoutparset 到 wim摸清 .p 文件的三个黑匣子拿到压缩包之后多数人的第一反应是直接在 MATLAB 里双击 example_channel_matrix.m回车。这个习惯不算错但我建议先花十分钟把文件清单过一遍因为这份代码的特点是文件多、依赖隐式、还有几个 .p 文件不能打开看。下面这张表是我按功能拆出来的分组方便你对照着查。分组文件作用主入口wim.m、wim_core.m、struct_generation.m、layout2link.m仿真调度与链路建立参数配置wimparset.m、layoutparset.m、arrayparset.m、ScenarioMapping.m模型、布点、阵列、场景天线方向图AntennaArray.m、AntennaResponse.m、3D_AntennaResponse.m、dipole.m、antenna_pol_vect.m天线响应与极化处理加密模块G_Calc1D_simple.p、G_Calc3D_simple.p、Aperture_Calc.p、interpbp.p增益计算、口径计算、宽带路径插值示例与文档example_channel_matrix.m、example_EADF_approx.m、example_syntetic_arrays.m、WIM2_ver064_notes.txt三个可运行实例与版本说明2.1 最小可运行链路先让 example_channel_matrix.m 出结果我一般不会先改任何参数而是先把 example_channel_matrix.m 跑通用 size() 把关键变量的维度全部打印出来确认当前 MATLAB 版本与这份代码兼容。以下是这个示例的核心调用逻辑。% 清理环境避免上一次运行留下的变量干扰 clear; close all; clc; % 布局参数小区半径、站点坐标、链路数等都由 layoutparset 提供默认值 layout layoutparset; % 模型参数载频、带宽、场景编号、极化方式等 [wim_par, link_par] wimparset; % 先跑单链路确认逻辑通 layout.num_links 1; layout.link_num 1; % 主入口data 返回信道数据simu_par 返回仿真过程参数 [data, simu_par, link_par] wim(layout, wim_par); % 打印信道矩阵维度一般能看到 nr x nt x n_paths 的结构 size(data.channel)这段代码的调用顺序是固定的layoutparset 返回布点与链路结构体小区半径、用户位置都在里面wimparset 返回信道模型参数包括路径损耗系数开关、阴影衰落开关等wim 执行真正的计算三个输出参数分别对应信道数据、仿真中间参数和链路参数。调试时我只动 layout.num_links其他保持默认这样出问题能排除“链路数过多导致计算量爆炸”的可能。跑完以后如果 size(data.channel) 能输出三维或四维矩阵说明主链路已经通了。附带说一下压缩包里还有一个 layoutparset.asv那是原作者编辑 layoutparset.m 时留下的自动保存备份和代码运行无关可以忽略。2.2 三个 .p 文件是黑匣子它们的输入顺序不能猜G_Calc1D_simple.p 和 G_Calc3D_simple.p 负责天线增益计算Aperture_Calc.p 负责口径计算interpbp.p 做宽带路径插值。这是 MATLAB 的 pcode 加密格式源码不可见调试时看到的报错往往很隐晦比如“输入参数的数目不足”或“输出参数太多”。最靠谱的做法不是猜参数顺序而是去读调用它们的 m 文件。# 在项目根目录搜一下哪些 .m 文件引用了 .p 文件 grep -n G_Calc3D_simple *.m grep -n interpbp *.m搜出的结果里能看到调用者的传参排列方式这就是 .p 文件实际期望的接口。WIM2_ver064_notes.txt 是版本说明对应 WIM2 ver 0.6.4跑不通的时候先打开它里面通常会记录版本边界和已知限制尤其是 MATLAB 版本与 p 文件编译格式的匹配问题。遇到 .p 报错时优先比对调用行的参数个数不要尝试把 .p 转回 .m这没有后悔药。提示在新版本 MATLAB 上如果 .p 文件加载失败先怀疑是压缩包里的 pcode 版本过旧换一台旧版本 MATLAB 验证最直接。2.3 struct_generation.m 和 layout2link.m 的作用struct_generation.m 负责把 wimparset 生成的参数结构体展开成每条链路独立的参数副本layout2link.m 负责把 layout 里的站点坐标映射到具体链路的收发端。这两步决定了你有没有“改一个参数影响所有链路”的问题。常见做法是改动参数后用 isequal 对比两条 link_par 的相关字段确认是否需要单独赋值。比如你在 wim_par 里改了场景编号两条链路都应该各自保留自己的传播场景副本而不是共享同一份指针。把这一层拖完后你会理解为什么很多人说 WIM2 难上手不是算法有多深而是文件依赖太多、黑匣子又多一上来就改参数很容易翻车。下一章把天线阵列和场景映射拆开讲这是 3D-MIMO 与传统 2D 信道模型最大的区别点。3. 3D 天线阵列与场景映射把 4×4 URA 配置和 LOS 概率接进信道模型3D-MIMO 的核心不只是“多加几根天线”而是三个维度同时起作用方位角、俯仰角、极化。WINNER II 默认阵列偏 2D 平面但做 3D 扩展时你要处理的就是这三者的联合。3D_AntennaResponse.m 是这套机制里的核心文件。3.1 arrayparset自定义均匀矩形阵列URA在 3D-MIMO 里最常见的阵列就是 URA行数控制垂直方向的波束个数列数控制水平方向的口径。下面是我常用的 4×4 双极化配置。% 自定义一个 4 行 x 4 列的均匀矩形阵列URA array arrayparset; % 载入默认阵列参数 array.num_ant 16; % 总阵元数 4 * 4 array.num_row 4; % 行数对应垂直维度 array.num_col 4; % 列数对应水平维度 array.element_spacing [0.5 0.5]; % 行、列方向阵元间距单位波长 array.pol VH; % 双极化垂直 水平阵元间距取 0.5 波长是常规选择小于 0.5 波长会增大阵元间相关性波束变宽大于 0.5 波长虽然能提高角度分辨率但会出现栅瓣方向图出现多余的峰值。垂直方向的行数越多俯仰角分辨率越好这是 3D-MIMO 相对 2D 方案在垂直维度上的额外自由度代价是天线响应矩阵的计算量成倍增长。arrayparset 输出的结构体最终会传给 3D_AntennaResponse.m由它生成每个阵元的 3D 方向图响应。3.2 ScenarioMapping.m 与 LOSprobability.m场景先于参数WINNER II 把传播环境划分成若干标准场景室内 A1、城区微蜂窝 B1、城区宏蜂窝 C1/C2、乡村 D1 等。ScenarioMapping.m 负责把场景名或场景编号映射到对应的系数表ScenParTables.m 里存的是这些系数。选错场景路径损耗和角度扩展都会跟着错所以第一步永远是确定场景。% 获取 B1 城区微蜂窝对应的场景编号 info ScenarioMapping(B1); % 生成距离-LOS 概率曲线距离单位是米 d 0:10:500; p_los LOSprobability(d, info); plot(d, p_los, LineWidth, 1.5); xlabel(距离 (m)); ylabel(LOS 概率); title(B1 场景 LOS 概率随距离变化);LOSprobability 的第二个参数可以直接传场景编号关键是和 ScenarioMapping 的输出保持一致。LOS 概率曲线决定了每条链路走 LOS 还是 NLOS 路径损耗公式B1 场景里 50 米内大概率是 LOS超过 200 米基本进入 NLOS。这个拐点直接改变路径损耗的斜率如果你发现仿真结果整体偏大或偏小先把布点距离和场景匹配关系查一遍。下表是我常用的场景选择参考。应用场景推荐 WINNER 场景典型站点间距室内热点A110–50 m城区微蜂窝B1100–500 m城区宏蜂窝C1 / C2500–1500 m乡村广覆盖D11–5 km3.3 角度生成与坐标转换fixedAods、fixedAoas 与 mycart2sphfixedAods.m 和 fixedAoas.m 用于固定出发角AoD与到达角AoA。做定向波束仿真或码本扫描时把目标角度填进这两个函数的参数里链路就不会随机调整角度方便你做波束对齐的对照实验。mycart2sph.m 与 mysph2cart.m 是坐标转换工具rotate_vector.m 负责把局部坐标系里的方向向量转到全局坐标系。% 把局部坐标系中的朝右向量绕 Z 轴旋转 30 度 v_local [1 0 0]; v_global rotate_vector(v_local, [0 0 1], 30*pi/180); % 将全局坐标向量转成球坐标输出方位角、俯仰角、半径 [az, el, r] mycart2sph(v_global(1), v_global(2), v_global(3));az 和 el 的单位是弧度拿到后要转成角度再和 3D_AntennaResponse.m 的方向图索引对齐。这类坐标转换是 3D 仿真里最容易出错的地方我踩过最典型的坑是俯仰角定义搞反有的模块以水平面为 0 度、向上为正有的以天顶为 0 度结果方向图主瓣直接朝着地面。所以每次新建工程我都会先用一个小脚本单独画一次方向图确认主瓣指向再接入整条链路。3.4 双极化天线响应与 antenna_pol_vect.mantenna_pol_vect.m 生成极化向量VH 双极化会生成两组正交通道信道矩阵的天线维度会翻倍。这个细节在容量计算时特别容易忽略你以为是 16×16 的系统实际上因为有双极化等效信道可能是 32×32。% 生成偶极子在指定方向的极化向量 pol_v antenna_pol_vect(theta, phi, V); pol_h antenna_pol_vect(theta, phi, H);每个极化向量包含三个分量对应 xyz 方向的场分量。双极化天线在 3D-MIMO 里几乎是标配因为极化分集能在不增加天线阵元数量的情况下提高链路稳定性但代价是信道矩阵维度变大、计算量增加提取 H 矩阵后做奇异值分解时你会发现奇异值分布被拉得更均匀这是极化分集的典型特征。4. WIM2 常见问题与避坑六个最容易让仿真中断的源头先说结论这个模型跑不通八成问题出在参数顺序和结构体字段而不是算法本身。下面六条是我在复现过程中遇到过的典型问题按现象、原因、解决三段写方便你对照排查。4.1 .p 文件报错“输入参数数目不足”现象wim 执行到一半控制台输出类似 “Error using G_Calc3D_simple / Too many input arguments”链路直接中断。 原因.p 文件是 pcode看不到函数头不同版本的 WIM2 里接口有差异调用时参数个数或顺序和原版不一致。 解决用 grep 在当前目录搜出所有引用该 .p 文件的 m 文件把调用行的参数列表原样复制到你的脚本不要凭记忆补参数。这一条是我见过最多的翻车点没有之一。4.2 路径损耗出现 NaN 或巨大负数现象pathloss.m 输出的路径损耗矩阵里出现 NaN或出现 -300 dB 之类的离谱值。 原因常见两种。一是站点和终端距离为 0对数距离公式里的 log10(0) 是负无穷二是传参时距离单位用错把米当成了千米。 解决给距离加下限站点间距小于 1 米时直接置为 1。同时检查 wim_par 里载频单位是否为 GHz别用 Hz 代入否则公式里的频率项会完全错位。4.3 3D 方向图全零或主瓣方向反了现象3D_AntennaResponse.m 返回的全 0 矩阵或者 example_syntetic_arrays.m 里画出的方向图主瓣指向地面。 原因俯仰角定义不统一。有的模块以水平面为 0 度、向上为正有的以天顶为 0 度从 mycart2sph 拿到弧度后没转角度就对齐方向图网格。 解决跑完整链路前先单独画一次方向图确认最大值落在期望的方位角和俯仰角再接进 wim。这个验证只需要一分钟但能省下半天排查时间。4.4 struct_generation 生成的字段和 wimparset 对不上现象用 wimparset 生成参数后直接给 layout 的某些字段赋值报字段不存在。 原因版本混用。网上下载的这份代码可能和你参照的 WIM2 基础版本不一致wimparset 输出的结构体字段有增删。 解决用 isfield 判断关键字段是否存在缺失的字段用点号直接赋值补上。以 wimparset.m 输出的实际结构体为准不要照抄网上旧帖子的字段名。4.5 链路数一多仿真时间成倍涨现象num_links 从 1 改成 8 之后运行时间从几十秒涨到十几分钟。 原因每条链路都独立完成角度生成、散射簇构建、路径损耗和天线响应计算多链路之间没有共享中间结果计算量接近线性甚至超线性增长。 解决调试阶段始终保持 num_links1正式批量前先用 num_links2 做一次时间基准测试。如果只是做容量对比可以分段跑每段 100 次用 parfor 做并行别一次性全塞进内存。4.6 MATLAB 新版本加载旧 p 文件报格式错误现象运行到 .p 文件时报 “Unable to read P-file” 或文件格式不识别。 原因p 文件的编译格式与 MATLAB 大版本绑定老版本编译出来的 p 文件在新版本上可能加载不了。 解决在项目文件夹里放一个版本说明标注这个包在哪一版 MATLAB 上验证过。如果新版本跑不了找一台旧版本 MATLAB或者换用同版本下重新编译的 p 文件。这个只能靠环境匹配没有代码层面的解决办法。提示压缩包里的 RP-121412.zip 是配套参考压缩包先跑通主链路再解压很多问题不是缺文件而是上面这些环境与参数问题。5. 把示例改造成自己的链路自定义阵列、H 矩阵提取与容量/BER 接法前面几章把坑扫了一遍这章开始真正干活。目标是把 example_channel_matrix.m 里写死的默认配置替换成你自己的天线阵列和场景并拿到能用于容量或误码率仿真的 H 矩阵。5.1 替换默认阵列从 example_syntetic_arrays.m 借用法改阵列之前先打开 example_syntetic_arrays.m 看一遍它专门演示天线阵列的结构配置。然后按下述方式做替换。% 在 wim 调用之前用自己的 array 覆盖默认阵列 array arrayparset; array.num_ant 16; array.num_row 4; % 垂直方向 4 层 array.num_col 4; % 水平方向 4 列 array.element_spacing [0.5 0.5]; % 把阵列信息写进 layoutWIM2 在主链路里按 layout 的阵列字段调用天线响应 layout.array array; [data, simu_par, link_par] wim(layout, wim_par);layout.array 这种赋值方式是我在多个版本里惯用的接法不同版本字段名可能有差异稳妥起见在源码里搜一下 array 关键字确认实际用的字段名。num_row 和 num_col 的乘积必须等于 num_ant不相等时 3D_AntennaResponse.m 会报维度不一致。改完阵列后建议先画一下方向图确认主瓣指向再跑主链路。5.2 从 data 里提取 H 矩阵并做归一化信道矩阵从 data 结构里取出来后第一个动作永远是打印维度确认顺序。% 打印 data.channel 的尺寸确认维度顺序 [nr, nt, n_paths] size(data.channel); % 取第一个时延抽头的所有天线得到链路冲激响应矩阵 H data.channel(:, :, 1); % 归一化让总的信道能量不随天线数变化方便对比容量 Hn H / norm(H, fro) * sqrt(nr * nt);squeeze 在这里没有用到因为直接索引就把单独维度压掉了。归一化的目的是让不同天线配置下比较容量时有一个公平基准如果不做归一化64 天线系统的容量天然高于 16 天线体现的是能量差异而不是复用增益差异。论坛里经常看到有人容量结果忽大忽小根因往往就是这一步漏了。5.3 容量计算SVD 与等功率注水拿到归一化 H 矩阵后容量计算就是标准的奇异值分解问题。% 容量计算等功率分配接收端已知 CSI snr_dB 20; snr 10^(snr_dB / 10); [~, S, ~] svd(Hn); Nt size(Hn, 2); % 等功率分配每根发射天线的功率是总功率的 1/Nt capacity sum(log2(1 snr / Nt * S.^2));svd 返回的奇异值 S 是列向量奇异值平方代表对应空间信道的增益。如果 Hn 秩不足奇异值里有 0对应项 log2(10) 为 0不贡献容量这正好反映了信道在空间维度上的“可用自由度”。把这段代码套进 SNR 循环就能画出 3D-MIMO 的容量曲线。5.4 BER 接入ZF 与 MMSE 的最小对比容量是理论上界工程上还要看检测算法能不能逼近它。以下是 ZF 和 MMSE 的最小接入代码。% 16QAM 调制ZF 与 MMSE 检测对比 M 16; s randi([0 M-1], Nt, 1); % 随机符号索引 x qammod(s, M, UnitAveragePower, true); noise (randn(nr, 1) 1j*randn(nr, 1)) / sqrt(2) * sqrt(1/snr); y Hn * x noise; % ZF 检测直接求逆 x_zf Hn \ y; s_zf qamdemod(x_zf, M, UnitAveragePower, true); % MMSE 检测加正则项防止噪声放大 x_mmse (Hn*Hn (1/snr)*eye(Nt)) \ (Hn * y); s_mmse qamdemod(x_mmse, M, UnitAveragePower, true); % 统计错误符号数 err_zf sum(s_zf ~ s); err_mmse sum(s_mmse ~ s);1/snr 就是噪声方差加在 MMSE 的正则项里。SNR 很高时 MMSE 趋近 ZF两者错误数会接近SNR 低时 MMSE 明显优于 ZF因为正则项抑制了噪声放大效应。这段代码可以直接塞进 Monte Carlo 循环统计 BER 曲线验证你改的阵列配置在高阶调制下是否还能正常工作。5.5 三个 example 的分工和阅读顺序example_channel_matrix.m 输出完整的信道矩阵适合直接接后续容量或检测模块example_syntetic_arrays.m 聚焦天线阵列结构适合验证方向图和阵列配置example_EADF_approx.m 做的是有效孔径分布函数近似研究天线口径时才用得上。我建议的阅读顺序是先用 example_channel_matrix.m 跑通主链路再借 example_syntetic_arrays.m 确认阵列最后才决定要不要碰 EADF三份示例分别对应“信道”“阵列”“口径”三个层次按需取用。6. 批量仿真的收尾技巧参数扫描与随机布点的落地顺序到这一步wim 的调用基本不再是玄学。我的习惯是不管是论文里的容量曲线还是工程里的覆盖统计都不贪多先把单链路调通再往上堆样本。6.1 参数扫描SNR 循环与信道多次平均snr_list 0:5:30; capacity_avg zeros(size(snr_list)); N_trial 100; for k 1:numel(snr_list) cap_sum 0; for trial 1:N_trial % 每次调用 wim 都会重新生成一次信道 [data, ~, ~] wim(layout, wim_par); H data.channel(:, :, 1); Hn H / norm(H, fro) * sqrt(size(H,1) * size(H,2)); S svd(Hn); snr 10^(snr_list(k) / 10); cap_sum cap_sum sum(log2(1 snr / size(H,2) * S.^2)); end capacity_avg(k) cap_sum / N_trial; end plot(snr_list, capacity_avg); xlabel(SNR (dB)); ylabel(平均容量 (bit/s/Hz));内层 trial 循环的用途是让随机性平均化WIM2 每次调用会重新生成角度和散射簇不平均的话容量曲线会抖动得很厉害。外层循环扫描 SNR 时不需要改 wim_par 里的噪声参数只改 snr 变量本身避免引入额外变量。6.2 随机布点用 StationDistXY 生成站点间距离固定坐标只能做单点验证做统计仿真时我会用圆形均匀布点代替手动坐标输入。% 基站固定原点10 个用户在半径 200m 的圆内均匀撒点 num_users 10; R 200 * sqrt(rand(num_users, 1)); % 均匀分布在圆内 theta 2 * pi * rand(num_users, 1); x R .* cos(theta); y R .* sin(theta); plot(x, y, o); title(随机布点); % 把坐标交给 StationDistXY 计算站点间距离矩阵 layout.station_dist StationDistXY(x, y);半径用 sqrt(rand) 而不是 rand否则点会往圆心集中布点密度不均匀。如果项目要求泊松点过程就先按每个子区域面积乘以密度算出期望点数再用 poissrnd 抽取实际点数在子区域内均匀撒点效果上比固定网格更贴近真实基站部署。6.3 最后一道验证先和自由空间损耗对一下再收工批量跑之前我会用一个固定距离把仿真路径损耗和自由空间损耗做一次对照。% 100 米2.6 GHz 下的自由空间损耗参考值 d_km 0.1; fc_mhz 2600; pl_fspl 20*log10(d_km) 20*log10(fc_mhz) 32.44; % 计算结果约为 80.7 dBWIM2 的 B1 LOS 路径损耗在 100 米、2.6 GHz 时会比自由空间损耗高出几 dB因为引入了额外的穿透和大气衰减项。如果仿真值和这个参考值差出 30 dB 以上基本可以断定场景选成了 NLOS或者 ScenParTables.m 的参数表加载出了问题。这个验证习惯我保持到现在不管换哪个版本先对一次路径损耗再对一次 LOS 概率曲线两项都通过才上批量仿真。从那次 .p 文件参数顺序搞反之后我每次拿到新代码包都是先跑最小链路、再验证基准值最后才动配置参数。希望帮到你。本文还有配套的精品资源点击获取