JP2016141300A - 鉄道車両の走行位置検出システム - Google Patents

鉄道車両の走行位置検出システム Download PDF

Info

Publication number
JP2016141300A
JP2016141300A JP2015019431A JP2015019431A JP2016141300A JP 2016141300 A JP2016141300 A JP 2016141300A JP 2015019431 A JP2015019431 A JP 2015019431A JP 2015019431 A JP2015019431 A JP 2015019431A JP 2016141300 A JP2016141300 A JP 2016141300A
Authority
JP
Japan
Prior art keywords
curve
curve information
candidate
data
information
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
JP2015019431A
Other languages
English (en)
Other versions
JP6509573B2 (ja
Inventor
健介 池本
Kensuke Ikemoto
健介 池本
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.)
Nippon Sharyo Ltd
Original Assignee
Nippon Sharyo Ltd
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 Nippon Sharyo Ltd filed Critical Nippon Sharyo Ltd
Priority to JP2015019431A priority Critical patent/JP6509573B2/ja
Publication of JP2016141300A publication Critical patent/JP2016141300A/ja
Application granted granted Critical
Publication of JP6509573B2 publication Critical patent/JP6509573B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Train Traffic Observation, Control, And Security (AREA)

Abstract

【課題】識別子情報を読み込んで、特定路線のうちの、鉄道車両が走行している特定路線区間を決定することにより、データ処理時間を短縮できる鉄道車両の走行位置検出システムを提供すること。【解決手段】鉄道車両は、線路近くに設置された識別子の識別子情報を読み込んで、特定路線のうちの、鉄道車両が走行している特定路線区間を決定する特定路線決定手段と、複数の路線網を鉄道車両が走行したときの実測走行データを記憶する走行データ記憶手段と、特定路線決定手段が決定した特定路線区間を含む特定路線内の実測走行データを、走行データ記憶手段から読み出し、現在走行中に検知している現在走行データと比較することにより、自車の走行位置を検出すること、を特徴とする。【選択図】図22

Description

