5G NR小区搜索中的PSS-SSS联合检测与polar译码仿真

简介：基于polar编解码的PSS-SSS联合检测仿真资源，面向通信方向本硕博教研与学习人群，用于掌握polar编解码算法及PSS-SSS联合检测流程。资源基于MATLAB实现，输出联合检测成功率、误码率BER与误块率BLER等关键指标，可作为算法验证与课程设计参考。压缩包共70个文件，包含56个m脚本、5个mat数据、4个bin信道文件、2个txt说明、1个mexw64可执行文件及1个avi操作录像，整体大小2.91MB，目录按发端、收端、信道等模块划分，结构清晰。已吸引288人学习下载。配套操作视频演示了从运行Runme.m到查看仿真结果的完整过程，同时提供Polar编解码、PSS/SSS检测、信道估计等子函数，方便读者分模块研读与二次开发。

1. 开机找小区这件事，PSS-SSS 联合检测为什么要拉着 polar 编解码一起算

5G UE 上电后第一件事是小区搜索：先在时频域里把 PSS 和 SSS 找出来，确定符号定时、频偏和物理小区 ID，才能去解 PBCH；而 PBCH 的 MIB 是 polar 码保护的，BER 和 BLER 是衡量这条链路最直接的体检指标。很多做物理层仿真的同事习惯把同步检测和信道译码分成两套代码维护，可真实链路上同步估计的残差会直接转化为 LLR 劣化，最后 BLER 曲线对不上谁都不认账。这篇笔记围绕“基于 polar 编解码的联合检测”拆讲：从 SSB 生成、PSS-SSS 联合检测，到 PBCH 解调和 CA-SCL 译码，把检测成功率、BER、BLER 三个指标的统计口径和常见翻车点一次说清。适合正在做 5G NR 小区搜索算法验证、物理层链路仿真，或者想把同步与译码串成完整链路的从业者。

2. 先立住协议框架：PSS-SSS 给 polar 译码喂了什么输入

2.1 同步信号不是“找峰值”那么简单

5G NR 的 SSB（SS/PBCH block）时域占 4 个 OFDM 符号，频域占 240 个子载波。PSS 放在第 1 个符号的中间 127 个子载波上，是长度 127 的 m 序列，只有 3 个候选循环移位，对应小区内 IDN_ID^2 ∈ {0,1,2}；SSS 放在第 3 个符号的同样位置，由两个 m 序列相乘得到，携带小区组 IDN_ID^1 ∈ {0,1,…,335}。物理小区 ID 按PCI = 3 × N_ID^1 + N_ID^2计算，所以同步检测不只是找相关峰、对齐符号定时，它本质上是一个 3×336 的联合假设检验：既要猜对定时，也要猜对 PCI。

接收端的常规做法是本地生成所有候选序列，把接收序列和每个候选做互相关，最大峰值对应的候选同时给出 N_ID 和定时。这里要区分两类相关：时域滑动相关适合在还没有 OFDM 符号边界时做定时捕获；频域相关则适合在定时已知后做精估计。PSS 检测通常先做时域相关，因为接收机此时连 FFT 窗口在哪都不知道。实现时用一次xcorr算完整相关函数，比逐点滑动循环快一个量级。

联合检测的第一步就在 PSS 和 SSS 之间发生：用 PSS 估计出的定时位置去约束 SSS 的频域抽取窗口，只有 PSS 定时可信的候选才进入 SSS 假设检验。单独看 PSS，检测成功率做到 99% 以上不难，难的是加上 SSS 联合校验后还能维持住。我习惯把成功率口径定成“PCI 完全正确且定时误差落在半个 CP 以内”，而不是“相关峰超过门限”——口径一松，仿真报告里的数字好看，但进了真实信道全是问题。

2.2 NR polar 编码链：56 比特怎么变成 864 个编码比特

PBCH 是 polar 码在 5G NR 里的标准样本。信息比特 A=32 比特（包括 MIB 主体以及 SFN 高位、半帧指示、SSB index 等），先按 TS 38.212 的 CRC24C 附加校验，得到 B=56 比特；送入 polar 编码器，母码长度 N=512。polar 编码的核心是信道极化：对一个 2 的幂长度的信道向量做 Kronecker 变换，一部分子信道容量趋于 1，另一部分趋于 0。编码器把信息比特放在可靠度最高的 K 个位置，其余位置放冻结比特。NR 的可靠信道排序是按规范里的 Q 函数近似量化表计算的，和经典论文里巴氏参数的排序有细微差别，自研编码器时这里最容易埋雷。

