CN111787500B - 一种基于能量优先的移动充电车辆多目标充电调度方法 - Google Patents

一种基于能量优先的移动充电车辆多目标充电调度方法 Download PDF

Info

Publication number
CN111787500B
CN111787500B CN202010390838.1A CN202010390838A CN111787500B CN 111787500 B CN111787500 B CN 111787500B CN 202010390838 A CN202010390838 A CN 202010390838A CN 111787500 B CN111787500 B CN 111787500B
Authority
CN
China
Prior art keywords
node
energy
cluster
charging
mcv
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.)
Active
Application number
CN202010390838.1A
Other languages
English (en)
Other versions
CN111787500A (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.)
Zhejiang University of Technology ZJUT
Original Assignee
Zhejiang University of Technology ZJUT
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 Zhejiang University of Technology ZJUT filed Critical Zhejiang University of Technology ZJUT
Priority to CN202010390838.1A priority Critical patent/CN111787500B/zh
Publication of CN111787500A publication Critical patent/CN111787500A/zh
Application granted granted Critical
Publication of CN111787500B publication Critical patent/CN111787500B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W4/00Services specially adapted for wireless communication networks; Facilities therefor
    • H04W4/30Services specially adapted for particular environments, situations or purposes
    • H04W4/38Services specially adapted for particular environments, situations or purposes for collecting sensor information
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/23Clustering techniques
    • G06F18/232Non-hierarchical techniques
    • G06F18/2321Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
    • G06F18/23213Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions with fixed number of clusters, e.g. K-means clustering
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/004Artificial life, i.e. computing arrangements simulating life
    • G06N3/006Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04QSELECTING
    • H04Q9/00Arrangements in telecontrol or telemetry systems for selectively calling a substation from a main station, in which substation desired apparatus is selected for applying a control signal thereto or for obtaining measured values therefrom
    • 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/04Communication route or path selection, e.g. power-based or shortest path routing based on wireless node resources
    • H04W40/10Communication route or path selection, e.g. power-based or shortest path routing based on wireless node resources based on available power or energy
    • 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/20Communication route or path selection, e.g. power-based or shortest path routing based on geographic position or location
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W40/00Communication routing or communication path finding
    • H04W40/24Connectivity information management, e.g. connectivity discovery or connectivity update
    • H04W40/32Connectivity information management, e.g. connectivity discovery or connectivity update for defining a routing cluster membership
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W84/00Network topologies
    • H04W84/18Self-organising networks, e.g. ad-hoc networks or sensor networks
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04QSELECTING
    • H04Q2209/00Arrangements in telecontrol or telemetry systems
    • H04Q2209/40Arrangements in telecontrol or telemetry systems using a wireless architecture
    • H04Q2209/43Arrangements in telecontrol or telemetry systems using a wireless architecture using wireless personal area networks [WPAN], e.g. 802.15, 802.15.1, 802.15.4, Bluetooth or ZigBee
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04QSELECTING
    • H04Q2209/00Arrangements in telecontrol or telemetry systems
    • H04Q2209/80Arrangements in the sub-station, i.e. sensing device
    • H04Q2209/88Providing power supply at the sub-station
    • 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

Abstract

一种基于能量优先的移动充电车辆多目标充电调度方法,首先,采用K‑Means‑PER算法对网络进行聚类分簇,使得能量需求紧急度高的节点优先成为聚类中心;其次,综合考虑节点的能量和距离双重权重,采用RED‑CHS机制选择通信簇头,使其在普通节点与基站之间中转信息;然后,根据虚拟力场理论,采用VFF‑AS算法确定簇内锚点的最佳位置,使其能够兼顾簇内各节点的充电需求;最后,对蚁群系统进行改进,采用PER‑ACS算法对MCV进行路径规划和调度,使得移动路线最优,并根据请求状况自适应调节MCV对节点的充电量,精确求解MCV在各个锚点的停靠时长。本发明的PER‑MTCS方法给出了MCV的最优调度,有效解决了WRSN节点的能量补给问题。

Description