本発明は、鉄道車両が走行中に自車の走行位置を検出する鉄道車両の走行位置検出システムに関するものである。
近年、鉄道車両が線路における曲線部を通過する場合、高速に且つ乗り心地良く通過するために、車体を左又は右に傾ける車体傾斜制御が実用化されている。この車体傾斜制御は、曲線部に進入する前に、的確な地点で曲率等の曲線情報に基づいて開始しなければならない。このため、鉄道車両は、走行中に線路における曲線部の正確な曲線情報を保持すると共に、自車の走行位置を検出すること(地点検知)が必要である。
ここで、自車の走行位置を検出する場合、例えば速度発電機を用いて走行した距離を測定している。しかし、走行中に車輪がレールに対して空転又は滑走する場合があるため、測定された実測距離は、正確な走行距離と一致せず、大きな誤差を含むものになる。このため、従来の地点検知では、線路に沿って任意の位置に配置された多数の地上子(識別子)の位置情報を利用して、実測距離を補正するようになっていた。
すなわち、鉄道車両は、地上子を通過する際にその地上子を検知して、地上子と車上子との通信によって検知した地上子の正確な位置情報を取得する。そして、取得した位置情報と実測距離とを比較して誤差を算出して、地上子を検知する度に実測距離を補正する。こうして、従来の地点検知では、自車の走行位置を、走行開始位置から直前に検知した地上子の位置までの距離と、その地上子の位置から測定した実測距離とを加算したものとしていた。
しかしながら、上述したように地上子の位置情報を利用する場合、仮に地上子を検知できなかったときには、鉄道車両は自車の走行位置を正確に検出できなくなる。このため、例えば、本出願人が提案している下記特許文献1では、地上子の検知を利用することなく、曲線部を通過したときに測定した実測曲線情報を利用して、どの曲線部を通過したかを特定し、特定した曲線部からの実測距離に基づいて自車の走行位置を検出するようになっている。
具体的には、検出したヨー角速度から算出した実測曲率と、予め記憶しているデータ曲率とを相関関数を用いて相関演算し、演算された相関値の最大値が設定値以上となった地点を特定して、自車の走行位置を検出するようになっている。模式的に説明すると、図21に示すように、検索区間内において、実測曲率とデータ曲率との相関演算を、演算区間の開始点をずらしながら連続的に複数回行う。そして、演算された相関値の最大値が例えば0.95以上になる地点を特定している。
また、鉄道車両が、線路近くに設置された識別子の識別子情報を読み込んで、特定路線のうちの、鉄道車両が走行している特定路線区間を決定する特定路線決定手段を有するものを本出願人は、台湾に出願し、公開されている(特許文献2)。
特開2007−15483号公報 台湾特許出願公開第201418088号公報
しかしながら、従来の鉄道車両の走行位置検出システムには、次のような問題があった。
(1)直前に通過した曲線部の形状と、鉄道車両が通過し得る全ての曲線部の中で何れかの曲線部の形状とが似ている場合がある。この場合、特許文献1の地点検知方法において、直前に曲線部を通過したときに測定した実測曲率と、通過し得る曲線部毎に予め記憶しているデータ曲率とを相関関数を用いて相関演算すると、設定値以上になる相関値の最大値が複数個存在して、実測曲率と各データ曲率とが一致する地点が非常に特定し難い。
この結果、鉄道車両は直前に通過した曲線部がどの曲線部なのかが明確に分からず、自車の走行位置を検出できない可能性があった。要するに、複数本の線路が並行して敷設された区間を走行する場合、片側通行が主流である状況においては、形状が類似した曲線部が検出され難くて、特定する曲線部を誤認することが少ない。しかしながら、単線双方向運用を行う状況においては、並行する複数の線路上に形状が類似した曲線部が存在し、それらを区別する必要がある。
(2)特に、複数路線網では、曲線部の個数が多くなるため、データの処理時間が長くなり、実用上問題となっていた。
そこで、本発明は、上記した課題を解決すべく、識別子情報を読み込んで、特定路線のうちの、鉄道車両が走行している特定路線区間を決定することにより、データ処理時間を短縮できる鉄道車両の走行位置検出システムを提供することを目的とする。
本発明に係る鉄道車両の走行位置検出システムは、次の構成を有し、次のような作用、効果を奏する。
(1)鉄道車両が、複数路線網から選択された複数の特定路線区間を含む特定路線を走行すること、鉄道車両が自車の走行位置を検出すること、を備える鉄道車両の走行位置検出システムにおいて、鉄道車両は、線路近くに設置された識別子の識別子情報を読み込んで、特定路線のうちの、鉄道車両が走行している特定路線区間を決定する特定路線決定手段と、複数の路線網を鉄道車両が走行したときの実測走行データを記憶する走行データ記憶手段と、特定路線決定手段が決定した特定路線区間を含む特定路線内の実測走行データを、走行データ記憶手段から読み出し、現在走行中に検知している現在走行データと比較することにより、自車の走行位置を検出すること、を特徴とする。
この場合には、複数路線網のうち、自車が走行している特定路線を決定し認識しているので、その特定路線内の実測走行データのみを読み出して、現在走行中に検知している現在走行データと比較すればよいため、処理を必要とするデータ量を大幅に低減でき、鉄道車両の傾斜制御をリアルタイムで適確に行うことができる。
(2)(1)に記載する鉄道車両の走行位置検出システムにおいて、鉄道車両が曲線部を通過したときに測定した実測曲線情報と予め記憶している多数の曲線部のデータ曲線情報とを比較して、どの曲線部を通過したかを特定し、特定した曲線部からの実測距離に基づいて自車の走行位置を検出すること、を特徴とする。
この場合には、複数路線網のうち、自車が走行している特定路線を決定し認識しているので、実測走行データのうち、その特定路線内のデータ曲線情報のみを読み出して、現在走行中に検知している現在走行データである実測曲線情報と比較すればよいため、処理を必要とするデータ量を大幅に低減でき、鉄道車両の傾斜制御をリアルタイムで適確に行うことができる。
(3)(2)に記載する鉄道車両の走行位置検出システムにおいて、曲線部を通過する毎に、測定した複数の実測値を含む通過曲線情報を作成する通過曲線情報作成部と、通過し得る全ての曲線部に対応していて、且つ前記データ曲線情報に基づいて前記実測値に対応するデータ値を含む候補曲線情報を作成する候補曲線情報作成部と、直前の通過曲線情報及び直前より前の通過曲線情報にそれぞれ対応する候補曲線情報の集まりを候補曲線情報群として複数個作成する候補曲線情報群作成部と、直前の通過曲線情報及び直前より前の通過曲線情報と各候補曲線情報群との一致度合いを求めるために、直前の通過曲線情報と候補曲線情報群の中で対応する候補曲線情報との一致度合いについて0から1までの範囲で大きい値ほど一致していることを表す単曲線一致率を演算すると共に、直前より前の通過曲線情報と候補曲線情報群の中で対応する候補曲線情報との一致度合いについて0から1までの範囲で大きい値ほど一致していることを表す単曲線一致率を演算して、各単曲線一致率を乗算して総合一致率を演算する一致率演算部と、前記候補曲線情報群毎にそれぞれ演算された総合一致率を比較して直前に通過した曲線部を特定する曲線特定部とを備えていることを特徴とする。
この場合には、一つの曲線部を通過したときに得られる通過曲線情報に基づいて一つの曲線部の一致度合いを比較するのではなく、複数の曲線部を通過したときに得られる複数の通過曲線情報に基づいて複数の曲線部の集まりの一致度合い(総合一致率)を比較している。すなわち、一つずつの曲線部の一致度合いを単曲線一致率として数値化し、各単曲線一致率を乗算することで総合一致率を演算しているため、正しい結果を示す総合一致率のみが大きな値になる。従って、直前に通過した曲線部の形状と、鉄道車両が通過し得る全ての曲線部の中で何れかの曲線部の形状とが似ていても、各総合一致率を比較することで、直前に通過した曲線部がどの曲線部なのかが的確に特定できる。
(4)(3)に記載された鉄道車両の走行位置検出システムにおいて、前記通過曲線情報作成部は、前記通過曲線情報の実測値として、通過する2つの曲線部間の曲線間長である実測曲線間長を含め、前記候補曲線情報作成部は、前記候補曲線情報のデータ値として、前記データ曲線情報に基づいて2つの曲線部間の曲線間長であるデータ曲線間長を含め、前記一致率演算部は、実測曲線間長とデータ曲線間長との一致度合いについて0から1までの範囲で大きい値ほど一致していることを表わす曲線間長一致率を演算し、且つ前記曲線間長一致率と前記各単曲線一致率とを乗算して総合一致率を演算することを特徴とする。
この場合には、曲線部間の曲線間長の一致度合いを曲線間長一致率として数値化して、曲線間長一致率と各単曲線一致率とを乗算して総合一致率を演算する。こうして、総合一致率を演算するためのパラメータとして曲線間長の一致度合いを含めることで、より正確な総合一致率を演算することができる。
(5)(3)又は(4)に記載された鉄道車両の走行位置検出システムにおいて、前記一致率演算部は、前記総合一致率を演算する前に、直前の通過曲線情報と各候補曲線情報との一致度合いについて0から1までの範囲で大きい値ほど一致していることを表わす単曲線一致率を演算して、前記候補曲線情報群作成部は、前記演算された単曲線一致率のうち所定値より大きい単曲線一致率が2つ以上存在する場合、各候補曲線情報群について前記所定値より大きい単曲線一致率が得られる原因になった候補曲線情報を含むように絞り込むことを特徴とする。
この場合には、候補曲線情報群を絞り込むことで、候補曲線情報群毎に演算する総合一致率の数が少なくなる。このため、総合一致率を演算するための計算量が少なくなり、計算負荷を軽減することができる。
(6)(3)乃至(5)の何れか1つに記載する鉄道車両の走行位置検出システムにおいて、前記候補曲線情報群毎に演算された総合一致率のうち基準値より大きい総合一致率が2つ以上存在する場合、前記通過曲線情報作成部は、それまでに用意していた通過曲線情報より1つ前の通過曲線情報を新たに用意し、前記候補曲線情報群作成部は、基準値より大きい総合一致率が得られる原因になった各候補曲線情報群に対して、前記新たに用意した通過曲線情報に対応する候補曲線情報を含めて、新しい各候補曲線情報群を作成し、前記一致率演算部は、前記新たに用意した通過曲線情報と前記新しい候補曲線情報群の中で対応する候補曲線情報との一致度合いについて0から1までの範囲で大きい値ほど一致していることを表す追加単曲線一致率を演算して、この追加単曲線一致率を前記総合一致率に乗算して新たに総合一致率を演算することを特徴とする。
この場合には、仮に直前に通過した曲線部及び直前より前に通過した曲線部と似ている複数の曲線部が存在していると、基準値より大きい総合一致率が2つ以上存在する可能性がある。このようなときに、単曲線一致率の数を増やして新たに総合一致率を演算することで、正しい結果を示す総合一致率のみが大きな値になる。こうして、繰り返し総合一致率を演算することで、直前に通過した曲線部を確実に特定することができる。
(7)(3)乃至(6)の何れか1つに記載する鉄道車両の走行位置検出システムにおいて、前記一致率演算部は、前記実測値と前記データ値との差から前記実測値又は前記データ値を除算した値を1から減算して部分一致率を算出し、前記実測値毎に算出した各部分一致率を乗算して、前記単曲線一致率を演算することを特徴とする。
(8)(7)に記載された鉄道車両の走行位置検出システムにおいて、個々の曲線部をパラメータの集合ととらえ、前記パラメータについて、0から1までの範囲で一致度合いを算出し、その乗算で前記単曲線一致率を演算することで個々の曲線部を評価することを特徴とする。
これらの場合には、一つの曲線部の一致度合いを意味する単曲線一致率を演算するために、相関演算せずに、曲線部のパラメータの各部分一致率を算出して、各部分一致率を乗算する。個々の曲線部をパラメータの集合ととらえ、曲線部を構成するパラメータについて、0.0〜1.0の範囲で一致度合いを算出し、その乗算で単一の曲線部としての一致度合いを算出することで、個々の曲線部を評価する。パラメータごとに、部分一致率に重みづけをすることで、変動の大きなパラメータの寄与度を調整することが可能となり、線路の状況に応じて変化しやすい値であっても、乗算の集合である単曲線一致率への影響を抑えることができる。従って、相関演算する場合のように、正確な結果を得るために何度か走行試験で取得した情報(データ曲率等)を加工して、常に更新する必要はなくて、各曲線部に対して予め記憶しているデータ曲線情報だけを利用して実施することができる。
(9)(7)又は(8)に記載された鉄道車両の走行位置検出システムにおいて、前記通過曲線情報作成部は、前記通過曲線情報の実測値として、実測曲率半径又は実測曲率の実測値を含め、前記候補曲線情報作成部は、前記候補曲線情報のデータ値として、データ曲率半径又はデータ曲率のデータ値を含めることを特徴とする。
(10)(9)に記載された鉄道車両の走行位置検出システムにおいて、前記通過曲線情報作成部は、前記通過曲線情報の実測値として、実測入口緩和曲線長と、実測円曲線長と、実測出口緩和曲線長とを含め、前記候補曲線情報作成部は、前記候補曲線情報のデータ値として、データ入口緩和曲線長と、データ円曲線長と、データ出口緩和曲線長とを含めることを特徴とする。
これらの場合には、走行中に測定する実測曲率半径又は実測曲率と、実測入口緩和曲線長と、実測円曲線長と、実測出口緩和曲線長は、その他に測定する値に比べて、真値に近い正確な値を測定できるものである。従って、各部分一致率を算出するための実測値として相応しく、より正確な単曲線一致率を演算することができる。
(11)(7)乃至(10)の何れか1つに記載する鉄道車両の走行位置検出システムにおいて、前記通過曲線情報作成部は、前記通過曲線情報の実測値として、実測曲線方向を含め、前記候補曲線情報作成部は、前記候補曲線情報のデータ値として、データ曲線方向を含めることを特徴とする。
この場合には、曲線部を特徴づける各パラメータに曲線方向を含めるため、曲線部を特定する際に、その曲線部の曲がっている方向が分かることになる。即ち、一つの曲線部であっても、例えば順行するときには右向きに曲がる曲線部であり、逆行するときには左向きに曲がる曲線部になるため、曲線部の曲がっている方向が分かることで、順行しているのか又は逆行しているのかが分かる。
(12)(2)に記載する鉄道車両の走行位置検出システムにおいて、曲線部を通過する毎に、測定した複数の実測値を含む通過曲線情報を作成する通過曲線情報作成部と、通過し得る全ての曲線部に対応していて、且つ前記データ曲線情報に基づいて前記実測値に対応するデータ値を含む候補曲線情報を作成する候補曲線情報作成部と、直前の通過曲線情報と各候補曲線情報との一致度合いを求めるために、前記実測値と前記データ値との差から前記実測値又は前記データ値を除算した値を1から減算して部分一致率を算出し、実測値毎に算出した各部分一致率を乗算して単曲線一致率を演算する一致率演算部と、前記候補曲線情報毎にそれぞれ演算された単曲線一致率を比較して直前に通過した曲線部を特定する曲線特定部とを備えていること、を特徴とする。
(13)(12)に記載された鉄道車両の走行位置検出システムにおいて、個々の曲線部をパラメータの集合ととらえ、前記パラメータについて、0から1までの範囲で一致度合いを算出し、その乗算で単曲線一致率を演算することで個々の曲線部を評価することを特徴とする。
これらの場合には、各候補曲線部と直前に通過した曲線部との一致度合いを求める際に、相関演算せずに、曲線部を特徴づける各パラメータの部分一致率を算出し、各部分一致率を乗算することで単曲線一致率を演算する。個々の曲線部をパラメータの集合ととらえ、曲線部を構成するパラメータについて、0.0〜1.0の範囲で一致度合いを算出し、その乗算で単一の曲線部としての一致度合いを算出することで、個々の曲線部を評価する。パラメータごとに、部分一致率に重みづけをすることで、変動の大きなパラメータの寄与度を調整することが可能となり、線路の状況に応じて変化しやすい値であっても、乗算の集合である単曲線一致率への影響を抑えることができる。このため、一つのパラメータが、鉄道車両や線路の状況に応じて変化し易い値であっても、乗算の集合である単曲線一致率にはほとんど影響はない。従って、相関演算する特許文献1の方法のように、正確な結果を得るために何度か走行試験で取得した情報を加工して、常にデータ曲率を更新する必要はない。すなわち、車体傾斜制御するために各曲線部に対して予め記憶しているデータ曲線情報だけを利用して実施することができ、事前にデータを準備する手間をほとんど無くすことができる。
本発明の鉄道車両の走行位置検出システムによれば、複数路線網のうち、自車が走行している特定路線を決定し認識しているので、その特定路線内の実測走行データのみを読み出して、現在走行中に検知している現在走行データと比較すればよいため、処理を必要とするデータ量を大幅に低減でき、鉄道車両の傾斜制御をリアルタイムで適確に行うことができる。
本実施形態の走行位置検出システムの構成を示した図である。 鉄道車両が自車の正確な走行位置を検出できない状況を説明するための図である。 本実施形態の路線状況を示した図である。 曲線部の構造を説明するための図である。 データベース処理部が各曲線部毎に記憶している情報値を模式的に示した図である。 曲線部を通過したときに検出されたヨー角速度を示した図である。 図6に示したヨー角速度が上昇する区間において、曲率比と曲率微分値を示した図である。 図6に示したヨー角速度が減少する区間において、曲率比と曲率微分値を示した図である。 本実施形態の走行位置検出システムの処理フローを示した図である。 単曲線一致率演算処理の処理フローを示した図である。 候補曲線情報作成部が作成した各候補曲線情報を模式的に示した図である。 本実施形態において演算された各単曲線一致率の値を示した図である。 候補曲線情報群作成処理の処理フローを示した図である。 候補曲線情報群作成部が用意する各候補曲線情報群を模式的に示した図である。 図14に示した各候補曲線情報群のうち絞り込まれた各候補曲線情報群を模式的に示した図である。 総合一致率演算処理の処理フローを示した図である。 本実施形態において演算された各総合一致率の値を示した図である。 新総合一致率演算処理の処理フローを示した図である。 相関演算する従来方法を説明するための図である。 単曲線一致率を演算する本実施形態の方法を説明するための図である。 特許文献1の地点検知方法において、実測曲率とデータ曲率との相関演算を説明するための模式図である。 本発明の走行位置検出システムの全体構成を示すフローチャートである。 別の実施形態の単曲線一致率を演算する本実施形態の方法を説明するための図である。
本発明に係る鉄道車両の走行位置検出システムの実施形態について、図面を参照しながら以下に説明する。図22に、走行位置検出システムの全体構成をフローチャートで示す。
始めに、地上子である識別子を読み込む(S1)。次に、複数の路線網の中から鉄道車両である自車が走行している特定路線区間を決定する(S2)。S2の特定路線区間の決定方法の詳細について、特許文献2を参照されたい。
次に、複数の路線網を鉄道車両が走行したときの実測走行データを記憶する走行データ記憶手段から、S2で決定された特定路線区間内の実測走行データを読み出す(S3)。
次に、現在の自車の現在走行データを取得する(S4)。次に、S3で読み出した特定路線区間内の実測走行データと、S4で取得した現在走行データとを比較して、自車の現在の走行位置を検出する(S5)。S5については、後で詳細に説明する。
本実施例によれば、鉄道車両が、複数の路線網から選択された特定路線区間を含む特定路線を走行すること、鉄道車両が自車の走行位置を検出すること、を備える鉄道車両の走行位置検出システムにおいて、鉄道車両は、線路近くに設置された識別子の識別子情報を読み込んで、特定路線のうちの、鉄道車両が走行している特定路線区間を決定する特定路線決定手段と、複数の路線網を鉄道車両が走行したときの実測走行データを記憶する走行データ記憶手段と、特定路線決定手段が決定した特定路線区間を含む特定路線内の実測走行データを、走行データ記憶手段から読み出し、現在走行中に検知している現在走行データと比較することにより、自車の走行位置を検出すること、を特徴とするので、複数路線網のうち、自車が走行している特定路線を決定し認識しているので、その特定路線内の実測走行データのみを読み出して、現在走行中に検知している現在走行データと比較すればよいため、処理を必要とするデータ量を大幅に低減でき、鉄道車両の傾斜制御をリアルタイムで適確に行うことができる。
次に、S5の一例について、図1〜図20を用いて詳細に説明する。
図1は、本実施形態の走行位置検出システムKSの構成を示した図である。走行位置検出システムKSは、図1に示すように、鉄道車両1に組み込まれたものであり、鉄道車両1は、車上子2と、速度発電機3と、ヨー角速度センサ4とを備えている。
走行位置検出システムKSは、実測距離演算部10と、走行速度演算部20と、データベース処理部30と、通過曲線情報作成部40と、候補曲線情報作成部50と、候補曲線情報群作成部50Aと、一致率演算部60と、曲線特定部70と、実測距離補正部80とを備えていて、周知のCPU、RAM、ROM、入出力インタフェース、及びこれらを接続するバスからなる電子制御装置によって実現されている。
鉄道車両1が走行する線路には、多数の地上子Tが配置されている。各地上子Tは、例えばトランスポンダ等によって構成され、線路に沿って間隔を隔てて任意の位置に設置されていて、固有の識別値を有している。車上子2は、鉄道車両1の車体の底部に取付けられていて、例えば、トランスポンダ等によって構成され、走行中に地上子Tと対向することで、地上子Tから送信される識別値を検知して取得するようになっている。取得された地上子Tの識別値は、データベース処理部30に送信される。
速度発電機3は、鉄道車両1の走行速度V(km/h)を測定するためのものであり、車輪の車軸に取付けられていて、車軸の回転に伴って発生するパルス又は正弦波の出力信号を走行速度演算部20に出力する。なお、速度発電機3の出力の方式は特に限定されるものではない。ヨー角速度センサ4は、車体に作用するヨー角の時間的変化をヨー角速度θ(rad/s)として検出する。検出されたヨー角速度θは、通過曲線情報作成部40に送信される。
実測距離演算部10は、走行速度演算部20で演算された走行速度Vと経過時間を乗算することで、鉄道車両1が走行した距離を実測距離S(m)として演算する。演算された実測距離Sは、通過曲線情報作成部40及び実測距離補正部80に送信される。走行速度演算部20は、速度発電機3から入力する出力信号をカウントすることで、周知の方法によって走行速度Vを演算する。演算された走行速度Vは、通過曲線情報作成部40に送信される。なお、実測距離Sを演算する方法は、上記した方法に限定されるものではなく適宜変更可能であり、例えば加速度センサを用いて演算しても良い。また、走行速度Vを演算する方法も適宜変更可能であり、例えば実測距離演算部10から入力する実測距離Sを微分して演算しても良い。
データベース処理部30は、地上子Tの正確な位置情報、即ち地上子Tの絶対位置や、地上子T間の絶対距離等を予め記憶している。データベース処理部30は、送信される地上子Tの識別値に基づいて、検知された地上子Tの位置情報を検索して、その地上子Tを特定する。こうして、鉄道車両1が地上子Tを検知した場合、地上子Tの位置情報が実測距離補正部80に送信され、実測距離補正部80は、地上子Tの位置情報に基づいて実測距離Sを補正して、現地点での自車の走行位置を特定する。そして、車体傾斜制御部90は、データベース処理部30から送信される曲線情報(カント、曲率半径、曲線方向等)と特定された自車の走行位置とに基づいて、車体を傾ける適切な地点を決定して、車体傾斜制御を適切に実行するようになっている。
ところで、鉄道車両1は、天候状態や電波状況が悪い場合、機器が故障した場合等に、地上子Tを検知できない事態が生じ得る。万一このような事態が生じた場合、地上子Tの正確な位置情報を取得できないため、実測距離Sを補正することができず、自車の正確な走行位置を検出することができない。特に、図2に示すように、鉄道車両が線路A1を走行して分岐駅を通過した後、仮に地上子Tを検知できなかった場合に、鉄道車両1は自車の正確な走行位置を検出できず、自身が線路A2又は線路A3のどちらを走行しているのかが分からなくなる。その結果、車体傾斜制御が誤った地点で実行されるおそれがある。
そこで、本実施形態の走行位置検出システムKSは、万一地上子Tが検知できなかった場合であっても、線路の曲線部を通過するときの実測値に基づいて、通過した曲線部を特定して、自車の走行位置を検出することができるように構成されている。以下、走行位置検出システムKSの各処理部30〜80の構成について、図3〜図18を参照しながら説明する。図3は、本実施形態の路線状況を示した図である。
鉄道車両の路線状況は各国に応じて事情が異なるが、本実施形態では図3に示した路線状況を例にして説明する。図3に示した路線状況では、2本の線路がおおよそ平行に並んで南北方向に延びていて、西側に配置された線路を「西線W」と呼び、東側に配置された線路を「東線E」と呼ぶことにする。そして、西線Wには3つの曲線部W1,W2,W3が形成され、東線Eにも3つの曲線部E1,E2,E3が形成されている。
鉄道車両1は、これら西線W及び東線Eの両方の線路において、北側及び南側の双方向に向かって走行するようになっている。即ち、鉄道車両1が走行する路線状況は、単線並列による双方向運用になっている。ここでは、北側に向かって走行する場合を「順行」と呼び、南側に向かって走行する場合を「逆行」と呼ぶことにする。また、曲線部W1と曲線部W2の間の直線部分と、曲線部E1と曲線部E2の間の直線部分との間には、渡り線C1,C2が配置されていて、鉄道車両1は渡り線C1,C2を通って西線Wから東線Eへ、又は東線Eから西線Wへ移ることができるようになっている。
次に、各曲線部E1,E2,E3,W1,W2,W3の構造について、図4を参照しながら、曲線部E1を例にして説明する。曲線部E1は、緩和曲線開始位置p1と円曲線開始位置p2と円曲線終了位置p3と緩和曲線終了位置p4とを有し、緩和曲線開始位置p1と円曲線開始位置p2との間が入口緩和曲線であり、円曲線開始位置p2と円曲線終了位置p3との間が円曲線であり、円曲線終了位置p3と緩和曲線終了位置p4との間が出口緩和曲線である。なお、入口及び出口の名称は、図4の右側に向かって走行する場合(順行)を基準にしている。また、円曲線は曲線部E1のうち曲率がほぼ最大で一定値になる部分であり、入口緩和曲線は直線部分から円曲線に向かって緩やかに曲がる部分であり、出口緩和曲線は円曲線から出口緩和曲線に向かって緩やかに曲がる部分である。
ここで、データベース処理部30が予め記憶している情報値について、図5を参照しながら説明する。図5は、データベース処理部30が各曲線部E1,E2,E3,W1,W2,W3毎に記憶している情報値を模式的に示した図である。先ず、「a」の文字は、円曲線の曲率半径を意味している。そして、「b」の文字は、入口緩和曲線の長さを意味していて、「c」の文字は、円曲線の長さを意味していて、「d」の文字は、出口緩和曲線の長さを意味している。
こうして、各曲線部E1,E2,E3,W1,W2,W3のそれぞれにおいて、緩和曲線開始位置p1、円曲線開始位置p2、円曲線終了位置p3、緩和曲線終了位置p4、曲率半径a、入口緩和曲線長b、円曲線長c、出口緩和曲線長dの正確な情報値を、鉄道車両1が事前の走行試験によって測定して、データベース処理部30が記憶している。例えば、曲線部E1に対応して記憶しているデータ曲線情報DE1は、「E1」の文字が付いていて、緩和曲線開始位置E1p1、円曲線開始位置E1p2、円曲線終了位置E1p3、緩和曲線終了位置E1p4、曲率半径E1a、入口緩和曲線長E1b、円曲線長E1c、出口緩和曲線長E1dである。なお、上記した各情報値は、車体傾斜制御を実行するためにデータベース処理部30が予め記憶しているものである。
続いて、鉄道車両1の走行状況について説明する。本実施形態では、図3に示すように、鉄道車両1が、東線Eを二点鎖線で示した位置から順行して、曲線部E1を通過して、渡り線C1を通って西線に移って、曲線部W2を通過した後、現時点では曲線部W2と曲線部W3の間の地点を走行している状況である。しかし、鉄道車両1は、地上子Tを検出できなくて、現地点では自身がどの曲線部を通過したかが分からないという状況で説明することにする。本実施形態では、データベース処理部30が地上子Tを検出できなかったことを判断すると、鉄道車両1は、地上子Tの位置情報に基づく走行位置の検出を中断し、本実施形態の走行位置検出システムKSを実行するようになっている。なお、後述する走行位置検出システムKSは、地上子による地点検知のバックアップのシステムとして常時稼働していて、地点検知を行う機器の故障を検知した場合又は地点検知による現在の走行位置の補正が失敗した場合に、地上子による地点検知から切り替わって実行しても良い。
走行位置検出システムKSの各処理部40〜80の説明に戻る。通過曲線情報作成部40は、随時ヨー角速度センサ4からヨー角速度θを入力されて、曲線部を通過する度に、通過した曲線部に対応する通過曲線情報Rを作成するものである。この通過曲線情報Rは、走行中に測定された値である実測値を複数含むものであり、具体的には、実測曲率半径az、実測入口緩和曲線長bz、実測円曲線長cz、実測出口緩和曲線長dz、実測曲線方向ez、実測曲線間長fzをそれぞれ実測値として含んでいる。ここで、実測曲線間長fzとは、2つの曲線部を通過したときに、一方の曲線部の緩和曲線終了位置p4から他方の曲線部の緩和曲線開始位置p1までの実測距離Sである。なお、上記した各実測値の最後の文字には、測定した実測値であることを意味する「z」を付して区別している。
本実施形態では、図3に示すように、直前に通過した曲線部(実際には、曲線部W2)の通過曲線情報を「R1」で表わし、各実測値を、実測曲率半径R1az、実測入口緩和曲線長R1bz、実測円曲線長R1cz、実測出口緩和曲線長R1dz、実測曲線方向R1ez、実測曲線間長R1fzで表わすことにする。実測曲線間長R1fzは、曲線部E1の緩和曲線終了位置p4から曲線部W2の入口緩和曲線開始位置p1までの実測距離である。
また、直前より1つ前に通過した曲線部(実際には、曲線部E1)の通過曲線情報を「R2」で表わし、各実測値を、実測曲率半径R2az、実測入口緩和曲線長R2bz、実測円曲線長R2cz、実測出口緩和曲線長R2dz、実測曲線方向R2ez、実測曲線間長R2fzで表わすことにする。また、直前より2つ前に通過した曲線部の通過曲線情報を「R3」で表わし、各実測値を、実測曲率半径R3az、実測入口緩和曲線長R3bz、実測円曲線長R3cz、実測出口緩和曲線長R3dz、実測曲線方向R3ez、実測曲線間長R3fzで表わすことにする。なお、通過曲線情報作成部40は、鉄道車両1が新たに曲線部を通過したときに測定した通過曲線情報Rを、直前の通過曲線情報R1として、それまでの通過曲線情報R1,R2,R3,・・・を一つ前の通過曲線情報R2,R3,・・・に更新する。
上述した各実測値の測定方法について説明する。図6は、曲線部を通過したときに検出されたヨー角速度θを示した図である。図6に示すように、ヨー角速度θは波形状で表わされ、曲線部の円曲線を通過中にヨー角速度θが最大値θ1となる。このため、このヨー角速度θ1を走行速度Vで除算することによって、円曲線の曲率が算出される。この円曲線の曲率は、曲線部の最大の曲率であるため、その曲率の逆数が実測曲率半径azになる。
ここで、図6に示したヨー角速度θが増加する区間Aにおいて、ヨー角速度θを走行速度Vで除算して曲率を算出し、その曲率と円曲線の曲率(最大の曲率)との曲率比αが、図7の実線で示されている。そして、図7において、波形状に示された曲率比αの時間微分が、曲率微分値βとして破線で示されている。
緩和曲線開始位置は、曲率微分値βが「0」から上昇し始めた位置であるため、本実施形態では、曲率微分値βが増加中で0.1になる点を求めて、図7のX1で示した位置を緩和曲線開始位置に決定している。また、円曲線開始位置は、曲率微分値βが「0」まで減少する直前の位置であるため、本実施形態では、曲率微分値が減少中で0.1になる点を求めて、図7のX2で示した位置を円曲線開始位置に決定している。なお、緩和曲線開始位置及び円曲線開始位置を決定するために、曲率微分値βの値が0.1になる位置を求めたが、例えば曲率微分値βの値が0.05になる位置を求めても良く、適宜変更可能である。
また、図6に示したヨー角速度θが減少する区間Bにおいて、ヨー角速度θを走行速度Vで除算して曲率を算出し、その曲率と円曲線の曲率(最大の曲率)との曲率比αが、図8の実線で示されている。そして、図8において、波形状に示された曲率比αの時間微分が、曲率微分値βとして破線で示されている。
円曲線終了位置は、曲率微分値が「0」から減少し始めた位置であるため、本実施形態では、曲率微分値βが減少中で−0.1になる点を求めて、図8のX3で示した点を円曲線終了位置に決定している。また、緩和曲線終了位置は、曲率微分値βが「0」まで上昇する直前の位置であるため、本実施形態では、曲率微分値βが増加中で−0.1になる点を求めて、図8のX4で示した位置を緩和曲線終了位置に決定している。なお、円曲線終了位置及び緩和曲線終了位置を求めるために、曲率微分値βの値が−0.1になる位置を求めたが、例えば曲率微分値βが−0.05になる位置を求めても良く、適宜変更可能である。
こうして、緩和曲線開始位置、円曲線開始位置、円曲線終了位置、緩和曲線終了位置が分かることによって、緩和曲線開始位置と円曲線開始位置との間の実測距離Sが実測入口緩和曲線長bzになり、円曲線開始位置と円曲線終了位置との間の実測距離Sが実測円曲線長czになり、円曲線終了位置と緩和曲線終了位置との間の実測距離Sが実測出口緩和曲線長dzになる。そして、曲線部を通過したときに検出するヨー角速度θの正負によって、右向きに曲がる曲線部又は左向きに曲がる曲線部が分かるため、実測曲線方向ezが分かる。なお、実測曲線方向ezは、例えば、右向きに曲がる曲線部の場合に「+1」に定義し、左向きに曲がる曲線部の場合に「−1」に定義すれば良い。また、通過した2つの曲線部において、一方の曲線部の緩和曲線終了位置と他方の曲線部の緩和曲線開始位置が分かるため、実測曲線間長fzが分かる。
上述したように、曲率(曲率微分値β)を利用して、緩和曲線開始位置、円曲線開始位置、円曲線終了位置、緩和曲線終了位置を特定する方法は、曲率がヨー角速度θを走行速度Vで除算したものであるため、曲率の値に走行速度Vの影響が出ないというメリットがある。一方、緩和曲線開始位置、円曲線開始位置、円曲線終了位置、緩和曲線終了位置を特定する別の方法として、ヨー角速度θをそのまま利用しても良い。この場合には、ヨー角速度θを走行速度Vで除算して曲率を算出しないため、計算量が減り、ソフトウェアの構造が簡単になるというメリットがある。なお、緩和曲線開始位置等を特定する方法は、ヨー角速度センサ4が検出したヨー角速度θを用いる方法以外にも適宜変更可能であり、例えば加速度センサが検出した鉄道車両1の横方向に作用する横加速度を用いても良い。
次に、走行位置検出システムKSの処理フローについて、図9を参照しながら説明する。図9に示すように、本実施形態において、走行位置検出システムKSが実行されると、単曲線一致率演算処理(ステップ100)、候補曲線情報群作成処理(ステップ200)、総合一致率演算処理(ステップ300)が順番に実行される。そして、後述するように、演算された総合一致率Jの値に応じて、走行位置検出システムKSを終了する、又は新総合一致率演算処理(ステップ400)を実行する、或いは再び単曲線一致率演算処理(ステップ100)を実行するようになっている。
<単曲線一致率演算処理>
図10に示すように、単曲線一致率演算処理では、先ず、候補曲線情報作成部50が、データベース処理部30に記憶された全ての曲線部、即ち鉄道車両1が通過し得る全ての曲線部に対応する候補曲線情報Kを作成する(ステップ110)。続いて、通過曲線情報作成部40が、直前に通過した曲線部の通過曲線情報R1(図3参照)を用意する(ステップ120)。なお、以下の説明では、候補曲線情報Kを作成する対象の曲線部が、路線に存在する全ての曲線部(曲線部W1,W2,W3,E1,E2,E3)としているが、直前に走行している路線が分かっている場合においてその路線から遷移可能な路線に存在する曲線部としても良い。
ここで、図11は、候補曲線情報作成部50が作成した候補曲線情報Kを模式的に示した図である。本実施形態では、図11に示すように、候補曲線情報作成部50は、各曲線部W1,W2,W3,E1,E2,E3に対応して、候補曲線情報KW1n,KW2n,KW3n,KW1s,KW2s,KW3s,KE1n,KE2n,KE3n,KE1s,KE2s,KE3sを作成する。先ず、これら候補曲線情報Kの種類について説明する。
各候補曲線情報Kでは、仮に鉄道車両1が西線Wを南側から順行(北側に向かって走行)した場合を想定すると、曲線部W1,W2,W3を順番に通過するため(図3参照)、曲線部W1,W2,W3の順番に、候補曲線情報KW1n,KW2n,KW3nを対応させている。なお、候補曲線情報Kの最後の文字には、北側に向かって走行した場合を意味するnを付して区別している。例えば、候補曲線情報KW3nは、曲線部W3を北側に向かって通過した場合の候補曲線情報Kを意味している。
また、仮に鉄道車両1が西線Wを北側から逆行(南側に向かって走行)した場合を想定すると、曲線部W3,W2,W1を順番に通過するため、曲線部W3,W2,W1の順番に、候補曲線情報KW1s,KW2s,KW3sを対応させている。なお、候補曲線情報Kの最後の文字には、南側に向かって走行した場合を意味するsを付して区別している。例えば、候補曲線情報KW3sは、曲線部W1を南側に向かって通過した場合の候補曲線情報Kを意味している。
東線Eでも同様に、仮に鉄道車両1が東線Eを南側から順行した場合を想定すると、曲線部E1,E2,E3を順番に通過するため、曲線部E1,E2,E3の順番に、候補曲線情報KE1n,KE2n,KE3nを対応させている。また、仮に鉄道車両1が東線Eを北側から逆行した場合を想定すると、曲線部E3,E2,E1を順番に通過するため、曲線部E3,E2,E1の順番に、候補曲線情報KE1s,KE2s,KE3sを対応させている。
次に、各候補曲線情報Kに含まれている各データ値について説明する。各データ値は、曲線部における曲率半径a、入口緩和曲線長b、円曲線長c、出口緩和曲線長d、曲線方向e、曲線間長fであるが、これらのデータ値は、曲線部を通過するときの方向(順行又は逆行)を考慮した値である。これは、例えば曲線部E1を順行する場合と逆行する場合とでは、曲率半径a、入口緩和曲線長b、円曲線長c、出口緩和曲線長d、曲線方向e、曲線間長fが異なるためである。本実施形態では、候補曲線情報Kに含まれるデータ値について、例えば「データ曲率半径」のように、「データ」を付けて区別することにする。
例えば、上述したように、候補曲線情報KW3nは、曲線部W3を北側に向かって通過した場合の候補曲線情報Kを意味しているため、候補曲線情報作成部50は、候補曲線情報KW3nの各データ値として、曲線部W3のデータ曲線情報DW3の各情報値(図5参照)に基づいて、データ曲率半径W3a−n、データ入口緩和曲線長W3b−n、データ円曲線長W3c−n、データ出口緩和曲線長W3d−n、データ曲線方向W3e−nを作成する。なお、データ曲線方向W3e−nは、曲線部W3を北側に向かって通過する場合の曲線方向として、例えば左向きに曲がることを意味する「−1」が入る。また、候補曲線情報作成部50は、データ曲線情報DW3とデータ曲線情報DW2に基づいて、曲線部W2の緩和曲線終了位置p4から曲線部W3の緩和曲線開始位置p1までの距離であるデータ曲線間長W3f−nを作成する。
同様に、例えば、上述したように、候補曲線情報KW3sは、曲線部W1を南側に向かって通過した場合の候補曲線情報Kを意味しているため、候補曲線情報作成部50は、候補曲線情報KW3sの各データ値として、曲線部W1のデータ曲線情報DW1の各情報値に基づいて、データ曲率半径W3a−s、データ入口緩和曲線長W3b−s、データ円曲線長W3c−s、データ出口緩和曲線長W3d−s、データ曲線方向W3e−sを作成する。また、候補曲線情報作成部50は、曲線部W1のデータ曲線情報DW1及び曲線部W2のデータ曲線情報DW2に基づいて、曲線部W2の緩和曲線終了位置p4から曲線部W1の緩和曲線開始位置p1までの距離であるデータ曲線間長W3f−sを作成する。こうして、候補曲線情報作成部50は、走行し得る全ての曲線部E1〜W3に対応する各候補曲線情報KW1n〜KE3sを作成する。
単曲線一致率演算処理の説明に戻る。図10に示すように、ステップ110及びステップ120の後に、一致率演算部60は、部分一致率Hを算出して(ステップ130)、単曲線一致率Iを演算する(ステップ140)。ここで、単曲線一致率Iとは、直前の通過曲線情報R1と各候補曲線情報KW1n〜KE3sとが、0から1までの範囲でどのくらい一致しているかを数値で表わしたものである。単曲線一致率Iが1に近い程、直前の通過曲線情報R1と対象となる候補曲線情報Kとが一致していることを意味する。この単曲線一致率Iを演算するために、直前の通過曲線情報R1の各実測値と、対象となる候補曲線情報Kの各データ値とによって、部分一致率Hが算出される。
部分一致率Hを求めるための演算式は、以下の数1で表わされる。
Figure 2016141300

