基于matlab的无线传感器网覆盖算法

本文相关程序如下: function [tot, act]=ran_rudun_dot (n, r) %input % n -the numbers of the dot, 1X1 % r -the radius of the dot, 1X1 %output % y -the numbers of the sleep dot, 1X1 %generate the rand dot %rand (’ state',0); dot=rand (n,2) *100; %xyn=round (sqrt (n)); % for i=0:(xyn-1) % dot(i*xyn+1 : (i+1)*xyn, 1)=100/xyn氺i ; % k=1 ; % for j = i*xyn+1: (i+1)氺xyn % dot(j,2)=100/xyn*k; % k:k+1; % end % end %to make a flag about being(1) or not (0) sleeping dot(:,3)=0; %n=xyn*xyn; for i=1:n if (dot (i, 3)==0) angle=[]; maxd=2*pi ; k=1 ; adj=[]; % begin: find the adjacent dot set and the repetitive dot for j=1:n d=norm ((dot (j,1:2) -dot (i,1:2))); if (j~=i && dot (j,3)==0 && d=:0) dot (i, 3)=1 ; break; end if(j~=i && dot (j ,3)==0 && (d〈=r) && d>0) adj (k,1)=j ; adj(k,2)=d; k=k+1; end end % begin: calculate the angle when the numbers of adjacent dot 〉=3 if size(adj, 1)>=3 for k=1:size (adj,1) a=atan ((dot (adj (k, 1),2)-dot (i,2)) /(dot (adj (k, 1),1)-dot (i,1))); b=acos (adj (k,2) / (2*r)); %a-b %a+b if (a-b〈0) %union (un ion ([0:0. 1 : (a+b) ],a+b), union ([ (2*pi +a_b) :0. 1 :(2*p i)],2*p i)) angIe=union(angIe,union (union ([0:0.1:(a+b) ],a+b),un ion ([ (2*p i +a-b) :0.1: (2*p i) ],2*pi))); else %union ([(a-b) :0.1: (a+b) ],a+b) ange=union(angle,union([(a-b):0.1:(a+b)],a+b)); end end %angIe diffn= [ang I e, 2*p i ] - [0,ang I e]; maxd=max (diffn); if maxd<=1.5 dot (i,3)=1 ; t=i ; end end % end : calculate end end tot=size (dot, 1); act=tot-sum (dot (:,3)); %begin: plot the initial dot and sleep dot % k=1 ; % sleep一dot=[]; % for i=1:n % if dot(1,3)^1 % sleep—dot (k,:)=dot (i,:); % k=k+1; % end % end plot (dot (:, 1),dot (:, 2),'.'); % if sum (dot (:,3))〉0 % hold on % plot (sleep—dot (:,1),sleep—dot (:,2),'+r'); % hold off % end % end : plot %begin:plot the active dot k=1 ; act_dot=[]; for i=1:n if dot(i, 3)==0 act_dot (k, :)=dot(i,:); k=k+1; end end %if sum (do (:, 3)) >0 hold on plot (act_dot (:, 1),act_dot (:,2),'+r'); title([‘ tot is ’,num2str (size (dot, 1)),', act is num2str (size (act_dot,1))]); hold off % size (act_dot, 1); %end %end: plot end