一种基于能量优先的移动充电车辆多目标充电调度方法
技术领域
本发明总体上归属于无线可充电传感器网络领域,具体涉及到网络节点聚类、无线能量传输、蚁群系统智能计算、路径规划以及移动调度等技术领域。
背景技术
传统的无线传感器网络(Wireless Sensor Network,WSN)是由大量低成本的微型传感器节点组成的多跳自组织网络,在军事侦查、健康医疗、环境监测、仓储管理和智能家居等领域有广阔的应用前景。但是,由于WSN节点体积小,其上搭载的电池容量有限,节点常常因为能量耗尽而死亡。依靠更换电池来延长网络寿命的方法往往操作不便、维护成本高、甚至有时根本不可行,这严重阻碍了WSN的进一步发展和应用。近年来,随着无线能量传输(Wireless Power Transfer,WPT)技术的快速发展,将能量补给装置放置在机器人或者移动设备上,以移动充电车辆(Mobile Charging Vehicle,MCV)的形式为节点提供稳定可持续的能量,已成为突破传统WSN节点能量瓶颈的重要方式。这种方式下,WSN节点能够按需、适时得到能量补充,从而形成了无线可充电传感器网络(Wireless Rechargeable SensorNetwork,WRSN)。与WSN相比,WRSN能够有效延长网络的生存周期,其应用的广度和深度得到了大幅度拓展与提升。
根据MCV对WRSN节点的充电方式不同,可将其分为一对一接触式充电和一对多非接触式充电两种形式。前者需要MCV移动至待充电节点处,以一对一接触的形式逐一为各个节点补充能量。然而,随着网络规模的扩大以及节点密度的增加,这种方式充电效率低,仅适合于节点较为稀疏的小规模网络。后者只需将MCV移动到合适的位置,就能够对其充电范围内的多个节点同时进行充电,这种批量并行操作的方式大大提高了网络的充电效率,因此可用于大规模无线传感器网络的能量补给。
在一对多非接触充电方式中,为了提高充电效率,可先对网络进行分簇,然后在每个簇内选择合适的停靠位置,即锚点(Anchor),使得MCV停留在锚点处对其簇内的多个节点同时进行无线充电。此外,为了提高通信效率,簇内的节点往往将数据发送给簇头,由簇头节点将数据转发给汇聚节点(Sink)。在这个过程中,有两个关键问题亟待解决:(1)如何对WRSN进行合理分簇,并确定通信簇头以及锚点的位置;(2)如何确定MCV对各个簇的充电次序以及充电停留时间,并调度MCV进行移动,使得其路线最优、充电效率最高。这两个因素均对WRSN的充电效果产生至关重要的影响。
发明内容
针对一对多非接触方式下对WRSN节点进行能量补给的问题,本发明提出了一种基于能量优先的MCV多目标充电调度方法——PER-MTCS(Priority of Energy Requirementbased Multiple Target Charging Scheduling)。该方法中,WRSN包括若干个带有可充电电池的同类传感器节点、一个携带有大容量可充电电池和无线充电装置的移动充电车辆MCV、一个用于给MCV补给能量的充电服务站、一个采集整个网络信息并调度MCV进行移动充电的基站;基站能够收集各个节点的数据信息、状态信息以及充电请求信息,并综合利用这些信息进行充电调度,最后将调度结果发送至MCV,由MCV执行具体的无线充电任务。整个系统的结构如图1所示。
本发明的技术方案为:
一种基于能量优先的MCV多目标充电调度方法,首先,采用基于能量优先的K-Means-PER(K-Means-Priority of Energy Requirement)方法对WRSN节点进行聚类,使得剩余能量越少、能量需求越紧急的节点优先成为聚类的中心,从而实现对WRSN网络的分簇;其次,采用RED-CHS(Residual Energy and Distance based Cluster Head Selection)算法来选择通信簇头,簇内其余节点通过簇头与基站进行信息交互;然后,采用VFF-AS(Virtual Force Field-Anchor Selection)算法来选择簇内锚点,使得在能够兼顾簇内各个节点能量需求的前提下,确定簇内最佳的锚点位置;最后,采用PER-ACS(Priority ofEnergy Requirement based Ant Colony System)算法对MCV进行最优路径规划与充电调度;
所述的K-Means-PER方法用于将整个WRSN网络节点聚类为不同的簇,以便提高其通信效率和MCV充电效率;
所述的RED-CHS算法用于选择簇内的通信簇头,使得其能够在簇内普通节点与网络基站之间中转信息;
所述的VFF-AS方法用于确定MCV在每个簇内的最佳停靠位置,使得MCV能够停留在该处为簇内的节点进行一对多非接触式充电;
所述的PER-ACS路径规划与充电调度方法用于求解MCV的最优移动路线,并确定在最优路径上每个站点的停靠时长,在该时长内MCV为相应簇内的节点进行无线充电。
进一步,所述的基于能量优先的K-Means-PER算法与传统的K-Means聚类算法不同,除了考虑节点之间的位置关系外,K-Means-PER算法更关注节点对能量的需求。由于能量需求与剩余能量有着密切联系,当前剩余能量越少,越需要补给能量,对充电需求越迫切,因此,K-Means-PER算法首先获取各个节点的剩余能量,依据其剩余能量多少对节点进行快速排序,剩余能量越少的节点在队列中的次序越靠前,充电优先级越高;然后取排序中的前K个节点作为初始的聚类中心;之后依据节点间的欧式距离进行聚类,最终获得K个以聚类中心为中心、半径不超过MCV有限充电范围的簇。
所述的通信簇头主要用于簇内节点与基站之间的信息中转,与普通通信方式相比,分簇机制通信效率更高,整个WRSN网络的通信开销更小。本发明中的RED-CHS通信簇头选择方法特征在于,它同时考虑了节点的能量状况以及节点与簇中心(即聚类中心)的位置关系双重因素,为节点设计了一个综合权重,作为它是否被选为簇头的度量依据。剩余能量越多、离簇中心距离越近的节点,权重越大,它当选为簇头的概率也越大。这是因为,一方面,在一个簇内,簇头由于要频繁地中转信息,故能量消耗很快,因此需要能量比较充足的节点担任簇头;另一方面,簇头负责中转簇内各个节点与基站之间的通信,从物理位置上讲,越是靠近簇中心位置的节点,越是能够兼顾簇内各个节点的通信需求。综上,那些剩余能量较多、同时又比较靠近簇中心的节点最合适作为簇头。
再进一步,在选择簇内锚点时,采用了VFF-AS方法,它借鉴了力学中的虚拟力场(Virtual Force Field)理论,将每个节点建模为质点,将其质量表示为该节点对能量需求的紧急程度;假设存在一个虚拟质点,根据牛顿万有引力定律,可以得到簇内每个节点对该虚拟质点的引力,进而求出虚拟质点所受的合力;合力为0的点即为平衡点,该点处能够兼顾到各个节点的引力需求,因此,该位置即为锚点的最佳位置。基于VFF-AS方法的锚点位置选择方法示意图如图3所示。
在规划MCV最优路径并调度其进行多目标充电时,采用了基于能量优先的蚁群系统方法PER-ACS,其过程为:(1)采用蚁群系统来解决MCV的路径规划问题,并对普通的蚁群系统进行改进,将锚点对能量需求的紧急程度整合到到蚁群系统的启发信息中,同时采用伪随机比例规则优化了蚂蚁从一个锚点到下一个锚点的转移决策,使得其能在先验知识和探索新路径之间做出权衡和调节,扩大了路径的搜索范围,避免了算法陷入局部最优,从而找到一条MCV移动充电的全局最优路径;(2)根据一轮调度中的充电请求数量动态调节MCV对簇内节点的充电量,并根据当前簇内节点的剩余能量信息和MCV的充电功率,精确地计算出MCV在各个锚点的停靠时长,在该时长内MCV对簇内的各个节点以一对多非接触方式进行无线充电。
本发明的有益效果为:充分考虑网络节点的能量需求,并结合聚类算法、蚁群系统和路径规划等技术,给出了MCV的最优调度,从而有效解决了WRSN节点的能量补给问题。
附图说明
图1是本发明的系统结构示意图;
图2是基于能量优先的MCV多目标充电调度方法PER-MTCS的实现流程图;
图3是基于虚拟力场的锚点位置选择方法VFF-AS的示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
现在将详细参考本发明的实施例,这些实施例的示例在附图中示出。文中“停靠点”和“锚点”在此用于方便描述,并且因此可以可交换地被使用,而且没有任何可区别的意义或功能。同样,文中“聚类中心”和“簇中心”也可以可交换地被使用,而且没有任何可区别的意义或功能。
图1给出了本发明的系统结构示意图,其核心组件包括可充电的传感器节点、网络基站、无线移动充电车辆和充电服务站四部分,各部分的功能如下:
(1)可充电的传感器节点(Rechargeable Sensors,RSs,简称节点)
整个WRSN网络由大量安装有可充电电池的传感器节点组成,每个节点感知并采集周围环境的信息,多个传感器以自组织的形式相互协作,共同完成对环境的监测任务。本发明中每个传感器节点同质同构,它们具有相同的计算能力和通信能力;每个节点所配备的可充电电池容量相等,最大容量均为Emax;MCV为节点充电时,若节点的能量达到最大容量Emax后,它不再接收能量;每个节点均能够通过GPS(Global Position System,全球定位系统)装置自行定位,且一旦部署,其位置就静止不动。
(2)网络基站(Base Station,BS,简称基站)
基站是整个WRSN网络的数据中心和控制中心,相当于网络的汇聚节点(Sink),有恒定的电源供电。它不仅能够以无线多跳形式收集各个节点采集的数据,进行分析和处理,还能获取各个节点的状态信息(例如:节点的位置坐标、节点当前的剩余能量等),并能接收节点的充电请求,然后调度MCV对节点进行一对多非接触式充电。
需要说明的是,每个普通的WRSN节点并不是直接与基站进行通信,而是首先对WRSN网络进行聚类,一个聚类形成一个簇,然后从中选出簇头,由簇头负责本簇内节点与基站之间的信息交互。
簇内节点向簇头发送的信息包括所采集的普通数据信息和充电请求信息等,簇头能对不同类型的信息进行区分和处理。具体地,在一段时间内,簇头可能会收到多个节点发来的充电请求,这时,它能够从中辨别出紧急度最高的请求,并将该请求转发到基站,由基站进行最终的充电调度。
(3)无线移动充电车辆(Mobile Charging Vehicle,MCV)
MCV是WRSN中的能量补给设备,它具有以下功能:
①携带有充足的电量(远远大于普通节点所需的充电量),具有无线充电能力,并且自身能够控制充电的开启与停止;
②具有无线通信能力,能够以无线方式接收基站发来的调度信息;
③具有移动能力,能够按照指定的调度路线在锚点间匀速移动;
④能够以一对多非接触方式为其充电半径内的所有节点进行无线充电;
初始运行时,系统根据各个节点的剩余能量进行聚类分簇,再根据虚拟力场理论确定簇内的最佳停靠点,即锚点。在WRSN运行过程中,若某节点有充电需求,则向自己所在簇的簇头发充电请求,簇头将此请求转发至基站。在一段时间内,基站根据所收到的充电请求状况,采用基于能量优先的蚁群系统进行调度,并将调度结果发送给MCV。MCV根据调度结果,沿着基站所规划的最优路径移动,依次停靠到各个锚点位置,为有充电需求的簇进行一对多无接触式充电。
(4)充电服务站(Charging Service Station,CSS)
充电服务站相当于一个后勤服务站,能够为MCV快速提供能量补给。随着MCV为其它节点进行能量输出,其自身的能量会逐渐减少,因此需要适时进行补给,以满足整个WRSN网络的充电需求。本发明中,当所有节点都没有充电需求时,MCV停留在充电服务站,为自身充满能量,然后处于待命状态;当接收到基站的调度信息后,它沿着最优路径依次为有能量需求的簇进行充电,最后再回到充电服务站补满能量,等待下一轮调度。
本发明公开了一种WRSN中基于能量优先的MCV多目标充电调度方法——PER-MTCS。图2给出了本发明实施例的流程图。PER-MTCS调度方法包括以下四步:(1)基于能量优先的K-Means-PER节点聚类分簇;(2)基于RED-CHS机制的通信簇头选择;(3)基于VFF-AS的最佳锚点确定方法;(4)基于PER-ACS的MCV路径规划与调度。各个步骤的具体实施方式如下。
步骤一:基于能量优先的K-Means-PER节点聚类分簇,过程如下:
步骤1.1:计算各节点的能量消耗速率
根据公式计算当前时刻t节点i的能量消耗速率/>其中,/>eij(t)和eiB(t)分别表示t时刻节点i接收1bit数据、向节点j发送1bit数据、向基站BS发送1bit数据所消耗的能量;dfki(t)、dfij(t)和dfiB(t)分别表示t时刻节点k到节点i、节点i到节点j、节点i到基站BS的数据流速率,单位均为bit/s;Ni-recv是节点i的邻居中向i发送数据的节点所形成的集合,即:节点i从这些邻居节点接收数据;Ni-send是节点i的邻居中从i发送数据所到节点形成的集合,即:节点i发送数据到这些邻居节点;/>为节点i接收数据时的总能量消耗率,/>为节点i发送数据时的总能量消耗率。由于WRSN节点收、发数据所消耗的能量远大于节点感知数据和自身运行所消耗的能量,故这里仅计算节点i收发数据时的能量消耗速率,而节点感知数据和自身运行所消耗的能量忽略不计。
步骤1.2:计算各节点的能量接收功率
根据公式计算节点i在t时刻的能量接收功率/>本发明中,MCV为能量发送器,即无线充电设备,节点为能量接收器,两者之间以无线形式传输能量。节点i的能量接收功率也即MCV对该节点的充电功率。MCV的充电半径为R,它与节点之间的欧式距离为d。当d≤R时,即节点在MCV的有效充电范围内,按照Friis无线信道传输模型/>为节点i进行充电;否则,当d>R时,节点超出了MCV的有效充电范围,MCV无法为其充电;Pemission为发送功率,Preceive为接收功率,Gs为发射天线增益,Gr为接收天线增益,ηrect为整流器的效率,λ为波长,Lp为极化损耗,d为发送器与接收器之间的距离,ξ为常数,用于调整短距离传输时的Friis自由空间传输方程。
步骤1.3:计算各节点的剩余能量
根据公式计算节点i在t时刻的剩余能量/>其中,Emax为节点可充电电池的最大容量,为节点i在t时刻所消耗的能量,根据(式1)求出;/>为节点i在t时刻所接收到的能量,根据(式2)求出。
步骤1.4:依据剩余能量值对节点进行排序
根据(式3)计算的剩余能量多少,对WRSN网络节点进行排序,剩余能量越少的节点次序越靠前。然后,选取前K个节点作为初始聚类中心。K为常数,在系统初始时设定。
设K个初始聚类中心所在的位置依次为当前轮聚类分簇时,K个簇中心位置依次为/>则分簇开始时,
步骤1.5:将每个节点划分到簇中
对于节点i,计算它与当前K个簇中心的距离,从中找出距离最近且小于R(R是MCV的有效充电半径)的那个簇中心,将其归为该簇。依次遍历所有节点,最终形成K个分簇。
步骤1.6:计算本轮分簇时各个簇的簇中心
对于第m个分簇Cm,按照公式计算其新的簇中心位置/>其中,|Cm|为Cm簇中节点的个数,(x,y)为Cm簇内节点的坐标。/>和/>分别表示Cm簇中所有节点横坐标与纵坐标的均值。对所有的簇进行以上操作,得到K个新的簇中心/>
步骤1.7:更新簇中心
对于每个簇,判断新的簇中心位置与上一轮的簇中心位置是否相同。若不同,用新的簇中心位置去替换旧的簇中心位置,即:/>否则,簇中心保持不变。
然后转到步骤1.5继续新一轮的分簇,直到所有簇的簇中心均不再变化为止。此时,每个簇的簇中心及其簇内的节点确定了下来。
步骤1.8:输出分簇结果
输出分簇结果C={C1,C2,…,Cm,…,CK},其中,Cm是第m个簇内的节点所形成的集合。
步骤二:基于RED-CHS机制的通信簇头选择,过程如下:
步骤2.1:计算簇内每个节点的选举权重
对于簇Cm内的节点i,根据公式计算其选举权重ωi。其中,/>为节点i的剩余能量,Emax为节点的满能量值,|Cm|为Cm簇中节点的个数,di-center和dj-center分别为节点i和节点j到簇中心的距离,/>是能量与距离选举因素权重系数,为常量,但值可根据需要进行调整;ε为节点选举权重调节因子,用来平衡节点的剩余能量比与距离比两者之间的差距,使其处于同一数量级水平。
步骤2.2:根据节点的选举权重确定通信簇头
按照选举权重大小为簇内每个节点进行降序排序,选择权重最大的节点作为该簇的通信簇头。这种机制使得那些剩余能量越多、越靠近簇中心的节点被选为通信簇头的概率越大。
步骤三:基于VFF-AS的最佳锚点确定方法,过程如下:
步骤3.1:计算簇内每个节点的质量
根据公式求得Cm簇内节点i的质量Mi,这里的质量Mi表示节点i对能量需求的紧急程度。其中,/>为节点i的剩余能量;为其能量消耗速率;ψ1和ψ2是剩余能量和能量消耗速率对能量需求紧急程度的影响权重,为常数,但值可根据需要进行调整;δ为节点质量调节因子,用来平衡剩余能量与能量消耗速率两者之间的差距,使其处于同一个数量级水平。剩余能量越少,能量消耗速率越快,则节点对能量需求的紧急程度越高,其所表示的质量就越大。
步骤3.2:计算簇内各节点对虚拟质点的引力
设簇Cm内存在一个虚拟质点Svirtual,其质量为Mvirtual,根据牛顿万有引力定律,按照公式求出该簇内节点i对虚拟质点Svirtual的引力/>其中,G是引力常数,Mi和Mvirtual分别表示节点i和虚拟质点Svirtual的质量,di-virtual是两者之间的距离,/>是节点i到虚拟质点Svirtual的单位向量。
步骤3.3:计算虚拟质点所受的合力,确定锚点的最佳位置
根据公式计算虚拟质点Svirtual所受的合力Fresultant。其中,|Cm|为Cm簇中节点的个数,/>为簇内节点i对虚拟质点Svirtual的引力。合力为0的点即为平衡点,该点处能够兼顾簇内各个节点对能量的需求,故该位置即为锚点的最佳位置。
基于虚拟力场的最佳锚点位置选择方法VFF-AS的示意图如图3所示。
步骤四:基于PER-ACS的MCV路径规划与充电调度
本发明中MCV路径规划与充电调度使用了基于能量优先的改进的蚁群系统来实现,其具体实施步骤如下:
步骤4.1:基站BS收集WRSN网络信息,并设置调度基本参数。
WRSN网络的每个节点定期更新自身信息,包括能量消耗速率、剩余能量、分簇信息和节点质量(表征该节点对能量需求的紧急程度)等。当节点的剩余能量低于设定的阈值时,它向通信簇头发送充电请求。簇头在一段时间内收到各个节点的充电请求后,从中找出质量最大的节点,将其充电请求转发给基站BS。BS定期收集各个通信簇头发来的充电请求,调度MCV执行充电任务。
路径规划与充电调度时,需设置以下基本参数:
K:WRSN网络被划分成的总簇数;
κ:本轮调度时基站收到簇头发来的充电请求个数;
Emax:节点的初始能量;
Ereq-th:节点的充电请求阈值;
Nant:蚁群系统的蚂蚁总数量;
Niter:蚁群系统的总迭代次数;
Ω:每个路段上的信息素浓度初始值;
ρ:全局信息素挥发因子;
μ:局部信息素挥发因子;
α:信息素启发式因子;
β:距离期望启发式因子;
γ:能量需求紧急程度期望启发式因子。
以上参数中,除了κ会随着节点的充电请求状况而动态变化外,其它参数均是常数,可根据需要设定。但是,一旦基站BS启动了路径规划过程,在同一轮调度中κ值保持不变。
除以上参数外,还需设置蚁群系统迭代计数变量Nc,其初值为1。
步骤4.2:初始部署蚂蚁
随机将Nant只蚂蚁部署在κ个有充电请求的锚点上。
步骤4.3:当迭代计数Nc≤Niter时,执行步骤4.4~4.7
步骤4.4:对于每一只蚂蚁ANTl(l=1,2,…,Nant),执行步骤4.5~4.7
步骤4.5:对于蚂蚁ANTl,选择当前的下一跳锚点
由于本轮调度时基站收到κ个簇头发来的充电需求,故MCV需要停靠κ个锚点为相应簇内的节点充电,这κ个锚点的编号为1,2,…,κ,对应的锚点集为{A1,A2,…,Aκ}。对于任一只蚂蚁ANTl,设其当前位于锚点n处,根据(式9)的伪随机比例规则确定它到下一个锚点l的转移决策。
(式9)中,τln(t)为t时刻锚点l到锚点n路段上的信息素浓度,初值为常数Ω,即:τln(0)=Ω,α为信息素启发式因子,表征信息素的重要程度;其中,dln(t)为t时刻锚点l到锚点n之间路段的长度,β为距离期望启发式因子,表征距离信息的重要程度;Mn-max(t)为t时刻锚点n所在簇中节点的最大质量,也即能量需求最迫切的节点的充电紧急程度,γ为能量需求紧急程度期望启发式因子,表征能量需求紧急程度的重要性;/>为t时刻蚂蚁ANT从锚点l到锚点n的转移概率;/>为蚂蚁ANTl下一步允许选择的锚点集合,在一轮路径搜索过程中,蚂蚁走过的站点被列入禁忌表,只有那些不在禁忌表中的尚未访问过的站点才可能作为下一站的后备访问点,即/>中的点,这样可以确保一轮路径搜索时蚂蚁不会重复访问任一路段。
(式9)采用了伪随机比例规则。q是在[0,1]区间均匀分布的随机变量,参数q0(0≤q0≤1)是当前选择最优移动方式的概率。当q≤q0时,蚂蚁根据先验知识选择目前最优路径;而当q>q0时,则按照概率探索新的路径。参数q0用来调节利用先验知识与探索新路径之间的相对重要性。
步骤4.6:局部信息素更新
对于锚点l到锚点n之间的路段,按照公式τln(t+T)=(1-μ)τln(t)+μτ0……(式10)进行局部信息素更新。其中,τln(t)为当前t时刻锚点l到锚点n路段上的信息素浓度,τln(t+T)为经过时间T后该路段上的信息素浓度;μ为局部信息素挥发因子,取μ=0.1时可以取得较好的搜索结果;其中κ为本轮调度的充电锚点个数,Lnn是根据最近邻原则计算而来的路径长度,即:对于κ个要充电的锚点,MCV随机从任意一个锚点出发,每次均选择离自己最近的锚点作为下一个要充电的目标,遍历所有κ个锚点得到的路径长度即为Lnn。局部更新规则会使蚂蚁访问过的路段上的信息素浓度逐渐减少,进而降低其它蚂蚁选择该路段的概率,提升蚂蚁探索未知路径的机率,从而避免了过早收敛于同一搜索路径。
步骤4.7:循环边界判定
判断蚂蚁ANTl是否访问完所有κ个锚点,形成一条完整的路径:
(1)若否,则转到步骤4.5处,继续寻找下一跳锚点,直到将所有κ个锚点访问完毕;
(2)若是,则蚂蚁ANTl已形成一条自己的访问路径,此时,转到步骤4.4处,继续下一只蚂蚁的路径搜索(l=l+1),直到所有Nant只蚂蚁都形成自己完整的访问路径。
步骤4.8:全局信息素更新
当前轮所有蚂蚁都形成自己的路径后,从中找出最优的一条路径对该路径上的信息素进行加强和更新。具体地,根据公式τln(t+T)=(1-ρ)τln(t)+ρΔτln(t+T)?…(式11)进行全局信息素更新。其中,τij(t)为当前t时刻锚点l到锚点n路段上的信息素浓度,τln(t+T)为经过时间T后该路段上的信息素浓度;ρ为全局信息素挥发因子;Δτln(t+T)为经过时间T后锚点l到锚点n路段上的信息素浓度增量,可根据公式求出。(式12)中,Nant为蚂蚁的总数;/>为第n只蚂蚁在ln路段上产生的信息素浓度增量,根据公式/>求出。(式13)中,/>为第Nc轮路径搜索时Nant只蚂蚁找到的当前最优路径,Q为常量。
步骤4.9:更新迭代次数Nc,判断是否路径搜索结束
若当前迭代计数Nc达到了总迭代次数Niter,则路径搜索结束,即为最终找到的最优路径,即/>否则,迭代计数更新,Nc=Nc+1,然后转到步骤4.3处继续执行,直到迭代次数达到Niter为止。
步骤4.10:输出最优的访问路径
输出最终的最优访问路径vj∈{1,2,...,κ}。其中,是本轮调度时MCV停靠的第j个锚点,该锚点在锚点集{A1,A2,…,Aκ}中的序号为vj
步骤4.11:计算MCV对节点的充电停止能量阈值
根据公式计算MCV对节点的充电停止能量阈值Estop-th。其中,Emax为节点的最大能量;Ereq-th为节点发送充电请求的能量阈值,一旦节点低于该值,就请求补充能量;K为所有锚点的总数,κ为本轮调度中有充电请求的锚点数量。由(式14)可得,当同一时刻请求充电的节点较多时,为了兼顾多节点的充电请求,MCV对各个节点的充电量较少;反之,当请求充电的节点较少时,MCV可为节点提供比较充裕的能量供给。这种能根据节点请求状况自适应调节充电量的方式能够更加灵活地满足整个WRSN网络的能量补给。
步骤4.12:确定MCV在最优路径上各个锚点的停靠时间
根据公式计算MCV停靠在锚点/>处为相应簇内节点充电的时间。其中,Estop-th为本轮调度时的充电停止能量阈值,/>为锚点/>所在的簇,/>为簇/>内的节点在时刻t剩余能量的最小值,/>为簇/>内的节点在时刻t接收功率的最小值。
按照4.1至4.12所描述的步骤,基站BS可以求出本轮有请求的κ个簇的最优访问路径LBest,然后按照该路径依次停靠在锚点在锚点/>处,MCV停靠的时长为/>该时长内MCV为锚点/>所在簇内的所有节点进行一对多无接触式充电。这即为基于PER-ACS的MCV路径规划与充电调度的全部过程。
综上,本发明中基于能量优先的MCV多目标充电调度方法PER-MTCS通过步骤一的K-Means-PER算法实现WRSN网络的聚类分簇,通过步骤二的RED-CHS机制来选择合适的通信簇头,通过步骤三的VFF-AS方法来选择簇内最佳的锚点位置,通过步骤四的PER-ACS算法来实现MCV的最优路线规划与充电调度。
以上所述的具体描述,对发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例,用于解释本发明,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (6)

