无线传感器网络的PEGASIS协议源代码

无线传感器网络PEGASIS协议在MATLAB上的仿真

function s = costofpegasis（rmax）
sum=0;

for r=1:1:rmax
%pegasis成chain与BS位置有关
%clear;
PS4=load（‘data‘）;
PS4=PS4.‘;

BSx=50;
BSy=300;

node=length（PS4）;

%按照到BS距离从近到远的顺序放置个节点坐标到PS4矩阵中，PS4（1:）中节点坐标是离BS最近的;
%PS4（i:）表示第i节点，具体是：
%PS4（i1）=x坐标，PS4（i2）=y坐标，PS4（i3）=i表示节点序号
for i=1:node
    for j=i:node
        if （PS4（i1）-BSx）^2+（PS4（i2）-BSy）^2>（PS4（j1）-BSx）^2+（PS4（j2）-BSy）^2
               s=PS4（i:）;
               PS4（i:）=PS4（j:）;
               PS4（j:）=s;
            end
        end
     PS4（i3）=i;
 end

 %d（ij）为第i节点到第j节点的距离值
for i=1:node
    for j=1:node
               if（i~=j）
                 d（ij）=sqrt（（PS4（i1）-PS4（j1））^2+（PS4（i2）-PS4（j2））^2）;
               else d（ij）=10^8;
               end
    end
end

% distance matrix

%从离BS最远节

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       3202  2009-05-15 18:04  pegasis\data

     文件       2122  2009-05-23 08:58  pegasis\costofpegasis.m

     文件         17  2009-05-24 16:21  pegasis\costofpegasis

     文件       2069  2009-05-23 08:53  pegasis\costofpegasis.asv

     目录          0  2009-05-16 09:34  pegasis

----------- ---------  ---------- -----  ----

                 7410                    5