CN106094833B - 一种基于混沌萤火虫算法的移动机器人路径规划方法 - Google Patents

一种基于混沌萤火虫算法的移动机器人路径规划方法 Download PDF

Info

Publication number
CN106094833B
CN106094833B CN201610569272.2A CN201610569272A CN106094833B CN 106094833 B CN106094833 B CN 106094833B CN 201610569272 A CN201610569272 A CN 201610569272A CN 106094833 B CN106094833 B CN 106094833B
Authority
CN
China
Prior art keywords
path
firefly
chaos
mobile robot
swarm algorithm
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
CN201610569272.2A
Other languages
English (en)
Other versions
CN106094833A (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.)
Chongqing University of Post and Telecommunications
Original Assignee
Chongqing University of Post and Telecommunications
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 Chongqing University of Post and Telecommunications filed Critical Chongqing University of Post and Telecommunications
Priority to CN201610569272.2A priority Critical patent/CN106094833B/zh
Publication of CN106094833A publication Critical patent/CN106094833A/zh
Application granted granted Critical
Publication of CN106094833B publication Critical patent/CN106094833B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/02Control of position or course in two dimensions
    • G05D1/021Control of position or course in two dimensions specially adapted to land vehicles
    • G05D1/0212Control of position or course in two dimensions specially adapted to land vehicles with means for defining a desired trajectory
    • G05D1/0217Control of position or course in two dimensions specially adapted to land vehicles with means for defining a desired trajectory in accordance with energy consumption, time reduction or distance reduction criteria

Landscapes

  • Engineering & Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Feedback Control In General (AREA)

Abstract

本发明请求保护一种基于混沌萤火虫算法的移动机器人路径规划方法,该方法包括步骤:S1、确定移动机器人系统代价函数;S2、对移动机器人系统进行初始化,得到所有可能路径;S3、利用混沌萤火虫算法对每一条路径进行迭代更新,包括对亮度和位置的更新;S4、对更新后的路径进行排序,得到局部最优路径;S5、迭代次数完成进行步骤S6,否则进行步骤S3;S6、利用最优调整策略调整当前局部最优路径;S7、输出全局最优路径。本发明能够使移动机器人在最短时间内寻找到通过狭窄区域的最优路径。

Description

