无线传感器网络的能量控制分簇发送数据包 - matlab算法设计

下载频道> 资源分类> matlab源码> 算法设计> 无线传感器网络的能量控制分簇发送数据包

标题：无线传感器网络的能量控制分簇发送数据包
分享到：

所属分类: 算法设计资源类型：文件大小: 192.66 KB 上传时间: 2016-01-25 19:14:56 下载次数: 5 资源积分：1分提供者: xiaopeng2 图像配准源码程序

内容：
无线传感器网络的能量控制分簇发送数据包，程序员在编程的过程中可以参考学习使用，希望对IT程序员有用，此源码程序简单易懂、方便阅读，有很好的学习价值！

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

% %

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

% %

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

clear;

close all;

clc;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 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(1,1)*xm;%坐标

XR(i)=S(i).xd;

S(i).yd=rand(1,1)*ym;

YR(i)=S(i).yd;

S(i).G=0;

%initially there are no cluster heads only nodes

S(i).type='N';%普通节点

temp_rnd0 = i;

%Random Election of Normal Nodes

%if (temp_rnd0>=m*n+1)

S(i).E=Eo;

S(i).ENERGY=0;

% plot(S(i).xd,S(i).yd,'o');

% hold on;

% end

%Random Election of Advanced Nodes

% if (temp_rnd0<m*n+1)

% S(i).E=Eo*(1+a)

% S(i).ENERGY=1;

% %plot(S(i).xd,S(i).yd,'+');

% %hold on;

% end

end

S(n+1).xd=sink.x;

S(n+1).yd=sink.y;

%plot(S(n+1).xd,S(n+1).yd,'x');

%First Iteration

figure(1);

%counter for CHs

countCHs=0;

%counter for CHs per round

rcountCHs=0;

cluster=1;

countCHs;

rcountCHs=rcountCHs+countCHs;

flag_first_dead=0;

for r=0:1:rmax %主循环,每次1轮

%Operation for epoch

if(mod(r, round(1/p) )==0)

for i=1:1:n

S(i).G=0;

S(i).cl=0;

end

hold off;

%Number of dead nodes

dead=0;

%Number of dead Advanced Nodes

dead_a=0;

%Number of dead Normal Nodes

dead_n=0;

%counter for bit transmitted to Bases Station and to Cluster Heads

packets_TO_BS=0;

packets_TO_CH=0;

%counter for bit transmitted to Bases Station and to Cluster Heads

%per round

PACKETS_TO_CH(r+1)=0;

PACKETS_TO_BS(r+1)=0;

figure(1);

for i=1:1:n

%checking if there is a dead node

if (S(i).E<=0)

% plot(S(i).xd,S(i).yd,'red .');

dead=dead+1;

end

if S(i).E>0

S(i).type='N';

end

%plot(S(n+1).xd,S(n+1).yd,'x');

if (dead == n)%节点全部死亡退出循环

break;

end

STATISTICS(r+1).DEAD=dead;

DEAD(r+1)=dead;

DEAD_N(r+1)=dead_n;

DEAD_A(r+1)=dead_a;

%When the first node dies

if (dead==1)

if(flag_first_dead==0)

first_dead=r

flag_first_dead=1;

end

countCHs=0;

cluster=1;

for i=1:1:n

if(S(i).E>0)

temp_rand=rand;

if ( (S(i).G)<=0) %如果该节点在候选集合中

%Election of Cluster Heads

if( temp_rand <= (p/(1-p*mod(r,round(1/p)))))

countCHs = countCHs+1;

S(i).type = 'C';

S(i).G = round(1/p)-1;

C(cluster).xd = S(i).xd;

C(cluster).yd = S(i).yd;

%plot(S(i).xd,S(i).yd,'k*');

distance=sqrt( (S(i).xd-(S(n+1).xd) )^2 + (S(i).yd-(S(n+1).yd) )^2 );%到sink的距离

C(cluster).distance = distance;

C(cluster).id = i;

X(cluster)=S(i).xd;

Y(cluster)=S(i).yd;

cluster=cluster+1;

%%广播自成为簇头

distanceBroad = sqrt(xm*xm+ym*ym);

if (distanceBroad > do)

S(i).E = S(i).E- ( ETX * ctrPacketLength + Emp* ctrPacketLength*( distanceBroad*distanceBroad*distanceBroad*distanceBroad ));%广播自成为簇头

else

S(i).E = S(i).E- ( ETX * ctrPacketLength + Efs * ctrPacketLength*( distanceBroad*distanceBroad));

end

%Calculation of Energy dissipated 簇头自己发送数据包能量消耗

distance;

if (distance > do)

S(i).E = S(i).E- ( (ETX+EDA)*(packetLength) + Emp * packetLength*( distance*distance*distance*distance ));

else

S(i).E = S(i).E- ( (ETX+EDA)*(packetLength) + Efs * packetLength*( distance * distance ));

end

packets_TO_BS = packets_TO_BS+1;

PACKETS_TO_BS(r+1) = packets_TO_BS;

end

STATISTICS(r+1).CLUSTERHEADS = cluster-1;%统计第r轮簇头数目,r是从0开始的,所以加1;cluster最后要-1,是应为上面的循环多加了1

CLUSTERHS(r+1)= cluster-1;

%Election of Associated Cluster Head for Normal Nodes

for i=1:1:n

if ( S(i).type=='N' && S(i).E>0 ) %普通节点

% min_dis = sqrt( (S(i).xd-S(n+1).xd)^2 + (S(i).yd-S(n+1).yd)^2 );%默认距离是到sink的距离

min_dis = INFINITY;

if(cluster -1 >= 1)%如果有簇头存在

min_dis_cluster = 1;

%加入最近的簇头

for c = 1:1:cluster - 1 %簇头数量一共是cluster - 1

%temp = min(min_dis,sqrt( (S(i).xd - C(c).xd)^2 + (S(i).yd - C(c).yd)^2 ) );

temp = sqrt( (S(i).xd - C(c).xd)^2 + (S(i).yd - C(c).yd)^2 );

if ( temp < min_dis )

min_dis = temp;

min_dis_cluster = c;

end

%接收簇头发来的广播的消耗

S(i).E = S(i).E - ETX * ctrPacketLength;

end

%Energy dissipated by associated Cluster Head普通节点发送数据包到簇头消耗,和加入消息

min_dis;

if (min_dis > do)

S(i).E = S(i).E - ( ETX*(ctrPacketLength) + Emp * ctrPacketLength*( min_dis * min_dis * min_dis * min_dis)); %向簇头发送加入控制消息

S(i).E = S(i).E - ( ETX*(packetLength) + Emp*packetLength*( min_dis * min_dis * min_dis * min_dis)); %向簇头数据包

else

S(i).E = S(i).E - ( ETX*(ctrPacketLength) + Efs*ctrPacketLength*( min_dis * min_dis)); %向簇头发送加入控制消息

S(i).E = S(i).E - ( ETX*(packetLength) + Efs*packetLength*( min_dis * min_dis)); %向簇头数据包

end

S(i).E = S(i).E - ETX*(ctrPacketLength); %接收簇头确认加入控制消息

%Energy dissipated %簇头接收簇成员数据包消耗能量,接收加入消息和和确认加入消息

if(min_dis > 0)

S(C(min_dis_cluster).id).E = S(C(min_dis_cluster).id).E - ( (ERX + EDA)*packetLength ); %接受簇成员发来的数据包

S(C(min_dis_cluster).id).E = S(C(min_dis_cluster).id).E - ERX *ctrPacketLength ; %接收加入消息

if (min_dis > do)%簇头向簇成员发送确认加入的消息

S(C(min_dis_cluster).id).E = S(C(min_dis_cluster).id).E - ( ETX*(ctrPacketLength) + Emp * ctrPacketLength*( min_dis * min_dis * min_dis * min_dis));

else

S(C(min_dis_cluster).id).E = S(C(min_dis_cluster).id).E - ( ETX*(ctrPacketLength) + Efs * ctrPacketLength*( min_dis * min_dis));

end

PACKETS_TO_CH(r+1) = n - dead - cluster + 1; %所有的非死亡的普通节点都发送数据包

