特表 2019-521311 低地球軌道衛星通信システムにおける位置推定

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】公表特許公報(A)

(11)【公表番号】特表2019-521311(P2019-521311A)

(43)【公表日】2019年7月25日

(54)【発明の名称】低地球軌道衛星通信システムにおける位置推定

(51)【国際特許分類】

   G01S 5/02 20100101AFI20190704BHJP

   G01S 19/48 20100101ALI20190704BHJP

【ＦＩ】

   G01S5/02 Z

   G01S19/48

【審査請求】未請求

【予備審査請求】未請求

【全頁数】36

(21)【出願番号】特願2018-554392(P2018-554392)

(86)(22)【出願日】2017年5月16日

(85)【翻訳文提出日】2018年10月15日

(86)【国際出願番号】AU2017000108

(87)【国際公開番号】WO2017197433

(87)【国際公開日】20171123

(31)【優先権主張番号】2016901913

(32)【優先日】2016年5月20日

(33)【優先権主張国】AU

(81)【指定国】 AP(BW,GH,GM,KE,LR,LS,MW,MZ,NA,RW,SD,SL,ST,SZ,TZ,UG,ZM,ZW),EA(AM,AZ,BY,KG,KZ,RU,TJ,TM),EP(AL,AT,BE,BG,CH,CY,CZ,DE,DK,EE,ES,FI,FR,GB,GR,HR,HU,IE,IS,IT,LT,LU,LV,MC,MK,MT,NL,NO,PL,PT,RO,RS,SE,SI,SK,SM,TR),OA(BF,BJ,CF,CG,CI,CM,GA,GN,GQ,GW,KM,ML,MR,NE,SN,TD,TG),AE,AG,AL,AM,AO,AT,AU,AZ,BA,BB,BG,BH,BN,BR,BW,BY,BZ,CA,CH,CL,CN,CO,CR,CU,CZ,DE,DJ,DK,DM,DO,DZ,EC,EE,EG,ES,FI,GB,GD,GE,GH,GM,GT,HN,HR,HU,ID,IL,IN,IR,IS,JP,KE,KG,KH,KN,KP,KR,KW,KZ,LA,LC,LK,LR,LS,LU,LY,MA,MD,ME,MG,MK,MN,MW,MX,MY,MZ,NA,NG,NI,NO,NZ,OM,PA,PE,PG,PH,PL,PT,QA,RO,RS,RU,RW,SA,SC,SD,SE,SG,SK,SL,SM,ST,SV,SY,TH,TJ,TM,TN,TR,TT,TZ

【公序良俗違反の表示】

（特許庁注：以下のものは登録商標）

１．ＺＩＧＢＥＥ

(71)【出願人】

【識別番号】518233372

【氏名又は名称】マイリオタ ピーティーワイ エルティーディー

【氏名又は名称原語表記】Ｍｙｒｉｏｔａ Ｐｔｙ Ｌｔｄ

(74)【代理人】

【識別番号】110000729

【氏名又は名称】特許業務法人 ユニアス国際特許事務所

(72)【発明者】

【氏名】ヘイリー、デイビッド、ビクター、ロウリー

(72)【発明者】

【氏名】グラント、アレキサンダー、ジェイムズ

(72)【発明者】

【氏名】マッキリアム、ロバート、ジョージ

(72)【発明者】

【氏名】ポロック、アンドレ

【テーマコード（参考）】

5J062

【Ｆターム（参考）】

5J062AA08

5J062CC07

5J062CC13

5J062CC15

5J062CC17

5J062DD02

5J062DD05

5J062FF01

(57)【要約】

【課題】地上局においてまたは衛星内で衛星通信における端末の位置を推定するための方法について説明する。
【解決手段】衛星位置の推定値とともに、端末から衛星に送信された信号の時間遅延、ドップラー、および／またはドップラー率の推定値を使用する。端末が安定した時間基準または周波数基準を有するならば、同期端末と非同期端末の両方とともに、ならびに、オン・ボード全世界位置基準とオン・ボード全世界タイミング基準とをもつ低地球軌道衛星と、これらをもたない低地球軌道衛星の両方とともに、使用可能である。加えて、衛星または端末のどちらかによって受信されたビーコン信号の時間遅延、ドップラー、またはドップラー率の推定値も使用することができる。これらのビーコン信号はＧＰＳ Ｌ１信号を含んでよい。標準的なＧＰＳ受信機とは対照的に、端末がＧＰＳエフェメリス・データを記憶することまたは連続的に動作することを必要としない。
【選択図】図２

【特許請求の範囲】

【請求項1】

衛星通信システム内の地上ベースの端末の位置の推定のための方法であって、
端末から離れたロケーションにおいて、該端末から衛星への送信の時間遅延、ドップラー周波数、またはドップラー周波数変化率のうちの１つまたは複数を、該送信に関する情報から推定するステップであって、該送信に関する該情報が、少なくとも該衛星によって受信された信号、または該衛星によって受信された信号の１つもしくは複数のデジタル・サンプルを含む、推定するステップと、
該送信の時間における該衛星の位置の推定値を受信するステップと、
該衛星の位置の該推定値と、該時間遅延の該推定値、該ドップラー周波数の該推定値、または該ドップラー周波数変化率の該推定値のうちの少なくとも１つとを使用して、該端末位置を推定するステップと
を含む方法。

【請求項2】

前記送信に関する前記情報が、衛星から地上局への１つまたは複数の送信を介して受信され、時間遅延、ドップラー周波数、またはドップラー周波数変化率のうちの１つまたは複数を推定するステップが、該地上局に動作可能に接続された通信受信機によって実行され、前記時間遅延、前記ドップラー周波数、または前記ドップラー周波数変化率が、前記送信におけるユーザ・データの推定値を回復するために該通信受信機によって復号プロセス中に推定される、請求項１記載の方法。

【請求項3】

前記送信に関する情報が、オン・ボード時間基準またはオン・ボード周波数基準に対する、前記送信が前記衛星によって受信されたタイム・スタンプをさらに含む、請求項１または２記載の方法。

【請求項4】

前記送信に関する情報が、前記端末における時間基準または周波数基準に対する、前記送信が前記端末によって送信されたタイム・スタンプをさらに含む、請求項１から３のいずれかに記載の方法。

【請求項5】

前記送信に関する情報が、前記端末時間基準または前記端末周波数基準に対して決定された所定の送信スケジュールに基づいた、前記端末による推定送信時間をさらに含む、請求項１から４のいずれかに記載の方法。

【請求項6】

前記端末が前記衛星に送信のシーケンスを送信し、前記衛星が、該送信のシーケンスに関する情報を地上局に送信し、遅延、ドップラー、またはドップラー率のうちの１つまたは複数を前記受信情報から推定するステップが、該送信のシーケンス内の各送信に対して実行され、前記端末位置を推定するステップが、前記各送信の時間における前記衛星の位置の前記推定値と、該送信のシーケンス内の各送信に対する前記推定値の各々とを使用する、請求項１から５のいずれかに記載の方法。

【請求項7】

前記端末位置を推定するステップが、前記送信のシーケンス内の送信間の推定値間の差に対して異なって実行される、請求項６記載の方法。

【請求項8】

前記端末が、安定したクロックを備え、いかなる全世界時間基準または全世界周波数基準にも同期されず、前記端末位置を推定する前記ステップが、全世界時間基準または全世界周波数基準に対する端末オフセットを推定するステップをさらに含む、請求項１から７のいずれかに記載の方法。

【請求項9】

前記衛星が、オン・ボード位置基準を含まず、前記端末位置の推定が、前記衛星位置を推定するステップをさらに含む、請求項１から８のいずれかに記載の方法。

【請求項10】

前記衛星において、地上ベースのソースまたは宇宙ベースのソースから１つまたは複数のビーコン信号を受信するステップであって、該ビーコン信号が、ビーコン位置と、基準に対する全世界タイミングとを含み、前記衛星が、該信号のサンプルを前記地上局に転送する、受信するステップをさらに含み、前記端末のロケーションを推定するステップが、各ビーコン信号に対するビーコン時間遅延、ビーコン・ドップラー周波数、またはビーコン・ドップラー周波数変化率のうちの少なくとも１つを推定するステップをさらに含む、請求項１から９のいずれか一項に記載の方法。

【請求項11】

前記端末において、地上ベースのソースまたは宇宙ベースのソースから１つまたは複数のビーコン信号を受信するステップであって、該ビーコン信号が、ビーコン位置と、基準に対する全世界タイミングとを含み、前記端末が、各ビーコン信号に対するビーコン・ドップラーまたはビーコン・ドップラー率のうちの１つまたは複数を推定し、前記衛星への１つまたは複数の送信内に該１つまたは複数の推定値を含む、受信するステップをさらに含み、前記端末のロケーションを推定するステップが、各ビーコン信号に対するビーコン時間遅延、ビーコン・ドップラー周波数、またはビーコン・ドップラー周波数変化率のうちの少なくとも１つを推定するステップをさらに含む、請求項１から９のいずれか一項に記載の方法。

【請求項12】

前記端末位置を推定するステップが、時間ｔにおける前記端末の位置ベクトルの仮説に基づく前記衛星の経時変化する位置ベクトルのための非線形最適化プロセスを使用して実行される、請求項１から１１のいずれかに記載の方法。

【請求項13】

前記衛星のフットプリントを推定するステップと、前記端末の前記位置ベクトルを該フットプリント内に存在するように拘束するステップとをさらに含む、請求項１２記載の方法。

【請求項14】

前記端末が、ＧＰＳ Ｌ１信号を受信し、Ｃ／Ａコードの「位相」（遅延）および前記ドップラー周波数の粗推定を実行し、前記衛星への１つまたは複数の送信内で該Ｃ／Ａコード位相の推定値および前記ドップラー周波数の推定値を提供するように構成され、前記端末位置を推定する前記ステップが、ＧＰＳアルマナックを受信するステップと、該Ｃ／Ａコードの受信された前記推定値および前記ドップラー周波数の受信された前記推定値を使用して、前記端末の前記位置を推定するステップとをさらに含む、請求項１から１３のいずれかに記載の方法。

【請求項15】

前記端末が、
信号を受信し、複素ベースバンド信号にダウンコンバートし、
該複素ベースバンド信号をサンプリングして、Ｗ個の複素サンプルのシーケンスを生じさせ、
該Ｗ個の複素サンプルの直列並列変換を実行し、該Ｗ個の複素サンプルにＭ個の複素正弦曲線を並列に乗算して、各々長さＮのサンプルのＭ個のシーケンスを生じさせ、使用される該Ｍ個の正弦曲線の周波数が、ドップラー周波数の予想範囲をカバーするように選択され、
サンプルの該Ｍ個のシーケンスの各々に対するゼロ・パディング離散フーリエ変換を実行して、Ｎ個の複素周波数領域サンプルのベクトルを生じさせ、
それぞれのＣ／ＡコードのＮ点ゼロ・パディング離散フーリエ変換である、該ベクトルの各々に長さＮのベクトルを要素ごとに乗算するコード位相推定器ごとのモジュールの各々に、該Ｍ個のベクトルを提供し、該乗算後の大きさが最大の要素に対応する関連付けられた前記コード遅延および前記ドップラー周波数を出力すること
によって、前記Ｃ／Ａコードの推定値および前記ドップラー周波数の推定値を取得する、請求項１４記載の方法。

【請求項16】

衛星通信システム内の地上ベースの端末の位置の推定のために地上局内に含まれ、または該地上局に動作可能に接続された通信受信機であって、
端末から衛星への送信の時間遅延、ドップラー周波数、またはドップラー周波数変化率のうちの１つまたは複数を、該衛星からの、地上局によって受信された該送信に関する情報から推定するように構成されたデコーダであって、該送信に関する該情報が、少なくとも該衛星によって受信された信号、または該衛星によって受信された信号の１つもしくは複数のデジタル・サンプルを含む、デコーダと、
該送信の時間における該衛星の位置の推定値を受信するためのインターフェースと、
該衛星の位置の該推定値と、該時間遅延の該推定値、該ドップラー周波数の該推定値、または該ドップラー周波数変化率の該推定値のうちの少なくとも１つとを使用して端末位置を推定するように構成されたプロセッサを備える位置推定モジュールと
を備える通信受信機。

