clear
%1.初始參數(shù)設定模塊
%.傳感器節(jié)點區(qū)域界限（單位?M）
xm=150;
ym=150;
%（1）匯聚節(jié)坐標給定
sink.x=0.5*xm;
sink.y=0.5*ym;
%區(qū)域內(nèi)傳器節(jié)數(shù)
n=250
%簇頭優(yōu)化比例（當選簇頭的概率）
p=0.05;
P=0.05;
INFINITY?=?999999999999999;
%packetLength=4000;
%能量模型（單位?焦）
%初始化能量模型
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;
%高能量節(jié)點超出一節(jié)點能量的百分比
%最大循環(huán)次數(shù)
rmax=2500
%算出參數(shù)?do
do=sqrt（Efs/Emp）;
Et=0;
%2.無線傳感器網(wǎng)絡模型產(chǎn)生模塊
%構(gòu)建無線傳感器網(wǎng)絡在區(qū)域內(nèi)均勻投放100個節(jié)點并畫出圖形
for?i=1:1:n
?????S1（i）.xd=rand（11）*xm;
????S2（i）.xd=S1（i）.xd;
????S3（i）.xd=S1（i）.xd;
????S4（i）.xd=S1（i）.xd;
????S5（i）.xd=S1（i）.xd;
????XR5（i）=S5（i）.xd;
????XR4（i）=S4（i）.xd;
????XR3（i）=S3（i）.xd;
????XR2（i）=S2（i）.xd;
????XR1（i）=S1（i）.xd;
????S1（i）.yd=rand（11）*ym;
????S2（i）.yd=S1（i）.yd;
????S3（i）.yd=S1（i）.yd;
????S4（i）.yd=S1（i）.yd;
????S5（i）.yd=S1（i）.yd;
????YR5（i）=S5（i）.yd;
????S5（i）.G=0;
????YR4（i）=S4（i）.yd;
????S4（i）.G=0;
????YR3（i）=S3（i）.yd;
????S3（i）.G=0;
????YR2（i）=S2（i）.yd;
????YR1（i）=S1（i）.yd;
????S1（i）.G=0;
????S2（i）.G=0;
????S1（i）.E=Eo;
????S2（i）.E=S1（i）.E;
????S3（i）.E=S1（i）.E;
????S4（i）.E=S1（i）.E;
????S5（i）.E=S1（i）.E;
????E3（i）=?S3（i）.E;
????E4（i）=?S4（i）.E;
????Et=Et+E3（i）;

????%initially?there?are?no?cluster?heads?only?nodes
????S1（i）.type=‘N‘;
????S2（i）.type=‘N‘;
????S3（i）.type=‘N‘;
????S4（i）.type=‘N‘;
????S5（i）.type=‘N‘;
end

