位置推定装置、位置推定方法及びプログラム
【課題】高い精度で移動局の位置を推定する。
【解決手段】位置が既知である複数の基地局２０と位置推定対象である移動局１０との一方向又は双方向通信により取得した所定の信号に基づき、複数の基地局２０と移動局１０との観測距離を計測する距離計測部３２と、距離計測部３２により計測された複数の基地局２０の観測距離に基づき、予め定められた移動局の推定初期位置４０の方向に各基地局２０の観測座標を決定する観測座標決定部３４と、観測座標決定部３４により決定された複数の基地局２０の観測座標から移動局１０の代表位置を移動局１０の推定位置５０として算出する推定位置算出部３６と、を有する位置推定装置を提供する。
【解決手段】位置が既知である複数の基地局２０と位置推定対象である移動局１０との一方向又は双方向通信により取得した所定の信号に基づき、複数の基地局２０と移動局１０との観測距離を計測する距離計測部３２と、距離計測部３２により計測された複数の基地局２０の観測距離に基づき、予め定められた移動局の推定初期位置４０の方向に各基地局２０の観測座標を決定する観測座標決定部３４と、観測座標決定部３４により決定された複数の基地局２０の観測座標から移動局１０の代表位置を移動局１０の推定位置５０として算出する推定位置算出部３６と、を有する位置推定装置を提供する。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、移動局の位置を推定する位置推定装置、位置推定方法及びプログラムに関する。より詳細には、本発明は、移動局と複数の基地局で通信を行い、その通信結果である無線信号の受発信時刻や、基地局側での受信時刻に基づいて移動局の位置を推定する位置推定装置、位置推定方法及びプログラムに関する。
【背景技術】
【０００２】
従来から、移動局と複数の基地局で通信を行い、その通信結果である無線信号や受発信時刻などに基づいて移動局の位置を推定する測位システムが提案されている。たとえば、特許文献１〜３には、複数の無線基地局と移動局との間で無線信号を送受信し、移動局の位置を推定するシステムが提案されている。
【０００３】
無線信号としては、電波、音波、光、磁気などを用いたものがある。それらを用いて位置を求める方式のうち主要なものを大別すると、信号の受信強度を用いるＲＳＳＩ（Ｒｅｃｅｉｖｅｄ Ｓｉｇｎａｌ Ｓｔｒｅｎｇｔｈ Ｉｎｄｉｃａｔｏｒ）方式（特許文献１参照）、信号の伝播時間を用いるＴＯＡ（Ｔｉｍｅ Ｏｆ Ａｒｒｉｖａｌ）方式（特許文献２参照）、信号の伝播時間差を用いるＴＤＯＡ（Ｔｉｍｅ Ｄｉｆｆｅｒｅｎｃｅ Ｏｆ Ａｒｒｉｖａｌ）方式、信号の受信方位を用いるＡＯＡ（Ａｎｇｌｅ Ｏｆ Ａｒｒｉｖａｌ）方式（特許文献３参照）の４方式が挙げられる。
【０００４】
ＲＳＳＩ方式は、特許文献１に示したように、信号の受信強度から距離を推定し、そこから３辺測量で位置を推定する方式である。この方式は、高精度に距離を推定することが難しい一方で、厳密かつ高精度な時刻管理が不要で、様々な機器に造りこみやすいといった特徴がある。受信強度から距離を推定することが難しいため、複数の受信信号の強弱のパターンの学習と識別により位置を推定するものもあるが、これには付近に学習済のエリアがなければ位置を推定できないといった課題がある。
【０００５】
ＴＯＡ方式は、信号の伝播時間から距離を計測し、そこから３辺測量で位置を推定する方式である。例えば、移動局から基地局に向かって信号を送信した際の信号の伝播時間をｔとし（通信方向が逆の場合や、双方向通信のものもある）、信号の伝播速度をｃとすると、基地局と移動局間の距離ｌは、ｌ＝ｃ×ｔで求めることができる。無線信号の一例として電波を媒体としたものの場合、一般に空気中の電波の伝播速度は一定（約３０万ｋｍ／秒）のため、時計の精度が十分に高ければ安定した位置精度を実現しやすい。例えば移動局と基地局の全ての時計が同期している場合、伝播時間ｔの計測は１回の通信で行なえる。他にも同期なしで距離を求める方法もあるが、その場合には、移動局から発信した電波が基地局に到着し、更に基地局から移動局にて受信した電波を検出する必要があるため処理の負荷が高くなる。
【０００６】
例えばＴＯＡ方式で用いられる測距の一方式であるＴＷＲ（Ｔｗｏ Ｗａｙ Ｒａｎｇｉｎｇ）方式では、移動局から基地局に信号を送出し、信号を受信した基地局は直ちに移動局に返信を行ない、移動局は、自局が信号を送出してから基地局からの信号を受信するまでの往復伝播時間から、基地局で折り返しに必要な処理時間を減算して２で割ることにより、片道にかかった伝播時間を求めることができる。いずれにしても、信号の伝播時間を計測するため、高精度な内部時計は必要である。たとえば電波を用いた場合、約１ｎｓの時刻誤差で約３０ｃｍの距離誤差が生じる。
【０００７】
ＴＤＯＡ方式は、ＴＯＡ方式と同様に信号の伝播時間をベースとしているが、基地局と移動局の距離を計測せずに、移動局からマルチキャスト発信した信号を同時に複数の基地局にて受信し、その受信時間差に基づいて位置を求める。２つの基地局で受信した信号の受信時刻差は、信号の伝播速度を掛け合わせると、移動局から見た距離差となる。２つの基地局から同一の距離差を発生し得る位置は双曲線を描くため、基地局３局以上で同時に信号を受信すると、複数の双曲線が得られ、それらの交点を求めることにより位置が決定する。この方式は、ＴＯＡ方式と比較し、移動局が信号を受信する必要がない点、１回の測位に必要な信号の送信回数が少ない（通常、移動局からの１回の信号発信で位置が確定する）といった利点がある一方、基地局間で時刻を厳密に同期しなければならず、また基地局に内蔵する時計には必ず誤差があるため、時計の精度にあわせて逐次時刻同期を行なう必要がある。
【０００８】
ＡＯＡ方式は、複数の基地局が移動局から信号を受信した際の信号の到達角度に基づき、三角測量を行なう方式である。信号の到達角度は、例えば電波の場合は、アレイアンテナ、音波の場合はマイクロフォンアレイなどを用い、ＭＵＳＩＣ（ＭＵｌｔｉｐｌｅ ＳＩｇｎａｌ Ｃｌａｓｓｉｆｉｃａｔｉｏｎ）法などを用いることで推定できる。距離的に非常に近いアンテナ（マイクロフォン）間の受信時刻差を用いるため、非常に高い精度で個々の受信信号をサンプリングする必要があるといった難しさがある。
【０００９】
このように、各方式には様々な状況により得手、不得手があり、すべてのケースにおいて常に最良となる方式は存在しない。特に、これら無線信号を用いる方式の共通の課題として、反射波による影響がある。
【００１０】
例えば特許文献２では、ＴＯＡ方式で位置を推定するシステムにおいて、反射波の影響を低減する方式が提案されている。通常、基地局と移動局の間に見通し（ＬＯＳ： Ｌｉｎｅ ｏｆ ｓｉｇｈｔ）があり、直線で結ばれる経路を通る直接波で通信が行なわれる場合には、直線距離を求めることが可能である。しかしながら、現実には経路途中の構造物などにより反射が発生することがあり、このような反射波を用いて距離を計測すると、直線距離よりも長い距離が得られてしまう。これは、基地局と移動局を最短で結ぶ経路がお互いを直線で結ぶ直接波の通る経路であり、それ以外の反射波の通る経路はすべて直線よりも長くなるためである。しかし、この関係を逆に利用することで、なるべく直接波の経路に近いものを選ぶことができる。例えば１回の通信の中で直接波と反射波が混在する合成波の場合は、一番初めに到達した成分が直接波による成分であり、この成分により伝播時間を求めることができれば、直線距離を算出することができる。また例えば複数回、信号の伝播時間を計測でき、その中に反射波と直接波による計測値が混在している場合は、最小の計測値を選ぶことで直接波による計測結果を選ぶことができ、それに基づき直線距離を得ることができる。また、この最小値を選ぶという処理は、たとえ直接波による計測結果がなく、すべて反射波によるものであったとしても、直接波と比較して最も差の少ない反射波を選ぶことになり、距離計測の誤差を低減するような直線距離に近い値を選ぶことができる。通信の受信品質をあわせて用いると、多重反射波などの品質の悪い通信を排除でき、より直接波に近いものを選ぶことも可能かもしれない。長期間にわたって継続して計測を行なう場合、移動局の基地局から見た平均移動速度を算出することで、よりもっともらしい位置を算出することも可能である。
【００１１】
特許文献３では、ＴＤＯＡ方式とＡＯＡ方式を組み合わせることで、位置精度を高めようとする技術が開示されている。ＡＯＡ方式は、角度の分解能が一定の場合、基地局と移動局の距離が離れるに従い、位置の分解能が低下する。一方で基地局と距離が十分に近い場合は、高い位置分解能を有する。かかる特徴により、信号の受信時刻が早いものが、移動局に最も近い基地局であり、位置分解能が高いとし、高い重み付けを行なう。このように時刻により重みを変えて加重平均をとることで、よりもっともらしい位置の推定を行なう。また、特許文献３では、受信時刻を得るため、ＴＤＯＡ方式の測位もあわせて行い、それらを重み付けによって加重平均を取ることにより、もっともらしい位置を計測する方式が提案されている。
【先行技術文献】
【特許文献】
【００１２】
【特許文献１】特開２００９−８５７８０号公報
【特許文献２】特開２００１−２７５１４８号公報
【特許文献３】特開２００７−１３５００号公報
【発明の概要】
【発明が解決しようとする課題】
【００１３】
このように、従来の位置推定方法では、基地局と移動局の距離を、電波等を利用して基地局毎に観測し、複数の基地局から得られる観測距離を用いて移動局の位置を推定していた。
【００１４】
しかし、上記従来手法では、例えば基地局と移動局の間に見通しがなく（このような状態を以下、ＮＬＯＳ（Ｎｏｎ−Ｌｉｎｅ−Ｏｆ−Ｓｉｇｈｔ）と呼ぶ。）、直接波が届かずに反射波で伝播経路が構成されるとき等、直接波による測定が十分に行なえない場合や、劣悪な伝播環境等の影響で異常な観測値を出力する基地局がある場合に、位置推定の精度が大きく劣化していた。
【００１５】
例えば特許文献２の方法では、ＮＬＯＳの状態にあったとしても、時間的な経過により移動局の移動や環境の変化などで、ＮＬＯＳの状態が一時的にでも解消し、直接波が到達すると仮定している。例えば屋外で移動局が道路上を継続的に移動しており、ビルなどの建築物の影響で一時的にＮＬＯＳ状態となっている場合は、時間的経過によってＬＯＳ状態（Ｌｉｎｅ ｏｆ ｓｉｇｈｔ：見通しがある状態）とＮＬＯＳ状態が交互に遷移し、長期的にみれば直接波による距離計測が可能な場合も多くある。このような場合では、文献２でもある程度正しく位置を推定可能である。しかし、たとえばオフィス空間などの屋内環境で人や物体を対象とした位置検知を行う場合では、室内の壁や棚などの構造物によりＮＬＯＳ状態になると、検知対象が継続的に移動しない場合が多いため、なかなかＮＬＯＳ状態が解消しないといった課題がある。特許文献２の方法では、このような課題に対処できない。たとえば電波の場合、電波暗室でない屋内など周囲に電波を良く反射する壁面のある環境などでは、直接波と反射波が複雑に混在することになるが、このような場合に電波の強度や雑音などの問題で、直接波を取り逃がす場合がよく起こり得る。このような場合、直接波を受信できる頻度が問題となり、頻度が低い場合には、十分な精度で位置測位が出来ないといった問題も生じていた。
【００１６】
また、特許文献３の方法では、特にＡＯＡ方式の分解能が高い基地局近くに移動局がある場合に反射波を拾ってしまうと、移動局が遠方にある場合よりも信号の入射角が大きく変わり、位置精度を大きく低下させる原因となっていた。
【００１７】
上記問題に鑑み、本発明は、ＮＬＯＳの状態の場合や直接波と反射波が混在して受信するような環境下であっても、高い精度で移動局の位置を推定することを目的とする。
【課題を解決するための手段】
【００１８】
上記課題を解決するために、本発明のある観点によれば、位置が既知である複数の基地局と位置推定対象である移動局との一方向又は双方向通信により取得した所定の信号に基づき、前記複数の基地局と前記移動局との観測距離を計測し、該移動局の位置を推定する位置推定装置であって、計測された前記複数の基地局の観測距離に基づき、予め定められた移動局の推定初期位置の方向に各基地局の観測座標を決定する観測座標決定部と、前記観測座標決定部により決定された前記複数の基地局の観測座標から移動局の代表位置を前記移動局の推定位置として算出する推定位置算出部と、を備えることを特徴とする位置推定装置が提供される。
【００１９】
かかる構成によれば、移動局の推定位置を予め適当な推定初期位置(座標)に仮決めし、観測した複数の基地局の観測距離に基づき、仮決めした推定初期位置（座標）の方向に各基地局の観測座標を決定する。
【００２０】
たとえば、各基地局と前記移動局の推定初期位置とを結ぶ直線上であって、計測された観測距離だけ各基地局から離れた位置に各基地局の観測座標を算出する。
【００２１】
このようにして、個々の基地局の観測距離に基づき推定初期位置の誤差成分を表すために観測座標を基地局毎に算出する。推定位置算出部は、複数の基地局の観測座標の集合から移動局の代表位置を算出し、その値を移動局の推定位置とする。
【００２２】
前記推定位置算出部により算出された移動局の推定位置と前記移動局の推定初期位置との変化量が予め定められた停止条件を満たすまで、前記距離計測部、前記観測座標決定部及び前記推定位置算出部の動作を繰り返す繰返し判定部をさらに有していてもよい。
【００２３】
前記繰返し判定部は、前記推定位置算出部により算出された移動局の推定位置と前記移動局の推定初期位置との変化量が予め定められた停止条件を満たしたとき、前記移動局の推定位置を前記移動局の位置と判定してもよい。
【００２４】
このようにして移動局の推定位置の変化量が十分に小さくなれば、処理が収束したとして最後に算出された移動局の推定位置を移動局の位置と判定する。これにより、異常な基地局が存在する場合やＮＬＯＳ状態の場合、又は直接波と反射波が混在して受信するような環境下であっても、少ない処理量で効率よく、かつ精度良く移動局の位置を推定することができる。
【００２５】
前記算出された移動局の推定位置を前記移動局の推定初期位置として再帰的に用いて、上記繰返し動作を繰り返してもよい。
【００２６】
前記推定位置算出部は、前記複数の基地局の観測座標の中央値を前記移動局の代表位置として算出してもよい。
【００２７】
前記推定位置算出部は、前記複数の基地局の観測座標の平均値を前記移動局の代表位置として算出してもよい。
【００２８】
前記推定位置算出部は、前記複数の基地局の観測座標の重み付き平均値を前記移動局の代表位置として算出してもよい。
【００２９】
前記推定位置算出部は、前記重み付き平均値の重みを前記観測距離の長さに対応して決定してもよい。
【００３０】
前記予め定められた移動局の推定初期位置は、前記複数の基地局に囲まれた範囲内の位置、又は前記複数の基地局から選択された一基地局から所定の範囲内の位置をデフォルト値として最初の繰返し処理に使用され、２回目以降の繰返し処理では、前回の繰返し処理で算出された移動局の推定位置を今回の繰返し処理の前記移動局の推定初期位置として再帰的に用いてもよい。
【００３１】
前記位置推定装置は、前記各基地局、前記移動局、前記複数の基地局及び前記移動局と通信可能なサーバ、又は前記移動局の位置を検知する位置検知装置のいずれかであってもよい。
【００３２】
また、上記課題を解決するために、本発明の別の観点によれば、位置が既知である複数の基地局と位置推定対象である移動局との一方向又は双方向通信により取得した所定の信号に基づき、距離計測部により前記複数の基地局と前記移動局との観測距離を計測し、該移動局の位置を推定する位置推定方法であって、計測された前記複数の基地局の観測距離に基づき、観測座標決定部により予め定められた移動局の推定初期位置の方向に各基地局の観測座標を決定する観測座標決定ステップと、前記観測座標決定ステップにより決定された前記複数の基地局の観測座標から移動局の代表位置を前記移動局の推定位置として算出する推定位置算出ステップと、を含むことを特徴とする位置推定方法が提供される。
【００３３】
また、上記課題を解決するために、本発明の別の観点によれば、位置が既知である複数の基地局と位置推定対象である移動局との一方向又は双方向通信により取得した所定の信号に基づき、距離計測部により前記複数の基地局と前記移動局との観測距離を計測し、該移動局の位置を推定する処理をコンピュータに実現させるためのプログラムであって、計測された前記複数の基地局の観測距離に基づき、観測座標決定部により予め定められた移動局の推定初期位置の方向に各基地局の観測座標を決定する観測座標決定処理と、前記観測座標決定処理により決定された前記複数の基地局の観測座標から移動局の代表位置を前記移動局の推定位置として算出する推定位置算出処理と、をコンピュータに実行させるためのプログラムが提供される。
【発明の効果】
【００３４】
以上説明したように本発明によれば、ＮＬＯＳ状態の場合や直接波と反射波が混在して受信するような環境下であっても、精度良く移動局の位置を推定することができる。
【図面の簡単な説明】
【００３５】
【図１】本発明の第１及び第２実施形態に係る位置推定システムの全体構成図である。
【図２】第１及び第２実施形態に係る位置検知装置の機能構成図である。
【図３】図３（ａ）〜図３（ｄ）は繰り返し回数が０〜３回、観測データが１個の場合の第１実施形態に係る位置推定の状態遷移図である。
【図４】図４（ａ）〜図４（ｄ）は繰り返し回数が０〜３回、観測データが複数個の場合の第１実施形態に係る位置推定の状態遷移図である。
【図５】第２実施形態に係る加重平均値の重みを説明するための図である。
【図６】ＴＯＡによる一般的な位置推定方法を説明するための図である。
【図７】観測誤差のヒストグラムを示した図である。
【図８】図８（ａ）は第１の基地局、図８（ｂ）は第１２の基地局での観測誤差のヒストグラムを示した図である。
【発明を実施するための形態】
【００３６】
以下に添付図面を参照しながら、本発明の好適な実施の形態について詳細に説明する。なお、本明細書及び図面において、実質的に同一の機能構成を有する構成要素については、同一の符号を付することにより重複説明を省略する。
【００３７】
まず、本発明の第１実施形態に係る位置推定システムについて説明する前に、本実施形態にて位置を推定する際の環境、及び位置推定に使用する観測データの性質について説明する。その後、本実施形態に係る位置推定システムによる観測データを用いた移動局の位置推定について説明する。
【００３８】
（位置推定が行われる環境）
屋内環境で便利なサービスを提供しようと考える場合、対象とする人や物の位置を正確に知る技術は重要である。その際、対象と一緒に持ち運べる無線センサ技術が注目されている。無線センサを用いた位置推定の実用例としては、カーナビ等で利用されるＧＰＳ（Ｇｌｏｂａｌ Ｐｏｓｉｔｉｏｎｉｎｇ Ｓｙｓｔｅｍ）がある。しかし、精度としては、推定位置に約１０ｍ以上の誤差が存在しており、屋内でのサービス提供への利用には適していない。
【００３９】
これに対して、超音波を用いて距離を測定する手法もある。この手法では、約数ｃｍ〜約十数ｃｍの精度を実現している。ただし、バッテリーの寿命が３日間から１週間と短く、伝播距離も短い。このように、屋内において高精度で安定した位置推定方法は、未だ確立されていないといえる。
【００４０】
本実施形態では、ＣＳＳ（Ｃｈｉｒｐ Ｓｐｒｅａｄ Ｓｐｅｃｔｒｕｍ）と呼ばれる無線技術を用いる。これは、屋内及び屋外のどちらでも利用可能であり、バッテリーの寿命は数年である。伝播距離も超音波より長い。精度は、約２ｍ〜５ｍである。しかしながら、統計的なアプローチを加えることにより、屋内や構内で利用可能な数十ｃｍ単位という高精度な位置推定の実現も可能である。
【００４１】
本実施形態では、基地局から観測ポイント（すなわち、移動局）までの距離をＴｏＡ（Ｔｉｍｅ ｏｆ Ａｒｒｉｖａｌ）方式でリアルタイムに測定する。ＴｏＡ方式では、移動局と１台又は複数台の基地局との間で、無線信号の伝播遅延時間に基づいて距離を測定する。ここでは、図６に示したように、移動局１０と３台の基地局２０ａ、２０ｂ、２０ｃとの間で、移動局と各基地局間の距離を観測する。時刻ｔ０に移動局１０から送られた無線信号が、時刻ｔ１に基地局２０ａに到達した場合を考える。この場合、移動局１０から基地局２０ａまでの観測距離ｄ１は、無線信号の伝達速度ｖを用いて次の式にて求められた距離であると測定される。
ｄ１＝ｖ（ｔ１−ｔ０）
【００４２】
移動局１０から基地局２０ｂ、２０ｃまでの観測距離ｄ２、ｄ３も同様にして次の式にて表される。
ｄ２＝ｖ（ｔ２−ｔ０）
ｄ３＝ｖ（ｔ３−ｔ０）
【００４３】
推定したい位置に移動局１０を配置し、基地局２０ａ、２０ｂ、２０ｃを実験対象となる部屋の指定位置に設置する。これにより、部屋の様々な場所から計った観測距離値を取得することができる。もし、どの基地局から得られた観測距離値にも誤差が含まれない場合、従来の位置推定方法では、三点測量の原理により移動局１０の位置を特定することができる。図６に示した移動局と各基地局との位置関係の場合、三点測量では、基地局２０ａ、２０ｂ、２０ｃの座標から観測距離を半径とする円を描き、すべての点が交わった位置を移動局１０の推定位置とする。
【００４４】
しかしながら、実際には観測距離には観測誤差が含まれるため、図６に示したように各基地局２０ａ、２０ｂ、２０ｃを中心とした半径ｄ１、ｄ２、ｄ３の円が１箇所で交わることはほとんどない。このような場合、交点の直線の式や重心、最小自乗法の考え方を用いて幾何学的に位置を求める方法も考案されているが、観測誤差の影響を考えていないため、屋内でも使える程の高精度は期待できない。
【００４５】
よって、本実施形態では、実際の観測データから観測値の性質を見出し、それに応じた位置推定方法を提供する。観測値の性質については、実験結果から以下のことがわかっている。
（１）観測誤差は正の偏りを示す。
（２）観測誤差は移動局と基地局との位置関係で決まる。
【００４６】
（１）については、各基地局から移動局までの観測距離と実際の距離の差を調べると、９９％以上の観測データにおいて実際の距離より観測距離の方が大きい値を持つという結果が出た。それを示したグラフが図７である。図７は、サンプル数Ｎが２９４２７４個の場合の観測距離から実際の距離を引いた観測距離の誤差ヒストグラムである。この結果を見ると、大部分の誤差が０ｍよりも大きくなっていることが分かる。観測データを算出する観測装置の内部の計算時間や、空気等の媒体が要因であると考えられる。
【００４７】
次に、大きな測定誤差が出現する原因について考察すると、ＮＬＯＳ(見通し外)による誤差が主要因ではないかと考えられる。たとえばオフィス空間などの屋内環境で人や物体を対象とした位置検知を行う場合では、室内の壁や棚などの構造物により、移動局から基地局までの直接波の経路が構成されない場合等が起こりえる。そのため、直接波の代わりに壁等に反射した反射波が観測されて観測誤差が大きく測定されてしまうことが考えられる。
【００４８】
また、図８（ａ）の第１の基地局における、サンプル数Ｎ＝１９８個の場合、及び図８（ｂ）の第１２の基地局、サンプル数Ｎ＝１５７個の場合の観測距離から実際の距離を引いた観測距離の誤差ヒストグラムを見ると、第１及び第１２の基地局に対して常に同じ程度の誤差が生じており、ＮＬＯＳによる誤差が定常的に発生していることが分かる。実験から、移動局がわずか１０ｃｍ以下だけ動いた場合でも、観測誤差が大きく変わることも確認されている。
【００４９】
＜第１実施形態＞
（位置推定システム）
さて、以上の観測環境や観測値の性質を考慮した本実施形態の位置推定について説明する。ここでは、図１に示した位置推定システム１００を例に挙げて説明する。位置推定システム１００には、位置推定対象である移動局１０と、位置が既知である複数の基地局２０と、位置検知装置３０とを有する。移動局１０、複数の基地局２０及び位置検知装置３０は、一方向又は双方向通信可能である。
【００５０】
複数の基地局２０は、例えば屋内の無線ＬＡＮの基地局や携帯電話の基地局等と同様に設置されている。複数の基地局２０の位置は、既知である。基地局２０は、各基地局の位置を知る方法があれば固定される必要はなく、移動局１０と同様に可動式の形態としてもよい。
【００５１】
移動局１０としては、ＰＣ（Ｐｅｒｓｏｎａｌ Ｃｏｍｐｕｔｅｒ）やセンサーネットワークのノード、携帯電話、ＰＤＡ（Ｐｅｒｓｏｎａｌ Ｄｉｇｉｔａｌ Ａｓｓｉｓｔａｎｔ）、アクティブタグなど、様々な形態への適用が可能である。ここでは、移動局が１局、基地局を８局の例をあげて説明するが、移動局は複数でも良く、また、基地局も２局以上であれば何局であっても良い。
【００５２】
各基地局２０と移動局１０との観測距離は、位置検知装置３０により測定される。なお、位置検知装置３０は、移動局の位置を推定する位置推定装置の一例であり、位置推定システム１００には、必ずしも位置検知装置３０が基地局２０や移動局１０と別体で存在していなくてもよい。つまり、位置検知装置３０は単独で存在していてもよいし、基地局２０に内蔵されていてもよいし、移動局１０に内蔵されていてもよい。位置検知装置３０は、複数の基地局２０及び移動局１０と通信可能な図示しないサーバに内蔵されていてもよい。したがって、位置推定装置は、位置検知装置３０として存在していてもよく、基地局２０として存在していてもよく、移動局１０として存在していてもよく、図示しないサーバとして存在していてもよい。
【００５３】
移動局の推定初期位置４０は、予め定められたデフォルト値であり、本実施形態では、移動局の推定初期位置４０は、各基地局２０の重心位置である。
【００５４】
（位置検知装置の機能構成）
図２に、本実施形態に係る位置検知装置３０の機能構成を示す。位置検知装置３０は、距離計測部３２、観測座標決定部３４、推定位置算出部３６及び繰返し判定部３８を有している。
【００５５】
距離計測部３２は、複数の基地局２０と移動局１０との一方向又は双方向通信により取得した所定の信号に基づき、複数の基地局２０と移動局１０との観測距離をそれぞれ計測する。所定の信号は、例えば、一方向又は双方向通信の通信結果である無線信号の受発信時刻や基地局側での受信時刻の情報を含み、各基地局２０と移動局１０との観測距離の計測にこれらの情報を利用してもよい。ただし、各基地局２０もしくは移動局１０にて受発信時刻を特定できれば、必ずしも時刻情報を前記所定の信号に含ませる必要はない。
【００５６】
観測座標決定部３４は、距離計測部３２により計測された各基地局２０の観測距離に基づき、予め定められた移動局の推定初期位置４０の方向に各基地局２０の観測座標を決定する。
【００５７】
具体的には、観測座標決定部３４は、各基地局２０と移動局１０の推定初期位置４０とを結ぶ直線上であって、計測された基地局毎の観測距離だけ各基地局２０から離れた位置に各基地局の観測座標（ａｉ、ｂｉ）を決定する。
【００５８】
たとえば、図３（ａ）に示したように、右上に配置された基地局２０と移動局１０の推定初期位置４０とを結ぶ直線ｌ上であって、計測された基地局２０の観測距離ｄ１だけ基地局２０から離れた位置に右上に配置された基地局２０の観測座標（ａ１、ｂ１）を定める。他の基地局２０についても同様にして、図３（ａ）にそれぞれ黒丸で示した観測座標が決定される。
【００５９】
推定位置算出部３６は、観測座標決定部３４により決定された複数の基地局２０の観測座標（ａｉ、ｂｉ）から移動局１０の代表位置を算出し、これを移動局１０の推定位置５０とする。本実施形態では、推定位置算出部３６は、８つの基地局２０の８つの観測座標の中央値を移動局の代表位置として算出し、移動局１０の推定位置５０とする。
【００６０】
繰返し判定部３８は、推定位置算出部３６により算出された移動局の推定位置５０と移動局の推定初期位置４０との変化量が予め定められた停止条件を満たすまで、ここで算出された移動局の推定位置５０を移動局の推定初期位置４０として置き換えることにより再帰的に用いて、観測座標決定部３４、前記推定位置算出部３６及び繰返し判定部３８の動作を繰り返す。
【００６１】
このように、予め定められた移動局の推定初期位置４０は、本実施形態では複数の基地局２０の重心位置がデフォルト値として設定され、最初の繰返し処理（ｓ＝０）に使用され（図３（ａ））、２回目以降の繰返し処理（図３（ｂ）〜図３（ｄ））では、前回の繰返し処理で算出された最も新しい移動局の推定位置５０を今回の繰返し処理の移動局の推定初期位置４０として再帰的に用いる。
【００６２】
繰返し判定部３８は、推定位置算出部３６により算出された移動局の推定位置５０と移動局の推定初期位置４０との変化量が予め定められた停止条件を満たしたとき、最新の移動局の推定位置５０を移動局の位置と判定する。停止条件については後述する。
【００６３】
なお、位置検知装置３０の各部への指令は、専用の制御デバイスあるいはプログラムを実行する図示しないＣＰＵにより実行される。推定初期位置や各種プログラムは、図示しないＲＯＭや不揮発性メモリに予め記憶されている。ＣＰＵが、これらのメモリから各プログラムやデータを読み出し実行することにより、距離計測部３２、観測座標決定部３４、推定位置算出部３６及び繰返し判定部３８の各機能が実行される。
【００６４】
（位置推定方法）
次に、本実施形態に係る位置検知装置３０において移動局１０の位置を推定する方法について説明する。前述したように、移動局の推定初期位置４０は、移動局１０の推定位置の初期値であり、本実施形態の位置推定方法が実行される前に予め決められている。本実施形態では、移動局の推定初期位置４０は、複数の基地局２０の重心位置であるが、これに限られず複数の基地局２０から選択された一基地局から所定の範囲内（例えば信頼性の高い基地局から１ｍ以内）の任意の位置であってもよいし、複数の基地局２０に囲まれた範囲Ａ内の任意の位置であってもよい。予め定められた移動局の推定初期位置４０は、後述する最初（１回目）の繰返し処理に使用される。
【００６５】
(入力／出力)
推定に用いる基地局の数をＫ、各基地局の二次元座標を（ｘｉ，ｙｉ）、ｉ＝１、・・・、Ｋ、各基地局からの観測距離ｄｉ、ｉ＝１、・・・、Ｋを入力とし、推定した位置の二次元座標（ｘ，ｙ）を出力とする。本実施形態に係る位置推定方法では、以下のステップ１〜ステップ４を逐次的に行う。
【００６６】
（ステップ１：初期値の設定）
繰返しカウンタＳを０に設定し、下記式（１）、式（２）に基づき、すべての基地局の座標からの重心位置（ｘ（０）、ｙ（０））を、移動局の推定初期位置４０として算出する。
【数１】
【００６７】
（ステップ２：観測座標の算出）
下記式（３）、式（４）に基づき、基地局の座標（ｘｉ、ｙｉ）から移動局の推定位置の座標（ｘ（ｓ）、ｙ（ｓ））の方向に観測距離ｄｉ分離れた位置にある観測座標（ａｉ、ｂｉ）を求める。
【数２】
【００６８】
（ステップ３：推定座標の更新）
下記式（５）、式（６）に基づき、求めた観測座標（ａｉ、ｂｉ）の中央値を算出し、移動局の推定位置５０とする。
【数３】
【００６９】
（ステップ４：収束判定）
算出された移動局の推定位置５０と移動局の推定初期位置４０との変化量が予め定められた停止条件を満たしたとき、算出された移動局の推定位置５０を移動局の位置と判定し、その推定位置５０の座標（ｘ、ｙ）＝（ｘ（Ｓ＋１）、ｙ（Ｓ＋１））として出力座標を生成し、繰返し処理を終了する。
【００７０】
停止条件としては、算出された移動局の推定位置５０と移動局の推定初期位置４０との変化量との差分や比が挙げられる。本実施形態では、算出された移動局の推定位置５０と移動局の推定初期位置４０との変化量が、次の式（７）で示した停止条件を満たしたとき、算出された移動局の推定位置５０を移動局の位置と判定し、繰返し処理を終了する。
＜停止条件＞
【数４】
なお、δは十分小さい数で、ここではδ＝０．０５ｍとする。
【００７１】
上記ステップ１〜４の実行に伴う状態変化を図３（ａ）〜図３（ｄ）を用いて視覚的に説明する。図３（ａ）〜図３（ｄ）は、繰返しカウンタＳ（ＳｔｅｐＮｕｍ）が０回〜３回の場合の移動局の推定位置５０の状態遷移を示す。
【００７２】
図３（ａ）に示したように、まず、ステップ１では、すべての基地局２０の重心位置がｓ＝０のときの移動局の推定初期位置４０となる。続いて、ステップ２により各基地局２０から推定初期位置４０の方向に観測距離分伸ばした点（図３では各黒丸）を求める。なお、ここでは１個の基地局につき１個の観測データ（１個の黒丸）を使った場合を示している。
【００７３】
次に、ステップ３では、すべての黒丸の座標について中央値を取り、この中央値を推定位置５０とする。ステップ４で「繰返し処理を行う（ステップ２に戻る）」と判定された場合、推定位置５０は、次時点（次の繰返し処理）の推定初期位置４０となる。つまり、図３（ａ）の推定位置５０が、図３（ｂ）の推定初期位置４０となる。
【００７４】
次に、ステップ４では、算出された移動局の推定位置５０と移動局の推定初期位置４０との距離を計算し、上記停止条件を満たすように十分短くなれば処理を終了する。上記停止条件を満たさない間は、繰返しカウンタＳに１を加算して（繰返しカウンタＳ＝１）、ステップ２に戻る。
【００７５】
次にステップ２〜４を実行した結果、図３（ｂ）に示したように、算出された移動局の推定位置５０と移動局の推定初期位置４０との距離は短くなるが、上記停止条件を満たさないため、さらに繰返しカウンタＳに１を加算して（繰返しカウンタＳ＝２）、ステップ２に戻る。
【００７６】
次にステップ２〜４を実行した結果、図３（ｃ）に示したように、算出された移動局の推定位置５０と移動局の推定初期位置４０との距離はさらに短くなるが、上記停止条件を満たさないため、さらに繰返しカウンタＳに１を加算して（繰返しカウンタＳ＝３）、ステップ２に戻る。
【００７７】
次にステップ２〜４を実行した結果、図３（ｄ）に示したように、算出された移動局の推定位置５０と移動局の推定初期位置４０との距離は十分に短くなり、上記停止条件を満たす。この結果、最後に算出された移動局の推定位置５０の座標（ｘ、ｙ）＝（ｘ（Ｓ＋１）、ｙ（Ｓ＋１））として出力座標を生成し、繰返し処理を終了する。
【００７８】
以上の説明では１個の基地局につき１個の観測データ（１個の黒丸）を使った場合について説明したが、図４（ａ）〜図４（ｄ）に示したように、１個の基地局につき複数個の観測データ（複数個の黒丸）を使った場合についても、本実施形態に係る位置推定方法を用いることができる。１個の基地局につき複数個の観測データを使う場合、１個の観測データを使う場合より、移動局の推定位置の精度を高めることができる。
【００７９】
以上に説明したように、本実施形態に係る位置検知装置３０によれば、移動局の推定位置を予め適当な推定初期位置(座標)に仮決めし、観測した複数の基地局２０の観測距離に基づき、仮決めした推定初期位置（座標）の方向に各基地局２０の観測座標を算出する。算出された観測座標の集合から移動局１０の代表位置を算出し、その値を移動局の推定位置５０とする。
【００８０】
このようにして求められた最新の移動局の推定位置５０の変化量が十分に小さくなれば、移動局の推定位置として収束したと判定する。以上の位置推定方法により、ＮＬＯＳ状態の場合や直接波と反射波が混在して受信するような環境下であっても、精度良く移動局１０の位置を推定することができる。
【００８１】
また、本実施形態によれば、三角測量をする必要も生じないため、少ない処理量で効率よく移動局１０の位置を推定することができる。
【００８２】
（中央値の効果）
特に、本実施形態では、複数の基地局２０の観測座標の中央値を移動局の代表位置（推定位置５０）として算出する。この効果について説明する。本実施形態では、移動局１０の観測座標を基地局数に応じて複数求め、その複数の観測座標の分布から推定座標（移動局の推定位置５０）を求め、さらに推定座標を用いて基地局毎に求める観測座標（ａｉ，ｂｉ）を更新する処理を繰返し、再帰的に推定座標に含まれる誤差の収束を図っている。
【００８３】
本実施形態では、複数の観測座標から推定座標を求める基準として中央値を求めることが容易にできる。そして、誤差を収束させる繰返し処理の過程で、常に基地局数分の観測座標サンプルが得られるため、複数の観測座標の中央値を用いて絶対的な誤差を低減させることができる。言い換えると、誤差の収束過程では、中央値は絶対誤差を最小化する基準ということができ、再帰的な誤差の収束過程に中央値を使えるということは、異常値等の影響を常に排除しながら誤差の収束を図ることができるということである。
【００８４】
一方、複数の観測座標から推定座標を求める基準として平均値を求め、求められた平均値を用いて誤差の収束を行うと、異常な値等すべてを含めた上での２乗誤差を最小化することとなり、最終的に得られる推定位置５０に、異常値が大きな影響を及ぼすこととなる。
【００８５】
また、収束前にどの基地局の値が異常値なのかを決めることは難しく、例えば別のやり方で複数の観測座標の分布から、例えば平均値と標準偏差を用いた分布を仮定してそこから大きく外れるものを異常値と定め、異常値以外の観測座標から平均値を求めて推定座標とするといった処理も考えられるが、誤差収束の過程では、前提としている誤差の分布が既知ではないため、誤ってただしいものと正しくないものを取り違えてしまう可能性が高い。
【００８６】
誤差収束の過程で、常に中央値が使えるということは、誤差増大の要因となりうる異常値除去等を用いず、すべての観測座標を用い、それらを繰り返し更新し続けることで、最終的に最も誤差の小さい値を抽出することができるようになる位置推定方法といえる。
【００８７】
たとえば、ＮＬＯＳ等による異常に大きな誤差の混入した観測距離が得られる状態においても、異常な誤差値による影響を排除しやすく（中央値を基準に使った場合、大きな観測誤差を返すデータの重みを、平均値よりも小さくすることができる）、精度良く推定位置を求めることができる。
【００８８】
このように、中央値を簡単に採用できるのも、移動局の推定位置を基地局数に応じて複数求めることができる本実施形態の特徴から得られるものである。
【００８９】
ただし、複数の基地局の観測座標の平均値を移動局の代表位置（移動局の推定位置）として算出してもよい。この場合、２乗誤差を最小化する方法として機能でき、また別のフィルタ演算を用いることも可能である。また、複数の基地局の観測座標の重心位置を移動局の代表位置（移動局の推定位置５０）として算出してもよい。
【００９０】
また、複数の基地局の観測座標の重み付き平均値を移動局の代表位置（移動局の推定位置５０）として算出してもよい。これについては、次の第２実施形態にて詳しく説明する。
【００９１】
＜第２実施形態＞
本発明の第２実施形態に係る移動局の位置推定方法ついて説明する。本実施形態では、上記ステップ１〜ステップ４までの基本処理は、第１実施形態と同様である。ステップ２で求めた観測座標（ａｉ、ｂｉ）について、ステップ３にて以下のとおり加重平均を行い、次の（ｓ＋１）時点の推定座標とする点で、中央値を用いる第１実施形態と異なる。よって、この相違点を中心に第２実施形態について説明する。
【００９２】
【数４】
ただし、
【００９３】
また、重みＷｉ（ｓ）は、時刻ｓにおける観測距離ｄｉにより求まる値であり、以下の式（８）で定義する。
Ｗｉ（ｓ）＝ｆ（ｄｉ） ・・・（８）
【００９４】
これは、例えば図５に示したような平均値０、分散が例えば全基地局の重心位置までの距離とした、ガウス分布の半分を２倍にしたようなものを用いてもよい。図５は、観測距離に応じた基地局の重み付けの一例を示す。これによれば、重み付き平均値の重みは、観測距離の長さに対応して決定される。
【００９５】
図５に示した重みによれば、観測距離が短いほど重みが大きく、観測距離が長いほど重みが小さい。ここで、遠くの基地局による観測距離ほど、反射などの影響を受けやすく、誤差を多く含みやすい。よって、本実施形態によれば、このように重みをつけることにより、観測距離が長い場合に含まれる誤差をより低減しやすくすることができる。この結果、移動局の位置推定の精度をより向上させることができる。
【００９６】
以上、各実施形態に係る位置推定によれば、ＮＬＯＳ状態の場合や直接波と反射波が混在して受信するような環境下であっても、精度良く移動局の位置を推定することができる。
【００９７】
（応用例）
最後に各実施形態の応用例について述べる。通常、移動局１０は複数存在することもあるが、このような場合、移動局１０は移動局間でユニークなＩＤを発信することにより、基地局２０はどの移動局１０からの発信であるかを認識することができる。これにより、基地局２０は、移動局１０毎にデータを収集することができる。
【００９８】
また、基地局２０と移動局１０の通信上の関係は、位置検知装置３０の機能を複数の基地局２０の少なくともいずれかに内蔵して、基地局２０が移動局１０との観測距離を収集し、基地局２０が移動局１０の位置を推定するようにしてもよいし、逆に移動局１０が基地局２０との観測距離を収集し、移動局１０が自機の位置を推定するようにしてもよい。
【００９９】
通信方法は、電波に限られず、例えば赤外線や超音波、画像処理等を用いて基地局２０と移動局１０間の観測距離を求めてもよい。
【０１００】
時刻ｓにおけるステップ１の移動局の推定初期位置４０は、時刻（ｓ−１）で求まる推定位置としてもよいし、推定位置にランダム誤差を加えた座標としてもよい。
【０１０１】
なお、上記実施形態において、各部の動作は互いに関連しており、互いの関連を考慮しながら、一連の動作及び一連の処理として置き換えることができる。これにより、位置推定装置の実施形態を、位置推定方法の実施形態とすることができる。また、位置推定装置の実施形態を、これらが有する機能をコンピュータに実現させるためのプログラムの実施形態とすることができる。
【０１０２】
これにより、位置が既知である複数の基地局と位置推定対象である移動局との一方向又は双方向通信により取得した所定の信号に基づき、距離計測部により前記複数の基地局と前記移動局との観測距離を計測し、該移動局の位置を推定する位置推定方法であって、計測された前記複数の基地局の観測距離に基づき、観測座標決定部により予め定められた移動局の推定初期位置の方向に各基地局の観測座標を決定する観測座標決定ステップと、前記観測座標決定ステップにより決定された前記複数の基地局の観測座標から移動局の代表位置を前記移動局の推定位置として算出する推定位置算出ステップと、を含むことを特徴とする位置推定方法を提供することができる。
【０１０３】
また、これにより、位置が既知である複数の基地局と位置推定対象である移動局との一方向又は双方向通信により取得した所定の信号に基づき、距離計測部により前記複数の基地局と前記移動局との観測距離を計測し、該移動局の位置を推定する処理をコンピュータに実現させるためのプログラムであって、計測された前記複数の基地局の観測距離に基づき、観測座標決定部により予め定められた移動局の推定初期位置の方向に各基地局の観測座標を決定する観測座標決定処理と、前記観測座標決定処理により決定された前記複数の基地局の観測座標から移動局の代表位置を前記移動局の推定位置として算出する推定位置算出処理と、をコンピュータに実行させるためのプログラムを提供することができる。
【０１０４】
以上、添付図面を参照しながら本発明の好適な実施形態について詳細に説明したが、本発明はかかる例に限定されない。本発明の属する技術の分野における通常の知識を有する者であれば、特許請求の範囲に記載された技術的思想の範疇内において、各種の変更例または修正例に想到し得ることは明らかであり、これらについても、当然に本発明の技術的範囲に属するものと了解される。
【符号の説明】
【０１０５】
１０ 移動局
２０ 基地局
３０ 位置検知装置
３２ 距離計測部
３４ 観測座標決定部
３６ 推定位置算出部
３８ 繰返し判定部
４０ 移動局の推定初期位置
５０ 移動局の推定位置
１００ 位置推定システム
【技術分野】
【０００１】
本発明は、移動局の位置を推定する位置推定装置、位置推定方法及びプログラムに関する。より詳細には、本発明は、移動局と複数の基地局で通信を行い、その通信結果である無線信号の受発信時刻や、基地局側での受信時刻に基づいて移動局の位置を推定する位置推定装置、位置推定方法及びプログラムに関する。
【背景技術】
【０００２】
従来から、移動局と複数の基地局で通信を行い、その通信結果である無線信号や受発信時刻などに基づいて移動局の位置を推定する測位システムが提案されている。たとえば、特許文献１〜３には、複数の無線基地局と移動局との間で無線信号を送受信し、移動局の位置を推定するシステムが提案されている。
【０００３】
無線信号としては、電波、音波、光、磁気などを用いたものがある。それらを用いて位置を求める方式のうち主要なものを大別すると、信号の受信強度を用いるＲＳＳＩ（Ｒｅｃｅｉｖｅｄ Ｓｉｇｎａｌ Ｓｔｒｅｎｇｔｈ Ｉｎｄｉｃａｔｏｒ）方式（特許文献１参照）、信号の伝播時間を用いるＴＯＡ（Ｔｉｍｅ Ｏｆ Ａｒｒｉｖａｌ）方式（特許文献２参照）、信号の伝播時間差を用いるＴＤＯＡ（Ｔｉｍｅ Ｄｉｆｆｅｒｅｎｃｅ Ｏｆ Ａｒｒｉｖａｌ）方式、信号の受信方位を用いるＡＯＡ（Ａｎｇｌｅ Ｏｆ Ａｒｒｉｖａｌ）方式（特許文献３参照）の４方式が挙げられる。
【０００４】
ＲＳＳＩ方式は、特許文献１に示したように、信号の受信強度から距離を推定し、そこから３辺測量で位置を推定する方式である。この方式は、高精度に距離を推定することが難しい一方で、厳密かつ高精度な時刻管理が不要で、様々な機器に造りこみやすいといった特徴がある。受信強度から距離を推定することが難しいため、複数の受信信号の強弱のパターンの学習と識別により位置を推定するものもあるが、これには付近に学習済のエリアがなければ位置を推定できないといった課題がある。
【０００５】
ＴＯＡ方式は、信号の伝播時間から距離を計測し、そこから３辺測量で位置を推定する方式である。例えば、移動局から基地局に向かって信号を送信した際の信号の伝播時間をｔとし（通信方向が逆の場合や、双方向通信のものもある）、信号の伝播速度をｃとすると、基地局と移動局間の距離ｌは、ｌ＝ｃ×ｔで求めることができる。無線信号の一例として電波を媒体としたものの場合、一般に空気中の電波の伝播速度は一定（約３０万ｋｍ／秒）のため、時計の精度が十分に高ければ安定した位置精度を実現しやすい。例えば移動局と基地局の全ての時計が同期している場合、伝播時間ｔの計測は１回の通信で行なえる。他にも同期なしで距離を求める方法もあるが、その場合には、移動局から発信した電波が基地局に到着し、更に基地局から移動局にて受信した電波を検出する必要があるため処理の負荷が高くなる。
【０００６】
例えばＴＯＡ方式で用いられる測距の一方式であるＴＷＲ（Ｔｗｏ Ｗａｙ Ｒａｎｇｉｎｇ）方式では、移動局から基地局に信号を送出し、信号を受信した基地局は直ちに移動局に返信を行ない、移動局は、自局が信号を送出してから基地局からの信号を受信するまでの往復伝播時間から、基地局で折り返しに必要な処理時間を減算して２で割ることにより、片道にかかった伝播時間を求めることができる。いずれにしても、信号の伝播時間を計測するため、高精度な内部時計は必要である。たとえば電波を用いた場合、約１ｎｓの時刻誤差で約３０ｃｍの距離誤差が生じる。
【０００７】
ＴＤＯＡ方式は、ＴＯＡ方式と同様に信号の伝播時間をベースとしているが、基地局と移動局の距離を計測せずに、移動局からマルチキャスト発信した信号を同時に複数の基地局にて受信し、その受信時間差に基づいて位置を求める。２つの基地局で受信した信号の受信時刻差は、信号の伝播速度を掛け合わせると、移動局から見た距離差となる。２つの基地局から同一の距離差を発生し得る位置は双曲線を描くため、基地局３局以上で同時に信号を受信すると、複数の双曲線が得られ、それらの交点を求めることにより位置が決定する。この方式は、ＴＯＡ方式と比較し、移動局が信号を受信する必要がない点、１回の測位に必要な信号の送信回数が少ない（通常、移動局からの１回の信号発信で位置が確定する）といった利点がある一方、基地局間で時刻を厳密に同期しなければならず、また基地局に内蔵する時計には必ず誤差があるため、時計の精度にあわせて逐次時刻同期を行なう必要がある。
【０００８】
ＡＯＡ方式は、複数の基地局が移動局から信号を受信した際の信号の到達角度に基づき、三角測量を行なう方式である。信号の到達角度は、例えば電波の場合は、アレイアンテナ、音波の場合はマイクロフォンアレイなどを用い、ＭＵＳＩＣ（ＭＵｌｔｉｐｌｅ ＳＩｇｎａｌ Ｃｌａｓｓｉｆｉｃａｔｉｏｎ）法などを用いることで推定できる。距離的に非常に近いアンテナ（マイクロフォン）間の受信時刻差を用いるため、非常に高い精度で個々の受信信号をサンプリングする必要があるといった難しさがある。
【０００９】
このように、各方式には様々な状況により得手、不得手があり、すべてのケースにおいて常に最良となる方式は存在しない。特に、これら無線信号を用いる方式の共通の課題として、反射波による影響がある。
【００１０】
例えば特許文献２では、ＴＯＡ方式で位置を推定するシステムにおいて、反射波の影響を低減する方式が提案されている。通常、基地局と移動局の間に見通し（ＬＯＳ： Ｌｉｎｅ ｏｆ ｓｉｇｈｔ）があり、直線で結ばれる経路を通る直接波で通信が行なわれる場合には、直線距離を求めることが可能である。しかしながら、現実には経路途中の構造物などにより反射が発生することがあり、このような反射波を用いて距離を計測すると、直線距離よりも長い距離が得られてしまう。これは、基地局と移動局を最短で結ぶ経路がお互いを直線で結ぶ直接波の通る経路であり、それ以外の反射波の通る経路はすべて直線よりも長くなるためである。しかし、この関係を逆に利用することで、なるべく直接波の経路に近いものを選ぶことができる。例えば１回の通信の中で直接波と反射波が混在する合成波の場合は、一番初めに到達した成分が直接波による成分であり、この成分により伝播時間を求めることができれば、直線距離を算出することができる。また例えば複数回、信号の伝播時間を計測でき、その中に反射波と直接波による計測値が混在している場合は、最小の計測値を選ぶことで直接波による計測結果を選ぶことができ、それに基づき直線距離を得ることができる。また、この最小値を選ぶという処理は、たとえ直接波による計測結果がなく、すべて反射波によるものであったとしても、直接波と比較して最も差の少ない反射波を選ぶことになり、距離計測の誤差を低減するような直線距離に近い値を選ぶことができる。通信の受信品質をあわせて用いると、多重反射波などの品質の悪い通信を排除でき、より直接波に近いものを選ぶことも可能かもしれない。長期間にわたって継続して計測を行なう場合、移動局の基地局から見た平均移動速度を算出することで、よりもっともらしい位置を算出することも可能である。
【００１１】
特許文献３では、ＴＤＯＡ方式とＡＯＡ方式を組み合わせることで、位置精度を高めようとする技術が開示されている。ＡＯＡ方式は、角度の分解能が一定の場合、基地局と移動局の距離が離れるに従い、位置の分解能が低下する。一方で基地局と距離が十分に近い場合は、高い位置分解能を有する。かかる特徴により、信号の受信時刻が早いものが、移動局に最も近い基地局であり、位置分解能が高いとし、高い重み付けを行なう。このように時刻により重みを変えて加重平均をとることで、よりもっともらしい位置の推定を行なう。また、特許文献３では、受信時刻を得るため、ＴＤＯＡ方式の測位もあわせて行い、それらを重み付けによって加重平均を取ることにより、もっともらしい位置を計測する方式が提案されている。
【先行技術文献】
【特許文献】
【００１２】
【特許文献１】特開２００９−８５７８０号公報
【特許文献２】特開２００１−２７５１４８号公報
【特許文献３】特開２００７−１３５００号公報
【発明の概要】
【発明が解決しようとする課題】
【００１３】
このように、従来の位置推定方法では、基地局と移動局の距離を、電波等を利用して基地局毎に観測し、複数の基地局から得られる観測距離を用いて移動局の位置を推定していた。
【００１４】
しかし、上記従来手法では、例えば基地局と移動局の間に見通しがなく（このような状態を以下、ＮＬＯＳ（Ｎｏｎ−Ｌｉｎｅ−Ｏｆ−Ｓｉｇｈｔ）と呼ぶ。）、直接波が届かずに反射波で伝播経路が構成されるとき等、直接波による測定が十分に行なえない場合や、劣悪な伝播環境等の影響で異常な観測値を出力する基地局がある場合に、位置推定の精度が大きく劣化していた。
【００１５】
例えば特許文献２の方法では、ＮＬＯＳの状態にあったとしても、時間的な経過により移動局の移動や環境の変化などで、ＮＬＯＳの状態が一時的にでも解消し、直接波が到達すると仮定している。例えば屋外で移動局が道路上を継続的に移動しており、ビルなどの建築物の影響で一時的にＮＬＯＳ状態となっている場合は、時間的経過によってＬＯＳ状態（Ｌｉｎｅ ｏｆ ｓｉｇｈｔ：見通しがある状態）とＮＬＯＳ状態が交互に遷移し、長期的にみれば直接波による距離計測が可能な場合も多くある。このような場合では、文献２でもある程度正しく位置を推定可能である。しかし、たとえばオフィス空間などの屋内環境で人や物体を対象とした位置検知を行う場合では、室内の壁や棚などの構造物によりＮＬＯＳ状態になると、検知対象が継続的に移動しない場合が多いため、なかなかＮＬＯＳ状態が解消しないといった課題がある。特許文献２の方法では、このような課題に対処できない。たとえば電波の場合、電波暗室でない屋内など周囲に電波を良く反射する壁面のある環境などでは、直接波と反射波が複雑に混在することになるが、このような場合に電波の強度や雑音などの問題で、直接波を取り逃がす場合がよく起こり得る。このような場合、直接波を受信できる頻度が問題となり、頻度が低い場合には、十分な精度で位置測位が出来ないといった問題も生じていた。
【００１６】
また、特許文献３の方法では、特にＡＯＡ方式の分解能が高い基地局近くに移動局がある場合に反射波を拾ってしまうと、移動局が遠方にある場合よりも信号の入射角が大きく変わり、位置精度を大きく低下させる原因となっていた。
【００１７】
上記問題に鑑み、本発明は、ＮＬＯＳの状態の場合や直接波と反射波が混在して受信するような環境下であっても、高い精度で移動局の位置を推定することを目的とする。
【課題を解決するための手段】
【００１８】
上記課題を解決するために、本発明のある観点によれば、位置が既知である複数の基地局と位置推定対象である移動局との一方向又は双方向通信により取得した所定の信号に基づき、前記複数の基地局と前記移動局との観測距離を計測し、該移動局の位置を推定する位置推定装置であって、計測された前記複数の基地局の観測距離に基づき、予め定められた移動局の推定初期位置の方向に各基地局の観測座標を決定する観測座標決定部と、前記観測座標決定部により決定された前記複数の基地局の観測座標から移動局の代表位置を前記移動局の推定位置として算出する推定位置算出部と、を備えることを特徴とする位置推定装置が提供される。
【００１９】
かかる構成によれば、移動局の推定位置を予め適当な推定初期位置(座標)に仮決めし、観測した複数の基地局の観測距離に基づき、仮決めした推定初期位置（座標）の方向に各基地局の観測座標を決定する。
【００２０】
たとえば、各基地局と前記移動局の推定初期位置とを結ぶ直線上であって、計測された観測距離だけ各基地局から離れた位置に各基地局の観測座標を算出する。
【００２１】
このようにして、個々の基地局の観測距離に基づき推定初期位置の誤差成分を表すために観測座標を基地局毎に算出する。推定位置算出部は、複数の基地局の観測座標の集合から移動局の代表位置を算出し、その値を移動局の推定位置とする。
【００２２】
前記推定位置算出部により算出された移動局の推定位置と前記移動局の推定初期位置との変化量が予め定められた停止条件を満たすまで、前記距離計測部、前記観測座標決定部及び前記推定位置算出部の動作を繰り返す繰返し判定部をさらに有していてもよい。
【００２３】
前記繰返し判定部は、前記推定位置算出部により算出された移動局の推定位置と前記移動局の推定初期位置との変化量が予め定められた停止条件を満たしたとき、前記移動局の推定位置を前記移動局の位置と判定してもよい。
【００２４】
このようにして移動局の推定位置の変化量が十分に小さくなれば、処理が収束したとして最後に算出された移動局の推定位置を移動局の位置と判定する。これにより、異常な基地局が存在する場合やＮＬＯＳ状態の場合、又は直接波と反射波が混在して受信するような環境下であっても、少ない処理量で効率よく、かつ精度良く移動局の位置を推定することができる。
【００２５】
前記算出された移動局の推定位置を前記移動局の推定初期位置として再帰的に用いて、上記繰返し動作を繰り返してもよい。
【００２６】
前記推定位置算出部は、前記複数の基地局の観測座標の中央値を前記移動局の代表位置として算出してもよい。
【００２７】
前記推定位置算出部は、前記複数の基地局の観測座標の平均値を前記移動局の代表位置として算出してもよい。
【００２８】
前記推定位置算出部は、前記複数の基地局の観測座標の重み付き平均値を前記移動局の代表位置として算出してもよい。
【００２９】
前記推定位置算出部は、前記重み付き平均値の重みを前記観測距離の長さに対応して決定してもよい。
【００３０】
前記予め定められた移動局の推定初期位置は、前記複数の基地局に囲まれた範囲内の位置、又は前記複数の基地局から選択された一基地局から所定の範囲内の位置をデフォルト値として最初の繰返し処理に使用され、２回目以降の繰返し処理では、前回の繰返し処理で算出された移動局の推定位置を今回の繰返し処理の前記移動局の推定初期位置として再帰的に用いてもよい。
【００３１】
前記位置推定装置は、前記各基地局、前記移動局、前記複数の基地局及び前記移動局と通信可能なサーバ、又は前記移動局の位置を検知する位置検知装置のいずれかであってもよい。
【００３２】
また、上記課題を解決するために、本発明の別の観点によれば、位置が既知である複数の基地局と位置推定対象である移動局との一方向又は双方向通信により取得した所定の信号に基づき、距離計測部により前記複数の基地局と前記移動局との観測距離を計測し、該移動局の位置を推定する位置推定方法であって、計測された前記複数の基地局の観測距離に基づき、観測座標決定部により予め定められた移動局の推定初期位置の方向に各基地局の観測座標を決定する観測座標決定ステップと、前記観測座標決定ステップにより決定された前記複数の基地局の観測座標から移動局の代表位置を前記移動局の推定位置として算出する推定位置算出ステップと、を含むことを特徴とする位置推定方法が提供される。
【００３３】
また、上記課題を解決するために、本発明の別の観点によれば、位置が既知である複数の基地局と位置推定対象である移動局との一方向又は双方向通信により取得した所定の信号に基づき、距離計測部により前記複数の基地局と前記移動局との観測距離を計測し、該移動局の位置を推定する処理をコンピュータに実現させるためのプログラムであって、計測された前記複数の基地局の観測距離に基づき、観測座標決定部により予め定められた移動局の推定初期位置の方向に各基地局の観測座標を決定する観測座標決定処理と、前記観測座標決定処理により決定された前記複数の基地局の観測座標から移動局の代表位置を前記移動局の推定位置として算出する推定位置算出処理と、をコンピュータに実行させるためのプログラムが提供される。
【発明の効果】
【００３４】
以上説明したように本発明によれば、ＮＬＯＳ状態の場合や直接波と反射波が混在して受信するような環境下であっても、精度良く移動局の位置を推定することができる。
【図面の簡単な説明】
【００３５】
【図１】本発明の第１及び第２実施形態に係る位置推定システムの全体構成図である。
【図２】第１及び第２実施形態に係る位置検知装置の機能構成図である。
【図３】図３（ａ）〜図３（ｄ）は繰り返し回数が０〜３回、観測データが１個の場合の第１実施形態に係る位置推定の状態遷移図である。
【図４】図４（ａ）〜図４（ｄ）は繰り返し回数が０〜３回、観測データが複数個の場合の第１実施形態に係る位置推定の状態遷移図である。
【図５】第２実施形態に係る加重平均値の重みを説明するための図である。
【図６】ＴＯＡによる一般的な位置推定方法を説明するための図である。
【図７】観測誤差のヒストグラムを示した図である。
【図８】図８（ａ）は第１の基地局、図８（ｂ）は第１２の基地局での観測誤差のヒストグラムを示した図である。
【発明を実施するための形態】
【００３６】
以下に添付図面を参照しながら、本発明の好適な実施の形態について詳細に説明する。なお、本明細書及び図面において、実質的に同一の機能構成を有する構成要素については、同一の符号を付することにより重複説明を省略する。
【００３７】
まず、本発明の第１実施形態に係る位置推定システムについて説明する前に、本実施形態にて位置を推定する際の環境、及び位置推定に使用する観測データの性質について説明する。その後、本実施形態に係る位置推定システムによる観測データを用いた移動局の位置推定について説明する。
【００３８】
（位置推定が行われる環境）
屋内環境で便利なサービスを提供しようと考える場合、対象とする人や物の位置を正確に知る技術は重要である。その際、対象と一緒に持ち運べる無線センサ技術が注目されている。無線センサを用いた位置推定の実用例としては、カーナビ等で利用されるＧＰＳ（Ｇｌｏｂａｌ Ｐｏｓｉｔｉｏｎｉｎｇ Ｓｙｓｔｅｍ）がある。しかし、精度としては、推定位置に約１０ｍ以上の誤差が存在しており、屋内でのサービス提供への利用には適していない。
【００３９】
これに対して、超音波を用いて距離を測定する手法もある。この手法では、約数ｃｍ〜約十数ｃｍの精度を実現している。ただし、バッテリーの寿命が３日間から１週間と短く、伝播距離も短い。このように、屋内において高精度で安定した位置推定方法は、未だ確立されていないといえる。
【００４０】
本実施形態では、ＣＳＳ（Ｃｈｉｒｐ Ｓｐｒｅａｄ Ｓｐｅｃｔｒｕｍ）と呼ばれる無線技術を用いる。これは、屋内及び屋外のどちらでも利用可能であり、バッテリーの寿命は数年である。伝播距離も超音波より長い。精度は、約２ｍ〜５ｍである。しかしながら、統計的なアプローチを加えることにより、屋内や構内で利用可能な数十ｃｍ単位という高精度な位置推定の実現も可能である。
【００４１】
本実施形態では、基地局から観測ポイント（すなわち、移動局）までの距離をＴｏＡ（Ｔｉｍｅ ｏｆ Ａｒｒｉｖａｌ）方式でリアルタイムに測定する。ＴｏＡ方式では、移動局と１台又は複数台の基地局との間で、無線信号の伝播遅延時間に基づいて距離を測定する。ここでは、図６に示したように、移動局１０と３台の基地局２０ａ、２０ｂ、２０ｃとの間で、移動局と各基地局間の距離を観測する。時刻ｔ０に移動局１０から送られた無線信号が、時刻ｔ１に基地局２０ａに到達した場合を考える。この場合、移動局１０から基地局２０ａまでの観測距離ｄ１は、無線信号の伝達速度ｖを用いて次の式にて求められた距離であると測定される。
ｄ１＝ｖ（ｔ１−ｔ０）
【００４２】
移動局１０から基地局２０ｂ、２０ｃまでの観測距離ｄ２、ｄ３も同様にして次の式にて表される。
ｄ２＝ｖ（ｔ２−ｔ０）
ｄ３＝ｖ（ｔ３−ｔ０）
【００４３】
推定したい位置に移動局１０を配置し、基地局２０ａ、２０ｂ、２０ｃを実験対象となる部屋の指定位置に設置する。これにより、部屋の様々な場所から計った観測距離値を取得することができる。もし、どの基地局から得られた観測距離値にも誤差が含まれない場合、従来の位置推定方法では、三点測量の原理により移動局１０の位置を特定することができる。図６に示した移動局と各基地局との位置関係の場合、三点測量では、基地局２０ａ、２０ｂ、２０ｃの座標から観測距離を半径とする円を描き、すべての点が交わった位置を移動局１０の推定位置とする。
【００４４】
しかしながら、実際には観測距離には観測誤差が含まれるため、図６に示したように各基地局２０ａ、２０ｂ、２０ｃを中心とした半径ｄ１、ｄ２、ｄ３の円が１箇所で交わることはほとんどない。このような場合、交点の直線の式や重心、最小自乗法の考え方を用いて幾何学的に位置を求める方法も考案されているが、観測誤差の影響を考えていないため、屋内でも使える程の高精度は期待できない。
【００４５】
よって、本実施形態では、実際の観測データから観測値の性質を見出し、それに応じた位置推定方法を提供する。観測値の性質については、実験結果から以下のことがわかっている。
（１）観測誤差は正の偏りを示す。
（２）観測誤差は移動局と基地局との位置関係で決まる。
【００４６】
（１）については、各基地局から移動局までの観測距離と実際の距離の差を調べると、９９％以上の観測データにおいて実際の距離より観測距離の方が大きい値を持つという結果が出た。それを示したグラフが図７である。図７は、サンプル数Ｎが２９４２７４個の場合の観測距離から実際の距離を引いた観測距離の誤差ヒストグラムである。この結果を見ると、大部分の誤差が０ｍよりも大きくなっていることが分かる。観測データを算出する観測装置の内部の計算時間や、空気等の媒体が要因であると考えられる。
【００４７】
次に、大きな測定誤差が出現する原因について考察すると、ＮＬＯＳ(見通し外)による誤差が主要因ではないかと考えられる。たとえばオフィス空間などの屋内環境で人や物体を対象とした位置検知を行う場合では、室内の壁や棚などの構造物により、移動局から基地局までの直接波の経路が構成されない場合等が起こりえる。そのため、直接波の代わりに壁等に反射した反射波が観測されて観測誤差が大きく測定されてしまうことが考えられる。
【００４８】
また、図８（ａ）の第１の基地局における、サンプル数Ｎ＝１９８個の場合、及び図８（ｂ）の第１２の基地局、サンプル数Ｎ＝１５７個の場合の観測距離から実際の距離を引いた観測距離の誤差ヒストグラムを見ると、第１及び第１２の基地局に対して常に同じ程度の誤差が生じており、ＮＬＯＳによる誤差が定常的に発生していることが分かる。実験から、移動局がわずか１０ｃｍ以下だけ動いた場合でも、観測誤差が大きく変わることも確認されている。
【００４９】
＜第１実施形態＞
（位置推定システム）
さて、以上の観測環境や観測値の性質を考慮した本実施形態の位置推定について説明する。ここでは、図１に示した位置推定システム１００を例に挙げて説明する。位置推定システム１００には、位置推定対象である移動局１０と、位置が既知である複数の基地局２０と、位置検知装置３０とを有する。移動局１０、複数の基地局２０及び位置検知装置３０は、一方向又は双方向通信可能である。
【００５０】
複数の基地局２０は、例えば屋内の無線ＬＡＮの基地局や携帯電話の基地局等と同様に設置されている。複数の基地局２０の位置は、既知である。基地局２０は、各基地局の位置を知る方法があれば固定される必要はなく、移動局１０と同様に可動式の形態としてもよい。
【００５１】
移動局１０としては、ＰＣ（Ｐｅｒｓｏｎａｌ Ｃｏｍｐｕｔｅｒ）やセンサーネットワークのノード、携帯電話、ＰＤＡ（Ｐｅｒｓｏｎａｌ Ｄｉｇｉｔａｌ Ａｓｓｉｓｔａｎｔ）、アクティブタグなど、様々な形態への適用が可能である。ここでは、移動局が１局、基地局を８局の例をあげて説明するが、移動局は複数でも良く、また、基地局も２局以上であれば何局であっても良い。
【００５２】
各基地局２０と移動局１０との観測距離は、位置検知装置３０により測定される。なお、位置検知装置３０は、移動局の位置を推定する位置推定装置の一例であり、位置推定システム１００には、必ずしも位置検知装置３０が基地局２０や移動局１０と別体で存在していなくてもよい。つまり、位置検知装置３０は単独で存在していてもよいし、基地局２０に内蔵されていてもよいし、移動局１０に内蔵されていてもよい。位置検知装置３０は、複数の基地局２０及び移動局１０と通信可能な図示しないサーバに内蔵されていてもよい。したがって、位置推定装置は、位置検知装置３０として存在していてもよく、基地局２０として存在していてもよく、移動局１０として存在していてもよく、図示しないサーバとして存在していてもよい。
【００５３】
移動局の推定初期位置４０は、予め定められたデフォルト値であり、本実施形態では、移動局の推定初期位置４０は、各基地局２０の重心位置である。
【００５４】
（位置検知装置の機能構成）
図２に、本実施形態に係る位置検知装置３０の機能構成を示す。位置検知装置３０は、距離計測部３２、観測座標決定部３４、推定位置算出部３６及び繰返し判定部３８を有している。
【００５５】
距離計測部３２は、複数の基地局２０と移動局１０との一方向又は双方向通信により取得した所定の信号に基づき、複数の基地局２０と移動局１０との観測距離をそれぞれ計測する。所定の信号は、例えば、一方向又は双方向通信の通信結果である無線信号の受発信時刻や基地局側での受信時刻の情報を含み、各基地局２０と移動局１０との観測距離の計測にこれらの情報を利用してもよい。ただし、各基地局２０もしくは移動局１０にて受発信時刻を特定できれば、必ずしも時刻情報を前記所定の信号に含ませる必要はない。
【００５６】
観測座標決定部３４は、距離計測部３２により計測された各基地局２０の観測距離に基づき、予め定められた移動局の推定初期位置４０の方向に各基地局２０の観測座標を決定する。
【００５７】
具体的には、観測座標決定部３４は、各基地局２０と移動局１０の推定初期位置４０とを結ぶ直線上であって、計測された基地局毎の観測距離だけ各基地局２０から離れた位置に各基地局の観測座標（ａｉ、ｂｉ）を決定する。
【００５８】
たとえば、図３（ａ）に示したように、右上に配置された基地局２０と移動局１０の推定初期位置４０とを結ぶ直線ｌ上であって、計測された基地局２０の観測距離ｄ１だけ基地局２０から離れた位置に右上に配置された基地局２０の観測座標（ａ１、ｂ１）を定める。他の基地局２０についても同様にして、図３（ａ）にそれぞれ黒丸で示した観測座標が決定される。
【００５９】
推定位置算出部３６は、観測座標決定部３４により決定された複数の基地局２０の観測座標（ａｉ、ｂｉ）から移動局１０の代表位置を算出し、これを移動局１０の推定位置５０とする。本実施形態では、推定位置算出部３６は、８つの基地局２０の８つの観測座標の中央値を移動局の代表位置として算出し、移動局１０の推定位置５０とする。
【００６０】
繰返し判定部３８は、推定位置算出部３６により算出された移動局の推定位置５０と移動局の推定初期位置４０との変化量が予め定められた停止条件を満たすまで、ここで算出された移動局の推定位置５０を移動局の推定初期位置４０として置き換えることにより再帰的に用いて、観測座標決定部３４、前記推定位置算出部３６及び繰返し判定部３８の動作を繰り返す。
【００６１】
このように、予め定められた移動局の推定初期位置４０は、本実施形態では複数の基地局２０の重心位置がデフォルト値として設定され、最初の繰返し処理（ｓ＝０）に使用され（図３（ａ））、２回目以降の繰返し処理（図３（ｂ）〜図３（ｄ））では、前回の繰返し処理で算出された最も新しい移動局の推定位置５０を今回の繰返し処理の移動局の推定初期位置４０として再帰的に用いる。
【００６２】
繰返し判定部３８は、推定位置算出部３６により算出された移動局の推定位置５０と移動局の推定初期位置４０との変化量が予め定められた停止条件を満たしたとき、最新の移動局の推定位置５０を移動局の位置と判定する。停止条件については後述する。
【００６３】
なお、位置検知装置３０の各部への指令は、専用の制御デバイスあるいはプログラムを実行する図示しないＣＰＵにより実行される。推定初期位置や各種プログラムは、図示しないＲＯＭや不揮発性メモリに予め記憶されている。ＣＰＵが、これらのメモリから各プログラムやデータを読み出し実行することにより、距離計測部３２、観測座標決定部３４、推定位置算出部３６及び繰返し判定部３８の各機能が実行される。
【００６４】
（位置推定方法）
次に、本実施形態に係る位置検知装置３０において移動局１０の位置を推定する方法について説明する。前述したように、移動局の推定初期位置４０は、移動局１０の推定位置の初期値であり、本実施形態の位置推定方法が実行される前に予め決められている。本実施形態では、移動局の推定初期位置４０は、複数の基地局２０の重心位置であるが、これに限られず複数の基地局２０から選択された一基地局から所定の範囲内（例えば信頼性の高い基地局から１ｍ以内）の任意の位置であってもよいし、複数の基地局２０に囲まれた範囲Ａ内の任意の位置であってもよい。予め定められた移動局の推定初期位置４０は、後述する最初（１回目）の繰返し処理に使用される。
【００６５】
(入力／出力)
推定に用いる基地局の数をＫ、各基地局の二次元座標を（ｘｉ，ｙｉ）、ｉ＝１、・・・、Ｋ、各基地局からの観測距離ｄｉ、ｉ＝１、・・・、Ｋを入力とし、推定した位置の二次元座標（ｘ，ｙ）を出力とする。本実施形態に係る位置推定方法では、以下のステップ１〜ステップ４を逐次的に行う。
【００６６】
（ステップ１：初期値の設定）
繰返しカウンタＳを０に設定し、下記式（１）、式（２）に基づき、すべての基地局の座標からの重心位置（ｘ（０）、ｙ（０））を、移動局の推定初期位置４０として算出する。
【数１】
【００６７】
（ステップ２：観測座標の算出）
下記式（３）、式（４）に基づき、基地局の座標（ｘｉ、ｙｉ）から移動局の推定位置の座標（ｘ（ｓ）、ｙ（ｓ））の方向に観測距離ｄｉ分離れた位置にある観測座標（ａｉ、ｂｉ）を求める。
【数２】
【００６８】
（ステップ３：推定座標の更新）
下記式（５）、式（６）に基づき、求めた観測座標（ａｉ、ｂｉ）の中央値を算出し、移動局の推定位置５０とする。
【数３】
【００６９】
（ステップ４：収束判定）
算出された移動局の推定位置５０と移動局の推定初期位置４０との変化量が予め定められた停止条件を満たしたとき、算出された移動局の推定位置５０を移動局の位置と判定し、その推定位置５０の座標（ｘ、ｙ）＝（ｘ（Ｓ＋１）、ｙ（Ｓ＋１））として出力座標を生成し、繰返し処理を終了する。
【００７０】
停止条件としては、算出された移動局の推定位置５０と移動局の推定初期位置４０との変化量との差分や比が挙げられる。本実施形態では、算出された移動局の推定位置５０と移動局の推定初期位置４０との変化量が、次の式（７）で示した停止条件を満たしたとき、算出された移動局の推定位置５０を移動局の位置と判定し、繰返し処理を終了する。
＜停止条件＞
【数４】
なお、δは十分小さい数で、ここではδ＝０．０５ｍとする。
【００７１】
上記ステップ１〜４の実行に伴う状態変化を図３（ａ）〜図３（ｄ）を用いて視覚的に説明する。図３（ａ）〜図３（ｄ）は、繰返しカウンタＳ（ＳｔｅｐＮｕｍ）が０回〜３回の場合の移動局の推定位置５０の状態遷移を示す。
【００７２】
図３（ａ）に示したように、まず、ステップ１では、すべての基地局２０の重心位置がｓ＝０のときの移動局の推定初期位置４０となる。続いて、ステップ２により各基地局２０から推定初期位置４０の方向に観測距離分伸ばした点（図３では各黒丸）を求める。なお、ここでは１個の基地局につき１個の観測データ（１個の黒丸）を使った場合を示している。
【００７３】
次に、ステップ３では、すべての黒丸の座標について中央値を取り、この中央値を推定位置５０とする。ステップ４で「繰返し処理を行う（ステップ２に戻る）」と判定された場合、推定位置５０は、次時点（次の繰返し処理）の推定初期位置４０となる。つまり、図３（ａ）の推定位置５０が、図３（ｂ）の推定初期位置４０となる。
【００７４】
次に、ステップ４では、算出された移動局の推定位置５０と移動局の推定初期位置４０との距離を計算し、上記停止条件を満たすように十分短くなれば処理を終了する。上記停止条件を満たさない間は、繰返しカウンタＳに１を加算して（繰返しカウンタＳ＝１）、ステップ２に戻る。
【００７５】
次にステップ２〜４を実行した結果、図３（ｂ）に示したように、算出された移動局の推定位置５０と移動局の推定初期位置４０との距離は短くなるが、上記停止条件を満たさないため、さらに繰返しカウンタＳに１を加算して（繰返しカウンタＳ＝２）、ステップ２に戻る。
【００７６】
次にステップ２〜４を実行した結果、図３（ｃ）に示したように、算出された移動局の推定位置５０と移動局の推定初期位置４０との距離はさらに短くなるが、上記停止条件を満たさないため、さらに繰返しカウンタＳに１を加算して（繰返しカウンタＳ＝３）、ステップ２に戻る。
【００７７】
次にステップ２〜４を実行した結果、図３（ｄ）に示したように、算出された移動局の推定位置５０と移動局の推定初期位置４０との距離は十分に短くなり、上記停止条件を満たす。この結果、最後に算出された移動局の推定位置５０の座標（ｘ、ｙ）＝（ｘ（Ｓ＋１）、ｙ（Ｓ＋１））として出力座標を生成し、繰返し処理を終了する。
【００７８】
以上の説明では１個の基地局につき１個の観測データ（１個の黒丸）を使った場合について説明したが、図４（ａ）〜図４（ｄ）に示したように、１個の基地局につき複数個の観測データ（複数個の黒丸）を使った場合についても、本実施形態に係る位置推定方法を用いることができる。１個の基地局につき複数個の観測データを使う場合、１個の観測データを使う場合より、移動局の推定位置の精度を高めることができる。
【００７９】
以上に説明したように、本実施形態に係る位置検知装置３０によれば、移動局の推定位置を予め適当な推定初期位置(座標)に仮決めし、観測した複数の基地局２０の観測距離に基づき、仮決めした推定初期位置（座標）の方向に各基地局２０の観測座標を算出する。算出された観測座標の集合から移動局１０の代表位置を算出し、その値を移動局の推定位置５０とする。
【００８０】
このようにして求められた最新の移動局の推定位置５０の変化量が十分に小さくなれば、移動局の推定位置として収束したと判定する。以上の位置推定方法により、ＮＬＯＳ状態の場合や直接波と反射波が混在して受信するような環境下であっても、精度良く移動局１０の位置を推定することができる。
【００８１】
また、本実施形態によれば、三角測量をする必要も生じないため、少ない処理量で効率よく移動局１０の位置を推定することができる。
【００８２】
（中央値の効果）
特に、本実施形態では、複数の基地局２０の観測座標の中央値を移動局の代表位置（推定位置５０）として算出する。この効果について説明する。本実施形態では、移動局１０の観測座標を基地局数に応じて複数求め、その複数の観測座標の分布から推定座標（移動局の推定位置５０）を求め、さらに推定座標を用いて基地局毎に求める観測座標（ａｉ，ｂｉ）を更新する処理を繰返し、再帰的に推定座標に含まれる誤差の収束を図っている。
【００８３】
本実施形態では、複数の観測座標から推定座標を求める基準として中央値を求めることが容易にできる。そして、誤差を収束させる繰返し処理の過程で、常に基地局数分の観測座標サンプルが得られるため、複数の観測座標の中央値を用いて絶対的な誤差を低減させることができる。言い換えると、誤差の収束過程では、中央値は絶対誤差を最小化する基準ということができ、再帰的な誤差の収束過程に中央値を使えるということは、異常値等の影響を常に排除しながら誤差の収束を図ることができるということである。
【００８４】
一方、複数の観測座標から推定座標を求める基準として平均値を求め、求められた平均値を用いて誤差の収束を行うと、異常な値等すべてを含めた上での２乗誤差を最小化することとなり、最終的に得られる推定位置５０に、異常値が大きな影響を及ぼすこととなる。
【００８５】
また、収束前にどの基地局の値が異常値なのかを決めることは難しく、例えば別のやり方で複数の観測座標の分布から、例えば平均値と標準偏差を用いた分布を仮定してそこから大きく外れるものを異常値と定め、異常値以外の観測座標から平均値を求めて推定座標とするといった処理も考えられるが、誤差収束の過程では、前提としている誤差の分布が既知ではないため、誤ってただしいものと正しくないものを取り違えてしまう可能性が高い。
【００８６】
誤差収束の過程で、常に中央値が使えるということは、誤差増大の要因となりうる異常値除去等を用いず、すべての観測座標を用い、それらを繰り返し更新し続けることで、最終的に最も誤差の小さい値を抽出することができるようになる位置推定方法といえる。
【００８７】
たとえば、ＮＬＯＳ等による異常に大きな誤差の混入した観測距離が得られる状態においても、異常な誤差値による影響を排除しやすく（中央値を基準に使った場合、大きな観測誤差を返すデータの重みを、平均値よりも小さくすることができる）、精度良く推定位置を求めることができる。
【００８８】
このように、中央値を簡単に採用できるのも、移動局の推定位置を基地局数に応じて複数求めることができる本実施形態の特徴から得られるものである。
【００８９】
ただし、複数の基地局の観測座標の平均値を移動局の代表位置（移動局の推定位置）として算出してもよい。この場合、２乗誤差を最小化する方法として機能でき、また別のフィルタ演算を用いることも可能である。また、複数の基地局の観測座標の重心位置を移動局の代表位置（移動局の推定位置５０）として算出してもよい。
【００９０】
また、複数の基地局の観測座標の重み付き平均値を移動局の代表位置（移動局の推定位置５０）として算出してもよい。これについては、次の第２実施形態にて詳しく説明する。
【００９１】
＜第２実施形態＞
本発明の第２実施形態に係る移動局の位置推定方法ついて説明する。本実施形態では、上記ステップ１〜ステップ４までの基本処理は、第１実施形態と同様である。ステップ２で求めた観測座標（ａｉ、ｂｉ）について、ステップ３にて以下のとおり加重平均を行い、次の（ｓ＋１）時点の推定座標とする点で、中央値を用いる第１実施形態と異なる。よって、この相違点を中心に第２実施形態について説明する。
【００９２】
【数４】
ただし、
【００９３】
また、重みＷｉ（ｓ）は、時刻ｓにおける観測距離ｄｉにより求まる値であり、以下の式（８）で定義する。
Ｗｉ（ｓ）＝ｆ（ｄｉ） ・・・（８）
【００９４】
これは、例えば図５に示したような平均値０、分散が例えば全基地局の重心位置までの距離とした、ガウス分布の半分を２倍にしたようなものを用いてもよい。図５は、観測距離に応じた基地局の重み付けの一例を示す。これによれば、重み付き平均値の重みは、観測距離の長さに対応して決定される。
【００９５】
図５に示した重みによれば、観測距離が短いほど重みが大きく、観測距離が長いほど重みが小さい。ここで、遠くの基地局による観測距離ほど、反射などの影響を受けやすく、誤差を多く含みやすい。よって、本実施形態によれば、このように重みをつけることにより、観測距離が長い場合に含まれる誤差をより低減しやすくすることができる。この結果、移動局の位置推定の精度をより向上させることができる。
【００９６】
以上、各実施形態に係る位置推定によれば、ＮＬＯＳ状態の場合や直接波と反射波が混在して受信するような環境下であっても、精度良く移動局の位置を推定することができる。
【００９７】
（応用例）
最後に各実施形態の応用例について述べる。通常、移動局１０は複数存在することもあるが、このような場合、移動局１０は移動局間でユニークなＩＤを発信することにより、基地局２０はどの移動局１０からの発信であるかを認識することができる。これにより、基地局２０は、移動局１０毎にデータを収集することができる。
【００９８】
また、基地局２０と移動局１０の通信上の関係は、位置検知装置３０の機能を複数の基地局２０の少なくともいずれかに内蔵して、基地局２０が移動局１０との観測距離を収集し、基地局２０が移動局１０の位置を推定するようにしてもよいし、逆に移動局１０が基地局２０との観測距離を収集し、移動局１０が自機の位置を推定するようにしてもよい。
【００９９】
通信方法は、電波に限られず、例えば赤外線や超音波、画像処理等を用いて基地局２０と移動局１０間の観測距離を求めてもよい。
【０１００】
時刻ｓにおけるステップ１の移動局の推定初期位置４０は、時刻（ｓ−１）で求まる推定位置としてもよいし、推定位置にランダム誤差を加えた座標としてもよい。
【０１０１】
なお、上記実施形態において、各部の動作は互いに関連しており、互いの関連を考慮しながら、一連の動作及び一連の処理として置き換えることができる。これにより、位置推定装置の実施形態を、位置推定方法の実施形態とすることができる。また、位置推定装置の実施形態を、これらが有する機能をコンピュータに実現させるためのプログラムの実施形態とすることができる。
【０１０２】
これにより、位置が既知である複数の基地局と位置推定対象である移動局との一方向又は双方向通信により取得した所定の信号に基づき、距離計測部により前記複数の基地局と前記移動局との観測距離を計測し、該移動局の位置を推定する位置推定方法であって、計測された前記複数の基地局の観測距離に基づき、観測座標決定部により予め定められた移動局の推定初期位置の方向に各基地局の観測座標を決定する観測座標決定ステップと、前記観測座標決定ステップにより決定された前記複数の基地局の観測座標から移動局の代表位置を前記移動局の推定位置として算出する推定位置算出ステップと、を含むことを特徴とする位置推定方法を提供することができる。
【０１０３】
また、これにより、位置が既知である複数の基地局と位置推定対象である移動局との一方向又は双方向通信により取得した所定の信号に基づき、距離計測部により前記複数の基地局と前記移動局との観測距離を計測し、該移動局の位置を推定する処理をコンピュータに実現させるためのプログラムであって、計測された前記複数の基地局の観測距離に基づき、観測座標決定部により予め定められた移動局の推定初期位置の方向に各基地局の観測座標を決定する観測座標決定処理と、前記観測座標決定処理により決定された前記複数の基地局の観測座標から移動局の代表位置を前記移動局の推定位置として算出する推定位置算出処理と、をコンピュータに実行させるためのプログラムを提供することができる。
【０１０４】
以上、添付図面を参照しながら本発明の好適な実施形態について詳細に説明したが、本発明はかかる例に限定されない。本発明の属する技術の分野における通常の知識を有する者であれば、特許請求の範囲に記載された技術的思想の範疇内において、各種の変更例または修正例に想到し得ることは明らかであり、これらについても、当然に本発明の技術的範囲に属するものと了解される。
【符号の説明】
【０１０５】
１０ 移動局
２０ 基地局
３０ 位置検知装置
３２ 距離計測部
３４ 観測座標決定部
３６ 推定位置算出部
３８ 繰返し判定部
４０ 移動局の推定初期位置
５０ 移動局の推定位置
１００ 位置推定システム
【特許請求の範囲】
【請求項１】
位置が既知である複数の基地局と位置推定対象である移動局との一方向又は双方向通信により取得した所定の信号に基づき、前記複数の基地局と前記移動局との観測距離を計測し、該移動局の位置を推定する位置推定装置であって、
計測された前記複数の基地局の観測距離に基づき、予め定められた移動局の推定初期位置の方向に各基地局の観測座標を決定する観測座標決定部と、
前記観測座標決定部により決定された前記複数の基地局の観測座標から移動局の代表位置を前記移動局の推定位置として算出する推定位置算出部と、を備えることを特徴とする位置推定装置。
【請求項２】
前記推定位置算出部により算出された移動局の推定位置と前記移動局の推定初期位置との変化量が予め定められた停止条件を満たすまで、前記観測座標決定部及び前記推定位置算出部の動作を繰り返す繰返し判定部をさらに備えることを特徴とする請求項１に記載の位置推定装置。
【請求項３】
前記繰返し判定部は、前記推定位置算出部により算出された移動局の推定位置と前記移動局の推定初期位置との変化量が予め定められた停止条件を満たしたとき、前記移動局の推定位置を前記移動局の位置と判定する請求項１又は２に記載の位置推定装置。
【請求項４】
前記算出された移動局の推定位置を前記移動局の推定初期位置として再帰的に用いて、上記繰返し動作を繰り返す請求項２または３に記載の位置推定装置。
【請求項５】
前記推定位置算出部は、前記複数の基地局の観測座標の中央値を前記移動局の代表位置として算出する請求項１〜４のいずれか一項に記載の位置推定装置。
【請求項６】
前記推定位置算出部は、前記複数の基地局の観測座標の平均値を前記移動局の代表位置として算出する請求項１〜４のいずれか一項に記載の位置推定装置。
【請求項７】
前記推定位置算出部は、前記複数の基地局の観測座標の重み付き平均値を前記移動局の代表位置として算出する請求項１〜４のいずれか一項に記載の位置推定装置。
【請求項８】
前記推定位置算出部は、前記重み付き平均値の重みを前記観測距離の長さに対応して決定する請求項７に記載の位置推定装置。
【請求項９】
前記予め定められた移動局の推定初期位置は、前記複数の基地局に囲まれた範囲内の位置、又は前記複数の基地局から選択された一基地局から所定の範囲内の位置をデフォルト値として最初の繰返し処理に使用され、２回目以降の繰返し処理では、前回の繰返し処理で算出された移動局の推定位置を今回の繰返し処理の前記移動局の推定初期位置として再帰的に用いる請求項１〜８のいずれか一項に記載の位置推定装置。
【請求項１０】
前記観測座標決定部は、前記各基地局と前記移動局の推定初期位置とを結ぶ直線上であって、前記計測された基地局毎の観測距離だけ各基地局から離れた位置に前記各基地局の観測座標を決定する請求項１〜９のいずれか一項に記載の位置推定装置。
【請求項１１】
前記位置推定装置は、前記各基地局、前記移動局、前記複数の基地局及び前記移動局と通信可能なサーバ、又は前記移動局の位置を検知する位置検知装置のいずれかであることを特徴とする請求項１〜１０のいずれか一項に記載の位置推定装置。
【請求項１２】
位置が既知である複数の基地局と位置推定対象である移動局との一方向又は双方向通信により取得した所定の信号に基づき、距離計測部により前記複数の基地局と前記移動局との観測距離を計測し、該移動局の位置を推定する位置推定方法であって、
計測された前記複数の基地局の観測距離に基づき、観測座標決定部により予め定められた移動局の推定初期位置の方向に各基地局の観測座標を決定する観測座標決定ステップと、
前記観測座標決定ステップにより決定された前記複数の基地局の観測座標から移動局の代表位置を前記移動局の推定位置として算出する推定位置算出ステップと、を含むことを特徴とする位置推定方法。
【請求項１３】
位置が既知である複数の基地局と位置推定対象である移動局との一方向又は双方向通信により取得した所定の信号に基づき、距離計測部により前記複数の基地局と前記移動局との観測距離を計測し、該移動局の位置を推定する処理をコンピュータに実現させるためのプログラムであって、
計測された前記複数の基地局の観測距離に基づき、観測座標決定部により予め定められた移動局の推定初期位置の方向に各基地局の観測座標を決定する観測座標決定処理と、
前記観測座標決定処理により決定された前記複数の基地局の観測座標から移動局の代表位置を前記移動局の推定位置として算出する推定位置算出処理と、をコンピュータに実行させるためのプログラム。
【請求項１】
位置が既知である複数の基地局と位置推定対象である移動局との一方向又は双方向通信により取得した所定の信号に基づき、前記複数の基地局と前記移動局との観測距離を計測し、該移動局の位置を推定する位置推定装置であって、
計測された前記複数の基地局の観測距離に基づき、予め定められた移動局の推定初期位置の方向に各基地局の観測座標を決定する観測座標決定部と、
前記観測座標決定部により決定された前記複数の基地局の観測座標から移動局の代表位置を前記移動局の推定位置として算出する推定位置算出部と、を備えることを特徴とする位置推定装置。
【請求項２】
前記推定位置算出部により算出された移動局の推定位置と前記移動局の推定初期位置との変化量が予め定められた停止条件を満たすまで、前記観測座標決定部及び前記推定位置算出部の動作を繰り返す繰返し判定部をさらに備えることを特徴とする請求項１に記載の位置推定装置。
【請求項３】
前記繰返し判定部は、前記推定位置算出部により算出された移動局の推定位置と前記移動局の推定初期位置との変化量が予め定められた停止条件を満たしたとき、前記移動局の推定位置を前記移動局の位置と判定する請求項１又は２に記載の位置推定装置。
【請求項４】
前記算出された移動局の推定位置を前記移動局の推定初期位置として再帰的に用いて、上記繰返し動作を繰り返す請求項２または３に記載の位置推定装置。
【請求項５】
前記推定位置算出部は、前記複数の基地局の観測座標の中央値を前記移動局の代表位置として算出する請求項１〜４のいずれか一項に記載の位置推定装置。
【請求項６】
前記推定位置算出部は、前記複数の基地局の観測座標の平均値を前記移動局の代表位置として算出する請求項１〜４のいずれか一項に記載の位置推定装置。
【請求項７】
前記推定位置算出部は、前記複数の基地局の観測座標の重み付き平均値を前記移動局の代表位置として算出する請求項１〜４のいずれか一項に記載の位置推定装置。
【請求項８】
前記推定位置算出部は、前記重み付き平均値の重みを前記観測距離の長さに対応して決定する請求項７に記載の位置推定装置。
【請求項９】
前記予め定められた移動局の推定初期位置は、前記複数の基地局に囲まれた範囲内の位置、又は前記複数の基地局から選択された一基地局から所定の範囲内の位置をデフォルト値として最初の繰返し処理に使用され、２回目以降の繰返し処理では、前回の繰返し処理で算出された移動局の推定位置を今回の繰返し処理の前記移動局の推定初期位置として再帰的に用いる請求項１〜８のいずれか一項に記載の位置推定装置。
【請求項１０】
前記観測座標決定部は、前記各基地局と前記移動局の推定初期位置とを結ぶ直線上であって、前記計測された基地局毎の観測距離だけ各基地局から離れた位置に前記各基地局の観測座標を決定する請求項１〜９のいずれか一項に記載の位置推定装置。
【請求項１１】
前記位置推定装置は、前記各基地局、前記移動局、前記複数の基地局及び前記移動局と通信可能なサーバ、又は前記移動局の位置を検知する位置検知装置のいずれかであることを特徴とする請求項１〜１０のいずれか一項に記載の位置推定装置。
【請求項１２】
位置が既知である複数の基地局と位置推定対象である移動局との一方向又は双方向通信により取得した所定の信号に基づき、距離計測部により前記複数の基地局と前記移動局との観測距離を計測し、該移動局の位置を推定する位置推定方法であって、
計測された前記複数の基地局の観測距離に基づき、観測座標決定部により予め定められた移動局の推定初期位置の方向に各基地局の観測座標を決定する観測座標決定ステップと、
前記観測座標決定ステップにより決定された前記複数の基地局の観測座標から移動局の代表位置を前記移動局の推定位置として算出する推定位置算出ステップと、を含むことを特徴とする位置推定方法。
【請求項１３】
位置が既知である複数の基地局と位置推定対象である移動局との一方向又は双方向通信により取得した所定の信号に基づき、距離計測部により前記複数の基地局と前記移動局との観測距離を計測し、該移動局の位置を推定する処理をコンピュータに実現させるためのプログラムであって、
計測された前記複数の基地局の観測距離に基づき、観測座標決定部により予め定められた移動局の推定初期位置の方向に各基地局の観測座標を決定する観測座標決定処理と、
前記観測座標決定処理により決定された前記複数の基地局の観測座標から移動局の代表位置を前記移動局の推定位置として算出する推定位置算出処理と、をコンピュータに実行させるためのプログラム。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【公開番号】特開２０１１−２１４９２０（Ｐ２０１１−２１４９２０Ａ）
【公開日】平成２３年１０月２７日（２０１１．１０．２７）
【国際特許分類】
【出願番号】特願２０１０−８１８１３（Ｐ２０１０−８１８１３）
【出願日】平成２２年３月３１日（２０１０．３．３１）
【新規性喪失の例外の表示】特許法第３０条第１項適用申請有り 発行者：日本計算機統計学会 刊行物名：日本計算機統計学会 第２３回シンポジウム論文集 発行年月日：２００９年１１月７日
【出願人】（００００００２９５）沖電気工業株式会社 (6,645)
【出願人】（５９９０１１６８７）学校法人 中央大学 (110)
【Ｆターム（参考）】
【公開日】平成２３年１０月２７日（２０１１．１０．２７）
【国際特許分類】
【出願日】平成２２年３月３１日（２０１０．３．３１）
【新規性喪失の例外の表示】特許法第３０条第１項適用申請有り 発行者：日本計算機統計学会 刊行物名：日本計算機統計学会 第２３回シンポジウム論文集 発行年月日：２００９年１１月７日
【出願人】（００００００２９５）沖電気工業株式会社 (6,645)
【出願人】（５９９０１１６８７）学校法人 中央大学 (110)
【Ｆターム（参考）】
[ Back to top ]