无线传感器网络实验指导书

无线传感器网络

实验指导书

信息工程学院

实验一质心算法

一、实验目得

掌握合并质心算法得基本思想;

学会利用MATLAＢ实现质心算法;

学会利用数学计算软件解决实际问题。

二、实验内容与原理

无需测距得定位技术不需要直接测量距离与角度信息。定位精度相对较低,不过可以满足某些应用得需要。

在计算几何学里多边形得几何中心称为质心,多边形顶点坐标得平均值就就是质心节点得坐标。

假设多边形定点位置得坐标向量表示为p i= （xｉ,y i）T,则这个多边形得质心坐标为:

例如,如果四边形ABCD得顶点坐标分别为(x1，ｙ１),（ｘ2, y2）,(ｘ3,y３) 与(x4,y4),则它得质心坐标计算如下:

这种方法得计算与实现都非常简单,根据网络得连通性确定出目标节点周围得信标参考节点,直接求解信标参考节点构成得多边形得质心。

锚点周期性地向临近节点广播分组信息,该信息包含了锚点得标识与位置。当未知结点接收到来自不同锚点得分组信息数量超过某一门限或在一定接收时间之后,就可以计算这些锚点所组成得多边形得质心,作为确定出自身位置。由于质心算法完全基于网络连通性,无需锚点与未知结点之间得协作与交互式通信协调,因而易于实现。

三、实验内容及步骤

该程序在Maｔlab环境下完成无线传感器中得质心算法得实现。在长为１00米得正方形区域，信标节点(锚点)为90个,随机生成50个网络节点。节点得通信距离为30米。

需完成:

分别画出不同通信半径,不同未知节点数目下得误差图，并讨论得到得结果

所用到得函数:

1、M = mｉn（A)返回A最小得元素、

如果A就是一个向量,然后min(A)返回A得最小元素、

如果A就是一个矩阵,然后ｍｉｎ(A)就是一个包含每一列得最小值得行向量。

2、ｒａnd

Ｘ＝rand返回一个单一均匀分布随机数在区间(0,1)。

X =ｒand(n)返回n--n矩阵得随机数字。

3、S =ｓum(Ａ）返回Ａ沿其大小不等于1 得第一个数组维度得元素得总与。

如果Ａ就是一个向量，sum(Ａ)可返回元素得总与。

如果Ａ就是一个矩阵,然后sum(A)返回一个行向量包含每个列得总与。

４、iｎf无穷大

此ＭＡTLAB 函数返回正无穷大得IEEE 算术表示。除以零与溢出等操作会生成无穷值,从而导致结果因太大而无法表示为传统得浮点值

５、zeｒos－创建全零数组

X= ｚeｒos返回标量0、

X ＝zerｏs(n) -由－n矩阵得零返回n、

6、ｐｌot（X，Y）画出Y随Ｘ变化得2D曲线。

plot(X,Y，o)用o描述(X,Y)这一点。

7、n= nｏｒｍ（ｖ)返回得2－范数或欧氏范数得向量v、

n = norm (v,p)返回向量范数定义得sｕm(abs(v)＾p)^(1/p),这里p就是任何正值，Ｉnf 或-Inf、

８、s = nuｍ2str（A)数值数组转换为字符数组输出,它表示得数字。输出格式取决于原始值得大小。num2sｔr就是用于标签与标题情节与数字值。

所用到得变量:

xｙ:均匀分布得信标节点位置矩阵

n:未知节点数量

SS：未知节点位置矩阵

dｍ:通信半径

cent:质心

ＭＭ:未知节点估计坐标矩阵

ｅ：估计位置与实际位置距离矩阵（误差)

四、源程序

ｃlear

ｃlc

%锚节点节点设置

fｏr i=1:１：10％1到10，步长就是1；画出锚点,前一个括号就是标号。

foｒj=１:1：10

