【MATLAB源码-第421期】基于MATLAB的SC-FDE单载波频域均衡系统仿真,包含MPSK、链路同步、信道估计。

操作环境:

MATLAB 2024a

1 、算法描述

摘要

单载波频域均衡技术兼具单载波信号峰均比较低和频域均衡计算复杂度较低的优势,在宽带多径环境下具有较高的工程应用价值。围绕单载波频域均衡通信体制,本文建立了一个包含帧同步、定时同步、载波同步、信道估计、频域均衡和误码率评估的完整多阶MPSK通信链路仿真平台,并在MATLAB R2024a环境下完成系统实现与性能分析。系统采用统一字序列与训练块联合构帧方式,在接收端通过归一化相关完成帧起始检测,通过多支路采样峰均比判决完成定时同步,通过重复训练序列相位差实现载波频偏估计,并结合已知训练块完成残余相位补偿。在信道处理部分,本文采用最小二乘初始估计、频域平滑和MMSE收缩相结合的信道估计方法,并配合频域最小均方误差均衡器恢复数据信号。针对QPSK、8PSK、16PSK和32PSK四种调制方式,本文从载波同步前星座、同步后均衡前星座、均衡后星座以及误码率曲线四个维度进行对比分析。仿真结果表明,所构建的接收机结构能够在存在定时偏移、载波频偏、固定初始相位和多径衰落的条件下稳定完成同步与均衡,且误码率随信噪比升高呈总体下降趋势。随着调制阶数提高,系统频谱效率提升,但对同步精度和信道估计误差更加敏感,均衡后星座聚集程度下降、误码率性能劣化。本文的研究结果能够为SC-FDE体制下多阶MPSK链路的参数设计、算法选取和系统级仿真提供参考。

关键词

单载波频域均衡;MPSK;帧同步;载波同步;MMSE信道估计;MATLAB仿真

1 引言

随着高速无线通信、宽带接入和复杂传播环境下可靠传输需求的不断增长,通信系统在提高频谱利用率的同时,必须兼顾多径衰落抑制、同步稳定性和接收端实现复杂度等关键问题。正交频分复用技术在频率选择性信道中应用广泛,但其较高的峰均比会对功率放大器线性度和射频前端效率提出更高要求。相比之下,单载波频域均衡技术在保持单载波低峰均比优势的同时,可以借助快速傅里叶变换将卷积均衡转化为频域逐点运算,因此在宽带通信、海洋通信、无人机载荷通信以及低功耗终端链路中得到持续关注[1-5]。

从接收机实现角度看,SC-FDE系统的性能不仅取决于均衡器结构,还高度依赖于前端同步与信道估计质量。帧同步偏差会导致有效数据块提取位置错误,定时偏差会引起采样失真与符号间干扰加剧,载波频偏和相位偏差会破坏星座点的相位稳定性,而信道估计误差则会直接影响频域均衡权值的准确性[6-10]。因此,构建一个同步、估计、均衡一体化的完整通信链路模型,对于理解SC-FDE系统的工作机理和性能边界具有重要意义。

另一方面,多阶MPSK调制在同一符号速率下能够提供不同水平的频谱效率。低阶调制如QPSK具有较大的相位间隔,抗噪声与抗估计误差能力较强;高阶调制如16PSK和32PSK能够提升单位带宽传输比特数,但其星座点间欧氏距离减小,对残余频偏、相位噪声、均衡误差和信道估计误差更为敏感。基于这一特点,在统一接收框架下比较不同MPSK调制的星座演化过程和误码率表现,具有较强的分析价值。

基于上述背景,本文围绕"基于SC-FDE单载波频域均衡的MPSK通信链路仿真"这一主题,建立包含发射成帧、信道作用、匹配滤波、帧同步、定时同步、载波同步、LS/MMSE信道估计、频域均衡以及误码率统计的完整MATLAB仿真系统,并对QPSK、8PSK、16PSK和32PSK四种调制进行统一测试与比较。本文工作的重点不在于给出抽象理论推导,而在于面向系统实现,验证完整接收链路中各功能模块对最终性能的影响,并在同一仿真平台下揭示调制阶数、同步精度与均衡效果之间的关系。

2 系统模型与关键算法

2.1 发射端帧结构与信号模型

