[go: up one dir, main page]

CN116261200A - An MPOLSR protocol optimization method based on link transmission quality for multi-UAV flight ad hoc network - Google Patents

An MPOLSR protocol optimization method based on link transmission quality for multi-UAV flight ad hoc network Download PDF

Info

Publication number
CN116261200A
CN116261200A CN202310326422.7A CN202310326422A CN116261200A CN 116261200 A CN116261200 A CN 116261200A CN 202310326422 A CN202310326422 A CN 202310326422A CN 116261200 A CN116261200 A CN 116261200A
Authority
CN
China
Prior art keywords
node
link
neighbor
link transmission
mpolsr
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202310326422.7A
Other languages
Chinese (zh)
Inventor
魏守明
沃滋植
马琳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin Institute of Technology Shenzhen
Original Assignee
Harbin Institute of Technology Shenzhen
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Harbin Institute of Technology Shenzhen filed Critical Harbin Institute of Technology Shenzhen
Priority to CN202310326422.7A priority Critical patent/CN116261200A/en
Publication of CN116261200A publication Critical patent/CN116261200A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W40/00Communication routing or communication path finding
    • H04W40/02Communication route or path selection, e.g. power-based or shortest path routing
    • H04W40/12Communication route or path selection, e.g. power-based or shortest path routing based on transmission quality or channel quality
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L45/00Routing or path finding of packets in data switching networks
    • H04L45/12Shortest path evaluation
    • H04L45/122Shortest path evaluation by minimising distances, e.g. by selecting a route with minimum of number of hops
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W40/00Communication routing or communication path finding
    • H04W40/02Communication route or path selection, e.g. power-based or shortest path routing
    • H04W40/22Communication route or path selection, e.g. power-based or shortest path routing using selective relaying for reaching a BTS [Base Transceiver Station] or an access point
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W84/00Network topologies
    • H04W84/18Self-organising networks, e.g. ad-hoc networks or sensor networks
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00Reducing energy consumption in communication networks
    • Y02D30/70Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Data Exchanges In Wide-Area Networks (AREA)

Abstract

本发明创造提供了一种用于多无人机飞行自组网中的基于链路传输质量的MPOLSR协议优化方法,基于多无人机之间进行信息交互的MPOLSR协议优化通过节点连接度与节点链路传输数期望进行MPR节点筛选,同时以链路传输期望作为链路权值指标进行最短路计算,实现MPOLSR协议优化。本发明创造用以解决MPR筛选机制单一、最短路径算法执行复杂度较高、未区分多路权值等问题,综合考虑节点连接度与节点链路传输数期望进行MPR节点筛选,同时更换最短路算法,并以链路传输期望作为路径权值进行区分,提高网络的有效性与可靠性。

Figure 202310326422

The present invention provides an MPOLSR protocol optimization method based on link transmission quality used in multi-UAV flight ad hoc networks. The MPOLSR protocol optimization based on information interaction between multi-UAVs is achieved through node connectivity and node connectivity. The link transmission number is expected to be screened by MPR nodes, and at the same time, the link transmission expectation is used as the link weight index to perform the shortest path calculation to realize the optimization of the MPOLSR protocol. The invention is created to solve the problems of single MPR screening mechanism, high execution complexity of the shortest path algorithm, and undifferentiated multi-path weights, etc., comprehensively considering node connectivity and node link transmission number expectations to perform MPR node screening, and replace the shortest path at the same time Algorithm, and the link transmission expectation is used as the path weight to distinguish, so as to improve the effectiveness and reliability of the network.

Figure 202310326422

Description

一种用于多无人机飞行自组网中的基于链路传输质量的 MPOLSR协议优化方法A MPOLSR protocol optimization method based on link transmission quality for multi-UAV flying ad hoc networks

技术领域Technical Field

本发明属于无人机通信领域,特别是涉及一种用于多无人机飞行自组网中的基于链路传输质量的MPOLSR协议优化方法。The present invention belongs to the field of unmanned aerial vehicle communication, and in particular relates to an MPOLSR protocol optimization method based on link transmission quality for use in a multi-unmanned aerial vehicle flight ad hoc network.

背景技术Background Art

多无人机之间进行信息交互需要飞行自组网(Flying Ad-hoc Networks,FANETs),这是移动自组织网络(MobileAd-hoc Networks,MANETs)在无人机领域的一种扩展应用,属于一种特殊意义上的MANETs。相比于其他组网的信息收集方式,多无人机自组织网络有着可以快速移动的极高灵活性,无需事先布置地面基站的极强机动性以及无视各种复杂恶劣地面环境的极佳适应性。因此,多无人机自组织网络在火情监控、抗震救灾、电力输送和目标跟踪等方面都获得了越来越广泛的应用,但其具有节点移动速度快、场景范围大、拓扑变化频繁的特点,因此设计一种高稳定性的路由协议,保证数据传输的有效性和可靠性是多无人机自组网的重要研究内容之一。Information exchange between multiple drones requires Flying Ad-hoc Networks (FANETs), which is an extended application of Mobile Ad-hoc Networks (MANETs) in the field of drones and belongs to a special kind of MANETs. Compared with other networking information collection methods, multi-drone self-organizing networks have extremely high flexibility in rapid movement, strong mobility without the need for pre-arranged ground base stations, and excellent adaptability to various complex and harsh ground environments. Therefore, multi-drone self-organizing networks have been increasingly widely used in fire monitoring, earthquake relief, power transmission, and target tracking. However, they have the characteristics of fast node movement, large scene range, and frequent topology changes. Therefore, designing a highly stable routing protocol to ensure the effectiveness and reliability of data transmission is one of the important research contents of multi-drone self-organizing networks.

MPOLSR(Multi-Path Optimized Link State Routing)协议相比其他主动路由协议,引入了多中继节点(Multi Point Relay,MPR)这一概念,只有被某节点选为MPR的点才能够转发该节点的拓扑控制信息,从而大幅度减少了洪泛机制的路由开销,具有比DSDV(Destination Sequenced Distance Vector)等主动式路由协议更高的性能。同时计算多条最短路径,有效提高了可靠性,可以保证无人机集群执行任务时的效率和质量。但该协议存在MPR筛选机制单一、最短路径算法执行复杂度较高、未区分多路权值等问题。Compared with other active routing protocols, the MPOLSR (Multi-Path Optimized Link State Routing) protocol introduces the concept of multi-point relay (MPR). Only the points selected as MPR by a node can forward the topology control information of the node, thereby greatly reducing the routing overhead of the flooding mechanism and having higher performance than active routing protocols such as DSDV (Destination Sequenced Distance Vector). Calculating multiple shortest paths at the same time effectively improves reliability and can ensure the efficiency and quality of drone clusters when performing tasks. However, the protocol has problems such as a single MPR screening mechanism, high execution complexity of the shortest path algorithm, and failure to distinguish between multiple path weights.

发明内容Summary of the invention

本发明提供一种用于多无人机飞行自组网中的基于链路传输质量的MPOLSR协议优化方法,用以解决MPR筛选机制单一、最短路径算法执行复杂度较高、未区分多路权值等问题。The present invention provides an MPOLSR protocol optimization method based on link transmission quality for a multi-UAV flying ad hoc network, which is used to solve the problems of single MPR screening mechanism, high execution complexity of the shortest path algorithm, and failure to distinguish multi-path weights.

本发明通过以下技术方案实现:The present invention is achieved through the following technical solutions:

一种用于多无人机飞行自组网中的基于链路传输质量的MPOLSR协议优化方法,所述MPOLSR协议优化方法包括,A MPOLSR protocol optimization method based on link transmission quality for a multi-UAV flight ad hoc network, the MPOLSR protocol optimization method comprising:

基于多无人机之间进行信息交互的MPOLSR协议优化通过节点连接度与节点链路传输数期望进行MPR节点筛选,同时采用SPFA最短路算法,并以链路传输期望作为链路权值指标进行最短路计算,实现MPOLSR协议优化。The MPOLSR protocol optimization based on information interaction between multiple UAVs performs MPR node screening through node connectivity and node link transmission number expectation. At the same time, the SPFA shortest path algorithm is adopted, and the link transmission expectation is used as the link weight indicator for shortest path calculation to achieve MPOLSR protocol optimization.

进一步的,所述MPR节点筛选为基于网络状态感知的MPR集选择,具体包括以下步骤,Further, the MPR node screening is an MPR set selection based on network state perception, which specifically includes the following steps:

步骤S1:首先建立一份备选的MPR节点列表M集,初始化为空集;Step S1: First, establish a candidate MPR node list M set, which is initialized as an empty set;

步骤S2:遍历节点1跳邻居集建立N1,遍历节点2跳邻居集建立严格2跳邻居集N2;Step S2: traverse the node 1-hop neighbor set to establish N1, and traverse the node 2-hop neighbor set to establish a strict 2-hop neighbor set N2;

步骤S3:判断MPR集合中的点能否到达中心节点的所有严格2跳邻居点,若能覆盖则算法结束,否则进入步骤S4;Step S3: Determine whether the points in the MPR set can reach all strictly 2-hop neighboring points of the central node. If they can be covered, the algorithm ends, otherwise it goes to step S4;

