WSN DV-Hop定位的Matlab仿真代码

基于DV-Hop的无线传感网络定位算法的matlab仿真代码。研究定位误差与锚节点个数和通信半径的关系。设定在100*100的区域内，总节点数为100，分别计算出通信半径为15,25,50时，定位误差随锚节点个数3~30改变时的变化情况。

% BorderLengh---- 正方形区域的边长，单位：m
% NodeSum-------- 网络节点总数
% AnchorSum------ 锚节点数
% UnodeSum------- 未知节点数=NodeSum-AnchorSum
% Nxy------------ 节点标号{Nxy} 3*NodeSum
% Anchor--------- 锚节点坐标矩阵 2*AnchorSum
% Unode---------- 未知节点坐标矩阵 2*（NodeSum-AnchorSum）
% R-------------- 节点通信距离
%{
    struct ANode  ---- 锚节点
    { ID
      X
      Y
      hops（NodeSum）         -- 与所有节点的跳数
      Md                    -- 该节点的平均每跳距离
      dis（AnchorSum）        -- 与所有锚节点的距离
    }

    struct UNode  ---- 未知节点
     { ID
       X
       Y
       hops（AnchorSum）      -- 与所有锚节点的跳数
       xx                   -- 估计坐标
       yy
       dis（AnchorSum1）     -- 该未知节点到各锚节点的距离
      }
%}
clear all; close all;
BorderLength=100;
NodeSum=100;
R=[152550];
ASum=[36912151821242730];
times=100;       % times 次测量，取平均定位误差

for rs=1:3     % R=152550
    
for as=1:10    % 横坐标AnchorSum 变化
    clear AnchorSum; clear UnodeSum; 
    AnchorSum=ASum（as）;
    UnodeSum=NodeSum-AnchorSum;
    clear error; clear Etemp; clear Error; clear badNo;
for T=1:times 
%111111111 初始化网络拓扑，随机生成节点的坐标 111111111111---------------
clear Nxy; clear Anchor; clear Unode;
Anchor=BorderLength.*（rand（2AnchorSum）-0.5）; % 产生随机坐标锚节点
Unode=BorderLength.*（rand（2UnodeSum）-0.5）;   
Nxy=[[1:NodeSum];Anchor Unode];           % 产生带逻辑号的节点坐标
% 画出节点分布图
%{
if T==1
   figure（1）
   plot（Nxy（21:AnchorSum）Nxy（31:AnchorSum）‘r*‘） %画出锚节点分布图
   hold on; %画出未知点分布
   plot（Nxy（2（AnchorSum+1）:NodeSum）Nxy（3（AnchorSum+1）:NodeSum）‘bo‘）
   xlim（[-BorderLength/2BorderLength/2]）;
   ylim（[-BorderLength/2BorderLength/2]）;
   title（‘* 锚节点   o 未知节点‘）
   grid on;
end
%}

%2222222222 计算节点间的最小跳数 2222222222222222----------
% 初始化节点间距离、跳数矩阵
clear Dnode; clear h;
for i=1:NodeSum
    for j=1:i
        Dnode（ij）=sqrt（（Nxy（2i）-Nxy（2j））^2+（Nxy（3i）-Nxy（3j））^2）; 
        Dnode（ji）=Dnode（ij）;         % 第i个节点与第j个节点间的相互距离
        if （Dnode（ij）<=R（rs））&&（Dnode（ij）>0）
            h（ij）=1; h（ji）=1;        % 通信范围内为1跳邻居节点
        elseif i==j
            h（ij）=0;                  % 自身节点跳距为0
        else h（ij）=inf; h（ji）=inf;   % 通信范围外初始跳数为无穷大
        end
    end
end       
% 最短路径算法计算节点间跳数 
for k=1:NodeSum
    for i=1:NodeSum
        for j=1:i
            if （h（ik）+h（kj）                h（ij）=h（ik）+h（kj）;
                h（ji）=h（ij）;
            end
        end
    end
end
%===============================================================
clear ANode; clear UNode;
for i=1:AnchorSum
    ANode（i）.ID=i;
    ANode（i）.X=Anchor（1i）;
    ANode（i）.Y=Anchor（2i）;
    ANode（i）.hops=h（i:）;               % -- 与所有节点的跳数
    ANode（i）.Md=0;                      % -- 平均每跳距离初始化为0
    ANode（i）.dis=Dnode（i1:AnchorSum）;  % -- 与所有锚节点的距离
end

for i=1:UnodeSum
    UNode（i）.ID=i;
    UNode（i）.X=Unode（1i）;
    UNode（i）.Y=Unode（2i）;
    UNode（i）.hops=h（AnchorSum+i1:AnchorSum）;  % -- 与所有锚节点的跳数
    UNode（i）.xx=0;