leach经典算法matlab代码

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%                                                                      %
% SEP: A Stable Election Protocol for clustered                        %
%      heterogeneous wireless sensor networks                          %
%                                                                      %
% （c） Georgios Smaragdakis                                             %
% WING group Computer Science Department Boston University           %
%                                                                      %
% You can find full documentation and related information at:          %
% http://csr.bu.edu/sep                                                %
%                                                                      %  
% To report your comment or any bug please send e-mail to:             %
% gsmaragd@cs.bu.edu                                                   %
%                                                                      %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%                                                                      %
% This is the LEACH [1] code we have used.                             %
% The same code can be used for FAIR if m=1                            %
%                                                                      %
% [1] W.R.Heinzelman A.P.Chandrakasan and H.Balakrishnan             %
%     “An application-specific protocol architecture for wireless      % 
%      microsensor networks“                                           % 
%     IEEE Transactions on Wireless Communications 1（4）:660-6702002  %
%                                                                      %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

clear;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% PARAMETERS %%%%%%%%%%%%%%%%%%%%%%%%%%%%

%Field Dimensions - x and y maximum （in meters）
xm=100;
ym=100;

%x and y Coordinates of the Sink
%sink.x=0.5*xm;
%sink.y=0.5*ym;
sink.x =0.5 * xm;  %modified by jhf
sink.y = ym + 50;  %modified by jhf

%Number of Nodes in the field
n=100;
n=20;  %modified by jhf 11_6

%Optimal Election Probability of a node
%to become cluster head
p=0.1;

%Energy Model （all values in Joules）
%Initial Energy 
Eo=0.5;
%Eelec=Etx=Erx
ETX=50*0.000000001;
ERX=50*0.000000001;
%Transmit Amplifier types
Efs=10*0.000000000001;
Emp=0.0013*0.000000000001;
%Data Aggregation Energy
EDA=5*0.000000001;

%Values for Hetereogeneity
%Percentage of nodes than are advanced
m=0.1;
%\alpha
a=1;

%maximum number of rounds
rmax=9999;


%%%%%%%%%%%%%%%%%%%%%%%%% END OF PARAMETERS %%%%%%%%%%%%%%%%%%%%%%%%

%Computation of do
do=sqrt（Efs/Emp）;

%Creation of the random Sensor Network
figure（1）;
for i=1:1:n
    S（i）.xd=rand（11）*xm;
    XR（i）=S（i）.xd;
    S（i）.yd=rand（11）*ym;
    YR（i）=S（i）.yd;
    S（i）.G=0;
    %initially there are no cluste