x（ｊ+（i－1）*10）＝(i-1）*10；

y(ｊ＋(i-1)*10）=(j-1)*10;

eｎd

end

fiｇure％出现图形界面

pｌoｔ(x,y,'k、＇); ％黑点

hｏld ｏn %继续画图

axｉs([0 1０00１00]）；

xｙ=[x;y]；把X，Y得坐标付给矩阵ＸY

xy;

hold on

xm=90;

ｙm=9０;

n=５0;％未知节点

foｒi＝１：1:n

Sx（i)=raｎｄ（1,1)＊xm;%产生一个一行一列得矩阵；依然就是0-1中任意一个值。

Sy(ｉ）=ranｄ(1，１）*ym;

ｐｌｏt(Sx(ｉ）,Sy(i),'r*');%红星

xlabel('x轴＇);

ｙlａｂel('y轴＇）;

hｏlｄon

end

ｄm=30;%通信半径

m＝１00;一共100个点;

ｆor ｊ＝１:１：n％未知节点循环;每一个未知节点都与在通讯范围内所有得锚点算一遍距离，放到一个矩阵中，通讯距离以外得点为零。

ＳＳ＝［Sx(j);Sy(j)];

k=0；%表示通信半径内得信标节点数

for ｉ=１:1:ｍ

d=nｏrｍ(（xy(:，ｉ)－SS），2);%２范数就就是求直线距离;不管行，第i 列;

if d＜=dm 就是否在范围内,取出信标节点,放到XX,yy中

ｘx（j，i)=ｘy（1，i);第１行第i列;%X坐标都给了XX

yy(j，i)=xy（２,i); 第2行第i列；

k=k+1；

else

xx(j,i)＝0;

ｙｙ（j，i)=0;

end

end

ｉf k~＝０%k不等于零

cent（：,j）=[sｕm（xx（ｊ，：))；sum(yy(ｊ，:))]/k；%第j行得所有列

else

ｃｅnｔ(:,j)＝0；令第j列得所有元素为零。

end

plｏｔ（cｅｎt（1,ｊ),cｅｎt(2,j),'o＇）;

ｈold on

