JP2014211760A - 経路修正装置 - Google Patents

経路修正装置 Download PDF

Info

Publication number
JP2014211760A
JP2014211760A JP2013087749A JP2013087749A JP2014211760A JP 2014211760 A JP2014211760 A JP 2014211760A JP 2013087749 A JP2013087749 A JP 2013087749A JP 2013087749 A JP2013087749 A JP 2013087749A JP 2014211760 A JP2014211760 A JP 2014211760A
Authority
JP
Japan
Prior art keywords
route
evaluation value
point interval
path
distance
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
JP2013087749A
Other languages
English (en)
Other versions
JP5986533B2 (ja
Inventor
亮暢 藤井
Akinobu Fujii
亮暢 藤井
豊 高岡
Yutaka Takaoka
豊 高岡
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.)
Toyota Motor Corp
Toyota Central R&D Labs Inc
Original Assignee
Toyota Motor Corp
Toyota Central R&D Labs Inc
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 Toyota Motor Corp, Toyota Central R&D Labs Inc filed Critical Toyota Motor Corp
Priority to JP2013087749A priority Critical patent/JP5986533B2/ja
Publication of JP2014211760A publication Critical patent/JP2014211760A/ja
Application granted granted Critical
Publication of JP5986533B2 publication Critical patent/JP5986533B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course or altitude of land, water, air, or space vehicles, e.g. automatic pilot
    • G05D1/02Control of position or course in two dimensions
    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G1/00Traffic control systems for road vehicles
    • G08G1/16Anti-collision systems