【請求項17】

前記送信に関する前記情報が、オン・ボード時間基準またはオン・ボード周波数基準に対する、前記送信が前記衛星によって受信されたタイム・スタンプをさらに含む、請求項１６記載の通信受信機。

【請求項18】

前記送信に関する前記情報が、前記端末における時間基準または周波数基準に対する、前記送信が前記端末によって送信されたタイム・スタンプをさらに含む、請求項１６または１７記載の通信受信機。

【請求項19】

前記送信に関する前記情報が、前記端末内に含まれる時間基準または周波数基準に対して決定された所定の送信スケジュールに基づいた、前記端末による推定送信時間をさらに含む、請求項１６から１８のいずれか一項に記載の通信受信機。

【請求項20】

前記端末が、前記衛星に送信のシーケンスを送信し、前記衛星が、該送信のシーケンスに関する情報を前記地上局に送信し、前記デコーダが、該送信のシーケンス内の各送信に対する遅延、ドップラー、またはドップラー率のうちの１つまたは複数を推定するように構成され、前記位置推定モジュールが、該送信のシーケンス内の各送信に対する前記各送信の時間における前記衛星の位置の前記推定値と、送信の前記時間遅延の前記推定値、前記ドップラー周波数の前記推定値、または前記ドップラー周波数変化率の前記推定値の各々とを使用する、請求項１６から１９のいずれか一項に記載の通信受信機。

【請求項21】

前記位置推定モジュールが、前記送信のシーケンス内の送信の前記時間遅延、前記ドップラー周波数、または前記ドップラー周波数変化率の推定値間の差を使用して、前記端末位置を異なるように推定するように構成される、請求項２０記載の通信受信機。

【請求項22】

前記位置推定モジュールが、全世界時間基準または全世界周波数基準に対する端末オフセットを推定し、前記端末の位置を推定する際に該端末オフセットの該推定値を使用するようにさらに構成される、請求項１６から２１のいずれか一項に記載の通信受信機。

【請求項23】

前記位置推定モジュールが、衛星位置を推定し、前記端末の位置を推定する際に前記衛星位置の推定値を使用するようにさらに構成される、請求項１６から２１のいずれか一項に記載の通信受信機。

【請求項24】

前記デコーダが、前記地上局を介して地上ベースのソースまたは宇宙ベースのソースから前記衛星によって受信された１つまたは複数のビーコン信号のサンプルから、ビーコン時間遅延、ビーコン・ドップラー周波数、またはビーコン・ドップラー周波数変化率のうちの少なくとも１つを推定するようにさらに構成され、
前記位置推定モジュールが、前記端末位置を推定する際に該ビーコン時間遅延、該ビーコン・ドップラー周波数、または該ビーコン・ドップラー周波数変化率のうちの該少なくとも１つを使用するようにさらに構成される、請求項１６から２３のいずれか一項に記載の通信受信機。

【請求項25】

前記インターフェースが、前記地上局を介して、地上ベースのソースまたは宇宙ベースのソースから前記端末によって受信された１つまたは複数のビーコン信号の各々に対するビーコン時間遅延、ビーコン・ドップラー周波数、またはビーコン・ドップラー周波数変化率のうちの１つまたは複数の推定値を受信するようにさらに構成され、
前記位置推定モジュールが、前記端末位置を推定する際に該ビーコン時間遅延、該ビーコン・ドップラー周波数、または該ビーコン・ドップラー周波数変化率のうちの該少なくとも１つを使用するようにさらに構成される、請求項１６から２３のいずれか一項に記載の通信受信機。

【請求項26】

前記位置推定モジュールが、時間ｔにおける前記端末の位置ベクトルの仮説に基づく前記衛星の経時変化する位置ベクトルのための非線形最適化プロセスを使用して前記端末位置を推定するようにさらに構成される、請求項１６から２５のいずれか一項に記載の通信受信機。

【請求項27】

前記位置推定モジュールが、前記衛星のフットプリントを推定するようにさらに構成され、前記端末位置を推定するステップが、前記衛星の該フットプリントの推定値を使用して、該フットプリント内にあるように前記端末の前記位置ベクトルを拘束する、請求項２６記載の通信受信機。

【請求項28】

前記インターフェースが、ＧＰＳアルマナックを受信するようにさらに構成され、前記地上局を介して、前記端末によって受信されたＧＰＳ Ｌ１信号からＣ／Ａコードおよびドップラー周波数の推定値を受信するようにさらに構成され、前記位置推定モジュールが、前記端末位置を推定する際に該Ｃ／Ａコードの前記推定値および該ドップラー周波数の前記推定値を使用するようにさらに構成される、請求項１６から２７のいずれか一項に記載の通信受信機。

【請求項29】

前記端末が、ＧＰＳ Ｌ１信号を受信し、前記Ｃ／Ａコードの「位相」（遅延）および前記ドップラー周波数の粗推定を実行し、前記衛星への１つまたは複数の送信内で前記Ｃ／Ａコード位相の前記推定値および前記ドップラー周波数の前記推定値を提供するように構成され、前記端末位置を推定する前記ステップが、ＧＰＳアルマナックを受信するステップと、前記Ｃ／Ａコードおよび前記ドップラー周波数の受信された前記推定値を使用して、前記端末の前記位置を推定するステップをさらに含む、請求項１６から２７のいずれか一項に記載の通信受信機。

【請求項30】

衛星通信システム内の端末であって、
地上ベースのソースまたは宇宙ベースのソースから１つまたは複数のビーコン信号を受信するように構成されたＲＦトランシーバ・フロント・エンドと、
該受信された１つまたは複数のビーコン信号からビーコン位置および基準に対する全世界タイミングを取得し、各ビーコン信号に対するビーコン・ドップラーまたはビーコン・ドップラー率のうちの１つまたは複数を推定するように構成されたデコーダと、
各ビーコン信号に対する該推定されたビーコン時間遅延、ビーコン・ドップラー周波数、またはビーコン・ドップラー周波数変化率のうちの少なくとも１つを送信するように構成された送信機と
を備える端末。

【請求項31】

ＧＰＳ Ｌ１信号を受信し、該受信されたＧＰＳ Ｌ１信号のＣ／Ａコードの「位相」（遅延）およびドップラー周波数の粗推定を実行するように構成された、衛星通信システム受信機内の端末であって、
ＲＦトランシーバ・フロント・エンドであって、
ＧＰＳ Ｌ１信号を受信し、複素ベースバンド信号にダウンコンバートし、
該複素ベースバンド信号をサンプリングして、Ｗ個の複素サンプルのシーケンスを生じさせる
ように構成されたＲＦトランシーバ・フロント・エンドと、
該Ｗ個の複素サンプルをＭ個のブロックに変換するように構成された直列並列変換器と、
各々長さＮのサンプルのＭ個のシーケンスを生じさせるために、Ｗ個の複素サンプルの各ブロックにＭ個の複素正弦曲線を並列に乗算することであって、使用される該Ｍ個の正弦曲線の周波数が、ドップラー周波数の予想範囲をカバーするように選択される、乗算することと、
該サンプルのＭ個のシーケンスの各々に対するゼロ・パディング離散フーリエ変換を実行して、Ｎ個の複素周波数領域サンプルのベクトルを生じさせることと、
衛星捕捉ユニットであって、
各々がＧＰＳ衛星のＣ／Ａコードを用いて構成された複数の衛星単位コード位相推定器であって、該複数の衛星単位コード位相推定器の各々が、Ｎ個の複素周波数領域サンプルのＭ個のベクトルのうちの１つを受信し、それぞれのＣ／ＡコードのＮ点ゼロ・パディング離散フーリエ変換である、長さＮのベクトルによる該ベクトルの各々の要素ごとの乗算を実行するように構成される、複数の衛星単位コード位相推定器と、
該乗算後の大きさが最大の要素に対応するコード遅延およびドップラー周波数を決定および出力するように構成されたピーク検出器と
を備える衛星捕捉ユニットと
を備える端末。

【請求項32】

地上局と、
複数の端末と、
複数の衛星であって、各衛星が、端末から信号を受信するように、および、該衛星によって受信された該信号または該衛星によって端末から受信された信号の１つもしくは複数のデジタル・サンプルのどちらかを該地上局に再送信するように構成される、複数の衛星と、
該地上局内に含まれる、または該地上局に動作可能に接続された、請求項１６から２９のいずれか一項に記載の通信受信機と
を備える衛星通信システム。

【請求項33】

各端末が請求項３０または３１に従って構成される、請求項３２記載の衛星通信システム。

【請求項34】

地上局と、
複数の端末と、
複数の衛星であって、各衛星が、端末から信号を受信するように、および、該衛星によって受信された該信号または該衛星によって端末から受信された信号の１つもしくは複数のデジタル・サンプルのどちらかを該地上局に再送信するように構成される、複数の衛星と、
該地上局内に含まれる、または該地上局に動作可能に接続された通信受信機であって、
端末から衛星への送信の時間遅延、ドップラー周波数、またはドップラー周波数変化率のうちの１つまたは複数を、該衛星からの、該地上局によって受信された該送信に関する情報から推定するように構成されたデコーダであって、該送信に関する該情報が、少なくとも該衛星によって受信された信号、または該衛星によって受信された信号の１つもしくは複数のデジタル・サンプルを含む、デコーダと、
該送信の時間における該衛星の位置の推定値を受信するためのインターフェースと、
該衛星の位置の該推定値と、該時間遅延、該ドップラー周波数、または該ドップラー周波数変化率の該推定値のうちの少なくとも１つとを使用して端末位置を推定するように構成されたプロセッサを備える位置推定モジュールと
を備える通信受信機と
を備える衛星通信システム。

【請求項35】

請求項１から１５のいずれか一項に記載の衛星通信システムにおける地上ベースの端末の位置の推定のための方法を実行するようにプロセスを構成するための命令を含む、プロセッサ可読媒体。

【発明の詳細な説明】

【技術分野】

【0001】

優先権書類
本出願は、２０１６年５月２０日に出願された「Ｐｏｓｉｔｉｏｎ Ｅｓｔｉｍａｔｉｏｎ ｉｎ ａ Ｌｏｗ Ｅａｒｔｈ Ｏｒｂｉｔ Ｓａｔｅｌｌｉｔｅ Ｃｏｍｍｕｎｉｃａｔｉｏｎｓ Ｓｙｓｔｅｍ」という名称のオーストラリア仮特許出願第２０１６９０１９１３号の優先権を主張するものであり、その内容は、その全体が参照により本明細書に組み込まれる。

【0002】

本開示は、衛星通信システムにおける地上ベースのトランシーバの位置の推定に関する。特定の形態では、本開示は、端末から離れて推定が実行されるシステムに関する。

【背景技術】

【0003】

本出願は、無線信号を使用してデータ１２を１つまたは複数の低地球軌道周回衛星２０に送信する、端末１０と呼ばれる地上ベースのトランシーバからなる衛星通信システム１を対象とする。これらの衛星は、信号２２を送信して、端末の無線信号を地上局３０へ中継し、地上局３０において、受信された信号２２が復号され、宛先に配信される。

【0004】

対象となるいくつかの適用例では、端末１０のロケーション（たとえば、緯度、経度および高度）の推定値を地上局３０において取得することが望ましい。例示的な適用例としては、家畜、産業用機械、出荷容器、パレット、および環境センサなどの、遠隔にあり、潜在的に可動性である、資産の追跡がある。これらは、資産またはセンサそれ自体がそれ自体の位置を知る必要はない適用例である。むしろ、中央管理システムが、すべての資産／センサのロケーションを決定することを望む。そのような適用例は、端末それ自体がそれ自体のロケーションの推定値を必要とするナビゲーション・アプリケーションと対比可能である。

