5G极化码解码工程实战：SC/SCL/CA-SCL选型与避坑指南

简介本资源面向5G通信物理层学习者与极化码算法研究者聚焦5G NR控制信道所采用的Polar码译码实现适合具备一定通信原理与MATLAB编程基础的中高级读者用于理解信道极化、编码构造与译码算法的完整链路。压缩包共19个文件约56KB以m脚本为主体辅以mat数据文件、cpp源码、mexw64动态库及pdb调试文件分别承担算法实现、信道与码字数据存储、编译加速与调试等用途。资源围绕极化码核心流程展开涵盖信息序列生成、极化变换矩阵构造、CRC校验、加噪与信道传输以及SC逐位取消译码和基于mex的快速译码实现可帮助读者对照代码理解可靠比特选取、似然比计算与低复杂度译码的工程落地方式。目前已有427人学习下载适合作为5G极化码译码入门与二次开发的参考素材。1. 从一份 polar_code_decode 压缩包说起5G 极化码解码到底难在哪如果你在搜索引擎里敲下polar_code_decode.rar、5G decode、5G极化码、Polar 5g、polar code这几个词大概率你手里已经拿到了一份别人打包好的极化码解码工程或者正准备自己从零搭一套。极化码Polar Code是 5G 控制信道 eMBB 场景下短码块的编码方案uplink 和 downlink 的控制信息都靠它兜底。它不像 LDPC 那样天生适合高吞吐极化码的看点在于码长可以做到很短比如 32、64、128在低信噪比下靠 SC 系列译码能逼近容量。但真正上手你会发现polar_code_decode这类工程最劝退的地方不是理论而是「SC 和 SCL 到底选哪个」「CRC 辅助怎么接」「LLR 怎么算」「冻结位表从哪来」这四件事。这篇笔记就按我实际复现一份极化码解码工程的顺序把选型、代码、参数和踩坑一次讲透适合已经懂一点通信原理、想动手跑通 5G 极化码解码的工程师。2. 5G 极化码解码的工程选型SC、SCL 还是 CA-SCL2.1 三种译码算法的取舍逻辑极化码的译码算法不是越多越好工程上真正会落地的就三条路SCSuccessive Cancellation、SCLSuccessive Cancellation List、CA-SCLCRC-Aided SCL。SC 是基础逐位判决复杂度 O(N log N)但短码下误块率不够看SCL 保留 L 条候选路径L 一般取 4 或 8复杂度翻 L 倍换来明显的纠错增益CA-SCL 在 SCL 基础上把 CRC 校验位拼在信息位后面用 CRC 从 L 条路径里挑一条最像的这是 5G 标准里控制信道实际采用的方案。选型时先问自己三个问题码长多少、目标误块率多少、算力预算多少。码长 N ≤ 128、目标 BLER 在 1e-3 量级CA-SCL L8 基本是标配如果只是做教学演示或者链路预算宽松SC 就够代码量少一半。我一般会先把 SC 跑通确认冻结位、LLR、比特反序这些基础环节没错再往上叠 SCL 的路径管理最后接 CRC。这样出问题的时候能快速定位是哪一层引入的。2.2 冻结位表与 5G 标准参数的对齐极化码的冻结位frozen bit集合决定了哪些位置放信息、哪些位置放已知的 0。5G 标准里冻结位不是随便选的它依赖一个可靠性序列reliability sequence标准给出了从 N32 到 N1024 的固定序列。工程上最常见的翻车就是自己用巴氏参数或者高斯近似算了一套冻结位结果和标准序列对不上解码端和编码端用的不是同一套误块率直接爆掉。参数常见取值说明码长 N32 / 64 / 128 / 256 / 512 / 10245G 控制信道常用 32~1024信息位 K由码率和 CRC 长度决定K N*R - CRC_lenCRC 长度6 / 11 / 24控制信道常用 6 或 11列表宽度 L4 / 8L 越大增益越高算力翻倍冻结位来源5G 标准可靠性序列不要自己算直接查表提示如果你拿到的polar_code_decode工程里冻结位是硬编码的数组先核对它和 5G 标准序列是否一致这是最高频的坑。2.3 从压缩包到可运行工程的最小步骤拿到一份polar_code_decode.rar之后不要急着编译。先解压看目录结构通常会有encoder、decoder、test三个部分。我的习惯是先找测试向量有没有一组已知的输入比特和对应的编码输出。如果有先跑编码确认编码输出对得上再跑解码。没有测试向量就自己造随机生成 K 个信息比特编码后加高斯白噪声再解码看能不能还原。# 解压后先看目录结构不要直接编译 unzip polar_code_decode.rar -d polar_decode cd polar_decode find . -maxdepth 2 -type f | head -50 # 找测试入口通常是 test 或 main 开头的文件 grep -rl main --include*.c --include*.cpp --include*.py .这段命令的逻辑是先摸清工程布局find限制深度避免输出爆炸grep找入口文件。参数上-maxdepth 2是防止递归太深head -50是防止文件太多刷屏。如果你看到的是 Python 工程入口通常是test_polar.py或main.pyC 工程则找main.c。找到入口后先读它的参数解析部分确认码长、码率、SNR 这些是从命令行传还是写死的。3. 用 Python 复现极化码 SC 解码的最小闭环3.1 编码端从信息比特到极化码字极化码编码的核心是极化变换本质是一个蝶形运算网络。把信息比特按冻结位表放到对应位置冻结位置 0然后做 log2(N) 层蝶形异或。下面是最小实现码长 N8信息位 K4冻结位用最简单的固定模式实际工程要换成 5G 标准序列。import numpy as np def polar_encode(info_bits, frozen_positions, N): info_bits: 信息比特长度 K frozen_positions: 冻结位索引集合 N: 码长 u np.zeros(N, dtypeint) info_idx [i for i in range(N) if i not in frozen_positions] for k, idx in enumerate(info_idx): u[idx] info_bits[k] # 蝶形极化变换 x u.copy() n int(np.log2(N)) for stage in range(n): step 1 (stage 1) half step 1 for start in range(0, N, step): for i in range(start, start half): x[i] x[i] ^ x[i half] return x N 8 frozen {0, 1, 2, 4} info [1, 0, 1, 1] codeword polar_encode(info, frozen, N) print(codeword:, codeword)逻辑说明u是长度为 N 的向量信息位填进去冻结位保持 0。蝶形变换按 stage 逐层做每层的步长是 2 的幂次half是半个步长异或操作把左半部分和右半部分合并。参数上frozen_positions用集合是为了 O(1) 查找N必须是 2 的幂。这段代码跑出来的codeword就是极化码字可以拿去加噪声。3.2 解码端SC 逐位判决与 LLR 递推SC 解码是逐位估计每一位的判决依赖前面所有位的估计结果。核心是 LLR 的递推计算分 f 函数和 g 函数两种。f 函数用于左子节点g 函数用于右子节点g 函数还要用到前面位的硬判决结果。def sc_decode(llr, frozen_positions, N): llr: 接收到的对数似然比长度 N frozen_positions: 冻结位索引 n int(np.log2(N)) # 存储每层的 LLR 和比特估计 L np.zeros((n 1, N)) B np.zeros((n 1, N), dtypeint) L[0] llr.copy() # 递归实现简化为迭代版本按位处理 u_hat np.zeros(N, dtypeint) for i in range(N): # 计算第 i 位的 LLR这里用简化递推 # 实际工程用树形递归这里展示逻辑 if i in frozen_positions: u_hat[i] 0 else: u_hat[i] 0 if L[0][i] 0 else 1 return u_hat # 模拟接收 LLR np.random.seed(42) tx 1 - 2 * codeword # BPSK 映射 snr_db 2.0 sigma 10 ** (-snr_db / 20) rx tx sigma * np.random.randn(N) llr 2 * rx / (sigma ** 2) u_hat sc_decode(llr, frozen, N) print(u_hat:, u_hat)逻辑说明llr是接收信号的软信息BPSK 映射下tx 1 - 2*codeword加高斯噪声后算 LLR。sc_decode里冻结位直接判 0信息位按 LLR 符号判决。参数上snr_db控制噪声强度sigma是噪声标准差。这段代码是简化版真实 SC 解码要用树形递归做 f/g 函数递推但逻辑框架就是这样。跑通之后你会看到u_hat和原始info在低 SNR 下会有差异这就是误码。3.3 加 CRC 与 SCL从能跑到能用SC 跑通只是第一步要能用必须上 SCL。SCL 的核心是维护 L 条路径每条路径有自己的 LLR 和比特估计在每一步分裂成两条按路径度量排序保留最好的 L。路径度量PM是累加的 LLR 绝对值PM 越小路径越可靠。def scl_decode(llr, frozen_positions, N, L4): SCL 解码L 为列表宽度 返回 L 条候选路径中 PM 最小的一条 paths [{pm: 0.0, u: np.zeros(N, dtypeint), llr: llr.copy()}] for i in range(N): new_paths [] for p in paths: if i in frozen_positions: # 冻结位不分裂 p[u][i] 0 new_paths.append(p) else: # 分裂成两条 for bit in [0, 1]: new_p {pm: p[pm], u: p[u].copy(), llr: p[llr].copy()} new_p[u][i] bit # PM 更新按 LLR 符号和判决是否一致累加 if (bit 0 and p[llr][i] 0) or (bit 1 and p[llr][i] 0): new_p[pm] abs(p[llr][i]) new_paths.append(new_p) # 按 PM 排序保留 L 条 new_paths.sort(keylambda x: x[pm]) paths new_paths[:L] return paths[0][u] u_hat_scl scl_decode(llr, frozen, N, L4) print(u_hat_scl:, u_hat_scl)逻辑说明paths是候选路径列表每条路径有 PM、比特估计和 LLR。冻结位不分裂信息位分裂成 0 和 1 两条PM 更新规则是如果判决比特和 LLR 符号不一致PM 增加该位的 LLR 绝对值。最后按 PM 排序保留 L 条。参数L控制列表宽度L4 是常见起点L8 增益更高但算力翻倍。CA-SCL 就是在最后用 CRC 从 L 条路径里挑而不是直接取 PM 最小的。4. 极化码解码工程避坑五个血泪教训4.1 冻结位表和标准对不上误块率直接爆现象解码端和编码端用同一套代码但误块率比理论值高一个数量级。原因冻结位集合是自己用巴氏参数算的和 5G 标准可靠性序列不一致。解决直接查 5G 标准给出的可靠性序列按序列从低到高选冻结位不要自己算。这个坑我踩过两次第一次以为是 LLR 算错了查了一天才发现是冻结位的问题。4.2 LLR 符号约定反了解码全错现象SC 解码输出全是 0 或者全是 1完全没纠错能力。原因LLR 的符号约定和判决规则不匹配。有的工程定义 LLR log(P(0)/P(1))有的定义反过来。解决先确认工程里 LLR 的定义判决规则bit 0 if LLR 0 else 1要和定义一致。测试方法不加噪声直接传编码后的 BPSK 信号LLR 应该是大正数或大负数解码应该 100% 正确。4.3 比特反序忘了做索引全乱现象编码输出和解码输入对不上但单独看编码和解码逻辑都没错。原因极化码的蝶形变换和标准里的比特反序bit reversal顺序不一致。解决确认工程里是否做了比特反序5G 标准里极化码的索引顺序是比特反序后的。如果编码端做了反序解码端也要做否则冻结位索引全错。4.4 SCL 路径度量溢出长码下崩溃现象码长 N1024 时 SCL 解码结果异常短码正常。原因PM 是累加的 LLR 绝对值长码下累加值可能溢出浮点精度或者整数范围。解决PM 用浮点数定期做归一化或者用对数域计算。另外路径排序时如果 PM 相等要有稳定的 tie-break 规则否则结果不可复现。4.5 CRC 校验位位置放错CA-SCL 失效现象CA-SCL 和 SCL 结果一样CRC 没起到筛选作用。原因CRC 校验位没有拼在信息位后面或者解码时没有把 CRC 位纳入路径筛选。解决编码时信息位 CRC 位一起做极化编码解码时从 L 条路径里用 CRC 校验选通过校验的那条。如果多条通过选 PM 最小的。5. 把极化码解码跑到可验证进阶技巧与自检习惯5.1 用已知测试向量做端到端自检工程上最可靠的自检方法是端到端测试向量。找一组已知的信息比特、编码输出、加噪后的 LLR、期望的解码输出跑一遍看是否一致。没有现成向量就自己造固定随机种子生成信息比特编码加固定 SNR 的噪声解码记录结果。下次改代码后跑同一组向量结果应该完全一致。def end_to_end_test(N64, K32, snr_db1.5, seed123): 端到端自检编码 - 加噪 - 解码 - 比对 np.random.seed(seed) frozen set(range(N - K)) # 简化冻结位实际用标准序列 info np.random.randint(0, 2, K) codeword polar_encode(info, frozen, N) tx 1 - 2 * codeword sigma 10 ** (-snr_db / 20) rx tx sigma * np.random.randn(N) llr 2 * rx / (sigma ** 2) u_hat scl_decode(llr, frozen, N, L8) info_hat [u_hat[i] for i in range(N) if i not in frozen] errors sum(1 for a, b in zip(info, info_hat) if a ! b) return errors # 跑 100 次统计平均误比特数 total_errors sum(end_to_end_test(seeds) for s in range(100)) print(average bit errors per block:, total_errors / 100)逻辑说明end_to_end_test封装了完整链路固定种子保证可复现。frozen用简化模式实际要换成标准序列。errors统计信息位里的错误数。跑 100 次取平均可以看到在当前 SNR 下的误比特率。参数snr_db控制信道质量seed控制随机性。这个自检脚本我一般会放在工程根目录每次改完解码逻辑先跑一遍。5.2 用 SNR-BLER 曲线验证解码器性能单点测试只能说明「能跑」要验证「跑得好」得画 SNR-BLER 曲线。固定码长和码率扫 SNR 从 -2dB 到 4dB每个点跑足够多的块比如 1000 块统计误块率。CA-SCL L8 的曲线应该明显优于 SC且在 1e-3 附近有陡降。SNR (dB)SC BLERSCL L4 BLERCA-SCL L8 BLER-1.00.450.320.250.00.280.150.081.00.120.040.012.00.030.0050.001注意上表是示意值实际曲线取决于码长、码率和冻结位表。如果你的曲线和这个趋势差太远先查冻结位和 LLR。5.3 我的自检习惯先固定种子再扫参数我自己的习惯是任何极化码解码工程第一步固定随机种子跑端到端测试确认结果可复现第二步扫 SNR 画曲线确认性能趋势合理第三步才去调 L 和 CRC 长度。这样出问题的时候能快速判断是逻辑错误还是参数问题。极化码解码的玄学不多大部分翻车都是冻结位、LLR 符号、比特反序这三个地方。把这三个地方的自检做扎实后面就是调参的体力活。希望帮到你。本文还有配套的精品资源点击获取