步骤S4:计算N1中所有节点的权值,包括节点连接度与节点链路传输期望;Step S4: Calculate the weights of all nodes in N1, including node connectivity and node link transmission expectations;

步骤S5:遍历之前计算的所有点,找出权值最大的点;Step S5: traverse all the points calculated previously and find the point with the largest weight;

步骤S6:将权值最大的节点加入M集中;Step S6: Add the node with the largest weight to the M set;

步骤S7:在N1中移除权值最大的节点;Step S7: remove the node with the largest weight in N1;

步骤S8:在N2中移除能够通过最大的节点到达的点,之后返回步骤S3。Step S8: Remove the points in N2 that can be reached through the largest node, and then return to step S3.

进一步的,所述步骤S2中的1跳邻居集包含四个字段:Furthermore, the 1-hop neighbor set in step S2 includes four fields:

N_neighbor_main_addr代表源节点的邻居节点ip地址,N_status代表源节点与邻居节点的链路对称关系,N_willingness代表邻居节点转发源节点的消息意愿程序,ETX_R为新增加的字段,代表源节点与邻居节点的链路传输期望的倒数,通过查询修改后的节点链路信息表获取值。N_neighbor_main_addr represents the IP address of the neighbor node of the source node, N_status represents the link symmetry relationship between the source node and the neighbor node, N_willingness represents the willingness of the neighbor node to forward the message of the source node, and ETX_R is a newly added field, which represents the inverse of the link transmission expectation between the source node and the neighbor node. The value is obtained by querying the modified node link information table.

进一步的,所述节点链路信息表包含七个字段,具体为,L_local_iface_addr表示本地节点的ip地址,L_neighbor_iface_addr表示本地节点的邻居节点的ip地址,L_SYM_time记载了源节点与邻居节点对称链路的到期时间,L_ASYM_time记载了源节点与邻居节点非对称链路的到期时间,L_time表示该链路信息的到期时间,LQ_I字段与ETX_R字段为新增加的字段,分别表示节点与邻居节点的修改后的正向链路质量与传输数期望的倒数。LQ_I字段通过查询本发明新定义的节点发送信息表与节点接收信息表获取,ETX_R字段通过获取本发明修改后的HELLO消息包后获取。Further, the node link information table includes seven fields, specifically, L_local_iface_addr represents the IP address of the local node, L_neighbor_iface_addr represents the IP address of the neighbor node of the local node, L_SYM_time records the expiration time of the symmetric link between the source node and the neighbor node, L_ASYM_time records the expiration time of the asymmetric link between the source node and the neighbor node, L_time represents the expiration time of the link information, and the LQ_I field and the ETX_R field are newly added fields, which respectively represent the modified forward link quality and the expected reciprocal of the number of transmissions between the node and the neighbor node. The LQ_I field is obtained by querying the node sending information table and the node receiving information table newly defined in the present invention, and the ETX_R field is obtained by obtaining the modified HELLO message packet of the present invention.

进一步的,所述节点发送信息表与节点接收信息表具体为,通过查询本发明新定义的节点发送信息表获取Sij,节点发送信息表包括N_neighbor_main_addr字段记录邻居ip、Send_cnt字段记录发给对应ip节点的消息数;与此同时,通过查询本发明新定义的节点接收信息表获取Rij,节点接收信息表包括N_neighbor_main_addr字段记录邻居ip、Recv_cnt字段记录收到对应ip节点的消息数;由此可以计算出LQ_I值。Further, the node sending information table and the node receiving information table are specifically as follows: S ij is obtained by querying the node sending information table newly defined in the present invention, the node sending information table includes an N_neighbor_main_addr field to record the neighbor IP, and a Send_cnt field to record the number of messages sent to the corresponding IP node; at the same time, R ij is obtained by querying the node receiving information table newly defined in the present invention, the node receiving information table includes an N_neighbor_main_addr field to record the neighbor IP, and a Recv_cnt field to record the number of messages received from the corresponding IP node; thereby, the LQ_I value can be calculated.

进一步的,所述修改后的HELLO消息包包括Reserved第一保留字段、Htime发送间隔、Willingness数据转发意愿、Linkcode链路编码信息、Reserved第二保留字段、LinkMessage Size链路状态的大小信息、NeighborInterfaceAddress邻居节点接口地址与新增加的LQ_I字段,表示产生该HELLO消息的源节点与接收到此HELLO消息的节点的LQ_I大小,源节点收到消息后即可得知自己与对应邻居节点的NLQ_I大小,从而通过前述方案计算出ETX_R,将结果填充入节点链路信息表中。Furthermore, the modified HELLO message packet includes the Reserved first reserved field, Htime sending interval, Willingness data forwarding intention, Linkcode link coding information, Reserved second reserved field, LinkMessage Size link state size information, NeighborInterfaceAddress neighbor node interface address and the newly added LQ_I field, which indicates the LQ_I size of the source node that generates the HELLO message and the node that receives the HELLO message. After receiving the message, the source node can know the NLQ_I size of itself and the corresponding neighbor node, thereby calculating ETX_R through the aforementioned scheme and filling the result into the node link information table.

进一步的,所述步骤S4中MPR节点权值指标的计算具体为,指标计算公式为:Furthermore, the calculation of the MPR node weight index in step S4 is specifically as follows:

Figure SMS_1
Figure SMS_1

其中ETX_R指的是链路传输期望ETX的倒数,D指的是节点连接度;α和β是权重指标,且满足α+β=1。Where ETX_R refers to the inverse of the expected ETX of link transmission, D refers to the node connectivity; α and β are weight indicators, and α+β=1.

进一步的,所述节点链路传输数期望计算具体为,Furthermore, the expected calculation of the node link transmission number is specifically as follows:

所述链路传输期望为,The link transmission expectation is:

Figure SMS_2
Figure SMS_2

Figure SMS_3
Figure SMS_3

Figure SMS_4
Figure SMS_4

Figure SMS_5
Figure SMS_5

Figure SMS_6
Figure SMS_6

其中Rji指的是节点j收到节点i的HELLO消息数,Sij指的是节点i发送给节点j的HELLO消息数;Rij指的是节点i收到节点j的HELLO消息数,Sji指的是节点j发送给节点i的HELLO消息数;LQ代表正向链路质量,NLQ代表反向链路质量,LQ_I代表修改后的正向链路质量,NLQ_I代表修改后的反向链路质量。Where R ji refers to the number of HELLO messages received by node j from node i, and S ij refers to the number of HELLO messages sent by node i to node j; R ij refers to the number of HELLO messages received by node i from node j, and S ji refers to the number of HELLO messages sent by node j to node i; LQ represents the forward link quality, NLQ represents the reverse link quality, LQ_I represents the modified forward link quality, and NLQ_I represents the modified reverse link quality.

进一步的,所述最短路为链路传输数期望,并采用SPFA算法的最短路优化计算,具体为,Furthermore, the shortest path is the expected number of link transmissions, and the shortest path optimization calculation using the SPFA algorithm is specifically as follows:

步骤D1:自定义一个node_map结构体,结构体中含有四个成员变量,分别是每个节点的地址、从源节点到该节点的上一跳地址、中心节点到该点的距离,以及一个布尔类型变量,表示该节点是否在计算队列中,通过遍历源节点的1跳邻居集、2跳邻居集与拓扑表进行填充;Step D1: Customize a node_map structure, which contains four member variables, namely the address of each node, the previous hop address from the source node to the node, the distance from the central node to the point, and a Boolean variable, indicating whether the node is in the calculation queue, which is filled by traversing the 1-hop neighbor set, 2-hop neighbor set and topology table of the source node;

步骤D2:建立一个etx_r_map二维数组用以存储每个节点到其他任意节点的链路权值大小,通过遍历源节点的1跳邻居集、2跳邻居集与拓扑表进行填充。链路权值的大小就是节点的链路传输期望;Step D2: Create a two-dimensional array etx_r_map to store the link weights from each node to any other node, and fill it by traversing the 1-hop neighbor set, 2-hop neighbor set and topology table of the source node. The link weight is the link transmission expectation of the node;

步骤D3:设置循环变量,初始化为0;Step D3: Set the loop variable and initialize it to 0;

步骤D4:判断循环变量的值是否达到了设定的多径路由数量,若达到则结束算法,反之进入步骤D5;Step D4: Determine whether the value of the loop variable has reached the set number of multipath routes. If so, the algorithm ends. Otherwise, it goes to step D5.

步骤D5:设定一个计算队列用于存储待更新的节点,同时将当前源节点放入计算队列中;Step D5: Set a calculation queue to store the nodes to be updated, and put the current source node into the calculation queue;

步骤D6:判断计算队列是否为空,若为空则进入步骤10,反之进入步骤D7;Step D6: Determine whether the calculation queue is empty, if it is empty, go to step 10, otherwise go to step D7;

步骤D7:取出队头节点a;Step D7: Take out the head node a;