【0005】

そのようなシステムでは、広範囲にわたる長期間の使用を可能にするために、端末のコスト、複雑さ、および電力要件を可能な限り低く保つことが望ましいことが多い。端末１０は、それらのロケーションを（ＬＥＯ衛星２０を介して）地上局３０に送信するために使用可能なＧＰＳ受信機を取り付け可能であるが、ＧＰＳモジュール／受信機を追加することによって、コストと電力要件の両方が増す。ＧＰＳモジュールは、位置を推定するために著しい計算リソースを必要とし、加えて、ＧＰＳモジュールは、最新の衛星エフェメリス・データを有さなければならない。これによって、一般的に、ＧＰＳモジュールが、周期的にウェイク・アップし、最新のＧＰＳアルマナックを獲得することが要求される。コールド・スタートでは、ＧＰＳ受信機は、端末に対する電力要件に著しい負荷をかける完全なアルマナックをダウンロードするのに１２分かかる。

【0006】

支援型ＧＰＳシステムは、最近のアルマナック情報を端末に送信することによって、および／または端末がＧＰＳ信号をサンプリングし、より能力の高い処理デバイスにサンプルを送信することを可能にすることによって、端末にかかる計算負荷または電力要件を減少させるために使用可能である。しかしながら、これには、端末が、端末とより能力の高い処理デバイスとの間の通信チャネルを通じてサンプルを送信するのに十分な能力をもつ通信リンクを確立することが可能であることが必要とされる。上記で説明されたような遠隔端末１０に関する適用例では、この要件によって、著しいコストおよび電力要件が端末１０に追加され得る。

【発明の概要】

【発明が解決しようとする課題】

【0007】

したがって、端末内でのＧＰＳモジュールの必要性を回避する、衛星通信システムにおける地上ベースの端末の位置の遠隔推定のためのシステムおよび方法を提供すること、または少なくとも既存のシステムおよび方法に有用な代替策を提供することが必要とされている。

【課題を解決するための手段】

【0008】

第１の態様によれば、衛星通信システム内の地上ベースの端末の位置の推定のための方法であって、
端末から離れたロケーションにおいて、端末から衛星への送信の時間遅延、ドップラー周波数、またはドップラー周波数変化率のうちの１つまたは複数を、送信に関する情報から推定することであって、この送信に関する該情報は、少なくとも衛星によって受信された信号、または衛星によって受信された信号の１つもしくは複数のデジタル・サンプルを含む、推定することと、
送信の時間における衛星の位置の推定値を受信することと、
衛星の位置の推定値と、時間遅延の推定値、ドップラー周波数の推定値、またはドップラー周波数変化率（遅延、ドップラー、またはドップラー率）の推定値のうちの少なくとも１つとを使用して、端末位置を推定することと
を含む方法が提供される。

【0009】

一形態では、情報は、衛星から地上局への１つまたは複数の送信を介して受信され、端末と衛星との間の送信に関する情報から遅延、ドップラー、またはドップラー率のうちの１つまたは複数を推定することは、地上局に動作可能に接続された通信受信機によって実行され、遅延、ドップラー、またはドップラー率は、送信におけるユーザ・データの推定値を回復するために通信受信機によって復号プロセス中に推定される。情報は、地上局に送信された送信信号の１つまたは複数のデジタル・サンプルであってよく、または、情報は、地上局に再送信される衛星（すなわち、ベント・パイプ型衛星）によって受信される再送信されたアナログ信号であってよい。

【0010】

一形態では、端末と衛星との間の送信に関する情報は、オン・ボード時間基準もしくはオン・ボード周波数基準に対する、送信が衛星によって受信されたタイム・スタンプ、および／または時間基準もしくは周波数基準に対する、送信が端末によって送信されたタイム・スタンプをさらに含む。一形態では、端末と衛星との間の送信に関する情報は、端末の時間基準または周波数基準に対して決定された所定の送信スケジュールに基づいた、端末による推定送信時間をさらに含む。

【0011】

一形態では、端末は衛星に送信のシーケンスを送信し、衛星は、送信のシーケンスに関する情報を地上局に送信し、遅延、ドップラー、またはドップラー率のうちの１つまたは複数を受信情報から推定することは、送信のシーケンス内の各送信に対して実行され、端末位置を推定することは、各送信の時間における衛星の位置の推定値と、送信のシーケンス内の各送信に対する推定値の各々とを使用する。

【0012】

さらなる形態では、端末位置を推定することは、送信のシーケンス内の送信間の推定値間の差に対して異なって実行される。これは、オン・ボード時間基準またはオン・ボード周波数基準がいかなる全世界時間基準もしくは全世界周波数基準にも同期されない、および／または衛星がオン・ボード位置基準を含まない、および／または端末がいかなる全世界時間基準または全世界周波数基準にも同期されない場合に、実行され得る。

【0013】

一形態では、端末は、安定したクロックを備え、いかなる全世界時間基準または全世界周波数基準にも同期されず、端末位置を推定することは、全世界時間基準または全世界周波数基準に対する端末オフセットを推定することをさらに含む。

【0014】

一形態では、衛星は、オン・ボード位置基準を含まず、端末位置の推定は、衛星位置を推定することをさらに含む。

【0015】

一形態では、方法は、
衛星において、地上ベースのソースまたは宇宙ベースのソースから１つまたは複数のビーコン信号を受信することであって、このビーコン信号は、ビーコン位置と、基準に対する全世界タイミングとを含み、衛星は、信号のサンプルを地上局に転送する、受信することをさらに含み、端末のロケーションを推定することは、各ビーコン信号に対するビーコン遅延、ビーコン・ドップラー、またはビーコン・ドップラー率のうちの少なくとも１つを推定することをさらに含む。この推定値は、衛星位置の推定値と、したがって端末位置とを改善するために使用可能である。

【0016】

一形態では、方法は、
端末において、地上ベースのソースまたは宇宙ベースのソースから１つまたは複数のビーコン信号を受信することであって、このビーコン信号は、ビーコン位置と、基準に対する全世界タイミングとを含み、端末は、各ビーコン信号に対するビーコン・ドップラーまたはビーコン・ドップラー率のうちの１つまたは複数を推定し、衛星への１つまたは複数の送信内に推定値を含む、受信することをさらに含み、端末のロケーションを推定することは、各ビーコン信号に対するビーコン遅延、ビーコン・ドップラー、またはビーコン・ドップラー率のうちの少なくとも１つを推定することをさらに含む。

【0017】

一形態では、位置推定は、時間ｔにおける端末の位置ベクトルの仮説に基づく衛星の経時変化する位置ベクトルのための非線形最適化プロセスを使用して実行される。

【0018】

さらなる形態では、方法は、衛星のフットプリントを推定することと、端末の位置ベクトルをこのフットプリント内に存在するように拘束することをさらに含む。

【0019】

一形態では、端末は、ＧＰＳ Ｌ１信号を受信し、Ｃ／Ａコードの「位相」（遅延）およびドップラー周波数の粗推定を実行し、衛星への１つまたは複数の送信内でＣ／Ａコードの推定値およびドップラー周波数の推定値を提供するように構成され、端末位置を推定することは、ＧＰＳアルマナックを受信することと、Ｃ／Ａコードの受信された推定値およびドップラー周波数の受信された推定値を使用して、端末の位置を推定することをさらに含む。

【0020】

さらなる形態では、端末は、
信号を受信し、複素ベースバンド信号にダウンコンバートし、
この複素ベースバンド信号をサンプリングして、Ｗ個の複素サンプルのシーケンスを生じさせ、
このＷ個の複素サンプルの直列並列変換を実行し、Ｗ個の複素サンプルにＭ個の複素正弦曲線を並列に乗算して、各々長さＮのサンプルのＭ個のシーケンスを生じさせ、使用されるＭ個の正弦曲線の周波数は、ドップラー周波数の予想範囲をカバーするように選択され、
サンプルのＭ個のシーケンスの各々に対するゼロ・パディング離散フーリエ変換を実行して、各ブロックに対するＮ個の複素周波数領域サンプルのベクトルを生じさせ、
それぞれのＣ／ＡコードのＮ点ゼロ・パディング離散フーリエ変換である、ベクトルの各々に長さＮのベクトルを要素ごとに乗算するコード位相推定器ごとのモジュールの各々に、Ｍ個のベクトルを提供し、乗算後の大きさが最大の要素に対応する関連付けられたコード遅延およびドップラー周波数を出力すること
によって、Ｃ／Ａコードの推定値およびドップラー周波数の推定値を取得する。

【0021】

第２の態様によれば、衛星通信システム内の地上ベースの端末の位置の推定のために地上局内に含まれ、またはこの地上局に動作可能に接続された通信受信機であって、
端末から衛星への送信の時間遅延、ドップラー周波数、またはドップラー周波数変化率のうちの１つまたは複数を、衛星からの、地上局によって受信された送信に関する情報から推定するように構成されたデコーダであって、この送信に関する情報は、少なくとも衛星によって受信された信号、または衛星によって受信された信号の１つもしくは複数のデジタル・サンプルを含む、デコーダと、
送信の時間における衛星の位置の推定値を受信するためのインターフェースと、
衛星の位置の推定値と、遅延の推定値、ドップラーの推定値、またはドップラー率の推定値のうちの少なくとも１つとを使用して端末位置を推定するように構成されたプロセッサを備える位置推定モジュールと
を備える通信受信機が提供される。

【0022】

さらなる形態では、通信受信機は、第１の態様の方法のさらなる形態を実行するようにさらに構成されてよい。通信システムおよび第１の態様の方法を実行するように構成された端末も提供されてよい。方法は、第１の態様の方法を実行するようにプロセスを構成するためにプロセッサ可読媒体内にも提供されてよい。

【0023】

本開示の実施形態について、添付の図面を参照しながら説明する。

【図面の簡単な説明】

【0024】

【図1】一実施形態による、低地球軌道衛星通信システムの概略図である。

【0025】

【図2】一実施形態による、ドップラー推定値および遅延推定値からの位置推定のためのシステムの概略図である。

【0026】

【図3】一実施形態による、ハイブリッド分散位置決定のためのモジュールの概略図である。

【0027】

【図4】一実施形態による、粗Ｃ／Ａコード位相／ドップラー推定モジュールの概略図である。

【0028】

【図5】一実施形態による、衛星通信システムにおける地上ベースの端末の位置の推定のための方法のフローチャートである。

【発明を実施するための形態】

【0029】

以下の説明では、同じ参照文字は、図全体を通して同じまたは対応する部分を指す。

【0030】

図１は、広範なシステム構成要素および信号経路を示す、一実施形態による、低地球軌道衛星通信システムの概略図である。衛星通信システム１は、無線信号を使用してデータ１２を１つまたは複数の低地球軌道周回衛星２０に送信する、端末１０と呼ばれる地上ベースのトランシーバからなる。これらの衛星は、受信信号２２が復号されて宛先に配信される地上局３０に端末の無線信号を中継するために、信号２２を送信する。図２は、本明細書において説明するようなドップラー推定値および遅延推定値からの位置推定を実行するためのさまざまなモジュールを示す、このシステムのより詳細な概略図である。

【0031】

システムおよび方法の実施形態は、低地球軌道周回（ＬＥＯ）衛星２０によって見られるように、飛行時間（端末から衛星までの送信経路長）、相対的ドップラー周波数（経路長の一次導関数）、および端末の送信信号１２のドップラー周波数（経路長の二次導関数）の変化率から得られた測定値から、ロケーション推定値を取得する。飛行時間とは、信号の送信と受信の間の時間遅延であり、したがって、遅延（または時間遅延）と呼ばれる。同様に、相対的ドップラー周波数およびドップラー周波数の変化率は、ドップラーおよびドップラー率と呼ばれる。方法の実施形態は、任意選択で、そのような測定値のシーケンスを使用可能であり、その場合、これらの測定値は、絶対的な時間基準または周波数基準を参照する必要はなく、特定の選定された基準（たとえば第１の測定値の時間および周波数）に対して異ならせることができる。

