CN101604166A - 一种基于粒子群优化算法的移动机器人路径规划方法 - Google Patents

一种基于粒子群优化算法的移动机器人路径规划方法 Download PDF

Info

Publication number
CN101604166A
CN101604166A CNA2009101006131A CN200910100613A CN101604166A CN 101604166 A CN101604166 A CN 101604166A CN A2009101006131 A CNA2009101006131 A CN A2009101006131A CN 200910100613 A CN200910100613 A CN 200910100613A CN 101604166 A CN101604166 A CN 101604166A
Authority
CN
China
Prior art keywords
path
robot
point
circle
radius
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.)
Granted
Application number
CNA2009101006131A
Other languages
English (en)
Other versions
CN101604166B (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.)
Haian Service Center for Transformation of Scientific Achievements
Original Assignee
Hangzhou Electronic Science and Technology University
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 Hangzhou Electronic Science and Technology University filed Critical Hangzhou Electronic Science and Technology University
Priority to CN2009101006131A priority Critical patent/CN101604166B/zh
Publication of CN101604166A publication Critical patent/CN101604166A/zh
Application granted granted Critical
Publication of CN101604166B publication Critical patent/CN101604166B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
  • Feedback Control In General (AREA)

Abstract

本发明涉及一种基于粒子群优化算法的移动机器人路径规划方法。目前移动机器人路径规划方法对于路径的选择比较单一、容易走入死胡同。本发明的具体步骤是:首先对环境进行建模;用多边形表示障碍物,用点表示机器人,用小方框表示机器人的起始位置,用十字表示目标点。机器人开始点标记为S,目标点标记为G;然后对机器人路径利用粒子群优化算法进行规划;最后对规划的路径进行深度优先搜索。本发明在极坐标粒子群中加入了深度优先搜索。用这种方法,机器人可以在大多数复杂环境中有效地找到一个无碰撞的路径,最终到达目标点。

Description