步骤D8:遍历节点a的etx_r_map中的所有邻边,对于a的邻居b,计算Wa+Wab。其中Wa表示源节点到a的距离,Wab表示点a到点b的距离;Step D8: Traverse all neighboring edges in etx_r_map of node a, and for a's neighbor b, calculate Wa+Wab, where Wa represents the distance from the source node to a, and Wab represents the distance from point a to point b;

步骤D9:若Wb<Wa+Wab,则在node_map中更新Wb=Wa+Wab,同时判断节点b是否在队列中,若不在则将b放入计算队列中,并将node_map中节点b的布尔类型变量置为真,之后进入步骤D6;Step D9: If Wb<Wa+Wab, then update Wb=Wa+Wab in node_map, and determine whether node b is in the queue. If not, put b into the calculation queue, and set the Boolean type variable of node b in node_map to true, and then go to step D6;

步骤D10:查看源节点到目的节点d的权值是否比预设最大值还要大,若是则表明这两点距离超过了预设值,默认不存在最短路径,算法结束;反之进入步骤D11;Step D10: Check whether the weight from the source node to the destination node d is greater than the preset maximum value. If so, it indicates that the distance between the two points exceeds the preset value. By default, there is no shortest path and the algorithm ends. Otherwise, go to step D11;

步骤D11:在etx_r_map中更新目的节点d与目的节点上一跳节点d_pre的链路权值为之前的两倍;Step D11: Update the link weight between the destination node d and the previous hop node d_pre of the destination node in etx_r_map to twice the previous value;

步骤D12:将d_pre信息赋值给当前节点目的节点d;Step D12: assign d_pre information to the current node destination node d;

步骤D13:判断节点d是否为源节点,若不是则进入步骤D11,反之进入步骤D14;Step D13: Determine whether node d is a source node, if not, proceed to step D11, otherwise proceed to step D14;

步骤D14:将除源节点之外所有点的权值都设置为无穷大,同时循环变量自增,并返回步骤D4。Step D14: Set the weights of all points except the source node to infinity, increment the loop variable, and return to step D4.

进一步的,所述节点间链路权值大小定义为ETX,通过查询本发明修改后的拓扑表获取;修改后的拓扑表包含五个字段:Furthermore, the link weight between the nodes is defined as ETX, which is obtained by querying the modified topology table of the present invention; the modified topology table includes five fields:

T_dest_addr表示目的节点ip地址、T_last_addr表示到达目的节点的上一跳节点ip地址、T_seq表示这条拓扑消息的序列号、T_time表示这条拓扑消息的过期时间、ETX_R为新增加的字段,表示传输数期望的倒数;拓扑表项通过节点接收到的TC消息进行填充;T_dest_addr indicates the IP address of the destination node, T_last_addr indicates the IP address of the previous hop node to the destination node, T_seq indicates the sequence number of this topology message, T_time indicates the expiration time of this topology message, and ETX_R is a newly added field, indicating the inverse of the expected number of transmissions; the topology table entry is filled by the TC message received by the node;

所述TC消息包括ANSN表示拓扑消息序列号、Reserved表示保留字段、AdvertisedNeighbor MainAddress表示选择TC消息产生者节点作为MPR点的节点地址、ETX_R字段为新增加的字段,表示TC消息产生者与对应Advertised Neighbor MainAddress字段ip地址的节点的链路传输期望的倒数,通过查询TC消息产生者的一跳邻居节点获取。The TC message includes ANSN indicating the topology message sequence number, Reserved indicating a reserved field, AdvertisedNeighbor MainAddress indicating the node address of the TC message generator node selected as the MPR point, and the ETX_R field is a newly added field indicating the inverse of the link transmission expectation between the TC message generator and the node corresponding to the IP address of the Advertised Neighbor MainAddress field, which is obtained by querying the one-hop neighbor node of the TC message generator.

本发明的有益效果是:The beneficial effects of the present invention are:

本发明改善了网络的分组投递率、端到端时延以及吞吐量。The present invention improves the packet delivery rate, end-to-end delay and throughput of the network.

本发明提高网络的有效性与可靠性。The present invention improves the effectiveness and reliability of the network.

本发明在保证中心节点选出的MPR点能够覆盖自己的所有严格2跳邻居点的大前提下,综合考虑了节点连接度与链路质量,使得选出的MPR点转发消息更为可靠了,同时减少了最短路计算时间,提高了分组投递率。Under the premise of ensuring that the MPR point selected by the central node can cover all its strictly 2-hop neighboring points, the present invention comprehensively considers the node connectivity and link quality, making the selected MPR point forwarding message more reliable, while reducing the shortest path calculation time and improving the packet delivery rate.

附图说明BRIEF DESCRIPTION OF THE DRAWINGS

图1是本发明改进的MPR节点筛选流程图。FIG1 is a flowchart of an improved MPR node screening process of the present invention.

图2是本发明改进的多路SPFA算法流程图。FIG. 2 is a flow chart of the improved multi-path SPFA algorithm of the present invention.

图3是本发明的端到端时延与节点移动速度的关系图。FIG3 is a diagram showing the relationship between the end-to-end delay and the node moving speed of the present invention.

图4是本发明的分组投递率与节点移动速度的关系图。FIG. 4 is a diagram showing the relationship between the packet delivery rate and the node moving speed of the present invention.

图5是本发明的吞吐量与节点移动速度的关系图。FIG. 5 is a graph showing the relationship between throughput and node moving speed according to the present invention.

图6是本发明的新增的节点发送信息表。FIG. 6 is a table of newly added node sending information of the present invention.

图7是本发明的新增的节点接收信息表。FIG. 7 is a newly added node reception information table of the present invention.

图8是本发明的修改的节点链路信息表。FIG. 8 is a modified node link information table of the present invention.

图9是本发明的修改的节点邻居表。FIG. 9 is a modified node neighbor table of the present invention.

图10是本发明的修改后的HELLO分组结构。FIG. 10 is a modified HELLO packet structure of the present invention.

图11是本发明的修改后的TC分组结构。FIG. 11 is a modified TC packet structure of the present invention.

图12是本发明的修改后的网络拓扑表结构。FIG. 12 is a modified network topology table structure of the present invention.

图13是本发明的计算框图。FIG13 is a calculation block diagram of the present invention.

具体实施方式DETAILED DESCRIPTION

下面将结合本发明实施例中的附图对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。The technical solutions in the embodiments of the present invention will be described clearly and completely below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

节点连接度的计算可以通过遍历节点邻居表与二跳表得到,节点链路传输质量的计算需要修改节点HELLO消息与相应表格,查表得到。The calculation of node connectivity can be obtained by traversing the node neighbor table and the two-hop table. The calculation of node link transmission quality requires modifying the node HELLO message and the corresponding table and looking it up.

图6为本发明新增的节点发送信息表,用于存储节点发送给对应ip地址其他节点的消息数。每当节点给邻居节点发送HELLO消息时,就在此表对应的N_neighbor_main_addr中让Send_cnt数量增加1。Figure 6 is a newly added node sending information table of the present invention, which is used to store the number of messages sent by the node to other nodes with corresponding IP addresses. Whenever a node sends a HELLO message to a neighbor node, the number of Send_cnt in the corresponding N_neighbor_main_addr in this table is increased by 1.

图7为本发明新增的节点接收信息表,用于存储节点接收到对应ip地址其他节点发送的消息数。每当节点接收到邻居节点发送的HELLO消息时,就在此表对应的N_neighbor_main_addr中让Recv_cnt数量增加1。Figure 7 is a newly added node receiving information table of the present invention, which is used to store the number of messages received by the node from other nodes with the corresponding IP address. Whenever a node receives a HELLO message sent by a neighboring node, the Recv_cnt number in the corresponding N_neighbor_main_addr in this table is increased by 1.

图8为本发明修改的节点链路信息表,用于记录节点的不同链路之间的具体信息。新增了LQ_I字段与ETX_R字段,用于存储节点的正向链路质量与传输期望。LQ_I与ETX_R字段的计算方式如前所述。图9为本发明修改的节点邻居表,用于记录当前节点的邻居节点相关信息。新增了ETX_R字段,用于存储节点的传输期望。ETX_R字段与链路信息表中一致。FIG8 is a node link information table modified by the present invention, which is used to record specific information between different links of a node. Newly added LQ_I field and ETX_R field are used to store the forward link quality and transmission expectation of the node. The calculation method of LQ_I and ETX_R fields is as described above. FIG9 is a node neighbor table modified by the present invention, which is used to record information related to neighbor nodes of the current node. Newly added ETX_R field is used to store the transmission expectation of the node. The ETX_R field is consistent with that in the link information table.

图10描述了修改后的HELLO分组结构,新增了LQ_I字段,指的是源节点与相应邻居节点的正向链路传输质量。每当节点需要发送HELLO消息时,查询自己的链路信息表,找到对应ip地址的LQ_I并填充。Figure 10 depicts the modified HELLO packet structure, with a new LQ_I field added, which refers to the forward link transmission quality between the source node and the corresponding neighbor node. Whenever a node needs to send a HELLO message, it queries its own link information table, finds the LQ_I of the corresponding IP address and fills it in.

