function[]=leach（）???
%x，y軸長(zhǎng)度
????xm=100;
????ym=100;
????%基站坐標(biāo)
????sink.x=0.5*xm;
????sink.y=0.5*ym;
????%number?of?sensors
????n=100;
????%簇頭所占比例
????p=0.1;
????%初始能量
????E0=0.02;
????ETX=50*0.000000000001;%傳輸能量，每bit
????ERX=50*0.000000000001;%接收能量，每bit
????Efs=10*0.000000000001;%耗散能量，每bit
????EDA=5*0.000000000001;%融合能耗，每bit
????%融合率
????cc=0.6;
????%最大輪數(shù)
????rmax=1000;
????CM=32;%控制信息大小
????DM=4000;%數(shù)據(jù)信息大小
????figure（1）;%顯示圖片
????for?i=1:1:n
????????S（i）.xd=rand（11）*xm;
????????S（i）.yd=rand（11）*ym;
????????S（i）.E=E0;
????????S（i）.G=0;%每一輪周期結(jié)束此變量為0
????????S（i）.type=‘N‘;%節(jié)點(diǎn)類型為普通節(jié)點(diǎn)????
????????plot（S（i）.xdS（i）.yd‘o‘）;
????????hold?on;
????end
????S（i+1）.xd=sink.x;
????S（i+1）.yd=sink.y;
????plot（S（i+1）.xdS（i+1）.yd‘x‘）;
????hold?on;
????flag_first_dead=0;%第一個(gè)死亡節(jié)點(diǎn)的標(biāo)志變量
????for?r=0:1:rmax
????????r+1
????????if（mod（rround（1/p）））
????????????for?i=1:1:n
????????????????S（i）.G=0;%輪數(shù)正好是一個(gè)周期的倍數(shù)設(shè)置S（i）.G=0
????????????end
????????end
????????hold?off;
????????cluster=0;%初始簇頭數(shù)為0
????????dead=0;%初始死亡節(jié)點(diǎn)數(shù)為0
????????figure（1）;
????????for?i=1:1:n
????????????if（S（i）.E<=0）%畫出能量小于0的節(jié)點(diǎn)
????????????????plot（S（i）.xdS（i）.yd‘red.‘）;
????????????????hold?on;
????????????????dead=dead+1;???
????????????????%將死亡節(jié)點(diǎn)標(biāo)記為紅色并且死亡數(shù)加一
????????????if（dead==1）
????????????????if（flag_first_dead==0）
????????????????????first_dead=r;
????????????????????save?Itest?first_dead;
????????????????????flag_first_dead=1;
????????????????end
????????????end
????????????else
????????????????S（i）.type=‘N‘;%畫出能量不小于0的節(jié)點(diǎn)
????????????????plot（S（i）.xdS（i）.yd‘o‘）;
????????????????hold?on;
????????????end
????????end
????????plot（S（n+1）.xdS（n+1）.yd‘x‘）;
????????Dead（r+1）=dead;
????????save?ItestDead（r+1）;
????????for?i=1:1:n
????????????if（S（i）.E>0）
????????????????if（S（i）.G<=0）
????????????????????temp_rand=rand;
????????????????????if（temp_rand<=p/（1-p*mod（rround（1/p））））%?temp_rand<=T（n）選舉成為簇頭節(jié)點(diǎn)
????????????????????????S（i）.type=‘C‘;
?????