一种基于混沌萤火虫算法的移动机器人路径规划方法
技术领域
本发明涉及移动机器人导航领域,特别是一种基于混沌萤火虫算法的移动机器人路径规划方法。
背景技术
路径规划是移动机器人导航中的关键技术之一,其目的是从起始位置到目标位置寻找一条最优的无碰撞路径。随着移动机器人应用范围的不断扩大,使其对路径规划技术的要求越来越高。一些新型人工智能(Artificial Intelligence,AI)技术逐渐应用于路径规划中,尤其是群集智能(Swarm Intelligence,SI)技术。例如,Garcia提出基于蚁群启发式的路径规划方法,并证明了该方法适合于静态和动态环境中的路径规划;Gong提出了一种基于多目标粒子群优化的路径规划方法,利用自适应变异操作提高可行路径的质量;Peng通过改进基于人工鱼群算法的路径规划方法来提高移动机器人在不同环境中路径规划的适应性。萤火虫算法(Firefly Algorithm,FA)与其他生物启发算法相比,它的理论和实现都相对简单,在路径规划中也有一定的应用,如:Li利用萤火虫算法和Bezier曲线来寻找最短无冲突的可行路径,但该算法仍然存在一些问题,如收敛速度慢,计算量大,稳定性差,易陷入局部最优等。因此,基于混沌萤火虫算法的移动机器人路径规划方法在狭窄区域内快速搜所到最优路径有非常重要的意义。本方法利用Lozi’s映射混沌序列调整萤火虫算法中的控制参数,避免了陷入局部最优,同时采用最优调整策略,对萤火虫算法进行高斯干扰,增强了该方法的搜索能力,并将该方法应用于移动机器人路径规划中,解决了移动机器人在狭窄区域内的路径规划问题。
发明内容
本发明旨在解决以上现有技术的问题。提出了一种收敛速度快,不易陷入局部最优的基于混沌萤火虫算法的移动机器人路径规划方法。本发明的技术方案如下:
一种基于混沌萤火虫算法的移动机器人路径规划方法,其包括以下步骤:
S1、获取移动机器人起始点到目标点的可行路径参数,确定移动机器人系统代价函数;S2、对移动机器人系统进行初始化,得到所有可能路径;S3、利用混沌萤火虫算法对步骤S2每一条路径进行迭代更新,包括对亮度和位置的更新;S4、对更新后的路径进行排序,得到局部最优路径;S5、步骤S3迭代次数完成进行步骤S6,否则继续进行步骤S3;S6、利用最优调整策略调整当前局部最优路径;S7、输出全局最优路径。
进一步的,所述步骤S1移动机器人系统代价函数表示为:
E=ω1EL2ED
其中ω12分别表示长度系数和冒险系数,且ω12=1,取值分别为ω1=0.75,ω2=0.25,EL表示移动机器人从起始点到目标点的可行路径距离,ED表示可行路径中各节点和环境中障碍物的代价函数;
其中(xi,yi,zi),表示可行路径中第i个节点的坐标,i=1…n-1,表示移动机器人起始点与可行路径第一节点间的距离,表示移动机器人目标点与可行路径最后一个节点间的距离,ΔLi表示可行路径中每两个相邻节点之间的距离,且
ED=C(k1+k2),其中C表示代价系数,且C>0,k1表示可行路径中节点个数,k2表示环境中障碍物个数。
进一步的,所述步骤S3利用混沌萤火虫算法对步骤S2每一条路径进行迭代更新包括:
S31、采用步骤S1得到的代价函数表示萤火虫的绝对亮度,代价函数值表示萤火虫所在位置处潜在的解,即每个萤火虫代表一条可行的路径,萤火虫的数量代表可行路径的数量,萤火虫的绝对亮度代表路径的质量;
S32、绝对亮度小的萤火虫向绝对亮度大的萤火虫移动,这里的绝对亮度是指萤火虫在当前位置的亮度;萤火虫相对亮度是指两个萤火虫之间相对的亮度,表示为其中I0表示最大亮度,分别进行萤火虫亮度和位置的更新;
萤火虫相对亮度与吸引力成正比,吸引力表示为:
萤火虫位置更新公式表示为:
xi=xiij(ij)(xi-xj)+α(t)(rand+0.5)
其中t为迭代次数,xi,xj为萤火虫i和萤火虫j所处的空间位置,β0为萤火虫的最大吸引力,rij为两个萤火虫之间的距离,rand表示为[0,1]之间的均匀分布;
γ(t)为Lozi’s映射混沌序列表示的光吸收系数,表示为:
α(t)为Lozi’s映射混沌序列表示的随机参数,表示为:
其中,t为迭代次数,a1,b1为控制γ(t)变化的控制参数,a2,b2为控制α(t)变化的控制参数。
进一步的,所述γ(1)=0.7(t=1),α(1)=1(t=1),a1=a2=1.7,b1=b2=0.5。
进一步的,所述步骤S4对更新后的路径进行排序,得到局部最优路径采用的排序方法包括比较法、冒泡法。
进一步的,所述步骤S6的最优调整策略为:
用高斯分布对每次迭代的最优路径的位置进行微小扰动,即:
其中,为高斯扰动后的最优位置,xbest为当前最优位置,η为控制参数,N(0,1)为高斯分布。
本发明的优点及有益效果如下:
本发明是一种基于混沌萤火虫算法的移动机器人路径规划方法,利用Lozi’s映射混沌序列调整萤火虫算法中的控制参数,避免了陷入局部最优,提高了收敛速度;采用最优调整策略,对萤火虫算法进行高斯干扰,增强了该算法的搜索能力。并将本方法应用于移动机器人路径规划中,使其在较短时间内规划出一条通过狭窄区域的最优路径。
附图说明
图1是本发明提供优选实施例基于混沌萤火虫算法的移动机器人路径规划方法流程图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、详细地描述。所描述的实施例仅仅是本发明的一部分实施例。
本发明的技术方案如下:
本发明提供了一种基于混沌萤火虫算法的移动机器人路径规划方法,其特征在于,包括以下步骤:
S1,确定移动机器人系统代价函数。在本实施方式中,移动机器人系统代价函数表示为:
E=ω1EL2ED
其中ω12分别表示长度系数和冒险系数,且ω12=1,取值分别为ω1=0.75,ω2=0.25,EL表示移动机器人从起始点到目标点的可行路径距离,ED表示可行路径中各节点和环境中障碍物的代价函数。
其中(xi,yi,zi),i=1…n-1,表示可行路径中第i个节点的坐标,表示移动机器人起始点与可行路径第一节点间的距离,表示移动机器人目标点与可行路径最后一个节点间的距离,ΔLi表示可行路径中每两个相邻节点之间的距离,且
ED=C(k1+k2),其中C表示代价系数,且C>0,k1表示可行路径中节点个数,k2表示环境中障碍物个数。
S2,对移动机器人系统进行初始化,得到所有可能路径。
S3,利用混沌萤火虫方法对每一条路径进行迭代更新,包括对亮度和位置的更新。在本实施方式中,混沌萤火虫方法为:
用代价函数表示萤火虫的绝对亮度,代价函数值表示萤火虫所在位置处潜在的解,即每个萤火虫代表一条可行的路径,萤火虫的数量代表可行路径的数量,萤火虫的绝对亮度代表路径的质量。
绝对亮度小的萤火虫向绝对亮度大的萤火虫移动,分别进行萤火虫亮度和位置的更新。萤火虫相对亮度是指两个萤火虫之间相对的亮度,表示为其中I0表示最大亮度,绝对亮度是指萤火虫在当前位置的亮度。
萤火虫相对亮度与吸引力成正比。吸引力表示为:
萤火虫位置更新公式表示为:
xi=xiij(ij)(xi-xj)+α(t)(rand+0.5)
其中t为迭代次数,xi,xj为萤火虫i和萤火虫j所处的空间位置,β0为萤火虫的最大吸引力,rij为两个萤火虫之间的距离,rand表示为[0,1]之间的均匀分布。
γ(t)为Lozi’s映射混沌序列表示的光吸收系数,表示为:
α(t)为Lozi’s映射混沌序列表示的随机参数,表示为:
其中,t为迭代次数,a1,b1为控制γ(t)变化的控制参数,a2,b2为控制α(t)变化的控制参数,且γ(1)=0.7(t=1),α(1)=1(t=1),a1=a2=1.7,b1=b2=0.5。
S4,对更新后的路径进行排序,得到局部最优路径采用的排序方法包括比较法、冒泡法。
S5,迭代次数完成进行步骤S6,否则进行步骤S3。
S6,利用最优调整策略调整当前局部最优路径。在本实施方式中,最优调整策略为:
用高斯分布对每次迭代的最优路径的位置进行微小扰动,即:
其中,为高斯扰动后的最优位置,xbest为当前最优位置,η为控制参数,N(0,1)为高斯分布。
S7,输出全局最优路径。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不一定指的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
以上这些实施例应理解为仅用于说明本发明而不用于限制本发明的保护范围。在阅读了本发明的记载的内容之后,技术人员可以对本发明作各种改动或修改,这些等效变化和修饰同样落入本发明权利要求所限定的范围。

