降低峰均比的PTS程序

clear all;
N=1000;
num=128;
d2=zeros（112*5+1）;
d8=zeros（112*5+1）;
d16=zeros（112*5+1）;
for i=1:N
   for j=1:2*num
      if randn<0
        a（j）=1;
      else
        a（j）=0;
      end
   end
   papr2=PTS_ideal（anum4*num2）;
   papr8=PTS_ideal（anum4*num8）;
   papr16=PTS_ideal（anum4*num16）;
   k=1;
   for pap0=0:1/5:12
      if papr2>10^（pap0/10）
         d2（k）=d2（k）+1;
      end
      if papr8>10^（pap0/10）
         d8（k）=d8（k）+1;
      end
      if papr16>