%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%??????????????????????????????????????????????????????????????????????%
%?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?=?200;
ym?=?200;

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

%Number?of?Nodes?in?the?field
n?=?200

%Optimal?Election?Probability?of?a?node
%to?become?cluster?head
p=0.05;
packetLength?=?4000;%數據包長度?
ctrPacketLength?=?100;%控制包長度
%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=1;
%\alpha
a=1;
INFINITY?=?999999999999999;
%maximum?number?of?rounds
rmax=2000

%%%%%%%%%%%%%%%%%%%%%%%%%?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）*y