Claims (5)

1.一种基于混沌萤火虫算法的移动机器人路径规划方法,其特征在于,包括以下步骤:
S1、获取移动机器人起始点到目标点的可行路径参数,确定移动机器人系统代价函数;S2、对移动机器人系统进行初始化,得到所有可能路径;S3、利用混沌萤火虫算法对步骤S2每一条路径进行迭代更新,包括对亮度和位置的更新;所述步骤S3利用混沌萤火虫算法对步骤S2每一条路径进行迭代更新包括:
S31、采用步骤S1得到的代价函数表示萤火虫的绝对亮度,代价函数值表示萤火虫所在位置处潜在的解,即每个萤火虫代表一条可行的路径,萤火虫的数量代表可行路径的数量,萤火虫的绝对亮度代表路径的质量;
S32、绝对亮度小的萤火虫向绝对亮度大的萤火虫移动,这里的绝对亮度是指萤火虫在当前位置的亮度;萤火虫相对亮度是指两个萤火虫之间相对的亮度,表示为其中I0表示最大亮度,分别进行萤火虫亮度和位置的更新;
萤火虫相对亮度与吸引力成正比,吸引力表示为:
萤火虫位置更新公式表示为:
xi=xiij(ij)(xi-xj)+α(t)(rand+0.5)
其中t为迭代次数,xi,xj为萤火虫i和萤火虫j所处的空间位置,β0为萤火虫的最大吸引力,rij为两个萤火虫之间的距离,rand表示为[0,1]之间的均匀分布;
γ(t)为Lozi’s映射混沌序列表示的光吸收系数,表示为:
α(t)为Lozi’s映射混沌序列表示的随机参数,表示为:
其中,t为迭代次数,a1,b1为控制γ(t)变化的控制参数,a2,b2为控制α(t)变化的控制参数;S4、对更新后的路径进行排序,得到局部最优路径;S5、步骤S3迭代次数完成进行步骤S6,否则继续进行步骤S3;S6、利用最优调整策略调整当前局部最优路径;S7、输出全局最优路径。
2.根据权利要求1所述的基于混沌萤火虫算法的移动机器人路径规划方法,其特征在于,所述步骤S1移动机器人系统代价函数表示为:
E=ω1EL2ED
其中ω12分别表示长度系数和冒险系数,且ω12=1,取值分别为ω1=0.75,ω2=0.25,EL表示移动机器人从起始点到目标点的可行路径距离,ED表示可行路径中各节点和环境中障碍物的代价函数;
其中(xi,yi,zi),表示可行路径中第i个节点的坐标,i=1…n-1,LSP1表示移动机器人起始点与可行路径第一节点间的距离,LPnE表示移动机器人目标点与可行路径最后一个节点间的距离,ΔLi表示可行路径中每两个相邻节点之间的距离,且
ED=C(k1+k2),其中C表示代价系数,且C>0,k1表示可行路径中节点个数,k2表示环境中障碍物个数。
3.根据权利要求1所述的基于混沌萤火虫算法的移动机器人路径规划方法,其特征在于,所述γ(1)=0.7(t=1),α(1)=1(t=1),a1=a2=1.7,b1=b2=0.5。
4.根据权利要求1所述的基于混沌萤火虫算法的移动机器人路径规划方法,其特征在于,所述步骤S4对更新后的路径进行排序,得到局部最优路径采用的排序方法包括比较法、冒泡法。
5.根据权利要求1所述的基于混沌萤火虫算法的移动机器人路径规划方法,其特征在于,所述步骤S6的最优调整策略为:
用高斯分布对每次迭代的最优路径的位置进行微小扰动,即:
其中,为高斯扰动后的最优位置,xbest为当前最优位置,η为控制参数,N(0,1)为高斯分布。
CN201610569272.2A 2016-07-19 2016-07-19 一种基于混沌萤火虫算法的移动机器人路径规划方法 Active CN106094833B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610569272.2A CN106094833B (zh) 2016-07-19 2016-07-19 一种基于混沌萤火虫算法的移动机器人路径规划方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610569272.2A CN106094833B (zh) 2016-07-19 2016-07-19 一种基于混沌萤火虫算法的移动机器人路径规划方法

