降低ofdm系统中papr改进pts算法仿真

clear all; clc; close all;

K = 128;                                                                    % SIZE OF FFT 
V = 4;                                                                       % NUMBER OF SELECTIONS
QPSK_Set  = [1 -1 j -j];
 

MAX_SYMBOLS  = 1e5;
PAPR_Orignal = zeros（1MAX_SYMBOLS）;

%  original
for nSymbol=1:MAX_SYMBOLS
    Index = randint（1Klength（QPSK_Set））+1;%1*128
    X= QPSK_Set（Index（1:））;                                               % Orignal Frequency domain signal
    x = ifft（X[]2）;         %1*128                                              % Time domain signal
    Signal_Power0 = abs（x.^2）;
    Peak_Power0   = max（Signal_Power0[]2）;
    Mean_Power0   = mean（Signal_Power02）;
    PAPR_Orignal（nSymbol） = 10*log10（Peak_Power0./Mean_Power0）;
end

%随机交织＋遍历搜索

PAPR_PTS1 = zeros（1MAX_SYMBOLS）;

for nSymbol=1:MAX_SYMBOLS
    Index1 = randint（1Klength（QPSK_Set））+1;
    X1 = QPSK_Set（Index1（1:））;                                               % Orignal Frequency domain signal
    l1=length（X1）;
    z1=zeros（Vl1）;
    pilot1=randperm（l1/2）;
    for k1=1:V/2 
        z1（k1pilot1（k1:V/2:l1/2））=X1（pilot1（k1:V/2:l1/2））;%随机
    end
    for k1=（V/2）+1:V
        z1（k1k1-2+l1/2:V/2:l1）=X1（k1-2+l1/2:V/2:l1）;%交织
    end
    a1= ifft（z1[]2）;

  
   b（1）=1; 
   b（2）=1;
   b（3）=1;
   b（4）=1;
 
   b=[b（1） b（2） b（3） b（4）];
   aa1=sum（a1.*repmat（b.‘1K））;
     
  
    Signal_P