一种用于多无人机飞行自组网中的基于链路传输质量的MPOLSR协议优化方法,所述MPOLSR协议优化方法包括,A MPOLSR protocol optimization method based on link transmission quality for a multi-UAV flight ad hoc network, the MPOLSR protocol optimization method comprising:

基于多无人机之间进行信息交互的MPOLSR协议优化通过节点连接度与节点链路传输数期望进行MPR节点筛选,同时采用SPFA最短路算法,并以链路传输期望作为链路权值指标计算最短路,实现MPOLSR协议优化。The MPOLSR protocol optimization based on information interaction between multiple UAVs performs MPR node screening through node connectivity and the expected number of node link transmissions. At the same time, the SPFA shortest path algorithm is adopted, and the link transmission expectation is used as the link weight indicator to calculate the shortest path to achieve MPOLSR protocol optimization.

进一步的,所述MPR节点筛选为基于网络状态感知的MPR集选择,具体包括以下步骤,Further, the MPR node screening is an MPR set selection based on network state perception, which specifically includes the following steps:

步骤S1:首先建立一份备选的MPR节点列表M集,初始化为空集;Step S1: First, establish a candidate MPR node list M set, which is initialized as an empty set;

步骤S2:遍历节点1跳邻居集建立N1,遍历节点2跳邻居集建立严格2跳邻居集N2(严格2跳指的是源节点只能通过2跳到达的点,剔除了其中的1跳邻居节点);Step S2: traverse the node 1-hop neighbor set to establish N1, and traverse the node 2-hop neighbor set to establish a strict 2-hop neighbor set N2 (strict 2-hop means that the source node can only be reached through 2 hops, excluding the 1-hop neighbor nodes);

步骤S3:判断MPR集合中的点能否到达中心节点的所有严格2跳邻居点(也就是N2是否为空),若能覆盖则算法结束,否则进入步骤S4;Step S3: Determine whether the points in the MPR set can reach all strictly 2-hop neighboring points of the central node (that is, whether N2 is empty). If so, the algorithm ends, otherwise proceed to step S4;

步骤S4:计算N1中所有节点的权值,包括节点连接度与节点链路传输期望;Step S4: Calculate the weights of all nodes in N1, including node connectivity and node link transmission expectations;

步骤S5:遍历之前计算的所有点,找出权值最大的点;Step S5: traverse all the points calculated previously and find the point with the largest weight;

步骤S6:将权值最大的节点加入M集中;Step S6: Add the node with the largest weight to the M set;

步骤S7:在N1中移除权值最大的节点;Step S7: remove the node with the largest weight in N1;

步骤S8:在N2中移除能够通过最大的节点到达的点,之后返回步骤S3。Step S8: Remove the points in N2 that can be reached through the largest node, and then return to step S3.

进一步的,所述步骤S2中的1跳邻居集包含四个字段:Furthermore, the 1-hop neighbor set in step S2 includes four fields:

N_neighbor_main_addr代表源节点的邻居节点ip地址,N_status代表源节点与邻居节点的链路对称关系,N_willingness代表邻居节点转发源节点的消息意愿程序,ETX_R为新增加的字段,代表源节点与邻居节点的链路传输期望的倒数,通过查询修改后的节点链路信息表获取值。N_neighbor_main_addr represents the IP address of the neighbor node of the source node, N_status represents the link symmetry relationship between the source node and the neighbor node, N_willingness represents the willingness of the neighbor node to forward the message of the source node, and ETX_R is a newly added field, which represents the inverse of the link transmission expectation between the source node and the neighbor node. The value is obtained by querying the modified node link information table.

进一步的,所述节点链路信息表包含七个字段,具体为,L_local_iface_addr表示本地节点的ip地址,L_neighbor_iface_addr表示本地节点的邻居节点的ip地址,L_SYM_time记载了源节点与邻居节点对称链路的到期时间,L_ASYM_time记载了源节点与邻居节点非对称链路的到期时间,L_time表示该链路信息的到期时间,LQ_I字段与ETX_R字段为新增加的字段,分别表示节点与邻居节点的修改后的正向链路质量与传输数期望的倒数。LQ_I字段通过查询本发明新定义的节点发送信息表与节点接收信息表获取,ETX_R字段通过获取本发明修改后的HELLO消息包后获取。Further, the node link information table includes seven fields, specifically, L_local_iface_addr represents the IP address of the local node, L_neighbor_iface_addr represents the IP address of the neighbor node of the local node, L_SYM_time records the expiration time of the symmetric link between the source node and the neighbor node, L_ASYM_time records the expiration time of the asymmetric link between the source node and the neighbor node, L_time represents the expiration time of the link information, and the LQ_I field and the ETX_R field are newly added fields, which respectively represent the modified forward link quality and the expected reciprocal of the number of transmissions between the node and the neighbor node. The LQ_I field is obtained by querying the node sending information table and the node receiving information table newly defined in the present invention, and the ETX_R field is obtained by obtaining the modified HELLO message packet of the present invention.

进一步的,所述节点发送信息表与节点接收信息表具体为,通过查询本发明新定义的节点发送信息表获取Sij,节点发送信息表包括N_neighbor_main_addr字段记录邻居ip、Send_cnt字段记录发给对应ip节点的消息数;与此同时,通过查询本发明新定义的节点接收信息表获取Rij,节点接收信息表包括N_neighbor_main_addr字段记录邻居ip、Recv_cnt字段记录收到对应ip节点的消息数;由此可以计算出LQ_I值。Further, the node sending information table and the node receiving information table are specifically as follows: S ij is obtained by querying the node sending information table newly defined in the present invention, the node sending information table includes an N_neighbor_main_addr field to record the neighbor IP, and a Send_cnt field to record the number of messages sent to the corresponding IP node; at the same time, R ij is obtained by querying the node receiving information table newly defined in the present invention, the node receiving information table includes an N_neighbor_main_addr field to record the neighbor IP, and a Recv_cnt field to record the number of messages received from the corresponding IP node; thereby, the LQ_I value can be calculated.

进一步的,所述修改后的HELLO消息包包括Reserved第一保留字段、Htime发送间隔、Willingness数据转发意愿、Linkcode链路编码信息、Reserved第二保留字段、LinkMessage Size链路状态的大小信息、NeighborInterfaceAddress邻居节点接口地址与新增加的LQ_I字段,表示产生该HELLO消息的源节点与接收到此HELLO消息的节点的LQ_I大小,源节点收到消息后即可得知自己与对应邻居节点的NLQ_I大小,从而通过前述方案计算出ETX_R,将结果填充入节点链路信息表中。Furthermore, the modified HELLO message packet includes the Reserved first reserved field, Htime sending interval, Willingness data forwarding intention, Linkcode link coding information, Reserved second reserved field, LinkMessage Size link state size information, NeighborInterfaceAddress neighbor node interface address and the newly added LQ_I field, which indicates the LQ_I size of the source node that generates the HELLO message and the node that receives the HELLO message. After receiving the message, the source node can know the NLQ_I size of itself and the corresponding neighbor node, thereby calculating ETX_R through the aforementioned scheme and filling the result into the node link information table.

进一步的,所述步骤S4中MPR节点权值指标的计算具体为,指标计算公式为:Furthermore, the calculation of the MPR node weight index in step S4 is specifically as follows:

Figure SMS_7
Figure SMS_7

其中ETX_R指的是链路传输期望ETX的倒数,D指的是节点连接度;α和β是权重指标,且满足α+β=1;权重指标具体的设置可以根据实际的不同应用场景作出对应调整。Among them, ETX_R refers to the inverse of the expected ETX of link transmission, and D refers to the node connectivity; α and β are weight indicators, and α+β=1; the specific setting of the weight indicator can be adjusted accordingly according to different actual application scenarios.

进一步的,所述节点链路传输数期望计算具体为,Furthermore, the expected calculation of the node link transmission number is specifically as follows:

所述链路传输期望为,The link transmission expectation is:

Figure SMS_8
Figure SMS_8

Figure SMS_9
Figure SMS_9

Figure SMS_10
Figure SMS_10

Figure SMS_11
Figure SMS_11

Figure SMS_12
Figure SMS_12

其中Rji指的是节点j收到节点i的HELLO消息数,Sij指的是节点i发送给节点j的HELLO消息数;Rij指的是节点i收到节点j的HELLO消息数,Sji指的是节点j发送给节点i的HELLO消息数;LQ代表正向链路质量,NLQ代表反向链路质量,LQ_I代表修改后的正向链路质量,NLQ_I代表修改后的反向链路质量。对于节点i来说,i可以获得的关于HELLO消息的数目只有“i发送给j的HELLO消息数”和“i收到j发送的HELLO消息数”,仅凭此数据是无法计算LQ、NLQ的,进一步无法计算ETX。因此,本发明对传输数期望的两个参量进行了一定程度的改变,将LQ和NLQ的分子互换,变成新的两个参量LQ_I和NLQ_I同时不影响传输数期望的计算。Wherein R ji refers to the number of HELLO messages received by node j from node i, and S ij refers to the number of HELLO messages sent by node i to node j; R ij refers to the number of HELLO messages received by node i from node j, and S ji refers to the number of HELLO messages sent by node j to node i; LQ represents the forward link quality, NLQ represents the reverse link quality, LQ_I represents the modified forward link quality, and NLQ_I represents the modified reverse link quality. For node i, the number of HELLO messages that i can obtain is only "the number of HELLO messages sent by i to j" and "the number of HELLO messages sent by i to j". It is impossible to calculate LQ and NLQ based on this data alone, and further it is impossible to calculate ETX. Therefore, the present invention makes a certain degree of change to the two parameters of the expected transmission number, and the numerators of LQ and NLQ are exchanged to become two new parameters LQ_I and NLQ_I, while not affecting the calculation of the expected transmission number.

MPOLSR原协议计算多条最短路采用的为Dijkstra算法,复杂度较高,且未区分不同链路的权值;本发明以链路传输期望作为各个链路权值进行区分,同时采用复杂度更低的SPFA算法进行实现。通过修改TC消息使每个节点都能感知到网络中所有链路的边权值,同时修改网络拓扑表,查表计算最短路径。The original MPOLSR protocol uses the Dijkstra algorithm to calculate multiple shortest paths, which is highly complex and does not distinguish the weights of different links. The present invention uses link transmission expectations as the weights of each link to distinguish, and uses the SPFA algorithm with lower complexity to implement it. By modifying the TC message, each node can perceive the edge weights of all links in the network, and at the same time modify the network topology table to look up the table to calculate the shortest path.

图11描述了修改后的TC分组结构,新增了ETX_R字段,指的是源节点与相应邻居节点的链路传输数期望的倒数。FIG11 illustrates the modified TC packet structure, with a newly added ETX_R field, which refers to the expected inverse of the link transmission number between the source node and the corresponding neighbor node.

图12描述了修改后的网络拓扑表结构,新增了ETX_R字段,记载了T_last_addr与T_dest_addr链路之间的传输期望的倒数。FIG12 illustrates the modified network topology table structure, with a newly added ETX_R field recording the expected inverse of the transmission between the T_last_addr and T_dest_addr links.

基于链路传输期望的最短路算法流程图如图2所示,包括如下步骤:The shortest path algorithm based on link transmission expectation is shown in FIG2 , and includes the following steps:

进一步的,所述最短路为链路传输数期望的最短路优化计算,具体为,Furthermore, the shortest path is an optimization calculation of the expected shortest path of the link transmission number, specifically,

步骤D1:自定义一个node_map结构体,结构体中含有四个成员变量,分别是每个节点的地址、从源节点到该节点的上一跳地址、中心节点到该点的距离,以及一个布尔类型变量,表示该节点是否在计算队列中,通过遍历源节点的1跳邻居集、2跳邻居集与拓扑表进行填充;Step D1: Customize a node_map structure, which contains four member variables, namely the address of each node, the previous hop address from the source node to the node, the distance from the central node to the point, and a Boolean variable, indicating whether the node is in the calculation queue, which is filled by traversing the 1-hop neighbor set, 2-hop neighbor set and topology table of the source node;

步骤D2:建立一个etx_r_map二维数组用以存储每个节点到其他任意节点的链路权值大小,通过遍历源节点的1跳邻居集、2跳邻居集与拓扑表进行填充。链路权值的大小就是节点的链路传输期望;Step D2: Create a two-dimensional array etx_r_map to store the link weights from each node to any other node, and fill it by traversing the 1-hop neighbor set, 2-hop neighbor set and topology table of the source node. The link weight is the link transmission expectation of the node;

步骤D3:设置循环变量,初始化为0;Step D3: Set the loop variable and initialize it to 0;

步骤D4:判断循环变量的值是否达到了设定的多径路由数量,若达到则结束算法,反之进入步骤D5;Step D4: Determine whether the value of the loop variable has reached the set number of multipath routes. If so, the algorithm ends. Otherwise, it goes to step D5.

步骤D5:设定一个计算队列用于存储待更新的节点,同时将当前源节点放入计算队列中;Step D5: Set a calculation queue to store the nodes to be updated, and put the current source node into the calculation queue;

步骤D6:判断计算队列是否为空,若为空则进入步骤10,反之进入步骤D7;Step D6: Determine whether the calculation queue is empty, if it is empty, go to step 10, otherwise go to step D7;

步骤D7:取出队头节点a;Step D7: Take out the head node a;

步骤D8:遍历节点a的etx_r_map中的所有邻边,对于a的邻居b,计算Wa+Wab。其中Wa表示源节点到a的距离,Wab表示点a到点b的距离;Step D8: Traverse all neighboring edges in etx_r_map of node a, and for a's neighbor b, calculate Wa+Wab, where Wa represents the distance from the source node to a, and Wab represents the distance from point a to point b;

步骤D9:若Wb<Wa+Wab,则在node_map中更新Wb=Wa+Wab,同时判断节点b是否在队列中,若不在则将b放入计算队列中,并将node_map中节点b的布尔类型变量置为真,其中Wb表示源节点到b的距离,之后进入步骤D6;Step D9: If Wb<Wa+Wab, then update Wb=Wa+Wab in node_map, and determine whether node b is in the queue. If not, put b into the calculation queue, and set the Boolean type variable of node b in node_map to true, where Wb represents the distance from the source node to b, and then go to step D6;

步骤D10:查看源节点到目的节点d的权值是否比预设最大值还要大,若是则表明这两点距离超过了预设值,默认不存在最短路径,算法结束;反之进入步骤D11;Step D10: Check whether the weight from the source node to the destination node d is greater than the preset maximum value. If so, it indicates that the distance between the two points exceeds the preset value. By default, there is no shortest path and the algorithm ends. Otherwise, go to step D11;

步骤D11:在etx_r_map中更新目的节点d与目的节点上一跳节点d_pre的链路权值为之前的两倍,这是为了防止计算重复路径;Step D11: Update the link weight between the destination node d and the previous hop node d_pre in etx_r_map to twice the previous value, in order to prevent the calculation of duplicate paths;

步骤D12:将d_pre信息赋值给当前节点目的节点d;Step D12: assign d_pre information to the current node destination node d;

步骤D13:判断节点d是否为源节点,若不是则进入步骤D11,反之进入步骤D14;Step D13: Determine whether node d is a source node, if not, proceed to step D11, otherwise proceed to step D14;

步骤D14:将除源节点之外所有点的权值都设置为无穷大,同时循环变量自增,并返回步骤D4。Step D14: Set the weights of all points except the source node to infinity, increment the loop variable, and return to step D4.

进一步的,所述节点间链路权值大小定义为ETX,通过查询本发明修改后的拓扑表获取;修改后的拓扑表包含五个字段:Furthermore, the link weight between the nodes is defined as ETX, which is obtained by querying the modified topology table of the present invention; the modified topology table includes five fields:

T_dest_addr表示目的节点ip地址、T_last_addr表示到达目的节点的上一跳节点ip地址、T_seq表示这条拓扑消息的序列号、T_time表示这条拓扑消息的过期时间、ETX_R为新增加的字段,表示传输数期望的倒数;拓扑表项通过节点接收到的TC消息进行填充;T_dest_addr indicates the IP address of the destination node, T_last_addr indicates the IP address of the previous hop node to the destination node, T_seq indicates the sequence number of this topology message, T_time indicates the expiration time of this topology message, and ETX_R is a newly added field, indicating the inverse of the expected number of transmissions; the topology table entry is filled by the TC message received by the node;

所述TC消息包括ANSN表示拓扑消息序列号、Reserved表示保留字段、AdvertisedNeighbor MainAddress表示选择TC消息产生者节点作为MPR点的节点地址、ETX_R字段为新增加的字段,表示TC消息产生者与对应Advertised Neighbor MainAddress字段ip地址的节点的链路传输期望的倒数,通过查询TC消息产生者的一跳邻居节点获取。The TC message includes ANSN indicating the topology message sequence number, Reserved indicating a reserved field, AdvertisedNeighbor MainAddress indicating the node address of the TC message generator node selected as the MPR point, and the ETX_R field is a newly added field indicating the inverse of the link transmission expectation between the TC message generator and the node corresponding to the IP address of the Advertised Neighbor MainAddress field, which is obtained by querying the one-hop neighbor node of the TC message generator.

将采用本次发明优化后的MPOLSR协议命名为AMPOLSR协议。为验证本发明所取得的效果,基于网络仿真软件NS2,实现了不同场景下MPOLSR、采用了Dijkstra算法计算最短路的AMPOLSR-Dijkstra与采用了SPFA算法优化后的AMPOLSR-SPFA这三种协议的仿真,对比了这三种路由协议的主要性能指标,包括分组投递率、平均端到端时延与网络吞吐量。本发明所采用的主要仿真参数见下表,仿真模型采用的是RPGM移动模型。The MPOLSR protocol optimized by this invention is named as AMPOLSR protocol. To verify the effect achieved by the present invention, based on the network simulation software NS2, the simulation of MPOLSR in different scenarios, AMPOLSR-Dijkstra using Dijkstra algorithm to calculate the shortest path, and AMPOLSR-SPFA optimized by SPFA algorithm are implemented, and the main performance indicators of these three routing protocols are compared, including packet delivery rate, average end-to-end delay and network throughput. The main simulation parameters used by the present invention are shown in the table below, and the simulation model uses the RPGM mobile model.