一种基于粒子群优化算法的移动机器人路径规划方法
技术领域
本发明计算机技术领域,具体是涉及一种基于粒子群优化算法的移动机器人路径规划方法。
背景技术
移动机器人是智能控制技术中的一个重要领域,已经被广泛应用在军事、工业、农业和教育等领域。路径规划是移动机器人系统的一个重要内容,它直接影响机器人完成任务的质量。粒子群优化算法(Particle SwarmOptimization,PSO)由Kennedy和Eberhart首先提出,是一种智能演化算法。把粒子群优化算法(Particle Swarm Optimization,PSO)用在路径规划中是一个新的尝试。通过改进的粒子群优化算法可以使机器人的行走路径变得光滑。在粒子群优化算法(PARTICLE SWARM OPTIMIZATION,PSO)的适应值函数中加入了障碍物躲避参数,使得机器人可以有效地避开障碍物。这些都是粒子群优化算法(Particle Swarm Optimization,PSO)在全局路径规划中的应用,改变适应值函数使得机器人在局部路径规划中可以成功地躲避障碍物而到达目标点。然而这个方法只能用在简单的环境中,不能帮助机器人躲避环境陷阱,如无路径必须返回。于是提出了一种基于极坐标粒子群(Polar Particle Swarm  Optimization,PPSO)的路径规划方法。虽然这个方法可以帮助机器人逃避环境陷阱,但它在障碍物太大的时候就会失效。
发明内容
本发明的目的在于针对现有技术的不足,提供一种基于粒子群优化算法的移动机器人路径规划方法。
本发明方法的具体步骤是:
步骤(1)设机器人的当前位置为Si,机器人传感器感应圆的半径为R,Si与目标点之间的距离为D;以Si为圆心做圆C,如果R≤D,则圆C的半径为R,如果R>D,则圆C的半径为D;以Si为射线出发点,做平行于x轴的射线,射线与圆C相交于Q点,并将此射线作为极坐标轴;
步骤(2)将Si点与Q点间线段进行K等分,得到长度为||SiQ||/K的K段线段,以各个等分点到Si点的距离为半径,以Si为圆心做圆,得到K个同心圆;
步骤(3)可行路径是一条从起点Si到圆C上点的若干线段组成的折线,用这些折线的端点即路径关键点的序列表示,机器人移动过的路径就是这些点的连线,定义为P={Si,P1,...,Pj,...,Pm},1≤j≤m,m<K;其中P1,...,Pj,...,Pm都是非障碍物点,即相邻的两点之间的连线上没有障碍物;
步骤(4)在当前非障碍物点上将传感器感应到的半径为R的圆等分为n个扇形区域{d1,d2,...,dj,dj+1,...,dn},每个扇形的角度为360/n度;
步骤(5)每次选择一个扇形区域作为探索区域:扇形区域dj被选择,用粒子群优化算法探索dj,如果在扇形区域dj中找到可行路径,移动机器人至下一个非障碍物点;如果找不到路径,探索下一个扇形区域dj+1;对每一个扇形区域探索利用PSO算法对机器人路径进行规划,最终找到无障碍物的路径;
步骤(6)如果在传感器感应到的半径为R的圆的范围内找不到可行路径,则返回到前一点,探索在前一点中没有探索过的扇形区域;把机器人每次走过的路径存储在一个链表中,以方便机器人搜索;
步骤(7)当圆C的半径为D时,目标点即在圆周C上,机器人利用PSO算法所做的移动即可到达目标点。当圆C的半径小于D时,机器人移动pm到后,将pm的极坐标值赋值给Si,重复步骤(1)到(7),直到找到一条无碰撞路径,到达目标点。
经过基于极坐标粒子群的深度优先搜索找到一个无碰撞的路径,最终到达目标点。
本发明在极坐标粒子群中加入了深度优先搜索。用这种方法,机器人可以在大多数复杂环境中有效地找到一个无碰撞的路径,最终到达目标点。
具体实施方式
采用深度优先搜索算法和粒子群优化算法相结合的方法,对移动机器人在未知环境信息下进行路径规划。旨在移动机器人能用粒子群优化算法在未知的复杂环境中找到有效路径。
该基于粒子群优化算法的移动机器人路径规划方法的具体步骤是:
步骤1:对环境进行建模。
用多边形表示障碍物,用点表示机器人,用小方框表示机器人的起始位置,用十字表示目标点。
步骤2:对机器人路径进行规划。
(1)假设机器人的当前位置为Si,以Si为圆点做圆C,如果R≤D,则C的半径为R,如果R>D,则C的半径为D。其中R表示机器人传感器感应圆的半径,D表示Si与目标点的距离。以Si为射线出发点,做平行于x轴的射线,射线与圆C相交于Q点,并将此射线作为极坐标轴。
(2)将SiQ点进行K等分,得到间隔长度为||SiQ||/K的K段线段,以各个等分点到原点的连线为半径,以Si为圆心做圆,得到一系列的同心圆。
(3)可行路径是一条从起点Si,到圆C上点的若干线段组成的折线,用这些折线的端点即路径关键点的序列表示。机器人移动过的路径就是这些点的连线,定义为P={S,P1,...,Pj,...,Pm},1≤j≤m。其中Pj是非障碍物点,和Pj相连的点的连线上没有障碍物。
(4)当圆C的半径为D时,目标点即在圆周C上,机器人利用PSO算法所做的移动即可到达目标点。当圆C的半径小于D时,机器人移动pm到后,将pm的极坐标值赋值给Si,重复步骤1)到4),直到找到一条无碰撞路径,到达目标点。
步骤3:基于极坐标粒子群的深度优先搜索。
(1)记录机器人移动过的点记为{s0,s1,s2,...si}。假设机器人的当前位置为Si,传感器感应到的范围为圆Ci,将圆Ci平分为n个扇形{d1,d2,...,dj,dj+1,...,dn},每个扇形的角度为360/n度。
(2)每次选择一个扇形作为探索区域。假设dj被选择,用粒子群优化算法PSO探索dj,如果在dj中找到可行路径,移动机器人;如果找不到路径,探索dj+1部分。对每一个扇形区域探索利用PSO算法对机器人路径进行规划,最终找到无障碍物的路径。
(3)如果在圆Ci范围内找不到可行路径,则返回si-1点,探索在Ci-1中没有探索过的区域(扇形)。把机器人每次走过的路径存储在一个链表中,以方便其搜索。
经过基于极坐标粒子群的深度优先搜索找到一个无碰撞的路径,最终到达目标点。