end

S(i).min_dis = min_dis;

S(i).min_dis_cluster = min_dis_cluster;

end

%hold on;

countCHs;

rcountCHs = rcountCHs + countCHs;

%Code for Voronoi Cells

%Unfortynately if there is a small

%number of cells, Matlab's voronoi

%procedure has some problems

%[vx,vy]=voronoi(X,Y);

%plot(X,Y,'r*',vx,vy,'b-');

% hold on;

% voronoi(X,Y);

% axis([0 xm 0 ym]);

end

x=1:1:r;

y=1:1:r;

z=1:1:r;

for i=1:r;

x(i)=i;

y(i) = n - STATISTICS(i).DEAD;

z(i)=CLUSTERHS(i);

end

%plot(x,y,'r',x,z,'b');

plot(x,y,'r');

hold on;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% STATISTICS %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% %

% DEAD : a rmax x 1 array of number of dead nodes/round

% DEAD_A : a rmax x 1 array of number of dead Advanced nodes/round

% DEAD_N : a rmax x 1 array of number of dead Normal nodes/round

% CLUSTERHS : a rmax x 1 array of number of Cluster Heads/round

% PACKETS_TO_BS : a rmax x 1 array of number packets send to Base Station/round

% PACKETS_TO_CH : a rmax x 1 array of number of packets send to ClusterHeads/round

% first_dead: the round where the first node died

% http://hi.baidu.com/benjevons/blog/item/b06533dfe0fc6b106327983a.html

文件列表（点击上边下载按钮，如果是垃圾文件请在下面评价差评或者投诉）:

图像配准源码程序/
图像配准源码程序/code/
图像配准源码程序/code/demons2d/
图像配准源码程序/code/demons2d/data/
图像配准源码程序/code/demons3d/
图像配准源码程序/code/demons3d/data/
图像配准源码程序/code/demons2d/data/heart-110.png
图像配准源码程序/code/demons2d/data/heart-64.png
图像配准源码程序/code/demons2d/data/lace1.png
图像配准源码程序/code/demons2d/data/lace2.png
图像配准源码程序/code/demons2d/data/lenag1.png
图像配准源码程序/code/demons2d/data/lenag2.png
图像配准源码程序/code/demons2d/data/lenag3.png
图像配准源码程序/code/demons2d/data/statue-rio-deformed.png
图像配准源码程序/code/demons2d/data/statue-rio.png
图像配准源码程序/code/demons2d/demons.m
图像配准源码程序/code/demons3d/compose.m
图像配准源码程序/code/demons3d/data/im1.mat
图像配准源码程序/code/demons3d/data/im2.mat
图像配准源码程序/code/demons3d/demons.m
图像配准源码程序/code/demons3d/expfield.m
图像配准源码程序/code/demons3d/findupdate.m
图像配准源码程序/code/demons3d/imagepad.m
图像配准源码程序/code/demons3d/imgaussian.m
图像配准源码程序/code/demons3d/iminterpolate.m
图像配准源码程序/code/demons3d/jacobian.m
图像配准源码程序/code/demons3d/register.m
图像配准源码程序/code/demons3d/resize.m
图像配准源码程序/code/demons3d/runme.m
图像配准源码程序/code/demons3d/showgrid.m
图像配准源码程序/code/demons3d/showimage.m
图像配准源码程序/code/demons3d/showvector.m

关键词：传感器 能量 无线

matlab源码下载排行

回到顶部

友情链接

matlab源码网站源码编程语言网络编程手机通讯编程 KUNAG 官网 KUNAG 电商更多

: 联系方式| 版权声明| 招聘信息| 广告服务| 银行汇款| 法律顾问| 兼职技术| 付款方式| 关于我们|; 网站客服程序员兼职招聘 411241577@qq.com; 沪ICP备19040327号-3; 公安备案号：沪公网安备 31011802003874号; 库纳格流体控制系统（上海）有限公司版权所有; Copyright © 1999-2014, GUSUCODE.COM, All Rights Reserved

关键词： 传感器 能量 无线

相关推荐

matlab源码下载排行

友情链接

关键词：传感器能量无线