JP5039063B2 - 携帯電子機器の状態を追跡する方法 - Google Patents

携帯電子機器の状態を追跡する方法 Download PDF

Info

Publication number
JP5039063B2
JP5039063B2 JP2008557597A JP2008557597A JP5039063B2 JP 5039063 B2 JP5039063 B2 JP 5039063B2 JP 2008557597 A JP2008557597 A JP 2008557597A JP 2008557597 A JP2008557597 A JP 2008557597A JP 5039063 B2 JP5039063 B2 JP 5039063B2
Authority
JP
Japan
Prior art keywords
probability value
lattice
grid
electronic device
state
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2008557597A
Other languages
English (en)
Other versions
JP2009529128A (ja
Inventor
ニーロ シロラ
シモ アリ‐ロイッティー
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nokia Oyj
Original Assignee
Nokia Oyj
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 Nokia Oyj filed Critical Nokia Oyj
Publication of JP2009529128A publication Critical patent/JP2009529128A/ja
Application granted granted Critical
Publication of JP5039063B2 publication Critical patent/JP5039063B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S5/00Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
    • G01S5/02Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using radio waves
    • G01S5/0294Trajectory determination or predictive filtering, e.g. target tracking or Kalman filtering
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W4/00Services specially adapted for wireless communication networks; Facilities therefor
    • H04W4/02Services making use of location information
    • H04W4/029Location-based management or tracking services

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Mobile Radio Communication Systems (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Description

本発明は、携帯電子機器の状態を追跡する方法と、その方法を実行するように構成される処理装置を備える携帯電子機器とに関する。
発明の背景
パーソナルな測位における課題の1つとして、測定源の数が少なく、測定源が大きな誤差を含むものである可能性があり、測定源の分布も適当ではない状況、具体的には屋内または市街地において、正確な位置情報を提供することが挙げられる。この課題を解決するには、これらの複数の測定源の融合から生じる非線形フィルタ方程式をうまく数値解析しなければならない。この場合、測定毎に最大量の情報を抽出することが有利である。
GPSなどの衛星ベースのシステムの動作は、屋内において高感度モードで使用されても予測不可能である。セルラーネットワーク、WLAN、またはブルートゥースなどのローカル無線ネットワークにより、ある程度の測位能力がもたらされるが、GPSに比べて精度が劣る。携帯電子機器のその他の可能な部品として、加速度計、気圧計、またはデジタルコンパスなどの搭載センサが挙げられる。
種々の測定源を組み合わせることは、異なる誤差特性、予想不可能な歪み、測定における定誤差、強い非線形性、複雑な時間依存性、およびデータの欠落により、困難である。一般的な方法で全事例をモデル化することは容易ではなく、また、モデルを正確に解くことは容易でないことは言うまでもない。正確なモデルを使用しても、一般的なカルマンフィルタおよびその非線形拡張は、突然失敗する可能性がある。
発明の摘要
本発明の第1の側面によると、携帯電子機器の状態を追跡する方法であって、
(i)携帯電子機器の状態を、複数のセルを含む格子(グリッド:grid)を用いて表すことであって、セルの各々が、1つ以上の状態変数により規定される状態空間における領域を表すと共に、携帯電子機器の状態が状態空間におけるその領域内にある確率値を有し、格子は、既定の閾値を超える確率値を有するセルのみを含むように境界付けられることと、
(ii)1つ以上の状態変数の値を示す測定信号を入手することと、
(iii)測定信号に基づいて格子の確率値を更新すること、および格子を再び境界付けること、
を、反復的に実行することを含む方法が提供される。
携帯電子機器の状態は、位置、速度、加速度、時計誤差を、単独でまたは組み合わせて含むがそれだけに限定されない状態変数を含んでもよい。
状態空間の次元の数は、状態変数の次元の総数に相当し、各セルは、状態空間の数に等しい次元の数を含む。例えば、状態空間およびセルは、3次元における位置および3次元における速度を表すために6次元であってもよい。
本発明は、全情報が、位置の解を一意にはもたらさない、または多数の解をもたらす測定ジオメトリから保持され、かつ事後分布の形状が、そのいかなるピークも降下させずに保持される手段を提供する。
本発明は、粒子フィルタ(particle filter)よりも高速に実行され、位置指紋方法(location fingerprint method)よりも一般的であり、完全な位置分布を伝搬させることから、不確定または多数解決のシステムを正確に表すことができる。
格子は、本方法の計算負荷に関する問題を緩和するために、一様な間隔を有する平行6面体格子であることが好ましい。
反復されるステップ(i)のうち1回以上は、測定信号を入手することと、測定信号に基づいて格子を形成することを含んでもよい。位置推定は、例えば、領域、領域差、または平面測定のうちのいずれかを含む測定信号から導かれてもよい。領域測定は、セルラー基地局、WLANまたはブルートゥース送受信機、および/または音響センサから得られてもよく、例えば、時間遅延、往復(round-trip)、または信号強度測定の形式であってもよい。位置を示すその他の種類の測定信号には、到来角測定、最大最小領域窓、基地セクター情報、ならびに搭載気圧計およびデジタルコンパスが含まれるが、それだけに限定されない。速度は、デルタレンジおよび方位測定を使用して測定可能である。加速度は、オンボードの加速度計を使用して測定可能である。
測定信号の詳細は、本発明にとって重要でないことが理解されたい。但し、該当する1つまたは複数の状態変数の1つまたは複数の値は、単独での使用またはその他の測定信号と併用する使用か否かに関わらず、測定信号を使用して入手または推定可能である。
反復されるステップ(i)のうち2回目以降の1回以上は、先に反復されたステップ(iii)において更新された格子を使用することを含んでもよい。
1回以上のステップ(i)は、(iv)更新された格子と、携帯電子機器の移動モデルとに基づいて、後続の反復における格子の確率値を予測して、予測格子を入手すること、を含んでもよい。
ステップ(iv)は、後続の反復におけるステップ(i)において、格子を再び境界付けること、および/または予測格子を使用することを含んでもよい。
格子を再び境界付けることは、既定の閾値未満の確率値を有するセルを除外するように格子の境界を移動すること、および/または既定の閾値を超える確率値を有するセルを含むように境界を移動すること、を含んでもよい。場合により、セルは、境界外に存在するが、その確率値はゼロである。いかなる過程の後であっても、そのようなセルの確率値が増加する場合に、境界は、既定の閾値を超える確率値を有するいかなるそのようなセルも含むように移動される。場合により、セルは、境界外には存在しない。確率値の変更をもたらすいかなる過程の後であっても、本発明は、境界外の暫定セルを規定することと、これらのセルの確率値を計算することと、既定の閾値を超える確率値を有するいかなるこのようなセルも含むように境界を移動すること、を含んでもよい。
移動モデルは、方法の計算負荷に関する問題を緩和するために、線形であることが好ましい。
反復的な実行の1回又はそれ以上において、更新された格子に基づいて、携帯電子機器の状態に関する期待値および分散を計算すること(v)を含んでもよい。
本発明の第2の側面によると、本発明の第1の側面の方法を実行するように構成される処理装置を備える携帯電子機器が提供される。
また、本発明は、上記の方法を実行するように構成されるコンピュータプログラムと、本発明の携帯電子機器が動作するシステムとを含む。本発明は、組み合わせが具体的に言及(または請求)されているか否かに関わらず、種々の組み合わせによるその1つ以上の側面および実施形態を包含する。
発明の詳細な説明
本発明をより容易に理解するために、添付の図面を参照して、ほんの一例として以下に説明される。以下の図面において、同一の数字は、その特徴が属する実施形態に関わらず、同一の特徴を言及している。
図6は、本発明に従う携帯電子機器の一部に関する概略図である。携帯電子機器は、メモリ102と、複数の送受信機モジュール104a〜dと通信するプロセッサ100とを備える。メモリ102は、後述する本発明の方法を実行するためのソフトウェアコード部分を含む。送受信機モジュール104a〜dは、GPSシステム、セルラー基地局、WLAN送信機、およびブルートゥース送信機からそれぞれ測定信号を受信することができる。送受信機モジュール104a〜dは、便宜上このように言及されるが、実際は、モジュール104a〜dによっては、送信する能力を含まないものもあれば、測定信号を受信するために送信することを要求されるものもあってもよい。
本発明の方法について、以下の用語を使用して説明していく。
Figure 0005039063
〔ステップ(i)〕
本ステップにおいて、複数のセルを含む格子を使用して、携帯電子機器の状態xkを示す。各セルは、状態空間における領域を示し、携帯電子機器の状態xkがその領域内にある確率値を有する。
例えば、状態xkは、3次元状態変数の位置rkやvkを含んでもよい。すなわち
Figure 0005039063
である。次に、各セルは、状態空間の6次元領域を表すために6次元である。しかしながら、説明目的のために、図1から図5は、状態xkが2次元位置rkを含む簡単な例を示す。
図1は、95%、66%、および50%の信頼領域を有する携帯電子機器の位置に関する事前分布(prior distribution )10を示す。事前分布10は、(i)以前に繰り返し実行した本方法から得られる事後分布(posterior distribution)、(ii)携帯電子機器の移動を予測するために、移動モデルを使用して事後分布を変更することによって得られる予測事前分布、および(iii)送受信機モジュール104a〜dのいずれかを介して得られる測定信号、のいくつかのソースのうちの1つから導かれてもよい。上記ソースのいくつかが利用可能である場合、組み合わせて使用可能である。
事前分布10は、サイズおよび形状が均一である多数の2次元セル14を含む事前格子12を使用して近似化される。各セル14は、地球表面における領域を表し、携帯電子機器がその領域内に位置する事前確率値を有する。
図1は、説明目的のためだけに、格子12における信頼領域により示される事前分布10を示すことを理解されたい。事前分布10は連続的であるように図1に示されるが、事前格子12の確率値により近似化される形式でのみ携帯電子機器のメモリ102に存在する。図1のセル14の相対的陰影は、確率値を表し、つまり、セル14の陰影が濃くなるほど、セル14の確率値が高くなることを示す。
事前格子12は境界16を含み、その境界内の全てのセル14が既定の閾値を超える事前確率値を有する。重要な領域Sを表すべく、事前分布の格子近似が、境界16によって定められる。領域Sは、
Figure 0005039063
の中で(単純に結合された)領域であり、その中の事前確率値は無視できない。
図1の例では、事前分布の99%が事前格子12に含まれるように閾値を設定した。しかしながら、閾値レベルは、設計上の選択の問題であることを理解されたい。閾値を低くすると、境界16内に含まれるセル14が増加するため、計算負荷の増加を伴うが、近似の精度が増す。反対に、閾値を高くすると、境界内に含まれるセル14が減少するため、近似の精度は低下するが、計算負荷は減少する。
境界16内のセル14だけが、携帯電子機器のメモリ102に書き込まれる。変形例において、メモリ102は、境界外のセルをも書き込んでもよいが、これらの事前確率値はゼロに設定される。いずれの場合においても、事前分布10の事前格子12により、近似領域は、計算し易くするように境界16によって範囲が定められる。
事前格子12におけるセル14の数は、設計上の選択の問題であり、計算負荷と精度をバランスさせねばならない。極端な選択すると、多数の小さいセルを生成することになる。この場合、本発明に従うセル14の中心の密度値に相当する密度値のみを使用する既知の点質量フィルタ(point-mass filter)の場合のように、近似表現は、確率値が準最適であったとしても正確なものに近づく。別の極端な選択をすると、少数の大きなセル14を使用することになる。従って、事前確率値を可能な限り正確に計算することは有利である。事前分布10の構造の多くは、大きいセル14で近似される場合に失われる。セル14は、事前分布10の最も微細な特徴よりもさらに小さくはすべきではない。
本ステップにおいて、時間インデックスkは、k=1に設定される。

〔ステップ(ii)〕
このステップにおいて、1つ以上の状態変数の値を示す送受信機モジュール104a〜dを介して測定信号を入手する。
図2は、送受信機モジュール104a〜dを介して得られる測定信号に基づき、携帯電子機器の推定位置を示す測定位置を示す測定尤度関数(measurement likelihood function)18を示す。図示されるように、測定尤度関数は、95%、66%、および50%信頼領域を含む。
図3は、事前分布10と同様の要領で、事前格子12を使用して近似化される測定尤度関数18を示す。図3において、携帯電子機器がそれぞれのセル14によって表される領域内に位置する確率を再び表す確率値は、測定尤度関数18のみに基づく。
図3において、測定尤度関数18は、事前分布10と同一の事前格子12を使用して近似化されるが、異なるサイズ、形状、および/または配向の格子を使用して測定尤度関数を表してもよいことを理解されたい。しかしながら、これにより、さらなる計算が必要になり得る。
測定信号には、位置推定を入手するための、例えば、領域、領域差、または平面測定のうちのいずれかが含まれてもよい。
真の位置をrと仮定して、位置sにおける局までの領域測定は、
Figure 0005039063
として記述可能である。関連する測定誤差νは、正規である必要はなく、現状に一致する経験的に決定される分布によって示される。
GPSシステムから得られる、バイアスが含まれる領域測定は、領域差として処理される。局のうちの1つが基準局として選択され、それに対する全ての差が計算される。基準局がs0にある場合、領域差測定は、
Figure 0005039063
となる。
最後に、平面測定は、
Figure 0005039063
である。ここでuは単位ベクトルである。
このステップ中、利用可能な全ての測定はベクトルykに格納され、対応する測定式は、ベクトル関数hk(x)に格納される。個の
例として、領域差測定がnd回、領域測定がnr回、および平面測定がnp回の場合を考える。測定ベクトルは、
Figure 0005039063
であり、測定モデルは、
Figure 0005039063
である。
全ての測定誤差が正規分布を有する場合、測定尤度関数18は、
Figure 0005039063
であり、ここでΣは測定誤差νの共分散行列である。

〔ステップ(iii)〕
本ステップにおいて、測定信号に基づいて事前格子12の事前確率値を更新することによって事後格子12'を生成し、事後格子12'の境界を再定義する。
図4は、事前分布10および測定尤度関数18から導かれた事後分布20を示す。事後分布20は、セル14'を有する事後格子12'によって表され、各セル14'は、事後分布20を近似化するために、事後確率値を有する。
事後格子12'の事後確率値は、各セル14における事前確率値を、上の尤度関数をセル14に亘って積分することによって求められる、セル14の全尤度で乗算することによって求められる。
Figure 0005039063
ここで
Figure 0005039063
である。
事後格子12'は、既定の閾値を超える確率値を有するセル14'のみを含むように新たな境界16'で再び境界付けられる。前境界16内のセル14が既定の閾値未満の事後確率値を有する場合、新しい境界16'は、そのようなセル14を除外するように配置される。

〔ステップ(iv)〕
本ステップはオプションである。本ステップにおいて、セル14"を有する予測格子12"は、現在の試行における事後格子12'および携帯電子機器の移動モデルに基づいて得られる。予測格子12"は、後続の試行における事前格子12を形成する。
図5は、移動モデル適用後の予測事前分布10'を示す。この移動モデルは、一定速度正規分布加速度モデル(constant-velocity-normally-distributed-acceleration model)またはシンガーモデル(singer model)などの、適切な任意のモデルである。予測事前分布10'は、本方法の後続の試行における事前分布10を形成する。
セル14"の予測事前確率値は、以下の式に従って、そのセル14"への移行に関する確率によって重み付けされた現在の試行における全確率値を総和することによって求められる。予測事前分布10'は、全確率値の計算後に正規化される。
Figure 0005039063
各セル14"の体積を
Figure 0005039063
で示すことによって、予測事前確率値
Figure 0005039063
は、セル14"上の予測pdfを積分することによって計算される。
Figure 0005039063
Figure 0005039063
をその格子近似に置き換えると、
Figure 0005039063
となり、ここで
Figure 0005039063
は、(k-1)番目の格子のj番目のセルからk番目の格子i番目のセルへの移行確率である。
計算の効率を上げるために、例えば
Figure 0005039063
などの線形移動モデルを使用する。予測格子12"は、移動モデルを事後格子12'に適用することによって形成される。旧格子のi番目のセルと、新格子のj番目のセルとの間の移行確率は、差(i - j)にのみ依存し、
Figure 0005039063
と記述可能である。予測事前確率値は、
Figure 0005039063
となり、これにより、d次元の離散線形畳み込みの計算が速くなる。
次に、移行確率
Figure 0005039063
は、
Figure 0005039063
となり、ここで
Figure 0005039063
である。これは、
Figure 0005039063
と簡略化することが可能である。
ここで
Figure 0005039063
は、修正された過程雑音(process noise)確率密度関数の省略表現である。具体的には、
Figure 0005039063
である場合、
Figure 0005039063
となる。
積分は、ハイパーボックスにおける単なる多正規確率であり、数値的に計算可能である。
wkが非ガウス型である場合、移行確率
Figure 0005039063
は、累積分布を使用して計算可能である。
移動モデルの適用後、既定の閾値を超える予測確率値を有する、前境界16'内外の全セル14"を含む境界16"になるように境界16'を移動する。
セル14"が前境界16'外に存在し、かつ事前および事後確率値がゼロであった場合、新境界16"は、既定の閾値を超える予測確率値を有するこれらのセル14"のいずれも含むように配置される。
前境界16'外にセルが存在しなかった変形例において、既存のセル14'のものと対応するサイズ、形状、および幾何学的配置を有する新しいセル14"が加えられ、新境界16"は、新しいセル14"を含むように配置される。新しいセル14"のサイズ、形状、および/または幾何学的配置は、既存の予測格子12"のものと一致する必要はないことを理解されたい。このような存在しないセルの確率値が、既定の閾値を上回るか否かに関しては、前境界16'外の多数の暫定セル14"を規定することによって、ならびに暫定セル14"の予測確率値を計算することによって、容易に判断可能である。このような暫定セル14"は、継続的に規定され、その確率値は、確率値が閾値未満である暫定セル14"のシェルを生成する点に達するまで計算される。次に、新境界16"は、閾値を超える確率値を有する暫定セル14"を全て含むように配置される。

〔ステップ(v)〕
本ステップにおいて、携帯電子機器の状態に関する期待値および分散を、以下の数式に従って、予測格子12"に基づき計算する。
Figure 0005039063
Figure 0005039063
最終的に、時間インデックスkは増分され、方法は、ステップ(i)から反復される。
本発明は、格子−質量アプローチ(grid-mass approach)を使用して、携帯電子機器の状態を追跡する方法を実行することを理解されたい。
本発明による方法の簡単な例を示す。 本発明による方法の簡単な例を示す。 本発明による方法の簡単な例を示す。 本発明による方法の簡単な例を示す。 本発明による方法の簡単な例を示す。 本発明に従う装置をしめる。

Claims (17)

  1. 携帯電子機器の状態を追跡する方法であって、
    (i)前記携帯電子機器の前記状態を、複数のセルを含む格子を用いて表すことであって、前記セルの各々は、1つ以上の状態変数により規定される状態空間における領域を表すと共に、前記携帯電子機器の前記状態が状態空間におけるその領域内にある第1の確率値を有し、前記格子は、既定の閾値を超える確率値を有するセルのみを含むように境界付けられることと、
    (ii)1つ以上の状態変数の値を示す測定信号を入手することと、
    (iii)第2の確率値を有するセルを含む更新された格子を入手するため、前記格子の前記第1の確率値を更新することおよび前記格子を再び境界付けることであって、前記更新された格子のセルの前記第2の確率値は、測定尤度関数をセルに亘って積分したもので前記第1の確率値を乗算することによって求められる、前記更新することおよび再び境界付けることと、
    (iv)後続の試行におけるステップ(i)で使用される予測格子を入手するため、前記更新された格子と、前記携帯電子機器の移動モデルとに基づいて、前記格子の第3の確率値を予測することであって、前記予測された格子のセルの前記第3の確率値は、そのセルへの移行に関する確率によって重み付けされた現在の試行における全ての前記第2の確率値を総和することによって求められる、前記予測することと、
    を反復的に実行することを含む、方法。
  2. 前記格子は、一様な間隔を有する平行6面体格子である、請求項1に記載の方法。
  3. 反復されるステップ(i)のうち1回以上は、測定信号を入手することと、前記測定信号に基づいて前記格子を形成することを含む、請求項1又は2に記載の方法。
  4. ステップ(iv)は、前記格子を再び境界付けることを含む、請求項1から3のいずれかに記載の方法。
  5. 前記移動モデルは線形である、請求項1から4のいずれかに記載の方法。
  6. 前記反復的な実行の1回又はそれ以上において、
    (v)前記予測された格子に基づいて、前記携帯電子機器の前記状態に関する期待値および分散を計算すること、
    を含む、請求項1から5のいずれかに記載の方法。
  7. 処理装置を備える携帯電子機器であって、前記処理装置は、
    (i)前記携帯電子機器の状態を、複数のセルを含む格子を用いて表し、前記セルの各々は、1つ以上の状態変数により規定される状態空間における領域を表すと共に、前記携帯電子機器の前記状態が状態空間におけるその領域内にある第1の確率値を有し、前記格子は、既定の閾値を超える確率値を有するセルのみを含むように境界付けられ、
    (ii)1つ以上の状態変数の値を示す測定信号を入手し、
    (iii)第2の確率値を有するセルを含む更新された格子を入手するため、前記格子の前記第1の確率値を更新し、前記格子を再び境界付け、前記更新された格子のセルの前記第2の確率値は、測定尤度関数をセルに亘って積分したもので前記第1の確率値を乗算することによって求められ、
    (iv)後続の試行におけるステップ(i)で使用される予測格子を入手するため、前記更新された格子と、前記携帯電子機器の移動モデルとに基づいて、前記格子の第3の確率値を予測し、前記予測された格子のセルの前記第3の確率値は、そのセルへの移行に関する確率によって重み付けされた現在の試行における全ての前記第2の確率値を総和することによって求められ
    前記(i)〜(iv)の処理を反復するように構成される、携帯電子機器。
  8. 携帯電子機器の状態を追跡する際に使用するコンピュータプログラムであって、
    (i)前記携帯電子機器の前記状態を、複数のセルを含む格子を用いて表すことであって、前記セルの各々は、1つ以上の状態変数により規定される状態空間における領域を表すと共に、前記携帯電子機器の前記状態が状態空間におけるその領域内にある第1の確率値を有し、前記格子は、既定の閾値を超える確率値を有するセルのみを含むように境界付けられることと、
    (ii)1つ以上の状態変数の値を示す測定信号を入手することと、
    (iii)第2の確率値を有するセルを含む更新された格子を入手するため、前記格子の前記第1の確率値を更新することおよび前記格子を再び境界付けることであって、前記更新された格子のセルの前記第2の確率値は、測定尤度関数をセルに亘って積分したもので前記第1の確率値を乗算することによって求められる、前記更新することおよび再び境界付けることと、
    (iv)後続の試行におけるステップ(i)で使用される予測格子を入手するため、前記更新された格子と、前記携帯電子機器の移動モデルとに基づいて、前記格子の第3の確率値を予測することであって、前記予測された格子のセルの前記第3の確率値は、そのセルへの移行に関する確率によって重み付けされた現在の試行における全ての前記第2の確率値を総和することによって求められる、前記予測することと、
    を、反復的に実行するように構成されるコンピュータコードを含む、コンピュータプログラム。
  9. 処理装置を具備する携帯電子機器を有するシステムであって、前記処理装置は、
    (i)前記携帯電子機器の状態を複数のセルを含む格子を用いて表し、前記セルの各々は、1つ以上の状態変数により規定される状態空間における領域を表すと共に、前記携帯電子機器の前記状態が状態空間におけるその領域内にある第1の確率値を有し、前記格子は、既定の閾値を超える確率値を有するセルのみを含むように境界付けられ、
    (ii)1つ以上の状態変数の値を示す測定信号を入手し、
    (iii)第2の確率値を有するセルを含む更新された格子を入手するため、前記格子の前記第1の確率値を更新し、前記格子を再び境界付け、前記更新された格子のセルの前記第2の確率値は、測定尤度関数をセルに亘って積分したもので前記第1の確率値を乗算することによって求められ、
    (iv)後続の試行におけるステップ(i)で使用される予測格子を入手するため、前記更新された格子と、前記携帯電子機器の移動モデルとに基づいて、前記格子の第3の確率値を予測し、前記予測された格子のセルの前記第3の確率値は、そのセルへの移行に関する確率によって重み付けされた現在の試行における全ての前記第2の確率値を総和することによって求められ
    前記(i)〜(iv)の処理を反復するように構成される、システム。
  10. 携帯電子機器の状態を複数のセルを含む格子を用いて表すことを反復的に行いうる手段であって、前記セルの各々は、1つ以上の状態変数により規定される状態空間における領域を表すと共に、前記携帯電子機器の前記状態が状態空間におけるその領域内にある第1の確率値を有し、前記格子は、既定の閾値を超える確率値を有するセルのみを含むように境界付けられる、手段と、
    1つ以上の状態変数の値を示す測定信号を入手することを反復的に行いうる手段と、
    第2の確率値を有するセルを含む更新された格子を入手するため、前記格子の前記第1の確率値を更新し、前記格子を再び境界付けることを反復的に行いうる手段であって、前記更新された格子のセルの前記第2の確率値は、測定尤度関数をセルに亘って積分したもので前記第1の確率値を乗算することによって求められる、前記更新および再び境界付けることを反復的に行いうる手段と、
    後続の試行における前記携帯電子機器の状態の表現で使用される予測格子を入手するため、前記更新された格子と、前記携帯電子機器の移動モデルとに基づいて、前記格子の第3の確率値を予測することを反復的に行いうる手段であって、前記予測された格子のセルの前記第3の確率値は、そのセルへの移行に関する確率によって重み付けされた現在の試行における全ての前記第2の確率値を総和することによって求められる、前記予測することを反復的に行いうる手段と、
    を備える、装置。
  11. (i)携帯電子機器の状態を複数のセルを含む格子を用いて表し、前記セルの各々は、1つ以上の状態変数により規定される状態空間における領域を表すと共に、前記携帯電子機器の前記状態が状態空間におけるその領域内にある第1の確率値を有し、前記格子は、既定の閾値を超える確率値を有するセルのみを含むように境界付けられ、
    (ii)1つ以上の状態変数の値を示す測定信号を入手し、
    (iii)第2の確率値を有するセルを含む更新された格子を入手するため、前記格子の前記第1の確率値を更新し、前記格子を再び境界付け、前記更新された格子のセルの前記第2の確率値は、測定尤度関数をセルに亘って積分したもので前記確率値を乗算することによって求められ、
    (iv)後続の試行におけるステップ(i)で使用される予測格子を入手するため、前記更新された格子と、前記携帯電子機器の移動モデルとに基づいて、前記格子の第3の確率値を予測し、前記予測された格子のセルの前記第3の確率値は、そのセルへの移行に関する確率によって重み付けされた現在の試行における全ての前記第2の確率値を総和することによって求められ
    前記(i)〜(iv)の処理を反復するように構成される装置。
  12. 前記装置は、一様な間隔を有する平行6面体格子である格子を使用して、携帯電子機器の状態を表すように構成される、請求項10又は11に記載の装置。
  13. 前記装置は、測定信号を入手することと、前記測定信号に基づいて前記格子を形成することを、1回又は複数回行うことにより、携帯電話機器の状態を表すように構成される、請求項10から12のいずれかに記載の装置。
  14. 前記装置は、前記第3の確率値を予測した後、前記格子を再び境界付けるように構成される、請求項10から13のいずれかに記載の装置。
  15. 前記装置は、携帯電子機器前記状態を表すために、線形移動モデルを使用するように構成される、請求項10から14のいずれかに記載の装置。
  16. 前記装置は、1回の試行又は複数回の反復試行により、前記予測された格子に基づいて、前記携帯電子機器の前記状態に関する期待値および分散を計算するように構成される、請求項10から15のいずれに記載の装置。
  17. 装置であって、
    (i)携帯電子機器の状態を複数のセルを含む格子を用いて表し、前記セルの各々は、1つ以上の状態変数により規定される状態空間における領域を表すと共に、前記携帯電子機器の前記状態が状態空間におけるその領域内にある第1の確率値を有し、前記格子は、既定の閾値を超える確率値を有するセルのみを含むように境界付けられ、
    (ii)1つ以上の状態変数の値を示す測定信号を入手し、
    (iii)第2の確率値を有するセルを含む更新された格子を入手するため、前記格子の前記第1の確率値を更新し、前記格子を再び境界付け、前記更新された格子のセルの前記第2の確率値は、測定尤度関数をセルに亘って積分したもので前記第1の確率値を乗算することによって求められ、
    (iv)後続の試行におけるステップ(i)で使用される予測格子を入手するため、前記更新された格子と、前記携帯電子機器の移動モデルとに基づいて、前記格子の第3の確率値を予測し、前記予測された格子のセルの前記第3の確率値は、そのセルへの移行に関する確率によって重み付けされた現在の試行における全ての前記第2の確率値を総和することによって求められ、
    前記(i)〜(iv)の処理を反復する手段を備える、装置。
JP2008557597A 2006-03-07 2006-03-07 携帯電子機器の状態を追跡する方法 Active JP5039063B2 (ja)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/EP2006/002272 WO2007101453A1 (en) 2006-03-07 2006-03-07 A method of tracking a state of a mobile electronic device

Publications (2)

Publication Number Publication Date
JP2009529128A JP2009529128A (ja) 2009-08-13
JP5039063B2 true JP5039063B2 (ja) 2012-10-03

Family

ID=37101900

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2008557597A Active JP5039063B2 (ja) 2006-03-07 2006-03-07 携帯電子機器の状態を追跡する方法

Country Status (7)

Country Link
US (1) US8213956B2 (ja)
EP (1) EP1991882B1 (ja)
JP (1) JP5039063B2 (ja)
KR (1) KR101106276B1 (ja)
CN (1) CN101438184B (ja)
BR (1) BRPI0621604B1 (ja)
WO (1) WO2007101453A1 (ja)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5094589B2 (ja) * 2008-06-25 2012-12-12 ヤフー株式会社 現在位置推定装置、方法及びシステム
FR2944603B1 (fr) * 2009-04-17 2014-07-11 Univ Troyes Technologie Systeme et procede de localisation de cible par un reseau d'emetteurs recepteurs
US20110274094A1 (en) * 2010-05-06 2011-11-10 Qualcomm Incorporated Methods and apparatus for estimating a geolocation of a wireless communications device
US9062978B2 (en) * 2011-05-31 2015-06-23 Massachusetts Institute Of Technology Tracking a body by nonlinear and non-Gaussian parametric filtering
FR2979900B1 (fr) 2011-09-12 2013-08-30 Eurocopter France Aeronef rapide a grande distance franchissable
CN103916954B (zh) * 2013-01-07 2017-11-03 华为技术有限公司 基于wlan的概率定位方法及定位装置
US20150181372A1 (en) * 2013-12-20 2015-06-25 Apple Inc. Location Accuracy Prediction
US10488505B2 (en) 2014-05-30 2019-11-26 The Boeing Company Positioning in indoor locations and other GPS-denied environments
US9929817B2 (en) 2014-11-19 2018-03-27 Oath Inc. System and method for 3D tracking for ad-hoc cross-device interaction
CN105699937B (zh) * 2014-11-27 2017-08-25 中国科学院沈阳自动化研究所 一种基于css定位技术的变电站高精度混合定位方法
US10444269B2 (en) 2017-05-26 2019-10-15 Honeywell International Inc. Apparatus and method for performing grid adaption in numerical solution of recursive bayesian estimators
US11378403B2 (en) 2019-07-26 2022-07-05 Honeywell International Inc. Apparatus and method for terrain aided navigation using inertial position

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4796191A (en) 1984-06-07 1989-01-03 Etak, Inc. Vehicle navigational system and method
EP0514887B1 (en) 1991-05-21 1997-04-16 Matsushita Electric Industrial Co., Ltd. Vehicle position detecting apparatus
US6889053B1 (en) * 1999-07-26 2005-05-03 Lucent Technologies Inc. Likelihood-based geolocation prediction algorithms for CDMA systems using pilot strength measurements
JP3943339B2 (ja) * 2000-02-25 2007-07-11 株式会社エヌ・ティ・ティ・ドコモ 移動通信システムにおける移動機の位置推定方法及びシステム
SE0004096D0 (sv) * 2000-11-08 2000-11-08 Nira Automotive Ab Positioning system
FI111901B (fi) * 2000-12-29 2003-09-30 Ekahau Oy Sijainnin arviointi langattomissa tietoliikenneverkoissa
US6542116B1 (en) * 2001-06-22 2003-04-01 Enuvis, Inc. Determining the spatio-temporal and kinematic parameters of a signal receiver and its clock by information fusion
FI113092B (fi) * 2002-05-31 2004-02-27 Ekahau Oy Paikannusepävarmuuden mittauksia ja niiden sovelluksia
CN100407852C (zh) * 2004-07-13 2008-07-30 中兴通讯股份有限公司 一种移动通讯中移动终端的定位方法

Also Published As

Publication number Publication date
JP2009529128A (ja) 2009-08-13
CN101438184B (zh) 2016-02-17
US8213956B2 (en) 2012-07-03
BRPI0621604A2 (pt) 2013-03-12
EP1991882B1 (en) 2016-11-02
BRPI0621604B1 (pt) 2019-02-19
EP1991882A1 (en) 2008-11-19
WO2007101453A1 (en) 2007-09-13
KR101106276B1 (ko) 2012-01-18
KR20080108521A (ko) 2008-12-15
CN101438184A (zh) 2009-05-20
US20090233619A1 (en) 2009-09-17

Similar Documents

Publication Publication Date Title
JP5039063B2 (ja) 携帯電子機器の状態を追跡する方法
Liu et al. Mercury: An infrastructure-free system for network localization and navigation
JP6285628B2 (ja) 移動端末及び位置推定方法並びにコンピュータ読出可能記録媒体
KR101912195B1 (ko) 목표의 실시간 위치 측정과 지도 구축 방법 및 장치
US9062978B2 (en) Tracking a body by nonlinear and non-Gaussian parametric filtering
KR101108751B1 (ko) 복수의 선택된 초기 위치 추정값을 사용하여 위치를결정하기 위한 방법 및 시스템
KR101634126B1 (ko) 위치 추정 방법, 장치 및 그 기록매체
KR101914550B1 (ko) 레이더의 표적 위치 추적 방법
US20150005000A1 (en) Method and apparatus for ascertaining a location of a personal portalbe wireless communication device
CN105424030A (zh) 基于无线指纹和mems传感器的融合导航装置和方法
Ibrahim et al. Inertial measurement unit based indoor localization for construction applications
CN111678513A (zh) 一种超宽带/惯导紧耦合的室内定位装置与系统
KR102172145B1 (ko) 추측 항법 시스템에서의 밀결합 측위 방법 및 그 장치
EP2982157B1 (en) Determining effects on communication signals
KR101045348B1 (ko) 직각 좌표계에서 도플러 속도를 이용한 표적 추적 방법
KR101462007B1 (ko) 자세 추정 장치 및 자세 추정 방법
CN115103299A (zh) 一种基于rfid的多传感器融合定位方法
KR20190113013A (ko) 비콘 신호와 핑거프린트 맵 기반의 실내 측위 방법 및 시스템
Choe et al. Point-mass filtering with boundary flow and its application to terrain referenced navigation
CN110426562B (zh) 基于分层搜索和距离空间投影的高精度闪电三维定位方法
Sung et al. TDoA based UGV localization using adaptive Kalman filter algorithm
CN104469782A (zh) 一种移动授权用户频谱检测与其地理位置估计方法与装置
Sirola et al. Local positioning with parallelepiped moving grid
Duan et al. Toward swarm robots tracking: A constrained Gaussian condensation filter method
CN114370879B (zh) 基于水下环境特征的auv鲁棒vbhiakf-slam导航方法

Legal Events

Date Code Title Description
A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20110829

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20111128

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20120309

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20120522

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

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20120706

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20150713

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

Ref document number: 5039063

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

S531 Written request for registration of change of domicile

Free format text: JAPANESE INTERMEDIATE CODE: R313531

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313113

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

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

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