WINNER II信道模型MATLAB实战：从链路预算翻车到MIMO波束赋形验证

简介：本资源面向无线通信仿真方向的研究生、工程师及科研人员，提供WINNER II信道模型的MATLAB实现代码，可视为3GPP SCM信道模型的增强版本，适用于MIMO系统建模、链路级仿真与信道特性分析等场景。压缩包共47个文件，约639KB，以41个.m脚本为核心，涵盖信道核心计算、天线阵列与极化响应、路径损耗、角度扩展、场景映射及参数配置等模块，另含4个.p文件、1份PDF说明文档与1个txt版本记录，便于理解代码结构与调用流程。目前已有1655人学习下载，说明其在相关领域具备一定参考价值。读者可借助配套文档与示例脚本，快速掌握WINNER II信道的参数设置、场景配置与信道矩阵生成方法，并在此基础上开展二次开发与算法验证。

1. 从一次链路预算翻车说起：WINNER II 信道模型到底能干什么

去年帮一个做室内覆盖的团队看仿真报告，他们用自由空间路损加瑞利衰落跑出来的覆盖图，和现场实测差了将近 12 dB。问题不在发射功率，也不在天线增益，而是信道模型选错了——自由空间加瑞利只适合窄带、低相关场景，而他们测的是 2.6 GHz 宽带 MIMO，多径时延扩展和空间相关性完全没建模。换成 WINNER II 之后，误差压到了 3 dB 以内。这就是 WINNER II 信道模型存在的意义：它不是教科书里那个理想化的瑞利或莱斯，而是一套面向宽带 MIMO 系统的几何随机信道模型，能同时刻画时延扩展、角度扩展、多普勒谱和极化，覆盖从室内热点到城市宏站再到农村宏站的完整场景谱系。

这套 MATLAB 代码包做的事情很具体：给定场景类型、载波频率、带宽、天线阵列配置和移动速度，它输出信道系数矩阵，你可以直接拿去做链路级仿真、波束赋形验证、MIMO 容量计算或者接收机算法测试。适合谁用？做 5G/6G 物理层仿真的研究生、做基站或终端天线评估的射频工程师、以及需要给系统级仿真喂真实信道数据的链路预算人员。如果你只用过rayleighchan或者comm.RicianChannel，那这套代码会让你第一次感受到“空间相关性”这四个字在 MIMO 里到底意味着什么。

2. WINNER II 的建模逻辑：为什么它不是“又一个瑞利信道”

2.1 几何随机模型的核心：把多径当成“簇”来处理

WINNER II 和传统统计信道模型最大的区别在于，它不直接假设包络服从某个分布，而是从几何出发：先随机撒一批散射体，再计算每条路径的时延、角度和复增益，最后叠加成信道冲激响应。具体来说，它把多径分成若干“簇”，每个簇内部又有多条子径。簇的时延服从某种分布，簇内子径的角度偏移服从另一组分布，所有这些参数都由场景决定。

这种建模方式带来的直接好处是：时延扩展和角度扩展是“自然涌现”的，而不是你手动指定的。比如在城市宏站场景下，基站高、终端低，散射体主要集中在终端周围，所以到达角扩展大、离开角扩展小。你不需要告诉模型“离开角扩展是 5 度”，它自己会算出来。这也是为什么 WINNER II 在 MIMO 波束赋形仿真里比rayleighchan靠谱得多——后者根本没有角度维度。

常见做法是：先确定场景（A1 室内热点、B1 城市微站、C1 城市宏站、D1 农村宏站等），然后查表拿到该场景下时延扩展、角度扩展、阴影衰落的标准差等大尺度参数，再用这些参数驱动几何随机过程。代码包里通常有一个winner2_core或类似名字的函数，输入一个配置结构体，输出信道系数。

2.2 场景参数表：选错场景比调错参数更致命

很多人拿到代码直接跑默认配置，结果发现仿真结果和实测对不上，回头查代码才发现默认场景是 A1 室内热点，而自己做的明明是 C1 城市宏站。WINNER II 的场景分类不是随便分的，每个场景对应一组固定的大尺度参数，这些参数来自大量实测数据的统计拟合。下面这张表是几个常用场景的关键参数对比，数值来自 WINNER II 最终报告的标准表格：