本文所建立的发射端采用块传输方式。首先,将随机二进制比特流按照每符号`k=log2(M)`比特进行分组,并映射为MPSK复符号序列。随后,系统构造由两段重复短统一字序列、一个带循环前缀的已知训练块和一个带循环前缀的数据块组成的发送帧。该构帧结构兼顾了不同功能需求:重复短序列用于载波频偏粗估计,训练块用于帧定位与信道估计,数据块用于实际有效信息传输。

设发送符号序列为`s(n)`,多径离散信道冲激响应为`h(l)`,载波频偏归一化值为`Δf`,初始相位为`φ0`,则接收基带离散信号可表示为:

`r(n)=exp[j(2πΔf n+φ0)] * Σ h(l)s(n-l) + w(n)`

其中,`w(n)`为加性高斯白噪声。由该模型可见,接收信号同时受到多径卷积、相位旋转与噪声扰动的共同影响。因此,接收端必须先完成同步和信道恢复,才能实现可靠判决。

在本文的仿真中,信道采用三径离散模型,主径和两条回波路径的复增益分别对应不同的幅度和相位,能够形成典型的频率选择性衰落特征。为增强可观测性,系统在发送端采用过采样构造离散波形,并在接收端进行匹配滤波和抽样支路选择,从而显式体现定时同步模块的作用。

2.2 帧同步与定时同步

帧同步的任务是从连续接收信号中准确找到一帧起始位置。若仅采用非归一化相关度量,在接收功率起伏明显或噪声较大时,相关峰值容易受到局部幅度变化影响,进而出现误判。基于此,本文采用统一字序列与接收信号之间的归一化相关度量:

`M(n)=|Σ r(n+i)u*(i)| / sqrt(Σ|r(n+i)|^2 · Σ|u(i)|^2)`

其中,`u(i)`为本地已知统一字序列。该度量能够削弱接收功率波动对相关峰定位的干扰,提升复杂环境下的帧同步鲁棒性。

考虑到系统使用`4`倍过采样,实际最佳采样时刻可能对应不同抽样相位。为此,接收端在四个采样支路上分别计算归一化相关度量,并采用"峰值与平均值之比"作为支路优选指标。其物理含义在于:若某一采样相位更接近最佳采样时刻,则该相位上的正确相关峰不仅更高,而且与背景相关底噪之间的区分更明显。因此,使用峰均比而非单纯峰值作为定时同步判据,可以在多支路比较时获得更稳定的效果。

完成支路选择后,系统记录最佳抽样相位及对应的帧起始位置,为后续载波同步、训练块提取和数据块截取提供依据。该方法结构清晰、实现直接,适合课程设计和系统仿真场景下的链路级建模。

2.3 载波同步与残余相位补偿

载波同步主要用于补偿发射端与接收端本振失配、平台运动或其他因素引起的频率偏差和相位旋转。本文的帧头中包含两段完全相同的短统一字序列,因此可利用它们的相位差实现频偏估计。设两段序列分别为`a(n)`和`b(n)`,则归一化频偏估计值可写为:

`Δf_hat = angle[Σ a*(n)b(n)] / (2πL0)`

其中,`L0`为短统一字序列长度。该方法利用重复结构直接估计相位增量,具有实现简单、开销较小的特点,适合作为载波频偏粗估计。

在得到`Δf_hat`后,接收端通过复指数乘法对整帧信号进行频偏校正。然而,粗频偏校正后通常仍会残留少量公共相位偏差,特别是在噪声、多径和信道估计误差共同存在时更为明显。为提高判决精度,本文进一步利用已知训练块在均衡后的输出与本地训练序列之间的相位差估计残余相位,并在数据块恢复阶段统一补偿。该两级处理结构能够显著改善均衡后星座点的聚集效果,使后续MPSK解调更稳定。

2.4 LS/MMSE信道估计与频域均衡

由于SC-FDE采用块传输方式,若循环前缀长度不小于信道有效长度,则数据块在去除循环前缀后可近似看作循环卷积,从而能够在频域内实现逐子载波均衡。本文首先在训练块上进行快速傅里叶变换,设训练块频域发送序列为`X(k)`,接收序列为`Y(k)`,则最小二乘信道估计为:

`H_LS(k)=Y(k)/X(k)`