Figure SMS_13
Figure SMS_13

Figure SMS_14
Figure SMS_14

图3描述了不同协议下端到端时延与节点移动速度的关系图。由图3可知:Figure 3 describes the relationship between end-to-end delay and node movement speed under different protocols. It can be seen from Figure 3 that:

平均端到端时延随着节点移动速度的增加呈现上升趋势,这是因为在节点移动速度较低时,网络拓扑相对静止,路由稳定,因此时延较低;但随着节点最大移动速度的增加,网络拓扑结构变化进一步加快,拓扑表更新不及时导致先前计算的最短路径可能并不是当前网络状况下的最短路径,使得信息传播时延加大。同时也可以看到修改后的协议比原协议的时延更低,SPFA算法的时延小于Dijkstra算法的时延,说明改进后的算法找到的最短路更佳,所以信息传递时延降低了,且SPFA算法因为复杂度低故时延低于Dijkstra算法。The average end-to-end delay shows an upward trend as the node movement speed increases. This is because when the node movement speed is low, the network topology is relatively static and the routing is stable, so the delay is low; but as the maximum node movement speed increases, the network topology structure changes further, and the topology table is not updated in time, resulting in the shortest path previously calculated may not be the shortest path under the current network conditions, which increases the information propagation delay. At the same time, it can be seen that the modified protocol has a lower delay than the original protocol, and the delay of the SPFA algorithm is less than the delay of the Dijkstra algorithm, indicating that the improved algorithm finds a better shortest path, so the information transmission delay is reduced, and the SPFA algorithm has a lower delay than the Dijkstra algorithm because of its low complexity.

图4描述了不同协议下分组投递率与节点移动速度的关系图。由图4可知:Figure 4 describes the relationship between packet delivery rate and node movement speed under different protocols. It can be seen from Figure 4 that:

分组投递率随着节点移动速度的增加呈现下降趋势,这是因为在节点移动速度较低时,网络拓扑相对静止,信息丢包率较少,但随着节点最大移动速度的增加,网络拓扑结构变化进一步加快,链路稳定性下降,因此分组投递率整体上呈现下降趋势。同时也可以看到修改后的协议比原协议的分组投递率要大,SPFA算法的投递率大于Dijkstra算法的投递率,说明该算法确实优化了分组投递率,且SPFA算法效果更佳,降低了更多的丢包率,提高了网络传输可靠性。The packet delivery rate shows a downward trend as the node movement speed increases. This is because when the node movement speed is low, the network topology is relatively static and the information packet loss rate is low. However, as the maximum node movement speed increases, the network topology structure changes further and the link stability decreases. Therefore, the packet delivery rate shows a downward trend overall. At the same time, it can be seen that the packet delivery rate of the modified protocol is larger than that of the original protocol. The delivery rate of the SPFA algorithm is greater than that of the Dijkstra algorithm, which shows that the algorithm does optimize the packet delivery rate, and the SPFA algorithm is more effective, reducing more packet loss rates and improving network transmission reliability.

图5描述了不同协议下吞吐量与节点移动速度的关系图。由图5可知:Figure 5 describes the relationship between throughput and node movement speed under different protocols. It can be seen from Figure 5 that:

网络吞吐量随着节点移动速度的增加大体呈现下降趋势,这是因为在节点移动速度较低时,网络拓扑相对静止,节点之间链路稳定,网络控制开销小,能够发送更多的实际消息包,因此网络吞吐量大。但随着节点最大移动速度的增加,网络拓扑结构变化进一步加快,链路不稳定性变大,信息丢包率加大,导致网络吞吐量降低。这一趋势与前述分组投递率一致。同时也可以看到修改后的协议比原协议的吞吐量要大,SPFA算法的吞吐量大于Dijkstra算法的吞吐量,说明传输了更多信息流。The network throughput generally shows a downward trend as the node movement speed increases. This is because when the node movement speed is low, the network topology is relatively static, the links between nodes are stable, the network control overhead is small, and more actual message packets can be sent, so the network throughput is large. However, as the maximum node movement speed increases, the network topology changes further accelerate, the link instability becomes greater, and the information packet loss rate increases, resulting in a decrease in network throughput. This trend is consistent with the aforementioned packet delivery rate. At the same time, it can be seen that the throughput of the modified protocol is larger than that of the original protocol. The throughput of the SPFA algorithm is greater than the throughput of the Dijkstra algorithm, indicating that more information flows are transmitted.

仿真结果表明:改进的MPOLSR路由协议通过优化MPR选择方法与最短路算法,不仅提高了网络的分组投递率与吞吐量,同时还降低了网络的端到端时延,即有效性与可靠性均得到了提升。因此,基于改进后的MPOLSR路由协议能够更好的应用于无人机飞行自组网中,满足任务需求。The simulation results show that the improved MPOLSR routing protocol not only improves the network's packet delivery rate and throughput by optimizing the MPR selection method and the shortest path algorithm, but also reduces the network's end-to-end delay, that is, the effectiveness and reliability are improved. Therefore, the improved MPOLSR routing protocol can be better applied to UAV flight ad hoc networks to meet mission requirements.

本发明采取了一种基于链路传输期望的最短路优化算法,其核心是采用SPFA算法,该算法的平均复杂度为O(km),其中k是一个很小的常数,m代表图论中边的数量。在最坏情况下复杂度为O(nm),其中n代表图论中点的数量。而原协议最短路算法采用Dijkstra算法,受限于协议的节点拓扑表结构,其执行复杂度为O(nm),故一般而言SPFA算法复杂度要较Dijkstra算法小,且该方法以链路传输期望区分了不同边的权值大小,执行此算法能够减小端到端时延,有效提高分组投递率与网络吞吐量。The present invention adopts a shortest path optimization algorithm based on link transmission expectation, the core of which is to adopt the SPFA algorithm, the average complexity of which is O(km), where k is a very small constant and m represents the number of edges in graph theory. In the worst case, the complexity is O(nm), where n represents the number of points in graph theory. The original protocol shortest path algorithm adopts the Dijkstra algorithm, which is limited by the node topology table structure of the protocol and has an execution complexity of O(nm). Therefore, in general, the complexity of the SPFA algorithm is smaller than that of the Dijkstra algorithm, and this method distinguishes the weights of different edges based on link transmission expectations. Executing this algorithm can reduce end-to-end delay and effectively improve packet delivery rate and network throughput.

以上所述仅为本发明创造的较佳实施例而已,并不用以限制本发明创造,凡在本发明创造的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明创造的保护范围之内。The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principle of the present invention should be included in the protection scope of the present invention.

Claims (10)

