(57)【要約】【課題】端末装置が使用するための衛星の拡張された衛星エフェメリス・データを生成する方法を提供すること。【解決手段】過去の軌道パラメータ・データ（例えば過去のＴＬＥから取られる）を使用して、軌道パラメータの将来の値の予測器を構築する。その後、軌道パラメータ予測モデルが端末装置に記憶され（または既存のモデルが更新されることができ）、所望のエポックに、予め記憶された軌道パラメータ予測器を使用してそのエポックに対応する衛星の軌道パラメータの推定値を作成する（すなわち合成または予測ＴＬＥ）。その後、予測された軌道パラメータ（予測ＴＬＥ）を標準的な軌道予測器への入力として使用して、例えば衛星通過時刻を判定するために、所望の時刻における衛星のロケーションの推定値を作成する。この方法は、数時間または数日である通常のエフェメリス・データの有効性と比べて、一年またはそれ以上にわたり有効な衛星エフェメリス・データの生成を可能にする。【選択図】図４
2020504820.tif 000002
【特許請求の範囲】
【請求項１】
衛星通信システムにおいて端末装置が使用するための衛星の拡張された衛星エフェメリス・データを生成する方法であって、該衛星は最小制御衛星であり、該方法は、
衛星の過去のエフェメリス・データを受信するステップであって、該衛星は最小制御衛星である、ステップと、
該受信した過去のエフェメリス・データを使用して、該衛星についての軌道パラメータのセットの時間的進展の推定値を生成するための動的軌道モデルを推定するステップと、
該動的軌道モデルを該端末装置によって記憶する、または該端末装置によって記憶されている動的軌道モデルを更新するステップであって、該動的軌道モデルは該端末装置の外部にある供給源から受信されるか、または該動的軌道モデルは該端末装置によって推定されるかのいずれかである、ステップと、
該端末装置により、該動的軌道モデルを使用してエポック時刻における合成軌道要素のセットを生成するステップと、
該端末装置が該合成軌道要素のセットを軌道予測器に提供して、予め定められた端末位置に対応する１つもしくは複数の衛星通過時刻の推定値、または所定の将来の時刻における１つもしくは複数の衛星ロケーションを取得するステップと、
を含む方法。
【請求項２】
前記動的軌道モデルを推定するステップが、前記端末装置から遠隔にあるコンピューティング装置によって行われ、前記動的軌道モデルが、前記端末装置に送信またはアップロードされる、請求項１記載の方法。
【請求項３】
前記軌道パラメータのセットは、軌道にある最小制御衛星の運動をモデル化するのに充分なパラメータの任意のセットであり、前記動的軌道モデルは、各軌道パラメータを直接または間接的にフィットさせるように構成される、請求項１または２記載の方法。
【請求項４】
前記軌道パラメータのセットは、ケプラー軌道パラメータ（ａ、ｅ、ｉ、Ω、ω、Ｍ）を含み、前記過去のエフェメリス・データは、時間期間ｔ＝（ｔ１，ｔ２…，ｔｎ）にわたる軌道パラメータ（ａ［ｔ］、ｅ［ｔ］、ｉ［ｔ］、Ω［ｔ］、ω［ｔ］、Ｍ［ｔ］）の値のセットを含み、動的軌道モデルを推定するステップは、軌道パラメータ（Ω、ω、Ｍ）についての前記過去のエフェメリス・データと、前記時間期間ｔにわたる前記衛星の平均運動とにモデルをフィットさせ、各軌道パラメータ（ａ、ｅ、ｉ）に対して時間期間ｔ内で最も最近の時間点（ｔｎ）に対応する前記軌道パラメータの値を使用することを含む、請求項３記載の方法。
【請求項５】
前記軌道パラメータのセットは、ケプラー軌道パラメータ（ａ、ｅ、ｉ、Ω、ω、Ｍ）を含み、前記過去のエフェメリス・データは、時間期間ｔ＝（ｔ１，ｔ２…，ｔｎ）にわたる軌道パラメータ（ａ［ｔ］、ｅ［ｔ］、ｉ［ｔ］、Ω［ｔ］、ω［ｔ］、Ｍ［ｔ］）の値のセットを含み、ここで前記ｔｎは時間期間ｔ内で最も最近の時間点の時刻であり、動的軌道モデルを推定するステップは、軌道パラメータ（ａ、ｅ、ｉ、Ω、ω、Ｍ）についての前記過去のエフェメリス・データと、前記時間期間ｔにわたる前記衛星の平均運動とにモデルをフィットさせることを含む、請求項３記載の方法。
【請求項６】
モデルをフィットさせるステップは、軌道パラメータ（ａ、ｉ、Ω、ω、Ｍ）については線形最小二乗モデルをフィットさせ、（ｅ）については非線形最小二乗モデルをフィットさせることを含む、請求項５記載の方法。
【請求項７】
前記時間点ｔｎは、最も最近の軌道パラメータ・データの時刻であり、前記時間期間は少なくとも、ｔｎに先行しかつｔｎを含んだ２週間を含む、請求項４から６のいずれか一項記載の方法。
【請求項８】
前記エポック時刻は、予め定められた対象の時間期間の直前に発生する、前記衛星の昇交点通過の時刻である、請求項１から７のいずれか一項記載の方法。
【請求項９】
時間期間ｔ＝（ｔ１，ｔ２…，ｔｎ）にわたる平均運動ρに対して線形フィットρ＝ｍρｔ＋ｃρを行って勾配パラメータｍρおよびオフセットｃρを取得し、
【数１】
2020504820.tif 000008
もしくはｍρ＝０の場合はＴ＝Ｋ／ρを使用して、最も最近の軌道パラメータ・データの時刻ｔｎからの、前記衛星が次のＫ回の完全な軌道を完了するための所要時間Ｔを求めるステップ、または、
時間の関数としての過去の軌道番号に対して線形フィットを行って完全な軌道の数Ｋの推定値を取得し、時間に対する該軌道番号の導関数から前記平均運動ρを推定し、該推定される平均運動ρを使用して、該最も最近の軌道パラメータ・データの時刻ｔｎからの、前記衛星が次のＫ回の完全な軌道を完了するための所要時間Ｔを求めるステップと、
前記エポック時刻Ｔｅｐｏｃｈ＝ｔｎ＋Ｔを求めるステップであって、ｔｎは該最も最近の軌道パラメータ・データの時刻である、ステップと、
をさらに含む、請求項８記載の方法。
【請求項１０】
モデルをフィットさせるステップが、ロバスト・フィッティング法を使用することを含む、請求項１から９のいずれか一項記載の方法。
【請求項１１】
線形モデルをフィットさせるための前記ロバスト・フィッティング法が、Ｔｈｅｉｌ−Ｓｅｎロバスト線形フィッティング法を使用することを含む、請求項１０記載の方法。
【請求項１２】
前記動的軌道モデルを使用してエポック時刻における合成軌道要素のセットを生成するステップが、該エポック時刻を使用して、外挿された軌道要素を計算することを含む、請求項１から１１のいずれか一項記載の方法。
【請求項１３】
前記予め定められた端末位置が、前記端末装置のメモリに記憶された位置である、請求項１から１２のいずれか一項記載の方法。
【請求項１４】
前記端末が衛星測位システム受信機をさらに備え、前記予め定められた端末位置が該衛星測位システム受信機から取得される、請求項１から１３のいずれか一項記載の方法。
【請求項１５】
衛星通信システムで使用するための端末装置であって、該衛星は最小制御衛星であり、該端末装置は、
少なくとも１つのアンテナと、
該少なくとも１つのアンテナを使用して該衛星にデータを送信するように少なくとも構成された通信モジュールと、
少なくとも１つのプロセッサおよびメモリと、を備え、該メモリは、請求項１から１４のいずれか一項記載の方法の端末装置ステップを該プロセッサに行わせる命令を備えている、端末装置。
【請求項１６】
前記端末装置が、前記通信モジュールを介して過去のエフェメリス・データを受信し、前記動的軌道モデルを推定するように構成されている、請求項１５記載の端末装置。
【請求項１７】
前記端末装置の位置推定値を前記軌道予測器に提供するように構成された衛星測位システム受信機をさらに備える、請求項１５記載の端末装置。
【請求項１８】
前記メモリは、前記端末装置の位置の推定値を記憶するように構成され、該位置を前記軌道予測器モジュールに提供するための命令を含む、請求項１５記載の端末装置。
【請求項１９】
衛星通過時間期間の推定値を受信し、推定される該衛星通過時間期間中に前記端末装置をウェイクアップさせるように構成されたアラーム・モジュールをさらに備える、請求項１５記載の端末装置。
【請求項２０】
請求項１から１４のいずれか一項記載の方法を実施するように構成された衛星通信システムにおけるコンピュータ装置と少なくとも１つの端末装置とを備えるシステムであって、該コンピューティング装置は、前記動的軌道モデルを推定し、前記動的軌道モデルを、請求項１の少なくとも最後の３つのステップを行う該少なくとも１つの端末装置に提供するように構成されている、システム。
【発明の詳細な説明】
【技術分野】
【０００１】
＜優先権文献＞
本願は、２０１６年１２月２２日に出願されたオーストラリア仮特許出願第２０１６９０５３１４号、名称「ＳＹＳＴＥＭ  ＡＮＤ  ＭＥＴＨＯＤ  ＦＯＲ  ＧＥＮＥＲＡＴＩＮＧ  ＥＸＴＥＮＤＥＤ  ＳＡＴＥＬＬＩＴＥ  ＥＰＨＥＭＥＲＩＳ  ＤＡＴＡ」に対する優先権を主張し、同出願の内容は参照により全体が本明細書に援用される。
【０００２】
本開示は、衛星軌道に関する。特定の形態では、本開示は、衛星エフェメリス（ｅｐｈｅｍｅｒｉｓ）データを生成する方法に関する。
【背景技術】
【０００３】
地球を周回する軌道にある人工衛星は、遠隔通信および測位の用途に使用されることがある。そのようなシステムでは、地球の表面上またはその近傍に位置するデバイスが、衛星からの信号を受信するか、または衛星に信号を送信する。これらの用途の多くでは、軌道にある衛星の精密な位置を知ることが必要である。例えば、全地球測位システム（ＧＰＳ）は、相対論的効果を考慮する必要があるほど精密にシステムの衛星の位置を判定できることに依拠する。いくつかの遠隔通信用途では、高度に指向性の（場合によっては操縦可能な）アンテナを衛星の方に向ける必要がある。
【０００４】
他の用途では、衛星の精密な位置は必要とされず、むしろ各衛星通過のタイミングを予測する必要がある。衛星通過とは、地上の観察者にとって衛星が有用である時間期間を言う。例えば、これは、衛星が観察者に対して任意の最小の仰角を超えている時間間隔であり得る。
【０００５】
多くの衛星通信および測位システムに対する主要な要件の１つは、地上のデバイスが衛星通過を予測できることである。ＧＰＳシステムの場合、各衛星の通過時刻を知ることにより、ＧＰＳ受信機は任意の特定の時に見るべき衛星を知ることができ、それによりＧＰＳ受信機は衛星を迅速に取得し、その衛星にロックすることができ、その結果、最初の位置算出までの時間が短縮する。すなわち、ＧＰＳ受信機は、この情報を使用して、ユーザ（またはユーザ・デバイス）に対して位置推定値を迅速に作成することができる。多くの補助ＧＰＳシステムは、ユーザ体感を改善する（すなわち待ち時間を短縮することによる）ために、最初の位置算出までの時間を可能な限り短く短縮する（例えば秒以下）ことに着目する。
【０００６】
通信システムでは、通過時刻を知ることにより、トランシーバ・デバイスは、どの通信と通信すべきか、およびどの時間に通信すべきかを知ることができる。そのようなシステムでは、最初の位置算出までの時間はそれほど重要ではないが、一般には電力を浪費するため、一般には最初の位置算出までの時間は数十分から数時間ではないことが望ましい。
【０００７】
通過時刻を推定するための既存の方法は軌道予測アルゴリズムに基づいており、これらのアルゴリズムは、通例は２行要素（ＴＬＥ（ｔｗｏ−ｌｉｎｅ  ｅｌｅｍｅｎｔ））であるエフェメリス・データ（または軌道要素）を使用して衛星の軌道をモデル化する。ＳＧＰ４アルゴリズム（Ｆｅｌｉｘ  Ｒ．Ｈｏｏｔｓ  ａｎｄ  Ｒｏｎａｌｄ  Ｌ．Ｒｏｅｈｒｉｃｈ，  “Ｍｏｄｅｌｓ  ｆｏｒ  ｐｒｏｐａｇａｔｉｏｎ  ｏｆ  ＮＯＲＡＤ  ｅｌｅｍｅｎｔ  ｓｅｔｓ”，  Ｔｅｃｈｎｉｃａｌ  Ｒｅｐｏｒｔ  ３，  Ｓｐａｃｅｔｒａｃｋ，  １９８８）は、ＴＬＥデータに使用するために設計された軌道予測アルゴリズムの一例である。
【０００８】
これらの方法は、衛星の位置および通過タイミングの非常に正確な予測を提供することができる。しかし、これらの方法はＴＬＥデータが古過ぎない場合にのみ有効であることは十分に立証されている。ＴＬＥデータに使用される標準的なＳＧＰ４予測器は、ＴＬＥデータが約一カ月以上経過すると、非常に不正確になる。この不正確さは、（ａ）軌道ドリフト、および（ｂ）予測器自体の不正確さ、という２つの主要な影響に起因し得る。
【０００９】
衛星の軌道パラメータは、時間と共に変化することがある（自然力に起因するか、または衛星の外部制御に起因するかのいずれか）。一部の衛星は厳密に制御された軌道にあり、「軌道上」にあり続けるために推進剤を費やす。例えば、中地球軌道（ＭＥＯ）にあるＧＰＳ衛星の軌道は、衛星を望ましい軌道に保つために、常に監視され、頻繁に調節される（例えば毎週から数カ月ごとの時間スケール）。
【００１０】
他の衛星は、無制御であるかまたは弱い制御のみを受けている場合がある。ここでは、そのような衛星をまとめて最小制御衛星と呼ぶことにする。これは、多くの場合、多量の推進剤を積んでおらず、ステーション維持の目的というよりは、衝突回避のために操縦動作を確保しておく、低地球軌道にある小型の衛星が該当する。
【００１１】
太陽同期軌道は、衛星が常に同じローカル太陽時に赤道上空（または実際には地球表面上の任意の地点）を通過するように設計される。これらの軌道は、地球の非球面性によって軌道が正確に一年に一回歳差運動するように、傾き（通例は近極）と高度との精密に選ばれた組み合わせを用いて設計される。実際には、これは厳密にはならず、無制御の軌道は意図される面からドリフトすることになる。
【００１２】
図１は、名目上は太陽同期軌道にある小型衛星についての昇交点の平均ローカル時刻（ＭＬＴＡＮ）（過去のＴＬＥデータから取得された）の曲線１をプロットしている。横軸は打ち上げからの日数を示している。縦軸は、ＵＴＣ時間で表された、昇交点の平均ローカル時刻を示している。軌道は太陽同期となるように設計されているが、その軌道面は、一年につきおよそ１時間ずつドリフトする。これらのようなシナリオに期限切れのＴＬＥデータを使用すると、結果として不正確な予測となる。図１の例では、期限切れになって１年が経過したＴＬＥでは、少なくとも１時間の誤差が生じることになる。そのため、衛星の獲得または衛星との通信を試みる地上端末は、衛星が視界になかった間に衛星の獲得を試みて１時間分の電力を浪費するか、または、衛星が当該メッセージを受信することを誤って予想したメッセージを送信する可能性があり、その場合、メッセージは失われる可能性があり、またはメッセージの再送信を試みて電力が浪費される可能性がある。しばしば、電力は衛星端末にとって貴重な資源である。
【００１３】
一方、図２は、制御された太陽同期軌道にある衛星についての昇交点の平均ローカル時刻（ＭＬＴＡＮ）の曲線をプロットしている。軌道面の変動は、約４年後に、図１における約４時間と比べてわずか数分である（両図のスケールの違いに注意されたい）。
【００１４】
軌道ドリフトを見越して、ＳＧＰ４予測アルゴリズムが「新鮮でない」ＴＬＥデータに使用されると、安定軌道の場合でも、アルゴリズム自体が不正確な結果を生ずることになる。例えば、図２に示す制御軌道衛星は、その軌道面を一カ月につき約４０秒変動させる。しかし、期限切れになって６カ月が経過したＴＬＥをＳＧＰ４に使用すると、通過タイミングの不正確さは１時間を越える結果となる。
【００１５】
衛星システムは現在、衛星の位置およびダイナミクスを定期的に観測し、更新されたＴＬＥデータを作成し（一日につき数回）、更新されたデータを、最新のＴＬＥを必要とする地上のデバイスに定期的に送信することにより、この問題を解決している。例えば、ＧＰＳシステムでは、衛星は４時間の間有効である新しいエフェメリスを２時間おきにブロードキャストし、すべての衛星が、そのフレーム構造の一部としてすべてのエフェメリス・データを含んでいる完全なアルマナック（ａｌｍａｎａｃ）を継続的に送信する（完全なアルマナックは衛星からダウンロードするのに１２．５分かかる）。
【発明の概要】
【発明が解決しようとする課題】
【００１６】
しかし、更新されたエフェメリス・データを送信するには、更新すべきデバイスへの通信リンク（または複数のリンク）の存在が前提となる。例えば、ＧＰＳシステムでは、エフェメリスはＧＰＳ衛星から取得されることがあり、または、補助ＧＰＳシステムでは、これは地上の送信機によって、または別のデータ・リンクを通じて提供されることがある。しかし、多くの衛星通信システムでは、デバイス（端末）への通信は、高費用であるか、利用不可能であるか、またはその他の理由で望ましくないことがある。そのため、いくつかの用途では、衛星通過タイミングの適度に正確な推定値を見つけることができる期間を拡張することに関心が持たれ得る。それらの用途には以下が含まれる。
１．）エフェメリス・データを定期的に更新することが難しいかまたは不可能であり得る、低帯域幅または断続的な通信リンクを備えるシステム、
２．）例えば一方向通信デバイスなど、エフェメリス・データを遠隔から更新する能力を持たないシステム、
３．）データ通信コスト、バッテリ使用、または他のコストのいずれかの点から、定期的にエフェメリス・データを更新するコストが高すぎるシステム、および
４．）デバイスに送信されるエフェメリス・データの通信妨害または電波妨害によって影響され得るシステム。
【００１７】
そのため、エフェメリス・データの更新が必要となるまでに端末デバイスが何カ月またはそれ以上にわたってエフェメリス・データを記憶もしくは生成できるように、拡張された衛星エフェメリス・データを生成するシステムおよび方法を提供するか、または少なくとも現在のシステムおよび方法の代替物を提供する必要性がある。
【課題を解決するための手段】
【００１８】
本明細書に記載される実施形態は、現在のＴＬＥおよびＳＧＰ４などの軌道予測器を使用して、現在可能なよりもはるかに長い時間枠にわたって衛星通過時刻を正確に予測するために使用することができる。
方法は、衛星の過去のエフェメリス・データが提供される衛星端末内で実施することができ、または、方法は、コンピューティング装置が動的軌道モデルの過去のエフェメリス・データを推定し、この動的軌道モデルが衛星端末に提供される通信システム内で分散させることができる。端末は、動的軌道モデルを使用して所望の時刻（すなわち、衛星への送信または衛星からの受信のためにウェイクアップ時刻をスケジュールしたい時）における合成軌道要素を生成し、その後これらの合成軌道要素（例えば合成ＴＬＥ）はＳＧＰ４のような標準的な軌道予測器などの軌道予測器に渡される。
【００１９】
第１の態様によれば、端末装置が使用するための衛星の拡張された衛星エフェメリス・データを生成する方法が提供され、該衛星は最小制御衛星であり、方法は、
衛星の過去のエフェメリス・データを受信するステップであって、衛星は最小制御衛星である、ステップと、
受信した過去のエフェメリス・データを使用して、衛星についての軌道パラメータのセットの時間的進展の推定値を生成するための動的軌道モデルを推定するステップと、
動的軌道モデルを端末装置に提供するステップであって、端末は動的軌道モデルを記憶する、または記憶された動的軌道モデルを更新する、ステップと、
端末装置により、動的軌道モデルを使用してエポック時刻における合成軌道要素のセットを生成するステップと、
端末装置が合成軌道要素のセットを軌道予測器に提供して、予め定められた端末位置に対応する１つもしくは複数の衛星通過時刻の推定値、または所定の将来の時刻における１つもしくは複数の衛星ロケーションを取得するステップと、を含む。
【００２０】
一形態では、動的軌道モデルを推定するステップは、端末装置から遠隔にあるコンピューティング装置によって行われ、動的軌道モデルは、端末装置に送信またはアップロードされる。
【００２１】
軌道パラメータのセットは、軌道にある最小制御衛星の運動をモデル化するのに充分なパラメータの任意のセットであり、動的軌道モデルは、各軌道パラメータを直接または間接的にフィットさせるように構成される。一形態では、軌道パラメータのセットは、ケプラー軌道パラメータ（ａ、ｅ、ｉ、Ω、ω、Ｍ）を含み、過去のエフェメリス・データは、時間期間ｔ＝（ｔ１，ｔ２…，ｔｎ）にわたる軌道パラメータ（ａ［ｔ］、ｅ［ｔ］、ｉ［ｔ］、Ω［ｔ］、ω［ｔ］、Ｍ［ｔ］）の値のセットを含み、動的軌道モデルを推定するステップは、軌道パラメータ（Ω、ω、Ｍ）についての過去のエフェメリス・データと、時間期間ｔにわたる衛星の平均運動とにモデルをフィットさせ、各軌道パラメータ（ａ、ｅ、ｉ）に対して時間期間ｔ内で最も最近の時間点（ｔｎ）に対応する軌道パラメータの値を使用することを含む。別の形態では、動的軌道モデルを推定するステップは、軌道パラメータ（ａ、ｅ、ｉ、Ω、ω、Ｍ）についての過去のエフェメリス・データと、時間期間ｔにわたる衛星の平均運動ρとにモデルをフィットさせることを含み、ここでｔｎは時間期間ｔ内で最も最近の時間点の時刻である。さらなる形態では、モデルをフィットさせるステップは、軌道パラメータ（ａ、ｉ、Ω、ω、Ｍ）については線形最小二乗モデルをフィットさせ、（ｅ）については非線形最小二乗モデルをフィットさせることを含む。一形態では、時間点ｔｎは、最も最近の軌道パラメータ・データの時刻であり、時間期間は少なくとも、ｔｎに先行しかつｔｎを含んだ２週間を含む。
【００２２】
一形態では、エポック時刻は、予め定められた対象の時間期間の直前に発生する、衛星の昇交点通過の時刻である。さらなる形態では、上記方法は、
時間期間ｔ＝（ｔ１，ｔ２…，ｔｎ）にわたる平均運動ρに対して線形フィットρ＝ｍρｔ＋ｃρを行って勾配パラメータｍρおよびオフセットｃρを取得し、
【数１】
2020504820.tif 000003
もしくはｍρ＝０の場合はＴ＝Ｋ／ρを使用して、最も最近の軌道パラメータ・データの時刻ｔｎからの、衛星が次のＫ回の完全な軌道を完了するための所要時間Ｔを求めるステップ、または、
時間の関数としての過去の軌道番号に対して線形フィットを行って完全な軌道の数Ｋの推定値を取得し、時間に対する軌道番号の導関数から平均運動ρを推定し、推定される平均運動ρを使用して、最も最近の軌道パラメータ・データの時刻ｔｎからの、衛星が次のＫ回の完全な軌道を完了するための所要時間Ｔを求めるステップと、をさらに含む。
所要時間Ｔが取得されると、Ｔｅｐｏｃｈ＝ｔｎ＋Ｔを使用してエポック時刻が求められる。
【００２３】
一形態では、モデルをフィットさせるステップが、ロバスト・フィッティング法を使用することを含む。さらなる形態では、線形モデルをフィットさせるためのロバスト・フィッティング法は、Ｔｈｅｉｌ−Ｓｅｎロバスト線形フィッティング法を使用することを含む。一形態では、動的軌道モデルを使用してエポック時刻における合成軌道要素のセットを生成するステップは、エポック時刻を使用して、外挿された軌道要素を計算することを含む。一形態では、予め定められた端末位置は、端末装置のメモリに記憶された位置である。一形態では、端末は衛星測位システム受信機をさらに備え、予め定められた端末位置は衛星測位システム受信機から取得される。
【００２４】
第２の態様によれば、衛星通信システムで使用するための端末装置が提供され、衛星は最小制御衛星であり、端末装置は、
少なくとも１つのアンテナと、
少なくとも１つのアンテナを使用して衛星にデータを送信するように少なくとも構成された通信モジュールと、
少なくとも１つのプロセッサおよびメモリと、を備え、メモリは、第１の態様の端末装置ステップをプロセッサに行わせる命令を備えている。
【００２５】
一形態では、端末装置は、通信モジュールを介して過去のエフェメリス・データを受信し、動的軌道モデルを推定するように構成されている。一形態では、端末装置は、端末装置の位置推定値を軌道予測器に提供するように構成された衛星測位システム受信機をさらに備える。一形態では、メモリは、端末装置の位置の推定値を記憶するように構成され、位置を軌道予測器モジュールに提供するための命令を含む。一形態では、端末装置は、衛星通過時間期間の推定値を受信し、推定される衛星通過時間期間中に端末装置をウェイクアップさせるように構成されたアラーム・モジュールをさらに備える。
【００２６】
第３の態様によれば、第１の態様の方法を実施するように構成された衛星通信システムにおけるコンピューティング装置と端末装置とを備えるシステムが提供され、該コンピューティング装置は、動的軌道モデルを推定し、動的軌道モデルを、第１の態様の方法の少なくとも最後の３つのステップを行う少なくとも１つの端末装置に提供するように構成されている。
【００２７】
本開示の実施形態を、以下で添付図面を参照して解説する。
【図面の簡単な説明】
【００２８】
【図１】名目上は太陽同期軌道にある小型衛星についての昇交点の平均ローカル時刻（ＭＬＴＡＮ）の曲線のプロットである。
【００２９】
【図２】制御された太陽同期軌道にある衛星についての昇交点の平均ローカル時刻（ＭＬＴＡＮ）の曲線のプロットである。
【００３０】
【図３】楕円軌道およびその楕円軌道を定義する６つのケプラー軌道要素の概略プロットである。
【００３１】
【図４】一実施形態に係る端末装置が使用するための衛星の拡張された衛星エフェメリス・データを生成する方法のフローチャートである。
【００３２】
【図５Ａ】一実施形態に係る衛星軌道の準主軸ａの過去の値のプロットである。
【００３３】
【図５Ｂ】一実施形態に係る衛星軌道の離心率ｅの過去の値のプロットである。
【００３４】
【図５Ｃ】一実施形態に係る衛星軌道の傾きｉの過去の値のプロットである。
【００３５】
【図５Ｄ】一実施形態に係る衛星軌道の赤経昇交点Ωの過去の値のプロットである。
【００３６】
【図５Ｅ】一実施形態に係る衛星軌道の近地点の引数ωの過去の値のプロットである。
【００３７】
【図５Ｆ】一実施形態に係る衛星軌道の平均近点離角Ｍの過去の値のプロットである。
【００３８】
【図６】一実施形態に係る衛星の平均運動のプロットである。
【００３９】
【図７】一実施形態に係るエポックの予測を実際のエポックと比較するプロットである。
【００４０】
【図８Ａ】一実施形態に係る衛星軌道の推定される準主軸ａの時間の誤差のプロットである。
【００４１】
【図８Ｂ】一実施形態に係る衛星軌道の推定される離心率ｅの時間の誤差のプロットである。
【００４２】
【図８Ｃ】一実施形態に係る衛星軌道の推定される傾きｉの時間の誤差のプロットである。
【００４３】
【図８Ｄ】一実施形態に係る衛星軌道の推定される赤経昇交点Ωの時間の誤差のプロットである。
【００４４】
【図８Ｅ】一実施形態に係る衛星軌道の推定される近地点の引数ωの時間の誤差のプロットである。
【００４５】
【図８Ｆ】一実施形態に係る衛星軌道の推定される平均近点離角Ｍの時間の誤差のプロットである。
【００４６】
【図９】一実施形態に係る端末装置の概略図である。
【発明を実施するための形態】
【００４７】
以下の説明では、同様の参照符号は、すべての図を通じて同様のまたは対応する部分を指し示す。
【００４８】
次いで、端末装置が使用するための衛星の拡張された衛星エフェメリス・データを生成する方法の実施形態について説明する。この新しい方法は、現在のＴＬＥおよびＳＧＰ４などの軌道予測器を使用して、現在可能なよりもはるかに長い時間枠にわたって衛星通過時刻を正確に予測するために使用することができる。この方法の実施形態は、無制御軌道または最小制御軌道にある衛星に特に適用可能である。大まかには、過去の軌道パラメータ・データ（例えば過去のＴＬＥから取られる）を使用して、軌道パラメータの将来の値の予測器を構築する。その後、軌道パラメータ予測モデルが端末装置に記憶され（または既存のモデルが更新されることができ）、所望のエポックに、予め記憶された軌道パラメータ予測器を使用してそのエポックに対応する衛星の軌道パラメータの推定値を作成する（すなわち合成または予測ＴＬＥ）。その後、予測された軌道パラメータを軌道予測器への入力として使用して、例えば衛星通過時刻または衛星がいつ端末に対して特定の位置に来るかを判定するために、衛星のロケーションの推定値を作成する。
【００４９】
表１にリストされ、図３に例示される６つのケプラー軌道要素によって楕円軌道３を記述することができる。この表は、以降で使用される表記および単位も紹介している。単位の選択は任意であり、例えば、角度量は、代わりにラジアンおよびメートル単位の距離で表されることも可能である。
【表１】
2020504820.tif 000004
【００５０】
ケプラー・モデルは理想化されており、地球軌道にある物体に関するいくつかの追加的な影響、例えば、大気抵抗、太陽および月への重力による引力を考慮に入れない。しかし、この実施形態では、この理想的ケプラー・モデルを使用すれば十分である。他の実施形態では、これらの追加的な影響を考慮に入れるモデルを使用することもできる。
【００５１】
数量ω、ｅ、Ω、ｉ、ａは、軌道楕円の幾何学的形状および地球に対するその向きを記述し、一方、Ｍは、軌道における衛星の位置を記述する（すなわち楕円上の点）。離心率ｅおよび準主軸ａは、楕円に関係する通常の数量である。さらに、ａは、主天体からの衛星の平均距離である。傾きｉは、地球の回転軸に対する軌道面の角度傾斜である。赤道軌道はｉ＝０度であり、極軌道はｉ＝９０度である。昇交点は、軌道を回っている物体が赤道面を南から北に交差する、軌道内の点である。赤経昇交点Ωは、春分点から測定した昇交点の角距離（経度）である。近地点の引数ωは、軌道の周りで昇交点から地球に最も近い軌道内の点まで測定される角距離である。
【００５２】
Ｔを軌道周期とし、ｎ＝３６０／Ｔを時間単位当たりの度単位の平均角運動とする。衛星が時刻ｔ０に近地点にあり、現在の時間がｔ０＋Δｔであるとする。すると、平均近点離角は以下となる。
Ｍ＝ｎ・Δｔ      式（１）
【００５３】
これは、周期Ｔの円軌道にある仮想上の衛星が同じ量の時間内に進むであろう、近地点から測定された角度である。あるいは、離心近点角Ｅを使用することができ、これら２つの数量が、
Ｍ＝Ｅ−ｅｓｉｎＥ      式（２）
を介して関係付けられ、円軌道の場合、ｅ＝０およびＭ＝Ｅである。
【００５４】
図４は、端末装置１００が使用するための衛星４の拡張された衛星エフェメリス・データを生成する方法のフローチャート１０である。この方法は、衛星４が最小制御される衛星通信システムで使用するのに特に適する。上記で定義したように、これは、例えば、６カ月に一度、一年に一度、数年に一度、緊急状況時のみに補正されるか、または全く補正が行われないなど、低頻度のみで、小さい量（すなわち小さいΔｖ）だけ補正される衛星である。多くのマイクロサット、ナノサット、およびキューブサットを含む多くの低地球軌道衛星はこの基準を満たす。図４を参照すると、どこで演算が行われるかによって定義される２つの主要な演算ブロック１１、２１がある。第１の主要コンポーネント１１は、衛星４の軌道パラメータ（例えばケプラー軌道パラメータ）の時間的進展の動的軌道モデル２０を推定するためのソフトウェア・モジュールである。下記で説明するように、推定またはフィッティング工程１４は、過去のエフェメリス・データ１２にモデルをフィットさせ、動的軌道モデル２０をまとめて定義するモデル・パラメータのセット１８を取得する。
【００５５】
動的軌道モデル１１のこの推定は、端末装置１００とは異なるコンピューティング装置（またはシステム）で行われる、オフラインまたはリモートの演算１６として行われてよい。例えば、これは、過去のエフェメリス・データ１２にアクセスできる、独立型コンピュータ、スマートフォン、クラウド・サーバ等であってよい。いくつかの実施形態では、このコンピューティング装置は、衛星通信システムの中央ハブ内に、または衛星通信システムの中央ハブとして働くクラウド・ベースのサーバ上に位置する。これは、継続的に、毎月など周期的に、または６カ月もしくは１２カ月ごとのようにより低頻度で、または要求に応じて行うことができる。生成されると、これは周期的に端末装置１００に提供される。これは、データ接続が利用可能な時は有線もしくは無線のリンクを通じて、もしくはリモートの更新手順の一部として提供されてよく、または例えばサービス時などに手動でアップロードされてもよい（例えばＵＳＢなどの物理的データ接続により）。
【００５６】
あるいは、過去のエフェメリス・データ１２が端末装置に供給されるかまたは端末装置によって取得される場合、端末装置は、この推定工程を時折デバイス上で行い、これを使用して記憶された動的軌道モデルを更新することができる。例えば、端末装置は、時折遠隔から更新され、サービスされ（例えば年に一度のサービス）てよく、または過去のエフェメリス・データを取得できるように時折データ接続へのアクセスを有してもよい。
【００５７】
動的軌道モデルが端末装置の外部にある供給源から生成された場合、動的軌道モデルは端末装置に送られるか、または提供される。動的軌道モデルが端末装置に対して利用可能になると、動的軌道モデルは端末装置によって記憶されるか、または、新たに計算されたモデル・パラメータに基づいて、以前に記憶された動的モデルが更新（もしくは置換）される。
【００５８】
第２の主要な演算ブロック２１は、動的軌道モデルが生成され１１、端末装置のメモリに記憶された後に、端末装置１００上で行われる。このブロックで、端末装置は、記憶された動的軌道モデル（またはより具体的にはモデル・パラメータのセット１８）を使用して、衛星通過時刻の予測値４０、および／または衛星位置４２、および／またはデバイスに対する衛星位置４４を取得する。この工程は、衛星のセット（またはコンステレーション）ごとに行うことができる。１つの実施形態では、端末装置はアラーム・モジュール４６を備え、アラーム・モジュール４６は、低電力のスリープ・モードから端末装置をウェイクさせて、衛星が端末装置の視界内にあると推定される時に衛星の１つにデータを送信するかまたは衛星の１つからデータを受信するために、この情報を使用して対象の１つまたは複数の衛星の通過時刻４０を判定する。
【００５９】
方法は大まかに、動的軌道モデル２０を使用してエポック時刻２８における合成軌道要素のセット３２を生成することを含む。このことは、内部クロック、外部クロック、またはＵＴＣに同期されたタイミング源などの時刻源２２から参照時刻２４を取得することを伴う。エポック演算２６を行って、合成軌道要素のセット３２を生成することが望まれる適切なエポック時刻を求める。通例、これは、衛星のロケーションが予測されることになる所望の時間よりも前の短い量の時間である。１つの実施形態では、エポック時刻は、参照時刻に対する昇交点の最も最近の時刻など、軌道モデルに依存するか、または軌道モデルに基づいて選択される。エポック時刻演算２６が行われ、エポック時刻２８が取得されると、軌道要素予測器モジュール３０が、動的軌道モデル２０についてのモデル・パラメータのセット１８を使用して、そのエポックに対応する軌道要素のセット３２を生成する。これらを、ローカルに生成された推定値であることを示すために合成軌道要素と呼ぶことにする（衛星または参照源によって送信された軌道要素に対して）。合成軌道要素３２は、ケプラー・パラメータのセット、または軌道予測器３４が使用するためのＴＬＥなど、予め定められた軌道要素のセットを備える。その後、合成軌道要素３２は、ＳＧＰ４、Ｊ２、または何らかの他の代替物などの軌道予測器３２に提供される。軌道予測器２４は、衛星４２の絶対ロケーション、および／または端末位置３８が既知である場合には端末装置に対する衛星位置４４、および／または端末の位置が既知である場合には１つまたは複数の通過時刻４０を予測するために使用される。端末装置３８デバイスのロケーションは位置源３６から取得され、位置源３６は、静止デバイスの場合には端末装置の既知の位置を記憶するメモリであってよく、または位置３８はＧＰＳなどの測位サービス３６から取得されてもよい。この方法の実施形態について以下でより詳細に解説する。
【００６０】
上記で概説したように、動的軌道モデル２０の推定は、軌道予測器が使用するための改善された（また長期間有効な）軌道要素を予測するために当該パラメータを使用することを可能にする軌道モデルを定義する軌道パラメータ１８の時間的進展をモデル化することに基づく。以下の実施形態では、軌道モデルはケプラー・パラメータとして表されるが、この選択は任意であり、便宜上採用されるに過ぎない。この軌道モデルは、通例は低地球軌道に見られる、無制御の太陽同期衛星などの無制御衛星、ならびにやはり通例は低地球軌道に見られる最小制御衛星に適する。モデルによって使用されるモデル軌道パラメータのセットは、軌道にある最小制御衛星の運動をモデル化するのに十分なパラメータの任意のセットである。他の実施形態では、より精緻化された軌道モデルと、座標系をモデル化するためのパラメータの他の選択とを使用することができる。同様に、記載される方法は、他の種類の軌道に対する用途を有し、それらの軌道への用途は、本明細書に記載される技法の単純な発展である。
【００６１】
時間期間Ｔにわたる過去の軌道パラメータ・データ１２を使用して、それらの軌道パラメータが時間と共にどのように進展したかと整合するモデルを開発またはフィットし、それを使用してこれらの軌道パラメータ１８の将来の値を予測することができる。図５Ａ〜５Ｆは、先の図１の例と同じ低地球軌道衛星についてのケプラー要素の過去の値を示す。横軸は打ち上げからの日数である。Ω、ω、およびＭに対応する角度は「アンラップ（ｕｎｗｒａｐｐｅｄ）」である。図５Ａは、準主軸ａの過去の値の曲線５１を示し、図５Ｂは、離心率ｅの過去の値の曲線５２を示し、図５Ｃは、傾きｉの過去の値の曲線５３を示し、図５Ｄは、赤経昇交点Ωの過去の値の曲線５４を示し、図５Ｅは、近地点の引数ωの過去の値の曲線５５を示し、図５Ｆは、平均近点離角Ｍの過去の値の曲線５６を示す。
【００６２】
これらの図から、Ｍは時間と共に線形に進展することが見て取れる（またこれは式（１）から予想される）。興味深いことに、ωおよびΩは両方とも時間と共に線形に進展するようにも見える。ａ、ｉ、およびｅはそれほど線形ではないが、それらは大きく変動してもいない。そのため、１つの実施形態では、軌道パラメータについてのモデルは以下のように生成される。ここでは、ａ、ｅ、ｉ、Ω、ω、Ｍの各値の最も最近の組がいくらかの時間期間（または時間の間隔）Ｔにわたって入手可能であるとする。ｔを時間変数とし、時間期間ｔにわたる時刻インスタンスｔ＝（ｔ１，ｔ２．．．，ｔｎ）についてのデータを持っているものとする。時刻インスタンスｔ１は、時間期間ｔ内で最も古い時間点であり、ｔｎは時間期間ｔの間で最も最近の時間点である。時刻インスタンスは、軌道パラメータの発行の時刻、または特定の軌道位置（例えば軌道の開始や、軌道内での参照点）などの任意の事象のエポック時刻であってよい。対象の時間期間ｔ内で最も最近の時間点ｔｎは、最も最近の軌道パラメータ・データの時刻（すなわちそれ以降軌道パラメータ・データが発行されていない）であってよく、または多少古いものであってもよいが、好ましくは過去１カ月以内または最も最近の軌道補正以降のものである。時間期間は、直前の２週間（ｔｎに先行し、かつｔｎを含む）、直前の１カ月にわたるか、または、軌道パラメータ・データの最も最近の発行（ｔｎ）から最も最近の軌道補正までの時間期間であってよい。いくつかの実施形態では、時間期間は、直近のｎ個のデータ点によって定義されてよい（すなわち、直近の１０、２０、または３０個の値のみを使用することができる）。これらの時刻インスタンスは、時間期間にわたって規則的な間隔である場合も、そうでない場合もある（またＴＬＥデータの場合は通例そうではない）ことに留意されたい。さらに、時刻インスタンスは、最初の時刻インスタンスと最後の時刻インスタンスとの間でデータが発行されたあらゆる時刻インスタンスである必要はない（すなわち時間期間は連続した時間点を備える必要がある）。例えば、１つおきのデータ点を使用することができ、または１日につき単一のデータ点がいくらかの時間期間にわたって使用されてもよい等である。ａ［ｔ］、ｅ［ｔ］、ｉ［ｔ］、Ω［ｔ］、ω［ｔ］、Ｍ［ｔ］のそれぞれを、所与の時刻ｔにおける軌道パラメータの各値とする。
【００６３】
そして、例えば、線形もしくは非線形の回帰もしくはフィッティング技術を使用するか、またはカルマン・フィルタ（および関連するカルマン予測器）もしくはベイズ推定器などの代替の推定技術を使用することにより、各時系列データ・セットに対してモデルをフィットさせる。本明細書の文脈では、モデルのフィッティングは、回帰技術およびそれに関係する推定技術の両方を伴う。
【００６４】
フィッティング技法は、対象の特定のパラメータに合わせて選択されるべきである。例えば、線形最小二乗を使用して、およそ線形に進展するパラメータについてのフィッティング・パラメータを取得することができる。あるいは、例えば二次または正弦など、より高次の多項式フィットまたは非線形の周期関数が適することもある。ノイズのあるデータに対して線形関数または二次関数などの多項式をフィットさせることは標準的な問題であり、応用することが可能なよく知られた解を有する。データによっては、非線形最小二乗などの非線形フィッティングを使用することが有利であり得る。ここでも、すでに利用可能となっているそのようなアルゴリズムの多数の選択肢がある。一部の場合では、適切なフィッティング関数の選択は、反復的な工程であり得る。例えば、線形フィットが剰余における構造または湾曲を示す場合、このことは非線形の二次または多項式フィットがより適することを示す可能性があり、その後、そのようなアプローチを使用してデータが再フィットされてよい。あるいは、高次多項式（例えば３次式）を使用してフィッティングが行われてよく、その後重要でない項が除去され、データが再フィッティングされてよい。フィッティング技法は、各軌道パラメータを直接フィットさせても、または間接的にフィットさせてもよい。例えば、モデル・フィッティングは、変数の変化、変数の関数を利用する、および／または座標の変換を利用するなどにより、ケプラー・パラメータに基づくモデルをフィットさせてよい。例えば、１つの実施形態では、平均近点離角Ｍを直接フィットさせるのではなく、平均近点離角の代わりに衛星の真の近点離角を使用／フィットし、真の近点離角から平均近点離角が導出される。同様に、平均近点離角の代わりに衛星の離心近点角を使用し／フィットさせ、ケプラーの方程式により離心近点角から平均近点離角が導出される。別の実施形態では、平均運動の代わりに衛星の軌道番号を使用し／フィットさせ、時間に対する軌道番号の導関数として平均運動が導出される。
【００６５】
また、ロバスト・フィッティング（またはロバスト回帰）法を使用して外れ値の影響を低めに考慮してもよい。取得された過去のデータの一部は、（ａ）ノイズがある／信頼性が低いか、または（ｂ）例えば小さな操縦に起因する外れ値を含むことがある。そのような場合には、外れ値を予め考慮するかまたは低めに考慮するロバスト・フィッティング（または回帰）アルゴリズム、例えば、反復的に重み付けが更新される最小二乗、ｌｅａｓｔ  ｔｒｉｍｍｅｄ  ｓｑｕａｒｅｓ、またはＴｈｅｉｌ−Ｓｅｎ法のようなロバスト技術を使用することが有益であり得る。Ｔｈｅｉｌ−Ｓｅｎ法は、サンプル点のすべてのペアを通るすべての線の勾配の中央値（厳密には、それぞれ異なるｘ座標をもつ点のペアからのみ定義された勾配の中央値）を選ぶことにより、サンプル点を通るように線をフィットさせる。ロバスト・フィッティング法の使用は、データの量をオフセットするためにも使用することができる（例えば、ロバスト・フィッティング法が使用される場合、ｎはより小さくてもよい）。場合によっては、例えば、多重線形回帰（ロバスト・フィッティング法の使用を含む）などの一般的な線形モデル化アプローチを使用して、パラメータの一部またはすべての合同フィッティングが行われてもよい。
【００６６】
図５の観察に基づき、表２は、一実施形態に係るフィッティング技法を概説する。
【表２】
2020504820.tif 000005
【００６７】
図５は、パラメータ（ｅ、ｉ、ａ）の変動が小さいことを示している。そのため、１つの実施形態では、回帰フィッティングはパラメータ（Ω、ω、Ｍ）のみに行われ、パラメータ（ｅ、ｉ、ａ）は最も最近の（すなわちｔｎの時の）過去の値に設定される。これは、例えば、これらの値がゆっくりと変化していたことを示す過去のデータの以前の分析に基づいて予め決定されてよく、または過去のデータの分析の一部として動的に行われてもよい。例えば、フィットされるパラメータの閾値が予め定められてよく、フィットされた値がこれらの事前に定義された変数よりも小さい場合、パラメータは最も最近の（すなわちｔｎの時の）過去の値に設定される。
【００６８】
ケプラー軌道パラメータの動的推定値を取得することに加えて、下記で説明されるエポック演算は、衛星の平均運動についてのモデルを必要とする。これは、ここでも、他のパラメータに関して上記で説明したのと同様にして、例えば過去のデータに対する線形最小二乗フィットを使用して、フィッティング（または回帰）によって取得されてよい。ρが平均運動である場合、１つの実施形態では、線形最小二乗フィットで、最良フィットの線ｍρｔ＋ｃρを記述するｍρ、ｃρが見つかる。１つの実施形態では、平均運動の代わりに衛星の軌道番号がフィットされ、平均運動は、時間に対する軌道番号の導関数として導出される。
【００６９】
動的軌道モデル２０は、フィッティング・ステップ１４中に求められたモデル・パラメータ１８からなる（平均運動を含む）。動的軌道モデルは、フィットされたパラメータの誤差（推定も行う場合）を含むこともあり得る。いくつかの実施形態では、動的軌道モデルによってフィットされる軌道パラメータのセットは、軌道にある最小制御衛星の運動をモデル化するのに十分なパラメータの任意のセットである。例えば、一部の軌道パラメータは、大幅には変動しない可能性があり、そのため固定値に設定され、モデルによってフィットされなくてよい。
【００７０】
動的軌道モデルが端末装置１００によって取得され、記憶されると、端末装置１００はこれを使用して取得することができる。方法は、エポック演算２６と、エポック３０における軌道パラメータの予測／生成と、その後の、エポックの前後の時刻における衛星ロケーションを予測する標準的な軌道予測器３４への入力としてのこれらの「合成」軌道要素（すなわち合成ＴＬＥ）の提供と、の３つの主要な処理ステップを含む。本発明ではＴＬＥを軌道要素の汎用的な記述子として使用しており、そのため、必ずしもＮＯＲＡＤの２行要素だけはなく、軌道パラメータの任意のセットに当てはまることに留意されたい。
【００７１】
端末装置は、十分に正確にＵＴＣと同期されている、一日のうちの時刻である参照時刻２４を提供する時刻源２２へのアクセスを有する。例えば、端末装置は、ＧＰＳモジュール、周期的にＵＴＣと同期される安定ローカル発振器、または端末装置１００の内部および外部の両方にある何らかの他のタイミング参照源から時刻を取得してよい。
【００７２】
大半の軌道予測器（ＳＧＰ４等）は、入力として、エポックと呼ばれる特定の時刻インスタンスにおける軌道および衛星位置を記述する軌道パラメータのセット（すなわちエフェメリス・データ）を受け取る。多くの場合、エポックが、昇交点の時刻など軌道上の特定の事象と位置合わせされることが利便である。通例は、ただし常にそうではないが、エポックを参照時刻の近傍に設定することを望む。そのため、１つの実施形態では、エポック時刻２８は、衛星通過予測の対象となる時間期間の直前に発生する昇交点通過の時刻となるように演算される。
【００７３】
ｔｎを最も最近の過去の軌道パラメータのセットのエポック時刻とし、ｍρおよびｃρを、過去の平均運動データρに対する最良の（ロバストな）線形フィットの勾配パラメータおよびオフセット・パラメータとする。図６は、図５Ａ〜５Ｆに示されるものと同じ衛星の平均運動のプロットを示す。実線６１は過去のデータ（３０日分の実際のＴＬＥデータ）を示し、破線６２は、フィットされたパラメータｃρ＝１４．２２６９９７８７回転／日、および一日当たりｍρ＝１．６５６６４×１０−６回転の二乗での、Ｔｈｅｉｌ−Ｓｅｎのロバスト線形フィットである。
【００７４】
ρが周回数／時間の単位で表されるとする（ｔと同じ単位での時間）。衛星が次のＫ回の完全な軌道を完了するための所要時間Ｔを推定したいとする。すると、以下が示されることができる。
【数２】
2020504820.tif 000006
【００７５】
同様の式を他のフィッティング・モデルに対して取得することが可能である。この積分を行い、Ｔについて解くと、
【数３】
2020504820.tif 000007
が得られ、ここでは非負の解を取る。ｍρ＝０の場合、ロピタルの定理を使用してＴ＝Ｋ／ρを回復できることに留意されたい。
【００７６】
式（４）を使用すると、所与の整数の軌道周回数Ｋに対応するｔｎへのオフセットである適切なＴを求めることは単純である。そのため、１つの実施形態では、勾配パラメータｍρおよびオフセットｃρを取得するために、時間期間ｔ＝（ｔ１，ｔ２．．．，ｔｎ）にわたる平均運動ρについてのρ＝ｍρｔ＋ｃρの線形フィットが行われ、これらはＴを求めるために式（４）に提供される。あるいは、時間の関数としての過去の軌道番号についての線形フィットを行い、このフィットを使用して直接Ｋを推定することにより、完全な軌道の数Ｋを推定することができる。次いで、平均運動ρの推定値が時間に対する軌道番号の導関数として取得され、次いでこれを使用して、例えば式（４）を使用して所要時間Ｔを求める。逆に、所与の目標時間オフセットＴ’については、式（４）を使用して、何らかの整数の軌道周回数に対応する最も近いＴ（必要に応じてＴ’の前または後）を容易に見つけることができる。
【００７７】
すると、エポックをＴｅｐｏｃｈ＝ｔｎ＋Ｔに設定することができる。例えば、勾配パラメータｍρおよびオフセットｃρを取得した後に、衛星が次のＫ回の完全な軌道を完了するための所要時間Ｔを、Ｔが対象のタイマ期間Ｔ’の開始の直前になるまで、ある範囲のＫの値に対して式（４）を使用して推定することができる。すなわち、ＴＫ≦Ｔ’＜ＴＫ＋１となるようにＫを探索することができ、その際にＴ（＝ＴＫ）を取得することができる。これは、初期値（例えば１または１０００）で開始し、次いで、時刻が対象の時間期間の開始の前であるか後であるかに応じて、Ｋを増分または減分することによって行うことができる。あるいは、式（３）においてＴ＝Ｔ’を使用することによって式（３）を使用することができ、次いでＫの値を最も近い整数に切り捨て（完全な軌道の数を得るため）、次いで式（４）を使用してＴを見つける。図７は、本発明の例示的衛星に対するこの工程の結果を示している。エポック時刻は、過去のエフェメリス・データから、図６で取得されたフィットを使用して演算された。実線７１は、予測されたエポック時刻と実際のＴＬＥエポック時刻との間の誤差を示している。これより、誤差は約３００日後にわずか１分であり、一年後の将来の誤差はわずかおよそ３分であることを見て取ることができる。この誤差の大きさは最終的な衛星通過予測の誤差を支配することになり、よってこれは非常に良好なパフォーマンスを示している。
【００７８】
次いで、軌道要素予測器３０を使用して、時刻インスタンスＴｅｐｏｃｈについての予測される軌道要素のセット（すなわちエフェメリス・データまたはＴＬＥ）を作成する。この予測器３０は、入力として、以前に演算されたＴｅｐｏｃｈ２８、ならびに上記で演算された軌道モデル・パラメータ３４を受け取り、外挿を使用して、所望の参照時間に対応する所望の軌道パラメータを取得する。例えば、本実施形態では以下を演算することになる。
ａ［Ｔｅｐｏｃｈ］＝ｍａＴｅｐｏｃｈ＋ｃａ      式（５）
ｉ［Ｔｅｐｏｃｈ］＝ｃｅｓｉｎ（νｅＴｅｐｏｃｈ＋θｉ）      式（６）
Ω［Ｔｅｐｏｃｈ］＝ｍΩＴｅｐｏｃｈ＋ｃΩ      式（７）
ω［Ｔｅｐｏｃｈ］＝ｍωＴｅｐｏｃｈ＋ｃω      式（８）
Ｍ［Ｔｅｐｏｃｈ］＝ｍＭＴｅｐｏｃｈ＋ｃＭ      式（９）
【００７９】
他のモデル（例えば、より高次の多項式）およびそれらに対応するフィットされたパラメータを使用して、同様の予測を見つけることができる。例えば、カルマン予測器が、カルマン・フィッティングの工程によって求められたパラメータと共に用いられてもよい。
【００８０】
図８Ａ〜８Ｆは、本発明の例示的衛星の軌道パラメータについての時間に伴う絶対誤差のプロットを示している。横軸は、最も最近の過去のデータ点以降の、日単位で測定された時間である。直前の３０日間分のデータを使用してフィットを作成した。図８Ａは、準主軸ａの絶対誤差の曲線８１を示し、図８Ｂは、離心率ｅの絶対誤差の曲線８２を示し、図８Ｃは、傾きｉの絶対誤差の曲線８３を示し、図８Ｄは、赤経昇交点Ωの絶対誤差の曲線８４を示し、図８Ｅは、近地点の引数ωの絶対誤差の曲線８８を示し、図８Ｆは、平均近点離角Ｍの絶対誤差の曲線８６を示している。
【００８１】
図８Ａ〜８Ｆは、単純な線形フィットを使用して取得された結果が充分なものであることを例示しているが、前述したようにｅ、ｉ、ω、およびＭに正弦関数をフィットさせることにより、改善された（すなわち低減した誤差）が取得されよう。また、ωおよびＭの誤差は相殺される傾向があることにも留意されたい（Ｍはωからのオフセットとして測定されるため）。これらの図は、３５０日前後に宇宙船に加えられた小さな調節の影響も明らかにしている。このことは、ＴＬＥデータが通例有効である通常の数時間から数日と比べて、エフェメリス・データの有効時間の大幅な延長を意味する。
【００８２】
これで更新されたまたは合成の軌道要素を取得したので、ＳＧＰ４またはＪ２（Ｓｈｒａｄｄｈａ  Ｇｕｐｔａ，  Ｍ．Ｘａｖｉｅｒ  Ｊａｍｅｓ  Ｒａｊ，  ａｎｄ  Ｒａｍ  Ｋｒｉｓｈａｎ  Ｓｈａｒｍａ．  Ｓｈｏｒｔ−ｔｅｒｍ  ｏｒｂｉｔ  ｐｒｅｄｉｃｔｉｏｎ  ｗｉｔｈ  Ｊ２  ａｎｄ  ｍｅａｎ  ｏｒｂｉｔａｌ  ｅｌｅｍｅｎｔｓ．  Ｉｎｔｅｒｎａｔｉｏｎａｌ  Ｊｏｕｒｎａｌ  ｏｆ  Ａｓｔｒｏｎｏｍｙ  ａｎｄ  Ａｓｔｒｏｐｈｙｓｉｃｓ，  １：１３５｛１４６、２０１１）などの軌道予測器モジュール３４が、軌道要素予測器３０によって作成された合成軌道要素３２を使用して、Ｔｅｐｏｃｈに近い時刻における衛星の絶対位置４２を推定する。合成軌道要素は、端末装置１００の視点からの衛星４の相対位置４４（例えば範囲、仰角、方位角）を求めるためにも使用することができる。この場合、端末装置の位置３８も必要とされ、位置３８は、ＧＰＳなどの位置参照源３６から提供されるか、または静止している端末装置１００の場合はメモリに記憶されてよい。軌道予測器３０は、端末装置１００から見た次の通過（または複数回の次の通過）の時刻を推定するためにも使用することができる。これには、この場合も、端末装置１００のおよその位置３８の知識が必要となる。
【００８３】
通過時刻を推定する際に、生じ得る誤差（例えば通過時刻プラスまたはマイナス誤差時間）の推定値を報告することも有用であり得る。例えば、デバイスがスリープからウェイクして衛星と通信するために予測通過時刻を使用する用途では、デバイスは少し早くウェイクして、少し長めにウェイク状態のままであることにより、通過時刻の不確定性を考慮に入れることができる。
【００８４】
上記で解説したように、エポック演算、軌道要素の予測、および実際の軌道予測は、端末装置上で行うことができる。図９は、一実施形態に係る端末装置１００の概略図である。端末装置は、１つまたは複数のアンテナ１１２を備えたＲＦフロントなどの通信モジュール１１０と、符号化および変調を含めて、データを準備し、そのデータを無線周波リンクを通じて衛星４に送信するための関連するハードウェアとを備える。いくつかの実施形態では、通信モジュール１１０は、衛星または他の供給源からデータを受信（および復号）するようにも構成される。通信モジュールは、ＲＦフロント・エンドと、追加的なハードウェア・コンポーネントおよび／またはソフトウェア・コンポーネントとを備える。受信信号は、ＲＦフロント・エンドによってベースバンドに変換され、メッセージ・データを抽出するための復調および復号などの処理のためにベースバンド受信機に送られる。同様に、端末送信機が、ソース・ビットを受け取り、それらを使用して、変調および符号化方式の適用、ならびに送信サブキャリアおよび電力の選択など、ユーザ端末のＲＦフロント・エンドを介して送信するためのベースバンド信号を生成する。通信モジュールは、独立型モジュールもしくは基板であってよく、または他のコンポーネントと一体化されてもよい。ソフトウェアによって定義される無線実装が使用されてよく、その場合は、ＲＦフロント・エンドが受信信号をアナログ−デジタル（ＡＤＣ）変換器に提供し、変換器はスペクトル・サンプルを信号プロセッサに提供する。
【００８５】
端末装置は、プロセッサ・モジュール１２０およびメモリ１３０も備える。メモリは、軌道モデル２０を記憶すること、エポック２８を演算または推定すること、軌道要素３２の予測を生成すること、およびそれらを軌道予測器に渡して通過時刻または衛星位置を取得することなど、本明細書に記載される方法をプロセッサに実施させるためのソフトウェア命令またはソフトウェア・モジュールを含む。メモリは、所望の時間に（例えば予測される衛星通過時刻の間に）端末をウェイクアップさせるアラーム・モジュールなど、他の機能に関するモジュールを記憶するためにも使用されてよい。電源、クロック、センサ・プラットフォーム等の他のコンポーネントも端末装置に含まれてよい。
【００８６】
メモリ１３０は、動的軌道モデル２０およびモデル・パラメータ１８を記憶するためにも使用される。端末装置は、過去のエフェメリス・データ１２を受け取り、更新された動的軌道モデル２０を取得（および記憶）するためのフィッティング工程を行う。あるいは、動的軌道モデル３０（および関連するモデル・パラメータ１８）は、端末装置に提供されてもよく、端末装置はモデルをメモリ１２０に記憶する。これは、以前に記憶された値を更新すること、または以前に記憶されたモデルを完全に置換する（もしくはモデルを初めて保存する）ことを伴い得る。
【００８７】
過去のエフェメリス・データ１２または軌道モデル２０は、通信モジュール１１０を介して提供されてよい。これは、スケジュールされた時刻に取得されてよく、その場合、デバイスはウェイクアップして送信を受信するか、またはデータの送信源の範囲内にある時に時機を見て取得されてもよい（例えば端末装置が可動性である場合）。あるいは、ＵＳＢインターフェースなどの物理インターフェースが提供されて、サービス時または保守時にこのデータがデバイスに物理的に転送（またはアップロード）されるようにしてもよい。端末は、位置および時刻の推定値を提供するため（すなわち図４のブロック２２および３６の両方）に使用できるＧＰＳ受信機１４０を有してよい。あるいは、端末装置は、通信モジュール１１０を介してタイミング情報を受信してもよく、または、端末装置は、例えばサービス時または保守時にＵＴＣと周期的に同期される安定した内蔵クロックを含んでよい。
【００８８】
本願は、端末装置が使用するための衛星の拡張された衛星エフェメリス・データを生成する方法およびシステムの実施形態を記載する。これらの実施形態は、最小制御衛星（無制御衛星を含む）に使用するのに特に適する。この方法は、現在のＴＬＥおよびＳＧＰ４などの軌道予測器を使用して現在可能なよりもはるかに長い時間枠、例えば数時間／数日から数カ月／一年またはそれ以上にわたって衛星通過時刻または衛星位置を正確に予測するために使用することができる。大まかには、過去の軌道パラメータ・データ（例えば過去のＴＬＥから取られる）を使用して、軌道パラメータの将来の値の予測器を構築する（すなわち動的軌道モデルを構築する）。これは、端末上かまたはオフラインのいずれかで行うことができ、モデルの結果が端末に提供される。
【００８９】
動的軌道パラメータ予測モデルが取得されると、モデルは端末装置上に記憶される。そして、所望のエポックに、事前に記憶された軌道パラメータ予測モデルを使用して、そのエポックに対応する衛星の軌道パラメータの推定値を作成する（すなわち合成または予測ＴＬＥ）。その後、これらの予測された軌道パラメータがＳＧＰ４などの軌道予測器への入力として使用されて、例えば衛星通過時刻またはいつ衛星が端末に対して特定の位置に来るかを判定するために、所望の時刻における衛星のロケーションの推定値を作成する。本明細書に示されるように、モデルが生成された時から数百日にわたって衛星ロケーションの正確な推定値を作成することができる。これにより、端末は衛星通過時刻を正確に推定できるため、低電力が分配される端末装置が効率的に動作することが可能となる。これにより、頻繁なエフェメリス・データを受信するためにウェイクアップする、または衛星通過を逃さないようにするために長期間にわたってウェイクアップすることなどによって、端末が不必要に電力を浪費しないことが保証される。
【００９０】
本明細書に記載される方法の実施形態は、衛星通過タイミングまたは衛星ロケーションの適度に正確な推定値を見つけることができる期間を拡張することに関心が持たれ得る種々の衛星通信システムで使用することができる。これらには、エフェメリス・データを定期的に更新することが難しいかまたは不可能であり得る、低帯域幅または断続的な通信リンクを備えるシステム、例えば一方向通信デバイスなど、エフェメリス・データを遠隔から更新する能力を持たないシステム、データ通信コスト、バッテリ使用、または他のコストのいずれかの点から、定期的にエフェメリス・データを更新するコストが高すぎるシステム、およびデバイスに送信されるエフェメリス・データの通信妨害または電波妨害によって影響され得るシステム、が含まれる。
【００９１】
当業者は、情報および信号は様々なテクノロジーおよび技術のいずれかを使用して表され得ることを理解されよう。例えば、上記の説明全体を通じて参照され得るデータ、命令、コマンド、情報、信号、ビット、シンボル、およびチップは、電圧、電流、電磁波、磁場もしくは磁性粒子、光学場もしくは光学粒子、またはそれらの任意の組み合わせによって表されてよい。
【００９２】
当業者はさらに、本明細書に開示される実施形態との関連で記載される様々な説明用の論理ブロック、モジュール、回路、およびアルゴリズム・ステップは、電子ハードウェア、コンピュータ・ソフトウェアもしくは命令、または両者の組み合わせによって実施されてよいことを認識されよう。ハードウェアとソフトウェアとの互換性を明確に例示するために、様々な説明用のコンポーネント、ブロック、モジュール、回路、およびステップについて、上記では概してそれらの機能性の点から説明した。そのような機能性がハードウェアとして実施されるかソフトウェアとして実施されるかは、特定の用途およびシステム全体に課される設計の制約に依存する。当業者は、記載される機能性を特定の用途ごとに異なる方式で実施してよいが、そのような実施の決定は、本発明の範囲からの逸脱を生じさせるものとは解釈すべきでない。
【００９３】
本明細書に開示される実施形態に関連して記載される方法またはアルゴリズムのステップは、直接ハードウェアとして、プロセッサによって実行されるソフトウェア・モジュールとして、またはそれら２つの組み合わせとして具現化されてよい。ハードウェア実装の場合、処理は、１つまたは複数の特定用途集積回路（ＡＳＩＣ）、デジタル信号プロセッサ（ＤＳＰ）、デジタル信号処理デバイス（ＤＳＰＤ）、プログラム可能論理デバイス（ＰＬＤ）、フィールド・プログラマブル・ゲートアレイ（ＦＰＧＡ）、プロセッサ、コントローラ、マイクロコントローラ、マイクロプロセッサ、もしくは本明細書に記載される機能を行うように設計された他の電子ユニット、またはそれらの組み合わせの中で実施されてよい。
【００９４】
いくつかの実施形態では、プロセッサ・モジュール１１０は、方法のステップのいくつかを行うように構成された１つまたは複数の中央演算処理装置（ＣＰＵ）を備える。同様に、端末装置に供給される軌道モデルを生成するためにコンピューティング装置が使用されてもよく、コンピューティング装置は１つまたは複数のＣＰＵを備えてよい。ＣＰＵは、入出力インターフェースと、演算論理ユニット（ＡＬＵ）と、制御ユニットと、入出力インターフェースを通じて入力装置および出力装置と通信状態にあるプログラム・カウンタ要素とを備えてよい。入出力インターフェースは、予め定められた通信プロトコル（例えばＢｌｕｅｔｏｏｔｈ、Ｚｉｇｂｅｅ、ＩＥＥＥ８０２．１５、ＩＥＥＥ８０２．１１、ＴＣＰ／ＩＰ、ＵＤＰ等）を使用して別のデバイス内の同等の通信モジュールと通信するためのネットワーク・インターフェースおよび／または通信モジュールを備えてよい。コンピューティング装置または端末装置は、単一のＣＰＵ（コア）もしくは複数のＣＰＵ（マルチ・コア）、または複数のプロセッサを備えてよい。コンピューティング装置または端末装置は、並列プロセッサ、ベクトル・プロセッサを使用してよく、または分散コンピューティング・デバイスであってもよい。メモリが、動作可能にプロセッサに結合され、ＲＡＭおよびＲＯＭコンポーネントを含んでよく、またデバイスまたはプロセッサ・モジュールの内部または外部に設けられてよい。メモリは、オペレーティング・システムおよび追加的なソフトウェア・モジュールまたは命令を記憶するために使用されてよい。プロセッサは、メモリに記憶されたソフトウェア・モジュールまたは命令をロードし、実行するように構成されてよい。
【００９５】
コンピュータ・プログラム、コンピュータ・コード、または命令とも呼ばれるソフトウェア・モジュールは、いくつかのソース・コードまたはオブジェクト・コード・セグメントまたは命令を含んでいてよく、ＲＡＭメモリ、フラッシュ・メモリ、ＲＯＭメモリ、ＥＰＲＯＭメモリ、レジスタ、ハード・ディスク、取り外し可能ディスク、ＣＤ−ＲＯＭ、ＤＶＤ−ＲＯＭ、Ｂｌｕ−ｒａｙディスク、または任意の他の形態のコンピュータ可読媒体などの任意のコンピュータ可読媒体に常駐してよい。いくつかの態様では、コンピュータ可読媒体は、非一時的なコンピュータ可読媒体（例えば有形の媒体）を備えてよい。また、他の態様には、コンピュータ可読媒体は、一時的なコンピュータ可読媒体（例えば信号）を備えてよい。上記の組み合わせもコンピュータ可読媒体の範囲に含められるべきである。別の態様では、コンピュータ可読媒体は、プロセッサと一体であってよい。プロセッサおよびコンピュータ可読媒体は、ＡＳＩＣまたは関係するデバイスに常駐してよい。ソフトウェア・コードはメモリ・ユニットに記憶されてよく、プロセッサはそれらを実行するように構成されてよい。メモリ・ユニットは、プロセッサの内部またはプロセッサの外部に実施されてよく、外部の場合、メモリ・ユニットは当技術分野で知られる様々な手段を介してプロセッサに動作可能に結合されることができる。
【００９６】
さらに、本明細書に記載される方法および技術を行うモジュールおよび／または他の適切な手段は、コンピューティング・デバイスによってダウンロードされるか、および／またはその他の形で取得できることを認識すべきである。例えば、そのようなデバイスをサーバに結合して、本明細書に記載される方法を行うための手段の移送を容易にすることができる。あるいは、本明細書に記載される様々な方法は、記憶手段（例えばＲＡＭ、ＲＯＭ、コンパクト・ディスク（ＣＤ）またはフロッピー・ディスクなどの物理記憶媒体等）を介して提供することができ、それによりコンピューティング・デバイスは、記憶手段をデバイスに結合するかまたは提供すると、様々な方法を取得することができる。さらに、本明細書に記載される方法および技術をデバイスに提供するための任意の他の適切な技術を利用することができる。
【００９７】
本明細書に開示される方法は、記載される方法を実現するための１つまたは複数のステップまたは動作を含む。これらの方法ステップおよび／または動作は、特許請求の範囲から逸脱することなく互いに入れ替えられてよい。換言すると、ステップまたは動作の特定の順序が指定されない限り、特定のステップおよび／または動作の順序および／または使用は、特許請求の範囲から逸脱することなく変更されてよい。
【００９８】
本明細書で使用される場合、用語「推定する（ｅｓｔｉｍａｔｉｎｇ）」または「判定する／求める（ｄｅｔｅｒｍｉｎｉｎｇ）」は、幅広い動作を包含する。例えば、「推定する（ｅｓｔｉｍａｔｉｎｇ）」または「判定する／求める（ｄｅｔｅｒｍｉｎｉｎｇ）」は、計算すること、演算すること、処理すること、導出すること、調査すること、参照すること（例えばテーブル、データベース、または別のデータ構造内の参照）、確定すること等を含み得る。また、「推定する（ｅｓｔｉｍａｔｉｎｇ）」または「判定する／求める（ｄｅｔｅｒｍｉｎｉｎｇ）」は、受信すること（例えば情報の受信）、アクセスすること（例えばメモリ内のデータへのアクセス）等を含み得る。また、「判定する／求める（ｄｅｔｅｒｍｉｎｉｎｇ）」は、解決すること、選択すること、選ぶこと、立証すること等を含み得る。
【００９９】
明細書および以下の特許請求の範囲の全体を通じて、文脈が別様に必要としない限り、単語「〜を備える（ｃｏｍｐｒｉｓｅ）」および「〜を含む（ｉｎｃｌｕｄｅ）」ならびに「〜を備える（ｃｏｍｐｒｉｓｉｎｇ）」および「〜を含む（ｉｎｃｌｕｄｉｎｇ）」などの変化形は、述べられる整数または整数の群の包含を示唆するが、他の整数または整数の群の排除は示唆しないものと理解されるものとする。
【０１００】
本明細書における従来技術への参照は、そのような従来技術が通常の一般的知識の一部をなすという任意形態の示唆の認知ではなく、またそのようにみなすべきではない。
【０１０１】
本開示は、その使用において記載される１つまたは複数の特定の用途に制約されないことが当業者により認識されよう。また、本開示は、本明細書に記載または図示される特定の要素および／または特徴に関してその好ましい実施形態に制約されることもない。本開示は、開示される１つまたは複数の実施形態に制限されず、以下の特許請求の範囲によって述べられ、定義される範囲から逸脱することなく、多数の再構成、変更、および置換が可能であることが認識されよう。本明細書で使用される場合、項目のリスト「のうちの少なくとも１つ」を言及する語句は、個々の単一の部材を含む、それらの項目の任意の組み合わせを指す。一例として、「ａ、ｂ、またはｃのうちの少なくとも１つ」は、ａ、ｂ、ｃ、ａ−ｂ、ａ−ｃ、ｂ−ｃ、およびａ−ｂ−ｃを網羅することが意図される。
【図１】
2020504820.tif 000009
【図２】
2020504820.tif 000010
【図３】
2020504820.tif 000011
【図４】
2020504820.tif 000012
【図５Ａ】
2020504820.tif 000013
【図５Ｂ】
2020504820.tif 000014
【図５Ｃ】
2020504820.tif 000015
【図５Ｄ】
2020504820.tif 000016
【図５Ｅ】
2020504820.tif 000017
【図５Ｆ】
2020504820.tif 000018
【図６】
2020504820.tif 000019
【図７】
2020504820.tif 000020
【図８Ａ】
2020504820.tif 000021
【図８Ｂ】
2020504820.tif 000022
【図８Ｃ】
2020504820.tif 000023
【図８Ｄ】
2020504820.tif 000024
【図８Ｅ】
2020504820.tif 000025
【図８Ｆ】
2020504820.tif 000026
【図９】
2020504820.tif 000027
【国際調査報告】
2020504820.tif 000028
2020504820.tif 000029
2020504820.tif 000030