ｐlot(［ｃｅnｔ（1,ｊ)Ｓx(ｊ)],[cent(2，ｊ）Sｙ(j）],＇R--') ；%估计位置与真实位置连接

ｔitle('Ceｎtroｉd'）；

hold on

MM=[cent(1,j);cｅｎt(2，j）];

ｅ（j)=nｏｒm((MM-SS),2）；%定位误差ｅnd

%figｕrｅ/ｄｍ

% axｉｓ([0 n ０1］)

％j=1:1：n

％plot(j,e(j) ,'－r、'）

% hold on

% tiｔlｅ（'Ceｎｔrｏid＇)

% E=suｍ(ｅ)/n

E=sum(e)/(ｎ＊dｍ）;

ｄisp(['定位误差=',num2ｓtｒ(E）])；%将结果转换成字符串输出

?实验

二

DV－

ｈop

算法

一、实

验

目

得

掌

握

DV-ho

ｐ算法

得基本思想;

学会利用MＡTLAB实现DV-hoｐ算法;

学会利用数学计算软件解决实际问题。

二、实验内容与原理

DV-Hｏp算法解决了低锚点密度引发得问题,它根据距离矢量路由协议得原理在全网范围内广播跳数与位置。

已知锚点Ｌ1与Ｌ２、L３之间得距离与跳数。L2计算得到校正值(即平均每跳距离)为(4０+75)/(2+５)=16、42m。假设传感器网络中得待定位节点A从L2获得校正值,则它与3个锚点之间得距离分别就是Ｌ１＝3×16、42,L２=２×１6、42,Ｌ3=3×1６、4２，然后使用多边测量法确定节点得位置。

图1 ＤVｈoｐ算法示意图

三、实验内容及步骤

该程序在Maｔｌａｂ环境下完成无线传感器中得DV－Hop算法得实现。在长为100米得正方形区域，随机生成1０0个网络节点,其中信标节点(锚点)为8个,未知坐标节点为9２个。节点得通信距离为50米。可算出最短路经算法计算节点间跳数、每个信标节点得校正值并可用用跳数估计距离。(100０-１000, ３00个节点,６0个信标节点)

所用到得函数：

参见上节:质心算法;

需完成：

1）运行程序，生成随机分布得节点;

2）更正程序中得错误,求每个信标节点得校正值；

3）继续运行程序,利用跳数估计待求节点得距离。

四、源程序

％~~~~～~~～~~ DＶ-Hop算法~~~~～~~～~~ ~~~~～～~~～~ ~~~~~~～~~~

% ＢorderＬengtｈ－－－--正方形区域得边长,单位:ｍ

% NoｄeＡｍount-－－--－-网络节点得个数

％BｅacoｎＡmouｎｔ---信标节点数

%Sｘy－--－--－－---－--用于存储节点得序号,横坐标，纵坐标得矩阵

%Beacon---------－信标节点坐标矩阵;BeaconＡｍounｔ*BeaconＡｍount

%ＵＮ------－-－--－-未知节点坐标矩阵;2*UNAmouｎt

% Diｓtanｃe-－---－未知节点到信标节点距离矩阵;２*BeaｃonAmount

%h--－-－--------－-节点间初始跳数矩阵

％X-－---－－-－------节点估计坐标初始矩阵,X=［x，y]'

% R------－--------－--节点得通信距离，一般为10－100m

ｃleａr,ｃlｏse aｌl;

BordeｒLenｇｔｈ=1０0;

NodｅＡｍｏunｔ＝1００;

ＢeaconAｍount＝8；

UＮAmoｕnt=NodｅＡmoｕnｔ-BｅａconAmｏunt;

Ｒ＝5０;

%D=zｅｒos（NodｅＡmount,NｏdeＡmoｕnt);%未知节电到信标节点距离初始矩阵;BeaconAmoｕnt行NodeAmｏunｔ列

ｈ=zeros(NoｄeAmoｕnｔ，NoｄｅAmoｕnｔ);%初始跳数为0;ＮｏｄeAmounｔ行NodeAmount列

X=zｅros(2，UNAmoｕnt);%节点估计坐标初始矩阵

％～～～~~~~~~在正方形区域内产生均匀分布得随机拓扑~~~～~~~~~～～~~~~~～~~~

C=BｏｒderLenｇth、*rand(2,ＮodｅAmｏｕｎｔ)；

%带逻辑号得节点坐标

Ｓｘy=[[1:ＮoｄeＡmounｔ]；C];％将1－NoｄeAmoｕnｔ与C上下合到一起,产生一个总矩阵。

Beaｃoｎ=[Ｓxy(2,1：ＢeaconAmｏuｎt）;Sxy(3，1:BeaconAmouｎt）];％信标节点坐标UN=[Sxy（2,(BeａcｏｎAmｏunt+１):NodeAｍｏunt）;Sxy(3,(BeacoｎAmoｕnｔ+１):N ｏdeＡmｏunｔ)];%未知节点坐标

%画出节点分布图

plｏt(Sｘｙ(2,1：BｅａcoｎAｍｏｕｎt),Sxy（3,1:BeaconＡmｏuｎt）,＇ｒ*',Sｘｙ(2,(ＢeacoｎAmｏunt+1):NoｄeＡmoｕnt)，Sxy(3,(ＢeａcｏnAｍount+1):NｏｄeＡｍｏunt），'k、')

ｘlｉｍ([０,BoｒdｅrLength])；

yliｍ([0,BorderLeｎgth]);

titｌｅ('＊红色信标节点、黑色未知节点')

%~~~~～～~~~~~~～～~~~~~~~～～~～～~~初始化节点间距离、跳数矩阵~~～~~~~～~~~~~～～~~~~~~~

for i＝1：NｏｄｅAｍoｕnt

ｆor j＝1:NodｅＡmｏuｎt

Dａｌl(ｉ,j)=((Ｓxy(2,i)-Ｓxy(2,j）)^２+(Sxy(3,ｉ）-Sxｙ(3,j）)^2）^０、5;%所有节点间相互距离

iｆ(Dalｌ(i，ｊ)<=R)&(Daｌl(i,j)>０）

h（ｉ,j）=1;％初始跳数矩阵

elseiｆi==j

h(i,j)=0;%自己到自己得距离

elsｅh(i，ｊ）=inｆ;

end

end

enｄ

%～~~~~～~~~~~最短路经算法计算节点间跳数~～~~～~～~～~~~~~～~~～~～

fｏｒk=1：NodeAmｏuｎｔ

ｆｏr i=１:ＮodeAmounｔ

for j＝1:NｏdeAmoｕnt

if h(i,k)+h（k,j）＜h(ｉ,j)%min（h(i,j）,h(i,k)＋h(k,j))

h(i,ｊ)=h(i,k)+h（ｋ，j);

enｄ

ｅnd

ｅnｄ

end

h

%～~~~~～~~~~~~~求每个信标节点得校正值~~~~～~~~~~～～~~~～~～~~~~~~~~ ｈ1＝ｈ(1：BｅacｏnAmouｎt，1:ＢeaｃｏnAmount);%从大矩阵中取出信标节点得跳数

D1=Dall(1：BeaconAmount,1:BeaconAmｏｕnt);%从大矩阵中取出信标节点得距离for i=1:BeaｃｏnAmount

dｈｏｐ(i，1）=sｕm(D1(ｉ,:））/ｓum（ｈ1（ｉ,:));%每个信标节点得平均每跳距离

end

D2=Ｄall(1:BeａcｏｎAmouｎt,(BｅaconAmount+１）:ＮodeAmount);%新标节点与未知节点得距离，逗号前面就是行得范围,后面就是列得范围。

for i=1:BｅaｃｏｎＡｍoｕnt

for j=１:UNＡmｏunｔ

if min(D２(:,ｊ)）==D２（i,j）

Dhop(1，j）=dhoｐ(i,j）;%未知节点从最近得信标获得校正值end

end

ｅnd

Dhop

%~～～~~~~~~~~~~~~~用跳数估计距离～~~～~~～~～~～~~~~~~～～

hoｐ1=h（1:ＢeaｃonＡmount,（BeaconAmoｕnt+1)：NｏdeAｍounｔ）％未知节点到信标跳数,BeacoｎAｍount行UNAmount列

foｒｉ=1:ＵNAｍouｎt

ｈop＝Dhop(1,i);%hｏp为从最近信标获得得校正值

Distaｎｃe(:,ｉ)=hop*hoｐ１(:,i）;%%Ｂｅacon行UN列；

end

% %～~～~～～～～~~~～～最小二乘法求未知点坐标~~~～~~~～~~~~~～~~~~～～～~～～~～～～~~~

d=Dｉstance;

ｆoｒi＝1:2

fｏｒj=1:（ＢeacｏnAmount-1)

a（i，ｊ）=Beacｏｎ(i,j)-Beaｃoｎ(i，ＢeａconAmount);

ｅnｄ

ｅｎｄ

%a为第一个信标节点与最后一个信号节点得坐标差值,保留最后一个信标节点

% ａ'

A=－２＊(a＇);

％d=d1';

for m=1:UＮAmouｎt

for ｉ=1:(BeaconAmｏunt-１)

B(i,1)=d(i，m）^2-ｄ(BeaconAmounｔ，m)^2－Beaｃon(1,i)^2+Beaｃon(１，BeaconAmount)^2-Beacon(２，i）^2+Beacｏn(2,BeacoｎAmouｎt)^2;

% 前九个信标节点到未知节点得距离平方

%最后信标节点到未知节点得距离平方

%前九个信标节点横坐标平方

% 最后信标节点到未知节点得横坐标平方

%

前九个信标节点纵坐标平方

％

最后信标节点纵坐标平方

ｅnｄ

X1=ｉnv(A＇*A)*A＇*Ｂ；%ｉｎｖ就是求逆矩阵

Ｘ(1，m)=Ｘ1(1,１);

X(2,m)=X１（2,1）;

end

%X为估计得未知节点位置

for i=1:UNＡmount

eｒroｒ(Ｚ,i)=(((X（1,ｉ）-UN(1,ｉ）)^2＋(X(2,i)-ＵN(2,i))^2)＾０、5); ％估计出得信号位置到实际位置间得距离

end

avｇｅｒrｏr(Ｚ)=suｍ(erｒor（Ｚ,:))/UＮAｍount;%平均错误

Ａccuracｙ（Z)=ａvgｅｒroｒ(Ｚ）/R; ％节点精确度

R=R+１０；

eｎd

%％～~~~~~~~～～~~～~~~～～~~~~~～画出节点精确度图～~~～~~~~~~~~

% Accｕracy

fiｇure；pｌot（20:１0：9０，Accuｒacy,'-o')

xlabｅｌ（'通信半径/ｍ');

yｌａbel('定位精度/１00%');

titｌe('节点精确度');

D Ｄ D D ＤU ＵU U U D D Ｄ D D U U U U ＵD D D D D U U U U U

U U U ＵU

U U U U U