LDPC编码在CDMA移动通信系统中的性能仿真与应用研究
项目介绍
本项目实现了一个基于MATLAB的综合仿真系统,专门用于评估低密度奇偶校验码(LDPC)在直接序列码分多址(DS-CDMA)移动通信系统中的性能表现。系统重点模拟了多用户并发环境下的物理层通信链路,通过引入准循环LDPC(QC-LDPC)构造和置信传播译码算法,深入探讨其在抑制多址干扰(MAI)以及克服加性高斯白噪声(AWGN)和瑞利(Rayleigh)衰落对通信质量影响的能力。
功能特性
- 准循环构造方案:支持通过基础矩阵和位移因子快速构造大规模校验矩阵,生成符合工程实际需求的QC-LDPC码。
- 多用户并发仿真:模拟4个用户共享信道的CDMA场景,利用正交Walsh码实现地址复用。
- 混合信道模型:系统完整实现了AWGN信道与多径平坦瑞利衰落信道的对比仿真。
- 高性能和积译码:采用在对数域运行的和积算法(Log-SPA),避免了乘法运算带来的数值溢出,并极大提高了译码效率。
- 多维统计分析:自动生成比特误码率(BER)随信噪比(SNR)变化的曲线,并动态定量分析译码过程中LLR(对数似然比)的收敛趋势。
实现逻辑与流程
本项目通过主运行模块串联了从信源编码到信宿解码的全过程,具体执行逻辑如下:
- 环境初始化与参数定义:设置信噪比范围(0-6 dB)、扩频因子(16)、LDPC码率(0.5)以及仿真规模。定义QC-LDPC的基础矩阵,设定子块大小为48,最终生成(288, 144)维度的系统码。
- 校验与生成矩阵构造:
- 根据循环移位参数构建满秩的校验矩阵H。
- 对H进行模二运算下的高斯消元,将其转换为标准阵形式。
- 从标准阵提取校验分量,构造出满足 $GH^T=0$ 的生成矩阵G,确保实现系统码编码(即信息位直接出现在码字前部)。
- 发射机信号处理:
- 随机产生多用户原始信息位。
- 执行LDPC编码,将信息位映射为具有纠错能力的码字。
- 进行BPSK映射,将二进制0/1转换为+1/-1电平信号。
- 使用Hadamard矩阵生成的正交Walsh码对每个用户的符号进行展宽,完成扩频。
- 将所有用户的扩频信号在时域进行线性叠加。
- 信道传输模拟:
- 在AWGN信道中,为合成信号添加特定功率的随机噪声。
- 在瑞利信道中,产生瑞利分布的衰减因子,并模拟多径平坦衰落过程中的信号畸变。
- 接收机信号处理:
-
解扩:利用预分配的各用户正交码,对接收信号做相关运算,抑制其他用户的多址干扰,恢复基带软信息。
-
似然计算:根据当前信道的噪声方差,将解扩后的信号转换为LDPC译码器所需的初始对数似然比(LLR)。
-
迭代译码:将LLR送入Log-SPA译码器,在校验节点与变量节点间传递置信信息,通过最大20次迭代逼近真实码字。
- 性能定量分析:统计不同信道下的误码数量,计算各信噪比点的BER,并记录译码过程中LLR平均置信度的增强过程。
关键函数与算法说明
1. 准循环构造算法
算法通过预定义的基础位移矩阵 B,对单位矩阵进行循环右移(circshift)操作。这种结构化的构造方式不仅降低了存储校验矩阵所需的内存,还极大简化了硬件实现中的地址索引逻辑。
2. 对数域和积算法(Log-SPA)
这是译码部分的核心。其核心思想是在对数域进行信息的更新更新:
- 校验节点更新:使用双曲正切函数(tanh)及其逆函数处理多个输入信息的乘积,计算校验方程对特定变量节点的约束置信度。
- 变量节点更新:对来自所有相邻校验节点的信息进行简单累加,并结合原始信道观测值,更新该节点的LLR估计。
- 早期停止准则:每一轮迭代后都会进行硬判决并校验 $H cdot x = 0$。若校验通过,说明码字已成功恢复,立即停止迭代以节省资源。
3. 简化高斯消元算法
该算法专门针对GF(2)域设计,通过行交换和行异或操作,将稀疏矩阵H转换为系统形式。这保证了编码器可以直接从生成矩阵G中提取出校验位,与原始信息位拼接完成编码。
4. CDMA扩频与解扩
系统利用Hadamard矩阵的行正交性。在发送端,每个比特被扩展为16个码片(Chip);在接收端,通过与相同掩码做內积并归一化,有效地从多用户热噪声混合体中提取出目标用户的能量。
系统要求
- 软件环境:MATLAB R2016b 或更高版本。
- 工具箱需求:无需特殊工具箱,代码由纯MATLAB语言编写,具备良好的跨平台兼容性。
- 硬件说明:由于涉及多用户多帧迭代运算,建议内存4GB以上。
使用方法
- 在MATLAB的工作路径中打开项目文件夹。
- 运行主程序函数。
- 观察命令行窗口输出的实时仿真进度(包含各SNR点下的BER结果)。
- 仿真结束后,程序将自动弹出两个可视化图表:
- 图1展示AWGN与瑞利信道下系统的误码率性能对比曲线。
- 图2展示在特定SNR下,LDPC译码算法在迭代过程中的收敛情况。