编码之后不是直接调制，还要做速率匹配：子块交织加比特选择，把 512 比特变成 E=864 比特，再 QPSK 调制到 432 个 RE 上。E 比 N 长，所以这里的速率匹配以循环重复为主，但重复不是简单复制，而是交织后按顺序选取。为什么 PBCH 要重复这么多？因为一个 SSB 里 PBCH 需要占满 432 个 QPSK 符号，而信息比特只有 32 个，polar 母码 512 也不够填。这个“交织+重复”的组合让同步估计误差非常致命：同一个信息比特的多个副本被拉开到不同的频域位置，信道频率选择性会同时影响所有副本，LLR 一错全错。

译码端对应 SC 或 CA-SCL。PBCH 自带 CRC24C，天然适合 CRC 辅助的列表译码。list 大小 L=8 是常见折中，L=16 能再捞回零点几 dB，但译码复杂度和运行时间接近翻倍。我一般先把 L=8 跑通，最后再确认一次 L=16 的增益值不值得。译码器工作前必须先做与发射端镜像的去速率匹配：把 864 个 LLR 按交织映射关系放回 512 个母码位置，重复位做 LLR 相加，冻结比特位置置 0，然后才进 polar 译码器。整个尺寸链我建议打印出来核对：A=32 → B=56 → N=512 → E=864 → 432 symbols，哪个环节对不上，后面全白跑。

2.3 联合检测到底“联合”了什么

联合检测的信息流是：PSS 给出粗定时和粗频偏 → SSS 在约束窗口内确认 PCI 并做精频偏 → DMRS 做信道估计 → PBCH 软解调 → polar 译码。同步估计的残余误差不会停留在同步模块内部，它会一路传进 LLR。比如残余频偏 Δf 会让相邻 OFDM 符号间产生相位旋转，DMRS 估计出的信道相位如果跟不上这种旋转，数据 RE 的 LLR 幅度会被整体压缩，CA-SCL 译码器把这些不可靠比特当可靠比特处理，BLER 从 10^-2 恶化到 10^-1 很常见。

所以联合检测的“联合”不是指某个模块同时吐出同步和译码结果，而是指同步质量参数必须能一路传递到底，并且在仿真里可以单独注入偏差来观察影响。我在工程里把残余频偏、定时偏移做成可配置输入，解调模块统一从结构体里取参数。这样既能跑“理想同步”的基线，也能跑“真实同步”的劣化曲线。仿真时三个指标分开度量：检测成功率找同步级的问题，BLER 找端到端的问题，BER 夹在中间用来分辨是解调坏了还是译码坏了。三条曲线一起看，才能定位到具体是哪一级掉了链子。

3. 仿真链路搭建：从 SSB 生成到 PSS-SSS 联合检测

3.1 生成 SSB 资源网格：频域 240 点、时域 4 符号

先固定仿真参数：子载波间隔 30 kHz，常规 CP，FFT 大小 4096（对应 100 MHz 带宽），SSB 的频域起点给一个固定偏移，不做栅格盲检。资源网格里符号 0 放 PSS，符号 1 和 3 放 PBCH 数据加 DMRS，符号 2 中间放 SSS、两侧放 PBCH。PBCH 数据加 DMRS 一共 576 个 RE，其中 DMRS 按每 4 个 RE 一个的密度插入，数据 RE 正好 432 个，对应 864 个编码比特。

% 生成 PSS 和 SSS 序列 ncellid = 42; % 随便挑一个 PCI nid1 = floor(ncellid / 3); % 小区组 ID，0~335 nid2 = mod(ncellid, 3); % 小区内 ID，0~2 pss = nrPSS(nid2); % 127 长 PSS m 序列 sss = nrSSS(nid1); % 127 长 SSS 序列 % 摆进 SSB 资源网格：240 子载波 x 4 符号 ssb_grid = zeros(240, 4); pss_idx = (240 - 127)/2 + 1 : (240 + 127)/2; % 中间 127 个子载波 ssb_grid(pss_idx, 1) = pss; ssb_grid(pss_idx, 3) = sss; % PBCH 和 DMRS 在符号 1/2/3 的剩余位置分别填充

nrPSS和nrSSS是 5G Toolbox 的序列生成函数，也可以按 TS 38.211 的 7.4.2 节自己写生成式，核心就是 m 序列的移位映射。PSS 的循环移位是43 × N_ID^2对 127 取模，SSS 的 m0/m1 由 N_ID^1 查表得到，不要自己编规律。PSS 和 SSS 的映射位置必须是 127 个点，且放在 240 个子载波正中间，这样才能保证接收端抽取时不需要知道 SSB 的频域绝对位置以外的额外信息。

