clc;
clear;
close all;
N=1:1:10;
pe=0.1;
for(N=1:1:50)
p(N)=0;
if(rem(N,2))
for(i=(N+1)/2:1:N)
p(N)=p(N)+nchoosek(N,i)*pe^i;
end
end
if(rem(N,2)==0)%%假设当1、0数量相同时,为随机情况
for(i=(N)/2+1:1:N)
p(N)=p(N)+nchoosek(N,i)*pe^i;
end
p(N)=p(N)+0.5*nchoosek(N,N/2)*pe^(N/2);
end
if(p(N)<=0.001)
N
break;
end
end
plot(1:1:N,p)