matlab分簇教程,leach分簇算法仿真(基于matlab)

leach分簇算法仿真(基于matlab)

2007-12-07 21:46

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% %

% 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://www.doczj.com/doc/37103e42be1e650e52ea9911.html/sep %

% %

% To report your comment or any bug please send e-mail to: %

% gsmaragd@http://www.doczj.com/doc/37103e42be1e650e52ea9911.html %

% %

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% %

% 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-670,2002 %

% %

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

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: