峰均比的计算与仿真

峰均比（Peak-to-Average Power Ratio, PAPR）是通信系统中一个重要的性能指标，特别是在多载波系统（如OFDM）中。PAPR定义为信号的最大峰值功率与平均功率的比值，通常用分贝（dB）表示。

一、峰均比的计算方法

1. 基本定义与公式

峰均比定义为峰值功率与平均功率的比值：

PAPR=max{∣x(t)∣²}/E{∣x(t)∣²}

其中：

x(t) 为信号波形

分子是信号瞬时功率的最大值

分母是信号平均功率（期望值）

电压域的峰均比（波峰因素 Crest Factor, CF）：

功率域PAPR与电压域CF的关系为：PAPR=CF²

2. 常见信号的峰均比计算

(1) 确定性周期信号

正弦波：CF=2≈1.414，PAPR=2（即3 dB）

方波：取决于占空比，50%占空比时与正弦波相同（PAPR=2）；25%占空比时PAPR=4（即6 dB）

恒包络调制信号：如理想FSK，PAPR = 0 dB

QAM信号：如16QAM

(2) 随机信号（高斯白噪声）

噪声的峰值是概率性的，通常取99.9%概率区间对应的幅度作为统计峰值：

PAPR≈6.6(约8.4 dB)

这是基于高斯分布6.6倍标准差（σ）对应的概率区间

(3) OFDM多载波信号

理论最大峰均比随子载波数N线性增长：

PAPRmax=10log10(N)(dB)

当N较大时，信号近似高斯分布，实际PAPR需用统计方法分析

统计量为CCDF

为了评估PAPR性能，通常绘制CCDF曲线，即PAPR超过某个门限值的概率。

CCDF曲线横坐标为PAPR值（dB），纵坐标为概率（PAPR>门限）。

在MATLAB中仿真PAPR的步骤如下：

生成信号 ：例如，生成OFDM信号或其他需要分析PAPR的信号。

计算每个符号的PAPR ：对于每个生成的信号符号，计算其瞬时PAPR。

统计多个符号 ：通过生成大量的信号符号，统计PAPR的分布，通常用互补累积分布函数（CCDF）表示。

OFDM信号生成 ：

生成随机的QAM符号。

进行IFFT变换得到时域信号。

可以添加循环前缀（但在PAPR计算中，循环前缀通常不包括在内，因为它不改变峰均比）。

PAPR计算 ：

对于每个OFDM符号（去掉循环前缀后），计算其峰值功率（即最大的瞬时功率）和平均功率。

然后计算该符号的PAPR。

% 参数设置

N = 1024;          % OFDM符号数

K = 256;           % 子载波数

M = 16;            % 调制阶数（如16-QAM）

cp_len = 16;       % 循环前缀长度

% 1. 信号生成

data = randi([0 M-1], K, N);          % 随机数据

mod_data = qammod(data, M);           % QAM调制

ifft_data = ifft(mod_data, K);        % IFFT变换

tx_signal = [ifft_data(end-cp_len+1:end,:); ifft_data]; % 加循环前缀

% 2. 峰均比计算

signal_power = abs(tx_signal).^2;                     % 瞬时功率

papr_dB = 10*log10(max(signal_power(:)) / mean(signal_power(:)));

% 3. 统计分析（CCDF）

num_sim = 10000;  % 仿真次数

papr_vals = zeros(1, num_sim);

for i = 1:num_sim

    % 重复上述过程...

    % 记录每次的PAPR值

end

ccdf = 1 - ecdf(papr_vals);  % 互补累积分布函数

 常见误区

恒包络调制（如QPSK）经滤波后会产生幅度波动，PAPR不再为0 dB

对于多载波系统，实际PAPR远小于理论最大值10log10(N)

峰均比（PAR，Peak-to-Average Ratio）是衡量信号质量和系统性能的一项不可或缺的指标。

峰均比计算需根据信号类型选择合适方法：

确定性信号：直接通过解析式计算

随机信号：采用统计方法，通过大量样本仿真获取CCDF曲线

多载波系统：理论估算结合仿真验证