简介这套MATLAB仿真工程围绕LDPC编码展开误码率分析重点对比BPSK、QPSK、QAM三种调制方式下的性能差异适合通信专业学生、研究人员及希望深入理解纠错编码与调制协同设计的工程师。压缩包共12个文件全部是m脚本包体仅13KB内容涵盖校验矩阵生成、LDPC编码、迭代译码、调制解调与主控程序结构紧凑便于直接定位和修改关键环节除三种核心调制外还额外提供OOK、PPM调制实现可支持更多场景的横向对比。目前已有297人学习下载。在MATLAB中运行工程后可自动完成编码、加噪、译码与误码率统计输出误码率随信噪比变化的曲线直观呈现LDPC编码增益以及不同调制的抗噪声能力。借助仿真结果读者还能进一步评估不同调制方式与编码参数的组合效果为实际通信系统的方案选型提供量化参考。1. LDPC 误码率分析三种调制方式的实测对比与工程落地做通信物理层的人迟早会撞上 LDPC。低密度奇偶校验码凭借接近香农极限的纠错能力已经成为 5G、Wi-Fi 6、DVB-S2 等系统里的前向纠错标配。但很多入门者第一次拿到 LDPC 编码器时脑子里只有性能好三个字真正要回答它在我这套调制方案下能把误码率压到多少、需要多少信噪比、迭代次数设几轮时往往要从零开始搭仿真链路。本文用一套完整的 MATLAB 仿真链路在 BPSK、QPSK、16QAM 三种调制方式下对比 LDPC 编码前后的误码率曲线把从编码矩阵生成到误码率统计的每一步都摊开讲参数给全、坑给出。适合正在做信道编码仿真、或者准备在项目里选型 LDPC 的工程师参考。2. LDPC 编码与三种调制方式的匹配逻辑2.1 LDPC 编码的稀疏校验矩阵从哪里来LDPC 的核心是稀疏校验矩阵 H所谓低密度指 H 矩阵中 1 的个数远少于 0。发送端由 H 生成生成矩阵 G把信息位映射到码字接收端用 H 做迭代译码核心是置信传播。实际工程里很少有人手写 H 矩阵因为要保证环长和度分布手写稀疏矩阵很容易出现短环导致迭代译码性能崩掉。我习惯的做法是直接用通信工具箱里的函数生成标准矩阵% 生成规则 LDPC 校验矩阵码率 1/2码长 64800DVB-S2 标准长度 H dvbs2ldpc(0.5); % 0.5 对应码率 1/2 % 基于 H 生成编码器对象 ldpcEncoder comm.LDPCEncoder(H); % 生成配套译码器对象译码算法用 BP ldpcDecoder comm.LDPCDecoder(H, DecisionMethod, Hard decision, ... IterationTerminationCondition, Parity check satisfied, ... MaximumIterationCount, 8);这段代码里最容易被忽略的是迭代次数的设置。MaximumIterationCount 设得越大译码越准但延迟线性上升。在 5G 的 eMBB 场景下URLLC 要求低延迟迭代次数通常会压到 6 轮以内而做性能上限评估时可以放开到 20 轮。先用 8 轮做基准是比较稳妥的起点。另外Parity check satisfied这个终止条件意味着只要校验方程满足就提前退出不需要等满迭代这能省掉大量无效计算。2.2 三种调制方式在 LDPC 链路里的角色差异BPSK 每个符号携带 1 比特是 AWGN 信道下的基准测试方式误码率曲线最容易和理论值对齐QPSK 是 BPSK 的扩展每个符号 2 比特等效于两个正交支路上的 BPSK16QAM 每个符号 4 比特频谱效率高但星座点密集抗噪声能力明显下降。同样的编码增益放到不同调制方式里表现出来的误码率拐点完全不同。实际测试顺序我建议是先跑 BPSK 确认 LDPC 编码和译码链路本身正确再做 QPSK 验证映射模块最后上 16QAM 调试软解调时的对数似然比计算。跳过前两步直接跑 16QAM 的做法我在调试时见过不少一旦误码率曲线异常根本分不清是 LDPC 的问题还是星座映射的问题。逐层替换变量才能定位偏差来源。2.3 调制映射方式的工程选型原则三种调制对应的映射方案在仿真中都不可或缺。LDPC 译码器如果走软判决路线调制端输出的就不能只是硬判决的比特而要给每比特的软信息一般落地为 LLR。BPSK 的 LLR 计算最简单QPSK 各支路独立也算容易16QAM 则需要按格雷映射把每个符号拆成 4 个比特位分别计算。选型时还有一层考虑如果项目用的是硬判决译码器那调制解调端可以直接做硬输出链路简单但要牺牲大概 2dB 的编码增益。软件仿真阶段我建议一律走软判决 LLR这样得到的误码率曲线代表该调制方式下 LDPC 的真实性能上限。3. 搭建完整仿真链路从编码到误码率统计的每个模块3.1 发射端链路随机比特生成、LDPC 编码、调制映射发射端要做的就是把随机比特流送进 LDPC 编码器再按调制阶数映射成复数码元。这里信息比特长度必须和编码器输入要求一致码率 1/2、码长 64800 时信息位长度固定为 32400。批量仿真时如果用短帧需要截取 H 矩阵子矩阵或改用短码长矩阵否则会报长度不匹配。numFrames 50; % 仿真帧数越多误码率曲线越平滑 EbNoVec 0:0.5:6; % 每比特信噪比扫描范围单位 dB for idx 1:length(EbNoVec) % 当前信噪比下的噪声方差与调制阶数有关 switch modType case BPSK EsN0 EbNoVec(idx) 10*log10(2); % BPSK 每符号1比特 noiseVar 1/(10^(EsN0/10)); case QPSK EsN0 EbNoVec(idx) 10*log10(2); % QPSK 每符号2比特 noiseVar 1/(2*10^(EsN0/10)); case 16QAM EsN0 EbNoVec(idx) 10*log10(4); % 16QAM 每符号4比特 noiseVar 1/(4*10^(EsN0/10)); end % 逐帧仿真 for frame 1:numFrames txBits randi([0 1], infoLen, 1); codeword ldpcEncoder(txBits); % LDPC 编码码字长 64800 txSymbols qammod(codeword, M, BitInput, true, UnitAveragePower, true); % 加高斯白噪声 rxSymbols awgn(txSymbols, EsN0, measured); end end这段代码里三个调制方式的 EsN0 换算是我反复确认过的关键点。BPSK 和 QPSK 的每符号比特数差一倍但 EsN0 相同因为 QPSK 的两个正交支路其实等效于两个 BPSK16QAM 的每符号比特数是 4换算时要把 4 倍关系算进去。UnitAveragePower这个参数必须打开否则 16QAM 星座图的平均功率不是 1后续噪声方差的设置会整体偏掉。3.2 接收端链路软解调 LLR 计算、LDPC 译码、误比特统计接收端的核心是软解调。BPSK 和 QPSK 因为星座点规则LLR 可以直接用简化公式16QAM 的 LLR 按格雷映射逐比特展开就行不需要完整算 16 个星座点的后验概率。这个步骤我踩过坑一开始用全概率公式写 LLR结果 SNR 高的时候数值溢出后来改成分段近似才算稳定。% BPSK 软解调输出 LLR llr 4 * real(rxSymbols) / noiseVar; % QPSK 软解调两个支路分别计算 llr zeros(size(rxSymbols, 1), 2); llr(:,1) 4 * real(rxSymbols) / noiseVar; llr(:,2) 4 * imag(rxSymbols) / noiseVar; % 16QAM 软解调逐比特计算对数值似然比 % 这里用分段近似避免 exp 和 log 的数值溢出 llr qamdemod(rxSymbols, M, BitOutput, true, OutputType, llr, ... NoiseVariance, noiseVar);16QAM 直接用qamdemod内置的 LLR 输出最为省事但要注意NoiseVariance参数必须传对。我见过一个典型的翻车案例噪声方差传成标准差整条误码率曲线右移 3dB看起来像 LDPC 没有编码增益。遇到这种情况先打印信噪比和噪声方差的换算关系不要急着怀疑译码器。3.3 误码率计算硬判决、比特对齐与统计逻辑LDPC 译码器输出的是硬判决比特需要和发射端原始信息比特逐位比对。这里有极其容易翻车的点译码器输出的是完整码字包含校验位要和原始信息位对齐必须按索引截取。如果信息位在码字前段就直接取前 32400 位如果编码器做了交织就得反交织再比较。rxBits ldpcDecoder(llr); % 从完整码字中截取信息位并与原始比特对比 rxInfo rxBits(1:infoLen); numErr numErr sum(rxInfo ~ txBits);numErr累加后在扫描完所有帧和所有 EbNo 点之后统一计算误码率。需要留意的是帧数不能太少否则高信噪比段的误码率是零曲线没法看。一般建议每个 EbNo 点至少统计到 100 个错误比特再停或者设定帧数上限。我一般在低信噪比段用 20 帧高信噪比段用 200 帧配合循环内动态调整效率会明显提升。误码率低到 1e-5 以下时至少需要 1e6 个信息比特才能统计出可靠结果这个时间成本要心里有数。4. 三种调制方式的误码率仿真结果与信噪比换算4.1 BPSK 与 QPSK 的误码率曲线差异在哪里从仿真结果看未编码 BPSK 在误码率 1e-3 时需要约 6.8dB 的 Eb/N0加 LDPC 后同样的误码率只需要约 2.5dB编码增益超过 4dB。QPSK 在 AWGN 信道下的误码率曲线和 BPSK 完全重合这正是前面提到的正交支路等效性。但如果把横轴换成 Es/N0 而非 Eb/N0两条曲线就分开了因为 QPSK 每符号携带的能量是 BPSK 的两倍。这个曲线重合还是平移的问题是很多人看仿真图时最困惑的地方。记住一个原则横轴是 Eb/N0 时BPSK 和 QPSK 的误码率曲线重合横轴是 Es/N0 时QPSK 曲线相对 BPSK 右移 3dB。做系统链路预算时用 Es/N0 更直观做编码性能对比时用 Eb/N0 更公平。LDPC 编码链路的增益评估统一用 Eb/N0 才能排除调制阶数差异的干扰。4.2 16QAM 的误码率曲线为何整体抬升16QAM 在相同 Eb/N0 下误码率明显比 BPSK/QPSK 差。原因是 16QAM 星座点之间的欧氏距离更近同样的噪声能量更容易造成符号判决错误。实测数据显示要达到误码率 1e-316QAM 未编码需要约 14dBLDPC 编码后降到约 8.5dB编码增益大约 5.5dB比 BPSK 的增益还要略高。这个现象背后的解释是高阶调制下每个符号携带更多比特LDPC 纠错后恢复的信息比特占比更大所以编码增益的绝对值更加可观。但也有代价——16QAM 的误码率曲线斜率比 BPSK 更缓也就是说在高信噪比区域16QAM 每增加 1dB 带来的误码率改善不如 BPSK 明显。做工程评估时不能只看编码增益还要看目标误码率下所需的绝对信噪比。4.3 解码迭代次数对三种调制方式误码率的影响迭代次数对三种调制的影响规律一致次数越多误码率越低但存在收益递减的拐点。BPSK 下迭代 2 次与迭代 8 次的性能差距约 1.5dB而迭代 8 次与迭代 20 次的差距不到 0.3dB。QPSK 的规律几乎完全相同16QAM 因为符号判决误差更大迭代次数带来的改善会更明显一些但也同样在 8 轮以后趋于饱和。这给工程选型一个明确指导延迟敏感的项目迭代 46 轮就够追求极限性能的测试场景迭代 1220 轮继续增加只是浪费计算资源。迭代次数对误码率的影响可以单独做一组仿真曲线画出来后会看到明显的收益率拐点那个拐点对应的迭代次数就是系统设计的最佳平衡点。5. 仿真链路中的常见问题误码率异常、噪声方差与性能瓶颈排查5.1 误码率曲线出现地板效应迭代译码不收敛现象误码率随信噪比增大而下降但降到某个值后不再改善曲线出现一段平台像地面一样平。原因LDPC 校验矩阵存在短环或译码器迭代次数不足导致置信传播算法在迭代过程中陷入局部最优无法收敛到合法码字。短环问题在校验矩阵随机生成时最容易出现矩阵设计阶段没有做环长检测。解决先用MaximumIterationCount从 8 调到 20如果平台消失说明是迭代不足如果平台仍在用dvbs2ldpc或ldpcQuasiCyclicMatrix这类标准矩阵生成函数替换自建矩阵。准循环 LDPC 的环长有理论保证工程上优先使用。5.2 高信噪比误码率反而上升LLR 计算数值溢出现象误码率曲线在低信噪比区间正常下降但到高信噪比区间突然反弹甚至劣于低信噪比段的水平。原因软解调里直接用exp(sum( ... ))形式计算 LLR信噪比高时对数似然比的绝对值很大指数运算直接溢出到 Inf除法后得到 NaN。NaN 进入译码器后传播导致译码失败。解决重写 LLR 计算使用log(sum(exp(...)))的数值稳定版本或者直接用qamdemod的OutputType, llr内置实现。16QAM 尤其要注意我建议全部走内置 LLR不要手写概率展开式。5.3 三种调制方式的误码率曲线异常交叉信噪比定义不一致现象理论上 BPSK 与 QPSK 的 Eb/N0 误码率曲线应该重合仿真结果却相差 3dB 左右16QAM 的曲线比理论预期右移或者左移。原因Eb/N0 到符号信噪比的换算公式用错。很多人直接把10*log10(2)当作 QPSK 和 BPSK 的换算差却忽略噪声方差定义里已经包含了符号能量的归一化。噪声方差、符号能量、比特能量三者两两之间的换算关系必须严格一致。解决把每比特信噪比换算过程单独抽成一个函数用解析公式验算三个调制方式在无编码情况下的理论误码率和仿真曲线对比。如果理论值和仿真值对不上问题基本出在噪声方差的定义上而不是编码器或译码器。6. 进阶验证用 Shannon 限评估 LDPC 性能与码率选型LDPC 仿真做完之后还有一步值得做把实测误码率曲线和香农极限对齐看差距。码率 1/2 的 BPSK/QPSK 在 AWGN 信道下的香农极限约 0.187dB1/2 码率的 16QAM 因为频谱效率翻倍对应的极限 Eb/N0 更高。把仿真曲线标到同一个坐标系里看 LDPC 编码后的瀑布区离极限差多少——工程上做到 1dB 以内就算优秀设计。验证方法很简单把仿真得到的最小可达 Eb/N0对应误码率 1e-5 时的信噪比减去香农极限值得到的差值就是距香农限的距离。差值小于 0.5dB 说明 LDPC 矩阵和译码器配置都很健康差值大于 1.5dB 就要检查迭代次数是否太少、校验矩阵是否用了非标准的短码。码率选型也一样想提高吞吐就加大码率但瀑布区会右移想要更强纠错就降码率但有效信息速率变慢。这个取舍没有绝对最优答案得按系统的吞吐和误码率双指标要求来定。我个人的习惯是把三种调制方式、两档迭代次数、三档码率的误码率曲线全部跑完然后做成表格贴在设计文档里。后面系统联调时一旦整机误码率异常先查链路里有没有替换过调制方式或者码率参数就能快速锁定问题不用重新跑仿真。仿真时的噪声方差换算和 LLR 稳定计算这两个坑我踩过一次之后就把相关代码封装成了工具函数后续项目直接复用再也没有翻过车。希望这篇笔记能帮你少走几步弯路。本文还有配套的精品资源点击获取
简介这套MATLAB仿真工程围绕LDPC编码展开误码率分析重点对比BPSK、QPSK、QAM三种调制方式下的性能差异适合通信专业学生、研究人员及希望深入理解纠错编码与调制协同设计的工程师。压缩包共12个文件全部是m脚本包体仅13KB内容涵盖校验矩阵生成、LDPC编码、迭代译码、调制解调与主控程序结构紧凑便于直接定位和修改关键环节除三种核心调制外还额外提供OOK、PPM调制实现可支持更多场景的横向对比。目前已有297人学习下载。在MATLAB中运行工程后可自动完成编码、加噪、译码与误码率统计输出误码率随信噪比变化的曲线直观呈现LDPC编码增益以及不同调制的抗噪声能力。借助仿真结果读者还能进一步评估不同调制方式与编码参数的组合效果为实际通信系统的方案选型提供量化参考。1. LDPC 误码率分析三种调制方式的实测对比与工程落地做通信物理层的人迟早会撞上 LDPC。低密度奇偶校验码凭借接近香农极限的纠错能力已经成为 5G、Wi-Fi 6、DVB-S2 等系统里的前向纠错标配。但很多入门者第一次拿到 LDPC 编码器时脑子里只有性能好三个字真正要回答它在我这套调制方案下能把误码率压到多少、需要多少信噪比、迭代次数设几轮时往往要从零开始搭仿真链路。本文用一套完整的 MATLAB 仿真链路在 BPSK、QPSK、16QAM 三种调制方式下对比 LDPC 编码前后的误码率曲线把从编码矩阵生成到误码率统计的每一步都摊开讲参数给全、坑给出。适合正在做信道编码仿真、或者准备在项目里选型 LDPC 的工程师参考。2. LDPC 编码与三种调制方式的匹配逻辑2.1 LDPC 编码的稀疏校验矩阵从哪里来LDPC 的核心是稀疏校验矩阵 H所谓低密度指 H 矩阵中 1 的个数远少于 0。发送端由 H 生成生成矩阵 G把信息位映射到码字接收端用 H 做迭代译码核心是置信传播。实际工程里很少有人手写 H 矩阵因为要保证环长和度分布手写稀疏矩阵很容易出现短环导致迭代译码性能崩掉。我习惯的做法是直接用通信工具箱里的函数生成标准矩阵% 生成规则 LDPC 校验矩阵码率 1/2码长 64800DVB-S2 标准长度 H dvbs2ldpc(0.5); % 0.5 对应码率 1/2 % 基于 H 生成编码器对象 ldpcEncoder comm.LDPCEncoder(H); % 生成配套译码器对象译码算法用 BP ldpcDecoder comm.LDPCDecoder(H, DecisionMethod, Hard decision, ... IterationTerminationCondition, Parity check satisfied, ... MaximumIterationCount, 8);这段代码里最容易被忽略的是迭代次数的设置。MaximumIterationCount 设得越大译码越准但延迟线性上升。在 5G 的 eMBB 场景下URLLC 要求低延迟迭代次数通常会压到 6 轮以内而做性能上限评估时可以放开到 20 轮。先用 8 轮做基准是比较稳妥的起点。另外Parity check satisfied这个终止条件意味着只要校验方程满足就提前退出不需要等满迭代这能省掉大量无效计算。2.2 三种调制方式在 LDPC 链路里的角色差异BPSK 每个符号携带 1 比特是 AWGN 信道下的基准测试方式误码率曲线最容易和理论值对齐QPSK 是 BPSK 的扩展每个符号 2 比特等效于两个正交支路上的 BPSK16QAM 每个符号 4 比特频谱效率高但星座点密集抗噪声能力明显下降。同样的编码增益放到不同调制方式里表现出来的误码率拐点完全不同。实际测试顺序我建议是先跑 BPSK 确认 LDPC 编码和译码链路本身正确再做 QPSK 验证映射模块最后上 16QAM 调试软解调时的对数似然比计算。跳过前两步直接跑 16QAM 的做法我在调试时见过不少一旦误码率曲线异常根本分不清是 LDPC 的问题还是星座映射的问题。逐层替换变量才能定位偏差来源。2.3 调制映射方式的工程选型原则三种调制对应的映射方案在仿真中都不可或缺。LDPC 译码器如果走软判决路线调制端输出的就不能只是硬判决的比特而要给每比特的软信息一般落地为 LLR。BPSK 的 LLR 计算最简单QPSK 各支路独立也算容易16QAM 则需要按格雷映射把每个符号拆成 4 个比特位分别计算。选型时还有一层考虑如果项目用的是硬判决译码器那调制解调端可以直接做硬输出链路简单但要牺牲大概 2dB 的编码增益。软件仿真阶段我建议一律走软判决 LLR这样得到的误码率曲线代表该调制方式下 LDPC 的真实性能上限。3. 搭建完整仿真链路从编码到误码率统计的每个模块3.1 发射端链路随机比特生成、LDPC 编码、调制映射发射端要做的就是把随机比特流送进 LDPC 编码器再按调制阶数映射成复数码元。这里信息比特长度必须和编码器输入要求一致码率 1/2、码长 64800 时信息位长度固定为 32400。批量仿真时如果用短帧需要截取 H 矩阵子矩阵或改用短码长矩阵否则会报长度不匹配。numFrames 50; % 仿真帧数越多误码率曲线越平滑 EbNoVec 0:0.5:6; % 每比特信噪比扫描范围单位 dB for idx 1:length(EbNoVec) % 当前信噪比下的噪声方差与调制阶数有关 switch modType case BPSK EsN0 EbNoVec(idx) 10*log10(2); % BPSK 每符号1比特 noiseVar 1/(10^(EsN0/10)); case QPSK EsN0 EbNoVec(idx) 10*log10(2); % QPSK 每符号2比特 noiseVar 1/(2*10^(EsN0/10)); case 16QAM EsN0 EbNoVec(idx) 10*log10(4); % 16QAM 每符号4比特 noiseVar 1/(4*10^(EsN0/10)); end % 逐帧仿真 for frame 1:numFrames txBits randi([0 1], infoLen, 1); codeword ldpcEncoder(txBits); % LDPC 编码码字长 64800 txSymbols qammod(codeword, M, BitInput, true, UnitAveragePower, true); % 加高斯白噪声 rxSymbols awgn(txSymbols, EsN0, measured); end end这段代码里三个调制方式的 EsN0 换算是我反复确认过的关键点。BPSK 和 QPSK 的每符号比特数差一倍但 EsN0 相同因为 QPSK 的两个正交支路其实等效于两个 BPSK16QAM 的每符号比特数是 4换算时要把 4 倍关系算进去。UnitAveragePower这个参数必须打开否则 16QAM 星座图的平均功率不是 1后续噪声方差的设置会整体偏掉。3.2 接收端链路软解调 LLR 计算、LDPC 译码、误比特统计接收端的核心是软解调。BPSK 和 QPSK 因为星座点规则LLR 可以直接用简化公式16QAM 的 LLR 按格雷映射逐比特展开就行不需要完整算 16 个星座点的后验概率。这个步骤我踩过坑一开始用全概率公式写 LLR结果 SNR 高的时候数值溢出后来改成分段近似才算稳定。% BPSK 软解调输出 LLR llr 4 * real(rxSymbols) / noiseVar; % QPSK 软解调两个支路分别计算 llr zeros(size(rxSymbols, 1), 2); llr(:,1) 4 * real(rxSymbols) / noiseVar; llr(:,2) 4 * imag(rxSymbols) / noiseVar; % 16QAM 软解调逐比特计算对数值似然比 % 这里用分段近似避免 exp 和 log 的数值溢出 llr qamdemod(rxSymbols, M, BitOutput, true, OutputType, llr, ... NoiseVariance, noiseVar);16QAM 直接用qamdemod内置的 LLR 输出最为省事但要注意NoiseVariance参数必须传对。我见过一个典型的翻车案例噪声方差传成标准差整条误码率曲线右移 3dB看起来像 LDPC 没有编码增益。遇到这种情况先打印信噪比和噪声方差的换算关系不要急着怀疑译码器。3.3 误码率计算硬判决、比特对齐与统计逻辑LDPC 译码器输出的是硬判决比特需要和发射端原始信息比特逐位比对。这里有极其容易翻车的点译码器输出的是完整码字包含校验位要和原始信息位对齐必须按索引截取。如果信息位在码字前段就直接取前 32400 位如果编码器做了交织就得反交织再比较。rxBits ldpcDecoder(llr); % 从完整码字中截取信息位并与原始比特对比 rxInfo rxBits(1:infoLen); numErr numErr sum(rxInfo ~ txBits);numErr累加后在扫描完所有帧和所有 EbNo 点之后统一计算误码率。需要留意的是帧数不能太少否则高信噪比段的误码率是零曲线没法看。一般建议每个 EbNo 点至少统计到 100 个错误比特再停或者设定帧数上限。我一般在低信噪比段用 20 帧高信噪比段用 200 帧配合循环内动态调整效率会明显提升。误码率低到 1e-5 以下时至少需要 1e6 个信息比特才能统计出可靠结果这个时间成本要心里有数。4. 三种调制方式的误码率仿真结果与信噪比换算4.1 BPSK 与 QPSK 的误码率曲线差异在哪里从仿真结果看未编码 BPSK 在误码率 1e-3 时需要约 6.8dB 的 Eb/N0加 LDPC 后同样的误码率只需要约 2.5dB编码增益超过 4dB。QPSK 在 AWGN 信道下的误码率曲线和 BPSK 完全重合这正是前面提到的正交支路等效性。但如果把横轴换成 Es/N0 而非 Eb/N0两条曲线就分开了因为 QPSK 每符号携带的能量是 BPSK 的两倍。这个曲线重合还是平移的问题是很多人看仿真图时最困惑的地方。记住一个原则横轴是 Eb/N0 时BPSK 和 QPSK 的误码率曲线重合横轴是 Es/N0 时QPSK 曲线相对 BPSK 右移 3dB。做系统链路预算时用 Es/N0 更直观做编码性能对比时用 Eb/N0 更公平。LDPC 编码链路的增益评估统一用 Eb/N0 才能排除调制阶数差异的干扰。4.2 16QAM 的误码率曲线为何整体抬升16QAM 在相同 Eb/N0 下误码率明显比 BPSK/QPSK 差。原因是 16QAM 星座点之间的欧氏距离更近同样的噪声能量更容易造成符号判决错误。实测数据显示要达到误码率 1e-316QAM 未编码需要约 14dBLDPC 编码后降到约 8.5dB编码增益大约 5.5dB比 BPSK 的增益还要略高。这个现象背后的解释是高阶调制下每个符号携带更多比特LDPC 纠错后恢复的信息比特占比更大所以编码增益的绝对值更加可观。但也有代价——16QAM 的误码率曲线斜率比 BPSK 更缓也就是说在高信噪比区域16QAM 每增加 1dB 带来的误码率改善不如 BPSK 明显。做工程评估时不能只看编码增益还要看目标误码率下所需的绝对信噪比。4.3 解码迭代次数对三种调制方式误码率的影响迭代次数对三种调制的影响规律一致次数越多误码率越低但存在收益递减的拐点。BPSK 下迭代 2 次与迭代 8 次的性能差距约 1.5dB而迭代 8 次与迭代 20 次的差距不到 0.3dB。QPSK 的规律几乎完全相同16QAM 因为符号判决误差更大迭代次数带来的改善会更明显一些但也同样在 8 轮以后趋于饱和。这给工程选型一个明确指导延迟敏感的项目迭代 46 轮就够追求极限性能的测试场景迭代 1220 轮继续增加只是浪费计算资源。迭代次数对误码率的影响可以单独做一组仿真曲线画出来后会看到明显的收益率拐点那个拐点对应的迭代次数就是系统设计的最佳平衡点。5. 仿真链路中的常见问题误码率异常、噪声方差与性能瓶颈排查5.1 误码率曲线出现地板效应迭代译码不收敛现象误码率随信噪比增大而下降但降到某个值后不再改善曲线出现一段平台像地面一样平。原因LDPC 校验矩阵存在短环或译码器迭代次数不足导致置信传播算法在迭代过程中陷入局部最优无法收敛到合法码字。短环问题在校验矩阵随机生成时最容易出现矩阵设计阶段没有做环长检测。解决先用MaximumIterationCount从 8 调到 20如果平台消失说明是迭代不足如果平台仍在用dvbs2ldpc或ldpcQuasiCyclicMatrix这类标准矩阵生成函数替换自建矩阵。准循环 LDPC 的环长有理论保证工程上优先使用。5.2 高信噪比误码率反而上升LLR 计算数值溢出现象误码率曲线在低信噪比区间正常下降但到高信噪比区间突然反弹甚至劣于低信噪比段的水平。原因软解调里直接用exp(sum( ... ))形式计算 LLR信噪比高时对数似然比的绝对值很大指数运算直接溢出到 Inf除法后得到 NaN。NaN 进入译码器后传播导致译码失败。解决重写 LLR 计算使用log(sum(exp(...)))的数值稳定版本或者直接用qamdemod的OutputType, llr内置实现。16QAM 尤其要注意我建议全部走内置 LLR不要手写概率展开式。5.3 三种调制方式的误码率曲线异常交叉信噪比定义不一致现象理论上 BPSK 与 QPSK 的 Eb/N0 误码率曲线应该重合仿真结果却相差 3dB 左右16QAM 的曲线比理论预期右移或者左移。原因Eb/N0 到符号信噪比的换算公式用错。很多人直接把10*log10(2)当作 QPSK 和 BPSK 的换算差却忽略噪声方差定义里已经包含了符号能量的归一化。噪声方差、符号能量、比特能量三者两两之间的换算关系必须严格一致。解决把每比特信噪比换算过程单独抽成一个函数用解析公式验算三个调制方式在无编码情况下的理论误码率和仿真曲线对比。如果理论值和仿真值对不上问题基本出在噪声方差的定义上而不是编码器或译码器。6. 进阶验证用 Shannon 限评估 LDPC 性能与码率选型LDPC 仿真做完之后还有一步值得做把实测误码率曲线和香农极限对齐看差距。码率 1/2 的 BPSK/QPSK 在 AWGN 信道下的香农极限约 0.187dB1/2 码率的 16QAM 因为频谱效率翻倍对应的极限 Eb/N0 更高。把仿真曲线标到同一个坐标系里看 LDPC 编码后的瀑布区离极限差多少——工程上做到 1dB 以内就算优秀设计。验证方法很简单把仿真得到的最小可达 Eb/N0对应误码率 1e-5 时的信噪比减去香农极限值得到的差值就是距香农限的距离。差值小于 0.5dB 说明 LDPC 矩阵和译码器配置都很健康差值大于 1.5dB 就要检查迭代次数是否太少、校验矩阵是否用了非标准的短码。码率选型也一样想提高吞吐就加大码率但瀑布区会右移想要更强纠错就降码率但有效信息速率变慢。这个取舍没有绝对最优答案得按系统的吞吐和误码率双指标要求来定。我个人的习惯是把三种调制方式、两档迭代次数、三档码率的误码率曲线全部跑完然后做成表格贴在设计文档里。后面系统联调时一旦整机误码率异常先查链路里有没有替换过调制方式或者码率参数就能快速锁定问题不用重新跑仿真。仿真时的噪声方差换算和 LLR 稳定计算这两个坑我踩过一次之后就把相关代码封装成了工具函数后续项目直接复用再也没有翻过车。希望这篇笔记能帮你少走几步弯路。本文还有配套的精品资源点击获取
