reliability1.m

来自「无线传感器网络的可靠路由算法」· M 代码 · 共 36 行

r=zeros(6,11);
for l=1:6
    x=0;
    for p=0:0.1:1
        x=x+1;
        suc=0;
        for i=1:1000
            count=0;
            count1=0;
            for j=1:8
                a=rand();
                if a>p
                   count=count+1;
                end
            end
            for j=8:8+2^(l-1)
                a=rand();
                if a>p
                   count1=count1+1;
                end
            end
            if count1+count>=8
               suc=suc+8;
            else
                suc=suc+count;
            end
        end
        r(l,x)=1-suc/8000;
    end
end
p=[0:0.1:1];
plot(p,r(1,:),'-o',p,r(2,:),'-d',p,r(3,:),'-+',p,r(4,:),'-*',p,r(5,:),'->',p,r(6,:),'-x');
grid;
xlabel('raw loss rate');
ylabel('final loss rate');
legend 1 2 4 8 16 64 127;