3.2 PSS 检测：用 xcorr 做滑动相关，先拿到定时和粗频偏

PSS 检测的输入是基带接收波形，输出是符号定时估计和粗频偏。最简单可靠的做法是时域滑动相关：本地准备 3 个 PSS 候选，与接收波形做互相关，找全局峰值。

% 接收信号 rx，已下变频到基带 pss_cands = {pss_0, pss_1, pss_2}; % 3 个 PSS 假设，由 3.1 序列循环移位得到 peak_val = 0; for cand = 1:3 [acor, lags] = xcorr(rx, pss_cands{cand}); % 一次算完整互相关 [v, idx] = max(abs(acor)); if v > peak_val peak_val = v; nid2_est = cand - 1; % PSS 假设号即 N_ID^2 估计 timing_est = lags(idx); % 符号定时，单位是采样点 phase_est = angle(acor(idx)); % 峰值相位，含频偏信息 end end % 粗频偏换算：delta_f = -phase_est / (2 * pi * Ts)

逻辑说明：xcorr返回所有时延点的相关值，省掉逐点循环。timing_est是峰值对应的时延，折算成采样点后就是 OFDM 符号起始位置；phase_est里混着频偏和定时偏差，需要按符号时长换算成频偏。nid2_est是三个 PSS 假设里胜出的那个。实际工程里，这一步之后要先做一次频偏补偿再去做 SSS，否则 SSS 的频域相关峰会明显变矮。

参数上要留意：PSS 相关归一化应该用接收信号能量而不是本地序列能量，否则噪声大的时段会出现假峰。频偏较大时（比如超过子载波间隔的一半），时域相关峰会展宽，定时估计会偏向一侧，这是正常现象，交给后面的 CP 相关或 SSS 精估计去修正。

3.3 SSS 假设检验与联合成功率统计口径

PSS 定时确定后，把接收波形按定时做 FFT 抽取符号 2 的频域数据，取中间 127 个子载波，和 336 个 SSS 假设做互相关，选最大峰。SSS 的互相关特性比 PSS 差一些，因为它是两个 m 序列相乘的复合序列，所以低 SNR 下更容易选错假设。

