clear;
clc;
%仿真CRC-16的漏檢概率
p='11000000000000101';
k=16;
n=k+16;
d='1111100000101010';
%d='111110000010101011111111';
D=bin2dec(d);
P_CRC=0;
num=0;


%加干擾
pb=0.1;  %誤碼率

for k1=1:n
    for k2=k1+1:n
%         for k3=k2+1:n
%         for k4=k3+1:n
%         for k5=k4+1:n
%         for k6=k5+1:n
%         for k7=k6+1:n
%         for k8=k7+1:n
        t=CRC_send(D,p);

        if t(k1)=='1'
            t(k1)='0';
        else
            t(k1)='1';
        end
        
        if t(k2)=='1'
            t(k2)='0';
        else
            t(k2)='1';
        end
        
%         if t(k3)=='1'
%             t(k3)='0';
%         else
%             t(k3)='1';
%         end
%         
%         if t(k4)=='1'
%             t(k4)='0';
%         else
%             t(k4)='1';
%         end
%         if t(k5)=='1'
%             t(k5)='0';
%         else
%             t(k5)='1';
%         end
%         if t(k6)=='1'
%             t(k6)='0';
%         else
%             t(k6)='1';
%         end
%         if t(k7)=='1'
%             t(k7)='0';
%         else
%             t(k7)='1';
%         end
%         if t(k8)=='1'
%             t(k8)='0';
%         else
%             t(k8)='1';
%         end
        cont=2;           %%%%%%%
        rt=t;
 
%接收檢錯(cuò)
        [q,l,rr,m]=mod2_div(rt,n,p,length(p));
        RR=bin2dec(rr);
        
            if RR==0
                num=num+1;
                P_CRC=P_CRC+(pb^cont)*(1-pb)^(n-cont);
            end
%         end
%         end   
%         end               %%%%%%%
%         end
%         end
%         end
    end
end

format short e;
P_CRC