Publications (2)

Publication Number Publication Date
CN106094833A CN106094833A (zh) 2016-11-09
CN106094833B true CN106094833B (zh) 2018-12-28

Family

ID=57220758

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610569272.2A Active CN106094833B (zh) 2016-07-19 2016-07-19 一种基于混沌萤火虫算法的移动机器人路径规划方法

Country Status (1)

Country Link
CN (1) CN106094833B (zh)

Families Citing this family (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106406320B (zh) * 2016-11-29 2019-08-20 重庆重智机器人研究院有限公司 机器人路径规划方法及规划路线的机器人
CN107423838A (zh) * 2017-04-16 2017-12-01 江西理工大学 基于混沌引力搜索的车辆路径规划方法
CN107219858B (zh) * 2017-04-26 2020-04-03 西北工业大学 一种改进萤火虫算法的多无人机协同耦合任务分配方法
CN108255178A (zh) * 2018-01-18 2018-07-06 山东理工大学 一种巡逻机器人全覆盖遍历路径规划算法的设计方法
CN109350510B (zh) * 2018-08-27 2021-06-22 重庆市中医院 一种自动火针治疗仪及其操作方法
CN109297498A (zh) * 2018-10-11 2019-02-01 南昌工程学院 一种基于改进多目标萤火虫算法的机器人路径规划方法
CN109542106A (zh) * 2019-01-04 2019-03-29 电子科技大学 一种移动机器人多约束条件下的路径规划方法
CN109613921B (zh) * 2019-01-08 2021-07-30 智慧航海(青岛)科技有限公司 基于快速移动萤火虫算法的无人船舶局部路径规划方法
US11353878B2 (en) * 2019-03-26 2022-06-07 Baidu Usa Llc Soft-boundary based path optimization for complex scenes for autonomous driving vehicles
CN110389587A (zh) * 2019-05-20 2019-10-29 长沙理工大学 一种目标点动态变化的机器人路径规划新方法
CN110823229B (zh) * 2019-12-16 2023-07-07 湖北工业大学 一种基于优化萤火虫算法的移动机器人路径规划方法系统
CN111552296B (zh) * 2020-05-14 2021-03-26 宁波智能装备研究院有限公司 一种基于弯曲柱坐标系的局部平滑轨迹规划方法
CN112101655B (zh) * 2020-09-11 2023-04-28 长沙理工大学 萤火虫算法在路径规划应用中的种群动态调节策略与实现
CN112731961A (zh) * 2020-12-08 2021-04-30 深圳供电局有限公司 路径规划方法、装置、设备及存储介质
CN112550497B (zh) * 2021-03-01 2021-06-01 苏州澳昆智能机器人技术有限公司 一种自动装车机器人用行走组件
CN115639828B (zh) * 2022-12-23 2023-03-31 河北科技大学 基于改进人工鱼群算法的机器人路径规划方法及装置

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102360214A (zh) * 2011-09-02 2012-02-22 哈尔滨工程大学 一种基于萤火虫算法的舰船路径规划方法
CN102768536A (zh) * 2012-07-20 2012-11-07 哈尔滨工程大学 一种基于多目标萤火虫算法的路径规划方法
CN103196449A (zh) * 2013-03-28 2013-07-10 哈尔滨工程大学 基于潮流潮汐预报信息的船舶航路规划方法
CN103968841A (zh) * 2014-06-03 2014-08-06 哈尔滨工程大学 一种基于改进萤火虫算法的auv三维航路规划方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2011055145A1 (en) * 2009-11-05 2011-05-12 Bae Systems Plc Generating a set of solutions to a multi-objective problem

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102360214A (zh) * 2011-09-02 2012-02-22 哈尔滨工程大学 一种基于萤火虫算法的舰船路径规划方法
CN102768536A (zh) * 2012-07-20 2012-11-07 哈尔滨工程大学 一种基于多目标萤火虫算法的路径规划方法
CN103196449A (zh) * 2013-03-28 2013-07-10 哈尔滨工程大学 基于潮流潮汐预报信息的船舶航路规划方法
CN103968841A (zh) * 2014-06-03 2014-08-06 哈尔滨工程大学 一种基于改进萤火虫算法的auv三维航路规划方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
A New Path Planning Method Based on Firefly Algorithm;Liu Chang 等;《2012 Fifth International Joint Conference on Computational Sciences and Optimization》;20120823;第775-778页 *
Chaotic Bee Swarm Optimization Algorithm for Path Planning of Mobile Robots;Jiann-Horng Lin 等;《Proceedings of the 10th WSEAS International Conference on EVOLUTIONARY COMPUTING》;20090325;第84-89页 *
基于改进萤火虫算法的动态自动聚集路径规划方法;刘鹏 等;《计算机应用研究》;20111130;第28卷(第11期);第4146-4149页 *
基于混沌理论的动态种群萤火虫算法;冯艳红 等;《计算机应用》;20130301;第796-799、805页 *
带高斯变异的人工萤火虫优化算法;莫愿斌 等;《计算机应用研究》;20130131;第30卷(第1期);第121-123页 *

Also Published As

Publication number Publication date
CN106094833A (zh) 2016-11-09

Similar Documents

Publication Publication Date Title
CN106094833B (zh) 一种基于混沌萤火虫算法的移动机器人路径规划方法
CN107368076B (zh) 一种智能环境下机器人运动路径深度学习控制规划方法
WO2016045615A1 (zh) 机器人静态路径规划方法
CN107607120B (zh) 基于改进修复式Anytime稀疏A*算法的无人机动态航迹规划方法
CN107402018B (zh) 一种基于连续帧的导盲仪组合路径规划方法
Ali et al. Cooperative path planning of multiple UAVs by using max–min ant colony optimization along with cauchy mutant operator
CN111024092A (zh) 一种多约束条件下智能飞行器航迹快速规划方法
CN105426992B (zh) 移动机器人旅行商优化方法
CN109542106A (zh) 一种移动机器人多约束条件下的路径规划方法
CN106525047A (zh) 一种基于floyd算法的无人机路径规划方法
CN109186619B (zh) 一种基于实时路况的智能导航算法
JP2009156617A (ja) ナビゲーション装置およびプログラム
CN109357678A (zh) 一种基于异质化鸽群优化算法的多无人机路径规划方法
CN106052692A (zh) 一种最短路径规划导航方法及系统
CN111381600A (zh) 一种基于粒子群算法的uuv路径规划方法
CN105869512B (zh) 多信息的混杂度量地图建图方法和装置
CN109947132A (zh) 基于rfid库存盘点的无人机三维路径规划方法
CN116700329A (zh) 一种基于改进萤火虫算法的无人机航迹规划方法
Yu et al. Balanced multi-region coverage path planning for unmanned aerial vehicles
CN108645505B (zh) 一种随机共振微弱信号检测方法
CN111797331B (zh) 基于群智感知的多目标多约束路线推荐方法
CN109145462A (zh) 一种基于力导引的异形区域拓扑布局算法
CN112484733B (zh) 一种基于拓扑图的强化学习室内导航方法
CN110187306A (zh) 一种应用于复杂室内空间的tdoa-pdr-map融合定位方法
Zhou et al. Routing optimization of intelligent vehicle in automated warehouse

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant