routing4.m

无线传感器网络的可靠路由算法

clear;
global neighbour;
global MR_reliability;
global ME_reliability;
global RE_reliability;
global link_reliability;
global link_energy;
global MR_parent;
global ME_parent;
global RE_parent;
global nei_num;
global MR_energy;
global ME_energy;
global RE_energy;
global flag;
MR_reliability=0.005*ones(100,1);
ME_reliability=ones(100,1);
RE_reliability=ones(100,1);
MR_energy=100*ones(100,1);
ME_energy=100*ones(100,1);
RE_energy=100*ones(100,1);
link_reliability=zeros(100,100);%0.999*ones(100,100);
link_energy=ones(100,100);
MR_parent=ones(100,1);
ME_parent=ones(100,1);
RE_parent=ones(100,1);
% flag=zeros(100,1);
N=100;
M=100;
neighbour=zeros(N,N*0.2);
for i=1:N
    x(i)=rand(1)*M;
    y(i)=rand(1)*M;
end
x(1)=M*0.95;y(1)=M*0.95;
for i=1:N
    nei_num(i)=1;
   for j=1:N
      if (x(i)-x(j))^2+(y(i)-y(j))^2<M^2*0.04 & i~=j
          neighbour(i,nei_num(i))=j;
          nei_num(i)=nei_num(i)+1;
      end
   end
end
n0=2;bn=1;r=2;rou=1/8/20;f=20;
MR_reliability(1)=1;ME_reliability(1)=1;link_reliability(1,1)=1;
for loop=1:20
    for i=2:100
        for m=1:nei_num(i)-1
            j=neighbour(i,m);
            d(i,j)=sqrt((x(i)-x(j))^2+(y(i)-y(j))^2);
            g(i,j)=10+10*n0*log(10/d(i,j));
            pk(i,j)=(1-1/2*exp(-10^(0.1*g(i,j))*bn/2/r))^(rou*8*f);
            link_reliability(i,j)=1-(1-pk(i,j))^3;
            if link_reliability(i,j)>0.9999
                link_reliability(i,j)=0.9999;
            end
            ctemp=link_reliability(i,j)*MR_reliability(j);
            if MR_reliability(i)<ctemp
               MR_reliability(i)=ctemp;
               MR_parent(i)=j;
            end
        end
    end
end
MR_parent(1)=1;ME_parent(1)=1;
plot(x,y,'o',x(1),y(1),'*');
for i=1:N
     cordinate_x(1)=x(i);
     cordinate_x(2)=x(MR_parent(i));
     cordinate_y(1)=y(i);
     cordinate_y(2)=y(MR_parent(i));
     line(cordinate_x,cordinate_y);     
     text(x(i)+2,y(i)+2,num2str(i));
end
hop=100*ones(N,1);hop(1)=0;k=0;computed=1;visited=zeros(N,1);l=0;
while computed<=N
    for i=1:N
        for j=1:nei_num(i)-1
             l=l+1;
             if hop(neighbour(i,j))==k & visited(i)==0
                hop(i)=k+1;
                visited(i)=1;
                computed=computed+1;
                if l>10000
                    error('graph not connected');
                end
             end
        end
    end     
    k=k+1;
end
RE_parent=ones(100,1);children=zeros(N,N*0.2);ch=ones(100,1);hop(1)=0;ME_reliability=ones(100,1);
for i=1:N
    k=i;
    for j=1:nei_num(i)-1
        etemp=neighbour(i,j);
        if hop(etemp)<hop(k) | (hop(etemp)==hop(k) & ch(etemp)<ch(k))
            k=etemp;
        end
    end
    if k~=i
       ME_parent(i)=k;
       childrenei_num(k,ch(k))=i;
       ch(k)=ch(k)+1;
    end
end
% plot(x,y,'o',x(1),y(1),'*');
% for i=1:N
%      cordinate_x(1)=x(i);
%      cordinate_x(2)=x(ME_parent(i));
%      cordinate_y(1)=y(i);
%      cordinate_y(2)=y(ME_parent(i));
%      line(cordinate_x,cordinate_y);     
%      text(x(i)+2,y(i)+2,num2str(i));
% end
% for h=0:20
%     for i=1:100
%         if hop(i)==h
%             j=ME_parent(i);
%             if i~=j
%                ME_reliability(i)=link_reliability(i,j)*ME_reliability(j);
%            else 
%                ME_reliability(i)=1;
%             end
%         end
%     end
% end   
% plot(x,y,'o',x(1),y(1),'*');
% for i=1:N
%      cordinate_x(1)=x(i);
%      cordinate_x(2)=x(RE_parent(i));
%      cordinate_y(1)=y(i);
%      cordinate_y(2)=y(RE_parent(i));
%      line(cordinate_x,cordinate_y);     
%      text(x(i)+2,y(i)+2,num2str(i));
% end
% axis([0,100,0,100]);
% i=[1:100];
% plot(i,MR_reliability,'-o',i,ME_reliability,'-o');
% legend reliability hops;
% xlabel('node ID');
% ylabel('probability');
% for i=1:10
%     count(i)=0;
%     count1(i)=0;
%     for j=1:100
%         if MR_reliability(j)>=(i-1)*0.1 && MR_reliability(j)<i*0.1
%             count(i)=count(i)+1;
%         end
%          if ME_reliability(j)>=(i-1)*0.1 && ME_reliability(j)<i*0.1
%             count1(i)=count1(i)+1;
%         end
%     end
% end
% i=[0:0.1:0.9];
%  plot(i,count,'-o',i,count1,'-o');
%  xlabel('probobility');
%  ylabel('node number');
e=ones(100,1);
for i=1:100
     j=i;
     MR_energy(i)=0;
     while j~=1
         MR_energy(i)=MR_energy(i)+link_energy(j,MR_parent(j))/e(j);
         j=MR_parent(j);
     end
end
for i=1:100
     j=i;
     ME_energy(i)=0;
     while j~=1
         ME_energy(i)=ME_energy(i)+link_energy(j,ME_parent(j))/e(j);
         ME_reliability(i)=ME_reliability(i)*link_reliability(j,ME_parent(j));
         j=ME_parent(j);
     end
end
% for i=1:100
%     a=found(i,0.85);
% end
% for i=1:100
%      j=i;
%      RE_energy(i)=0;
%      while j~=1
%          RE_energy(i)=RE_energy(i)+link_energy(j,RE_parent(j))/e(j);
%          RE_reliability(i)=RE_reliability(i)*link_reliability(j,RE_parent(j));
%          j=RE_parent(j);
%      end
% end
%  i=[1:100];
% plot(i,MR_energy,'-o',i,ME_energy,'-o',i,RE_energy,'-o');
%  plot(i,MR_reliability,'-o',i,ME_reliability,'-o',i,RE_reliability,'-o');
%  legend MR_energy ME_energy RE_energy
% %  for i=1:10
%     count(i)=0;
%     count1(i)=0;
%     for j=1:100
%         if MR_energy(j)>=(i-1)*4 && MR_energy(j)<i*4
%             count(i)=count(i)+1;
%         end
%          if ME_energy(j)>=(i-1)*4 && ME_energy(j)<i*4
%             count1(i)=count1(i)+1;
%         end
%     end
%  end
% i=[0:4:36];
% plot(i,count,'-o',i,count1,'-o');
% legend reliability hops