LS估计实现简单,但在中低信噪比下容易受到噪声扰动。为减弱估计抖动,本文对LS结果进行频域滑窗平滑,得到较为稳定的平滑信道谱。进一步地,利用估计误差功率构造MMSE型收缩因子,对LS结果进行幅度压缩,以降低噪声主导频点对均衡器权值的放大作用。该策略兼顾了估计准确性与计算复杂度,适合课程设计级别的系统仿真实现。

在此基础上,本文构造频域最小均方误差均衡器:

`G(k)=H*(k)/( |H(k)|^2 + λ )`

其中,`H(k)`为经平滑和MMSE修正后的信道估计,`λ`为等效噪声功率项。对接收数据块进行FFT后,与均衡器权值逐点相乘,再通过IFFT变换回时域即可恢复数据符号。最后,对均衡后的符号执行残余相位补偿和MPSK硬判决,完成比特恢复。

上述处理流程表明,SC-FDE的核心思想并非回避多径影响,而是利用循环前缀将时域卷积问题转移到频域,在同步与信道估计的支撑下,以较低复杂度完成块级均衡。这也是本文系统设计的理论基础。

3 仿真平台与参数设置

3.1 仿真平台

本文采用MATLAB R2024a作为仿真平台,建立发射机、信道与接收机一体化的系统级模型。之所以选择MATLAB,一方面是因为其复数信号处理、向量化运算和频域分析功能完善,便于快速搭建SC-FDE链路;另一方面,MATLAB具备较成熟的绘图与数据处理能力,适合输出星座图、误码率曲线等结果用于对比分析。

3.2 关键参数

系统主要参数设置如下:数据块长度`N=256`,循环前缀长度`CP=32`,过采样倍数`sps=4`,短统一字序列长度`L0=64`。接收端载波频偏设为`0.003`,固定初始相位为`0.35 rad`,帧前加入长度为`160`的定时偏移。信道采用三径离散模型,其复增益分别对应主径、次径和第三径,可形成较明显的频率选择性衰落。误码率曲线的扫描范围为`0~30 dB`,步长为`2 dB`,每个信噪比点采用`3000`次蒙特卡洛统计,以保证结果具有一定稳定性。

在调制方式方面,系统分别测试QPSK、8PSK、16PSK和32PSK四种MPSK体制。不同调制方式共享相同的同步和均衡结构,仅在每符号比特数和相位映射方式上存在差异。该统一框架能够较好地反映调制阶数变化对接收机性能的影响。

3.3 程序流程

整个仿真流程可以概括为以下几个步骤。首先,发射端生成随机比特并完成MPSK映射,然后插入统一字序列、训练块和数据块,形成完整发送帧。其次,信号经过多径信道、频偏旋转、固定相位偏移和加性白噪声后,进入接收端。接收端先执行匹配滤波与多支路抽样,再通过归一化相关完成定时支路选择与帧同步。随后,系统利用重复短统一字序列估计载波频偏并进行校正,从训练块获得LS/MMSE信道估计,构建频域均衡器,对数据块进行恢复。最后,通过MPSK判决输出比特并统计误码率,同时保存不同处理阶段的星座分布用于可视化分析。

为全面评估系统行为,本文对每种调制方式均输出四类结果:载波同步前星座图、同步后均衡前星座图、均衡后星座图以及误码率曲线。这样的结果组织方式能够直观体现同步、信道估计和频域均衡对信号恢复过程的逐级作用。

4 仿真结果与性能分析

4.1 星座图演化分析

从载波同步前星座图可以看到,接收符号在载波频偏和相位旋转影响下呈现明显的环状扩散分布。对于QPSK和8PSK,这种现象尤其直观:理想星座点位于固定相位位置,而实际接收点由于相位随时间不断漂移,被拉伸为近似圆环或厚环结构。随着调制阶数提高,理想相位点数量增多,环状扩散对相邻相位簇的覆盖也更加明显。

完成载波同步后,频偏引起的整体旋转得到抑制,星座图由连续环状逐渐向若干模糊的簇状结构收缩。但此时由于多径引起的幅相失真尚未完全消除,星座点仍围绕理想位置呈现一定范围的离散扩散,特别是在8PSK、16PSK和32PSK条件下,相邻簇之间的分界并不十分清晰。这一现象说明仅靠同步处理并不足以恢复信号,还必须依赖后续的信道估计和均衡过程。