这张表里最容易被忽略的是阴影衰落标准差。C1 和 D1 都是 8 dB，但 D1 的时延扩展是 C1 的近两倍。如果你做的是农村覆盖，用 C1 的参数会导致时延扩展偏小，均衡器设计过于乐观，实测误码率会打脸。我一般会建议：场景选完之后，把时延扩展和角度扩展的分布类型也确认一遍，WINNER II 里不同场景用的分布不一样，有的是对数正态，有的是均匀，代码里通常用randn或rand配合变换实现。

2.3 从配置结构体到信道系数：一次完整的调用链

代码包的核心接口通常是一个配置结构体，字段包括场景名、载波频率、带宽、发射天线数、接收天线数、天线间距、移动速度、仿真时长等。下面是一个典型的调用示例，我把它拆成三段来看：

% 第一段：构建配置结构体 cfg = struct(); cfg.scenario = 'C1'; % 场景：城市宏站 cfg.fc = 2.6e9; % 载波频率 2.6 GHz cfg.bw = 20e6; % 带宽 20 MHz cfg.Nt = 4; % 发射天线数 cfg.Nr = 4; % 接收天线数 cfg.dt = 0.5; % 发射天线间距（波长倍数） cfg.dr = 0.5; % 接收天线间距 cfg.v = 30; % 移动速度 km/h cfg.T = 0.01; % 仿真时长 秒 cfg.fs = 30.72e6; % 采样率

这段配置里，dt和dr用波长倍数表示，0.5 表示半波长间距，这是 MIMO 里最常见的配置。v是移动速度，直接影响多普勒频移。fs是采样率，通常要大于等于带宽，这里 30.72 MHz 是 LTE 的典型采样率。配置结构体建好之后，调用核心函数：

% 第二段：生成信道系数 [H, delay, angle] = winner2_core(cfg); % H 维度：Nr x Nt x Npath x Ntime % delay：每条路径的时延 % angle：每条路径的离开角和到达角

H是四维矩阵，四个维度分别是接收天线、发射天线、路径数和时间采样点。这个维度顺序很关键，后面做波束赋形的时候如果搞反了，得到的波束图会完全不对。delay和angle可以用来验证模型是否合理，比如 C1 场景下最大时延不应该超过几微秒，到达角扩展应该在几十度量级。

% 第三段：验证信道统计特性 H_flat = squeeze(H(:,:,1,:)); % 取第一条路径 R = H_flat * H_flat' / size(H_flat,3); % 接收相关矩阵 eig_val = eig(R); disp(['条件数：', num2str(cond(R))]);

这段验证代码做的是：取第一条路径的信道系数，计算接收相关矩阵，然后看条件数。条件数过大说明空间相关性太强，可能是天线间距设得太小或者场景参数不对。我一般会要求条件数在 10 以内，超过 20 就要回头检查配置。

3. 把代码跑起来：从环境准备到第一个信道系数

3.1 环境依赖与文件结构：别急着改代码，先看依赖

这套代码包通常包含几个核心文件：主函数、场景参数表、天线阵列生成函数、以及一些辅助工具。在 MATLAB 里跑之前，先确认你的版本支持randn的种子控制，因为 WINNER II 的几何随机过程依赖随机数，不同 MATLAB 版本的随机数生成器默认算法可能不同，会导致结果不可复现。我一般会在脚本开头加一行rng(42)，固定种子。

文件结构上，常见做法是把所有.m文件放在同一个目录下，然后在脚本里addpath进去。不要试图把函数文件改成脚本文件，WINNER II 的代码里有很多嵌套函数和持久变量，改结构容易出玄学问题。另外，如果你用的是 MATLAB R2016b 之前的版本，注意string类型和struct的字段名大小写敏感问题，场景名'C1'和'c1'在代码里可能被当成两个不同的场景。

提示：跑之前先用which winner2_core确认函数在路径里，如果返回空，说明addpath没生效。

3.2 最小可运行示例：四发四收城市宏站

下面这个脚本是我自己常用的最小验证模板，直接复制就能跑。它做了一件完整的事：配置 C1 场景、生成信道、计算容量、画图。

