PTS matlab 代碼源

clc;clear all;close all;

N=128;%FFT_size
v=4;%number of disjoint subblocks
M=4;% M-psk
len=100;%number of symbol
papr_0=zeros（1len）;

%% original
for n=1:len
   msg=randi（[0 M-1]1N）;
   tx=pskmod（msgM）;
   x=ifft（tx[]2）;
   power0=abs（x.^2）;
   peak_power0=max（power0）;
   mean_power0=mean（power0）;
   papr_0（n）=10*log10（peak_power0./mean_power0）;
end

%%  PTS（綟だ澄）

phase_factor=[1 -1];% exp（2*j*pi*（w-1））
W=2;%す砛跑て计秖

Choose = [1 1 1 1; 1 1 1 2; 1 1 2 1; 1 2 1 1; 2 1 1 1;...
          1 1 2 2; 1 2 1 2; 1 2 2 1; 2 2 1 1; 2 1 2 1; 2 1 1 2;...      
          2 2 2 1; 2 2 1 2; 2 1 2 2; 1 2 2 2; 2 2 2 2];% b
Choose_Len = 16;% W^v 2^4 

papr_pts=zeros（1len）;
for n=1:len
   msg_p=randi（[0 M-1]1N）;
   tx_p=pskmod（msg_pM）;
   m=zeros（vlength（tx_p））;%V*N
   sub=N/v;
   for k=1:v %partition
       m（k:）=[zeros（1（k-1）*sub）tx_p（（k-1）*sub+1:k*sub）zeros（1（N-（k*sub）））]
   end
   xp=ifft（m[]2）; 
   
   min_ini=10;%砞﹚papr程﹍
   %碝т程ㄎ
   for k2=1:Cho