1.一种用于多无人机飞行自组网中的基于链路传输质量的MPOLSR协议优化方法,其特征在于,所述MPOLSR协议优化方法包括,1. a kind of MPOLSR protocol optimization method based on link transmission quality in many unmanned aerial vehicle flight ad hoc networks, it is characterized in that, described MPOLSR protocol optimization method comprises, 基于多无人机之间进行信息交互的MPOLSR协议优化通过节点连接度与节点链路传输数期望进行MPR节点筛选,同时选取最短路并以链路传输期望作为MPR节点权值指标,实现MPOLSR协议优化。The MPOLSR protocol optimization based on information interaction between multiple UAVs screens MPR nodes through node connectivity and node link transmission expectations, and at the same time selects the shortest path and uses link transmission expectations as the MPR node weight index to implement the MPOLSR protocol optimization. 2.根据权利要求1所述一种用于多无人机飞行自组网中的基于链路传输质量的MPOLSR协议优化方法,其特征在于,所述MPR节点筛选为基于网络状态感知的MPR集选择,具体包括以下步骤,2. according to claim 1, a kind of MPOLSR protocol optimization method based on link transmission quality for multi-unmanned aerial vehicle flight ad hoc network, it is characterized in that, described MPR node is screened as the MPR set based on network state awareness selection, specifically including the following steps, 步骤S1:首先建立一份备选的MPR节点列表M集,初始化为空集;Step S1: First, create a candidate MPR node list M set, which is initialized as an empty set; 步骤S2:遍历节点1跳邻居集建立N1,遍历节点2跳邻居集建立严格2跳邻居集N2;Step S2: Traversing the node 1-hop neighbor set to establish N1, traversing the node 2-hop neighbor set to establish a strict 2-hop neighbor set N2; 步骤S3:判断MPR集合中的点能否到达中心节点的所有严格2跳邻居点,若能覆盖则算法结束,否则进入步骤S4;Step S3: Determine whether the points in the MPR set can reach all the strict 2-hop neighbor points of the central node. If they can be covered, the algorithm ends, otherwise, enter step S4; 步骤S4:计算N1中所有节点的权值,包括节点连接度与节点链路传输期望;Step S4: Calculate the weights of all nodes in N1, including node connectivity and node link transmission expectations; 步骤S5:遍历之前计算的所有点,找出权值最大的点;Step S5: Traversing all the previously calculated points to find out the point with the largest weight; 步骤S6:将权值最大的节点加入M集中;Step S6: Add the node with the largest weight into the M set; 步骤S7:在N1中移除权值最大的节点;Step S7: remove the node with the largest weight in N1; 步骤S8:在N2中移除能够通过最大的节点到达的点,之后返回步骤S3。Step S8: Remove the point that can be reached through the largest node in N2, and then return to step S3. 3.根据权利要求2所述一种用于多无人机飞行自组网中的基于链路传输质量的MPOLSR协议优化方法,其特征在于,所述步骤S2中的1跳邻居集包含四个字段:N_neighbor_main_addr代表源节点的邻居节点ip地址,N_status代表源节点与邻居节点的链路对称关系,N_willingness代表邻居节点转发源节点的消息意愿程序,ETX_R为新增加的字段,代表源节点与邻居节点的链路传输期望的倒数,通过查询修改后的节点链路信息表获取值。3. according to claim 2, a kind of MPOLSR protocol optimization method based on link transmission quality for multi-unmanned aerial vehicle flight ad hoc network, it is characterized in that, the 1-hop neighbor set in the described step S2 comprises four Field: N_neighbor_main_addr represents the ip address of the neighbor node of the source node, N_status represents the link symmetry relationship between the source node and the neighbor node, N_willingness represents the message willingness program of the neighbor node to forward the source node, ETX_R is a newly added field, representing the source node and the neighbor node The reciprocal of the link transmission expectation of , and obtain the value by querying the modified node link information table. 4.根据权利要求3所述一种用于多无人机飞行自组网中的基于链路传输质量的MPOLSR协议优化方法,其特征在于,所述节点链路信息表包含七个字段,具体为,L_local_iface_addr表示本地节点的ip地址,L_neighbor_iface_addr表示本地节点的邻居节点的ip地址,L_SYM_time记载了源节点与邻居节点对称链路的到期时间,L_ASYM_time记载了源节点与邻居节点非对称链路的到期时间,L_time表示该链路信息的到期时间,LQ_I字段与ETX_R字段为新增加的字段,分别表示节点与邻居节点的修改后的正向链路质量与传输数期望的倒数。LQ_I字段通过查询本发明新定义的节点发送信息表与节点接收信息表获取,ETX_R字段通过获取本发明修改后的HELLO消息包后获取。4. according to claim 3 a kind of MPOLSR protocol optimization method based on link transmission quality in the multi-unmanned aerial vehicle flight ad hoc network, it is characterized in that, described node link information table comprises seven fields, specifically L_local_iface_addr represents the ip address of the local node, L_neighbor_iface_addr represents the ip address of the neighbor node of the local node, L_SYM_time records the expiration time of the symmetric link between the source node and the neighbor node, L_ASYM_time records the time of the asymmetric link between the source node and the neighbor node The expiration time, L_time indicates the expiration time of the link information, the LQ_I field and the ETX_R field are newly added fields, which respectively represent the reciprocal of the modified forward link quality and transmission number expectation of the node and the neighbor node. The LQ_I field is obtained by querying the node transmission information table and the node reception information table newly defined in the present invention, and the ETX_R field is obtained by obtaining the modified HELLO message package of the present invention. 5.根据权利要求4所述一种用于多无人机飞行自组网中的基于链路传输质量的MPOLSR协议优化方法,其特征在于,所述节点发送信息表与节点接收信息表具体为,通过查询本发明新定义的节点发送信息表获取Sij,节点发送信息表包括N_neighbor_main_addr字段记录邻居ip、Send_cnt字段记录发给对应ip节点的消息数;与此同时,通过查询本发明新定义的节点接收信息表获取Rij,节点接收信息表包括N_neighbor_main_addr字段记录邻居ip、Recv_cnt字段记录收到对应ip节点的消息数;由此可以计算出LQ_I值。5. according to claim 4, a kind of MPOLSR protocol optimization method based on link transmission quality for multi-unmanned aerial vehicle flying ad hoc network, it is characterized in that, described node sends information table and node receives information table specifically as , obtain S ij by querying the node sending information table newly defined in the present invention, the node sending information table includes the N_neighbor_main_addr field record neighbor ip, the Send_cnt field record sending to the number of messages corresponding to the ip node; meanwhile, by querying the new definition of the present invention The node receiving information table obtains R ij , and the node receiving information table includes the N_neighbor_main_addr field to record the neighbor ip, and the Recv_cnt field to record the number of messages received from the corresponding ip node; thus the LQ_I value can be calculated. 6.根据权利要求4所述一种用于多无人机飞行自组网中的基于链路传输质量的MPOLSR协议优化方法,其特征在于,所述修改后的HELLO消息包包括Reserved第一保留字段、Htime发送间隔、Willingness数据转发意愿、Linkcode链路编码信息、Reserved第二保留字段、LinkMessage Size链路状态的大小信息、NeighborInterfaceAddress邻居节点接口地址与新增加的LQ_I字段,表示产生该HELLO消息的源节点与接收到此HELLO消息的节点的LQ_I大小,源节点收到消息后即可得知自己与对应邻居节点的NLQ_I大小,从而通过前述方案计算出ETX_R,将结果填充入节点链路信息表中。6. according to claim 4, a kind of MPOLSR protocol optimization method based on link transmission quality for multi-unmanned aerial vehicle flight ad hoc network, it is characterized in that, the HELLO message package after the modification includes Reserved first reserved Field, Htime sending interval, Willingness data forwarding willingness, Linkcode link encoding information, Reserved second reserved field, LinkMessage Size link state size information, NeighborInterfaceAddress neighbor node interface address and the newly added LQ_I field, indicating the source of the HELLO message The LQ_I size of the source node and the node that received the HELLO message. After receiving the message, the source node can know the NLQ_I size of itself and the corresponding neighbor node, so as to calculate ETX_R through the above scheme, and fill the result into the node link information table middle. 7.根据权利要求2所述一种用于多无人机飞行自组网中的基于链路传输质量的MPOLSR协议优化方法,其特征在于,所述步骤S4中MPR节点权值指标的计算具体为,指标计算公式为:7. according to claim 2, a kind of MPOLSR protocol optimization method based on link transmission quality for multi-unmanned aerial vehicle flight ad hoc network, it is characterized in that, the calculation of MPR node weight index in the described step S4 is specific The index calculation formula is:
Figure FDA0004153376030000021
Figure FDA0004153376030000021
其中ETX_R指的是链路传输期望ETX的倒数,D指的是节点连接度;α和β是权重指标,且满足α+β=1。Where ETX_R refers to the reciprocal of link transmission expectation ETX, D refers to node connectivity; α and β are weight indexes, and satisfy α+β=1.
8.根据权利要求7所述一种用于多无人机飞行自组网中的基于链路传输质量的MPOLSR协议优化方法,其特征在于,所述节点链路传输数期望计算具体为,8. according to claim 7, a kind of MPOLSR protocol optimization method based on link transmission quality for multi-unmanned aerial vehicle flight ad hoc network, it is characterized in that, the expected calculation of the number of node link transmissions is specifically, 所述链路传输期望为,The link transmission expectation is,
Figure FDA0004153376030000031
Figure FDA0004153376030000031
Figure FDA0004153376030000032
Figure FDA0004153376030000032
Figure FDA0004153376030000033
Figure FDA0004153376030000033
Figure FDA0004153376030000034
Figure FDA0004153376030000034
Figure FDA0004153376030000035
Figure FDA0004153376030000035
其中Rji指的是节点j收到节点i的HELLO消息数,Sij指的是节点i发送给节点j的HELLO消息数;Rij指的是节点i收到节点j的HELLO消息数,Sji指的是节点j发送给节点i的HELLO消息数;LQ代表正向链路质量,NLQ代表反向链路质量,LQ_I代表修改后的正向链路质量,NLQ_I代表修改后的反向链路质量。Among them, R ji refers to the number of HELLO messages that node j receives from node i, S ij refers to the number of HELLO messages that node i sends to node j; R ij refers to the number of HELLO messages that node i receives from node j, S ji refers to the number of HELLO messages sent by node j to node i; LQ stands for forward link quality, NLQ stands for reverse link quality, LQ_I stands for modified forward link quality, NLQ_I stands for modified reverse link quality road quality.
9.根据权利要求1所述一种用于多无人机飞行自组网中的基于链路传输质量的MPOLSR协议优化方法,其特征在于,所述最短路为链路传输数期望的最短路优化计算,具体为,9. according to claim 1, a kind of MPOLSR protocol optimization method based on link transmission quality for multi-unmanned aerial vehicle flight ad hoc network, it is characterized in that, described shortest path is the expected shortest path of link transmission number The optimization calculation, specifically, 步骤D1:自定义一个node_map结构体,结构体中含有四个成员变量,分别是每个节点的地址、从源节点到该节点的上一跳地址、中心节点到该点的距离,以及一个布尔类型变量,表示该节点是否在计算队列中,通过遍历源节点的1跳邻居集、2跳邻居集与拓扑表进行填充;Step D1: Customize a node_map structure, which contains four member variables, which are the address of each node, the address of the last hop from the source node to the node, the distance from the center node to the point, and a Boolean The type variable indicates whether the node is in the calculation queue, and is filled by traversing the 1-hop neighbor set, 2-hop neighbor set and topology table of the source node; 步骤D2:建立一个etx_r_map二维数组用以存储每个节点到其他任意节点的链路权值大小,通过遍历源节点的1跳邻居集、2跳邻居集与拓扑表进行填充。链路权值的大小就是节点的链路传输期望;Step D2: Create a etx_r_map two-dimensional array to store the link weight value from each node to any other node, and fill it by traversing the 1-hop neighbor set, 2-hop neighbor set and topology table of the source node. The size of the link weight is the link transmission expectation of the node; 步骤D3:设置循环变量,初始化为0;Step D3: Set the loop variable and initialize it to 0; 步骤D4:判断循环变量的值是否达到了设定的多径路由数量,若达到则结束算法,反之进入步骤D5;Step D4: Determine whether the value of the loop variable has reached the set number of multipath routes, and if so, end the algorithm, otherwise enter step D5; 步骤D5:设定一个计算队列用于存储待更新的节点,同时将当前源节点放入计算队列中;Step D5: Set a calculation queue to store the nodes to be updated, and put the current source node into the calculation queue; 步骤D6:判断计算队列是否为空,若为空则进入步骤10,反之进入步骤D7;Step D6: Determine whether the calculation queue is empty, if it is empty, go to step 10, otherwise go to step D7; 步骤D7:取出队头节点a;Step D7: Take out the queue head node a; 步骤D8:遍历节点a的etx_r_map中的所有邻边,对于a的邻居b,计算Wa+Wab。其中Wa表示源节点到a的距离,Wab表示点a到点b的距离;Step D8: Traverse all adjacent edges in the etx_r_map of node a, and calculate Wa+Wab for neighbor b of a. Where Wa represents the distance from the source node to a, and Wab represents the distance from point a to point b; 步骤D9:若Wb<Wa+Wab,则在node_map中更新Wb=Wa+Wab,同时判断节点b是否在队列中,若不在则将b放入计算队列中,并将node_map中节点b的布尔类型变量置为真,之后进入步骤D6;Step D9: If Wb<Wa+Wab, then update Wb=Wa+Wab in node_map, and judge whether node b is in the queue at the same time, if not, put b into the calculation queue, and set the Boolean type of node b in node_map The variable is set to true, and then enter step D6; 步骤D10:查看源节点到目的节点d的权值是否比预设最大值还要大,若是则表明这两点距离超过了预设值,默认不存在最短路径,算法结束;反之进入步骤D11;Step D10: Check whether the weight value from the source node to the destination node d is greater than the preset maximum value, if so, it indicates that the distance between these two points exceeds the preset value, there is no shortest path by default, and the algorithm ends; otherwise, enter step D11; 步骤D11:在etx_r_map中更新目的节点d与目的节点上一跳节点d_pre的链路权值为之前的两倍;Step D11: In etx_r_map, update the link weight between the destination node d and the last hop node d_pre of the destination node to twice the previous value; 步骤D12:将d_pre信息赋值给当前节点目的节点d;Step D12: Assign d_pre information to the destination node d of the current node; 步骤D13:判断节点d是否为源节点,若不是则进入步骤D11,反之进入步骤D14;Step D13: Determine whether node d is a source node, if not, go to step D11, otherwise go to step D14; 步骤D14:将除源节点之外所有点的权值都设置为无穷大,同时循环变量自增,并返回步骤D4。Step D14: Set the weights of all points except the source node to infinity, meanwhile the loop variable is incremented, and return to step D4. 10.根据权利要求9所述一种用于多无人机飞行自组网中的基于链路传输质量的MPOLSR协议优化方法,其特征在于,所述节点间链路权值大小定义为ETX,通过查询本发明修改后的拓扑表获取;修改后的拓扑表包含五个字段:10. according to claim 9, a kind of MPOLSR protocol optimization method based on link transmission quality for multi-unmanned aerial vehicle flight ad hoc network, is characterized in that, the link weight value size between described nodes is defined as ETX, Obtain by querying the modified topology table of the present invention; the modified topology table includes five fields: T_dest_addr表示目的节点ip地址、T_last_addr表示到达目的节点的上一跳节点ip地址、T_seq表示这条拓扑消息的序列号、T_time表示这条拓扑消息的过期时间、ETX_R为新增加的字段,表示传输数期望的倒数;拓扑表项通过节点接收到的TC消息进行填充;T_dest_addr indicates the ip address of the destination node, T_last_addr indicates the ip address of the last hop node to reach the destination node, T_seq indicates the sequence number of this topology message, T_time indicates the expiration time of this topology message, ETX_R is a newly added field, indicating the number of transmissions The expected reciprocal; the topology entry is filled by the TC message received by the node; 所述TC消息包括ANSN表示拓扑消息序列号、Reserved表示保留字段、AdvertisedNeighbor MainAddress表示选择TC消息产生者节点作为MPR点的节点地址、ETX_R字段为新增加的字段,表示TC消息产生者与对应Advertised Neighbor MainAddress字段ip地址的节点的链路传输期望的倒数,通过查询TC消息产生者的一跳邻居节点获取。The TC message includes ANSN representing the sequence number of the topology message, Reserved representing the reserved field, AdvertisedNeighbor MainAddress representing the node address for selecting the TC message generator node as the MPR point, and the ETX_R field being a newly added field, indicating that the TC message generator and the corresponding Advertised Neighbor The reciprocal of the link transmission expectation of the node with the ip address in the MainAddress field is obtained by querying the one-hop neighbor node of the TC message generator.
CN202310326422.7A 2023-03-30 2023-03-30 An MPOLSR protocol optimization method based on link transmission quality for multi-UAV flight ad hoc network Pending CN116261200A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310326422.7A CN116261200A (en) 2023-03-30 2023-03-30 An MPOLSR protocol optimization method based on link transmission quality for multi-UAV flight ad hoc network

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310326422.7A CN116261200A (en) 2023-03-30 2023-03-30 An MPOLSR protocol optimization method based on link transmission quality for multi-UAV flight ad hoc network