% winner2_minimal.m clear; clc; rng(42); % 配置 cfg.scenario = 'C1'; cfg.fc = 2.6e9; cfg.bw = 20e6; cfg.Nt = 4; cfg.Nr = 4; cfg.dt = 0.5; cfg.dr = 0.5; cfg.v = 30; cfg.T = 0.01; cfg.fs = 30.72e6; % 生成信道 [H, delay, angle] = winner2_core(cfg); % 计算遍历容量 SNR_dB = 0:2:20; C = zeros(size(SNR_dB)); for k = 1:length(SNR_dB) snr = 10^(SNR_dB(k)/10); for t = 1:size(H,4) Ht = H(:,:,:,t); Hsum = sum(Ht, 3); % 合并所有路径 C(k) = C(k) + log2(det(eye(cfg.Nr) + ... snr/cfg.Nt * Hsum * Hsum')); end C(k) = C(k) / size(H,4); end % 画图 figure; plot(SNR_dB, C, 'b-o', 'LineWidth', 1.5); xlabel('SNR (dB)'); ylabel('容量 (bps/Hz)'); title('C1 场景 4x4 MIMO 遍历容量'); grid on;

这段代码的逻辑说明：H的第四维是时间采样，每个采样点对应一个信道实现。sum(Ht, 3)是把所有路径的贡献合并，得到等效的窄带信道矩阵。容量公式用的是标准的 MIMO 容量表达式，det里的eye(cfg.Nr)是接收端噪声协方差归一化后的结果。参数SNR_dB从 0 到 20，步长 2，覆盖了典型的蜂窝通信工作区间。

跑完之后你会看到一条单调上升的曲线，20 dB 时容量大概在 15-20 bps/Hz 之间。如果曲线明显偏低，检查cfg.Nt和cfg.Nr是否和H的维度匹配；如果曲线抖动特别大，说明时间采样点太少，把cfg.T加大或者cfg.fs提高。

3.3 参数怎么改：频率、带宽、天线数的联动关系

改参数不是随便填数字，有几个联动关系必须遵守。第一，cfg.fs必须大于等于cfg.bw，否则时域采样会欠采样，信道系数会出现混叠。第二，cfg.dt和cfg.dr用波长倍数表示，如果你知道实际天线间距（比如 0.1 米），要先换算：dt = d / (c / fc)，其中c是光速。第三，cfg.v影响多普勒频移，最大多普勒频移fd = v * fc / c，如果fd接近fs/2，会出现多普勒混叠，这时候要么降低v，要么提高fs。

我一般会先算一遍fd，确认它远小于fs/2再跑仿真。比如v=30 km/h、fc=2.6 GHz，fd大约是 72 Hz，而fs=30.72 MHz，完全没问题。但如果做高铁场景，v=300 km/h，fd就到 720 Hz 了，虽然还是远小于fs/2，但多普勒谱的形状会变得很重要，WINNER II 里有多普勒谱的配置选项，默认是经典谱，高铁场景可能需要改成平坦谱。

4. 避坑与排查：那些让仿真结果对不上实测的细节

4.1 现象：容量曲线在高 SNR 段饱和，原因：天线相关性没配对

有人跑完 4x4 容量发现 20 dB 时只有 10 bps/Hz，远低于理论值。查了半天代码没报错，最后发现是cfg.dt和cfg.dr都设成了 0.1，天线间距太小，空间相关性接近 1，信道矩阵接近秩 1，容量自然上不去。解决方法是把间距改到 0.5 波长以上，或者显式检查相关矩阵的条件数。我一般会在生成信道后加一行cond(Hsum * Hsum')，条件数超过 100 就说明相关性太强。

4.2 现象：时延扩展明显偏小，原因：场景参数表被覆盖

WINNER II 的场景参数表通常是一个独立的.mat文件或函数，里面存着每个场景的时延扩展、角度扩展等。如果你在脚本里不小心定义了一个同名变量，或者修改了参数表但没保存，模型会用一个默认值。更隐蔽的情况是：参数表里的时延扩展是“均值”，而模型实际用的是对数正态分布，如果你只看均值不看分布，可能会觉得时延扩展“应该更大”。解决方法是打印delay的最大值和标准差，和场景参数表对照。

4.3 现象：不同 MATLAB 版本结果不一致，原因：随机数生成器变了

MATLAB R2017a 之后默认随机数生成器从twister变成了simdTwister，虽然都是梅森旋转，但种子相同的情况下产生的序列不同。如果你的代码依赖randn生成几何随机散射体，换版本后结果会变。解决方法是显式指定生成器：rng(42, 'twister')。这个坑我在两个不同版本的 MATLAB 上复现过，血泪经验是：任何依赖随机数的仿真，脚本开头必须固定生成器和种子。

4.4 现象：多普勒谱形状不对，原因：移动方向没设对

WINNER II 默认假设移动方向是随机的，但如果你做的是定向覆盖场景（比如沿公路的基站），移动方向应该固定。代码里通常有一个cfg.mobility_direction或类似字段，默认是'random'。改成'fixed'并指定角度后，多普勒谱会从经典 U 型变成偏移的 U 型。这个参数在高速场景下影响很大，低速场景下不明显，容易被忽略。

4.5 现象：跑大规模 MIMO 时内存爆了，原因：H 矩阵维度没控制

H是四维矩阵，维度是Nr x Nt x Npath x Ntime。如果你做 64x64 大规模 MIMO，Npath有几十条，Ntime有几千个采样点，内存占用会到 GB 级别。解决方法是分块生成：每次只生成一个时间窗的信道，处理完再生成下一个。代码包里通常有winner2_core_chunk或类似的分块接口，如果没有，可以自己写循环，每次改cfg.T和起始时间。

5. 进阶用法：用信道系数验证波束赋形算法

5.1 从 H 矩阵提取空间协方差矩阵

波束赋形的核心是空间协方差矩阵。对于 WINNER II 生成的信道，你可以按路径分别计算协方差，也可以合并所有路径。我一般会先看第一条路径的协方差，因为第一条路径通常是直射或最强径，对波束方向影响最大。下面这段代码演示了如何提取并可视化：

% 提取第一条路径的接收协方差 H1 = squeeze(H(:,:,1,:)); % Nr x Nt x Ntime Rr = zeros(cfg.Nr, cfg.Nr); for t = 1:size(H1,3) h = H1(:,:,t); Rr = Rr + h * h'; end Rr = Rr / size(H1,3); % 特征分解 [V, D] = eig(Rr); [~, idx] = sort(diag(D), 'descend'); V = V(:, idx); % 画波束图 theta = -90:0.5:90; w = V(:,1); % 主特征向量 pattern = zeros(size(theta)); for k = 1:length(theta) a = exp(1j * pi * (0:cfg.Nr-1)' * sind(theta(k))); pattern(k) = abs(w' * a)^2; end figure; plot(theta, 10*log10(pattern/max(pattern)), 'LineWidth', 1.5); xlabel('角度 (度)'); ylabel('归一化增益 (dB)'); title('基于 WINNER II 信道的主特征波束'); grid on;

这段代码的逻辑：Rr是接收协方差矩阵，eig分解后取最大特征值对应的特征向量作为波束权重。a是均匀线阵的导向矢量，sind(theta)里的theta是角度。画出来的波束图应该有一个明显的主瓣，指向散射体集中的方向。如果主瓣不明显或者出现多个等高的瓣，说明信道相关性太低，可能是天线间距太大或者散射体太分散。

5.2 用多普勒谱验证移动性建模

WINNER II 的多普勒谱可以通过对时间维做 FFT 得到。下面这段代码计算并画出多普勒谱，用来验证移动速度设置是否正确：

% 多普勒谱分析 h = squeeze(H(1,1,1,:)); % 取第一根天线的第一条路径 L = length(h); Hf = fftshift(fft(h .* hanning(L), 1024)); f = linspace(-cfg.fs/2, cfg.fs/2, 1024); figure; plot(f, 20*log10(abs(Hf)/max(abs(Hf))), 'LineWidth', 1.2); xlabel('频率 (Hz)'); ylabel('归一化功率 (dB)'); title('多普勒谱'); xlim([-200 200]); grid on;

跑完之后你会看到一个 U 型谱，两个峰分别在正负最大多普勒频移处。如果谱是平的，说明多普勒建模没生效，检查cfg.v是否为零或者移动方向设置有问题。如果峰的位置和fd = v * fc / c对不上，检查fc和v的单位，v是 km/h 还是 m/s，代码里通常有换算，但不同版本的换算系数可能不一样。

5.3 一个我常犯的错误：忘了归一化

WINNER II 生成的信道系数默认是归一化的，总功率为 1。但如果你手动合并路径或者做波束赋形，很容易破坏归一化，导致 SNR 定义混乱。我一般会在生成信道后加一行H = H / sqrt(mean(abs(H(:)).^2))，强制归一化。这个习惯是从一次翻车经历来的：当时做容量对比，两个场景的容量差了 3 dB，查了一天才发现是其中一个场景的信道功率没归一化。从那以后我每次生成信道都强制走一遍归一化检查，mean(abs(H(:)).^2)必须在 1 附近，偏差超过 0.1 就要查原因。

希望帮到你。

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