1.一种基于能量优先的移动充电车辆多目标充电调度方法,其特征在于,首先,采用K-Means-PER方法对WRSN网络进行聚类分簇,使得能量需求紧急度高的节点优先成为聚类中心;其次,综合考虑节点能量和距离信息双重权重,采用RED-CHS机制选择通信簇头,使其在普通节点与基站之间中转信息;然后,根据虚拟力场理论,采用VFF-AS算法确定簇内锚点的最佳位置,使其能够兼顾簇内各节点的充电需求;最后,对蚁群系统进行改进,采用PER-ACS算法对MCV进行路径规划和调度,使得移动路线最优,并根据请求状况自适应调节MCV对节点的充电量,精确求解MCV在各个锚点的停靠时长;
所述的K-Means-PER方法用于将整个WRSN网络节点聚类为K个不同的簇,以提高其通信效率和MCV充电效率;
所述的RED-CHS机制用于选择簇内的通信簇头,使得其能够在簇内普通节点与网络基站之间中转信息;
所述的VFF-AS算法用于确定MCV在每个簇内的最佳停靠位置,使得MCV能够停留在该处为簇内的节点进行一对多非接触式充电;
所述的PER-ACS路径规划与充电调度方法用于求解MCV的最优移动路线,并确定在最优路径上每个站点的停靠时长,在该时长内MCV为相应簇内的节点进行无线充电;
K-Means-PER聚类分簇方法中,除了考虑节点之间的位置关系与距离信息外,更关注节点对能量的需求;它根据剩余能量多少对节点进行快速排序,使得剩余能量越少的节点充电优先级越高;然后取排序中前K个节点作为初始聚类中心,之后依据节点间的欧氏距离进行聚类,最终获得K个以聚类中心为中心、半径不超过MCV有限充电范围的簇;
所述的RED-CHS通信簇头选择机制中,同时考虑了节点的剩余能量状况以及节点与簇中心的欧氏距离双重因素,为节点设计了一个综合权重,作为它是否被选为簇头的度量依据;依据该综合权重,剩余能量越多、离簇中心距离越近的节点当选为簇头的概率越大,这样使得所选出的簇头不仅有充足的能量胜任簇内节点与基站之间的信息转发,而且位置上也比较靠近簇中心,从而能够同时兼顾簇内各个节点的通信需求;
所述的VFF-AS簇内锚点最佳位置确定方法中,采用了力学中的虚拟力场理论,将WRSN节点建模为质点,将节点质量表示为该节点对能量需求的紧急程度;假设在簇内存在一个虚拟质点,根据牛顿万有引力定律,可以求解簇内每个节点对该虚拟质点的引力,进而得出虚拟质点所受的合力;利用力学理论与方法求解虚拟质点的受力平衡点,由于该平衡点处能够兼顾到各个节点的引力需求,因此该位置即为锚点的最佳位置;
所述的PER-ACS路径规划与充电调度方法中,采用了改进的蚁群系统,将锚点对能量需求的紧急程度整合到到蚁群系统的启发信息中,同时采用伪随机比例规则优化了蚂蚁从一个锚点到下一个锚点的转移决策,使得其能在先验知识和探索新路径之间做出权衡和调节,扩大了路径的搜索范围,避免了算法陷入局部最优,从而找到一条MCV移动充电的最优路径;根据一轮调度中的充电请求数量自适应调节MCV对簇内节点的充电量,并根据当前簇内节点的剩余能量信息和MCV的充电功率,精确地求解MCV在各个锚点的停靠时长,在该时长内MCV对簇内的各个节点以一对多非接触方式进行无线充电。
2.如权利要求1所述的一种基于能量优先的移动充电车辆多目标充电调度方法,其特征在于,节点的剩余能量由其能量消耗速率和能量接收功率综合计算而来;对于任意节点i,其能量消耗速率、能量接收功率和剩余能量分别按照(式1)、(式2)和(式3)的方法求出:
其中,为当前时刻t节点i的能量消耗速率,/>eij(t)和eiB(t)分别表示t时刻节点i接收1bit数据、向节点j发送1bit数据、向基站BS发送1bit数据所消耗的能量;dfki(t)、dfij(t)和dfiB(t)分别表示t时刻节点k到节点i、节点i到节点j、节点i到基站BS的数据流速率,单位均为bit/s;Ni-recv是节点i的邻居中向i发送数据的节点所形成的集合,即:节点i从这些邻居节点接收数据;Ni-send是节点i的邻居中从i发送数据所到节点形成的集合,即:节点i发送数据到这些邻居节点;/>为节点i接收数据时的总能量消耗率,/>为节点i发送数据时的总能量消耗率;/>为节点i在t时刻的能量接收功率,Pemission为MCV的发送功率,Gs为发射天线增益,Gr为接收天线增益,ηrect为整流器的效率,λ为波长,Lp为极化损耗,d为发送器与接收器之间的距离,R为MCV的充电半径,ξ为常数,用于调整短距离传输时的Friis自由空间传输方程;/>为节点i在t时刻的剩余能量,Emax为节点可充电电池的最大容量。
3.如权利要求1所述的一种基于能量优先的移动充电车辆多目标充电调度方法,其特征在于,节点的综合权重采用的方法求出,其中,/>和/>分别是剩余能量比权重与距离比权重,ε可调节这双重因素的放大比例,使其处于同一数量级水平;/>为节点i的剩余能量,Emax为节点可充电电池的最大容量,|Cm|为Cm簇中节点的个数,Cm是第m个簇内的节点所形成的集合,m∈{1,2,...,K},di-center和dj-center分别为节点i和节点j到簇中心的距离。
4.如权利要求1所述的一种基于能量优先的移动充电车辆多目标充电调度方法,其特征在于,所述的VFF-AS算法中,节点的质量采用的方法求得,其中,ψ1和ψ2分别是剩余能量和能量消耗速率对能量需求紧急程度的影响权重,节点质量调节因子δ能够调整两者之间的放大比例,这种方法可使剩余能量越少、能量消耗率越快的节点对充电需求的紧急度越高;/>为节点i的剩余能量;/>为其能量消耗速率。
5.如权利要求1所述的一种基于能量优先的移动充电车辆多目标充电调度方法,其特征在于,改进的蚁群系统中,蚂蚁从一个锚点到下一个锚点的选择时采用(式6)的转移决策,并在蚂蚁寻径时采用局部信息素更新和全局信息素更新两种更新机制,避免了算法过早收敛于同一条路径;
其中,κ为本轮调度的充电锚点个数,τln(t)为t时刻锚点l到锚点n路段上的信息素浓度,初值为常数Ω,即:τln(0)=Ω,α为信息素启发式因子,表征信息素的重要程度;其中,dln(t)为t时刻锚点l到锚点n之间路段的长度,β为距离期望启发式因子,表征距离信息的重要程度;Mn-max(t)为t时刻锚点n所在簇中节点的最大质量,也即能量需求最迫切的节点的充电紧急程度,γ为能量需求紧急程度期望启发式因子,表征能量需求紧急程度的重要性;/>为t时刻蚂蚁ANT从锚点l到锚点n的转移概率;allowedANT为蚂蚁ANT下一步允许选择的锚点集合,在一轮路径搜索过程中,蚂蚁走过的站点被列入禁忌表,只有那些不在禁忌表中的尚未访问过的站点才可能作为下一站的后备访问点,即allowedANT中的点,这样可以确保一轮路径搜索时蚂蚁不会重复访问任一路段;q是在[0,1]区间均匀分布的随机变量,参数q0是当前选择最优移动方式的概率,0≤q0≤1;当q≤q0时,蚂蚁根据先验知识选择目前最优路径;而当q>q0时,则按照概率/>探索新的路径,参数q0用来调节利用先验知识与探索新路径之间的相对重要性。
6.如权利要求1所述的一种基于能量优先的移动充电车辆多目标充电调度方法,其特征在于,在动态调整MCV对簇内节点的充电量时,基站BS能够根据当前一轮调度中所收到的锚点充电请求数、总的锚点数、节点的充电请求阈值以及节点电池的最大容量,采用的方法求得MCV为本轮请求锚点的充电停止阈值,其中,Estop-th为MCV对节点的充电停止能量阈值,Emax为节点可充电电池的最大容量;Ereq-th为节点发送充电请求的能量阈值,一旦节点低于该值,就请求补充能量;K为所有锚点的总数,κ为本轮调度中有充电请求的锚点数量;然后利用计算结果,并结合锚点所在簇内各节点的剩余能量和充电功率,采用/>方法精确求解MCV在最优路径上各个锚点的停靠时间,其中,/>为MCV停靠在锚点/>处为相应簇内节点充电的时间,Estop-th为本轮调度时的充电停止能量阈值,/>为锚点/>所在的簇,vn∈{1,2,...,κ},/>为簇/>内的节点在时刻t剩余能量的最小值,为簇/>内的节点在时刻t接收功率的最小值,使得MCV的充电效率最大。
CN202010390838.1A 2020-05-11 2020-05-11 一种基于能量优先的移动充电车辆多目标充电调度方法 Active CN111787500B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010390838.1A CN111787500B (zh) 2020-05-11 2020-05-11 一种基于能量优先的移动充电车辆多目标充电调度方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010390838.1A CN111787500B (zh) 2020-05-11 2020-05-11 一种基于能量优先的移动充电车辆多目标充电调度方法