【0032】

端末と衛星との間の（またはこれらからの）送信に関する情報は、（衛星において受信される）送信信号の１つまたは複数のデジタル・サンプルであってよく、この情報は、次いで、地上局に送信される。または、情報は、地上局に再送信される衛星（すなわち、ベント・パイプ型衛星）によって受信される再送信されたアナログ信号であってもよい。情報（または信号、またはサンプル）は、次いで、遅延、ドップラー、および／またはドップラー率の推定を含む推定プロセス（以下で詳細に説明する）を実行する通信受信機に提供される。

【0033】

端末と衛星との間の送信に関する情報は、オン・ボード時間基準またはオン・ボード周波数基準に対する、送信が衛星によって受信された時間のタイム・スタンプも含んでよい。以下で説明するように、この時間基準または周波数基準は、全世界基準に同期可能であり、非同期だが安定した時間基準または周波数基準であってもよい。同様に、情報は、端末における時間基準または周波数基準に対して送信が端末によって送信されたタイム・スタンプを含んでよい（すなわち、端末は、送信するときに、このタイム・スタンプを追加する）。同じく、これは、全世界基準に同期されてもよいし、非同期であってもよい。一実施形態では、端末は、その送信にタイム・スタンプを与えず、代わりに、端末の時間基準または周波数基準に対して決定された所定の送信スケジュールに基づいて送信する。たとえば、これは、毎秒であってもよいし、何らかの他の固定時間間隔（たとえば１００ｍｓ、２５０ｍｓ、５００ｍｓ、２ｓ、５ｓ）であってもよい。この場合、端末スケジュールは、地上局における位置推定モジュールによって知られている。タイム・スタンプ情報は、送信中に埋め込まれ、ひとたび送信が復号および回復されると利用可能になってもよいし、または衛星から地上局への追加送信すなわち側（ｓｉｄｅ）送信として提供されてもよい。本明細書の文脈では、端末と衛星との間の送信に関する情報は、信号または信号のサンプルと、この信号またはサンプルから回復可能である、または信号とともに送信される、メタデータとを含む。

【0034】

地上局３０（または地上局３０と通信するサイト）における位置推定器４０は、地上局３０において計算される端末の遅延／ドップラー／ドップラー率の推定値を、端末のロケーションに関する特定の仮説に対して予想されるであろう対応する仮説化された量と比較することによって、機能する。これらの量は、衛星の位置２０を知れば決定可能であり、その正確な推定値は、地上局に利用可能（別のソースから取得可能、または利用可能／取得可能データに基づいて推定可能、のどちらか）である。ロケーション推定値は、推定された量から仮説化された量の間の最小総合誤差をもつ仮説である。一般的には、遅延、ドップラー、およびドップラー率の推定は、これらの量における誤差の推定も含む。この誤差（または不確定度）は、位置推定器にも渡され、誤差に基づいて推定値を適切に重み付けすることなどによって、位置推定プロセスにおいて使用可能である。遅延、ドップラー、またはドップラー率のうちの１つのみが必要とされるが、性能は、一般的には、２つの値または３つの値すべての推定値が提供される場合に改善される。

【0035】

より一般には、遅延、ドップラー、またはドップラー率の推定が衛星において実行されてよいことが諒解されるであろう。さらに、方法全体は、十分な計算リソースおよび電力リソースを有するならば、衛星において実行されてもよい（遅延、ドップラー、またはドップラー率の推定、次いで、この情報を使用した位置推定）。これによって、一実施形態による衛星通信システムにおける地上ベースの端末の位置の推定のための方法のフローチャート５００である図５に示される一般的な方法がもたらされる。この方法は、
ステップ５１０において、端末から離れたロケーションにおいて、端末から衛星への送信の時間遅延、ドップラー周波数、またはドップラー周波数変化率のうちの１つまたは複数を、送信に関する情報から推定することと、
ステップ５２０において、送信の時間における衛星の位置の推定値を受信することと、
ステップ５３０において、衛星の位置の推定値と、時間遅延、ドップラー周波数、またはドップラー周波数変化率の推定値のうちの少なくとも１つとを使用して、端末位置を推定することとを含む。送信に関する情報は、衛星によって受信された信号、または衛星によって受信された信号の１つもしくは複数のデジタル・サンプルを少なくとも含み、送信タイム・スタンプまたは受信タイム・スタンプまたは推定送信時間などのさらなる情報を含んでよい。したがって、この方法は、衛星内で、衛星と地上局における位置推定モジュールとの間で共同で、もっぱら地上局内で（すなわち、地上局は、遅延／ドップラー／ドップラー率の推定を実行するための通信受信機と、位置推定器とを含む）、または地上局から離れて配置されているが（たとえば通信リンクを経由して）地上局に動作可能に接続された位置推定モジュールに遅延の推定値、ドップラーの推定値、もしくはドップラー率の推定値を提供する地上局（通信受信機を含む）間で、実施可能である。

【0036】

方法の実施形態は、加えて、端末によって受信された他のビーコン信号５０を任意選択で機会主義的に（ｏｐｐｏｒｔｕｎｉｓｔｉｃａｌｌｙ）活用することができる。これらのビーコン信号５０は、ビーコン衛星５１からのビーコン信号５２であってもよいし、全地球測位システム（ＧＰＳ）衛星５３からのＧＰＳ Ｌ１信号５４であってもよいし、または地上ビーコン５５からのビーコン信号５６であってもよい。本発明者らは、すべてのそのような信号５０を「ビーコン」と呼ぶものとする。そのようなビーコン信号５０が役立つために、
信号のソースのロケーションが、地上局では既知である、または推定可能であるが、端末では必ずしもそうではないこと、および
端末が、受信信号から遅延、ドップラー、またはドップラー変化率のうちの１つまたは複数を推定することができること
が必要とされる。

【0037】

この場合、端末１０が、ビーコン信号（または複数のビーコン信号）のための遅延／ドップラー／ドップラー率推定値のうちの１つまたは複数を計算し、または別の方法で推定し、次いで、そのようなビーコン信号５０の推定値をそのデータ送信１２の一部として低地球軌道衛星２０に送信する。次いで、位置推定器が、端末により計算された推定値を対応する仮説値と比較して、端末位置推定値を改良する。このプロセスについては、以下のセクション２でより詳細に説明する。

【0038】

いくつかの例では、低地球軌道衛星のロケーションは、地上局では演繹的に正確に知られていない場合がある。そのような場合、地上ビーコンは、任意選択で、地上局の衛星位置の推定値を改善するために活用される。この場合、低地球軌道衛星は、加えて、地上ビーコン信号をサンプリングし、これらのサンプルを地上局に中継する。次いで、地上局が、これらの信号の遅延／ドップラー／ドップラー率の推定値を使用して、衛星の位置を推定する、またはその衛星の位置の推定値を改善することができる。地上ビーコンは、信頼された既知のロケーションを有する、端末の一般集団の特殊な例であるかもしれない（たとえば、地上ビーコンは、自らの正確な位置を得て、それを地上局に送信する標準的なＧＰＳモジュールを搭載してよい）。あるいは、ビーコンは、既知のロケーションから送信された、機会の信号であるかもしれない（たとえば、地上ラジオ放送またはテレビジョン放送）。

【0039】

提案される測位方法は、遠隔追跡アプリケーションのためのいくつかの利点を有する。
１．この方法は、端末におけるＧＰＳモジュールを必要としない。これによって、端末は組み立てるためのコストが低くなり、はるかに少ない電力を使用するようになり、より長いバッテリ寿命をもたらす。
２．この方法は、端末が、ＧＰＳ衛星、ビーコン衛星、または地上ビーコンの位置の知識を取得することを要求しない。たとえば、コールド・スタート・モードでは、ＧＰＳ受信機が、完全なアルマナックをダウンロードするのに１２分かかる。本明細書において説明する方法の実施形態は、このプロセスを回避する。
３．この方法は、主に、地上局において実施され、地上局では、より多くの計算リソースが利用可能であり、電力使用はそれほど重要ではなく、ＧＰＳ衛星および他の衛星のためのエフェメリス・データは、たとえばインターネット・データベースへの直接接続を介して、利用可能である。

【0040】

本明細書において説明する方法は、非同期端末と同期端末の両方に適用可能である。
非同期：端末は、いかなる全世界時間基準または全世界周波数基準に同期されない。しかしながら、端末は、それ自体のローカルな時間基準および周波数基準を得ることができる、それ自体の十分に安定したクロックを有する。これは、端末の時間および周波数は、パケットによってあまり大きく変動しないことを意味する。
時間および／または周波数同期：端末は、全世界基準に正確に同期されている安定した時間基準および／または周波数基準を有する。

【0041】

本明細書において説明する方法はまた、オン・ボード全世界位置基準およびオン・ボード全世界タイミング基準の有無にかかわらず、低地球軌道衛星に独立して適用可能である。
位置認識：低地球軌道衛星は、宇宙におけるその位置の正確な知識および正確な全世界タイミング基準（たとえばオン・ボードＧＰＳペイロードから取得された）を有する。
位置未認識：低地球軌道衛星は、オン・ボード位置基準を有さない。加えて、低地球軌道衛星は、正確な全世界タイミング基準を有さないことがある。

【0042】

２．０ 位置推定

【0043】

２．１ 詳細なシステム・モデル

【0044】

図２は、一実施形態による、ドップラー推定値および遅延推定値からの位置推定のためのシステム２の詳細な概略図である。このシステムは、ユーザ・データを含むデータ１２のパケットのシーケンスを低地球軌道内の衛星ペイロード２４を介して地上局３０に送信する遠隔ユーザ端末１０からなる。

【0045】

時間ｔにおいて、端末と衛星との間の距離ｄ（ｔ）は
ｄ（ｔ）＝｜｜ｘ_ｔ（ｔ）−ｘ_ｓ（ｔ）｜｜ 式１
であり、ここで、ｘ_ｔ（ｔ）およびｘ_ｓ（ｔ）はそれぞれ、時間ｔにおける端末および衛星の位置である。これらの位置は、たとえば、地球中心慣性（ＥＣＩ）座標で表現可能であり、｜｜・｜｜は標準的なユークリッド距離測度である。しかしながら、他の類似の座標系も使用可能である。

【0046】

図２を参照すると、端末１０は、たとえば端末（図示せず）と統合されたもしくはこれに接続されたセンサ・モジュールからのユーザ・データ１１２、および／またはメモリに記憶され得る、一意の端末識別子番号などの識別データ１１４を受信する。パケット・フォーマッティング・モジュール１１０は、ユーザ・データ１１２および識別データ１１４を受け取り、メッセージ・ビット１１６を生成する。送信機ＩＱ波形生成器１２０は、メッセージ・ビット１１６およびスロット・タイミング情報１２２を受け取って、送信機ＩＱサンプル１２４を生成する。デジタル・アナログ（ＤＡＣ）変換器１３０は、アナログ信号をＲＦトランシーバ・フロント・エンド１４０に提供し、ＲＦトランシーバ・フロント・エンド１４０は、アンテナ１４２を介して信号１２を送信する。ＲＦトランシーバ・フロント・エンド１４０は、広帯域トランシーバ・フロント・エンドであってよい。

【0047】

低地球軌道衛星２０は、端末１０から送信された信号１２を受信するＲＦフロント・エンド２５を有する衛星ペイロード２４を含む。送信１２の飛行時間τ（ｔ）は、
τ（ｔ）＝ｄ（ｔ）／ｃ 式２
であり、ここで、ｃは光のスピードである。相対速度は、ドップラー周波数を直接的に与える

【数1】

であり、ω（ｔ）は、

【数2】

であり、ここで、ω_ｃは送信の中心周波数である。同様に、本発明者らは、ドップラー率を