上記した数1で算出される部分一致率Hは、0から1までの範囲でどのくらい実測値とデータ値とが一致しているかを数値で表わしたものであり、部分一致率Hが大きい程、実測値とデータ値とが一致していることを意味する。上記した数1の場合、実測値がデータ値の2倍より大きい値であると、部分一致率Hが0より小さい値になる。従って、本実施形態では、部分一致率Hが0より小さい場合には、「0」と定義するようになっている。
例えば、直前の通過曲線情報R1と候補曲線情報KW1nとの単曲線一致率IW1nを算出するためには、直前の通過曲線情報R1の各実測値(実測曲率半径R1az、実測入口緩和曲線長R1bz、実測円曲線長R1cz、実測出口緩和曲線長R1dz、実測曲線方向R1ez)と、候補曲線情報KW1nの各データ値(データ曲率半径W1a−n、データ入口緩和曲線長W1b−n、データ円曲線長W1c−n、データ出口緩和曲線長W1d−n、データ曲線方向W1e−n)とによって、各部分一致率H(Ha,Hb,Hc,Hd,He)が算出される。
すなわち、実測曲率半径R1azとデータ曲率半径W1a−nとの部分一致率Haは、上記した数1の実測値に実測曲率半径R1azが代入され、上記した数1のデータ値にデータ曲率半径W1a−nが代入されて、算出される。同様に、実測入口緩和曲線長R1bzとデータ入口緩和曲線長W1b−nとの部分一致率Hb、実測円曲線長R1czとデータ円曲線長W1c−nとの部分一致率Hc、実測出口緩和曲線長R1dzとデータ出口緩和曲線長W1d−nとの部分一致率Hd、実測曲線方向R1ezとデータ曲線方向W1e−nとの部分一致率Heが算出される。なお、実測曲線方向R1ezとデータ曲線方向W1e−nとの部分一致率Heは、例えば、上記した数1を用いずに、同じ値(同じ曲線方向)である場合には「1」になり、異なる値(逆の曲線方向)である場合には「−1」になるように設定しても良い。そして、これら各部分一致率Ha,Hb,Hc,Hd,Heが全て乗算されることで、直前の通過曲線情報R1と候補曲線情報KW1nとの単曲線一致率IW1n(図12参照)が演算される。なお、単曲線一致率Iは1つの曲線部の一致度合いを表すものであるため、単曲線一致率IW1nを演算する際に、直前の通過曲線情報R1に含まれる実測曲線間長R1fz及び候補曲線情報KW1nに含まれるデータ曲線間長W1f−nが用いられることはない。
すなわち、一つの曲線部の一致度合いを意味する単曲線一致率IW1nを演算するために、相関演算せずに、曲線部のパラメータの各部分一致率Ha,Hb,Hc,Hd,Heを算出して、各部分一致率Ha,Hb,Hc,Hd,Heを乗算する。個々の曲線部をパラメータの集合ととらえ、曲線部を構成するパラメータについて、0.0〜1.0の範囲で一致度合いを算出し、その乗算で単曲線一致率IW1nを演算することで、個々の曲線部を評価する。パラメータごとに、部分一致率Hに重みづけをすることで、変動の大きなパラメータの寄与度を調整することが可能となり、線路の状況に応じて変化しやすい値であっても、乗算の集合である単曲線一致率Iへの影響を抑えることができる。従って、相関演算する場合のように、正確な結果を得るために何度か走行試験で取得した情報(データ曲率等)を加工して、常に更新する必要はなくて、各曲線部に対して予め記憶しているデータ曲線情報Dだけを利用して実施することができる。
以下、同様に、図11に示す全ての候補曲線情報KW1n〜KE3sに対して、単曲線一致率IW1n〜IE3sがそれぞれ演算される。ここで、図12は、本実施形態において演算された単曲線一致率IW1n〜IE3sの値を示した図である。図12に示すように、例えば、単曲線一致率IW1nは0.1であるため、直前の通過曲線情報R1と候補曲線情報KW1nとがほとんど似ていないことを意味している。一方、単曲線一致率IW2nは0.95であるため、直前の通過曲線情報R1と候補曲線情報KW2nとが非常に似ていることを意味し、単曲線一致率IE2nも0.9であるため、直前の通過曲線情報R1と候補曲線情報KE2nとが非常に似ていることを意味している。また、例えば、単曲線一致率IW1sは0であり、直前の通過曲線情報R1と候補曲線情報KW1sとが全く似ていないことを意味しているが、これは曲線方向が逆であるためである。
単曲線一致率演算処理の説明に戻る。図10に示すように、ステップ140の後、曲線特定部70は、各単曲線一致率IW1n〜IE3nが所定値Ia以上であるか否かを判断して、所定値Ia以上の単曲線一致率Iのみを残す(ステップ150)。所定値Iaは、直前の通過曲線情報R1と各候補曲線情報KW1n〜KE3sとが一致していることを判断するための基準になる値であり、1に極めて近い値に設定されている。本実施形態では、所定値Iaは0.9に設定されているが、適宜変更可能である。このため、図11に示すように、単曲線一致率IW2n及び単曲線一致率IE2nのみが0.9以上であるため、ステップ150において、曲線特定部70は、単曲線一致率IW2n及び単曲線一致率IE2nを残す。
次に、ステップ160において、曲線特定部70は、残った単曲線一致率Iが有るか否かを判断する。本実施形態では、単曲線一致率IW2n及び単曲線一致率IE2nが残っているため、「Yes」と判断されて、ステップ170に進む。続いて、ステップ170において、曲線特定部70は、残った単曲線一致率Iが2つ以上有るか否かを判断する。本実施形態では、「Yes」と判断されて、単曲線一致率演算処理が終了し、候補曲線情報群作成処理が実行される。
ところで、上述したように残った単曲線一致率Iが2つ以上存在する場合とは、直前の通過曲線情報R1が各候補曲線情報KW1n〜KE3sのうち2つ以上の候補曲線情報Kと非常に似ている、即ち直前に通過した曲線部が全ての曲線部の中で2つ以上の曲線部と非常に似ている場合を意味している。具体的に、本実施形態では、直前に通過した曲線部が、順行(北側に向かって走行)する場合の曲線部W2、及び順行する場合の曲線部E2と非常に似ている場合を意味している。
このように、曲線部W2の形状と曲線部E2と形状とが似ていると、直前に通過した曲線部が、曲線部W2又は曲線部E2のどちらなのかが分からない事態が生じる。要するに、複数本の線路が並行して敷設された区間を走行する場合、例えば日本のように片側通行が主流である状況においては、形状が類似した曲線部が検出され難くて、特定する曲線部を誤認することが少ない。しかしながら、本実施形態で説明するように単線双方向運用を行う状況においては、並行する複数の線路上に形状が類似した曲線部が存在し、それらを区別する必要がある。このため、本実施形態の走行位置検出システムKSでは、直前に通過した曲線部の一致度合いを演算するだけでなく、直線より前に通過した曲線部の一致度合いを演算することで、直前に通過した曲線部がどの曲線部なのかを特定できるように構成されている。
ここで、単曲線一致率演算処理のステップ180、ステップ190、ステップ195について説明しておく。本実施形態と異なる場合として、仮に、各単曲線一致率IW1n〜IE3sが全て所定値Iaより小さい結果が得られた場合、ステップ160において、残った単曲線一致率Iが無いため、「No」と判断されて、ステップ180に進むことになる。これは、直前の通過曲線情報R1が各候補曲線情報KW1n〜KE3sの全てに似ていない、即ち直前に通過した曲線部が全ての曲線部の中で何れも似ていないという結果が得られた場合の処理である。
上記した仮の場合では、ステップ180において、曲線特定部70が直前に通過した曲線部を特定できないため、鉄道車両1が次の曲線部を通過するまで走行位置検出システムKSの実行を待機する。これは、直前に通過した曲線部を誤って特定することを回避するためである。従って、ステップ180に進んだ段階では、次の曲線部を通過するまで待機して、鉄道車両1が新たに曲線部を通過した段階でステップ120に戻る。そして、新たに用意する直前の通過曲線情報R1に基づいて上述したステップ120〜160の処理が実行されるようになっている。
また、本実施形態と異なる場合として、仮に、各単曲線一致率IW2nのみが所定値Iaより大きい結果が得られた場合について説明する。この場合、ステップ170において、残った単曲線一致率Iが2つ以上ではないため、「No」と判断されて、ステップ190に進むことになる。これは、直前の通過曲線情報R1が候補曲線情報KW2nのみに非常に似ていて、即ち直前に通過した曲線部が順行(北側に向かって走行)する場合の曲線部W2にのみ非常に似ている場合の処理である。上記した仮の場合では、ステップ190において、曲線特定部70は、鉄道車両1自身が北側に向かって走行して、直前に曲線部W2を通過したことを特定することになる。
その後、ステップ191において、実測距離補正部80は、上述したように特定された曲線部W2に基づいて、実測距離Sを補正する。即ち、鉄道車両1は、直前に曲線部W2を北側に向かって走行したことが分かると共に、その曲線部W2の緩和曲線終了位置p4(図4参照)が分かるため、実測距離補正部80が現地点での実測距離Sを、曲線部W2の緩和曲線終了位置p4から測定された実測距離Sを用いて補正する。こうして、鉄道車両1は、直前に通過した曲線部を特定することで、自身の走行位置を検出できるようになっている。ステップ191の後、ステップ192に進み、走行位置検出システムKSが終了する。
<候補曲線情報群作成処理>
しかしながら、上述したように、本実施形態では、直前に通過した曲線部が曲線部W2又は曲線部E2のどちらなのかが分からなかったため、単曲線一致率演算処理後に、候補曲線情報群作成処理が実行される。先ず、図13に示すように、ステップ210において、通過曲線情報作成部40は、直前より1つ前に通過した曲線部の通過曲線情報R2(図3参照)を用意する。次に、ステップ220において、候補曲線情報群作成部50Aは、連続して通過する2つの曲線部に対応する2つの候補曲線情報Kの組み合わせを、各候補曲線情報群Gとして用意する。
ここで、図14は、候補曲線情報群作成部50Aが用意する各候補曲線情報群Gを模式的に示した図である。直前より前に通過した曲線部と、直前に通過した曲線部との組み合わせは、曲線部W1→W2、曲線部W2→W3、曲線部W3→W2、曲線部W2→W1、曲線部E1→E2、曲線部E2→E3、曲線部E3→E2、曲線部E2→E1が考えられる。更に、渡り線C1,C2(図3参照)を通過する場合を含めると、曲線部W1→E2、曲線部W1→E1(順行から逆行)、曲線部W2→E1、曲線部W2→E2(逆行から順行)、曲線部E1→W2、曲線部E1→W1(順行から逆行)、曲線部E2→W1、曲線部E2→W2(逆行から順行)が考えられる。
そして、図14に示すように、上述した曲線部の組み合わせに各候補曲線情報Kが対応していて、2つの候補曲線情報Kの集まりが候補曲線情報群Gになっている。例えば、候補曲線情報群G1(W1→W2)は、直前より前に通過した曲線部が曲線部W1で、直前に通過した曲線部が曲線部W2である場合を想定して、候補曲線情報KW1nと候補曲線情報KW2nを含んでいる。以下、同様に、候補曲線情報群G2(W2→W3)〜G16(E2→W2)が用意される。
また、本実施形態では、各候補曲線情報群G1(W1→W2)〜G16(E2→W2)には、それぞれデータ曲線間長Fが含まれるようになっている。例えば、候補曲線情報群G1(W1→W2)は、曲線部W1を順行して通過した後に、曲線部W2を順行して通過した場合を想定しているため、曲線部W1の緩和曲線終了位置p4から曲線部W2の緩和曲線開始位置p1までのデータ曲線間長F1を含んでいる。なお、このデータ曲線間長F1は、候補曲線情報KW2nに含まれるデータ曲線間長W2f−nと同じものである。
また、例えば、候補曲線情報群G9(W1→E2)は、曲線部W1を順行して通過した後に、渡り線C2を通って東線Eに移り、曲線部E2を順行して通過した場合を想定しているため、曲線部W1の緩和曲線終了位置p4から曲線部E2の緩和曲線開始位置p1までのデータ曲線間長F9を含んでいる。なお、このデータ曲線間長F9は、曲線部W1の緩和曲線終了位置p4から渡り線C2の入口位置までの長さと、渡り線C2の長さと、渡り線C2の出口位置から曲線部E2の緩和曲線開始位置p1までの長さに基づいて作成されたものである。
候補曲線情報群作成処理の説明に戻る。本実施形態では、ステップ230において、候補曲線情報群作成部50Aが、各候補曲線情報群G1(W1→W2)〜G16(E2→W2)を絞り込むようになっている。具体的には、単曲線一致率演算処理で、直前に通過した曲線部が順行する場合の曲線部W2又は曲線部E2のどちらかであることが分かっているため、候補曲線情報群Gは、各候補曲線情報群G1(W1→W2)〜G16(E2→W2)の中で、直前に通過した曲線部に対応する候補曲線情報Kが候補曲線情報KW2n又は候補曲線情報KE2nである候補曲線情報群Gに絞り込む。この結果、図14に示す各候補曲線情報群G1(W1→W2)〜G16(E2→W2)は、図15に示す各候補曲線情報群G1(W1→W2),G5(E1→E2),G9(W1→E2),G12(W2→E2),G13(E1→W2),G16(E2→W2)に絞り込まれる。こうして、各候補曲線情報群Gが絞り込まれた後に、候補曲線情報群作成処理が終了し、総合一致率演算処理が実行される。
<総合一致率演算処理>
総合一致率演算処理では、先ず、図16に示すように、ステップ310において、一致率演算部60は、絞り込まれた候補曲線情報群G毎に、直前より前の曲線部の一致度合いを意味する単曲線一致率IIを演算する。この単曲線一致率IIは、直前より前の通過曲線情報R2と、対象になる候補曲線情報Kとが、0から1までの範囲でどのくらい一致しているかを数値で表わしたものである。単曲線一致率IIが1に近い程、直前より前の通過曲線情報R2と、対象になる候補曲線情報Kとが一致していることを意味する。本実施形態において、対象になる候補曲線情報Kは候補曲線情報KW1n,KE1n,KW2s,KE2sであるため(図15参照)、単曲線一致率IIを演算するために、上記した数1を用いて、直前の通過曲線情報R2の各実測値と各候補曲線情報のKW1n,KE1n,KW2s,KE2sの各データ値とによって、部分一致率Hが演算される。
例えば、直前より前の通過曲線情報R2と候補曲線情報KW1nとの単曲線一致率IIW1nを演算するためには、直前の通過曲線情報R2の各実測値(実測曲率半径R2az、実測入口緩和曲線長R2bz、実測円曲線長R2cz、実測出口緩和曲線長R2dz、実測曲線方向R2ez)と、候補曲線情報KW1nの各データ値(データ曲率半径W1a−n、データ入口緩和曲線長W1b−n、データ円曲線長W1c−n、データ出口緩和曲線長W1d−n、データ曲線方向W1e−n)とによって、各部分一致率Hが演算される。そして、これら各部分一致率Hが全て乗算されることで、単曲線一致率IIW1n(図17参照)が演算される。
以下、同様に、候補曲線情報のKE1n,KW2s,KE2sに対して、単曲線一致率IIE1n,IIW2s,IIE2sが演算される。ここで、図17は、本実施形態において演算された単曲線一致率IIW1n,IIE1n,IIW2s,IIE2sの値を示した図である。図17に示すように、例えば、単曲線一致率IIE1nは0.95であるため、直前より前の通過曲線情報R2と候補曲線情報KE1nとが非常に似ていて、直前より前に通過した曲線部が順行する場合の曲線部E1に非常に似ていることを意味している。また、例えば、単曲線一致率IIW2s,IIE2sは0であるため、直前より前の通過曲線情報R2と候補曲線情報KW2s,KE2sとが全く似ていないことを意味しているが、これは曲線方向が逆であるためである。
総合一致率演算処理の説明に戻る。本実施形態では、ステップ320において、一致率演算部60は、通過した2つの曲線部間の長さの一致度合いを意味する曲線間長一致率Lを演算するようになっている。この曲線間長一致率Lは、直前の通過曲線情報R1の実測曲線間長R1fz(図3参照)と、各候補曲線情報群Gのデータ曲線間長F1,F5,F9,F12,F13,F16とが、0から1までの範囲でどのくらい一致しているかを数値で表わしたものである。曲線間長一致率Lが1に近い程、2つの曲線部を通過したときに測定した実測距離S(実測曲線間長R1fz)と、候補として考える2つの曲線部間のデータ曲線間長Fとが一致していることを意味する。
例えば、候補曲線情報群G1の曲線間長一致率L1は、上記した数1を用いて、実測曲線間長R1fzを実測値に代入し且つデータ曲線間長F1をデータ値に代入することで、演算される。なお、上記した数1の部分一致率Hが、曲線間長一致率L1に代わる。以下、同様に、候補曲線情報群G5,G9,G12,G13,G16に対して、曲線間長一致率L5,L9,L12,L13,L16が演算される。ここで、図17には、本実施形態で演算された曲線間長一致率L1,L5,L9,L12,L13,L16の各値が示されている。例えば、曲線間長一致率L12は0.3であるため、実測曲線間長R1fzが候補として考える2つの曲線部の間のデータ曲線間長F12とほとんど一致していないことを意味している。また、例えば、曲線間長一致率L13は0.95であるため、実測曲線間長R1fzが候補として考える2つの曲線部の間のデータ曲線間長F13とほぼ一致していることを意味している。
こうして、単曲線一致率II及び曲線間長一致率Lを演算した後、ステップ330において、一致率演算部60は、各候補曲線情報群G1,G5,G9,G12,G13,G16毎に、総合一致率Jを演算する。総合一致率Jとは、直前の通過曲線情報R1と直前より前の通過曲線情報R2と実測曲線間長R1fzが、各候補曲線情報群Gに対して、0から1までの範囲でどのくらい一致しているかを数値で表わしたものである。総合一致率Jが1に近い程、複数の曲線部を通過したときに測定した各実測値の集まりが、候補曲線情報群Gに含まれる各データ値の集まりと一致していることを意味する。この総合一致率Jは、単曲線一致率IIと単曲線一致率Iと曲線間長一致率Lとが乗算されることによって、演算される。
例えば、図17に示すように、候補曲線情報群G1の総合一致率J1は、単曲線一致率IW1nと単曲線一致率IIW1nと曲線間長一致率L1とが乗算されることによって、演算される。以下、同様に、候補曲線情報群G5,G9,G12,G13,G16に対して、総合一致率J5,J9,J12,J13,J16が演算される。ここで、図17には、本実施形態で演算された総合一致率J1,J5,J9,J12,J13,J16の各値が示されている。本実施形態では、候補曲線情報群G13の総合一致率J13は0.86であり、他の総合一致率J1,J5,J9,J12,J13,J16よりも大きく且つ1に近い値である。こうして、総合一致率Jが演算された後、総合一致率演算処理が終了し、図9に示すステップ510に進む。
図9に示すように、ステップ510において、曲線特定部70が、各総合一致率J1,J5,J9,J12,J13,J16が、基準値Ja以上であるか否かを判断して、基準値Ja以上の総合一致率Jのみを残す。基準値Jaは、直前の通過曲線情報R1と直前より前の通過曲線情報R2と実測曲線間長R1fzとが、候補曲線情報群Gに対して一致していることを判断するための基準になる値であり、1に近い値に設定されている。本実施形態では、基準値Jaが0.85に設定されているが、適宜変更可能である。このため、ステップ510において、総合一致率J13のみが0.85以上であるため、曲線特定部70は、総合一致率J13のみを残す。
次に、ステップ520において、残った総合一致率J13が有るか否かを判断する。本実施形態では、総合一致率J13のみが残っているため、「Yes」と判断されて、ステップ530に進む。続いて、ステップ530において、曲線特定部70は、残った総合一致率Jが1つであるか否かを判断する。本実施形態では、総合一致率J13のみが残っているため、「Yes」と判断されて、ステップ540に進む。ステップ540において、曲線特定部70は、直前の通過曲線情報R1と直前より前の通過曲線情報R2と実測曲線間長R1fzとが、候補曲線情報群G13に一致していることを最終的に決定して、直前に通過した曲線部を特定する。こうして、鉄道車両1は自身が曲線部E1を順行して通過した後、直前に曲線部W2を順行して通過したことを特定することになる。
その後、ステップ550において、実測距離補正部80は、上述したように特定された曲線部W2に基づいて、実測距離Sを補正する。即ち、鉄道車両1は、直前に曲線部W2を北側に向かって走行したことが分かると共に、その曲線部W2の緩和曲線終了位置p4(図4参照)が分かるため、実測距離補正部80が現地点での実測距離Sを、曲線部W2の緩和曲線終了位置p4から測定された実測距離Sを用いて補正する。こうして、鉄道車両1は、直前に通過した曲線部を特定することで、自身の走行位置を検出することができる。ステップ550の後、ステップ560に進み、走行位置検出システムKSが終了する。なお、実測距離Sの補正は、曲線部を通過する毎に、上述した単曲線一致率演算処理(ステップ100)、候補曲線情報群作成処理(ステップ200)、総合一致率演算処理(ステップ300)が実行されることで行われる。つまり、実測距離補正部80が実測距離Sを補正するタイミングは、鉄道車両1が新たに曲線部を通過するときである。
走行位置が検出された後には、実測距離補正部80が現地点での走行位置を車体傾斜制御部90に送信し、データベース処理部30が次に通過する曲線部の曲線情報(カント、曲率半径、曲線方向、変曲点の位置等)を車体傾斜制御部90に送信する。これにより、車体傾斜制御部90は、送信された現地点での走行位置及び曲線情報に基づいて、次に通過する曲線部で車体を傾ける適切な地点を決定して、車体傾斜制御を適切に実行する。なお、実測距離Sは、実測距離補正部80から車体傾斜制御部90へ常時送られていて、車体傾斜制御部90は常時実行されている。
ここで、図9に示す走行位置検出システムKSのステップ570、及びステップ400のサブルーチンについて説明しておく。本実施形態と異なる場合として、仮に、各総合一致率J1,J5,J9,J12,J13,J16が全て基準値Jaより小さい結果が得られた場合、ステップ520において、残った総合一致率Jが無いため、「No」と判断されて、ステップ570に進むことになる。これは、直前に通過した曲線部及び直前より前に通過した曲線部が、候補の中で何れにも似ていないという結果が得られた場合の処理である。
上記した場合では、ステップ570において、曲線特定部70が直前に通過した曲線部を特定できないため、鉄道車両1が次の曲線部を通過するまで走行位置検出システムKSの実行を待機する。これは、直前に通過した曲線部を誤って特定することを回避するためである。従って、ステップ570に進んだ段階では、次の曲線部を通過するまで待機して、鉄道車両1が新たに曲線部を通過した段階でステップ100に戻る。そして、新たに用意する直前の通過曲線情報R1、直前より前の通過曲線情報R2、実測曲線間長R1fzに基づいて、走行位置検出システムKSの処理が実行されるようになっている。
また、本実施形態と異なる場合として、仮に、総合一致率J9及び総合一致率J13が基準値Jaより大きい結果が得られた場合(以下、「仮定の場合」と呼ぶ)、ステップ530において、残った総合一致率Jが2つであるため、「No」と判断されて、ステップ400のサブルーチンに進み、新総合一致率演算処理が実行されるようになっている(ステップ400)。これは、直前に通過した曲線部及び直前より前に通過した曲線部が、候補の中で似ているものが2つ以上存在するという結果が得られた場合の処理である。
このため、図18に示すように、新総合一致率演算処理が実行されると、先ず、ステップ410において、通過曲線情報作成部40が、それまでに用意していた通過曲線情報Rより1つ前の通過曲線情報Rを新たに用意する。このため、仮定の場合では、直前より前の前の通過曲線情報R3が新たに用意される。なお、以下の説明では、基本的に上述した処理(ステップ200及びステップ300)の繰り返しであるため、簡略して説明する。
次に、ステップ420において、候補曲線情報群作成部50Aは、基準値Jaより大きい総合一致率Jの原因になった各候補曲線情報群Gに対して、新たに用意した通過曲線情報Rに対応する候補曲線情報Kを含めて、新しい各候補曲線情報群Gnを作成する。このため、仮定の場合では、候補曲線情報群G9及び候補曲線情報群G13に対して、直前より前の前に通過したと考えられる曲線部の候補曲線情報Kを含めて、新しい各候補曲線情報群Gn9,Gn13,・・・を作成する。
続いて、ステップ430において、一致率演算部60は、新たに用意した通過曲線情報Rと、新しい各候補曲線情報群Gnの中で対応する候補曲線情報Kとの追加単曲線一致率Inを演算する。このため、仮定の場合では、直前より前の前の通過曲線情報R3の各実測値と、新しい各候補曲線情報群Gn9,Gn13,・・・の中で対応する候補曲線情報Kの各データ値とを、上記した数1に代入して、各部分一致率Hを算出して、追加単曲線一致率Inを演算する。
また、ステップ440において、一致率演算部60は、新たに通過曲線情報Rを用意したことよって考えられる実測曲線間長fnと新たに候補曲線情報Kを含めたことによって考えられるデータ曲線間長Fnとの追加曲線間長一致率Lnを演算する。このため、仮定の場合では、直前より前の前に通過した曲線部から直前より前に通過した曲線部までの実測距離S(実測曲線間長R2fz)と、新たに候補曲線情報Kを含めたことによって考えられるデータ曲線間長Fnとを、上記した数1に代入して、追加曲線間長一致率Lnを演算する。
最後に、ステップ450において、一致率演算部60は、新しい候補曲線情報群Gn毎に、それまでの総合一致率Jに対して、追加単曲線一致率In及び追加曲線間長一致率Lnを乗算して、新たに総合一致率Jを演算する。こうして、新総合一致率演算処理では、より多くの実測値及びデータ値を用いて、通過した複数の曲線部の一致度合いとして新たに総合一致率Jを演算するようになっている。このため、仮定の場合のように、総合一致率J9及び総合一致率J13が基準値Jaより大きい結果が得られた場合であっても、新総合一致率演算処理を実行することで、正しい結果を示す総合一致率Jのみが大きな値になるように演算することができる。
そして、新総合一致率演算処理が終わると、図9に示すステップ510に戻り、新たに演算された総合一致率Jが基準値Ja以上であるか否かが判断されて、基準値Ja以上の総合一致率Jのみが残される。以下同様に、ステップ520及びステップ530に進み、残った総合一致率Jが1つであれば、ステップ540において、直前に通過した曲線部を特定することができる。一方、仮に未だ基準値Ja以上の総合一致率Jが2つ以上存在する場合には、ステップ400のサブルーチンに進み、同様に新総合一致率演算処理が実行されることになる。従って、基準値Ja以上の総合一致率Jが2つ以上存在しても、新総合一致率演算処理を繰り返し実行することで、基準値Ja以上の総合一致率Jを1つだけ残すことができ、直前に通過した曲線部を特定することができる。
本実施形態の作用効果について説明する。
本実施形態によれば、一つの曲線部を通過したときに得られる通過曲線情報R1に基づいて一つの曲線部の一致度合い(単曲線一致率I)を比較するのではなく、複数の曲線部を通過したときに得られる複数の通過曲線情報R1,R2に基づいて複数の曲線部の集まりの一致度合い(総合一致率J)を比較している。即ち、一つずつの曲線部の一致度合いを単曲線一致率I,IIとして数値化し、更に曲線部間の曲線間長の一致度合いを曲線間長一致率Lとして数値化して、各単曲線一致率I,IIと曲線間長一致率Lとを乗算して総合一致率Jを演算している。このため、実際に直前に通過した曲線部W2とそれ以外の曲線部E2との形状とが似ていても、正しい結果を示す総合一致率J13のみを大きな値で表わすことができる。従って、直前に通過した曲線部がどの曲線部なのかが的確に特定できて、鉄道車両は自身の走行位置を検出することができる。
また、本実施形態によれば、上述したように、総合一致率Jを演算するためのパラメータとして、一つずつの曲線部の一致度合いだけでなく、曲線間長の一致度合いも含めることで、より正確な総合一致率Jを演算している。更に、本実施形態では、通過曲線情報Rの実測値として、実測曲率半径azと、実測入口緩和曲線長bzと、実測円曲線長czと、実測出口緩和曲線長dzとを用いている。これは、上記した各実測値が、その他の測定値に比べて、走行中に真値に近い正確な値を測定できるためである。従って、各部分一致率Hを算出するための実測値として相応しく、より正確な単曲線一致率I,IIを演算することができる。
また、本実施形態によれば、通過曲線情報Rの実測値として、実測曲線方向ezが含まれ、候補曲線情報のデータ値として、データ曲線方向eが含まれている。こうして、曲線部を特徴づけるパラメータに曲線方向が含まれているため、曲線部を特定する際に、その曲線部の曲がっている方向が分かる。即ち、一つの曲線部であっても、例えば順行するときには右向きに曲がる曲線部であり、逆行するときには左向きに曲がる曲線部になるため、曲線部の曲がっている方向が分かることで、順行しているのか又は逆行しているのかが分かることになる。
また、本実施形態によれば、候補曲線情報群作成部50Aは、通過する2つの曲線部の全ての組み合わせに対応する各候補曲線情報群G1〜G16(図14参照)を、単曲線一致率Iを演算した結果に基づいて、各候補曲線情報群G1,G5,G9,G12,G13,G16に絞り込む。これにより、候補曲線情報群G毎に演算する総合一致率Jの数が少なくなる。この結果、総合一致率Jを演算するための計算量を少なくすることができて、計算負荷を軽減することができる。
ところで、実測値(通過曲線情報R)とデータ値(候補曲線情報K)とによって、曲線部の一致度合いを0から1までの値で数値化する方法として、本実施形態のように上記した数1を用いて部分一致率H(単曲線一致率I)を演算する方法の他、相関関数を用いて相関演算する方法がある。そこで、相関演算する方法と、上記した数1を用いて演算する方法とを比較して説明する。なお、相関演算する方法の説明では、実測曲率とデータ曲率との一致度合いを求める場合を例にする。
先ず、図19は、相関演算する方法を説明するための図であり、或る曲線部のデータ曲率が波形状に実線で示され、或る曲線部を通過したときに測定された実測曲率が波形状に破線で示されている。図19に示すように、相関演算する方法では、データ曲率である実線と実測曲率である破線との一致度合いを求めるために、実線と破線の差である斜線で示した面積が最小になる地点を探すことになる。言い換えると、斜線で示した面積が最小になる地点を探すために、データ曲率と実測曲率を逐次読み込んで相関関数を用いて相関演算している。
このため、相関演算する方法では、計算量が非常に多くて、計算負荷が大きい。また、相関演算を行うためのデータ曲率として、波形を表わす非常に多くのデータの集合体を、予め走行試験で準備しておく必要がある。そして、準備しておく波形状のデータ曲率及び走行中に測定される波形状の実測曲率は、鉄道車両や線路の状況に応じて変化し易いため、正確な結果を得るために準備しなければならないデータ曲率は、何度か走行試験で取得した情報を加工して、常に更新していく仕組みが必要になる。従って、事前にデータ曲率を準備するための手間が大きい。
一方、図20は、上記した数1を用いて演算する方法を説明するための図であり、或る曲線部のデータ曲率が波形状に実線で示され、或る曲線部を通過したときに測定された実測曲率が破線で示されている。上記した数1を用いて演算する方法は、概念的には以下のように説明することができる。先ず、図20に示すように、線路の曲線部は入口緩和曲線と円曲線と出口緩和曲線とで構成されることに基づいて、実測曲率である破線を、入口緩和曲線と円曲線と出口緩和曲線とで直線状に区分けされた二点鎖線で理想化している。破線を直線状の二点鎖線で理想化できると考えたのは、円曲線での曲率の変化幅Z1は、曲線部全体での曲率の変化幅Zに比べて十分小さく、実測曲率である破線と理想化された二点鎖線との間の誤差は無視できるほど小さいためである。
こうして、上記した数1を用いて演算する方法では、実測曲率を理想化して、実線と二点鎖線との一致度合いを求めるようになっている。このとき、単曲線一致率Iを演算するために、上記した数1によって、曲率半径、入口緩和曲線長、円曲線長、出口緩和曲線長、曲線方向のそれぞれ部分一致率Ha,Hb,Hc,Hd,Heを算出し、且つそれぞれの部分一致率Ha,Hb,Hc,Hd,Heを乗算するだけであるため、計算量が少なくて、計算負荷が小さい。
そして、実線と一点鎖線を比較する方法が、曲線部のパラメータの各部分一致率Hを乗算する方法であるため、データ曲率(データ曲率半径)等の一つのパラメータにおいて、鉄道車両や線路の状況に応じて変化し易い値であっても、乗算の集合である単曲線一致率Iにはほとんど影響はない。従って、相関演算する方法のように、正確な結果を得るために何度か走行試験で取得した情報を加工して、常にデータ曲率を更新する必要はない。即ち、車体傾斜制御するために各曲線部に対して予め記憶しているデータ曲線情報(図5参照)だけを利用して実施することができ、事前にデータを準備する手間をほとんど無くすことができる。
また、上述したように、単曲線一致率Iは任意に選択したパラメータの各部分一致率を乗算するだけで演算されるものであるため、曲線部を特徴づけるために相応しいパラメータを新たに追加したい場合には、追加したパラメータの部分一致率を算出して乗算するだけで、対応することができる。従って、上記した数1を用いて単曲線一致率Iを演算する本実施形態の方法は、柔軟性が高い方法である。
以上、本発明に係る鉄道車両の走行位置検出システムの実施形態について説明したが、本発明はこの実施形態に限定されることはなく、その趣旨を逸脱しない範囲で様々な変更が可能である。
例えば、本実施形態では、部分一致率Hを求めるために上記した数1を用いたが、上記した数1に換えて、変形実施形態として下記の数2を用いても良い。
Figure 2016141300
本実施形態のように、上記した数1を用いて部分一致率Hを演算する場合には、各部分一致率Hを算出する際の分母が全てデータ値で統一されるため、各部分一致率Hを同じ基準で数値化することができる。一方、上記した数2の各部分一致率Hは、分母についてデータ値又は実測値のどちらかを用いて算出される。このため、上記した数1の場合には、分母を全てデータ値に統一することで、単曲線一致率Iをデータ値に対して実測値がどの程度類似しているかを正規化した数値で評価できる。一方、上記した数2の場合には、上記した数1を用いる場合のように、実測値がデータ値2倍より大きい値であるときに部分一致率Hを「0」に定義するという制約を無くすことができ、実測値の全ての値を対象にすることができる。
また、本実施形態では、通過曲線情報Rの実測値として実測曲率半径azを用い、候補曲線情報Kのデータ値としてデータ曲率半径aを用いたが、実測曲率半径azに換えて実測曲率を用い、データ曲率半径aに換えてデータ曲率を用いても良い。また、本実施形態では、通過曲線情報Rの各実測値、及び候補曲線情報Kの各データ値のパラメータとして、曲率半径、入口緩和曲線長、円曲線長、出口緩和曲線長を用いたが、パラメータは上記したものに限定されるものではない。即ち、本発明は、曲線部を特徴づける任意のパラメータについてデータ値(理想値)と実測値との一致度合いを照合することを特徴とするものであり、パラメータは主成分分析等によって任意に選択することができる。
また、本実施形態では、図9のステップ530及びステップ540に示すように、残った総合一致率Jが1つのみである場合に、直前に通過した曲線部を特定するように設定したが、直前に通過した曲線部を特定する方法は、適宜変更可能である。例えば、残った総合一致率Jが任意の指定数個である場合に、運転士が残った総合一致率Jの値を視覚的に認識することができ、その総合一致率Jの値に応じて直前に通過した曲線部を特定しても良い。
また、本実施形態では、分かり易く説明するために、図3に示した単線並列による双方向運用を例にしたが、本実施形態の走行位置検出システムKSは、あらゆる路線状況に対して適用することができる。従って、例えば、駅間が一本の線路でつながっている路線状況、又は駅間が三本以上の線路でつながっている路線状況に対して適用することができる。また、鉄道車両が線路に対して双方向に走行する路線状況だけでなく、線路に対して一方向にのみ走行する路線状況でも当然適用できる。この場合には、曲線部を特徴づけるパラメータに曲線方向を含めなくても良い。
次に、S5の別の実施形態について説明する。内容的には、前述した実施の形態と同じ部分が多いので、相違する点のみ詳細に説明する。図1において、候補曲線情報群作成部50Aを備えておらず、候補曲線情報作成部50は、直接一致率演算部60に接続している。 また、図10において、S170がなく、S160から直接S190に繋がっている。
図23は、単曲線一致率Iを演算する本実施形態の方法を説明するための図であり、或る曲線部のデータ曲率が波形状に実線で示され、或る曲線部を通過したときに測定された実測曲率が破線で示されている。単曲線一致率Iを演算する本実施形態の方法は、概念的には以下のように説明することができる。
先ず、図23に示すように、線路の曲線部は入口緩和曲線と円曲線と出口緩和曲線とで構成されることに基づいて、実測曲率である破線を、入口緩和曲線と円曲線と出口緩和曲線とで直線状に区分けされた二点鎖線で理想化している。破線を直線状の二点鎖線で理想化できると考えたのは、円曲線での曲率の変化幅Z1は、曲線部全体での曲率の変化幅Zに比べて十分小さく、実測曲率である破線と理想化された二点鎖線との間の誤差は無視できるほど小さいためである。
こうして、単曲線一致率Iを演算する本実施形態の方法では、実測曲率を理想化して、実線と二点鎖線との一致度合いを求めるようになっている。このとき、単曲線一致率Iを演算するためには、例えば、数式1によって、曲率半径、入口緩和曲線長、円曲線長、出口緩和曲線長、曲線方向のそれぞれ部分一致率Ha,Hb,Hc,Hd,Heを算出し、且つそれぞれの部分一致率Ha,Hb,Hc,Hd,Heを乗算するだけであるため、計算量が少なくて、計算負荷が小さい。
そして、実線と二点鎖線を比較する方法が、曲線部のパラメータの各部分一致率Hを乗算する方法であるため、データ曲率(データ曲率半径)等の一つのパラメータにおいて、鉄道車両や線路の状況に応じて変化し易い値であっても、乗算の集合である単曲線一致率Iにはほとんど影響はない。従って、相関演算する特許文献1の方法のように、正確な結果を得るために何度か走行試験で取得した情報を加工して、常にデータ曲率を更新する必要はない。即ち、車体傾斜制御するために各曲線部に対して予め記憶しているデータ曲線情報だけを利用して実施することができ、事前にデータを準備する手間をほとんど無くすことができる。
1 鉄道車両
2 車上子
3 速度発電機
4 ヨー角速度センサ
10 実測距離演算部
20 走行速度演算部
30 データベース処理部
40 通過曲線情報作成部
50 候補曲線情報作成部
50A 候補曲線情報群作成部
60 一致率演算部
70 曲線特定部
80 実測距離補正部
90 車体傾斜制御部
KS 走行位置検出システム
D データ曲線情報
R 通過曲線情報
K 候補曲線情報
G 候補曲線情報群
H 部分一致率
I,II 単曲線一致率
In 追加単曲線間一致率
L 曲線間長一致率
Ln 追加曲線間長一致率
J 総合一致率
a,az データ曲率半径,実測曲率半径
b,bz データ入口緩和曲線長,実測入口緩和曲線長
c,cz データ円曲線長,実測円曲線長
d,dz データ出口緩和曲線長,実測出口緩和曲線長
e,ez データ曲線方向,実測曲線方向
f,fz データ曲線間長,実測曲線間長