Publications (2)

Publication Number Publication Date
CN111787500A CN111787500A (zh) 2020-10-16
CN111787500B true CN111787500B (zh) 2023-07-25

Family

ID=72753496

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010390838.1A Active CN111787500B (zh) 2020-05-11 2020-05-11 一种基于能量优先的移动充电车辆多目标充电调度方法

Country Status (1)

Country Link
CN (1) CN111787500B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112261704B (zh) * 2020-10-23 2022-12-09 广州大学 一种基于可充电无线传感器网络的重复博弈路由方法
CN112907153B (zh) * 2021-01-15 2022-11-01 中原工学院 一种混合场景下考虑用户多种需求的电动汽车调度方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109495945A (zh) * 2018-12-18 2019-03-19 河海大学常州校区 一种WSNs中基于聚类分簇的协同充电方法
CN109511150A (zh) * 2018-06-20 2019-03-22 河海大学常州校区 一种基于单对多充电技术的移动充电车路径规划算法
CN110061538A (zh) * 2019-03-22 2019-07-26 中山大学 一种wsn节点智能分簇及移动充电设备路径规划方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9122693B2 (en) * 2010-11-30 2015-09-01 Nokia Technologies Oy Method and apparatus for determining contextually relevant geographical locations

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109511150A (zh) * 2018-06-20 2019-03-22 河海大学常州校区 一种基于单对多充电技术的移动充电车路径规划算法
CN109495945A (zh) * 2018-12-18 2019-03-19 河海大学常州校区 一种WSNs中基于聚类分簇的协同充电方法
CN110061538A (zh) * 2019-03-22 2019-07-26 中山大学 一种wsn节点智能分簇及移动充电设备路径规划方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
"An Energy-aware Anchor Point Selection Strategy for Wireless Rechargeable Sensor Networks";Fang Liu and Shuaishuai Lu;《2019 IEEE/CIC International Conference on Communications in China(ICCC)》;第422-426页 *
"无线可充电传感器网络中充电规划及其可调度性研究";胡诚;《中国博士学位论文全文数据库 计算机应用技术》;全文 *
"无线可充电传感器网络中最优充电策略的研究";梅焕;《中国优秀硕士学位论文全文数据库 电子科学与技术》;全文 *