【数3】

と定義することができる。

【0048】

以下で詳述するように、スケーリング定数に達するまで、基礎をなす対象となる量はｄ（ｔ）と、その一次時間導関数および二次時間導関数ｄ’（ｔ）およびｄ’’（ｔ）である。座標の選択は恣意的である。以下で、遅延、ドップラー、およびドップラー率を参照するが、方法は、距離、速度、および加速度、またはこれらの量の他の任意の可逆変換を使用しても、等しく機能するであろう。したがって、本明細書の文脈では、遅延、ドップラー、およびドップラー率の言及は、これらの量の可逆変換から取得された同等の量を包含することができる。

【0049】

端末によって送信されるパケットのシーケンスは、ｊによってインデックスが付けられる。一般性を失うことなく、ｊ＝０，１，・・・，Ｊ−１であると仮定し、ここで、Ｊは、衛星通過中に端末によって送信されたパケットの数である。衛星ペイロード２４は、（ドップラー・シフトを可能にする端末送信を捕らえるのに十分に広い帯域幅を使用する）ＲＦフロント・エンド２５と地上局３０への送信のためにダウンリンク・トランシーバ２８にスペクトル・サンプル２７を提供するアナログ・デジタル（ＡＤＣ）変換器２６とを介して、関連のある無線周波数帯域をサンプリングする。この実施形態では、端末と衛星との間の送信に関する情報は、少なくともこれらのスペクトル・サンプルを含む。加えて、情報は、追加されてもよいし地上局に別個に送信されてもよい、タイム・スタンプなどの他の任意のメタデータを含んでよい。

【0050】

衛星が、（たとえば全地球測位システムから取得される）正確な全世界時間基準へのアクセスを有する場合、衛星は、これらのデジタル化されたサンプル２７にタイム・スタンプを与え、ダウンリンク・トランシーバ２８は、異なる周波数上で異なる帯域幅であってよいデジタル通信リンク２２を使用して、それらを地上局３０に中継する。地上局３０は、これらのデジタル・サンプルを処理のために通信受信機３２に転送する。通信受信機３２は、地上局３０と同じ場所に配置されてもよいし、されなくてもよい。他の実施形態では、端末信号１２を地上局３０に中継することができる他の衛星ペイロード構成が使用可能である。衛星は、地上局が衛星２０の観測視野角の範囲内にある場合、リアル・タイムで信号１２を中継することができる、または、地上局が衛星２０の観測視野角の範囲内にあるとき、受信信号１２を後で記憶および転送することができることに留意されたい。

【0051】

通信受信機３２は、受信されたスペクトル・サンプルから端末パケットのシーケンスを復号する。復号プロセスの主要な目標は、ユーザ・データ１１２の推定値を取得または回復させることである。パケットｊに対する復号プロセスの一部として、通信受信機は、ドップラー・シフトω_０ｊ、ドップラー率ｖ_０ｊ、および遅延τ_０ｊ３４の関連付けられた推定値も生じさせる。これらの推定値３４はすべて、正確で安定した時間基準を基準とし、この時間基準は、衛星、地上局、または位置プロセッサに利用可能な同じ基準であってもよいし、同じ基準でなくてもよい（これについては、以下で詳述する）。

【0052】

下付き文字０は、これらの推定値３４を低地球軌道衛星を用いて識別するために使用され、この下付き文字は、後で、他の衛星に関連のある量を区別するために使用される（すなわち、これは、利用可能なビーコン衛星および／または地上ビーコンのリストに後でインデックスを付ける）。

【0053】

推定値のシーケンス｛ω_０ｊ，ｖ_０ｊ，τ_０ｊ｝、ｊ＝０，１，・・・，Ｊ−１（３４とラベルが付けられる）は、位置推定モジュール４０に転送され、位置推定モジュール４０は、通信受信機３２と同じ場所に配置されてもよいし、されなくてもよい。本明細書の文脈では、ターム（ｔｅｒｍ）・モジュールは、方法を実施するように構成された、特定用途向け集積回路（ＡＳＩＣ）、デジタル信号プロセッサ（ＤＳＰ）、デジタル信号処理デバイス（ＤＳＰＤ）、プログラマブル論理デバイス（ＰＬＤ）、もしくはフィールド・プログラマブル・ゲート・アレイ（ＦＰＧＡ）、または他の支持構成要素（電源、データ・バス、通信ボード、またはチップなど）とともに、方法を実施するために実行可能な命令を記憶するプロセッサおよび関連メモリ（オン・ボードまたは動作可能にリンクされる）をもつプロセッサ・ボードなどのハードウェアを包含することに留意されたい。ターム・モジュールとプロセッサは互換的に使用され、その場合、プロセッサは、プロセッサ・ユニットならびに関連するハードウェアおよびソフトウェアとして広義に解釈されるべきである。

【0054】

位置推定プロセッサ４０は、低地球軌道衛星に関する利用可能なエフェメリス・データ（軌道要素）４１を有する。これらは、たとえば、所与のエポック（時点）に関する軌道要素を符号化するごく最近に公表された２行軌道要素（ＴＬＥ）から取得される。シーケンス｛ω_０ｊ，ｖ_０ｊ，τ_０ｊ｝３４およびエフェメリス・データ４１を使用して、位置推定プロセッサ４０は、各パケットｊ＝０，１，・・・，Ｊ−１の送信の時間における端末の位置の位置推定値４２のシーケンスを生じさせる。

【0055】

時間同期動作および周波数同期動作の場合、位置推定プロセッサ４０は、ω_０ｊ，ｖ_０ｊ，τ_０ｊ３４が「絶対的」である、すなわち、既知の基準を参照しているということを活用することができる。

【0056】

非同期動作（端末または衛星のどちらかが正確な全世界時間基準を有さない）の場合、位置推定プロセッサは、代わりに、異なるパケットからのこれらの推定値の安定性を活用することができる、すなわち、異なるインデックスｊをもつパケット間で時間推定値および／または周波数推定値を比較することができる。たとえば、プロセッサは、以下の量を使用して、パケットｊ＝０に対して異なって動作することができる。

【数4】

【0057】

基準パケットｊ＝０の選択肢は恣意的であり、ｊ＝０，１，・・・，Ｊ−１の他の任意の選択肢であってよい。

【0058】

２．２ 非線形最適化を介した位置推定

【0059】

衛星エフェメリス・データ４１から、位置推定器４０は、少なくとも低地球軌道衛星の経時変化する位置ベクトルｘ_ｓ（ｔ）の近似を決定することができる。ｔ_ｊは、パケットｊ＝０，１，・・・，Ｊの送信時間とする。ｘ_ｔ（ｔ）は、（ｘ_ｓ（ｔ）と同じ座標系内の）時間ｔにおける端末の位置ベクトルに対する仮説とする。簡単にするために、さしあたり、時間間隔（ｔ_０，ｔ_Ｊ−１）に対するｘ_ｔ（ｔ）≒ｘ_ｔを仮定する。

【0060】

所与のｘ_ｓ（ｔ）の場合、衛星の「フットプリント」が決定可能である。これは、（障害がない場合）衛星が見える、地球上の位置Ｘ（ｔ）のセットである。これは、時間ｔにおける端末位置に対する選択肢の実現可能なセットを表し、すなわち、ｘ_ｔ∈Ｘ（ｔ）であることがわかっている。このフットプリントは、衛星位置、タイミング、および衛星との通信に必要とされる最小仰角における小さい不確定度について説明するように容易に修正可能である。Ｘ＝∩_ｊＸ（ｔ_ｊ）とする。

【0061】

次に、端末が同期端末である一実施形態による位置推定のための方法について説明する。次いで、セクション２．２．２において、この実施形態が、端末が非同期端末である一実施形態に一般化される。

【0062】

低地球軌道衛星２０によって受信可能な地上ビーコン信号または宇宙ベースのビーコン信号５８は、低地球軌道衛星の位置推定値を改善することによって性能を改善するために使用可能である。これについては、セクション３．１で説明する。ビーコン信号５０は、端末によっても受信され、位置推定値４２を改善するために使用可能である。これについては、セクション３．２で説明する。

【0063】

２．２．１ 同期端末

【0064】

まず第一に、端末１０が時間および周波数同期を有し、したがって、推定値ω_０ｊおよびｖ_０ｊがω_ｃに関し、τ_０ｊがパケットｊの送信の正確な時間ｔ_ｊに関する（すなわち、飛行時間の推定値τ（ｔ_ｊ）である）と仮定する。

【0065】

時間ｔにおける所与の仮説端末位置ｘ_ｔおよび衛星位置ｘ_ｓ（ｔ）の場合、
ｄ（ｔ｜ｘ_ｔ）＝｜｜ｘ_ｔ−ｘ_ｓ（ｔ）｜｜_２ 式９
を、時間ｔにおける端末と衛星との間の経時変化する距離に関する仮説とする。同様に、ω（ｔ｜ｘ_ｔ）、ｖ（ｔ｜ｘ_ｔ）、およびτ（ｔ｜ｘ_ｔ）を定義する。

【0066】

誤差関数
Δω_０ｊ（ｘ_ｔ）＝ω（ｔ_ｊ｜ｘ_ｔ）−ω_０ｊ 式１０
Δｖ_０ｊ（ｘ_ｔ）＝ｖ（ｔ_ｊ｜ｘ_ｔ）−ｖ_０ｊ 式１１
Δτ_０ｊ（ｘ_ｔ）＝τ（ｔ_ｊ｜ｘ_ｔ）−τ_０ｊ 式１２
を定義する。

【0067】

さらに、総合誤差を

【数5】

と定義し、ここで、α_０ｊ、β_０ｊ、γ_０ｊは、全世界誤差関数への寄与を重み付けするために使用される非負の定数である。これらの定数は、システム性能を調整するために使用可能であり、実験によって前もって決定可能であり、値をより同等の値に変換するための単位係数（ｕｎｉｔ ｆａｃｔｏｒ）の変更を含んでよい。これらの定数は、それらのパケットに関して報告された品質尺度により、異なるパケットからの推定値の寄与を重み付けするためにも使用可能である（たとえば、より低い信号対雑音比をもつパケットからの推定値にあまり注目しないことができる）。τ、ω、およびｖの異なる推定器は、これらの推定値の品質尺度も報告してよく、これらは、重み付け係数を選ぶために使用されてよい。

【0068】

位置推定プロセッサは、以下の非線形最適化問題を解くことによって、端末位置の推定値

【数6】

を生じさせる。

【数7】

【0069】

非線形最小二乗最適化問題の解決のための多くの既知の手法がある。一例が、Ｌｅｖｅｎｂｅｒｇ−Ｍａｒｑｕａｒｄｔアルゴリズムである。そのようなアルゴリズムは反復的に機能し、誤差関数の勾配を考慮することによって最初の推測を改善する。最初の推測のための良好な選択肢は、時間ｔにおける衛星の地上軌跡位置である。本発明者らは、領域Ｘにおける全世界的な解を必要とするので、この空間を（たとえばグリッディングによって）離散化し、Ｌｅｖｅｎｂｅｒｇ−Ｍａｒｑｕａｒｄｔなどのローカル最適化器を走らせ、各グリッド点を用いて初期化することができる。全世界的な解の候補は、結果として生じる最小のΔ（ｘ_ｔ）をもつものである。

【0070】

他の実施形態では、推定は、他の関数または最適化方法を使用して実行されることがある。たとえば、線形最小二乗法が使用されてもよいし、または第１の推定値として線形方法を使用し、この推定値をさらに改良するために非線形方法を使用する（または、その逆）ハイブリッド解決策が使用されてもよい。同様に、たとえば関数を定義し、その関数を最小または最大にすることによって誤差を最小にしようとする反復的方法を含む他の方法がある。いくつかの場合では、推定プロセスは、最適な解決策であることが保証されることが必要であるが、所定の閾値より小さい誤差、または連続した反復もしくは反復のシーケンス間の変化が所定の閾値よりも小さいことなどの、精度または安定性に基づいた停止基準を満たすによって定義されてよい。

