GPSR路由算法的MATLAB实现

GPSR基于节点地理位置路由信息,采用贪婪策略和右手准则的结合在邻居节点中选择下一跳节点进行数据转发。节点在进行路由选择时,只需知道自己、邻居和目标节点的地理位置信息,无需维护全局网络的链路状态,这在很大程度上降低了网络开销。同时,节点在数据转发时总是选择比自己近的且到目标节点距离最近的邻居节点作为下一跳。这不仅减少了所选路径的跳数,而且缩短了数据传输的路径总长度,整体降低了WSNs的能量消耗,能有效提高路径QoS时延性能和延长网络生命周期。

主程序代码如下:

%% GPSR路由算法改进

%%所有传感器节点分配消息发送到目的地

clear all;close all ;clc;%清除变量

randn('state', sum(100*clock));

%% ------------------参数设置开始------------------------------

area_length = 100;%区域

area_wide = 100;%区域

Range = 20;%节点最大传递范围

node_energy=100;%节点能源

node_down=0.4;%节点能量衰减系数

num_node = 60; %节点数

node_near=0.3;%节点紧致度,也即节点间距离不能小于0.3*Range

sent_number=20;

%% ------------------参数设置结束------------------------------

%% --------------------坐标设置开始----------------------------

%% 计算出目标节点坐标、源节点坐标、各路由节点几何位置

node=node_position(area_length,area_wide,num_node,node_energy,Range,node_near);

% 计算距离表

Dtable=distance_table(node);

%计算拓扑结构表

Ttable=topo_table(node,Range,Dtable);

x=1:num_node;

Ttable(x).topolopy;

RNG_Ttable=RNG_topology(node,Dtable,Ttable);

%% ----------------------坐标设置结束--------------------------

%% -------------------计算距离表和拓扑关系----------------------

%% -------------------贪婪算法开始-------------

route_length=zeros(sent_number,1);

for k1=1:sent_number

%随机生成源节点、目标节点

s1=randperm(num_node);

v1=find(s1==1);%源节点

v2=find(s1==2);%目标节点

%确定源节点、目标节点

for i=1:num_node

node(i).class=0;

if i==v1

node(i).class=1;

number_source=i;

end

if i==v2

node(i).class=2;

number_dem=i;

end

end

%计算距离

distance_s2d=Dtable(number_source,number_dem);%距离初始化

route=[number_source];%路径初始化

Ttable01=Ttable;%临时拓扑赋值,防止搞坏原始拓扑集

iter_number=0;

%% GPRS主函数

while distance_s2d~=0

%% 判断死链

if isempty(route)

disp('能量耗尽,无法产生路由路径');

break;

end

curr_node=route(end);%路径当前最后一个节点

curr_topolopy=Ttable01(curr_node).topolopy;%获取当前拓扑相邻集合

a03=numel(find(curr_topolopy==number_dem));

if a03>0%存在

route=[route number_dem];

break;

end

%% ------------------删除curr_topolopy中属于route的节点开始-----

n=length(route);%计算route长度

for k=1:n

v01=find(curr_topolopy==route(k));%route(k)是不是curr_topolopy的成员

if numel(v01)>0%是成员

curr_topolopy(v01)=[];%删除该元素

end

end

%% ------------------删除curr_topolopy中属于route的节点结束-----

%% ------------------删除curr_topolopy中属于比route(end)离目的节点距离远的开始-----

distance_route_end=Dtable(route(end),number_dem);

distance_topolopy=Dtable(curr_topolopy,number_dem);

v4=find(distance_topolopy<distance_route_end);%找出距离小于当前节点的节点

curr_topolopy=curr_topolopy(v4);%更新拓扑相邻集合

%% ------------------删除curr_topolopy中属于比route(end)离目的节点距离远的结束-----

if ~isempty(curr_topolopy)%拓扑相邻集合非空就找下

curr_distance=Dtable(curr_topolopy,number_dem);%获取当前拓扑相邻集合与目的地节点距离集合

[dis01,index01]=min(curr_distance);%排序

next_node=curr_topolopy(index01(1));

route=[route next_node];

distance_s2d=Dtable(number_source,route(end));%重新计算距离

else%如果是空的,进入 模式

s4=RNG_Ttable(route(end),:);

v4=find(s4==1);

if numel(v4)>0

s5=randperm(numel(v4));

v5=find(s5==1);

next_node=v4(v5);

route=[route next_node];

distance_s2d=Dtable(number_source,route(end));%重新计算距离

else

%% 从倒数第二个路径点的拓扑相邻集合删除当前这节点

if length(route)>=2

a01=Ttable01(route(end-1)).topolopy;%上一个节点

a01(find(a01==curr_node))=[];

Ttable01(route(end-1)).topolopy=a01;

%% 从当前路径中删除当前路径点

route(end)=[];%删掉最后一个

else

disp('能量耗尽,无法产生路由路径');

break;

end

end

end

iter_number=iter_number+1;

end

disp('route of GPRS')

route

%% 更新节点能量

[k31,k32]=size(route);

for i=1:k32

b01=node_down*node_energy;

node(route(i)).energy=node(route(i)).energy-b01;

end

%% 更新拓扑表

Ttable=topo_table(node,Range,Dtable);

%% 更新RNG_topology

% RNG_Ttable=RNG_topology(node,Dtable,Ttable);%这个太慢,因为有多层嵌套循环,计算距离,运算量太大

RNG_Ttable=RNG_topology_opt(node,RNG_Ttable);

%绘图

s='route of GPRS';

%绘制路由图

draw_GPRS_route(node,route,Range,s);

route_length(k1)=length(route);

end

%% -------------------贪婪算法结束-------------

route_length

程序结果:

参考文献:

1.改进的GPSR模型及其仿真分析,吴三斌,王小明,杨 涛,付 红

相关推荐
不去幼儿园20 分钟前
【MARL】深入理解多智能体近端策略优化(MAPPO)算法与调参
人工智能·python·算法·机器学习·强化学习
Mr_Xuhhh22 分钟前
重生之我在学环境变量
linux·运维·服务器·前端·chrome·算法
盼海1 小时前
排序算法(五)--归并排序
数据结构·算法·排序算法
网易独家音乐人Mike Zhou5 小时前
【卡尔曼滤波】数据预测Prediction观测器的理论推导及应用 C语言、Python实现(Kalman Filter)
c语言·python·单片机·物联网·算法·嵌入式·iot
Swift社区8 小时前
LeetCode - #139 单词拆分
算法·leetcode·职场和发展
Kent_J_Truman9 小时前
greater<>() 、less<>()及运算符 < 重载在排序和堆中的使用
算法
IT 青年9 小时前
数据结构 (1)基本概念和术语
数据结构·算法
Dong雨10 小时前
力扣hot100-->栈/单调栈
算法·leetcode·职场和发展
SoraLuna10 小时前
「Mac玩转仓颉内测版24」基础篇4 - 浮点类型详解
开发语言·算法·macos·cangjie
liujjjiyun10 小时前
小R的随机播放顺序
数据结构·c++·算法