无线传感器网络能量模型

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

% Mainframe of Comparision of LDEE-DGA LEACH and Multi-path energy aware
% routing algorithm.

%Author: SunHongyu

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

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

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

% We assume that the three algorithm are in the same situatuons:
% rx and ry represents the distribution region of sensor nodes;
% R represents the Radius of each sensor node
% Einit represents the initial energy of each sensor nodes 
% Eelect represents the energy consumption of /bit/m;
% EDA represents the energy consumotion of data aggregation;
% Etx and Erx represent the energy consumption of kbit packets
% Eamp represents the energy consumption of kbit*average hops of routing
% algorithm.
% N represents the number of the nodes in the whole distributuon.
% Ecal represents the energy consumption of calculation.
clear all;

rx=100;
ry=100;

R=20;
Einit=0.5;
Eelect=0.00000001;
EDA=0.0000000000001;
k=1024;
A_hops=rx/R;

Etx=k*Eelect;
Erx=k*Eelect;
Eamp=k*A_hops*Eelect;
Ecal=Eelect/30; % 3000Ecal=100Eelect

N=100;
Rate=0.1;

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

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

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

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%Distributions%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% This part discribe the simulation of the distribution of the three
% algorithm

% Figure（1） represents the network model of Multi-path energy aware routing
% algorithm.Thusthe plant network model with only one model

% Figure（2） reprensents the network model of leach algorithm Thusthe
% hierarchy network model.

% LDEE-DGA this paper proposed focus on the plant network model thus the
% LDEE-DGA is under the sutuations as shown in figure 1

%%%%%%%%%%%%%%%%%%%%%%%%%%Distribution%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


figure（1）;

Sink.x=0.5*rx;
Sink.y=0.5*ry;

plot（Sink.xSink.y‘r*‘）;
hold on;

for i=1:1:N
    Sensor（i）.x=rand（11）*rx;
    Sensor（i）.y=rand（11）*ry;
    plot（Sensor（i）.xSensor（i）.y‘bo‘）;
    hold on;
end


figure（2）;

Sink.x=0.5*rx;
Sink.y=0.5*ry;

plot（Sink.xSink.y‘r*‘）;
hold on;

for i=1:1:N
    Sensor（i）.x=rand（11）*rx;
    Sensor（i）.y=rand（11）*ry;
    Sensor（i）.E=Einit;
end 
for i=1:1:N*Rate;
    plot（Sensor（i）.xSensor（i）.y‘b*‘）;
    hold on;
end
for i=N*Rate+1:1:N
    plot（Sensor（i）.xSensor（i）.y‘bo‘）;
    hold on;
end


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

%%%%%%%%%%%%%%%%%%%%%%%%%%EnergyComparision%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% This part‘s function is to prove the energy consumption of LDEE-DGA-a is
% lower than Multi-path energy a