Claims (1)

1、一种基于粒子群优化算法的移动机器人路径规划方法,其特征在于该方法的步骤包括:
步骤(1)设机器人的当前位置为Si,机器人传感器感应圆的半径为R,Si与目标点之间的距离为D;以Si为圆心做圆C,如果R≤D,则圆C的半径为R,如果R>D,则圆C的半径为D;以Si为射线出发点,做平行于x轴的射线,射线与圆C相交于Q点,并将此射线作为极坐标轴;
步骤(2)将Si点与Q点间线段进行K等分,得到长度为||SiQ||/K的K段线段,以各个等分点到Si点的距离为半径,以Si为圆心做圆,得到K个同心圆;
步骤(3)可行路径是一条从起点Si到圆C上点的若干线段组成的折线,用这些折线的端点即路径关键点的序列表示,机器人移动过的路径就是这些点的连线,定义为P={Si,P1,...,Pj,...,Pm},1≤j≤m,m<K;其中P1,...,Pj,...,Pm都是非障碍物点,即相邻的两点之间的连线上没有障碍物;
步骤(4)在当前非障碍物点上将传感器感应到的半径为R的圆等分为n个扇形区域{d1,d2,...,dj,dj+1,...,dn},每个扇形的角度为360/n度;
步骤(5)每次选择一个扇形区域作为探索区域:扇形区域dj被选择,用粒子群优化算法探索dj,如果在扇形区域dj中找到可行路径,移动机器人至下一个非障碍物点;如果找不到路径,探索下一个扇形区域dj+1;对每一个扇形区域探索利用PSO算法对机器人路径进行规划,最终找到无障碍物的路径;
步骤(6)如果在传感器感应到的半径为R的圆的范围内找不到可行路径,则返回到前一点,探索在前一点中没有探索过的扇形区域;把机器人每次走过的路径存储在一个链表中;
步骤(7)当圆C的半径为D时,目标点即在圆周C上,机器人利用PSO算法所做的移动即可到达目标点;当圆C的半径小于D时,机器人移动pm到后,将pm的极坐标值赋值给Si,重复步骤(1)到(7),直到找到一条无碰撞路径,到达目标点。
CN2009101006131A 2009-07-10 2009-07-10 一种基于粒子群优化算法的移动机器人路径规划方法 Expired - Fee Related CN101604166B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2009101006131A CN101604166B (zh) 2009-07-10 2009-07-10 一种基于粒子群优化算法的移动机器人路径规划方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2009101006131A CN101604166B (zh) 2009-07-10 2009-07-10 一种基于粒子群优化算法的移动机器人路径规划方法

Publications (2)

Publication Number Publication Date
CN101604166A true CN101604166A (zh) 2009-12-16
CN101604166B CN101604166B (zh) 2011-01-12

Family

ID=41469927

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2009101006131A Expired - Fee Related CN101604166B (zh) 2009-07-10 2009-07-10 一种基于粒子群优化算法的移动机器人路径规划方法

Country Status (1)

Country Link
CN (1) CN101604166B (zh)