% 在 PSS 定时基础上抽取符号 2 的频域向量 sym2 = fft(rx_waveform(timing_est + sym_start(2) : ... timing_est + sym_start(2) + N_fft - 1)); sym2 = sym2(ssb_freq_range); % 取中间 127 个子载波 sss_metrics = zeros(1, 336); for nid1 = 0:335 sss_local = nrSSS(nid1); sss_metrics(nid1 + 1) = abs(sym2.' * conj(sss_local)) / ... (norm(sym2) * norm(sss_local) + eps); end [~, sss_idx] = max(sss_metrics); nid1_est = sss_idx - 1; pci_est = 3 * nid1_est + nid2_est;

联合检测成功率的统计口径要这样定义：pci_est == pci_true且abs(timing_est - timing_true)小于容限（比如 CP 长度的一半）。只有 PCI 猜对而定时差很多，SSB 后面的 PBCH 照样解不出来。统计脚本里把每个块的真实 PCI、真实定时和估计值全部存下来，最后统一汇总：

success = (pci_est == pci_true) & (abs(timing_est - timing_true) < tol); detect_rate = mean(success) * 100; % 联合检测成功率，单位 %

容限tol必须按采样率折算。仿真里加了频偏后，PSS 定时估计会有系统性偏移，如果发现检测成功率在低 SNR 下掉得异常快，先检查频偏补偿是否在 SSS 检测之前生效。另外一个工程技巧是：单个 SSB 检测成功率不够时，可以做跨帧的多次检测再投票，但那是后话，先把单次链路跑通。

4. PBCH 解调与 polar 译码：BER 和 BLER 怎么算才不虚

4.1 DMRS 信道估计与 QPSK 软解调 LLR

PBCH 的 432 个 QPSK 数据符号散布在符号 1、符号 2 两侧、符号 3 中，DMRS 按每 4 个 RE 一个的密度插入。信道估计的常见做法是先取 DMRS 位置的接收值与本地参考序列做 LS 估计，再插值到数据位置。

% DMRS 位置提取与 LS 信道估计 dmrs_rx = rx_grid(pbch_dmrs_indices); % DMRS 位置由 N_ID 决定 dmrs_ref = nrPBCHDMRS(ncellid); % 本地 DMRS 参考 hest = dmrs_rx .* conj(dmrs_ref) ./ (abs(dmrs_ref).^2 + eps); % LS 估计 hest_data = interp1(..., hest, pbch_data_indices); % 插值到数据 RE

逻辑说明：nrPBCHDMRS生成与 PCI 相关的参考序列，dmrs_rx是资源网格上对应位置的接收值。LS 估计得到每个 DMRS RE 上的信道响应，再插值铺到数据 RE 上。频率选择性信道下线性插值基本够用；多径时延扩展大时可以换样条插值，代价是复杂度上升。噪声方差noise_var必须和 SNR 统计口径保持一致，否则 LLR 幅度失真，polar 译码性能会整体偏差。

QPSK 软解调把 432 个符号拆成 864 个比特的 LLR，近似做法是把 QPSK 看成两路正交的 BPSK：

% 均衡后的符号按实部虚部分别计算 LLR eq = rx_data .* conj(hest_data) ./ (abs(hest_data).^2 + noise_var); llr = [real(eq); imag(eq)] * (2 / (sqrt(2) * noise_var)); % 结果是一个 864 长的 LLR 向量，对应 864 个编码比特

这里2 / (sqrt(2) * noise_var)是 QPSK 星座点距离和噪声方差的折算系数。有人图省事把噪声项直接填 1，曲线看起来也能跑，但换 SNR 口径后完全不可复现，属于典型的血泪经验。

4.2 CA-SCL 译码器配置：list 大小与去速率匹配

LLR 向量拿到后，先做去速率匹配，再进 CA-SCL 译码器。自研译码器时，我通常把 list 大小和 CRC 校验开关做成可配置参数。

% CA-SCL 译码：L=8，CRC24C 辅助路径选择 [info_bits, crc_ok] = ca_scl_polar_decode(... llr_dematted, ... % 去速率匹配后的 512 长 LLR n_polar=512, ... % 母码长度 k_info=32, ... % 信息比特数 crc_len=24, ... % CRC24C list_size=8); % 候选路径数

逻辑说明：译码器内部按 polar 因子图做 SC 展开，维护list_size条活跃路径，每扩展一层按路径度量排序剪枝，最终输出时用 CRC24C 从剩余路径中挑选合法的一条。crc_ok直接作为这个码块是否译码成功的判据，BLER 的分子就靠它。注意llr_dematted不是原始 LLR，而是经过解交织、重复位合并、冻结位置零之后的 512 长向量。

list 大小从 8 加到 16，BLER 通常能下移 0.2~0.5 dB，但译码时间接近翻倍。仿真预算有限时先在 L=8 上跑通，最后再决定是否提升到 16。

4.3 BER vs BLER：统计口径和蒙特卡洛点数

仿真输出三条曲线：PSS-SSS 联合检测成功率、BER、BLER。检测成功率按“PCI 和定时同时正确”的比例统计；BER 按译码输出的 32 位信息比特与原始 32 位逐位对比；BLER 按crc_ok == 0的比例统计。BER 和 BLER 不能互相换算：一个块错 1 个比特和错 10 个比特，对 BLER 的贡献完全一样。

蒙特卡洛点数决定曲线可信度。想看 BLER=10^-2，至少需要跑出 100 个错误块，也就是约 1 万个块；每个块做一次 4096 点 FFT 加 CA-SCL 译码，仿真时间会相当可观。常见做法是低 SNR 段少跑，高 SNR 段设置最大块数上限，按累计错误块数自适应停止。BER 的最小统计量建议比 BLER 更大，因为一个错误块只贡献一个 BLER 错误计数，但可能贡献几十个比特错误。BER 曲线通常比 BLER 曲线更快收敛，别被假象骗了。

5. 参数调优与避坑：三条曲线对不上时先查这五处

5.1 PSS 相关峰出现等间隔旁瓣，定时估计周期性跳变

现象：检测成功率在中等 SNR 下无故下降，打印定时误差发现估计值与真实值之差总在正负一个 OFDM 符号附近跳。原因：PSS 序列的循环移位自相关旁瓣不低，加上 CP 结构，在没做频偏预校正时，相关峰会周期性出现在符号边界上，最大峰值选到了旁瓣。解决：在 PSS 相关之前先用 CP 自相关做一次粗定时，把搜索窗口缩到 CP 范围内，再做 PSS 精相关；或者在频域做 PSS 相关时对序列做加窗处理压制旁瓣。

5.2 残余频偏让星座整体旋转，polar 译码雪崩

现象：PSS-SSS 联合检测成功率很高，98% 以上，但 PBCH 的 BLER 在低 SNR 下明显高于同参数的理论曲线，甚至出现误码平台。原因：PSS 粗频偏估计精度不够，残余频偏使 QPSK 星座缓慢旋转，DMRS 估计出的信道跟不上符号间相位变化，LLR 幅度被压缩。解决：在 SSS 确认之后用 SSS 或 DMRS 做一次精细频偏估计，对 PBCH 符号做逐符号相位补偿；仿真里把“理想频偏补偿”和“实际频偏补偿”做成两条 BLER 曲线，两条线的差就是频偏估计模块的优化空间。这个问题最隐蔽，因为它不影响检测成功率，只影响端到端 BLER。

5.3 SNR 口径不统一，曲线整体平移 2~3 dB

现象：同一套链路今天跑出的 BLER 曲线和上周差 2 dB，或者和标准参考曲线差 3 dB。原因：SNR 定义口径混乱——有人按信号总功率/噪声总功率，有人按每个 RE 的 Es/N0，还有人按 E_b/N0。QPSK 下 Es/N0 和 E_b/N0 差 3 dB，如果再把 864 个编码比特折算成 32 个信息比特，差距更大。解决：仿真脚本里固定用 Es/N0（每数据 RE 符号能量对噪声功率谱密度），在配置文件里显式打印 SNR 换算关系，确保每条曲线都能被别人复现。

5.4 蒙特卡洛点数不够，BLER 低段乱跳

现象：BLER 在 10^-3 附近出现强烈的锯齿，同样参数跑两遍结果差一倍。原因：高 SNR 下错误块是稀有事件，固定块数仿真可能只碰到一两个错误块，统计方差极大。解决：按累计错误块数自适应停止，至少累计 100 个错误块；对 10^-3 量级目标，预先估算需要的最大块数并设max_blocks上限，避免无限跑下去。检测成功率统计同理，想看 99% 的检测成功率至少跑 1000 次独立 SSB 检测。

5.5 信息位集合不按标准来，高 SNR 出现错误平台

现象：BLER 曲线在高 SNR 下不再往下走，稳定在 10^-4~10^-5 的平台，怎么加大 list 都没用。原因：自研 polar 编码器用的信息位集合是经典巴氏参数排序或旧版本计算，与 NR 标准指定的 Q 函数量化表不完全一致，个别本应冻结的比特被当成信息位，这些位置恰好是极化质量最差的，产生不可约的错误。解决：用 5G Toolbox 的编码器做交叉验证，把信息位集合导出成文本文件固化下来，每次跑仿真前校验一次，避免静默错误。这个坑一旦踩进去，排查周期通常以周为单位。

6. 代码怎么组织才能稳定复现三条曲线：工程结构与验证顺序

6.1 模块划分与主循环

仿真工程按发射、同步、解调译码、统计四个方向拆。主循环是 SNR 外层、块内层，块循环适合并行加速：

parfor blk = 1:max_blocks [det_ok, ber, bler] = sim_one_block(snr_esn0, cfg); results(blk, :) = [det_ok, ber, bler]; end

每个块独立生成 SSB、加信道、走完整接收链路，块之间不共享状态，这是并行仿真的前提。配置项集中在cfg结构体里，包括频偏、定时偏移、信道类型、list 大小、统计容限。打印每个阶段的尺寸链A→B→N→E→symbols，排障时能省一半时间。

6.2 先验证译码器，再逐级打开同步链路

建议按“理想同步 → 半联合 → 全联合”的顺序解锁。第一步，关闭频偏和定时偏移，直接把发射端的 PBCH 符号对齐到接收端，单独验证 polar 译码器和统计脚本，BER/BLER 应当和理想曲线一致。第二步，只加频偏不加定时偏移，验证频偏估计与补偿模块。第三步，全链路跑通，对比三条曲线。每一步留一个 golden 输出文件，防止后续改动把已验证的结果弄坏。

6.3 操作视频对应的演示路径

代码操作视频按这个演示路径录制：先跑单次仿真脚本，打印每个阶段的尺寸链和中间估计值；再跑蒙特卡洛脚本，展示三条曲线随 SNR 从 -5 dB 到 3 dB 的变化；最后演示如何把残余频偏注入解调模块，观察 BLER 劣化曲线。视频里出现的参数和这里的cfg结构体一一对应，照着改就能复现。

我做这个仿真最深刻的一次教训是：检测成功率做到 99.5% 之后得意了很久，结果发现 PBCH 一直没解扰，BER 卡在 50% 左右。那周我天天在查 polar 译码器，最后发现是扰码的问题。从那以后我给自己定了个规矩：所有模块先在理想条件下单测，再逐级打开非理想因素，每一步留 golden 输出。这个习惯帮我少走了很多弯路，希望帮到你。

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