Also Published As

Publication number Publication date
CN111787500A (zh) 2020-10-16

Similar Documents

Publication Publication Date Title
CN107835499B (zh) 一种WSNs中基于分簇和能量中继的移动充电方法
CN109495945B (zh) 一种WSNs中基于聚类分簇的协同充电方法
Han et al. A multicharger cooperative energy provision algorithm based on density clustering in the industrial Internet of Things
Xu et al. Minimizing the longest charge delay of multiple mobile chargers for wireless rechargeable sensor networks by charging multiple sensors simultaneously
CN111787500B (zh) 一种基于能量优先的移动充电车辆多目标充电调度方法
CN110418434A (zh) 一种无人机辅助的无线传感器网络充电方法与装置
CN107623901B (zh) 一种WRSNs中联合数据收集及能量补给方法
Xu et al. On-demand energy replenishment for sensor networks via wireless energy transfer
CN109451556A (zh) 基于uav对无线传感网充电的方法
Yu et al. Charging strategy and scheduling algorithm for directional wireless power transfer in WRSNs
Chelbi et al. A new hybrid routing protocol for wireless sensor networks
Liu et al. Maximizing sensor lifetime via multi-node partial-charging on sensors
Lei et al. An energy urgency priority based mobile charging scheme in Wireless Rechargeable Sensor Network
Zhao et al. Hybrid scheduling strategy of multiple mobile charging vehicles in wireless rechargeable sensor networks
CN110662175B (zh) 一种基于无线可充电传感器网络的移动车速度控制方法
CN110049500B (zh) 基于模拟退火算法的无线可充电传感网中uav能量补偿方法
CN112702688A (zh) 结合能量补充和数据收集的移动小车规划方法
Singh et al. An efficient approach for wireless rechargeable sensor networks for vehicle charging
Rao et al. Joint wireless charging and data collection using mobile element for rechargeable wsns
Prakash et al. Mixed Linear Programming for Charging Vehicle Scheduling in Large-Scale Rechargeable WSNs
CN110248330A (zh) 一种基于中继充电模型最大化充电小车休息时间调度方法
Yi et al. Multi-Task Transfer Deep Reinforcement Learning for Timely Data Collection in Rechargeable-UAV-aided IoT Networks
Zhao et al. Directional charging-based scheduling strategy for multiple mobile chargers in wireless rechargeable sensor networks
Fan et al. Exploiting an unmanned aerial vehicle for energy efficient data collection in wireless sensor networks
Rahaman et al. An efficient charging scheduling scheme to enhance the wireless rechargeable sensor networks' lifespan

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
GR01 Patent grant
GR01 Patent grant