Cited By (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102129249A (zh) * 2011-01-10 2011-07-20 中国矿业大学 一种危险源环境下的机器人全局路径规划方法
CN102169347A (zh) * 2011-03-08 2011-08-31 浙江工业大学 基于协作协进化和多种群遗传算法的多机器人路径规划系统
CN102175245A (zh) * 2011-01-28 2011-09-07 哈尔滨工程大学 一种基于海流历史统计信息的水下潜器路径规划方法
CN102854880A (zh) * 2012-10-08 2013-01-02 中国矿业大学 面向混合地形区域不确定环境的机器人全局路径规划方法
CN103336526A (zh) * 2013-06-20 2013-10-02 苏州经贸职业技术学院 基于协同进化粒子群滚动优化的机器人路径规划方法
CN103576685A (zh) * 2013-11-14 2014-02-12 哈尔滨工程大学 母船回收时uuv的路径确定方法
CN104050390A (zh) * 2014-06-30 2014-09-17 西南交通大学 一种基于可变维粒子群膜算法的移动机器人路径规划方法
CN104063745A (zh) * 2014-04-21 2014-09-24 河海大学 基于改进粒子群算法的多路径规划方法
CN104517297A (zh) * 2013-09-28 2015-04-15 沈阳新松机器人自动化股份有限公司 一种基于粒子群优化的机器人标定方法
CN104765371A (zh) * 2015-04-22 2015-07-08 福州大学 基于滚动窗的深度搜索与模糊控制融合的路径规划方法
CN105717929A (zh) * 2016-04-29 2016-06-29 中国人民解放军国防科学技术大学 一种多分辨率障碍物环境下移动机器人混合路径规划方法
CN106897800A (zh) * 2017-02-28 2017-06-27 东方网力科技股份有限公司 一种目标的运动范围预测方法及装置
CN107148064A (zh) * 2017-05-15 2017-09-08 南京邮电大学 一种基于粒子群的车载网络路由路径优化方法
CN107450563A (zh) * 2017-09-21 2017-12-08 景德镇陶瓷大学 一种基于多子群的自适应信息反馈粒子群机器人路径选择方法及系统
CN107843252A (zh) * 2017-10-18 2018-03-27 歌尔股份有限公司 导航路径优化方法、装置及电子设备
CN108133471A (zh) * 2016-11-30 2018-06-08 天津职业技术师范大学 一种自然光照条件下基于人工蜂群算法的农业移动机器人导航路径提取方法及装置
WO2018176596A1 (zh) * 2017-03-31 2018-10-04 深圳市靖洲科技有限公司 基于权重改进粒子群算法的无人自行车路径规划方法
CN108827314A (zh) * 2018-08-16 2018-11-16 电子科技大学 一种路径规划方法
CN108917770A (zh) * 2018-07-25 2018-11-30 重庆交通大学 一种工业机器人路径搜索优化方法
CN109394251A (zh) * 2017-08-15 2019-03-01 西门子保健有限责任公司 用于运行x射线设备的方法和具有铰接臂的x射线设备
CN109581987A (zh) * 2018-12-29 2019-04-05 广东飞库科技有限公司 一种基于粒子群算法的agv调度路径规划方法及系统
CN109782779A (zh) * 2019-03-19 2019-05-21 电子科技大学 基于种群超启发式算法的洋流环境下auv路径规划方法
CN113043286A (zh) * 2020-12-25 2021-06-29 杭州电子科技大学 一种多自由度机械臂实时避障路径规划系统及方法
CN113137969A (zh) * 2021-05-13 2021-07-20 洛阳尚奇机器人科技有限公司 一种移动机器人的局部路径优化方法
CN114237244A (zh) * 2021-12-15 2022-03-25 东南大学 一种基于光线跟踪的移动机器人路径规划方法

Cited By (39)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102129249A (zh) * 2011-01-10 2011-07-20 中国矿业大学 一种危险源环境下的机器人全局路径规划方法
CN102175245A (zh) * 2011-01-28 2011-09-07 哈尔滨工程大学 一种基于海流历史统计信息的水下潜器路径规划方法
CN102169347A (zh) * 2011-03-08 2011-08-31 浙江工业大学 基于协作协进化和多种群遗传算法的多机器人路径规划系统
CN102854880A (zh) * 2012-10-08 2013-01-02 中国矿业大学 面向混合地形区域不确定环境的机器人全局路径规划方法
CN102854880B (zh) * 2012-10-08 2014-12-31 中国矿业大学 面向混合地形区域不确定环境的机器人全局路径规划方法
CN103336526B (zh) * 2013-06-20 2015-08-05 苏州经贸职业技术学院 基于协同进化粒子群滚动优化的机器人路径规划方法
CN103336526A (zh) * 2013-06-20 2013-10-02 苏州经贸职业技术学院 基于协同进化粒子群滚动优化的机器人路径规划方法
CN104517297A (zh) * 2013-09-28 2015-04-15 沈阳新松机器人自动化股份有限公司 一种基于粒子群优化的机器人标定方法
CN103576685A (zh) * 2013-11-14 2014-02-12 哈尔滨工程大学 母船回收时uuv的路径确定方法
CN104063745A (zh) * 2014-04-21 2014-09-24 河海大学 基于改进粒子群算法的多路径规划方法
CN104063745B (zh) * 2014-04-21 2017-05-03 河海大学 基于改进粒子群算法的多路径规划方法
CN104050390A (zh) * 2014-06-30 2014-09-17 西南交通大学 一种基于可变维粒子群膜算法的移动机器人路径规划方法
CN104050390B (zh) * 2014-06-30 2017-05-17 西南交通大学 一种基于可变维粒子群膜算法的移动机器人路径规划方法
CN104765371A (zh) * 2015-04-22 2015-07-08 福州大学 基于滚动窗的深度搜索与模糊控制融合的路径规划方法
CN105717929A (zh) * 2016-04-29 2016-06-29 中国人民解放军国防科学技术大学 一种多分辨率障碍物环境下移动机器人混合路径规划方法
CN105717929B (zh) * 2016-04-29 2018-06-15 中国人民解放军国防科学技术大学 一种多分辨率障碍物环境下移动机器人混合路径规划方法
CN108133471A (zh) * 2016-11-30 2018-06-08 天津职业技术师范大学 一种自然光照条件下基于人工蜂群算法的农业移动机器人导航路径提取方法及装置
CN108133471B (zh) * 2016-11-30 2021-09-17 天津职业技术师范大学 一种基于人工蜂群算法的机器人导航路径提取方法及装置
CN106897800A (zh) * 2017-02-28 2017-06-27 东方网力科技股份有限公司 一种目标的运动范围预测方法及装置
WO2018176596A1 (zh) * 2017-03-31 2018-10-04 深圳市靖洲科技有限公司 基于权重改进粒子群算法的无人自行车路径规划方法
CN107148064B (zh) * 2017-05-15 2020-06-23 南京邮电大学 一种基于粒子群的车载网络路由路径优化方法
CN107148064A (zh) * 2017-05-15 2017-09-08 南京邮电大学 一种基于粒子群的车载网络路由路径优化方法
CN109394251A (zh) * 2017-08-15 2019-03-01 西门子保健有限责任公司 用于运行x射线设备的方法和具有铰接臂的x射线设备
CN107450563A (zh) * 2017-09-21 2017-12-08 景德镇陶瓷大学 一种基于多子群的自适应信息反馈粒子群机器人路径选择方法及系统
CN107450563B (zh) * 2017-09-21 2020-08-25 景德镇陶瓷大学 基于多子群的自适应信息反馈粒子群机器人路径选择方法
CN107843252A (zh) * 2017-10-18 2018-03-27 歌尔股份有限公司 导航路径优化方法、装置及电子设备
CN107843252B (zh) * 2017-10-18 2020-10-09 歌尔股份有限公司 导航路径优化方法、装置及电子设备
CN108917770A (zh) * 2018-07-25 2018-11-30 重庆交通大学 一种工业机器人路径搜索优化方法
CN108827314A (zh) * 2018-08-16 2018-11-16 电子科技大学 一种路径规划方法
CN108827314B (zh) * 2018-08-16 2020-06-16 电子科技大学 一种路径规划方法
CN109581987A (zh) * 2018-12-29 2019-04-05 广东飞库科技有限公司 一种基于粒子群算法的agv调度路径规划方法及系统
CN109782779B (zh) * 2019-03-19 2020-11-06 电子科技大学 基于种群超启发式算法的洋流环境下auv路径规划方法
CN109782779A (zh) * 2019-03-19 2019-05-21 电子科技大学 基于种群超启发式算法的洋流环境下auv路径规划方法
CN113043286A (zh) * 2020-12-25 2021-06-29 杭州电子科技大学 一种多自由度机械臂实时避障路径规划系统及方法
CN113043286B (zh) * 2020-12-25 2022-05-03 杭州电子科技大学 一种多自由度机械臂实时避障路径规划系统及方法
CN113137969A (zh) * 2021-05-13 2021-07-20 洛阳尚奇机器人科技有限公司 一种移动机器人的局部路径优化方法
CN113137969B (zh) * 2021-05-13 2023-08-04 洛阳尚奇机器人科技有限公司 一种移动机器人的局部路径优化方法
CN114237244A (zh) * 2021-12-15 2022-03-25 东南大学 一种基于光线跟踪的移动机器人路径规划方法
CN114237244B (zh) * 2021-12-15 2023-11-03 东南大学 一种基于光线跟踪的移动机器人路径规划方法

Also Published As

Publication number Publication date
CN101604166B (zh) 2011-01-12

Similar Documents

Publication Publication Date Title
CN101604166B (zh) 一种基于粒子群优化算法的移动机器人路径规划方法
CN106444769B (zh) 一种室内移动机器人增量式环境信息采样的最优路径规划方法
CN104298239B (zh) 一种室内移动机器人增强地图学习路径规划方法
CN107894773A (zh) 一种移动机器人的导航方法、系统及相关装置
CN111024092B (zh) 一种多约束条件下智能飞行器航迹快速规划方法
CN109726866A (zh) 基于q学习神经网络的无人船路径规划方法
Fang et al. Autonomous robotic exploration based on frontier point optimization and multistep path planning
Shi et al. A local obstacle avoidance method for mobile robots in partially known environment
CN106873599A (zh) 基于蚁群算法和极坐标变换的无人自行车路径规划方法
CN109933067A (zh) 一种基于遗传算法和粒子群算法的无人艇避碰方法
Zhang et al. Online path planning for UAV using an improved differential evolution algorithm
CN113625716B (zh) 一种多智能体动态路径规划方法
CN102541057A (zh) 一种基于激光测距仪的移动机器人避障方法
CN112925315A (zh) 一种基于改进蚁群算法和a*算法的履带车路径规划方法
CN112947594B (zh) 一种面向无人机的航迹规划方法
Ma et al. A two-level path planning method for on-road autonomous driving
Chen et al. An improved path planning method based on artificial potential field for a mobile robot
CN103869824A (zh) 基于生物触角模型的多机器人水下目标搜寻方法及装置
CN109445440A (zh) 基于传感器融合与改进q学习算法的动态避障方法
CN109189074A (zh) 一种用于仓储环境的室内自主建图方法
CN104765371A (zh) 基于滚动窗的深度搜索与模糊控制融合的路径规划方法
Gyenes et al. Motion planning for mobile robots using the safety velocity obstacles method
CN112857370A (zh) 一种基于时序信息建模的机器人无地图导航方法
CN113219990B (zh) 基于自适应邻域与转向代价的机器人路径规划方法
CN114564008A (zh) 基于改进A-Star算法的移动机器人路径规划方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
ASS Succession or assignment of patent right

Owner name: HAIAN SERVICE CENTER FOR TRANSFORMATION OF SCIENTI

Free format text: FORMER OWNER: HANGZHOU ELECTRONIC SCIENCE AND TECHNOLOGY UNIV

Effective date: 20131021

C41 Transfer of patent application or patent right or utility model
COR Change of bibliographic data

Free format text: CORRECT: ADDRESS; FROM: 310018 HANGZHOU, ZHEJIANG PROVINCE TO: 226600 NANTONG, JIANGSU PROVINCE

TR01 Transfer of patent right

Effective date of registration: 20131021

Address after: 226600 No. 106 middle Yangtze Road, Haian County, Nantong, Jiangsu

Patentee after: Haian Service Center for Transformation of Scientific Achievements

Address before: Hangzhou City, Zhejiang province 310018 Jianggan District Xiasha Higher Education Park No. 2 street

Patentee before: Hangzhou Electronic Science and Technology Univ

CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20110112

Termination date: 20180710

CF01 Termination of patent right due to non-payment of annual fee