峰值功率比

% SLM--PAPR Simulation

clear all; clc; close all;

K = 128;                                                                    % SIZE OF FFT 
N = 6;                                                                      % NUMBER OF SELECTIONS
QPSK_Set  = [1 -1 i -i];
Phase_Set = [1 -1];

MAX_SYMBOLS  = 1e4;
PAPR_Orignal = zeros（1MAX_SYMBOLS）;
PAPR_SLM     = zeros（3MAX_SYMBOLS）;

X     = zeros（NK）;
Index = zeros（NK）;
for nSymbol=1:MAX_SYMBOLS
    Index（1:）   = randint（1Klength（QPSK_Set））+1;
    Index（2:N:） = randint（N-1Klength（Phase_Set））+1;
    
    X（1:） = QPSK_Set（Index（1:））;                                          % Orignal Frequency domain signal
    Phase_Rot = Phase_Set（Index（2:N:））;
    X（2:N:） = repmat（X（1:）N-11）.*Phase_Rot;                             % Phase roated Frequency domain signal
    
    x = ifft（X[]2）;