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

【MATLAB源码-第442期】基于MATLAB的OFDM系统PAPR抑制算法仿真及限幅压扩SLM、PTS与TR性能对比

news 2026/5/21 18:56:26

操作环境：

MATLAB 2024a

1、算法描述

摘要
正交频分复用技术能够把高速数据流分解到多个正交子载波上传输，因此在宽带通信系统中具有较高的频谱利用率和较强的抗频率选择性衰落能力。公开资料显示，OFDM 已经用于 DAB、DVB、WLAN、WiMAX、第四代和第五代移动通信等系统，5G NR 物理层标准也以 OFDM 作为重要波形基础。但是，OFDM 多个子载波在时域叠加后容易产生较高峰均功率比。较高 PAPR 会提高功率放大器线性工作要求，降低功放效率，并可能带来带内失真和邻道泄漏问题。因此，PAPR 抑制是 OFDM 发射链路设计中的重要问题。本文基于 MATLAB 对原始 OFDM、限幅法、压扩法、SLM 选择映射、PTS 部分传输序列和 TR 保留音法进行建模仿真。系统采用 16QAM 调制、过采样 IFFT、循环前缀和 AWGN 信道，分别从 PAPR 互补累积分布、误码率、星座图、平均 PAPR、算法复杂度和功率谱六个角度展开对比。仿真结果能够直观看出不同 PAPR 抑制方法在峰值控制能力、信号失真、误码率影响和计算复杂度方面的差异。整体来看，限幅和压扩方法结构简单，但更容易引入失真；SLM 和 PTS 具有较好的无失真抑峰能力，但需要边信息和较高计算量；TR 方法通过保留子载波构造抵消信号，能够在不直接削弱数据子载波的情况下抑制峰值，但会牺牲部分频域资源。本文仿真模型结构清晰，指标完整，适合用于 OFDM PAPR 抑制算法的学习、验证与横向比较。

关键词： OFDM；峰均功率比；PAPR 抑制；16QAM；限幅法；压扩法；SLM；PTS；TR 保留音法；误码率分析

1. 引言
OFDM 的核心思想是把一个高速传输任务分配到多个低速正交子载波上完成。每个子载波之间保持频域正交关系，因此子载波频谱可以相互重叠而不产生传统意义上的严重互扰。这样做可以提高频谱利用率，也可以把复杂的宽带频率选择性信道转化为多个相对简单的窄带子信道。对于宽带无线通信、移动通信和多载波传输系统而言，这种结构具有明显工程价值。

但是，OFDM 的优势和问题是绑定在一起的。由于发送端需要把多个子载波信号叠加到同一个时域波形中，当多个子载波相位在某些采样点附近同向叠加时，时域瞬时幅度会显著升高。这个现象表现为峰均功率比增大。PAPR 越高，发射机功率放大器越难在高效率状态下保持线性输出。如果功放为了避免削顶而预留很大功率回退，系统能量效率会下降；如果功放没有足够回退，输出波形会出现非线性压缩，导致误码率变差和频谱扩展。近年的 PAPR 研究仍然持续关注限幅滤波、SLM、PTS、TR 等方法，说明该问题在 OFDM 及其扩展波形中仍有现实意义。

本文围绕 OFDM PAPR 抑制问题建立 MATLAB 仿真系统。系统不是只展示单一算法，而是把原始 OFDM、限幅、压扩、SLM、PTS 和 TR 放在同一参数条件下对比。这样做的价值在于可以直接观察不同方案的取舍关系。单纯看 PAPR 曲线并不够，因为有些方法降低了峰值，却损伤了数据符号；单纯看 BER 也不够，因为有些方法误码率较好，但计算复杂度过高或需要额外边信息。因此，本文采用多指标评价方式，使仿真结果更接近真实通信系统中的工程判断逻辑。

2. 系统模型
本文仿真系统采用典型 OFDM 基带链路。发送端首先产生随机二进制比特流，然后通过 16QAM 调制形成复数星座符号。调制符号被映射到有效数据子载波位置，并避开直流子载波。为了提高 PAPR 统计的准确性，系统采用过采样方式扩展 IFFT 点数，使时域峰值观察更加充分。随后，不同 PAPR 抑制算法分别作用在 OFDM 频域或时域信号上。最后，系统添加循环前缀，并通过 AWGN 信道传输。

