overall1.m

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

o=zeros(9,7);p=0.2628;r=zeros(9,7);
for l=1:9
    x=0;
    for k=1:7
        x=x+1;
        suc=0;
        s=0;
        for i=1:1000
            d=0;
            left=8;
            for h=1:5
                   v(h)=0;
                   count=0;
                   count1=0;
                   m=0;
                   while v(h)==0 && m<k
                         for j=count+1:left
                             s=s+1;
                             a=rand();
                             if a>p
                                count=count+1;
                             end
                         end
                         for j=left+count1:left+(l-1)*left/8
                             s=s+1;
                             a=rand();
                             if a>p
                                count1=count1+1;
                             end
                         end
                         if count1+count>=8
                            v(h)=1;
                         else
                             m=m+1;
                         end
                   end
                   if m>=k
                       left=count;
                   end
            end
            suc=suc+left;
        end
        o(l,x)=s/5/suc;
        r(l,x)=suc/8000;
    end
end
plot(o(:,1),r(:,1),'-o',o(:,2),r(:,2),'-*',o(:,3),r(:,3),'-+',o(:,4),r(:,4),'-d',o(:,5),r(:,5),'-x',o(:,6),r(:,6),'-p',o(:,7),r(:,7),'-h');%,o(8,:),r(8,:),'-o',o(9,:),r(9,:),'-o')
axis([1,6,0,1]);
legend 0 1 2 3 4 5 5+RF;
xlabel('overhead');
ylabel('reliability');
grid;