在频域均衡和残余相位补偿之后,星座聚集效果显著增强。QPSK均衡后四个星座点能够较稳定地收敛到理想位置附近,点云分布紧凑,表明在当前信道与噪声条件下系统能够获得较好的恢复效果。8PSK均衡后八个簇也能够清晰分开,虽然点云半径和角向扩散仍略大于QPSK,但已具备较好的判决边界。进入16PSK和32PSK后,均衡后星座仍然沿理想圆周附近分布,但点簇之间的相位间隔明显缩小,因此残余估计误差更容易引起判决翻转。该结果与高阶MPSK调制对相位精度更敏感的理论认识是一致的。

4.2 误码率曲线分析

误码率曲线表明,四种调制方式在`0~30 dB`信噪比范围内均呈现随信噪比提高而总体下降的趋势,这说明所建立的同步、信道估计与均衡链路能够在不同噪声水平下保持基本稳定。QPSK曲线下降最快,在中高信噪比区间内表现出最优的可靠性;8PSK次之;16PSK和32PSK的误码率明显高于前两者,且高信噪比下的下降速度相对缓慢。

从仿真结果的量级看,在`30 dB`附近,QPSK误码率已进入`10^-4`数量级,8PSK约处于`10^-3`数量级,16PSK约为`10^-2`数量级,而32PSK仍维持在`10^-2`至`10^-1`之间。造成这一差异的主要原因在于:随着调制阶数增加,相邻相位点之间的角距离减小,同样大小的残余相位误差、频偏估计偏差、信道估计误差或噪声扰动都更容易跨越判决边界,进而引发误码。

值得注意的是,高阶MPSK系统的误码性能不仅受噪声主导,也与同步精度紧密相关。对于32PSK而言,即使均衡后星座整体仍沿理想圆周排列,只要存在较小的公共相位误差或局部相位抖动,相邻点之间就可能发生混淆。因此,在高阶相位调制条件下,仅提升信噪比并不能完全弥补同步和估计误差,接收机往往需要更高精度的频偏跟踪、相位恢复或判决反馈机制。

4.3 同步与信道估计模块的作用分析

本文仿真结果还表明,前端同步模块对后续均衡性能具有基础性影响。若帧同步位置出现偏差,则训练块和数据块的截取边界将发生错位,频域均衡器权值也随之失效;若定时支路选择不正确,即便帧起始位置正确,接收信号仍可能因采样失真而导致星座扩散加剧。因此,采用归一化相关与峰均比联合判决的同步策略,是保证链路稳定工作的前提。

在信道估计方面,单纯使用LS估计虽然能够给出快速的频域信道响应,但在中低信噪比条件下估计曲线容易出现明显波动,直接用于构造均衡器时,部分频点会对噪声过度放大。本文通过频域平滑和MMSE型收缩对LS结果进行修正,使信道估计谱更加平滑,均衡器权值更稳定,从而改善了均衡后星座的集中程度和误码率表现。这说明在课程设计规模的仿真平台中,将"LS初估 + 平滑 + MMSE修正"作为组合策略,能够在计算量和性能之间取得较为合理的平衡。

4.4 系统特点与局限性

从整体实现效果看,本文系统具有以下特点。其一,发射端构帧结构清晰,统一字序列、训练块和数据块职责明确,便于实现链路分段处理。其二,接收端算法之间衔接紧密,帧同步、定时同步、载波同步、信道估计和均衡各模块构成完整闭环。其三,在统一接收框架下比较多阶MPSK调制,能够较直观地反映调制阶数变化带来的性能代价。

与此同时,本文模型也存在一定局限。首先,所用信道为块内近似静止的三径模型,尚未引入快时变信道和更复杂的多普勒扩展;其次,载波同步采用的是重复训练序列粗估计加训练块残余相位补偿的结构,对于更复杂的相位噪声场景仍有改进空间;再次,系统当前以硬判决误码率为评价指标,尚未引入信道编码、软判决译码或迭代均衡机制。因此,后续工作可以考虑将SC-FDE与信道编码、判决反馈、深度学习辅助估计或自适应跟踪方法结合,以进一步提升高阶MPSK条件下的鲁棒性[1][4][5][9]。

5 结论