接收端首先去除循环前缀，然后通过 FFT 恢复频域子载波。对于 SLM 和 PTS 等需要相位补偿的方法，接收端利用发送端保存的相位边信息对数据符号进行恢复。对于压扩方法，接收端执行对应的反压扩处理。完成频域恢复后，系统提取数据子载波并进行 16QAM 解调，将接收比特与原始发送比特进行比较，从而计算误码率。

该模型包含两个核心仿真流程。第一个流程用于 PAPR 统计，系统在大量随机 OFDM 帧上计算每种算法的 PAPR，并进一步绘制 CCDF 曲线、平均 PAPR 柱状图和复杂度对比图。第二个流程用于误码率分析，系统在多个 SNR 点上重复发送随机帧，统计不同算法的 BER，并在指定 SNR 下保存接收星座点。这样的设计能够同时覆盖峰值统计、传输可靠性和星座失真情况。

3. PAPR 抑制算法原理
原始 OFDM 不进行额外峰值处理，直接把频域符号通过 IFFT 转换为时域信号。它作为基准方案，用于观察未抑制时的 PAPR 水平和误码性能。原始方案结构最简单，不需要边信息，也不会引入额外失真，但其峰值较高，对功放线性度要求最大。

限幅法是最直接的 PAPR 抑制方法。它设置一个幅度门限，当时域信号幅度超过门限时，直接把超出部分压到门限附近。限幅法实现简单，计算量低，也不需要复杂边信息。但是限幅本质上属于非线性削峰，会直接改变原始 OFDM 时域波形。这样虽然能降低峰值，却会引入带内失真和频谱扩展，因此误码率和功率谱可能受到影响。

压扩法通过非线性幅度变换压缩大幅度样点，同时尽量保留小幅度样点的可辨识度。相比硬限幅，压扩通常更平滑，削峰过程没有限幅那么突兀。系统中采用 mu 律压扩思想，并在接收端执行反压扩。该方法能够降低峰值，但仍然属于非线性处理，因此在噪声环境下反压扩可能放大部分扰动，误码性能不一定总是优于原始 OFDM。

SLM 选择映射属于概率类无失真 PAPR 抑制方法。它通过多组相位扰动序列生成多个候选 OFDM 符号，然后从中选择 PAPR 最低的候选符号发送。SLM 的优点是不会像限幅那样直接削弱时域峰值，因此理论上对数据星座的破坏更小。缺点是发送端需要多次 IFFT 或等效候选计算，同时接收端需要知道所选相位序列，否则无法正确恢复原始符号。

PTS 部分传输序列方法把数据子载波划分成多个子块，每个子块分别变换到时域后乘以不同相位因子，再通过搜索选择合适的相位组合，使合成后的 OFDM 符号 PAPR 降低。PTS 通常具有较好的抑峰效果，但搜索组合会随着子块数和相位集合扩大而迅速增加。已有综述研究也把 PTS 作为 OFDM PAPR 抑制中的重要方向进行系统讨论。

TR 保留音法在频域中预留一部分子载波，不承载用户数据，而是用于生成峰值抵消信号。系统通过多次迭代检测时域峰值，并把抵消分量投影到保留音位置。TR 的优点是数据子载波不被直接削顶，缺点是需要牺牲部分频域资源。近年仍有研究使用范数优化、流形优化等方法改进 TR 的峰值控制能力和资源利用效率。

4. 仿真设计
本文 MATLAB 工程采用主函数与函数文件分离结构。主函数 main.m 负责清理环境、添加路径、创建图像文件夹、读取参数、调用仿真函数和输出结果。所有功能函数放在 function 文件夹中，包括参数配置、调制解调、OFDM 发送接收、PAPR 计算、不同算法实现和绘图函数。

参数设置方面，系统采用 128 个基础子载波，并设置 4 倍过采样。有效子载波避开直流位置，其中一部分子载波被设置为 TR 保留音，其余子载波用于承载 16QAM 数据。PAPR 统计帧数设置为 1500 帧，误码率统计帧数设置为每个 SNR 点 180 帧，SNR 范围为 0 dB 到 18 dB。SLM 使用 8 组候选相位序列，PTS 使用 4 个子块和四相位集合，TR 设置 6 次迭代。以上参数能够在仿真耗时和结果稳定性之间取得相对平衡。