Claims (13)

  1. 鉄道車両が、複数路線網から選択された複数の特定路線区間を含む特定路線を走行すること、前記鉄道車両が自車の走行位置を検出すること、を備える鉄道車両の走行位置検出システムにおいて、
    前記鉄道車両は、線路近くに設置された識別子の識別子情報を読み込んで、前記特定路線のうちの、前記鉄道車両が走行している前記特定路線区間を決定する特定路線決定手段と、
    前記複数の路線網を前記鉄道車両が走行したときの実測走行データを記憶する走行データ記憶手段と、
    前記特定路線決定手段が決定した前記特定路線区間を含む前記特定路線内の前記実測走行データを、前記走行データ記憶手段から読み出し、現在走行中に検知している現在走行データと比較することにより、自車の走行位置を検出すること、
    を特徴とする鉄道車両の走行位置検出システム。
  2. 請求項1に記載する鉄道車両の走行位置検出システムにおいて、
    鉄道車両が曲線部を通過したときに測定した実測曲線情報と予め記憶している多数の曲線部のデータ曲線情報とを比較して、どの曲線部を通過したかを特定し、特定した曲線部からの実測距離に基づいて自車の走行位置を検出すること、
    を特徴とする鉄道車両の走行位置検出システム。
  3. 請求項2に記載する鉄道車両の走行位置検出システムにおいて、
    曲線部を通過する毎に、測定した複数の実測値を含む通過曲線情報を作成する通過曲線情報作成部と、
    通過し得る全ての曲線部に対応していて、且つ前記データ曲線情報に基づいて前記実測値に対応するデータ値を含む候補曲線情報を作成する候補曲線情報作成部と、
    直前の通過曲線情報及び直前より前の通過曲線情報にそれぞれ対応する候補曲線情報の集まりを候補曲線情報群として複数個作成する候補曲線情報群作成部と、
    直前の通過曲線情報及び直前より前の通過曲線情報と各候補曲線情報群との一致度合いを求めるために、直前の通過曲線情報と候補曲線情報群の中で対応する候補曲線情報との一致度合いについて0から1までの範囲で大きい値ほど一致していることを表す単曲線一致率を演算すると共に、直前より前の通過曲線情報と候補曲線情報群の中で対応する候補曲線情報との一致度合いについて0から1までの範囲で大きい値ほど一致していることを表す単曲線一致率を演算して、各単曲線一致率を乗算して総合一致率を演算する一致率演算部と、
    前記候補曲線情報群毎にそれぞれ演算された総合一致率を比較して直前に通過した曲線部を特定する曲線特定部とを備えていること、
    を特徴とする鉄道車両の走行位置検出システム。
  4. 請求項3に記載された鉄道車両の走行位置検出システムにおいて、
    前記通過曲線情報作成部は、前記通過曲線情報の実測値として、通過する2つの曲線部間の曲線間長である実測曲線間長を含め、
    前記候補曲線情報作成部は、前記候補曲線情報のデータ値として、前記データ曲線情報に基づいて2つの曲線部間の曲線間長であるデータ曲線間長を含め、
    前記一致率演算部は、実測曲線間長とデータ曲線間長との一致度合いについて0から1までの範囲で大きい値ほど一致していることを表わす曲線間長一致率を演算し、且つ前記曲線間長一致率と前記各単曲線一致率とを乗算して総合一致率を演算すること、
    を特徴とする鉄道車両の走行位置検出システム。
  5. 請求項3又は請求項4に記載された鉄道車両の走行位置検出システムにおいて、
    前記一致率演算部は、前記総合一致率を演算する前に、直前の通過曲線情報と各候補曲線情報との一致度合いについて0から1までの範囲で大きい値ほど一致していることを表わす単曲線一致率を演算して、
    前記候補曲線情報群作成部は、前記演算された単曲線一致率のうち所定値より大きい単曲線一致率が2つ以上存在する場合、各候補曲線情報群について前記所定値より大きい単曲線一致率が得られる原因になった候補曲線情報を含むように絞り込むこと、
    を特徴とする鉄道車両の走行位置検出システム。
  6. 請求項3乃至請求項5の何れか1つに記載する鉄道車両の走行位置検出システムにおいて、
    前記候補曲線情報群毎に演算された総合一致率のうち基準値より大きい総合一致率が2つ以上存在する場合、
    前記通過曲線情報作成部は、それまでに用意していた通過曲線情報より1つ前の通過曲線情報を新たに用意し、
    前記候補曲線情報群作成部は、基準値より大きい総合一致率が得られる原因になった各候補曲線情報群に対して、前記新たに用意した通過曲線情報に対応する候補曲線情報を含めて、新しい各候補曲線情報群を作成し、
    前記一致率演算部は、前記新たに用意した通過曲線情報と前記新しい候補曲線情報群の中で対応する候補曲線情報との一致度合いについて0から1までの範囲で大きい値ほど一致していることを表す追加単曲線一致率を演算して、この追加単曲線一致率を前記総合一致率に乗算して新たに総合一致率を演算すること、
    を特徴とする鉄道車両の走行位置検出システム。
  7. 請求項3乃至請求項6の何れか1つに記載する鉄道車両の走行位置検出システムにおいて、
    前記一致率演算部は、前記実測値と前記データ値との差から前記実測値又は前記データ値を除算した値を1から減算して部分一致率を算出し、前記実測値毎に算出した各部分一致率を乗算して、前記単曲線一致率を演算すること、
    を特徴とする鉄道車両の走行位置検出システム。
  8. 請求項7に記載された鉄道車両の走行位置検出システムにおいて、
    個々の曲線部をパラメータの集合ととらえ、前記パラメータについて、0から1までの範囲で一致度合いを算出し、その乗算で前記単曲線一致率を演算することで個々の曲線部を評価すること、
    を特徴とする鉄道車両の走行位置検出システム。
  9. 請求項7又は請求項8に記載された鉄道車両の走行位置検出システムにおいて、
    前記通過曲線情報作成部は、前記通過曲線情報の実測値として、実測曲率半径又は実測曲率の実測値を含め、
    前記候補曲線情報作成部は、前記候補曲線情報のデータ値として、データ曲率半径又はデータ曲率のデータ値を含めること、
    を特徴とする鉄道車両の走行位置検出システム。
  10. 請求項9に記載された鉄道車両の走行位置検出システムにおいて、
    前記通過曲線情報作成部は、前記通過曲線情報の実測値として、実測入口緩和曲線長と、実測円曲線長と、実測出口緩和曲線長とを含め、
    前記候補曲線情報作成部は、前記候補曲線情報のデータ値として、データ入口緩和曲線長と、データ円曲線長と、データ出口緩和曲線長とを含めること、
    を特徴とする鉄道車両の走行位置検出システム。
  11. 請求項7乃至請求項10の何れか1つに記載する鉄道車両の走行位置検出システムにおいて、
    前記通過曲線情報作成部は、前記通過曲線情報の実測値として、実測曲線方向を含め、
    前記候補曲線情報作成部は、前記候補曲線情報のデータ値として、データ曲線方向を含めること、
    を特徴とする鉄道車両の走行位置検出システム。
  12. 請求項2に記載する鉄道車両の走行位置検出システムにおいて、
    曲線部を通過する毎に、測定した複数の実測値を含む通過曲線情報を作成する通過曲線情報作成部と、
    通過し得る全ての曲線部に対応していて、且つ前記データ曲線情報に基づいて前記実測値に対応するデータ値を含む候補曲線情報を作成する候補曲線情報作成部と、
    直前の通過曲線情報と各候補曲線情報との一致度合いを求めるために、前記実測値と前記データ値との差から前記実測値又は前記データ値を除算した値を1から減算して部分一致率を算出し、実測値毎に算出した各部分一致率を乗算して単曲線一致率を演算する一致率演算部と、
    前記候補曲線情報毎にそれぞれ演算された単曲線一致率を比較して直前に通過した曲線部を特定する曲線特定部とを備えていること、
    を特徴とする鉄道車両の走行位置検出システム。
  13. 請求項12に記載された鉄道車両の走行位置検出システムにおいて、
    個々の曲線部をパラメータの集合ととらえ、前記パラメータについて、0から1までの範囲で一致度合いを算出し、その乗算で前記単曲線一致率を演算することで個々の曲線部を評価すること、
    を特徴とする鉄道車両の走行位置検出システム。
JP2015019431A 2015-02-03 2015-02-03 鉄道車両の走行位置検出システム Active JP6509573B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2015019431A JP6509573B2 (ja) 2015-02-03 2015-02-03 鉄道車両の走行位置検出システム

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2015019431A JP6509573B2 (ja) 2015-02-03 2015-02-03 鉄道車両の走行位置検出システム

Publications (2)

Publication Number Publication Date
JP2016141300A true JP2016141300A (ja) 2016-08-08
JP6509573B2 JP6509573B2 (ja) 2019-05-08

Family

ID=56568281

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2015019431A Active JP6509573B2 (ja) 2015-02-03 2015-02-03 鉄道車両の走行位置検出システム

Country Status (1)

Country Link
JP (1) JP6509573B2 (ja)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106541963A (zh) * 2016-10-12 2017-03-29 同济大学 基于特征点触发的获取轨道曲线线路信息的系统及方法
JP2018095146A (ja) * 2016-12-15 2018-06-21 日本車輌製造株式会社 鉄道車両
WO2019198633A1 (ja) * 2018-04-09 2019-10-17 日本製鉄株式会社 鉄道車両の現在位置情報補正方法
WO2020129423A1 (ja) * 2018-12-20 2020-06-25 株式会社日立製作所 位置検出装置及び方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000168552A (ja) * 1998-12-01 2000-06-20 Nippon Sharyo Seizo Kaisha Ltd 鉄道車両の走行位置検知方法
JP2007015483A (ja) * 2005-07-06 2007-01-25 Odakyu Dentetsu Kk 鉄道車両の走行位置検出方法
JP2008143390A (ja) * 2006-12-12 2008-06-26 Yoshihiro Suda 姿勢制御装置
JP2014097704A (ja) * 2012-11-13 2014-05-29 Nippon Sharyo Seizo Kaisha Ltd 鉄道車両の走行距離検出システム

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000168552A (ja) * 1998-12-01 2000-06-20 Nippon Sharyo Seizo Kaisha Ltd 鉄道車両の走行位置検知方法
JP2007015483A (ja) * 2005-07-06 2007-01-25 Odakyu Dentetsu Kk 鉄道車両の走行位置検出方法
JP2008143390A (ja) * 2006-12-12 2008-06-26 Yoshihiro Suda 姿勢制御装置
JP2014097704A (ja) * 2012-11-13 2014-05-29 Nippon Sharyo Seizo Kaisha Ltd 鉄道車両の走行距離検出システム

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106541963A (zh) * 2016-10-12 2017-03-29 同济大学 基于特征点触发的获取轨道曲线线路信息的系统及方法
CN106541963B (zh) * 2016-10-12 2019-04-16 同济大学 基于特征点触发的获取轨道曲线线路信息的系统及方法
JP2018095146A (ja) * 2016-12-15 2018-06-21 日本車輌製造株式会社 鉄道車両
WO2019198633A1 (ja) * 2018-04-09 2019-10-17 日本製鉄株式会社 鉄道車両の現在位置情報補正方法
JPWO2019198633A1 (ja) * 2018-04-09 2021-02-12 日本製鉄株式会社 鉄道車両の現在位置情報補正方法
WO2020129423A1 (ja) * 2018-12-20 2020-06-25 株式会社日立製作所 位置検出装置及び方法

