当前位置：首页 > news >正文

▲基于OFDM+QPSK的通信链路matlab性能仿真,包含LDPC,Schmidl-Cox频偏估计和MMSE信道估计

news 2026/6/9 3:42:35

✅1.本系统整体构架

✨2.各个模块基本原理

2.1 LDPC信道编码

2.2 QPSK调制

2.3 OFDM调制

2.4 Schmidl-Cox频偏与符号同步算法

2.5 OFDM解调制

2.6 MMSE最小均方误差信道估计

2.7 LDPC译码

💡3.仿真结果

👇4.完整程序下载

✅1.本系统整体构架

整个程序，我们采用如下的流程图实现：

该通信链路分为发送端、无线信道、接收端三大部分。

发送端信号流向为：原始二进制比特流→LDPC信道编码→比特交织→QPSK星座调制→组帧插入训练序列→OFDM 调制→射频发射；

信号经过多径衰落信道、载波频偏、高斯白噪声叠加后到达接收端；

接收端流程为：射频接收→下变频→Schmidl-Cox算法完成符号同步与载波频偏估计补偿→OFDM 解调制→MMSE信道估计与信道均衡→QPSK解调→解交织→LDPC译码→恢复原始比特。

整套方案结合了LDPC强纠错能力、QPSK频谱效率与抗干扰特性、OFDM抗多径衰落优势，同时依靠同步与信道估计算法解决无线传输中的频偏、信道畸变问题，是宽带无线通信的经典组合方案。

✨2.各个模块基本原理

2.1 LDPC信道编码

LDPC是一类稀疏校验矩阵的线性分组码，凭借逼近香农极限的纠错性能，成为高速通信主流编码方案。设原始信息比特序列为u=[u0,u1,…,uk−1]，码长为n，信息位长度为k，校验位长度为 m=n−k，LDPC码的核心是奇偶校验方程。定义稀疏校验矩阵H，维度为m×n，矩阵中“1”的数量远少于“0”，满足线性分组码约束关系：

其中c=[c0,c1,…,cn−1]为编码后的码字序列。LDPC编码通过生成矩阵G完成映射，码字满足c=uG。工程实现中常采用系统码结构，将码字分为信息位与校验位两部分，仅对校验位进行迭代求解。编码完成后，连续的码字比特流送入后续调制模块，LDPC通过冗余校验比特对抗信道噪声与干扰，大幅降低误码率。

2.2 QPSK调制

QPSK是正交相移调制，将连续两个二进制比特映射为一个复数星座符号，具备频谱利用率高、抗衰落能力强的特点。将LDPC输出的串行比特流按每2比特分组，记分组比特为(b1,b2)，采用格雷映射规则完成星座映射。

QPSK归一化星座点对应的复数表达式为：

该公式将二进制比特转换为复基带符号，四个星座点均匀分布在单位圆上。所有QPSK复符号组成数据符号流，完成基带调制后进入OFDM组帧环节。

2.3 OFDM调制

OFDM核心思想是将高速串行数据流拆分为多路低速子载波数据流，利用子载波正交性对抗多径衰落，同时插入循环前缀（CP）消除符号间干扰。设单帧OFDM系统子载波总数为 N，有效数据子载波数为Nd，QPSK调制后的复符号映射到对应子载波，得到频域序列X=[X0,X1,…,XN−1]。

完成CP添加后，在每帧数据起始位置插入Schmidl-Cox 训练序列（同步前导），用于接收端同步与频偏估计，最终组帧完成的基带信号经上变频后送入无线信道。

2.4 Schmidl-Cox频偏与符号同步算法

Schmidl-Cox是OFDM系统经典的联合符号同步与载波频偏估计算法，依托帧内重复结构的训练序列实现同步，分为符号定时同步和载波频偏估计两部分。训练序列设计为前后两段完全相同的时域序列，设接收端滑动窗口内前半段采样为r(n)，后半段采样为:

当滑动窗口对准训练序列位置时，M(d)出现尖锐峰值，以此确定OFDM符号起始位置，完成符号定时同步。在同步完成后，利用训练序列的相位差估计载波频偏。由频偏带来的相位旋转满足 Δϕ=πε，对互相关项P(d)取相位，可得频偏估计值：

2.5 OFDM解调制

频偏补偿与符号同步完成后，首先剔除接收信号中的循环前缀，保留长度为N的有效时域OFDM符号r^(n)。通过快速傅里叶变换将时域信号转换回频域，完成OFDM解调：

2.6 MMSE最小均方误差信道估计

MMSE信道估计是在最小化均方误差准则下求解信道频域响应，性能优于LS估计，可有效抑制噪声干扰。系统利用帧内导频符号作为已知参考，设导频位置发射符号为Xp(k)，导频位置接收频域符号为Yp(k)。

2.7 LDPC译码

最后进入LDPC译码环节，主流采用置信传播算法（BP），基于稀疏校验矩阵进行迭代译码。利用接收比特的对数似然信息，在变量节点与校验节点之间迭代传递置信度，不断修正比特判决结果，迭代收敛后输出最终译码比特u^，完成整个通信链路的数据恢复。

💡3.仿真结果

仿真参数如下：

MMSE均衡 + LDPC译码 %系统参数配置 NFFT = 64; % FFT 点数(子载波总数) CPlen = 16; % 循环前缀长度 bitsPerSym = 2; % QPSK 每符号 2 bit Fs = 20e6; % 采样率 20 MHz df = Fs / NFFT; % 子载波间隔 = 312.5 kHz %802.11a dataIdx_log = [-26:-22, -20:-8, -6:-1, 1:6, 8:20, 22:26]; % 48 数据 pilotIdx_log = [-21 -7 7 21]; % 4 导频 allIdx_log = sort([dataIdx_log pilotIdx_log]); % 52 有效 %逻辑频点 log2bin = @(x) mod(x, NFFT) + 1; dataBins = log2bin(dataIdx_log); pilotBins = log2bin(pilotIdx_log); activeBins = log2bin(allIdx_log); nData = numel(dataBins); nPilot = numel(pilotBins); nActive = numel(activeBins); pilotVal = ones(nPilot,1); %LDPC编/译码器 %DVB-S.2 LDPC, 码率1/2 H = dvbs2ldpc(1/2); cfgEnc = ldpcEncoderConfig(H); cfgDec = ldpcDecoderConfig(H); K = cfgEnc.NumInformationBits; % 32400 Ncode = cfgEnc.BlockLength; % 64800 R = K / Ncode; % 1/2 maxIter= 30; % LDPC 最大迭代 numOFDM= Ncode / bitsPerSym / nData;% 675 个数据 OFDM 符号 assert(mod(Ncode/bitsPerSym, nData)==0,'码长需与子载波数匹配'); %信道与仿真参数 EbN0_dB = 0:1:15; % Eb/N0 扫描范围 numFrames = 500; % 每点帧数(码字数) cfoTrue = 1.05; % 归一化载波频偏 %多径瑞利信道(时域抽头), 指数功率延迟谱 --- chDelays = [0 1 3]; % 抽头延迟(样点) chPdB = [0 -8 -12]; % 各抽头平均功率(dB) chPow = 10.^(chPdB/10); chPow = chPow/sum(chPow); Lch = max(chDelays)+1; % 信道长度

仿真结果如下图所示：