(57)【要約】  加重最小二乗最適化を組み込んだ無線位置検出システムの結果を改良する方法は、最大事後（ＭＡＰ)確率メトリックを用いて、特定の測位問題（例えば、ＵＴＤＯＡ)の特性を組み込むことによって、加重最小二乗法を一般化する。加重最小二乗法は、通例、ＴＤＯＡ、ならびにＴＤＯＡ／ＡＯＡおよびＴＤＯＡ／ＧＰＳ混成システムを含む、関連のある位置検出システムによって用いられる。組み込まれる特性には、ＴＤＯＡ誤差についての経験的情報、および他のネットワーク・エレメントに対する移動体の位置の確率分布が含まれる。ＴＤＯＡ誤差分布および事前移動体位置をモデリングする技法を提供する。この新しい確率分布モデルを用いて、ＭＡＰ判断メトリックを計算する方法を提供する。フィールド・データを用いる検査は、この方法が既存の加重最小二乗法に対してかなりの改良をもたらすことを示す。【選択図】      図４ａ
2012514210.tif 000002
【特許請求の範囲】
【請求項１】
誤差検出方法であって、
フィールド・データを得るステップ（４０２）であって、前記フィールド・データが、信号相関モデルにおいて用いる基準線または位置依存値を有する、ステップと、
（１）信号相関モデルおよび関連する測定パラメータ（４０５）、（２）相関マトリクス規則（４０４）、ならびに（３）事前位置のモデル（４０６）を得るために、前記フィールド・データを分析するステップ（４０３）と、
前記測定の推定した可変性に基づいて、測定値に対する重みを計算するステップ（４０７）と、
共分散マトリクス（４０９）を発生するために、前記相関マトリクス規則と共に前記重みを用い、逆共分散マトリクス（４１０）を計算するステップと、
最大の事後（ＭＡＰ）メトリックを有する位置を発見するために、地理的領域にわたって反復サーチを行うステップ（４１５）と、
停止条件に達したと判定するステップ（４１９）と、
最大のＭＡＰメトリックを有する地理的位置を位置解として報告するステップと、
を含む、誤差検出方法。
【請求項２】
請求項１記載の方法であって、更に、スキューに対する測定パラメータと、前記測定に対する相関値との間のマッピングを規定する表を作成するステップを含む、方法。
【請求項３】
請求項１記載の方法であって、更に、スキューに対する測定パラメータと、基準線の数との間のマッピングを規定する表を作成するステップを含む、方法。
【請求項４】
請求項１記載の方法において、更に、位置検出受信機を外部アンテナにリンクする異なる受信ポート間の相関を分析し、その用途に対する相関値および規則を与えるステップを含む、方法。
【請求項５】
請求項１記載の方法において、前記反復サーチが、各反復において地理的サーチ空間分解能を低下させ、内挿補間によって新たな検査地点を発生する、分解能ループを含む、方法。
【請求項６】
請求項５記載の方法において、前記サーチが、実行前に、前の反復の最小誤差地点を中心に置き直され、現在の地理的領域サーチ空間内にある検査地点が、個々にサーチされ、検査地点毎に、ＭＡＰメトリックを計算する、方法。
【請求項７】
請求項１記載の方法において、前記反復サーチが、前記共分散マトリクス、誤差モデル、および測定パラメータ表を用いるＭＡＰメトリック計算を含む、方法。
【発明の詳細な説明】
【技術分野】
【０００１】
関連出願に対する相互引用
[0001]  本願は、２００８年１２月３０日に出願された米国特許出願第１２／３４６，１４６号の優先権を主張する。この特許出願をここで引用したことにより、その内容全体が本願にも含まれるものとする。
技術分野
[0002]  本願は、一般的には、ワイヤレス位置検出の分野、即ち、ワイヤレス・デバイスの位置を推定するシステムおよび方法に関し、更に特定すれば、一般化した誤差分布を用いる方法に関する。
【従来技術】
【０００２】
[0003]  連邦通信委員会がＰＳＡＰ−レベル位置検出精度命令に踏み出したため、異なる位置検出技術の方法を改良することが必要となった。本明細書において記載する主題は、通信および位置検出技術の分野に関する。これは、汎地球測位システム（ＧＰＳ）、アップリンク到達時間差（ＵＴＤＯＡ）、および 高度前進リンク三角測量（ＡＦＬＴ）というような位置検出技術の精度を高める手段を提供する。
【０００３】
[0004]  位置推定の一般的な手法の１つは、時間差、疑似距離、または電力レベルというような測定量から、加重最小二乗解を求めることである。加重最小二乗解は、入力誤差が独立しておりそして正規分布(Gaussian)（J. CafferyのWireless Location in CDMA Cellular Radio System”（ＣＤＭＡセルラ無線システムにおけるワイヤレス位置検出）, Boston-London: Kluwer Academic Publishers, 2000を参照のこと）であるときに最大尤度（ＭＬ）に達することが知られている。しかし、実際に遭遇するより一般的な条件の下では、これを行うことができない。例えば、ＴＤＯＡ誤差は、予測したマルチパス遅延プロファイル(multipath delay profile)の前端に対して正となる傾向がある。以下で説明するように、完全でない先端の検出、および非見通し線（ＮＬＯＳ：non-line-of-sight）伝搬というような、様々な要因がこれら正の誤差に寄与する。その結果、基準線毎(per-baseline)の誤差分布にスキューが生じる。このスキューが、基本的な加重最小二乗法の精度を低下させる。対照的に、本明細書において記載する方法は、このスキューの知識を利用して改良した結果を得る。更に、多くの場合これらの誤差間における相関は求めることができる。例えば、同じセクタにおいて別々のマルチパス成分を受信することができ、１つのサイトにＮＬＯＳ条件が存在することもあり、そして基準信号によって共通誤差が混入されることもあり得る。これらの相関は、以下で説明するように、最大事後確率（ＭＡＰ：maximum a posteriori）アルゴリズムに組み込むことができる。また、このフレームワークは、位置検出解における事前移動体位置分布の推定値を組み込むために用いることもできる。
【０００４】
[0005]  ＵＴＤＯＡは、任意のタイプの移動局（ＭＳ）から送信された任意の信号を任意の基地局において受信してＵＴＤＯＡ測定値を求めることを可能にする、ネットワーク・ベース技術である。図１に示すように、規準基地局が、各協同基地局とほぼ同時に、受信信号を測定する。
【０００５】
[0006]  図１は、測位のために移動体デバイスに利用可能な、または移動体デバイスから入手可能な信号の理想化したモデルを示し、本発明は、測位推定値（位置検出試行とも呼ばれる）の精度を高めるために用いることができる。また、この図は、ワイヤレス位置検出のためのシステム・コンポーネントも特定する。図１において、米国のＮａｖｓｔａｒ汎地球測位システム（ＧＰＳ）のようなグローバル・ナビゲーション衛星システム１０１（ＧＮＳＳ）が、鮮明な符号分割多元接続（ＣＤＭＡ）スペクトル拡散信号１０７をブロードキャストする。この信号は、緯度、経度、高度、および速度のＴＤＯＡ位置推定のために、空間に配備された移動体和やレス・デバイス１０２によって用いられる。位置計算のための衛星信号１０７を受信するための移動体デバイス１０２が装備されていない場合、ＴＤＯＡまたは到達時刻（ＴＯＡ）計算を用いるアップリンクおよびダウンリンク地上ワイヤレス技法を用いて、位置推定値を求めることができる。ダウンリンク（ネットワーク・ベース送信機からデバイスへ）ＴＤＯＡまたはＴＯＡ技法を用いる地上広域ワイヤレス位置検出技法には、高度前進リンク三角測量（ＦＬＴ）［ＩＳ−９５、ＩＳ−２０００］、改良到達時間差（Ｅ−ＯＴＤ）［ＧＳＭ］および観測到達時間差（ＯＴＤＯＡ）［ＵＭＴＳ］、ならびに分散ビーコン技法が含まれる。地上ダウンリンク技法は、移動体デバイス１０２がネットワーク・ベース送信機１０３、１０４からのダウンリンク無線信号１０８を測定し、次いで無線リンク（１つまたは複数）１０９、バックホール設備１１３、およびワイヤレス通信ネットワーク１１０を用いて、緯度、経度、そして場合によっては高度に変換するために、収集無線測定値を位置判定エンティティ（ＰＤＥ）１０６に伝えることを必要とする。
【０００６】
[0007]  アップリンク（デバイスからネットワーク・ベース受信機）ＴＤＯＡまたはＴＯＡ技法を用いる地上広域ワイヤレス位置検出技法には、Ｕ−ＴＤＯＡ、Ｕ−ＴＤＯＡ／到達角度（ＡｏＡ）混成、およびＵ−ＴＤＯＡ／補助ＧＰＳが含まれる。Ｕ−ＴＤＯＡおよび混成は、現在ＣＤＭＡ［ＩＳ−９５、ＩＳ−２０００］、ＧＳＭ、ＵＭＴＳ、ＷｉＭＡＸ（８０２．１６ｅ／ｍおよび８０２．２０）において機能しており、来るべき長期発展（ＬＴＥ：Long-Term-Evolution）ＯＦＤＭベース・ワイヤレス無線アクセス・ネットワーク（ＲＡＮ）のために概念的に機能する。地上アップリンク技法は、移動体デバイス１０２の送信１０９を、ネットワーク・ベース受信機（この場合、セル・サイト１０３、１０４内に一緒に配置されている）が測定することを必要とする。次いで、測定データはバックホール１１１によって、緯度、経度、速度、そして場合によっては高度に変換するために、位置判定エンティティ（ＰＤＥ）１０６に伝えられる。前述のワイヤレス位置検出技法とは関係なく、無線信号飛行時間(time-of-flight)の判定は、デバイス１０２の実際の位置の精度高い判定に対する鍵となる。図１では、加算的または減算的干渉による信号反射、回折、および減衰の実世界における影響は、示されていない。
【０００７】
[0008]  図１のシステムでは、基準基地局における受信信号を、協同基地局における受信信号と相関付けることによって、到達時間差が判定される。協同局は、それらのＴＤＯＡ測定値を位置判定エンティティ（ＰＤＥ）に送り、この位置判定エンティティにおいて、位置検出解が求められる。しかしながら、加法性ノイズおよび信号レベル変動によって、測定値に対する減損が生ずる可能性がある。これらの減損は、協同基地局において移動体信号の存在を検出する感度に影響を及ぼす可能性がある。推定に対する他の障害も、協同局見通し線（ＬＯＳ）経路遅延を検出する能力に影響を及ぼす。
【０００８】
[0009]  図２ａ、図２ｂ、図２ｃ、および図２ｄは、建造物のような物体がどのように直接経路を遮断して、アップリンク、ダウンリンク、ＧＮＳＳ、および混成ＧＮＳＳ／アップリンク・システム（ＧＮＳＳは、グローバル・ナビゲーション衛星システムを表す）を含む異なる位置検出環境において、非見通し線障害が生ずるかについて示す。建造物周囲を通過する回折経路が、非常に減衰されたまたは完全に遮断された直接経路よりも遅れて、受信機に到達する。加えて、障害物からの反射も散乱を発生する可能性があり、この散乱が異なる経路の到達時間の分散を生ずる。図２ａには、アップリンク・ワイヤレス位置検出システムの一例が図示されている。移動体デバイス１０２が信号１０９を送信する。場合によっては、基準受信機２０３に対するように、無線信号が直接受信されることもある（見通し線、即ち、ＬＯＳの場合）。しかし、他の受信機１０４は回折信号２０２または反射信号２０３を受信するかもしれない。各場合において、元のアップリンク信号１０９も受信されるかもしれず、あるいは障害物２０１によって遮断、減衰、または遅延された元の信号を有するかもしれない。
【０００９】
[0010]  図３は、最初の到達の検出を困難にし、ＴＤＯＡ誤差のスキューの原因となる障害を示す。図３における参照番号は、次のように用いられている。
【００１０】
３０３＝送信時刻
３０４＝検出しきい値
３０５＝見通し線（ＬＯＳ）飛行時間
３０６＝遅れ時間
３０７＝報告されるＴＯＡまたはＴＤＯＡの基準
３０８＝遅延拡散
３０９＝失われた信号成分
[0011]  図３は、振幅３０２対時間３０１のプロット３００における、マルチ・パス劣化信号の到達時刻を示す。時点３０３において信号が送信され、この信号は３０５として示す潜在的な直接経路見通し線を有する。最も早い信号成分の到達は検出されない。何故なら、これらは検出しきい値３０４未満のレベルで到達するからである。検出しきい値３０４は、過剰な偽警報を回避するように維持しなければならない。最も早い到達の検出イベントを逸したことにより、報告されるＴＯＡまたはＴＤＯＡが、望まれるＬＯＳ  ＴＯＡまたはＴＤＯＡよりも大きくなる。この例では、しきい値３０７よりも大きな最初の信号は、本当の最初に到達した信号成分から３０６の遅延を生ずる。加えて、最も早く到達するマルチパス成分が、ＮＬＯＳ遅延を生ずるＮＬＯＳ伝搬のために、予期したよりも遅く到達する可能性がある。これによっても、報告されるＴＤＯＡがＬＯＳ  ＴＯＤＡよりも大きくなる。これらの要因は、ＴＤＯＡ測定値と、測位アルゴリズムによってサーチまたは計算されるＬＯＳ  ＴＤＯＡとの間における誤差にスキューが生じる。本明細書において記載する本発明の解決策における測位判断は、これらの要因によって生ずる誤差のスキュー、および誤差分布の非正規分布形状の双方を利用する。
【００１１】
［0012］”Bayesian-update based location prediction method for CDMA Systems”（CDMAシステムのためのベイズ更新に基づく位置予測方法）と題するK. Chang et al.,の２００３年５月１３日付け米国特許第６，５６４，０６５号に記載されている方法は、シミュレーションを用いた事後電力分布から下した位置判断によって、CMDAパイロット・チャネル測定値から電力レベルを予測するように思われる。”Method of precise position determination”（正確な位置判定方法）と題するE. Freiの１９９３年１０月１２日付けの米国特許第５，２５２，９８２号に記載されている方法は、事後ＲＭＡ誤差を用いるＧＰＳ位置検出解に対する位相の曖昧さを反復的に求める加重最小二乗法を用いて、ガウス誤差を憶測するように思われる。
【発明の概要】
【発明が解決しようとする課題】
【００１２】
[0013]  加重最小二乗最適化を組み込んだ無線位置検出システムの結果を改良する方法は、最大事後（ＭＡＰ)確率メトリックを用いて、特定の測位問題（例えば、ＵＴＤＯＡ)の特性を組み込むことによって、加重最小二乗法を一般化する。先に論じたように、ＷＬＳ方法は、通例、ＴＤＯＡ、ならびにＴＤＯＡ／ＡＯＡおよびＴＤＯＡ／ＧＰＳ混成システムを含む、関連のある位置検出システムによって用いられる。組み込まれる特性には、ＴＤＯＡ誤差についての経験的情報、および他のネットワーク・エレメントに対する移動体の位置の確率分布が含まれる。ＴＤＯＡ誤差分布および事前移動体位置をモデリングする技法を提供する。この新しい確率分布モデルを用いて、ＭＡＰ判断メトリックを計算する方法を提供する。
【課題を解決するための手段】
【００１３】
[0014]  例示的な実施形態は、誤差検出方法を提供する。この方法は、フィールド・データを得るステップであって、前記フィールド・データが、前記信号相関モデルにおいて用いる基準線または位置依存値を有する、ステップと、（１）信号相関モデルおよび関連する測定パラメータ、（２）相関マトリクス規則、ならびに（３）事前位置のモデルを得るために、前記フィールド・データを分析するステップと、前記測定の推定した可変性に基づいて、測定値に対する重みを計算するステップと、共分散マトリクスを発生するために、前記相関マトリクス規則と共に前記重みを用い、逆共分散マトリクスを計算するステップと、最大の事後（ＭＡＰ）メトリックを有する位置を発見するために、地理的領域にわたって反復サーチを行うステップと、停止条件に達したと判定するステップと、最大のＭＡＰメトリックを有する地理的位置を位置解として報告するステップとを備えている。
【００１４】
[0015]  本明細書において記載する方法は、以下の様々な主要な改革を含むが、必ずしもこれらに限定される訳ではない。
【００１５】
[0016]  分析事前分布(Analytical a priori distribution)：実際の位置を与える経験的データを用いて、基準タワーからの正規化した距離の位置検出解に対する分布を求めて、サーチエリアにおけるタワーに対する事前位置の全体的な形状をモデル化する。事前位置分布の形状を近似する指数分布が示され、その分散が、経験的データから計算される。
【００１６】
[0017]  分析ＴＤＯＡ誤差分布：二重指数分布モデルを一般化して、指数にスキューおよび任意の累乗(power)を組み込む。経験的データから、モデル・パラメータを推定する。
【００１７】
[0018]  マルチパス／ＮＬＯＳ誤差指標：ＴＤＯＡ誤差分布の主要な指標には、基準線の数、予測されたマルチパス補正（観察された信号パラメータおよび／またはローカルＲＦ環境の知識に基づいて）、および各基準線のＴＤＯＡ相関が含まれる。経験的データを分析し、条件付き誤差分布を発生することによって、これらの指標からモデル・パラメータを求める方法を提供する。基準線毎に、スキューのようなモデル／パラメータをこれらの指標から計算する。
【００１８】
[0019]  ＴＤＯＡ誤差相関：前述の分析ＴＤＯＡ誤差分布を有する基準線間において相関付けした誤差に対する事後確率を計算する方法を提供する。これらの相関は、対応する合同誤差確率分布(joint error probability distribution)によって、ＭＡＰアルゴリズムに組み込まれる。
【００１９】
[0020]  共通バイアス軽減方法：更に一般的な分布を用いると、測定値に存在する可能性がある共通バイアスに対する分析解を求めるのが難しくなる。バイアスの除去方法を、種々の複雑さ−性能トレードオフと共に提供する。
【００２０】
[0021]  反復調節：以上の方法を応用する反復手順を制作する。この手順は、初期化動作と、残余値の推定を含む。  
【００２１】
[0022]  本発明の技術の他の特徴について、以下で説明する。
【図面の簡単な説明】
【００２２】
[0023]  添付図面は以下を含む。
【図１】図１は、測位ネットワークの図である。
【図２ａ】図２ａは、ＬＯＳ経路遅延推定に対する障害を示す。
【図２ｂ】図２ｂは、ＬＯＳ経路遅延推定に対する障害を示す。
【図２ｃ】図２ｃは、ＬＯＳ経路遅延推定に対する障害を示す。
【図２ｄ】図２ｄは、ＬＯＳ経路遅延推定に対する障害を示す。
【図３】図３は、測定誤差スキューの原因の図。
【図４ａ】図４ａは、ＭＡＰ誤差分布方法の成分である。
【図４ｂ】図４ｂは、ＭＡＰ誤差分布方法の成分である。
【図５】図５は、誤差分布モデリング・プロセス。
【図６】図６は、事前分布データ分析の論理フロー。
【図７】図７は、事前位置分散と経験的モデルとの比較である。
【図８】図8は、条件付き誤差分布データ分析である。
【図９】図９は、誤差分布対加重最小二乗法の総体的サンプル(sample overall)を示す。
【図１０】図１０は、サンプル全体誤差分布対新しい粗いモデルを示す。
【図１１】図１１は、スキュー率の相関に対するサンプル依存性を示す。
【図１２】図１２は、サンプル全体条件付き誤差分布対新たなモデル（小さなスキュー）の関係を示す。
【図１３】図１３は、サンプル条件付き誤差分布対新たなモデル（大きなスキュー）の関係を示す。
【図１４】図１４は、ＭＡＰ判断メトリック計算の論理フローである。
【図１５】図１５は、二重指数仮定（ｐ＝１、ξ＝１）に対するバイアス計算である。
【図１６】図１６は、条件付き誤差計算である。
【発明を実施するための形態】
【００２３】
[0040]  図４ａおよび図４ｂは、ＭＡＰ誤差検出方法の例示的な実施態様のコンポーネントを示す。図示のように、ＭＡＰプロセスはステップ４０１において開始する。フィールド・データ４０２を分析して（４０３）、１組の信号相関規則およびモデルを求める。これらのモデルおよび関連付けられた測定パラメータを、モデルにおいて用いられる基準線または位置(location)（または、場所(position)、ここで「位置」および「場所」という用語は本明細書では相互交換可能に用いられるものとする）依存値を有することができるフィールド・データから展開する(develop)。例えば、誤差のスキューは、低相関ＵＴＤＯＡ測定値の方が大きいこともある。このため、表４０５を作成することができる。この表は、スキューに対するモデル・パラメータと測定に対する相関値との間のマッピングを規定する。同様に、事前位置についてモデルおよび表を計算する（４０６）。また、フィールド・データ分析プロセスが、位置検出受信機(location receiver)（例えば、位置測定ユニット（ＬＭＵ）または信号収集システム（ＳＣＳ））を外部アンテナにリンクする異なる受信機ポート間の相関を分析し、それらの用途に合わせて相関値および規則を与える。例えば、ＮＬＯＳ効果のために、同じ位置（共通サイト・ポート）にあるポートにおいて、誤差の小さな相関がある場合がある。一旦フィールド・データを分析したなら、推定した測定の可変性に基づいて、測定値に対して重みを計算する。次いで、このこれらの重みを相関マトリクス規則と共に用いて、ポート毎の共分散を発生し、４０９において反転させる。
【００２４】
[0041]  図４ｂに示すように、次に、地理的領域にわたって反復サーチを実行する。その目的は、ＭＡＰメトリックが最も大きい位置を発見することである。一旦反復サーチを開始したなら（４１１）、分解能ループ(resolution loop)に入る（４１２）。このループにおいて、反復毎に地理的サーチ空間分解能を低下させて、内挿補間によって新たな検査地点を発生する。このサーチは、実行前に、前の反復の最小誤差地点を中心に置き直すことができる。現在の地理的領域における検査地点毎に、ＭＡＰメトリックを計算し（４１５）、メトリックが最も小さい地点を選択することによって、現在の地理的領域をサーチする（４１４）。ＭＡＰメトリックの計算には、共分散マトリクス４０４、誤差モデル、および測定パラメータ表４０５を用いる。サーチ空間内に未だ他にも検査地点が存在する場合、処理ロジック４１７は４１４に戻ってサーチを再開する。現在の分解能によって与えられる現在のサーチ空間の中に検査していない検査地点がない場合、ＭＡＰプロセスは、予め設定されている分解能限度集合(resolution limit set)をチェックする（４１８）。未だ最も高い分解能に達していない場合、プロセスはステップ４１２に戻る。それ以外の場合、ＭＡＰプロセスは、１つ以上の停止規準に到達したか否かチェックする（４１９）。停止規準を満たす場合、ＭＡＰプロセスは終了し（４２０）、メトリックが最も大きい地理的位置が位置検出解を与える。
【００２５】
[0042]  本発明の解決策の目的は、誤差の事後確率をモデル化し、この確率を最大にする位置検出解を求めることである。ベイズの定理から（A. Papoulis, Probability Random Variables, and Stochastic Processes（確率ランダム変数、および確率的プロセス）、 MacGraw Hill Inc., New York, NY, 1984を参照のこと）、ランダム位置ベクトルＬの条件付き確率密度関数は、Ｎ個の測定誤差ｅに関して、以下のように与えられる。
【００２６】
【数１】
2012514210.tif 000003
【００２７】
ここで、
【００２８】
【数２】
2012514210.tif 000004
【００２９】
は、ランダム位置ベクトル、
【００３０】
【数３】
2012514210.tif 000005
【００３１】
は、ＴＯＤＡ誤差＋共通バイアス、
そして、
τｉ（ｘ、ｙ、ｚ）は、ｉ番目の基準線に対する、地点ｘ、ｙ、ｚにおけるＬＯＳ  ＴＤＯＡである。
【００３２】
【数４】
2012514210.tif 000006
【００３３】
は、ｉ番目のＴＤＯＡ基準線測定値、
Ｂは、測定値の中に存在し得る共通バイアスである。
【００３４】
[0043]  計算を簡略化するために、（１）の対数を最大にすることができる。何故なら、（１）を最大にする位置は、（１）の対数を最大にする位置でもあるからである。（１）の自然対数は、
【００３５】
【数５】
2012514210.tif 000007
【００３６】
となる。
【００３７】
[0044]  最後の項は位置には依存しないので、異なる位置を考慮するときには、これは一定となるため、無視することができる。この結果、全ての位置において最大化すべき以下の関数が残る。
【００３８】
【数６】
2012514210.tif 000008
【００３９】
[0045]  最初の項は、事後誤差確率密度の対数であり、２番目の項は事前確率密度の対数である。
誤差分布モデリング・プロセス
[0046]  誤差分布モデリング・プロセスを図５に示す。最初に、フィールド・データ４０２を分析して、種々の基準線および位置特定測定パラメータの事前位置分布および誤差分布に対する影響を判定する。次に、異なる基準線に対する誤差間の相関を分析して、相関値、およびこれらの相関を得るための規則を決定する。誤差分布全体に当てはまる(fit)誤差分布の粗いモデルを求める（５０４）。フィールド・データの事前分布に基づいて、事前分布のモデルを求める。次に、種々の位置特定および基準線特定測定パラメータに基づいて、修正可能にするために、誤差分布モデルをリファインする。最後に、相関値、およびこれらの値を適用するための関連規則に基づいて、相関マトリクスを発生する。
事前分布
[0047]  適した事前分布を求めるロジックを図６に示す。位置サーチエリアは、位置特定データによっては、非常に異なる寸法を有する可能性があり、位置検出の問題も非常に異なるネットワーク構成を有することがあるので、サーチエリア内にある実際の位置の関数であるモデル・パラメータの分布を求めることが望ましい。サーチエリア内にある位置は、事前位置である。式（５）において用いるための事前位置確率を決定するために計算し用いることができる事前位置の関数であるモデル・パラメータを求めることが望ましい。最初に、フィールド・データを分析して、候補パラメータを決定する。例えば、基準タワーからの距離についての統計が、事前位置を表すことがある。あるいは、サービング・タワーからの距離を考慮してもよい。これらの距離を最大サーチ範囲に正規化するとよい。パラメータ毎に、範囲およびビン・サイズを選択しなければならず、フィールド・データにおける実際の（真の）位置に基づいて更新された事前分布を表すヒストグラムも選択しなければならない。パラメータ計算と共に、各潜在的分布を事前パラメータ・データベースに格納する。
【００４０】
[0048]  一旦潜在的事前分布を種々のモデル・パラメータについて計算したなら、図５に示したようにモデルを選択する。このモデルは、フィールド・エリアからの分布を表し、そして計算上効率的であるように、モデルを選択する。基準タワーからサーチ領域のエッジまでの最大距離に正規化された、基準タワーから候補位置の距離となるように、モデル・パラメータ例を選択する。これによって、以下のような、三次元ランダム・ベクトルＬから、基準タワーからの正規化ランダム距離Ｄへの変換が得られる。
【００４１】
【数７】
2012514210.tif 000009
【００４２】
ここで、
xref、yref、zrefは、基準タワーの位置座標である。
【００４３】
Rmaxは、基準タワーからサーチ領域のエッジまでの最大距離である。
【００４４】
以下のように、指数となるようにモデル例を選択する。
【００４５】
【数８】
2012514210.tif 000010
【００４６】
ここで、λａ＝１１が、フィールド・データに当てはまるように選択された。図７では、このモデルを、正確な一致を示すフィールド・データと比較する。
誤差分布
[0049]  また、誤差分布のモデルを得るためにも、フィールド・データを分析する。図８は、事前データ分析におけると同様の一連のステップを示す。この図では、測定パラメータを分析して、誤差分布においてどのパラメータが大きな変化の原因となるのか判定する。測定パラメータの例には、次のものが含まれる。
【００４７】
基準線毎のＵＴＤＯＡ相関、
基準線毎のマルチパス相関係数、
位置毎の測定数。
【００４８】
[0050]  条件付き統計および全体的統計を蓄積するために、これらのパラメータの範囲およびビン・サイズを決定する。
【００４９】
[0051]  誤差分布全体のサンプルが、図９に示されており、加重最小二乗法によって想定される正規分布(Gaussian distribution)と比較される。この図から、ガウス仮定は分布全体の形状を適正にモデル化していないことは明白である。分布全体の右へのスキューは明白である。
【００５０】
[0052]  この分布全体は、図５に示したような粗い誤差モデルへの入力を与える。ｉ番目の限界誤差(marginal error)の誤差モデル例を、次のようになるように決定する。
【００５１】
【数９】
2012514210.tif 000011
【００５２】
ここで、
piは、０よりも大きい任意の指数型ベキ(exponential power)であるモデル・パラメータである。
【００５３】
riは、分布のスキューを示す正の比率であるモデル・パラメータである。
【００５４】
σiは、ｉ番目の基準線に対する標準偏差である。
【００５５】
[0053]  ｋおよびＡの値は、所与のriおよびpiについて以下の条件、
【００５６】
【数１０】
2012514210.tif 000012
【００５７】
を満足するように選択される。正規分布では、ri＝１、pi＝２、そしてｋ＝１／２となる。二重指数分布またはラプラス分布では、ri＝１、pi＝１、そしてｋ＝√２となる。
【００５８】
[0054]  粗いモデリング・ステップでは、フィールド・データと一致する式（８）におけるモデル・パラメータの値を計算する。図１０は、pi＝１．１およびri＝１．１として、粗いモデルを適用する効果を示す。この図は、図９に示したガウス仮定よりも、遙かに正確にフィールド・データと一致することを示す。
【００５９】
[0055]  図５における「細かい誤差モデル決定」ブロック５０６への入力として、条件付き誤差分布が用いられる。条件付き誤差に対して、測定パラメータ・ビン毎に、誤差の平均および標準偏差に基づいて、スキュー率を計算する。ビンの中心は、誤差分布の条件付き値となる。
【００６０】
[0056]  以下のようにして、条件付き分布の平均および標準偏差に関するスキューを求めることができる。条件付き誤差分布を二重指数として近似する場合、指数におけるスケーリング係数は、
【００６１】
【数１１】
2012514210.tif 000013
【００６２】
となる。ここで、σは条件付き分布の標準偏差である。
【００６３】
[0057]  riを推定するために、２つの別のスケーリングした指数分布を考慮し、これらの１つを０を中心として折り返す(flip)。分布の双方の成分をスケーリングして、１／２まで積分する。その結果、条件付き分布の平均ｍを、スケーリング係数について、次のように置くことができる。
【００６４】
【数１２】
2012514210.tif 000014
【００６５】
ここで、λLおよびλRは、それぞれ、０の左側および右側における別の指数分布成分の指数におけるスケーリング係数である。尚、スキューの全てはλに対するλRの変化によるものであると仮定すると、λL≒λという仮定が可能になる。λRについて式（１０）を解き、λL≒λとして（９）を用いると、次の式が得られる。
【００６６】
【数１３】
2012514210.tif 000015
【００６７】
[0058]  次に、（９）および（１１）からスキュー率は次のようになる。
【００６８】
【数１４】
2012514210.tif 000016
【００６９】
[0059]  次に、条件付き誤差分布からのσおよびｍの値を用いて、式（１２）を用いてriを計算することができる。
【００７０】
[0060]  ＵＴＤＯＡ相関の関数としてのスキュー率の一例を図１１に示す。条件付き分散毎に、式（１２）において平均ｍおよび標準偏差σを用いて、riを計算する。図では、相関が低い場合には、スキューが増大することが明白である。サンプル数およびビン・サイズの関数である可変性を平均化するために表の参照の代わりに用いることができる線型データ当てはめを示す。「細かい誤差モデル決定」ブロック５０６（図５）の結果、ＵＴＤＯＡ相関のような測定パラメータ値の、スキューのようなモデル・パラメータへのマッピングが得られる。（８）におけるモデル・パラメータpiに対する同様の調節も、種々の測定パラメータ値の関数として行うことができる。
【００７１】
[0061]  細かいモデル調節の結果例を図１２に示す。ここでは、相関が高い方に小さなスキューがある。図１３は、相関が低い方に大きなスキューがある例を示す。これらの図は、誤差分布モデルにおいて、図１０に対して更なる改良が達成されることを示す。
相関マトリクス  
[0062]  また、図５に示したような異なるポートについて誤差間に存在する相関を分析するためにも、フィールド・データが用いられる。誤差間の相関は、種々の規則または条件に基づいて計算される。２つのポート間における相関を計算する規則例には、次のことが含まれる。
【００７２】
同じセクタにあるポート間に固定相関を適用する。
【００７３】
同じサイトにあるポート間に固定相関を適用する。
【００７４】
協同ポートと基準ポートとの間に固定相関を適用する。
【００７５】
[0063]  規則毎に、フィールド・データ統計から、誤差に対する正規化相関値または相関係数を計算する（A. Papoulis, Probability Random Variables, and Stochastic Processes（確率ランダム変数、および確率的プロセス）、 MacGraw Hill Inc., New York, NY, 1984を参照のこと）。相関値および規則は、「共分散マトリクス入力および反転実行」ブロック４０９（図４ａ）に入力を与える。相関マトリクス発生江ブロックでは、相関値のポート毎のマトリクスを計算する。複数の規則が１対のポートに適用される場合、このマトリクスでは最も大きな相関値が用いられる。
加重および分散の計算
[0064]  各基準線に対する加重は、クラメール・ラオ境界（R. MacDonough, A. Whalen, Detection of Signals in Noise, 2nd Ed., Academic Press, San Diego, CA, 1995を参照のこと）からのＲＭＳ誤差に基づく。ＡＷＧＮ(加法的ホワイト・ガウス・ノイズ）におけるＴＤＯＡ  ＲＭＳ誤差の下限は、次のようになる。
【００７６】
【数１５】
2012514210.tif 000017
【００７７】
ここで、Ｗは信号帯域幅、Ｔはコヒーレント積分長、そしてρiはｉ番目の基準線の相関である。ＡＷＧＮにおける平均誤差は０に近いので、誤差の標準偏差はほぼＲＭＳ誤差となる。重みは、ＲＭＳ誤差の二乗に対して１であり、以下のような理論的加重が与えられる。
【００７８】
【数１６】
2012514210.tif 000018
【００７９】
[0065]  これらの加重動作例は、図４に示したような「重み計算」ブロック４０７におけるフィールド・データ分析後に行われる。ＲＭＳ誤差を更に増大させるマルチパスのような劣化を考慮する他の効果を含めてもよい。
共分散マトリクスの計算
[0066]  合同誤差密度を用いて判定を下すために、共分散マトリクスが必要になることがある。共分散マトリクスＣは、以下のように計算される、ｉ番目のポートとｊ番目のポートとの間における共分散のポート毎のマトリクスである。
【００８０】
【数１７】
2012514210.tif 000019
【００８１】
ここで、βijは、相関マトリクスからのｉ番目のポートとｊ番目のポートとの間における相関係数である。
【００８２】
[0067]  あるいは、ポート間における相関レベルが小さすぎると見なされる場合には、計算効率のために、このステップを迂回してもよい。判断基準の一例は、βijの少なくとも１つが相関しきい値を超過する場合に、共分散マトリクスを用いることである。このしきい値を超過しない場合、独立誤差分析を用いるために、フラグをセットする。
ＭＡＰ判断メトリック計算
[0068]  合同誤差密度を用いたＭＡＰ判断計算は、相関付けられたＵＴＤＯＡ誤差に対して更なる一般化を用いる。合同ガウス誤差から開始して、事後確率は次のようになる。
【００８３】
【数１８】
2012514210.tif 000020
【００８４】
ここで、Ｇは一定である。個々のＵＴＤＯＡ誤差に関して、
【００８５】
【数１９】
2012514210.tif 000021
【００８６】
ここで、dijはＣ-１のエレメントである。限界誤差確率密度についての式（８）におけるモデルを想定すると、合同密度に対して以下の一般化が式（１７）に行われる。
【００８７】
【数２０】
2012514210.tif 000022
【００８８】
ここで、
【００８９】
【数２１】
2012514210.tif 000023
【００９０】
【数２２】
2012514210.tif 000024
【００９１】
そして、
【００９２】
【数２３】
2012514210.tif 000025
【００９３】
（１８）および（７）の自然対数を（５）に代入すると、次の式が得られる。
【００９４】
【数２４】
2012514210.tif 000026
【００９５】
[0069]  目的は（２２）を最大にするｘ、ｙ、ｚおよびＢを求めることであるので、ｘ、ｙ、ｚおよびＢに依存しない項を無視することができ、次の式が得られる。
【００９６】
【数２５】
2012514210.tif 000027
【００９７】
[0070]  計算効率のために、式（２３）を−ｋで除算してもよく、
【００９８】
【数２６】
2012514210.tif 000028
【００９９】
を最大にするｘ、ｙ、ｚは、以下のように求められる。
【０１００】
【数２７】
2012514210.tif 000029
【０１０１】
[0071]  次いで、最小にする判断メトリックは、
【０１０２】
【数２８】
2012514210.tif 000030
【０１０３】
となる。
【０１０４】
[0072]  基準線間に低い（〜０）相互相関がある位置では、共分散マトリクスは対角線状となる。ＵＴＤＯＡ誤差間に独立を仮定すると、（２５）を次のように簡略化することができる。
【０１０５】
【数２９】
2012514210.tif 000031
【０１０６】
[0073]  基準線毎に予め計算されている重みに関して、メトリックは、
【０１０７】
【数３０】
2012514210.tif 000032
【０１０８】
となる。
【０１０９】
[0074]  図１４は、各ｘ、ｙ、ｚ値におけるＭＡＰ判断メトリック計算の論理フローを示す。考慮するｘ、ｙ、ｚについて、誤差サンプルを計算する（１４０１）。一般に、Ｍを最小にする式（２７）における共通バイアス「Ｂ」についての分析解は、計算するのが難しい。したがって、計算するのが比較的容易な２つのサンプル点におけるバイアス（pi＝１、ri＝１、およびpi＝２、ri＝１）を発見し、組み合わせて近似値を得る。誤差に対してガウス統計を仮定して、最小バイアスを計算し（１４０２）、次いで、二重指数統計を仮定した計算を行う（１４０３）。次に、２つのバイアス点を組み合わせることによって、用いるべきバイアスを求める（１４０４）。この組み合わせえは、単にガウス指数バイアスおよび二重指数バイアスの平均を取ることによって行うことができる。あるいは、任意のpiおよびriに対するバイアスを、複雑さを犠牲にして、可能な全てのバイアスにわたるサーチによって発見することもできる。これを行う場合、２つのバイアス・サンプルの平均を、サーチの開始点として用いる。あるいは、オフラインでサーチを行い、２つのバイアス・サンプルに対するバイアス結果と比較することもできる。この場合、オフライン分析を犠牲にして、２つのバイアス点からの平均百分率偏差を、組み合わせにおいて用いることもできる。サンプル・バイアスを組み合わせた後、この組み合わせバイアス１４０５、共分散マトリクス１４０７、ならびにフィールド・データから既に展開されている誤差モデルおよびパラメータ表を用いて、条件付き誤差寄与を判定する（１４０５）。次に、計算した条件付き誤差寄与１４０５、および事前位置について既に展開されているモデル１４０９を用いて、ＭＡＰメトリックに対する事前寄与を求める（１４０６）。種々の地理的マップ分解能および反復によるｘ、ｙ、ｚに対する図１４におけるメトリック計算の最小値が、図４に示したように、最終解を与える。
ガウス・バイアス  
[0075]  （２７）においてpi＝２およびri＝１と設定することによって、ガウス・バイアスを求める。尚、ri＝１について、
【０１１０】
【数３１】
2012514210.tif 000033
【０１１１】
が、
【０１１２】
【数３２】
2012514210.tif 000034
【０１１３】
を与えることに注意すること。
ここで、Nは基準線の数であり、
【０１１４】
【数３３】
2012514210.tif 000035
【０１１５】
はバイアス無し誤差(unbiased error)である。
【０１１６】
[0076]  Ｂに対する（２８）の導関数を０に設定し、Ｂについて解くことによって、バイアスに対する最小解を求めると、次の式が得られる。
【０１１７】
【数３４】
2012514210.tif 000036
【０１１８】
式（２９）は、誤差分布が正規分布であるときのバイアスを与える。
指数バイアス
[0077]  指数バイアスは、（２７）においてpi＝１およびri＝１に設定することによって求められ、次の式が得られる。
【０１１９】
【数３５】
2012514210.tif 000037
【０１２０】
[0078]  再度、Ｂに対する（２８）の導関数を０に等しく設定し、Ｂについて解くことによって、バイアスに対する最小解を求める。Ｂに対する各項の導関数は、次の通りである。
【０１２１】
【数３６】
2012514210.tif 000038
【０１２２】
ここでＵ（ｘ）は、ユニット・ステップ関数である（A. Oppenhem and A. Willsky, Signals and Systems（信号およびシステム）, Prentice-Hall, Inc., Englewood Cliffs, NJ, 1983を参照のこと）。和の導関数を０に設定すると、次の式が得られる。
【０１２３】
【数３７】
2012514210.tif 000039
【０１２４】
[0079]  Ｂの関数としての（３２）における各項は、値Δτiに達するまでは−√Wiとなり、次いでＢｇａΔτiよりも大きくなると、１ステップで√Wiになる。これらの不連続のために、Ｂには正確な解がない。しかしながら、近似解を与えるＢの値を発見することができる。
【０１２５】
[0080]  図１５における解は、増大するΔτiにしたがって（３２）における加法を命令する(order)。次いで、（３２）をできるだけ０に近づけるＢの値を求める。これが生ずるのは、ｋ番目のステップ遷移が発生して、負項の合計が正項の合計にほぼ等しくなるＢの値のときである（即ち、
【０１２６】
【数３８】
2012514210.tif 000040
【０１２７】
ここで、Ｎは命令された加算における基準線の総数である）。図では、重みおよびサンプル・アレイが入力され、並び替えられる。停止条件を与えるために、
【０１２８】
【数３９】
2012514210.tif 000041
【０１２９】
となるしきい値を計算する。最も小さい遷移点を有する項から最も大きな遷移点を有する項まで順番に、項を蓄積する。しきい値に達した時点で、奇数個の項がある場合、値Δτｋが返され、それ以外の場合、ｋ番目の項の遷移点と、前の項の遷移点との平均を取る。
メトリック計算
[0081]  （２７）における最初の項を、図１６におけるステップにしたがって計算する。基準線毎に、ＵＴＤＯＡ相関、基準線の数、およびマルチパス・パラメータというような測定パラメータを決定する。これらの測定パラメータは、参照表からまたはパラメータ当てはめモデルを用いた直接計算によって誤差モデル・パラメータpiおよびriを決定するために用いられる。スキューに応じて、誤差を調節する。最後に、誤差が有意な相関を有するときには（２５）を用いて、また誤差が有意な相関を有していないときには（２７）を用いて、メトリックにおける合計を計算する。相関の有意性は、図４における共分散マトリクス計算の一部として判定される。最後に、事前確率を考慮に入れるために、（６）を用いて、（２７）および（２５）における最後の項を計算する。
サンプル結果
[0082]  以下の表は、加重最小二乗アルゴリズムに対するサンプルの改良を示す。約４６，０００個の位置測定値を用いて、加重最小二乗アルゴリズムおよび上述のアルゴリズムを用いて、測位誤差の分布を編集した。上述のモデルのパラメータは、３２，０００箇所の別の訓練データ集合を用いて選択した。この表は、９５番目および６７番目の百分位数において、それぞれ、約２０メートルおよび２メートルの改良を示す。平均誤差は、約１５メートルだけ改善した。
【０１３０】
【表１】
2012514210.tif 000042
【０１３１】
結論
[0083]  本発明、および以下の特許請求の範囲の保護範囲は、以上に記載した詳細には全く限定されない。ワイヤレス位置検出の当業者であれば、本明細書において開示した発明概念から逸脱することなく、例示的な実施形態に対して種々の変更を行うことができよう。
【図１】
2012514210.tif 000043
【図２ａ】
2012514210.tif 000044
【図２ｂ】
2012514210.tif 000045
【図２ｃ】
2012514210.tif 000046
【図２ｄ】
2012514210.tif 000047
【図３】
2012514210.tif 000048
【図４ａ】
2012514210.tif 000049
【図４ｂ】
2012514210.tif 000050
【図５】
2012514210.tif 000051
【図６】
2012514210.tif 000052
【図７】
2012514210.tif 000053
【図８】
2012514210.tif 000054
【図９】
2012514210.tif 000055
【図１０】
2012514210.tif 000056
【図１１】
2012514210.tif 000057
【図１２】
2012514210.tif 000058
【図１３】
2012514210.tif 000059
【図１４】
2012514210.tif 000060
【図１５】
2012514210.tif 000061
【図１６】
2012514210.tif 000062
【国際調査報告】
2012514210.tif 000063