绘图输出包括六类结果。第一类是 PAPR CCDF 曲线，用于观察不同阈值下 PAPR 超过该阈值的概率。第二类是 BER 曲线，用于评价算法在噪声信道下的可靠性。第三类是星座图，用于观察非线性处理或噪声对接收符号分布的影响。第四类是复杂度对比图，用理论运算量和实测单帧耗时共同反映算法代价。第五类是平均 PAPR 柱状图，用于展示每种方法的平均削峰能力和离散程度。第六类是功率谱图，用于观察不同方法是否对频谱形态产生明显影响。

5. 仿真结果分析
从 PAPR CCDF 曲线来看，原始 OFDM 通常具有最高的 PAPR 分布尾部。这符合 OFDM 多子载波时域叠加的基本特性。限幅法、压扩法、SLM、PTS 和 TR 都能够不同程度降低 PAPR 超限概率。CCDF 曲线越靠左，表示在相同概率下所需承受的 PAPR 阈值越低，也说明抑峰效果越明显。

从平均 PAPR 柱状图来看，各算法的整体排序能够进一步验证 CCDF 曲线的趋势。平均 PAPR 不是唯一评价指标，但它可以快速反映不同算法在大量随机帧上的平均削峰能力。如果某种方法平均值较低且标准差较小，说明该方法不仅有较强抑峰效果，而且对不同随机帧具有更稳定的表现。

从 BER 曲线来看，不同算法的传输可靠性差异不能简单等同于 PAPR 抑制能力。限幅和压扩可能降低峰值，但它们对时域波形进行了非线性处理，因此在接收端可能出现星座点扩散和判决误差增加。SLM 和 PTS 主要通过相位候选选择降低峰值，理论上不直接破坏数据符号，因此在边信息正确恢复的条件下，误码率一般更容易保持接近原始 OFDM。TR 方法不直接修改数据子载波，但保留音的抵消信号会改变整体时域波形，同时牺牲部分频域资源，因此其 BER 和 PAPR 表现需要结合参数共同判断。

从星座图来看，星座点越集中，说明接收端恢复质量越好。若某种方法的星座点明显发散，说明该方法可能引入了更强失真或噪声放大。限幅类方法在高削峰强度下容易造成星座扩散；压扩类方法在反压扩后也可能对噪声敏感。SLM 和 PTS 如果边信息准确，星座结构通常更清晰。TR 方法的星座质量与保留音数量、迭代次数和抵消步长有关。

从复杂度图来看，原始 OFDM、限幅和压扩的结构较简单，适合低复杂度场景。SLM 需要生成多个候选符号，其计算量与候选序列数量直接相关。PTS 需要搜索多个相位组合，复杂度随子块数和候选相位数量增长。TR 方法需要多次迭代检测峰值并更新保留音，因此计算代价也高于简单限幅和压扩。工程应用中不能只追求最低 PAPR，还必须考虑硬件资源、实时性、边信息开销和频谱约束。

从功率谱图来看，限幅等非线性处理可能带来频谱形态变化。虽然本仿真主要是基带算法对比，但功率谱观察仍然有价值。实际发射机中，频谱扩展会关系到邻道泄漏和射频滤波压力。TR、SLM 和 PTS 等方法通常比硬削峰更利于控制频谱副作用，但具体表现仍然依赖参数设置。

6. 工程意义
该仿真工程的价值不在于把某一种算法包装成绝对最优，而在于建立一个公平的横向比较框架。PAPR 抑制本身就是典型的多目标问题。一个方法可能削峰很强，但误码率变差；另一个方法可能误码率保持较好，但计算复杂度高；还有一些方法需要额外边信息或牺牲子载波资源。工程设计必须承认这些代价，而不是只看单一曲线。

对于实际通信系统，PAPR 抑制方案需要结合功率放大器线性区、发射功率目标、带宽配置、调制阶数、接收端边信息恢复能力和硬件算力共同选择。移动终端更关注功放效率和实现复杂度，基站侧可能有更强算力但也需要满足频谱模板要求。若系统更看重简单实现，可以优先考虑限幅和压扩；若系统更看重信号质量，可以考虑 SLM 或 PTS；若系统允许预留部分频域资源，可以考虑 TR 及其优化变体。