【0071】

端末位置に関する他の任意の情報がない場合に、解決不能の曖昧さがある可能性に留意されたい。これは、同じ最小Δ（ｘ_ｔ）を有する、Ｘに属する２つの端末位置候補があることを意味する。いくつかの場合には、そのような曖昧さは、演繹的に既知であり得る端末位置についての追加情報を用いることによって解決可能である。たとえば、端末のおおよその位置は、すでに、たとえば以前に取得された位置フィックスから知られている場合があり、この情報は、不正確な曖昧さを除外するために使用可能である。

【0072】

前述のように、簡単にするために、端末は動かないと仮定された。しかしながら、いくつかの実施形態では、この仮定が緩和されることがある。その場合、式９〜式１３におけるｘ_ｔをｘ_ｔ（ｔ）で置き換え、そのため、誤差が各送信に対して推定され、最適化式１４は、固定時間と速度成分における位置ベクトルまたは位置について解く。

【0073】

２．２．２ 非同期端末

【0074】

端末が、既知の基準に対する時間および／または周波数の同期がない（ただし、対象となる持続時間にわたって別の方法で安定している）とき、式１４に類似した最適化問題を作ることが依然として可能である。必要とされる主な修正は、異なるように動作し、式６〜式８において定義された量

【数8】

に関する式１０〜式１２に類似した誤差関数を定義することである。

【0075】

３．０ ビーコン信号

【0076】

システムは、任意選択で、地上無線送信機によって、または衛星によって送信された、ビーコン信号５０を活用することができる。これらのビーコンが既知の（または推定可能な）位置を有することが必要とされる。例としては、
・ＧＰＳ衛星５３によって送信されたＬ１信号５４、または、他の全地球測位システムからの信号
・地上ＧＰＳ補強システムによって送信された信号
・既知の信頼された位置およびタイミングをもつ端末（たとえば、ＧＰＳ受信機を取り付けられた端末）からの送信
・デジタル・オーディオ・ブロードキャスト無線信号などの、既知の周波数および既知のロケーションをもつ地上ビーコン５５からの地上ブロードキャスト５６
・他の衛星コンステレーション５１から送信されたビーコン５２
がある。

【0077】

セクション３．１は、低地球軌道衛星によって受信されたそのような信号をどのように使用するかについて説明する。セクション３．２は、端末によって受信されたそのような信号をどのように活用するかについて説明する。ビーコン信号はｋ＝１，２，・・・，Ｋ−１によってインデックスが付けられるものとし、ビーコン位置はｂ_ｋとする。

【0078】

３．１ 衛星によって受信されたビーコン

【0079】

端末１０の正確な位置推定値を取得するために、位置推定器４０が衛星２０の位置の正確な知識を有することが重要である。ＴＬＥデータから得られた推定値は、実際には、十分に正確でないことがある。性能を改善するために、衛星は、加えて、宇宙ベースのソースまたは地上ソース５１、５３、５５から送信されたビーコン信号５０をサンプリングすることができる。一例は、低地球軌道衛星２０がＧＰＳ受信機ペイロードを持っており、ＧＰＳ信号を受信してその位置を非常に正確に判断し、この情報を地上局３０に中継することである。これは、そのようなＧＰＳペイロードが存在しない場合があり、または大量の電力を使用する場合があり、および利用可能な帯域幅のうちのいくらかを使用することがあるので、いくつかの例では、望ましくない場合がある。

【0080】

ＧＰＳ信号を受信および復号する代わりに、加えて、衛星は、図１に示すように、地上ビーコン端末５５（または他の任意のビーコン）からの信号５８をサンプリングすることができる。

【0081】

地上ビーコン端末５５は、自らの正確な位置および全世界タイミングを（送信し得る他の任意のデータに加えて）送信する。加えて、衛星２０が、これらの信号をサンプリングし、これらのサンプルを、位置推定器４０への提供のために地上局３０に転送する。位置推定器４０は、ビーコン信号ｋ＝０，１，・・・，Ｋ−１の各々に対して遅延τ_ｋ、ドップラーω_ｋ、およびドップラー率ｖ_ｋの推定値を取得する。これらは、さらに、各ビーコンによって送信された特定のパケットによってインデックス付与可能であるが（たとえばω_ｋｊ）、これは、明確にするため省略されている。

【0082】

セクション２．２．１における発展を一般化すると、所与の仮説端末位置ｘ_ｔおよび仮説衛星位置ｘ_ｓ（ｔ）に対して、ω（ｔ｜ｘ_ｔ，ｘ_ｓ（ｔ））、ｖ（ｔ｜ｘ_ｔ，ｘ_ｓ（ｔ））および、τ（ｔ｜ｘ_ｔ，ｘ_ｓ（ｔ））を、端末および衛星位置の特定の仮説選択肢に対応する時間ｔにおけるドップラー、ドップラー率、および遅延と定義する。同様に、ω（ｔ｜ｂ_ｋ，ｘ_ｓ（ｔ））、ｖ（ｔ｜ｂ_ｋ，ｘ_ｓ（ｔ））、およびτ（ｔ｜ｂ_ｋ，ｘ_ｓ（ｔ））が、ビーコンｋの既知の位置および衛星の仮説位置に対する仮説量であると定義する。

【0083】

時間ｔ_ｊにおいて受信されるパケットｊ＝０，１，・・・，Ｊ−１に対する誤差関数を、
Δω_０ｊ（ｘ_ｔ，ｘ_ｓ（ｔ_ｊ））＝ω（ｔ_ｊ｜ｘ_ｔ，ｘ_ｓ（ｔ_ｊ））−ω_０ｊ 式１５
Δｖ_０ｊ（ｘ_ｔ，ｘ_ｓ（ｔ_ｊ））＝ｖ（ｔ_ｊ｜ｘ_ｔ，ｘ_ｓ（ｔ_ｊ））−ｖ_０ｊ 式１６
Δτ_０ｊ（ｘ_ｔ，ｘ_ｓ（ｔ_ｊ））＝τ（ｔ_ｊ｜ｘ_ｔ，ｘ_ｓ（ｔ_ｊ））−τ_０ｊ 式１７
と定義し、時間ｔ_ｋにおいて観測されるビーコンを、
Δω_ｋ（ｂ_ｋ，ｘ_ｓ（ｔ_ｋ））＝ω（ｔ_ｋ｜ｘ_ｔ，ｘ_ｓ（ｔ_ｋ））−ω_ｋ 式１８
Δｖ_ｋ（ｂ_ｋ，ｘ_ｓ（ｔ_ｋ））＝ｖ（ｔ_ｋ｜ｂ，ｘ_ｓ（ｔ_ｋ））−ｖ_ｋ 式１９
Δτ_ｋ（ｂ_ｋ，ｘ_ｓ（ｔ_ｋ））＝τ（ｔ_ｋ｜ｂ_ｋ，ｘ_ｓ（ｔ_ｋ））−τ_ｋ 式２０
と定義する。

【0084】

総合誤差を、

【数9】

と定義し、ここで、以前と同様に、α_ｋ、βｋ、γ_ｋは、全世界誤差関数に対するビーコン誤差の寄与を重み付けするために使用される追加の非負の定数である。

【0085】

Ｓを、衛星２０の可能な位置の範囲（たとえば、現在ＴＬＥによって予測される軌道についての不確実性楕円として取得される）とする。位置推定モジュール４０は、以下の非線形最適化問題

【数10】

を解くことによって、端末位置の推定値

【数11】

を生じさせる。

【0086】

式２２の副産物は、低地球軌道衛星２０に対する正確な位置推定値である。非同期動作への拡張は、セクション２．２．２において説明したものと同様のやり方で達成され、誤差関数のさまざまなバージョンに対して動作する。方法は、低地球軌道衛星２０が正確な全世界時間基準を有さない状況にさらに拡張可能である。この場合、最適化は、衛星タイミング基準ｔに不確定度を含むように拡張される。Ｔは、このタイミング基準

【数12】

における不確定度とする。

【0087】

３．２ 端末によって受信されたビーコン

【0088】

端末によって受信されたビーコンも利用可能である。これらの実施形態では、端末は、信号の時間遅延τ_ｋ、ドップラーω_ｋ、およびドップラー率ｖ_ｋのうちの少なくとも１つをビーコンｋから推定するために、およびこれらの推定値を、低地球軌道衛星２０を介して地上局３０に送信するために必要とされる。端末においてτ_ｋ、ω_ｋ、および／またはｖ_ｋの推定を実行することによって、これは、支援型ＧＰＳシステムと同様に信号のサンプルを地上局に送信する必要性を不要にし、端末が地上局に返送する必要があるデータの量を減少させる。推定は、受信されたビーコン信号を処理および復号するデコーダを使用して実行されてよい。

【0089】

セクション３．１において説明した拡張と同様に、次いで、位置推定器が、これらの追加推定値に対する項を含むように誤差関数（式１４）を拡張し、ビーコンの既知の位置および端末ｘ_ｔの仮説化された位置が与えられれば仮説値と比較することによって、これらの推定値をその非線形最適化に組み込む。この方法は、セクション３．１において説明した方法に加えて使用可能である。唯一の違いは、τ_ｋ、ω_ｋ、およびｖ_ｋ、およびビーコンの識別情報（またはその位置）が、地上局において計算されるのではなく、端末によって計算され、地上局に送信されることである。

【0090】

セクション３．３では、衛星のＧＰＳコンステレーションによって送信されたＬ１信号から（差動）遅延およびドップラー推定値を取得するための特定の低複雑度方法について説明する。

【0091】

ＧＰＳ信号は、長さ１０２３のゴールド符号であり、２つのｍ−シーケンスから構築されるｃｏａｒｓｅ ａｃｑｕｉｓｉｔｉｏｎ Ｃ／Ａコードを備え、衛星固有である。このコードは、正確に１ｍｓの周期を有する。暗号化された精度コード（Ｐ（Ｙ））は１０．２３Ｍｃｈｉｐ／ｓであり、１週間の周期を有する（約２．３５×１０^４チップ。ナビゲーション・データは、５０ビット／ｓでＣ／Ａコード上に変調される（ビット周期２０ｍｓ、すなわち、ビットあたり２０のＣ／Ａコード反復）。ナビゲーション・データは、３０ビットのテレメトリ・データと、３０ビットのハンドオーバ・データと、ＧＰＳ日付および時間と、エフェメリス・データと、アルマナック・データとを含む。テレメトリ・データおよびハンドオーバ・データは、６秒（３００ビット）ごとに反復される。Ｃ／ＡコードおよびＰ（Ｙ）コードは、Ｐ（Ｙ）に対する３ｄＢの減衰を用いて直交変調される。

【0092】

ＧＰＳ受信機用語では、擬似距離とは、（サブフレームの開始に対する）衛星における送信の時間と受信の時間との絶対時間差である。標準的なＧＰＳ受信機処理では、受信機は、これらの擬似距離を取得するために絶対時間基準を有すると仮定される。

【0093】

標準的受信機は、
ｃｏａｒｓｅ ａｃｑｕｉｓｉｔｉｏｎを取得し（どのＣ／Ａコードが存在するかを検出し、おそらく、その検索では、以前に受信されたアルマナック・データによって支援される）、
コード追跡ループに入り、
Ｃ／Ａコード内で搬送される５０ｂｐｓナビゲーション・データを復調し、
各サブフレーム内のハンド・オーバー・ワード（ＨＯＷ）から復調されたデータから衛星において時間を決定し、
ビュー内の各衛星のための擬似距離を計算する。これによって、受信機が、安定してＧＰＳ時間に同期されているクロックを有することが必要とされる。

【0094】

