leach路由协议的MATLAB代码

NodeNums = 100; % the num of node 
AreaR = 100 ;   % the area of simulate
NodeTranR=10;   %  the transit Radius
Elec=50 * 10^（-9）; %
Eamp=100*10^（-12）; 
Bx=50;  % The Postion of baseation
By=175;
MaxInteral =700; % the leach simulate time
Pch=0.05;  % the desired percentage of cluster heads 
InitEn=0.5;  % the init energy of all node
Tr=30;  
 TDMA=100;
Kbit=2000; % the bits of a node transmiting a packet every time
BandWitch = 1*10.^（6）; %  Channel Bandwitch
TOS_LOCAL_ADDRESS = 0; 
for i=1:（MaxInteral）
    AliveNode（i）=NodeNums;
    AmountData（i）=0;
end
sym alldata;
alldata=0;
LAECH = zeros（1MaxInteral）;
LAENO = zeros（1MaxInteral）;  
for i=1:1:NodeNums 
    EnNode（i）=InitEn;  % the init energy of all node
    StateNode（i）=1;    % the State of all node 1: alive 0:dead
    ClusterHeads（i）=0;  % the Set of Cluster Head 1: cluster head 0 :node
    
    Rounds=0;  % the round
end
Threshold=0;    % the threshold of node becoming a cluster-head

    Node.x=