本文围绕SC-FDE单载波频域均衡多阶MPSK通信链路,完成了包含帧同步、定时同步、载波同步、LS/MMSE信道估计、频域均衡和误码率评估的完整MATLAB系统仿真。仿真结果表明,采用"归一化相关帧同步 + 多支路峰均比定时选择 + 重复序列频偏估计 + 训练块残余相位补偿 + LS/MMSE信道估计 + MMSE频域均衡"的接收结构,能够在存在多径、频偏、相位偏差和噪声的条件下有效恢复发送数据。QPSK在可靠性方面表现最佳,8PSK在性能与频谱效率之间取得折中,16PSK和32PSK则体现出更高频谱效率下对同步与估计精度的更高要求。总体而言,本文构建的仿真平台较为完整地反映了SC-FDE系统中各关键模块的协同关系,可为后续开展更复杂信道模型、迭代接收机设计以及高阶调制鲁棒性研究提供基础。

参考文献

1\] Yun S M, Ryu Y B, Song D H, et al. Narrowband Interference Cancellation Using a Fine Frequency Shift in Single-Carrier Frequency Domain Equalization (SC-FDE) Systems\[J\]. Electronics, 2025, 14(4): 798. DOI: 10.3390/electronics14040798. \[2\] Marey M, Mostafa H. Estimation of IQI for AF Cooperative Single-Carrier Frequency Domain Equalization Systems Using Channel Decoder Feedback\[J\]. Electronics, 2023, 12(4): 863. DOI: 10.3390/electronics12040863. \[3\] Han R, Jia N, Guo Z, et al. Prefiltered Single-Carrier Frequency-Domain Equalization for Binary CPM over Shallow Water Acoustic Channel\[J\]. Sensors, 2022, 22(10): 3821. DOI: 10.3390/s22103821. \[4\] Lin Z, Wang L, Peng C, et al. Direct Adaptive Multi-Resampling Turbo Equalizer for Underwater Acoustic Single-Carrier Communication\[J\]. Journal of Marine Science and Engineering, 2024, 12(8): 1271. DOI: 10.3390/jmse12081271. \[5\] Zhao H, Yao K, Xiang D, et al. Learning to Equalize for Single-Carrier Underwater Acoustic Communications\[J\]. Journal of Marine Science and Engineering, 2025, 13(11): 2209. DOI: 10.3390/jmse13112209. \[6\] Qu F, Li Z, Zhang M, et al. A Carrier-Based Gardner Timing Synchronization Algorithm for BPSK Signal in Maritime Communication\[J\]. Journal of Marine Science and Engineering, 2023, 11(4): 829. DOI: 10.3390/jmse11040829. \[7\] Zhang S, Wang L, Liu K, et al. Low Complexity Parallel Carrier Frequency Offset Estimation Based on Time-Tagged QPSK Partitioning for Coherent Free-Space Optical Communication\[J\]. Photonics, 2024, 11(9): 885. DOI: 10.3390/photonics11090885. \[8\] Liu Z, Gong K, Sun P, et al. OQPSK Synchronization Parameter Estimation Based on Burst Signal Detection\[J\]. Electronics, 2021, 10(1): 69. DOI: 10.3390/electronics10010069. \[9\] Tarafder P, Chun C, Ullah A, et al. Channel Estimation in 5G-and-Beyond Wireless Communication: A Comprehensive Survey\[J\]. Electronics, 2025, 14(4): 750. DOI: 10.3390/electronics14040750. \[10\] Martins J, Conceição F, Gomes M, et al. Joint Channel Estimation and Synchronization Techniques for Time-Interleaved Block-Windowed Burst OFDM\[J\]. Applied Sciences, 2021, 11(10): 4403. DOI: 10.3390/app11104403. ## ****2** **、仿真结果演示**** ![](https://i-blog.csdnimg.cn/direct/ab3a42a0f2364d97abe9b7d0281ff14c.png)![](https://i-blog.csdnimg.cn/direct/6cc553805e9746d28d19eb0575061d26.png)![](https://i-blog.csdnimg.cn/direct/4a1908e393bd4b519e43742f2ea6314e.png)![](https://i-blog.csdnimg.cn/direct/acfc8e39f62948edac41cad1a8903f63.png) ## ****3** **、关键代码展示**** 略 ## ****4** **、MATLAB 源码获取**** V 点击下方名片关注公众号获取