衛星ｊのための擬似距離は、２つの構成要素、すなわち、（ａ）ミリ秒単位で離散化された、フレームの開始において測定された時間差（Ｃ／Ａコード期間）および（ｂ）ミリ秒未満の解像度である、ＰＬＬからのコード位相、の合計である。システムでは、ローカルで安定し、いかなる外部ソースにも同期されていない受信機クロックのみが必要とされる。ローカルで安定するとは、信号サンプルの１つのベクトルを処理する持続時間にわたって十分な正確さをもつ、それ自体の時間を維持することを意味する。絶対擬似距離を決定する代わりに、本発明者らは、差動擬似距離を計算する。差動擬似距離は、各衛星に対する関連するサブフレーム境界と最も早い衛星との相対遅延である。

【0095】

３．３ ビーコンとしてのＧＰＳ Ｌ１信号の使用

【0096】

このセクションでは、ビーコン信号５０の特定のインスタンス、すなわち、ＧＰＳ衛星５３によって送信されたＬ１信号５４について説明する。受信されたＬ１信号５４から必要とされる遅延およびドップラー推定値を取得するための位置推定方法の一実施形態について説明する。ＧＰＳ信号から位置を取得するための従来の方法に関する１つの問題は、かなりの計算リソースと、あまり古くないＧＰＳエフェメリス・データを必要とすることである。この実施形態では、そのような要件は、有利には回避可能である。

【0097】

端末におけるかなりの計算リソースという第１の問題に対処する１つの方策は、端末においてＧＰＳ信号をサンプリングし、これらのサンプルを、より能力の高い中央処理デバイスに送信することである。このアイデアは、セルラー式無線適用例において使用され、支援型ＧＰＳとして知られている。しかしながら、この手法は、端末と中央処理ユニットとの間における十分に有能な通信チャネルの存在を仮定するものである。地理的に分散された端末をもつ多くの適用例では、そのような通信チャネルまたは通信リンクは、インフラストラクチャの欠如により利用可能でない、もしくは実現可能でない、または、必要とされるリンクを確立もしくは維持するために必要とされるさらなるコスト、電力、もしくは帯域幅により禁止する、のどちらかである。

【0098】

本実施形態では、この第１の問題は、端末が端末から中央の位置推定プロセッサ４０に少数のビットを送信することのみが必要であるように、端末とより能力の高い中央の受信機との間の計算を注意深く分割することによって対処される。第２の問題に関して、ＧＰＳ衛星エフェメリス・データは、通常、「アルマナック」の送信を介して最新に保たれる。これは、ＧＰＳ受信機モジュールにとって標準的な慣例である。これによって、デバイスがこれらの信号を受信することができ、デバイスが、アルマナックを受信するのに十分に長い間、オンにされ得る（または、スリープから引き出され得る）ことが必要とされる。

【0099】

本明細書において説明する方法は、これらの要件の両方を回避する。有利には、他のソフトウェアにより定義されたＧＰＳ受信機とは異なって、この方法は、端末においてアルマナック（またはエフェメリス）・データを必要としない。その代わりに、位置推定が中央プロセッサで完了され、中央プロセッサは、標準的な通信リンクを介して（たとえばインターネットを介して）アルマナックを備えることができる。

【0100】

図３は、一実施形態による、ハイブリッド分散位置決定のためのモジュールの概略図３を示す。図２における対応する部分は、それに応じてラベルが付される。加えて、端末１０は、ＧＰＳアンテナ１４４と、コード位相およびドップラー推定器モジュール１６０にＧＰＳ ＩＱサンプルを提供するＡＤＣ１５０にベースバンド・アナログ信号１４６を提供する広帯域ＲＦトランシーバ・フロント・エンド１４０とを備える。これによって、パケット形成モジュール１１０にコード位相およびドップラー・ビットを提供する圧縮モジュール１７０によって圧縮されるコード位相およびドップラー推定値

【数13】

１６２が生成される。次いで、コード位相およびドップラー推定値

【数14】

１６２の推定値が、復号動作中に通信受信機３２によって取得され、位置推定器４０に提供される。

【0101】

３．３．１ 粗遅延およびドップラー推定

【0102】

図４は、一実施形態によりＣ／Ａコードの「位相」（遅延）およびドップラー周波数の粗推定を実行する、端末１０内の粗Ｃ／Ａコード位相／ドップラー推定モジュール４００のより詳細な概略図である。この実施形態では、直列ドップラー／並列位相手法が推定に使用されるが、他の実施形態では、直列ドップラー／並列位相手法は、その逆に行われることも可能である。

【0103】

共通フロント・エンド４１０は、信号４０２を受信し、複素ベースバンド４１２にダウンコンバートし４１１、４１２をサンプリングして、Ｗ個の複素サンプルのシーケンス４１４を生じさせる。この実施形態では、Ｗ個のサンプルは直列並列変換器４１４に渡され、次いで、フロント・エンドが、これらのサンプルに（並列に）Ｍ個の複素正弦曲線

【数15】

４１５を乗算して、サンプルのＭ個のシーケンスｘ_ｍ、ｍ＝１，２，・・・，Ｍ、４１６、各々の長さＮ４１６を生じさせる。各列に使用されるＭ個の正弦曲線の周波数は、予想ドップラー周波数の範囲をカバーするように選択される。固定された受信機のための一般的な選択肢は、５００Ｈｚ刻みで−５ｋＨｚから５ｋＨｚの範囲に及ぶことであろう。これにより、Ｍ＝２１が得られる。モバイル受信機は、この範囲を±１０ＫＨｚおよびＭ＝４１に拡張する必要があるであろう。次いで、フロント・エンドは、Ｗ個の時間領域サンプルのこれらのブロックの各々のゼロ・パディング離散フーリエ変換４１７を行って、Ｎ個の複素周波数領域サンプル

【数16】

４２１のベクトルを生じさせる。

【0104】

ベクトルｙ_１，ｙ_２，・・・，ｙ_Ｍ４２１が、衛星単位コード位相推定器モジュール４２０の各々に提供される。衛星ｊのためのコード位相推定器は、ｙ_ｍの各々に長さＮのベクトル

【数17】

を要素ごとに乗算し４２２、これは、Ｃ／ＡコードｊのＮ点ゼロ・パディング離散フーリエ変換４２３である。分かりやすくするため、これは、図４では、周波数ビンｍおよび衛星ｊに対してのみ示されている。この演算の結果は、ｚ_ｍｊ４２４とラベルが付されている。

【0105】

ピーク検出ブロック４２５は、複素ベクトルｚ_ｍｊの大きさが最大の要素のロケーションに対応するコード遅延（およびドップラー周波数ω_ｊ＝２πｆ_ｍ）を出力する。これは、隣接要素の補間を使用しても改善可能である。ピーク検出ブロック４２５は、このピークの（補間された）大きさａ_ｊも出力し、これを閾値と比較して、衛星ｊからの信号が存在するかどうかを判断することが可能である。あるいは、最大ピークと２番目に大きいピークの大きさの比が、Ｔｕｒｋｍａｎ−Ｗａｌｋｅｒ検定を使用して検査可能である。

【0106】

前述の方法は、すべての衛星捕捉ユニットによって共有される、以下のフロント・エンド演算を必要とする。
・ｘ_１，ｘ_２，・・・，ｘ_Ｍを生じさせるための、２つの長さのＷ個の複素ベクトルのＭ個の要素ごとの乗算。これは、Ｏ（ＭＷ）の演算を必要とする。
・ｙ_１，ｙ_２，・・・，ｙ_Ｍを生じさせる、長さＮのＭ回の離散フーリエ変換。これは、Ｏ（ＭＮｌｏｇ_２Ｎ）の演算である。

【0107】

加えて、各衛星ｊ＝１，２，・・・，Ｊおよび各周波数ｍ＝１，２，・・・，Ｍに対して、ベクトルｚ_ｊｍを形成するためのＮ個の複素乗算が存在し、Ｎ個の複素数の大きさが、ピークを見つけるために比較される。Ｊ個の衛星では、これは、Ｍ（Ｗ＋Ｎｌｏｇ_２Ｎ＋２ＪＮ）の複素数演算に達する。

【0108】

たとえば、固定された受信機およびＭ＝２１、Ｗ＝Ｎ＝１０２４、およびＪ＝２４（「コールド・スタート」に相当する）では、約１．２７×１０^６回の演算がある。１００ＭＨｚで動作し、サイクルごとに単一の命令を仮定した汎用プロセッサの場合、これは、約１２．７ミリ秒のプロセッサ時間に対応する。これは、ＧＰＳアルマナックをダウンロードするのに必要とされる１２分およびアルマナック・データに基づいてＧＰＳ信号を実際に捕捉するさらなる処理時間と比較される。

【0109】

３．３．２ 共同改良

【0110】

標準的なＧＰＳ受信機は、位相ロック・ループ（ＰＬＬ）を使用して、サブ・チップ・レベルでコード位相を連続的に追跡する。これは、ＧＰＳ衛星からの信号の連続的な受信を必要とし、安定したロックを確立するのに、ある程度の時間がかかる。これは、受信機が短いウィンドウの信号サンプルへのアクセスのみを有する場合には、可能でない。

【0111】

その代わり、本明細書において説明する方法では、検出された衛星に関する粗遅延とドップラー推定値１４２が共同で改良される。一実施形態では、これは、Ｌｅｖｅｎｂｅｒｇ−Ｍａｒｑｕａｒｄｔなどの非線形最小二乗法を使用して実施可能である。これは、受信信号と検出された衛星のＣ／Ａコードの時間シフト・バージョンおよび周波数シフト・バージョンの重ね合わせとして形成された仮説信号との間の誤差ノルムを最小にすることによって機能する。これらの個々の信号は、それらの複素利得の最小二乗推定値によって重み付けされる。この方法は、加えてドップラー率を仮説信号に適用させることによって、ドップラー率の推定値も生じさせることができる。

【0112】

ｉ∈Ｉは、粗捕捉中に検出されたＧＰＳ衛星のセットにインデックスを付けるものとする。ｃ_ｉ（ｔ）は、ＧＰＳ衛星ｉからのＣ／Ａコードとする。次いで、共同改良ステップは、

【数18】

として定義される誤差信号ε（ｔ）に対して動作し、ここで、ｔ_ｉ、ω_ｉ、ｖ_ｉは、ＧＰＳ衛星ｉに対する、仮説により改良される遅延、ドップラー、およびドップラー率である。加えて、ａ_ｉは、信号ｉに対する複素利得である。

【0113】

この手法の利点は、短いウィンドウのサンプルのみを使用して動作することができ、そのため、従来のＧＰＳ受信機とは異なり、連続的な動作を必要としないことである。これは、端末のかなりの電力節約を表す。第２の重要な利点は、この手法が、遅延、ドップラー、およびドップラー率を共に推定することである。これは、ＧＰＳによって使用されるゴールド符号間の小さな（しかし、重要な）相互相関を考慮に入れる。これは、観察窓サイズも適切に説明する。

【0114】

当業者は、情報および信号が、さまざまな技術および技法のいずれかを使用して表されてよいことを理解するであろう。たとえば、上記の説明全体を通じて参照され得るデータ、命令、コマンド、情報、信号、ビット、シンボル、およびチップは、電圧、電流、電磁波、磁場もしくは磁性粒子、光場もしくは光粒子、またはそれらの任意の組み合わせによって表されてよい。

【0115】

当業者は、本明細書で開示される実施形態に関連して説明したさまざまな例示的な論理ブロック、モジュール、回路、およびアルゴリズム・ステップは、電子ハードウェア、コンピュータ・ソフトウェアもしくは命令、または両方の組み合わせとして実装されてよいことをさらに諒解するであろう。ハードウェアとソフトウェアのこの互換性を明確に示すために、さまざまな例示的な構成要素、ブロック、モジュール、回路、およびステップについて、それらの機能に関して上記で概括的に説明してきた。そのような機能をハードウェアとして実装するか、ソフトウェアとして実装するかは、特定の適用例および全体的なシステムに課された設計上の制約に依存する。当業者は、説明した機能を、各特定の適用例に対してさまざまな方策で実装してよいが、そのような実装決定は、本発明の範囲からの逸脱を引き起こすと解釈されるべきではない。

