temp1.asv

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

clear;
N=100;
M=100;
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;
plot(x,y,'o',x(1),y(1),'*');
neighbor=zeros(N,N*0.2);
for i=1:N
    n(i)=1;
   for j=1:N
      if (x(i)-x(j))^2+(y(i)-y(j))^2<M^2*0.04 & i~=j
          neighbor(i,n(i))=j;
          n(i)=n(i)+1;
      end
   end
end
hop=N*ones(N,1);hop(1)=0;k=0;computed=1;visited=zeros(N,1);
while computed<=N
    for i=1:N
        for j=1:n(i)-1
             if hop(neighbor(i,j))==k && visited(i)==0
                hop(i)=k+1;
                visited(i)=1;
                computed=computed+1;
             end
        end
    end     
    k=k+1;
end
parent=ones(N,1);children=zeros(N,N*0.2);c=ones(N,1);hop(1)=0;
for i=1:N
    k=i;
    for j=1:n(i)-1
        temp=neighbor(i,j);
        if hop(temp)<hop(k) | (hop(temp)==hop(k) & c(temp)<c(k))
            k=temp;
        end
    end
    if k~=i
       parent(i)=k;
       children(k,c(k))=i;
       c(k)=c(k)+1;
    end
end
 parent(1)=1;
 for i=1:N
     cordinate_x(1)=x(i);cordinate_x(2)=x(parent(i));cordinate_y(1)=y(i);cordinate_y(2)=y(parent(i));
     line(cordinate_x,cordinate_y);
 end
t=zeros(N,1);
t(1)=1;
computed=1;
token=1;
forward=ones(N,1);
assigned=zeros(N,1);assigned(1)=1;used=zeros(N,1);
max_t=1;
while computed<=N |used(
    if used(token)==0
    for i=1:c(token)-1
        temp=children(token,i);
        if assigned(temp)==0
            color=1;find=0;
            compare=parent(temp);
            while find==0
                for j=1:n(compare)-1
                    find=1;
                    if color==t(neighbor(compare,j))|color==t(compare)
                        color=color+1;
                        find=0;
                        break;
                    end
                end
            end
            t(temp)=color;
            computed=computed+1;
            assigned(temp)=1;
        end
    end
    token=children(token,forward(token));
    if c(token)==1
        token=parent(token);
    end
    if forward(token)<c(token)-1
       forward(token)=forward(token)+1;
    else
        used(token)=1;
        token=parent(token);
    end
    end
end
for i=1:N    
    text(x(i)+2,y(i)+2,num2str(t(i)));
end