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>