7. 系统特点
本文系统具有较完整的算法覆盖范围。工程不仅包含原始 OFDM，还同时加入限幅法、压扩法、SLM、PTS 和 TR 六类方案，能够覆盖失真类、概率选择类、相位优化类和保留音类 PAPR 抑制思路。系统输出结果也比较完整，既有 PAPR CCDF 曲线，也有 BER、星座图、复杂度、平均 PAPR 和频谱图，能够从多个角度支撑结论。

本文系统结构清楚。主函数只负责整体流程调度，具体算法被封装到独立函数中。这样的结构便于后续扩展，例如替换调制方式、增加信道模型、改变子载波数量、增加功放非线性模型或加入新的 PAPR 抑制算法。代码中的关键位置已经加入中文注释，便于理解每个函数的作用和主要处理逻辑。

本文系统的评价方式较接近工程判断。它没有只关注 PAPR 曲线，而是把 BER 和复杂度一起纳入分析。这样可以避免“PAPR 降得越低越好”的片面判断。对于真实系统而言，通信质量、实现成本、实时性和频谱约束必须一起考虑。该工程能够帮助使用者理解不同算法之间的真实取舍关系。

8. 结论
本文基于 MATLAB 建立了 OFDM 系统 PAPR 抑制算法仿真模型，对原始 OFDM、限幅法、压扩法、SLM、PTS 和 TR 保留音法进行了统一参数下的对比分析。仿真系统采用 16QAM 调制、过采样 IFFT、循环前缀和 AWGN 信道，并从 PAPR CCDF、误码率、星座图、复杂度、平均 PAPR 和功率谱六个方面输出结果。

分析表明，原始 OFDM 结构简单但 PAPR 较高。限幅和压扩方法实现简单，能够快速降低峰值，但可能引入非线性失真。SLM 和 PTS 能够通过候选相位选择降低 PAPR，信号质量保持较好，但需要额外计算量和边信息支持。TR 方法利用保留音生成峰值抵消信号，能够避免直接削顶数据子载波，但需要牺牲部分频域资源，并且性能依赖迭代参数和保留音位置。

整体而言，本文仿真结果说明，OFDM PAPR 抑制没有单一绝对最优方案。不同方法适合不同工程条件。若系统关注低复杂度，可以选择限幅或压扩；若系统关注较低失真，可以选择 SLM 或 PTS；若系统允许保留部分子载波资源，可以考虑 TR 方法。该工程能够为 OFDM PAPR 抑制算法学习、性能验证和方案比较提供较完整的 MATLAB 仿真基础。

参考文献
[1] ETSI. ETSI TS 138 211 V18.2.0: 5G; NR; Physical Channels and Modulation (3GPP TS 38.211 Version 18.2.0 Release 18). ETSI, 2024.

[2] Kiambi, Stephen, Elijah Mwangi, and George Kamucha. “Reducing PAPR of OFDM Signals Using a Tone Reservation Method Based on l∞-Norm Minimization.” Journal of Electrical Systems and Information Technology, vol. 9, article 12, 2022.

[3] Al-Jawhar, Yasir Amer, et al. “A Review of Partial Transmit Sequence for PAPR Reduction in the OFDM Systems.” IEEE Access, vol. 7, 2019, pp. 18021–18041.

[4] Gökceli, Selahattin, et al. “Frequency-Selective PAPR Reduction for OFDM.” arXiv, 2019.

[5] Gökceli, Selahattin, et al. “PAPR Reduction with Mixed-Numerology OFDM.” arXiv, 2019.

[6] Liu, Maolin, Wei Chen, Jialong Xu, and Bo Ai. “A Comprehensive Study of PAPR Reduction Techniques for Deep Joint Source Channel Coding in OFDM Systems.” arXiv, 2023.

[7] Rexhepi, Getuar, et al. “Tone Reservation-Based PAPR Reduction Using Manifold Optimization for OFDM-ISAC Systems.” arXiv, 2024.

[8] Rexhepi, Getuar, Kengo Ando, and Giuseppe Thadeu Freitas de Abreu. “PAPR-Optimized OFDM Design for Opportunistic Communications and Sensing.” arXiv, 2024.

[9] MathWorks. “OFDM Modulation Using MATLAB.” MathWorks Documentation, 2024.

[10] MathWorks. “Downlink OFDM Modulation.” 5G Toolbox Documentation, 2024.