Also Published As

Publication number Publication date
JP6509573B2 (ja) 2019-05-08

Similar Documents

Publication Publication Date Title
JP6509573B2 (ja) 鉄道車両の走行位置検出システム
CN107403482B (zh) 一种确定车辆行驶里程数的方法、装置及系统
CN106104656B (zh) 地图信息生成系统、方法及程序
CN102481952B (zh) 用于稳健地识别车辆的车道变换的方法和控制设备
US7756634B2 (en) Navigation device and current-position-calculation method
US9714836B2 (en) Vehicle position identifying system and vehicle position identifying method
US11603114B2 (en) Apparatus for generation of data of travel path inside intersection, program for generation of data of travel path inside intersection, and storage medium
CN108107883B (zh) 一种基于磁带导引agv的多传感器信息融合定位方法
CN111656420B (zh) 交叉点内的行驶轨道数据生成装置、以及存储介质
CN105928529A (zh) 一种多证据融合地图匹配算法
JP5730068B2 (ja) 鉄道車両の位置検出装置および位置検出方法
JP6203607B2 (ja) 鉄道車両の走行位置検出システム
CN111856521A (zh) 数据处理方法、装置、电子设备和存储介质
CN111024079B (zh) 一种根据多个位置点与路线进行匹配的方法和系统
KR20130057893A (ko) 클로소이드 곡선을 이용한 주행경로 생성방법 및 이를 수행하기 위한 프로그램이 기록된 컴퓨터로 읽을 수 있는 기록매체
MX2022002423A (es) Metodo y vehiculo de medicion para determinar una posicion real de una via.
CN104422449A (zh) 一种车辆导航的方法、装置
CN106796126A (zh) 用于实时地提供动态测量值的动态误差值的方法和系统
JP2008134187A (ja) 車載用ナビゲーション装置
WO2013111421A1 (ja) 勾配情報取得方法、記憶媒体、勾配情報取得装置およびプログラム
JP2007015483A (ja) 鉄道車両の走行位置検出方法
CN112750219A (zh) 车辆出行费用计算方法及系统
JP2015093605A (ja) 鉄道車両の走行位置検出システム
KR100862721B1 (ko) 도시철도차량용 차륜경 측정 시스템 및 그 방법
CN110440826B (zh) 一种位置姿态测量系统里程计数据处理方法

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20180115

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20180822

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20180828

A601 Written request for extension of time

Free format text: JAPANESE INTERMEDIATE CODE: A601

Effective date: 20181001

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20181119

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: 20190402

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20190403

R150 Certificate of patent or registration of utility model

Ref document number: 6509573

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250