S1（n+1）.xd=sink.x;
S1（n+1）.yd=sink.y;
S2（n+1）.xd=sink.x;
S2（n+1）.yd=sink.y;
%3.網(wǎng)絡運行模塊
%簇頭節(jié)點數(shù)
countCHs1=0;
cluster1=1;%此定義的目的僅僅是給定一個1開始的下標參數(shù)，真正的簇頭數(shù)應該還減去1
flag_first_dead1=0;
flag_all_dead1=0;
%死亡節(jié)點數(shù)
dead1=0;
first_dead1=0;
all_dead1=0;
%活動節(jié)點數(shù)
allive1=n;
%counter?for?bit?transmitted?to?bases?Station?and?to?Cluster?Heads
packets_TO_BS1=0;
packets_TO_CH1=0;
%（1）循環(huán)模式設定
for?r=0:1:rmax?????%該?for?循環(huán)將下面的所有程序包括在內(nèi)，直到最后一?end?才結(jié)束循環(huán)
????r
??%每過一個輪轉(zhuǎn)周期（本程序為10次）使各節(jié)點的S（i）.G參數(shù)（該參數(shù)用于后面的簇選舉，在該輪轉(zhuǎn)周期內(nèi)已當選過簇頭的節(jié)點不能再當選）恢復為零
??if（mod（r?round（1/p）?）==0）
????for?i=1:1:n
????????S1（i）.G=0;
????????S1（i）.cl=0;
????end
??end
%（2）死亡節(jié)點檢查模塊
dead1=0;
for?i=1:1:n
????%檢查有無死亡節(jié)點
????if?（S1（i）.E<=0）
????????dead1=dead1+1;?
????????%（3）第一個死亡節(jié)點的產(chǎn)生時間（用輪次表示）
????????%第一個節(jié)點死亡時間
????????if?（dead1==1）
???????????if（flag_first_dead1==0）
??????????????first_dead1=r;
??????????????flag_first_dead1=1;
???????????end
????????end
??????
????????if（dead1==n）
???????????if（flag_all_dead1==0）
??????????????all_dead1=r;
??????????????flag_all_dead1=1;
???????????end
????????end
????end
????if?S1（i）.E>0
????????S1（i）.type=‘N‘;
????end
end
STATISTICS.DEAD1（r+1）=dead1;
STATISTICS.ALLIVE1（r+1）=allive1-dead1;
%（4）簇頭選舉模塊
countCHs1=0;
cluster1=1;
for?i=1:1:n
???if（S1（i）.E>0）
???temp_rand=rand;?????
???if?（?（S1（i）.G）<=0）??
???????%簇頭的選舉，當選的簇頭會把各種相關(guān)信存入下面程序所給定的變量中
????????if（temp_rand<=?（p/（1-p*mod（rround（1/p）））））
????????????countCHs1=countCHs1+1;
????????????packets_TO_BS1=packets_TO_BS1+1;
????????????PACKETS_TO_BS1（r+1）=packets

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件??????21684??2014-10-21?15:36??route\akskflhalsfhaslf.m

?????文件??????10454??2007-05-25?19:47??route\bp2.m

?????文件??????19001??2014-09-11?03:01??route\bp4.m

?????文件??????12515??2007-06-04?19:02??route\DEEC1VSDEEC2.m

?????文件???????6637??2014-09-13?00:53??route\DEEC2.m

?????文件??????13674??2014-10-19?19:39??route\duibi0000.m

?????文件???????7854??2014-09-15?17:06??route\duo1.fig

?????文件???????4542??2008-05-21?07:04??route\duotiao.m

?????文件???????9771??2014-10-15?14:58??route\EAMMH.m

?????文件????????313??2007-08-23?18:28??route\energy1.m

?????文件???????1418??2007-08-24?10:04??route\energybalancemodel.m

?????文件??????20215??2014-09-10?10:25??route\fencu.m

?????文件??????20473??2007-06-04?20:15??route\fps.m

?????文件???????8071??2014-10-18?20:38??route\gaijin000.m

?????文件???????8164??2014-10-21?10:53??route\gaijin001.m

?????文件???????4548??2010-06-26?15:03??route\ga_wsn.m

?????文件??????25539??2014-09-11?09:22??route\GJ.m

?????文件??????18776??2007-04-09?14:39??route\heed.m

?????文件???????2181??2010-12-02?10:27??route\hopbyhop.m

?????文件??????12481??2014-10-23?15:24??route\hunhe1.m

?????文件???????8244??2014-08-31?11:18??route\LEACH00.m

?????文件???????6764??2014-10-21?17:20??route\leach1.m

?????文件???????5955??2014-10-10?19:16??route\leachm1.m

?????文件??????10974??2014-08-31?11:18??route\leachm1vsleachm2.m

?????文件???????6402??2007-05-25?11:24??route\leachm2.m

?????文件???????6118??2007-05-28?22:11??route\leachs1.m

?????文件??????12137??2014-09-03?08:13??route\leachs1vsleachs2.m

?????文件???????6346??2014-10-10?19:10??route\leachs2.m

?????文件??????20471??2012-04-14?21:56??route\leachvsdeec.m

?????文件???????9012??2014-10-14?09:32??route\LEACH_ENERGY2.m

............此處省略49個文件信息