Landscapes

  • Navigation (AREA)
  • Traffic Control Systems (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Automation & Control Theory (AREA)

Abstract

【課題】経路を複数の離散した経路点で記述したときに、経路点間の距離が適切な距離とされた経路データに修正する。【解決手段】 経路修正装置は、入力された経路データから、経路データで記述される経路を評価したトータル評価値を計算するトータル評価関数が極値をとるときの経路データを探索する。経路データは、経路を複数の離散した経路点で記述するデータを含んでいる。経路修正装置は、隣接する経路点の経路点間隔に関する経路点間隔評価値を計算する処理と、移動体と障害物の離反距離に関する離反距離評価値を計算する処理と、少なくとも経路点間隔評価値と離反距離評価値とからトータル評価値を計算する処理を実行する。そして、経路点間隔が予め設定された設定範囲内から外れると、無限となる経路点間隔評価値を算出する。【選択図】図1

Description

本明細書では、移動体が移動するべき経路の初期値が与えられたときに、より好ましい経路に修正する技術を開示する。
移動体が移動するべき初期経路が与えられたときに、その初期経路をより好ましい経路に修正する技術が開発されている(例えば、非特許文献1)。非特許文献1には、経路を記述する経路データ(各種パラメータ)から、その経路の好ましさの程度を示す評価値を計算する評価関数を求め、その評価値を極小化する経路データ(各種パラメータ)を求める。これによって、初期経路がより好ましい経路に修正される。
車輪移動ロボットの形状と運動拘束を考慮した経路平滑化法,長坂直樹,津坂祐司,第17回ロボティクスシンポジア予稿集,pp.169−174,2012
非特許文献1では、経路を複数の離散した経路点で記述し、各経路点を修正してより好ましい経路データを計算する。この際、移動体と障害物との間の離反距離が大きくなるように経路を修正するため、障害物の近傍の経路点では、その経路点が障害物から離れるように修正される。このため、障害物の近傍では経路点が疎に分布し、それ以外の範囲では経路点が密に分布する経路データに修正されることがある。その結果、移動体と障害物との干渉を判定する精度が低下するという問題が生じる。
本明細書では、経路を複数の離散した経路点で記述したときに、経路点間の距離が適切な距離とされた経路データに修正することができる技術を開示する。
本明細書に開示する第1の経路修正装置は、経路データを入力して記憶する装置と、入力された経路データから、経路データで記述される経路を評価したトータル評価値を計算するトータル評価関数が極値をとるときの経路データを探索する探索装置を備えている。経路データは、経路を複数の離散した経路点で記述するデータを含んでいる。探索装置は、隣接する経路点の経路点間隔に関する経路点間隔評価値を計算する処理と、移動体と障害物の離反距離に関する離反距離評価値を計算する処理と、少なくとも経路点間隔評価値と離反距離評価値とからトータル評価値を計算する処理と、を実行する。探索装置は、経路点間隔が予め設定された設定範囲内から外れると、無限となる経路点間隔評価値を算出する。
この経路修正装置では、隣接する経路点間の距離を評価する経路点間隔評価値を計算する。そして、経路点間隔が予め設定された設定範囲内から外れると、無限となる経路点間隔評価値を算出する。このため、経路点間隔が予め設定された設定範囲内となるように、経路データが修正される。その結果、経路点間の距離を適切な距離としながら、初期経路をより好適な経路データに修正することができる。
本明細書が開示する第2の経路修正装置は、経路データを入力して記憶する装置と、入力された経路データから、経路データで記述される経路を評価したトータル評価値を計算するトータル評価関数が極値をとるときの経路データを探索する探索装置を備えている。経路データは、経路を複数の離散した経路点で記述するデータを含んでいる。探索装置は、隣接する経路点の経路点間隔に関する経路点間隔評価値を計算する処理と、移動体と障害物の離反距離に関する離反距離評価値を計算する処理と、少なくとも経路点間隔評価値と離反距離評価値とからトータル評価値を計算する処理を実行する。探索装置は、さらに計算途上の経路に対して、移動体と障害物との距離を算出し、その算出した距離が予め設定された設定長さを越えるときは、経路点間隔を短くするように評価値を算出する第1経路点間隔評価関数により経路点間隔評価値を算出し、その算出した距離が予め設定された設定長さ以下となるときは、経路点間隔が予め設定された設定範囲内から外れると、無限となる評価値を算出する第2経路点間隔評価関数により経路点間隔評価値を算出する。
この経路修正装置でも、移動体と障害物との距離が近いときには、経路点間隔が予め設定された設定範囲内となるように経路データが修正される。したがって、移動体と障害物との距離が近いときに、経路点間の距離が短くなり過ぎることが防止できる。このため、経路点間の距離を適切な距離としながら、初期経路をより好適な経路データに修正することができる。
実施例の経路修正装置のシステム構成を示す図である。 経路記述パラメータ、経路方向、移動体方向、障害物、経路点間隔に関する評価関数等の意味内容を示す図である。 経路方向と移動体方向の偏差に関する評価関数の意味内容を示す図である。 障害物までの離反距離に関する意味内容を示す図である。 移動体と障害物との離反距離に関する離反距離評価関数を示す図である。 移動体と障害物とが閾値以上離れているときの第1経路間隔評価関数を示す図である。 移動体と障害物とが閾値より接近しているときの第2経路間隔評価関数を示す図である。 経路修正装置によって初期経路を修正するときの処理手順を示すフローチャートである。 初期経路を例示する図である。 修正後の経路を例示する図である。 修正後の経路における経路点の間隔をグラフ化した図である。
以下に説明する実施例の主要な特徴を列記しておく。なお、以下に記載する技術要素は、それぞれ独立した技術要素であって、単独であるいは各種の組合せによって技術的有用性を発揮するものであり、出願時請求項記載の組合せに限定されるものではない。
(特徴1) 上記の経路修正装置では、経路データは、移動体の方向を記述するデータを含んでいてもよい。この場合に、探索装置は、さらに、経路方向と移動体方向の間の偏差に関する偏差評価値を計算する処理を実行してもよい。そして、トータル評価値を計算する処理では、経路点間隔評価値と離反距離評価値と偏差評価値とからトータル評価値を計算してもよい。このような構成によると、移動体の方向(向き)まで考慮して、より好適な経路に修正することができる。
経路データの修正装置の実施例を説明する。最初に用語の意味を明らかにしておく。
経路データとは、経路点(移動体の基準点が通過する点)のx座標と、経路点のy座標と、経路点における移動体方向(x軸の正方向から移動体に固定されている基準線が指向している方向までの反時計回転方向の角θ)が、経路点の順序に従って並んでいるデータをいう。経路データは、式(1)で示され、Piは式(2)で示される。xi,yi,θiを、本明細書では、経路記述パラメータという。i=1は出発点であり、i=Nは目標点である。x1、y1,θ1,xN、yN,θNは与えられている。
Figure 2014211760
Figure 2014211760
初期経路とは、修正装置に入力される、修正処理前の経路データで記述される経路をいう。初期経路データは、RRTのような確率的手法で求めることもできるし、A探索のようなグリッド探索法で求めることができる。本実施例の修正装置は、初期経路データの生成プロセスに限定されないで適用することができる。初期経路は、粗い精度で探索された経路であり、移動体がその初期経路に沿って移動すると、移動体が障害物に接触することがある。あるいは、経路方向が急激に変化することがある。
ここで、移動体は有限の大きさを備えている。障害物は有限の大きさを備えている。凸包とは、物体を包含する凸多角形であり、最短距離計算法とは、離れて存在する2個の凸包のそれぞれにおける最近接点を求め、その距離(凸包間の最短距離)を計算する手法である。最短距離計算法には、GJKアルゴリズムなどの既知の数学手法を利用することができる。
また、ポテンシャル増大開始距離とは、離反距離がそれ以下となると、離反距離評価値を計算するポテンシャル関数がゼロでなくなる距離をいう。許容距離とは、移動体と障害物の間に確保する最小離反距離をいう。許容距離は、ポテンシャル増大開始距離に等しくしてもよいし、異なっていてもよい。また、許容偏差とは、先行経路点から後続経路点を指す方向と、x軸の正方向がなす角をδとしたときに、先行経路点における移動体方向とδの偏差、ならびに、後続経路点における移動体方向とδの偏差に関して許容される最大偏差をいう。
(修正装置の全体構成)
図1は、入力装置2から初期経路を入力し、その初期経路を好ましい経路に修正し、修正後経路を出力装置22から出力する経路修正装置のシステム構成を示している。入力装置2から入力した初期経路は、修正途上の経路記述パラメータ記憶装置4に記憶される。後記するように、経路記述パラメータ記憶装置4の記憶内容は、計算の進行に伴って更新されていく。経路記述パラメータ記憶装置4は、修正途上の経路記述パラメータを記憶している装置ということができる。
経路修正装置は、移動体と障害物との最短距離(離反距離)を判定する装置5と、経路点間隔に関する評価関数を複数記憶している装置6と、移動体と障害物との最短距離判定に基づいて経路点間隔評価関数記憶装置6から経路点間隔に関する評価関数を選択する装置7と、移動体から障害物までの離反距離に関する評価関数を記憶している装置8と、経路方向と移動体方向の偏差に関する評価関数を記憶している装置10と、経路点間隔評価値と偏差評価値と離反距離評価値を合計したトータル評価値を計算するトータル評価関数記憶装置12を備えている。経路記述パラメータ記憶装置4に記憶されている経路記述パラメータをトータル評価関数に代入することで、トータル評価値を計算することができる。図1は、システム構成を理解しやすくするために、トータル評価関数記憶装置12とは別に、経路点間隔評価関数記憶装置6と経路点間隔評価関数選択装置7と離反距離評価関数記憶装置8と偏差評価関数記憶装置10が図示されているが、これらの装置6,7,8,10は、トータル評価関数記憶装置12に内蔵されていてもよい。
経路修正装置は、トータル評価関数を経路記述パラメータで偏微分した偏微分関数記憶装置14を備えている。すなわち、トータル評価関数をxiで偏微分した偏微分関数と、トータル評価関数をyiで偏微分した偏微分関数と、トータル評価関数をθiで偏微分した偏微分関数を記憶している。ここで、i=2,3,・・N―1である。i=1は出発点であり、i=Nは目標点であるから、x1、y1,θ1,xN、yN,θNは固定されている。偏微分関数は、中間の経路点に関して偏微分する。経路記述パラメータ記憶装置4に記憶されている経路記述パラメータを、トータル評価関数をxiで偏微分した偏微分関数と、トータル評価関数をyiで偏微分した偏微分関数と、トータル評価関数をθiで偏微分した偏微分関数に代入することで、3種類の偏微分係数を計算することができる。図1では、微分関数が一つのブロック14に図示されているが、xiの偏微分関数と、yiの偏微分関数と、θiの偏微分関数に分けて図示してもよい。
経路記述パラメータを変数とするトータル評価関数が与えられ、経路記述パラメータを変数とする偏微分関数(トータル評価関数を経路記述パラメータで偏微分した関数)が与えられれば、既知の数学的手法で、トータル評価関数が極値をとる際の経路記述パラメータの値を決定することができる。本実施例では、好ましい経路ほど小さな評価値をとるトータル評価関数を採用する。経路修正装置は、トータル評価関数が最小値をとる経路記述パラメータの値を決定することによって、初期経路を最も好ましい経路に修正する。本実施例では、好ましい経路に修正する速度と安定性にすぐれている準ニュートン法を採用する。ただし、準ニュートン法に限られず、最急降下法、ニュートン法、共役勾配法などの勾配法を用いることもできる。
多変数関数が最小値をとる条件を探査する勾配法では、偏微分関数を利用して、関数値を最小化させる変数の修正方向を決定する。多変数関数の場合には、各変数の修正量の比率を計算する。ここでは、それを探索方向ベクトルという。経路修正装置は、探索方向ベクトルを計算する装置16を備えている。
多変数関数が最小値をとる条件を探査する勾配法では、修正量の大きさが問題となる。修正量が過大であれば、修正方向は正しくても、最小化点を超えた範囲で関数値(トータル評価値)を計算することになり、トータル評価値が増大するという計算結果になってしまう。修正量が過小であれば、修正によって得られる関数値の減少量が微小となり、探索速度が遅くなる。本実施例の経路修正装置では、複数種類のステップ幅(修正量を決めるパラメータ)を順に用意しておき、関数値が減少する制約下で最大の修正量を採用するステップ幅設定装置18を備えている。
経路修正装置では、探索が進行して最小化位置に接近すると、探索方向ベクトル計算装置16で計算された方向に、ステップ幅設定装置18で設定された修正量だけ修正しても、評価値が変化しなくなる。最小化点では、偏微分係数がゼロであることに対応する。本実施例では、探索方向ベクトル計算装置16で計算された方向にステップ幅設定装置18で設定された修正量だけ修正することで関数値が変化(減少)するうちは、修正後の経路記述パラメータによって、経路記述パラメータ記憶装置4の記憶内容を更新する。更新した経路記述パラメータを用いて最小化条件探索処理を繰り返す。経路記述パラメータ記憶装置4の記憶内容は、計算途上の経路記述パラメータに更新されていく。探索方向ベクトル計算装置16で計算された方向にステップ幅設定装置18で設定された修正量だけ修正しても評価値が変化しなくなると、最小化点が検索されたとして処理を終了する。その段階の経路データ(経路記述パラメータ)が、出力装置22から出力される。
(経路データ)
図2において、参照番号26は移動体を示し、出力装置22から出力された修正後経路に沿って移動する。参照番号24は移動体26を包含する凸包であり、図1の修正装置では、凸包24と障害物間の最短距離を求めて離反距離を計算する。参照番号28は移動体26に固定されている基準点であり、参照番号30は移動体26に固定されている基準線である。
経路は、データPi=(xi,yi,θi)(ここでi=1,2,・・N)で定義される。移動体26は、基準点28が(xi,yi)の座標点を辿るように移動する。また、移動体26は、基準線30がθiの方向を向く車体姿勢をとる。本明細書では、角度をx軸の正方向から反時計回転方向に計測した角度で示す。参照記号Ojは、移動体26の走行を妨げる障害物である。本明細書では、出発点(x1,y1)から目標点(xN,yN)までの経路を決定する際に考慮する必要がある障害物がM個存在するものとする。
(初期経路)
図1の経路修正装置に入力される初期経路は、RRTのような確率的手法あるいはA探索のようなグリッド探索法で求めた大局的経路である。データPi=(xi,yi,θi)のN個の並びで記述されている。A探索のようなグリッド探索法で経路を探索すると、移動体方向θiは計算されない。この場合、経路方向と移動体方向が等しいとすることができる。あるいは別の手法で、移動体方向を決定してもよい。
初期経路のうち、出発点のデータP1=(x1,y1,θ1)と目的点のデータPN=(xN,yN,θN)は、確定的なものであり、図1の修正装置で修正することはない。ただし、経路点のデータPi(i=2,・・,N−1)は最善なものでない可能性があり、図1の修正装置で修正する。初期経路には下記の問題が内在されていることがある。
(1)移動体と障害物が干渉する(重なりあう)経路である可能性がある。
(2)移動体と障害物が重なり合わないまでも非常に接近する経路である可能性がある。
(3)経路方向が急激に変化する経路である可能性がある。
(4)経路方向と移動体方向が不一致であって、移動体が追従できない経路である可能性がある。図1の修正装置は、上記の問題を解消する側に修正した経路を計算する。その修正計算のために、修正装置は下記の3種類の評価関数を利用する。
(経路点間隔評価関数)
出発点から目標点に至る経路の全長(各経路点の間隔(先行経路点と後続経路点の距離)の総和)は短いことが好ましい。ただし、出発点から目標点に至る経路に、経路点の間隔が長い領域が含まれると、その領域において障害物との干渉判定の精度が低下する。その結果、障害物との干渉(衝突)が生じる経路が最終的な経路として出力される可能性が生じる。したがって、障害物との干渉が問題となる領域においては、経路点の間隔が長くなり過ぎることを抑制する必要がある。本実施例では、移動体と障害物との最短距離(離反距離)に応じて、2種類の経路点間隔評価関数を利用する。すなわち、修正途上の経路において、移動体と障害物との離反距離が設定値(正の値)以上であると判定される場合は、経路点間隔が短い経路に修正される第1経路点間隔評価関数を利用する。一方、修正途上の経路において、移動体と障害物との離反距離が設定値未満であると判定される場合は、経路点間隔が予め設定された範囲内となるように修正される第2経路点間隔評価関数を利用する。なお、移動体と障害物との距離判定については後述する。
第1経路点間隔評価関数を図6に示す。係数Kaは、トータル評価値に占める経路点間隔評価値の重みを示し、距離liは、座標点(xi,yi)から座標点(xi+1,yi+1)までの距離である。第1経路点間隔評価関数は全長(座標点(xi,yi)から座標点(xi+1,yi+1)までの距離)が短いほど小さな評価値が計算される。したがって、第1経路点評価関数は、経路点間隔を短くする方向に経路を修正する。
第2経路点間隔評価関数を図7に示す。係数Kaと距離liとは、第1経路点間隔評価関数と同一である。係数Kは、第2経路点間隔評価関数の第2項の重みを示す係数であり、係数Kは、第2経路点間隔評価関数の第3項の重みを示す係数である。距離l,lhigh(>l)は、予め設定された距離である。図7から明らかなように、第2経路点間隔評価関数は、3つの項(評価関数)の和となる。具体的には、第1項(1/2・Ka(li−l)は、経路点間隔(距離li)が距離lとなるときに最小値(ゼロ)となる評価値を算出し、経路点間隔(距離li)が距離lより離れるほど大きな評価値を算出する。第2項(K・log0.5li)は、経路点間隔(距離li)が距離ゼロとなると無限大となる評価値を算出し、経路点間隔(距離li)が大きくなるほど小さな評価値を算出する。第3項(K・log0.5{−(li−lhigh)}は、経路点間隔(距離li)が距離lhighとなると無限大となる評価値を算出し、経路点間隔(距離li)が小さくなるほど小さな評価値を算出する。したがって、第2経路点間隔評価関数は、経路点間隔(距離li(座標点(xi,yi)から座標点(xi+1,yi+1)までの距離))が距離lとなるときに最小値となり、経路点間隔がゼロ又は距離lhighに近づくほど大きくな値となり、経路点間隔がゼロ又は距離lhighのときに無限大となる。このため、第2経路点間隔評価関数は、経路点間隔が予め設定された範囲(0<l<lhigh)となるように経路を修正する。より詳細には、経路点間隔が予め設定された距離lとなるように経路を修正する。
なお、距離lは、移動体と障害物との干渉判定の精度が低下しないような距離に設定される。例えば、距離lは、初期経路長を(N−1)で割った値とすることができる。この場合、第2経路点間隔評価関数は、初期経路長に近づくように経路を修正することとなる。また、第1経路点間隔評価関数と第2経路点間隔評価関数とを切換える設定値は、適宜設定することができる。例えば、設定値としては、上記の距離lを用いてもよい。この場合、移動体と障害物との最短距離(離反距離)が距離l以上となると、経路点間隔が短くなるように経路を修正し、一方、移動体と障害物との最短距離(離反距離)が距離l未満となると、経路点間隔が距離lとなるように経路を修正する。これによって、移動体と障害物との最短距離(離反距離)が距離l未満となる領域においては、経路点間隔が距離lとなるように経路が修正される。
経路修正装置は、i=1からN−1までの各距離liについて上記の経路点間隔評価関数を用いて評価値を算出し、それら評価値の総和を取れば、経路全長に関する評価値が求められる。具体的には、各距離liについて移動体と障害物との最短距離が設定値(例えば、距離l)以上となるか否かを判定し、移動体と障害物との最短距離が設定値以上となる場合は第1経路点間隔評価関数によって評価値を算出し、移動体と障害物との最短距離が設定値未満となる場合は第2経路点間隔評価関数によって評価値を算出する。そして、各距離liについて算出した評価値の総和を取ることで、経路全長に関する評価値が得られる。これによって、経路点間隔が適切な範囲としながら、全長が短い経路に修正される。もっとも他の評価関数をも利用するために、経路点間隔のみを考慮して修正されるものでないが、経路点間隔をも考慮した修正計算が実行されるということはできる。
(最短距離(離反距離))
ここで、移動体と障害物との最短距離(離反距離)について図4を参照して説明する。図4は、移動体が経路点Piに移動した状態を示し、近傍にj番目の障害物が位置している状態を示している。移動体26の凸包24の存在範囲は、座標値(xi,yi)と移動体方向θiから計算することができる。障害物は、動かない凸包32の存在範囲を記述するデータで記述されている。凸包24と凸包32の位置と姿勢が与えられれば、GJKアルゴリズムなどの既知の数学手法によって、凸包24の輪郭上にあって凸包32に最も近い点Ri,jと、凸包32の輪郭上にあって凸包24に最も近い点Oi,jを計算し、その距離lei,jを計算することができる。図示の場合は、最近接点が頂点である場合を例示しているが、辺上に最近接点がある場合にも対応することができる。図4に示す式には、近接点Oi,jの座標値と、近接点Ri,jの座標値と、最短距離(離反距離)lei,jの値を示している。したがって、最短距離(離反距離)lei,jが設定値(例えば、距離l)以上となるか否かによって、選択する経路点間隔評価関数が変化することとなる。
なお、移動体と障害物が接近して重なりあう場合(すなわち、干渉する場合)、GJKアルゴリズムでは正の値を持つlei,jが計算される。この場合のlei,jは、移動体と障害物が重なりあっている距離であり、両者を離反させるのに必要な移動距離に等しい。本実施例では、移動体と障害物が接近して重なりあう場合には、lei,jが負の値を持つ距離だとして扱う。したがって、最短距離(離反距離)lei,jが0以下となると、移動体と障害物とが干渉していると判断することができる。したがって、最短距離(離反距離)lei,jが0以下の場合、移動体と障害物との最短距離が設定値(例えば、距離l)以下となり、第2経路間隔評価関数が選択される。その結果、経路が長くなるように修正され、移動体と障害物との干渉が解消される。
なお、移動体と障害物との距離判定は、経路上の各経路点Pi(i=2,・・,N−1)において行われる。また、経路上の各経路点Pi(i=2,・・,N−1)において、M個の障害物のそれぞれについて最短距離(離反距離)lei,jが算出される。
(離反距離評価関数)
上述した最短距離(離反距離)lei,jは、移動体と障害物とが干渉する確率を評価するための離反距離評価関数としても利用される。すなわち、離反距離lei,jが大きいほど移動体と障害物とが干渉する確率は低くなり、離反距離lei,jが小さいほど移動体と障害物とが干渉する確率は高くなる。したがって、離反距離を評価する離反距離評価関数によって、移動体と障害物とが干渉する確率を評価することができる。
図5は、離反距離評価関数Uei,jを示している。係数Keは、離反距離評価関数の重み係数である。離反距離lei,jがポテンシャル増大開始距離le0以上であれば、離反距離評価値はゼロである。離反距離lei,jがポテンシャル増大開始距離le0以下となると、離反距離評価値は増加し始める。移動体と障害物が接近して離反距離lei,jがゼロに近づくほど、離反距離評価値は大きくなる。ただし、移動体と障害物が接近して離反距離lei,jがゼロとなっても、離反距離評価値は無限大に発散せず、有限の値が計算される関数形を採用している。
なお、実際には、移動体と障害物が重なりあう現象は生じない。しかしながら、経路を修正する計算の途上では、移動体と障害物とが重なりあう経路に関する計算が必要とされることがある。例えば、移動体と障害物が重なりあう初期経路を、重なり合わない経路に修正する場合には、移動体と障害物が重なりあう経路に関する計算が必要とされる。本実施例では、移動体と障害物が重なりあう場合も有限の離反距離評価値が計算される離反距離評価関数を用いるので、修正計算の途上で計算不能に陥ることがない。また、障害物はM個存在するため、移動体が経路点Piにあるときの離反距離評価値は、各障害物に対する離反距離評価値の和として求められる。
(偏差評価関数)
本実施例では、移動体26は、移動方向(経路方向)と独立して車体方向(姿勢方向、移動体方向)を変える機構を備えていない。このため、経路点Pi=(xi,yi)を辿って移動する際の移動体方向θiが経路方向δiに一致することは保証されていない。その問題を解消するために、経路修正装置には経路方向δiと移動体方向θiを一致させるように修正する評価関数が導入されており、移動体がPi(xi,yi)を辿って移動すると、移動体方向θiが経路方向δiによく一致する関係が得られるようにしている。移動方向から独立して車体方向を変える機構を備えていない移動体26の経路を修正する際に、経路方向δiと独立して移動体方向θiを定義する経路記述データを採用しても、問題が生じないようにしている。
図3の(b)は、経路点Piから経路点Pi+1に至る経路方向(δi)と、移動体方向の偏差を示している。実際には、経路点Piにおける移動体方向θiと、経路点Pi+1における移動体方向θi+1は相違している。そこで、2種類の偏差評価関数を導入している。
第1偏差評価関数は、図3の(3)に示すように、経路方向δiと移動体方向θiの偏差φiを自乗した値を計算する。係数Kbは、トータル評価値に占める第1偏差評価値の重みを示している。図3の(4)に示すように、i=1からN−1までの総和を取れば、第1偏差評価値を経路全長に亘って累積した第1偏差評価値Ubを計算することができる。偏差φiが小さいほど評価値が小さくなる。評価値を最小化する方向に修正すると、偏差φiが小さくなる経路に修正される。
第2偏差評価関数は、図3の(5)に示すように、経路方向δiと移動体方向θi+1の偏差ηiを自乗した値を計算する。係数Kcは、トータル評価値に占める第2偏差評価値の重みを示している。図3の(6)に示すように、i=1からN−1までの総和を取れば、第2偏差評価値を経路全長に亘って累積した第2偏差評価値Ucを計算することができる。偏差ηiが小さいほど評価値が小さくなる。評価値を最小化する方向に修正すると、偏差ηiが小さくなる経路に修正される。
移動体方向θを考慮すると、障害物までの離反距離を増大させられることがある。移動体方向θを加味して経路を修正すると、移動体方向θを振りながら移動することで移動体が障害物の間を縫って移動する経路に修正することができる。この場合には、移動体方向θによく一致する経路方向δとなる経路点位置に修正されることになる。
図3の(a)と(b)とを比較すると明らかなように、経路点Piで第1評価関数と第2評価関数を計算することになる。その際には、共通して移動体方向θiを利用する。第1評価関数は、移動体方向θiと経路方向δiの偏差φiが小さくなる方向に修正する。第2評価関数は、移動体方向θiと経路方向δi-1の偏差ηi-1が小さくなる方向に修正する。共通の移動体方向θiに対する偏差が小さくなる側に修正すると、経路方向δi-1と経路方向δiの差も小さくなっていく。経路方向δi-1と経路方向δiの差が大きいことは、経路が急激に方向を変えることを意味する。経路方向δi-1と経路方向δiの差を小さくすることは、緩やかに方向を変える経路に修正することに相当する。1つの経路点Piに対して、その経路点Piにおける移動体方向θiと直前の経路方向δi-1の偏差ηi-1に対応する第2偏差評価関数と、その経路点Piにおける移動体方向θiと直後の経路方向δiの偏差φiに対応する第1偏差評価関数を用いて評価値を計算するようにすると、急激に方向を変える経路から緩やかに方向を変える経路に修正される。経路の平滑化処理が実施される。
(トータル評価関数)
本実施例では、下記式で示されるトータル評価関数を利用する。
トータル評価関数
=経路点間隔評価関数+偏差評価関数+離反距離評価関数
=経路点間隔評価関数+第1偏差評価関数+第2偏差評価関数+離反距離評価関数
すなわち、下記の式3を用いて、経路の修正処理を実行する。
Figure 2014211760
本実施例では、トータル評価関数を経路記述パラメータxi,yi,θiの各々で偏微分した偏微分関数を利用して、トータル評価値を最小化する経路記述パラメータxi,yi,θiを探索する。トータル評価関数=経路点間隔評価関数+偏差評価関数+離反距離評価関数であるため、トータル評価関数の偏微分関数=経路点間隔評価関数の偏微分関数+偏差評価関数の偏微分関数+離反距離評価関数の偏微分関数の関係が成立する。
(経路点間隔評価関数の偏微分関数)
経路点間隔評価関数Uaの偏微分関数を求めると下記となる。上述したように、経路点間隔評価関数は、移動体と障害物との最短距離が設定値以上となるか否かで2種類の経路点間隔評価関数を用いる。具体的には、移動体と障害物との最短距離が設定値以上となるときの第1経路点間隔評価関数の偏微分関数は式4となり、移動体と障害物との最短距離が設定値未満となるときの第2経路点間隔評価関数の偏微分関数は式5となる。
Figure 2014211760
Figure 2014211760
(離反距離評価関数の偏微分関数)
離反距離評価関数Ueの偏微分関数を求めると下記となる。
Figure 2014211760
(偏差評価関数の偏微分関数)
図3の(4)(6)に示す第1偏差評価関数Ubと第2偏差評価関数Ucの偏微分関数から偏差評価関数Ub+Ucの偏微分関数を求めると下記となる。
Figure 2014211760
(評価関数の最小化過程)
関数とその関数の偏微分関数が判明していれば、既知の数学的手法で、関数を最小化する変数の値を決定することができる。その問題は、最急降下法、ニュートン法、準ニュートン法、共役勾配法などの勾配法で解くことができる。本実施例では、計算の収束速度と安定性のバランスが優れている準ニュートン法を用いる。
準ニュートン法などの勾配法では、計算途上の変数値(複数個の変数を扱うからベクトルとなる)における関数値から、関数値を減少させる変数の修正方向を特定する。偏微分関数から変数の変化に起因して関数値が変化する量を計算することができるため、関数値を減少させる変数の修正方向を特定することができる。次に、修正量を仮定し、変数群を修正方向へ修正量だけ修正した場合の関数値を計算する。関数値を減少させる修正方向を用いることから、関数値は減少するはずである。しかしながら、修正量が大きすぎれば、関数が最小値を取る位置を越えて修正することになり、関数値が増大してしまう。そこで、関数値が減少するという制約を満たす修正量の範囲内で修正する。勾配法では、
(1)現在の変数値の並び(ベクトル)よりも関数値が減少する修正方向を計算し、
(2)関数値が減少するという制約を満たす修正量で変数値の並びを修正方向に修正し、
(2)で修正された変数値の並びを(1)の変数値の並びとして(1)と(2)を繰り返す。
この結果、変数値の並びは、関数値を最小化するものに近づいていく。
関数値を最小化する変数値に近づくと、偏微分係数がゼロなって、変数値を修正方向に修正量だけ修正しても関数値は変化しなくなる。このことから、関数値を最小化する変数の並びに修正されたことがわかる。
準ニュートン法では、修正方向dkが下記の式で計算される。なお添え字kは、前記(1)(2)の修正処理をすでにk-1回繰り返しており、これからk回目の修正処理を実施することを示している。
Figure 2014211760
上記式の各項目は下記に示されるものである。
Figure 2014211760
Figure 2014211760
本実施例では、修正量を下記の表に従って変えていく。
Figure 2014211760
本実施では、式8で決定された修正方向に、式8で決定された修正量に「1.0」を乗じた修正量だけ修正し、関数値を計算する。関数値が増加すれば、「1.0」の乗数では過大であるとし、乗数を「0.1」に変更して再度関数値を計算する。それでも関数値が増加すれば、乗数を「0.01」に変更する。それでも関数値が増加すれば、乗数を「0.001」に変更する。
最小値に近づくにつれ、修正前のトータル評価関数の値と修正後のトータル評価関数の値の差の絶対値は減少していく。本実施例では、その差の絶対値が、閾値(例えば、10−6)を下回ったときに、トータル評価値を最小化する経路記述パラメータの並びに修正されたとし、修正処理を終了する。そうして修正された経路記述パラメータの並びを出力装置22から出力する。
なお、計算された移動体方向θと実際の移動体方向が一致するようにするためには、経路方向δと移動体方向θの偏差に関して許容偏差角を設定しておいてもよい。下記式を加えて修正処理を進めれば、経路方向δと移動体方向θの偏差が許容偏差角η0,φ0以下であることが保証されている経路記述データに修正される。経路方向δと移動体方向θの偏差が小さい経路が計算されていれば、計算された移動体方向θと実際の移動体方向はよく一致する。
Figure 2014211760
移動体と障害物の離反距離に関して最短値を設定して修正作業を進めることもできる。下記式を加えて修正処理を進めれば、移動体と障害物の最短離反距離が許容距離minlen以上であることが保証されている経路記述データに修正される。許容距離minlenは、上述したポテンシャル増大開始距離le0に等しくてもよいし、異なっていてもよい。
Figure 2014211760
また、式11,12の両方を用いて経路修正処理を実施してもよい。
次に、上述した経路修正装置で初期経路を最適な経路に修正する際に実行される処理について、図8のフローチャートに基づいて説明する。経路修正装置は、まず、初期経路(修正途上においては修正途上の経路)を入力する(S10)。初期経路又は修正途上経路が入力されると、次に、経路点Piに移動体を配置したときに、移動体と障害物との最短距離が設定値以上となるか否かを判定する(S12)。具体的には、最短距離(離反距離)lei,jが設定値(距離l)以上となるか否かで判定する。移動体と障害物との最短距離が設定値以上となる場合(S12でYES)は、第1経路点間隔評価関数を用いて評価値を算出する(S16)。一方、移動体と障害物との最短距離が設定値未満の場合(S12でNO)は、第2経路点間隔評価関数を用いて評価値を算出する(S14)。そして、移動体と障害物の離反距離に関する評価値を離反距離評価関数を用いて算出し(S17)、経路方向と移動体の方向(車両方向)の偏差を評価するための評価値を偏差評価関数を用いて算出する(S18)。
S14,16,18によって各評価値が算出されると、それらの評価値からトータル評価値を算出する(S20)。次に、経路修正装置は、トータル評価値を最小化するために、経路点間隔評価関数と偏差評価関数と離反距離評価関数を偏微分する(S22)。次に、上述した方法で探索方向ベクトルを計算し(S24)、乗数(ステップ幅)を設定し(S26)、計算した探索方向ベクトルと設定した乗数に基づいてトータル評価関数を最小化する(S28)。次に、最小化されたトータル評価関数で算出されるトータル評価値の変化量(前回の処理で算出されたトータル評価値との差分)が閾値以下となるか否かを判定する(S30)。トータル評価値の変化量が閾値以下でない場合(S30でNO)は、S10に戻って、S10からの処理を繰り返す。これによって、修正された経路に基づいて、S10からの処理が繰り返される。一方、トータル評価値の変化量が閾値以下となる場合(S30でYES)は、初期経路が最適な経路に修正されたとして、その修正後の経路を出力して、処理を終了する(S32)。
上述した経路修正装置によって初期経路を修正した例を説明する。図9は、初期経路の一例を示している。図9に示す初期経路では、移動体は2つの障害物の間を通過し、その後、直角に進行方向を変えている。図10は、経路修正装置から出力された経路を示している。出力された経路では、移動体は、障害物を通過した後、穏やかに進行方向を変え、目標点に目標姿勢で到達している。図11は、修正後の経路(図10)における、隣接する経路点間の距離(経路点Pと経路点Pi+1の距離)を示すグラフである。グラフの縦軸は経路点間の距離(経路点間隔)であり、横軸が経路点間隔の番号i(すなわち、経路点Pと経路点Pi+1の距離)である。なお、第1経路点間隔評価関数のみを用いて経路を修正したときの比較例を併せて示す。また、計算条件は、l=0.05,lhigh=0.15とした。図11から明らかなように、経路点間隔は、障害物近傍では設定値l(0.05)よりずれ、障害物から離れると設定値lとなった。ただし、本実施例では、最小の経路点間隔が0.035mとなり、最大経路点間隔が0.068mとなったのに対して、比較例では、最小の経路点間隔が0.0003mとなり、最大経路点間隔が0.089mとなった。本実施例によって、経路点間隔が小さくなり過ぎることを防止できると共に、経路点間隔が大きくなり過ぎることを防止することができる。
上述したように、本実施例の経路修正装置では、移動体と障害物との最短距離が設定値以上となるときは第1経路点間隔評価関数を用いて評価値を算出し、移動体と障害物との最短距離が設定値未満となるときは第2経路点間隔評価関数を用いて評価値を算出する。このため、経路点間隔が大きくなり過ぎることが防止され、移動体と障害物の干渉の有無を判定する処理の精度が低下することを抑制することができる。また、経路点間隔が小さくなり過ぎることが防止され、曲率半径の小さな急カーブに経路が修正されることを防止することができる。
以上、本発明の具体例を詳細に説明したが、これらは例示にすぎず、特許請求の範囲を限定するものではない。特許請求の範囲に記載の技術には、以上に例示した具体例を様々に変形、変更したものが含まれる。
例えば、上述した実施例では、移動体と障害物との最短距離に応じて2種類の経路点間隔評価関数を用いたが、移動体と障害物との最短距離に関係なく、第2経路点間隔評価関数のみを用いて評価値を算出してもよい。すなわち、移動体と障害物との最短距離が設定値以上となる場合にも、第2経路点間隔評価関数を用いて評価値を算出してもよい。このような構成としても、経路点間隔を適切な範囲とすることができる。
また、本明細書または図面に説明した技術要素は、単独であるいは各種の組合せによって技術的有用性を発揮するものであり、出願時請求項記載の組合せに限定されるものではない。また、本明細書または図面に例示した技術は複数目的を同時に達成するものであり、そのうちの一つの目的を達成すること自体で技術的有用性を持つものである。
2:初期経路入力装置
4:修正途上の経路記述パラメータの記憶装置
5:移動体と障害物との距離を判定する判定装置
6:経路点間隔に関する評価関数の記憶装置
8:障害物までの離反距離に関する評価関数の記憶装置
10:経路方向と移動体方向の偏差に関する評価関数の記憶装置
12:トータル評価関数の記憶装置
14:経路記述パラメータについてのトータル評価関数の微分関数の記憶装置
16:探索方向ベクトルの計算装置
18:ステップ幅の設定装置
20:トータル評価値の変化量の判定装置
22:修正後経路出力装置

Claims (3)

  1. 経路データを入力して記憶する装置と、
    入力された経路データから、経路データで記述される経路を評価したトータル評価値を計算するトータル評価関数が極値をとるときの経路データを探索する探索装置を備えており、
    経路データは、経路を複数の離散した経路点で記述するデータを含んでおり、
    探索装置は、
    隣接する経路点の経路点間隔に関する経路点間隔評価値を計算する処理と、
    移動体と障害物の離反距離に関する離反距離評価値を計算する処理と、
    少なくとも経路点間隔評価値と離反距離評価値とからトータル評価値を計算する処理と、を実行し、
    探索装置は、経路点間隔が予め設定された設定範囲内から外れると、無限となる経路点間隔評価値を算出することを特徴とする経路データの修正装置。
  2. 経路データを入力して記憶する装置と、
    入力された経路データから、経路データで記述される経路を評価したトータル評価値を計算するトータル評価関数が極値をとるときの経路データを探索する探索装置を備えており、
    経路データは、経路を複数の離散した経路点で記述するデータを含んでおり、
    探索装置は、
    隣接する経路点の経路点間隔に関する経路点間隔評価値を計算する処理と、
    移動体と障害物の離反距離に関する離反距離評価値を計算する処理と、
    少なくとも経路点間隔評価値と離反距離評価値とからトータル評価値を計算する処理と、を実行し、
    探索装置は、さらに計算途上の経路に対して、移動体と障害物との距離を算出し、
    その算出した距離が予め設定された設定長さを越えるときは、経路点間隔を短くするように評価値を算出する第1経路点間隔評価関数により経路点間隔評価値を算出し、
    その算出した距離が予め設定された設定長さ以下となるときは、経路点間隔が予め設定された設定範囲内から外れると、無限となる評価値を算出する第2経路点間隔評価関数により経路点間隔評価値を算出することを特徴とする経路データの修正装置。
  3. 経路データは、移動体の方向を記述するデータを含んでおり、
    探索装置が、さらに、経路方向と移動体方向の間の偏差に関する偏差評価値を計算する処理を実行し、
    トータル評価値を計算する処理では、経路点間隔評価値と離反距離評価値と偏差評価値とからトータル評価値を計算することを特徴とする請求項1又は2に記載の経路データの修正装置。
JP2013087749A 2013-04-18 2013-04-18 経路修正装置 Active JP5986533B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2013087749A JP5986533B2 (ja) 2013-04-18 2013-04-18 経路修正装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2013087749A JP5986533B2 (ja) 2013-04-18 2013-04-18 経路修正装置

Publications (2)

Publication Number Publication Date
JP2014211760A true JP2014211760A (ja) 2014-11-13
JP5986533B2 JP5986533B2 (ja) 2016-09-06

Family

ID=51931475

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2013087749A Active JP5986533B2 (ja) 2013-04-18 2013-04-18 経路修正装置

Country Status (1)

Country Link
JP (1) JP5986533B2 (ja)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2017151687A (ja) * 2016-02-24 2017-08-31 本田技研工業株式会社 移動体の経路計画生成装置
JP2019108124A (ja) * 2016-03-15 2019-07-04 本田技研工業株式会社 車両制御装置、車両制御方法、および車両制御プログラム
CN112381425A (zh) * 2020-11-18 2021-02-19 天津市城市规划设计研究总院有限公司 适用于农业用地生产潜力评价的用地连片度指标计算方法
CN112712733A (zh) * 2020-12-23 2021-04-27 交通运输部公路科学研究所 基于车路协同的碰撞预警方法、系统和路侧单元

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004325083A (ja) * 2003-04-21 2004-11-18 I Transport Lab Co Ltd 経路決定装置、経路決定方法および経路決定プログラム
JP2006195969A (ja) * 2004-12-14 2006-07-27 Honda Motor Co Ltd 自律移動ロボットの移動経路生成装置
JP2010127646A (ja) * 2008-11-25 2010-06-10 Toshiba Corp 経路を選択する装置、方法およびプログラム
JP2012190405A (ja) * 2011-03-14 2012-10-04 Toyota Motor Corp 経路情報修正装置、軌道計画装置、及びロボット
JP2013019683A (ja) * 2011-07-07 2013-01-31 Hitachi Ltd 経路作成システム、経路作成方法及びプログラム

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004325083A (ja) * 2003-04-21 2004-11-18 I Transport Lab Co Ltd 経路決定装置、経路決定方法および経路決定プログラム
JP2006195969A (ja) * 2004-12-14 2006-07-27 Honda Motor Co Ltd 自律移動ロボットの移動経路生成装置
JP2010127646A (ja) * 2008-11-25 2010-06-10 Toshiba Corp 経路を選択する装置、方法およびプログラム
JP2012190405A (ja) * 2011-03-14 2012-10-04 Toyota Motor Corp 経路情報修正装置、軌道計画装置、及びロボット
JP2013019683A (ja) * 2011-07-07 2013-01-31 Hitachi Ltd 経路作成システム、経路作成方法及びプログラム

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2017151687A (ja) * 2016-02-24 2017-08-31 本田技研工業株式会社 移動体の経路計画生成装置
JP2019108124A (ja) * 2016-03-15 2019-07-04 本田技研工業株式会社 車両制御装置、車両制御方法、および車両制御プログラム
CN112381425A (zh) * 2020-11-18 2021-02-19 天津市城市规划设计研究总院有限公司 适用于农业用地生产潜力评价的用地连片度指标计算方法
CN112712733A (zh) * 2020-12-23 2021-04-27 交通运输部公路科学研究所 基于车路协同的碰撞预警方法、系统和路侧单元

Also Published As

Publication number Publication date
JP5986533B2 (ja) 2016-09-06

Similar Documents

Publication Publication Date Title
JP5986533B2 (ja) 経路修正装置
WO2017173990A1 (zh) 一种机器人避障中的最短路径规划方法
JP5993784B2 (ja) 経路修正装置
Yang Anytime synchronized-biased-greedy rapidly-exploring random tree path planning in two dimensional complex environments
KR101585504B1 (ko) 자율 이동 차량의 경로 생성 방법 및 경로 생성 장치
CN112683278B (zh) 一种基于改进a*算法和贝塞尔曲线的全局路径规划方法
CN109491389A (zh) 一种具有速度约束的机器人轨迹跟踪方法
KR101764653B1 (ko) 이동체의 경로계획 장치 및 그 계획 방법
CN112327876B (zh) 一种基于终距指数的机器人路径规划方法
Song et al. T $^\star $: Time-Optimal Risk-Aware Motion Planning for Curvature-Constrained Vehicles
US20220203534A1 (en) Path planning method and biped robot using the same
EP3316065A1 (en) Method for determining a path
US20210269030A1 (en) Vehicle Control Device, Vehicle Control Method, and Vehicle Control System
CN108803592B (zh) 用于避免机器人碰撞的方法、装置和系统
JP5923973B2 (ja) 経路修正装置
JP2008040702A (ja) 近似点群データからの近似曲線生成プログラム及び方法
Do et al. Narrow passage path planning using fast marching method and support vector machine
JP6131587B2 (ja) 経路修正装置
JP2017146893A (ja) 自己位置推定方法
Lin et al. Humanoid navigation planning in large unstructured environments using traversability-based segmentation
US11579587B2 (en) Automatic program-correction device, automatic program-correction method, and automatic path-generation device
CN110749332B (zh) Rs曲线的曲率优化方法、装置、计算机设备及存储介质
CN114939872B (zh) 基于MIRRT*-Connect算法的智能仓储冗余机械臂动态避障运动规划方法
JP6941039B2 (ja) 軌道計画装置、軌道計画方法、およびプログラム
JP2006155559A (ja) 経路計画装置

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20150611

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20160413

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20160426

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20160607

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20160802

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20160805

R150 Certificate of patent or registration of utility model

Ref document number: 5986533

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250