【0116】

本明細書で開示される実施形態に関連して説明した方法またはアルゴリズムのステップは、ハードウェアにおいて直接的に実施されてもよいし、プロセッサによって実行されるソフトウェア・モジュールにおいて実施されてもよいし、またはこれら２つの組み合わせで実施されてもよい。ハードウェア実装形態では、処理は、１つまたは複数の特定用途向け集積回路（ＡＳＩＣ）、デジタル信号プロセッサ（ＤＳＰ）、デジタル信号処理デバイス（ＤＳＰＤ）、プログラマブル論理デバイス（ＰＬＤ）、フィールド・プログラマブル・ゲート・アレイ（ＦＰＧＡ）、プロセッサ、コントローラ、マイクロコントローラ、マイクロプロセッサ、または本明細書において説明する機能を実行するように設計された他の電子のユニット、またはそれらの組み合わせの中で実装されてよい。いくつかの実施形態では、１つまたは複数の中央処理ユニット（ＣＰＵ）を備えるコンピューティング・デバイスは、位置の推定などの方法のステップのうちのいくつかを実行するために使用されることがある。ＣＰＵは、入出力インターフェースと、算術論理演算ユニット（ＡＬＵ）と、入出力インターフェースを通して入力デバイスおよび出力デバイスと通信する制御ユニットおよびプログラム・カウンタ要素とを備えてよい。入出力インターフェースは、あらかじめ定義された通信プロトコル（たとえば、Ｂｌｕｅｔｏｏｔｈ（登録商標）、Ｚｉｇｂｅｅ、ＩＥＥＥ８０２．１５、ＩＥＥＥ８０２．１１、ＴＣＰ／ＩＰ、ＵＤＰなど）を使用する別のデバイス内の等価な通信モジュールと通信するために、ネットワーク・インターフェースおよび／または通信モジュールを備えてよい。コンピューティング・デバイスは、単一のＣＰＵ（コア）または複数のＣＰＵ（複数のコア）を備えてもよいし、複数のプロセッサを備えてもよい。コンピューティング・デバイスは、並列プロセッサ、ベクトル・プロセッサを使用してもよいし、分散型コンピューティング・デバイスであってもよい。メモリは、プロセッサに動作可能に結合され、ＲＡＭ構成要素とＲＯＭ構成要素とを含んでよく、デバイスの中に設けられてもよいし、デバイスの外部に設けられてもよい。メモリは、オペレーティング・システムおよび追加ソフトウェア・モジュールまたは命令を記憶するために使用されてよい。プロセッサは、メモリ内に記憶されたソフトウェア・モジュールまたは命令をロードおよび実行するように構成されてよい。

【0117】

コンピュータ・プログラム、コンピュータ・コード、または命令としても知られるソフトウェア・モジュールは、いくつかのソース・コード・セグメントまたはオブジェクト・コード・セグメントまたは命令を含んでよく、ＲＡＭメモリ、フラッシュメモリ、ＲＯＭメモリ、ＥＰＲＯＭメモリ、レジスタ、バード・ディスク、リムーバブル・ディスク、ＣＤ−ＲＯＭ、ＤＶＤ−ＲＯＭ、Ｂｌｕ−ｒａｙ（登録商標）ディスク、または他の任意の形態のコンピュータ可読媒体などの、任意のコンピュータ可読媒体内に常駐してよい。いくつかの態様では、コンピュータ可読媒体は、非一時的なコンピュータ可読媒体（たとえば、有形媒体）を含んでよい。さらに、他の態様の場合、コンピュータ可読媒体は、一時的なコンピュータ可読媒体（たとえば、信号）を含んでよい。上記の組み合わせも、コンピュータ可読媒体の範囲内に含まれるべきである。別の態様では、コンピュータ可読媒体は、プロセッサと一体化されてもよい。プロセッサおよびコンピュータ可読媒体は、ＡＳＩＣまたは関連デバイス内にあってよい。ソフトウェア・コードはメモリ・ユニット内に記憶されてよく、プロセッサは、それらを実行するように構成されてよい。メモリ・ユニットは、プロセッサの中で実装されてもよいし、プロセッサの外部で実装されてもよく、その場合、メモリ・ユニットは、当技術分野で知られているようにさまざまな手段を介してプロセッサに通信可能に結合可能である。

【0118】

さらに、本明細書において説明する方法および技法を実行するためのモジュールおよび／または他の適切な手段は、コンピューティング・デバイスによってダウンロードされるおよび／または別の方法で取得されることが可能であることを諒解されたい。たとえば、そのようなデバイスは、本明細書において説明する方法を実行するための手段の移転を容易にするためにサーバに結合可能である。あるいは、本明細書において説明するさまざまな方法は、記憶手段（たとえば、ＲＡＭ、ＲＯＭ、コンパクト・ディスク（ＣＤ）またはフロッピー・ディスクなどの物理的記憶媒体など）を介して提供可能であり、したがって、コンピューティング・デバイスは、記憶手段をデバイスに結合または提供するときにさまざまな方法を取得することができる。さらに、本明細書において説明する方法および技法をデバイスに提供するための他の任意の適切な技法が利用可能である。

【0119】

本明細書で開示される方法は、説明された方法を達成するための１つまたは複数のステップまたはアクションを含む。方法ステップおよび／またはアクションは、特許請求の範囲から逸脱することなく、互いと置き換えられてよい。言い換えれば、ステップまたはアクションの特定の順序が指定されない限り、特定のステップおよび／またはアクションの順序および／または使用は、特許請求の範囲から逸脱することなく、修正されてよい。本明細書で使用されるとき、項目のリストのうちの「少なくとも１つ」を参照する句は、単一のメンバを含む、それらの項目の任意の組み合わせを指す。一例として、「ａ、ｂ、またはｃのうちの少なくとも１つ」は、ａ、ｂ、ｃ、ａ−ｂ、ａ−ｃ、ｂ−ｃ、およびａ−ｂ−ｃを包含することが意図されている。

【0120】

本明細書で使用されるとき、「推定する」または「決定する」という用語は、多種多様のアクションを包含する。たとえば、「推定する」または「決定する」は、算出する、計算する、処理する、導出する、調査する、調べる（たとえば、テーブル、データベースまたは別のデータ構造を調べる）、確かめるなどを含んでよい。また、「推定する」または「決定する」は、受信する（たとえば、情報を受信する）、アクセスする（たとえば、メモリ内のデータにアクセスする）などを含んでよい。また、「決定する」は、解く、選ぶ、選択する、確立するなどを含んでよい。

【0121】

上記の説明は、低地球軌道衛星に言及する。しかしながら、ペイロードを持っており長時間にわたってホバリングまたは飛行することができる、高高度バルーン、飛行船、または固定翼長航続時間無人航空機などの高高度プラットフォーム局を理解されたい。したがって、本明細書の文脈では、衛星は、ＬＥＯ衛星と高高度プラットフォームの両方を包含すべきものである。

【0122】

本明細書において説明するシステム、方法、および装置の実施形態は、端末自体においてではなく、地上局において、端末位置の推定値（たとえば、緯度、経度、および高度）を生じさせるために使用可能である。これは、ナビゲーション・アプリケーションと対比可能であり、ナビゲーション・アプリケーションでは、端末それ自体がそれ自体のロケーションの推定値を必要とする。

【0123】

１つの適用例は、家畜、産業用機械、出荷容器、パレット、および環境センサなどの、遠隔にあり、潜在的に可動性である、資産の追跡がある。これらは、資産またはセンサそれ自体がそれ自体の位置を知る必要はない適用例である。むしろ、中央管理システムが、すべての資産／センサのロケーションを決定することを望む。

【0124】

別の適用例は、その位置を別個に報告していてよい、資産のロケーションの独立した検証を提供している。たとえば、自動識別システムを使用する船は、自らの位置（オン・ボードＧＰＳモジュールから導出されてよい）をブロードキャストする。提案された方法は、セキュリティおよび／またはロバストネスの理由で、これらの報告された位置を独立して検証するために使用可能である。同様に、システムは、船団追跡システムにおける輸送手段の位置を検証するために使用可能である。

【0125】

提案される測位方法は、遠隔追跡アプリケーションのためのいくつかの利点を有する。
１．この方法は、端末におけるＧＰＳモジュールを必要としない。これによって、端末は組み立てるためのコストが低くなり、はるかに少ない電力を使用するようになり、バッテリ寿命の延長につながる。
２．この方法は、端末が、ＧＰＳ衛星、ビーコン衛星、または地上ビーコンの位置の知識を取得することを要求しない。たとえば、コールド・スタート・モードでは、ＧＰＳ受信機が、完全なアルマナックをダウンロードするのに１２分かかる。本明細書において説明する方法の実施形態は、このプロセスを回避する。
３．この方法は、主に、地上局において実施され、地上局では、より多くの計算リソースが利用可能であり、電力使用はそれほど重要ではなく、ＧＰＳ衛星および他の衛星のためのエフェメリス・データは、たとえばインターネット・データベースへの直接接続を介して、利用可能である。

【0126】

上記で説明したように、地上局３０（または地上局３０と通信するサイト）における位置推定器４０は、地上局３０において計算される端末の時間遅延／ドップラー／ドップラー率の推定値を、端末のロケーションに関する特定の仮説に対して予想されるであろう対応する仮説化された量と比較することによって、機能する。これらの量は、衛星の位置２０を知れば決定可能であり、その正確な推定値は、地上局に利用可能（別のソースから取得可能、または利用可能／取得可能データに基づいて推定可能、のどちらか）である。ロケーション推定値は、推定された量から仮説化された量の間の最小総合誤差をもつ仮説である。方法の実施形態は、加えて、任意選択で、端末によって受信された他のビーコン信号５０を機会主義的に活用することができる。この場合、端末１０が、ビーコン信号（または複数のビーコン信号）のための遅延／ドップラー／ドップラー率推定値のうちの１つまたは複数を計算し、または別の方法で推定し、次いで、そのようなビーコン信号５０の推定値をそのデータ送信１２の一部として低地球軌道衛星２０に送信する。本明細書において説明する方法は、非同期端末と同期端末の両方に適用可能である。本明細書において説明する方法はまた、オン・ボード全世界位置基準およびオン・ボード全世界タイミング基準の有無にかかわらず、低地球軌道衛星に独立して適用可能である。さらに、いくつかの実施形態では、遅延、ドップラー、またはドップラー率の推定は、衛星において実行されてよく、地上局におけるもしくは地上局に接続された位置推定モジュール内での位置推定、または方法全体は、衛星内で実行されてもよい（すなわち、位置推定モジュールは衛星内にある）。

【0127】

明細書および以下の特許請求の範囲の全体を通じて、文脈がそれ以外を必要としない限り、「備える（ｃｏｍｐｒｉｓｅ）」および「含む（ｉｎｃｌｕｄｅ）」という語ならびに「備える（ｃｏｍｐｒｉｓｉｎｇ）」および「含む（ｉｎｃｌｕｄｉｎｇ）」などの変形は、述べられた整数または整数のグループの包含を暗示するが、他の任意の整数または整数のグループの除外は暗示しないと理解される。

【0128】

本明細書における任意の従来技術への参照は、従来技術が共通一般知識の一部を形成する、任意の形式の提案の承認でなく、そのようにみなすべきでない。

【0129】

本開示は、その使用に関して、説明した特定の１つまたは複数の適用例に制限されないことが、当業者によって諒解されよう。本開示は、その好ましい実施形態では、本明細書において説明または図示した特定の要素および／または特徴に関して制限されない。本開示は、開示する１つまたは複数の実施形態に限定されず、以下の特許請求の範囲によって記載および定義された範囲から逸脱することなく、多数の再配列、修正、および置換が可能であることが諒解されるであろう。

【図1】

【図2】

【図3】

【図4】

【図5】

【国際調査報告】