Publications (1)

Publication Number Publication Date
CN116261200A true CN116261200A (en) 2023-06-13

Family

ID=86684457

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310326422.7A Pending CN116261200A (en) 2023-03-30 2023-03-30 An MPOLSR protocol optimization method based on link transmission quality for multi-UAV flight ad hoc network

Country Status (1)

Country Link
CN (1) CN116261200A (en)

Similar Documents

Publication Publication Date Title
Zheng et al. A mobility and load aware OLSR routing protocol for UAV mobile ad-hoc networks
CN113765544B (en) Method, system and device for multi-path concurrent transmission of dual-mode heterogeneous field area network based on HPLC and RF
Gharavi et al. Traffic scheduling technique for smart grid advanced metering applications
CN109639588B (en) Network congestion control routing method for aviation cluster
CN111132236B (en) MPR node selection method for multi-UAV self-organizing network based on improved OLSR protocol
Alzubi et al. Intelligent and dynamic neighbourhood entry lifetime for position-based routing protocol using fuzzy logic controller
Liu et al. Advanced fast recovery OLSR protocol for UAV swarms in the presence of topological change
CN114641050A (en) Energy-saving routing protocol method for unmanned aerial vehicle ad hoc network load sensing
Mule et al. Proactive source routing protocol for opportunistic data forwarding in MANETs
CN116456419B (en) Multi-unmanned aerial vehicle data return route method for energy balance high-reliability transmission
Quy et al. An Adaptive On-demand Routing Protocol With QoS Support for urban-MANETs.
CN116261200A (en) An MPOLSR protocol optimization method based on link transmission quality for multi-UAV flight ad hoc network
Verma et al. Stable energy proficient and load balancing based QoS routing in mobile Ad-Hoc networks: Mobile software based approach
Ahmad et al. Enhanced AODV route Discovery and Route Establishment for QOS provision for real time transmission in MANET
Huang et al. A routing algorithm based on cross-layer power control in wireless ad hoc networks
CN116782330A (en) Method and medium for selecting path in SD-UANET by considering load and link change
Abbas et al. Path diminution in node-disjoint multipath routing for mobile ad hoc networks is unavoidable with single route discovery
Bunu et al. Optimizing the Parameters of Relay Selection Model in OLSRv2 D2D Network
Vyas et al. Comparative analysis of routing protocols in MANETS
Selvakumar et al. Sustain route by tabu and amplifing qos with distributed scheduling in WMN
Abualola et al. Qos-olsr 2.0: A quality-of-service optimized link state routing protocol for mesh networks
Wang et al. A multi-point two-hop neighbor relay based OLSR routing protocol design method for self-organized multi-agent communication
CN115865775B (en) Unmanned aerial vehicle network rapid route recovery method based on OLSR
Hu et al. Performance analysis of routing protocols in uav ad hoc networks
Fu et al. Research of QoS routing algorithm in Ad Hoc networks based on reinforcement learning

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination