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