dv-hop算法定位误差与半径的关系.docVIP

%第二个程序 定位误差与半径关系 节点分布图 %~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~??DV-Hop算法??(定位误差与通信半径关系)~~~~~~~~~~~~~~~~~~~~~~~~%~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~% BorderLength-----正方形区域的边长，单位：m% NodeAmount-------网络节点的个数% BeaconAmount---信标节点数% Sxy--------------用于存储节点的序号，横坐标，纵坐标的矩阵%Beacon----------信标节点坐标矩阵;BeaconAmount*BeaconAmount%UN-------------未知节点坐标矩阵;2*UNAmount% Distance------未知节点到信标节点距离矩阵;2*BeaconAmount%h---------------节点间初始跳数矩阵%X---------------节点估计坐标初始矩阵,X=[x,y]% R------------------节点的通信距离，一般为10-100mclear,close all;BorderLength=100;NodeAmount=100;BeaconAmount=10;UNAmount=NodeAmount-BeaconAmount;R=20;% D=zeros(NodeAmount,NodeAmount);%未知节电到信标节点距离初始矩阵；BeaconAmount行NodeAmount列 for p=1:500 h=zeros(NodeAmount,NodeAmount);%初始跳数为0；BeaconAmount行NodeAmount列X=zeros(2,UNAmount);%节点估计坐标初始矩阵Y=5; %%运算次数，R=R*Y%~~~~~~~~~~~~~~~~~~~~~~~~~~~~~在正方形区域内产生均匀分布的随机拓扑~~~~~~~~~~~~~~~~~~~~C=BorderLength.*rand(2,NodeAmount);%带逻辑号的节点坐标Sxy=[[1:NodeAmount];C];%Sxy(2,1:BeaconAmount)=[0 100 0 100 50]; % 25 25 75 75];%Sxy(3,1:BeaconAmount)=[0 0 100 100 50];Beacon=[Sxy(2,1:BeaconAmount);Sxy(3,1:BeaconAmount)];%信标节点坐标UN=[Sxy(2,(BeaconAmount+1):NodeAmount);Sxy(3,(BeaconAmount+1):NodeAmount)];%未知节点坐标%画出节点分布图plot(Sxy(2,1:BeaconAmount),Sxy(3,1:BeaconAmount),r*,Sxy(2,(BeaconAmount+1):NodeAmount),Sxy(3,(BeaconAmount+1):NodeAmount),k.)xlim([0,BorderLength]);ylim([0,BorderLength]);title(* 红色信标节点 . 黑色未知节点)%~~~~~~~~~~~~~~~~~~~~~~~~~~~~初始化节点间距离、跳数矩阵~~~~~~~~~~~~~~~~~~~~~~%Accuracy=zeros(1,Y);for Z=1:Y? ? Dall=zeros(NodeAmount,NodeAmount);? ?for i=1:NodeAmount? ???for j=1:NodeAmount? ?? ???Dall(i,j)=((Sxy(2,i)-Sxy(2,j))^2+(Sxy(3,i)-Sxy(3,j))^2)^0.5;%所有节点间相互距离? ?? ?? ?if (Dall(i,j)=R)(Dall(i,j)0)? ?? ?? ?? ?h(i,j)=1;%初始跳数矩阵? ?? ?? ?elseif i==j? ?? ?? ?? ?h(i,j)=0;? ?? ?? ?else h(i,j)=inf;? ?? ?? ?end? ???end? ?end%~~~~~~~~~~~~~~~~~~~~~~~~~最短路经算法计算节点间跳数~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~? ?for k=1:NodeAmount? ???for i=1:NodeAmoun