基于AWGN信道的LDPC编解码仿真系统

资 源 简 介

该项目是一个专门用于MATLAB环境的LDPC（低密度奇偶校验码）编解码仿真方案，全面实现从信源产生到信道传输及终端译码的完整链路。系统首要功能是生成符合特定规则的稀疏校验矩阵，包括随机构造与规则构造等多种模式，确保矩阵具备低密度特性。在编码阶段，程序利用生成的矩阵对原始信息比特流进行高效编码处理；信道模型模拟了典型的加性高斯白噪声（AWGN）环境，支持自定义信噪比（SNR）参数。项目的核心在于译码环节，实现了基于置信传播（Belief Propagation）原理的和积算法（SPA），通过在校验节点与变量节点间迭代传递概率信息来纠正传输错误。程序还内置了对数域和积算法（Log-SPA）以提升计算效率并避免数值下溢风险。用户可以通过该程序自动运行多组仿真实验，获取不同信噪比下的误比特率（BER）性能曲线，从而验证LDPC码在接近香农极限方面的优异表现。该程序结构清晰，注释详尽，适用于无线通信、卫星通信研究以及数字编码课程的实验教学。

详 情 说 明

基于AWGN信道的LDPC低密度奇偶校验码MATLAB仿真系统

本项目提供了一个完整的LDPC（低密度奇偶校验码）编解码链路仿真环境，旨在通过MATLAB实现从矩阵构造、信号编码、信道传输到迭代译码的全过程。系统专注于展示LDPC码在加性高斯白噪声（AWGN）信道下的纠错性能，并提供了详细的性能评估可视化工具。

项目介绍

该系统通过程序化方式实现了LDPC码的经典模型。项目不仅涵盖了规则校验矩阵的构造算法，还集成了高效的对数域和积算法（Log-SPA）进行译码。用户可以通过调整预设的参数（如码长、码率、信噪比范围及迭代次数），观察系统在不同环境下的误比特率（BER）表现。系统输出包括校验矩阵的稀疏性视图、BER随信噪比变化的曲线图，以及译码过程中的收敛情况分析图。

功能特性

1. 规则矩阵构造 采用Gallager构造法生成规则LDPC校验矩阵。系统首先生成基础分块矩阵，满足预设的列重（wc=3）和行重（wr=6），随后通过随机列置换操作确保矩阵的随机性与稀疏性。

2. 矩阵系统化处理 为了便于编码，系统利用模2算术下的高斯消元法将原始校验矩阵转换为系统码形式 [P | I]。基于此形式，程序能够自动派生出对应的生成矩阵 G，从而简化了从信息位到码字的映射过程。

3. 对数域和积译码算法 (Log-SPA) 译码核心采用了数值稳定的对数域和积算法。相比于概率域的SPA算法，该实现利用对数似然比（LLR）进行计算，将大量的乘法运算转换为加法运算，有效避免了浮点数下溢问题，并显著提升了计算效率。

4. 自动化仿真与评估 系统支持多信噪比点的自动循环测试。每一信噪比点均采用蒙特卡罗方法，直到达到设定的最小错误计数或仿真帧数上限。最终生成标准的BER曲线，并提供单帧译码过程中的校验综合征（Syndrome）收敛分析。

实现逻辑与流程

第一阶段：初始化与矩阵准备 设定码长（N=504）和信息长度（K=252），计算校验位长度。通过Gallager法构造校验矩阵 H。对 H 进行高斯消元变换，确保其右侧为单位阵，并据此构造生成矩阵 G。

第二阶段：信道传输模拟 生成随机二进制信息序列，利用生成矩阵 G 完成编码。采用BPSK调制方案将 0/1 序列映射为 +1/-1。根据当前的信噪比计算噪声标准差，在信号中叠加 AWGN 噪声。

第三阶段：迭代译码处理 计算接收端信道对数似然比 L_c。在译码循环中，交替进行校验节点更新（利用 tanh 和 atanh 函数处理概率传递）和变量节点更新。每一轮迭代后，利用当前的判定结果进行综合征检测（H * c' = 0），若满足条件则提前终止迭代。

第四阶段：数据收集与绘图 对比原始信息位与译码后的信息位，统计误码数。在完成所有信噪比点测试后，绘制半对数坐标下的 BER 性能曲线。

关键算法细节

校验节点更新（Check Node Update）： 在对数域中，校验节点向变量节点传递的信息通过多个变量节点传入信息的 tanh 函数乘积计算得出。程序中加入了截断处理（1e-15），以防止 atanh 函数在处理极值时出现数值溢出。

变量节点更新（Variable Node Update）： 变量节点负责汇总所有连接的校验节点反馈的信息以及初始化的信道信息。通过简单的加减运算即可完成信息反馈的更新，体现了对数域计算的便捷性。

提前停止准则（Early Stopping Criterion）： 通过实时计算校验综合征，如果当前解满足 H*c' = 0，说明已找到合法码字，程序会立即退出当前帧的迭代，节省仿真时间。

系统要求

• MATLAB R2016b 或更高版本。
• 无需额外工具箱支持，核心逻辑基于标准 MATLAB 矩阵运算实现。

使用方法

1. 打开 MATLAB 软件。
2. 将相关代码脚本载入编辑器。
3. 点击运行按钮或在命令行输入主程序指令。
4. 程序将依次弹出三个窗口：

* 校验矩阵稀疏特性图：展示矩阵 H 中非零元素的分布。 * LDPC码性能曲线：展示不同信噪比下的误比特率。 * 译码迭代过程收敛分析：展示最后一组信噪比下，校验不匹配数随迭代次数下降的过程。

1. 在命令行窗口可以实时查看当前仿真的进度、总比特数、错误数及实时 BER 数值。