JPH1020010A - 衛星追尾方式 - Google Patents
衛星追尾方式Info
- Publication number
- JPH1020010A JPH1020010A JP16921696A JP16921696A JPH1020010A JP H1020010 A JPH1020010 A JP H1020010A JP 16921696 A JP16921696 A JP 16921696A JP 16921696 A JP16921696 A JP 16921696A JP H1020010 A JPH1020010 A JP H1020010A
- Authority
- JP
- Japan
- Prior art keywords
- satellite
- function
- coefficient
- antenna
- stored
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Landscapes
- Radio Relay Systems (AREA)
Abstract
(57)【要約】
【課題】 衛星の追尾精度を向上させる衛星追尾方式を
提供する。 【解決手段】 衛星軌跡情報記憶部11に蓄積されてい
る衛星軌跡情報を用いて、回帰分析処理部12は、線形
回帰分析により衛星の軌道関数を求める。ここで求めた
衛星の軌道関数の回帰係数は、回帰係数記憶部13に順
次記憶される。回帰係数評価部14は、回帰係数記憶部
13に記憶された回帰係数を平均化して最適評価する。
衛星位置・指令角計算処理部15は、最適評価された回
帰係数から衛星の予測位置およびアンテナ部17に与え
る指令角(方位角および仰角)を求める。アンテナ制御
部16は、衛星位置・指令角計算処理部15で求められ
た指令角にアンテナが向くようにアンテナ部17の駆動
を制御する。また、アンテナ制御部16は、一定時間毎
にステップトラック方式によりアンテナ部17を駆動し
て衛星を捕捉し、新たな衛星軌跡情報として衛星軌跡情
報記憶部11に順次記憶する。
提供する。 【解決手段】 衛星軌跡情報記憶部11に蓄積されてい
る衛星軌跡情報を用いて、回帰分析処理部12は、線形
回帰分析により衛星の軌道関数を求める。ここで求めた
衛星の軌道関数の回帰係数は、回帰係数記憶部13に順
次記憶される。回帰係数評価部14は、回帰係数記憶部
13に記憶された回帰係数を平均化して最適評価する。
衛星位置・指令角計算処理部15は、最適評価された回
帰係数から衛星の予測位置およびアンテナ部17に与え
る指令角(方位角および仰角)を求める。アンテナ制御
部16は、衛星位置・指令角計算処理部15で求められ
た指令角にアンテナが向くようにアンテナ部17の駆動
を制御する。また、アンテナ制御部16は、一定時間毎
にステップトラック方式によりアンテナ部17を駆動し
て衛星を捕捉し、新たな衛星軌跡情報として衛星軌跡情
報記憶部11に順次記憶する。
Description
【0001】
【発明の属する技術分野】本発明は、衛星追尾方式に関
し、特に静止衛星の自動追尾システムにおいて過去の衛
星軌跡情報を利用して現在の衛星位置を予測する衛星追
尾方式に関するものである。
し、特に静止衛星の自動追尾システムにおいて過去の衛
星軌跡情報を利用して現在の衛星位置を予測する衛星追
尾方式に関するものである。
【0002】
【従来の技術】従来の静止衛星の衛星追尾方式は、静止
衛星から伝搬されてくる電波の受信レベルを監視しなが
らアンテナを少しずつステップ状に動かして、受信レベ
ルが最も高くなるような位置を探すステップトラック方
式が主流であった。
衛星から伝搬されてくる電波の受信レベルを監視しなが
らアンテナを少しずつステップ状に動かして、受信レベ
ルが最も高くなるような位置を探すステップトラック方
式が主流であった。
【0003】しかし、ステップトラック方式は、受信波
を利用する方式であるため、受信レベルの低下や受信レ
ベルに大きな変動が起こると衛星追尾が中断する。この
衛星追尾の中断中に、静止衛星が大きく移動し、受信レ
ベルが復帰しないと自動追尾が不可能になるという問題
点があった。
を利用する方式であるため、受信レベルの低下や受信レ
ベルに大きな変動が起こると衛星追尾が中断する。この
衛星追尾の中断中に、静止衛星が大きく移動し、受信レ
ベルが復帰しないと自動追尾が不可能になるという問題
点があった。
【0004】そこで、バックアップ用としてメモリ追尾
方式あるいはIESS−412プログラム追尾方式が用
いられてきた。
方式あるいはIESS−412プログラム追尾方式が用
いられてきた。
【0005】メモリ追尾方式は、静止衛星の動きが1日
周期であることを利用して1周期前の角度情報またはこ
の1周期前の角度情報に2周期前の角度情報を加えて静
止衛星の位置を計算するものである。しかし、毎日同時
刻に受信レベルの変動が発生するシンチレーション現象
の場合には1周期前や2周期前の角度情報が取り込めな
い。従って、この場合には、メモリ追尾方式では対応で
きないという問題点があった。。
周期であることを利用して1周期前の角度情報またはこ
の1周期前の角度情報に2周期前の角度情報を加えて静
止衛星の位置を計算するものである。しかし、毎日同時
刻に受信レベルの変動が発生するシンチレーション現象
の場合には1周期前や2周期前の角度情報が取り込めな
い。従って、この場合には、メモリ追尾方式では対応で
きないという問題点があった。。
【0006】IESS−412プログラム追尾方式は、
受信波を利用せずに、静止衛星の11軌道要素を用いて
静止衛星の位置を計算する。しかし、この方式では定期
的に衛星運用センターから11軌道要素を入手しなけれ
ばならないという問題点があった。
受信波を利用せずに、静止衛星の11軌道要素を用いて
静止衛星の位置を計算する。しかし、この方式では定期
的に衛星運用センターから11軌道要素を入手しなけれ
ばならないという問題点があった。
【0007】このような問題点を解消するために、本願
出願人は特願平5−68986にて最小二乗法による線
形回帰分析を行って静止衛星の現在位置を予測する衛星
追尾方式を提案している。以下、この衛星追尾方式の概
要について説明する。
出願人は特願平5−68986にて最小二乗法による線
形回帰分析を行って静止衛星の現在位置を予測する衛星
追尾方式を提案している。以下、この衛星追尾方式の概
要について説明する。
【0008】図4は、この衛星追尾方式の機能ブロック
図である。図5は、この衛星追尾方式のフローチャート
である。
図である。図5は、この衛星追尾方式のフローチャート
である。
【0009】まず、過去の衛星軌跡情報が衛星軌跡情報
記憶部21に少なくとも1周期分記憶されているかどう
かを判定する(S21)。衛星軌跡情報が1周期分まで
記憶されていなければ、アンテナ制御部26はステップ
トラック方式によりアンテナ部27を駆動して衛星の位
置の捕捉を行い(S27)、その位置を衛星軌跡情報と
して衛星軌跡情報記憶部21に記憶する(M21)。一
方、衛星軌跡情報が1周期分以上記憶されているとき
は、衛星軌跡情報記憶部21から一定期間分(例えば、
1周期分)の衛星軌跡情報を読み出す(S22)。そし
て、回帰分析処理部22は最小二乗法による線形回帰分
析を用いて衛星の軌道関数を求める(S23)。
記憶部21に少なくとも1周期分記憶されているかどう
かを判定する(S21)。衛星軌跡情報が1周期分まで
記憶されていなければ、アンテナ制御部26はステップ
トラック方式によりアンテナ部27を駆動して衛星の位
置の捕捉を行い(S27)、その位置を衛星軌跡情報と
して衛星軌跡情報記憶部21に記憶する(M21)。一
方、衛星軌跡情報が1周期分以上記憶されているとき
は、衛星軌跡情報記憶部21から一定期間分(例えば、
1周期分)の衛星軌跡情報を読み出す(S22)。そし
て、回帰分析処理部22は最小二乗法による線形回帰分
析を用いて衛星の軌道関数を求める(S23)。
【0010】次に、新しい衛星軌跡情報を記憶する時刻
(例えば、1時間毎)となっているかどうかを判定する
(S24)。記憶時刻となっていなければ、衛星位置・
指令角計算処理部25は、ステップS23で求めた軌道
関数を使用して現在の衛星位置を計算し、アンテナ部2
7に与える指令角(方位角および仰角)を計算する(S
25)。ここで求めた指令角をアンテナ制御部26に入
力し、アンテナ制御部26からの指示によりアンテナ部
27を駆動する(S26)。
(例えば、1時間毎)となっているかどうかを判定する
(S24)。記憶時刻となっていなければ、衛星位置・
指令角計算処理部25は、ステップS23で求めた軌道
関数を使用して現在の衛星位置を計算し、アンテナ部2
7に与える指令角(方位角および仰角)を計算する(S
25)。ここで求めた指令角をアンテナ制御部26に入
力し、アンテナ制御部26からの指示によりアンテナ部
27を駆動する(S26)。
【0011】ステップS24の処理で新しい衛星軌跡情
報の記憶時刻となっていると判断したときは、アンテナ
制御部26はステップトラック方式によりアンテナ部2
7を駆動して衛星の位置の捕捉を行い(S27)、その
位置を新たな衛星軌跡情報として衛星軌跡情報記憶部2
1に追加して記憶する(M21)。
報の記憶時刻となっていると判断したときは、アンテナ
制御部26はステップトラック方式によりアンテナ部2
7を駆動して衛星の位置の捕捉を行い(S27)、その
位置を新たな衛星軌跡情報として衛星軌跡情報記憶部2
1に追加して記憶する(M21)。
【0012】以上の動作を繰り返し、常に最新の衛星軌
跡情報を使用することにより、衛星の予測追尾を行うも
のである。
跡情報を使用することにより、衛星の予測追尾を行うも
のである。
【0013】上述の衛星の予測追尾に使用している軌道
関数は、地上局から見た衛星の動きを多項式で近似した
関数である。ところが、実際には静止衛星は静止軌道上
を正確に周回するということはできない。従って、蓄積
された衛星軌跡情報には、この誤差が必然的に含まれて
いる。また、地上局の位置によって見た目上の衛星の軌
跡は様々に変化する。このため、ある特定の位置の地上
局では衛星の軌道関数に本来は係数が0となるべきであ
るが、前述の誤差のために実際にはある値を持っている
不必要な項が含まれてしまっている。この不必要な項を
含んだ軌道関数を以後の衛星の追尾に適用すると、かな
りの大きさで誤差の影響が生じ、予測追尾の精度を悪化
させてしまうという問題点があった。
関数は、地上局から見た衛星の動きを多項式で近似した
関数である。ところが、実際には静止衛星は静止軌道上
を正確に周回するということはできない。従って、蓄積
された衛星軌跡情報には、この誤差が必然的に含まれて
いる。また、地上局の位置によって見た目上の衛星の軌
跡は様々に変化する。このため、ある特定の位置の地上
局では衛星の軌道関数に本来は係数が0となるべきであ
るが、前述の誤差のために実際にはある値を持っている
不必要な項が含まれてしまっている。この不必要な項を
含んだ軌道関数を以後の衛星の追尾に適用すると、かな
りの大きさで誤差の影響が生じ、予測追尾の精度を悪化
させてしまうという問題点があった。
【0014】なお、本願に関連し、運動している物体の
追尾精度を高める技術として、特開昭58−10987
0に記載の目標追尾予測システムと、特開平1−182
776に記載の速度予測回路がある。
追尾精度を高める技術として、特開昭58−10987
0に記載の目標追尾予測システムと、特開平1−182
776に記載の速度予測回路がある。
【0015】まず、特開昭58−109870に記載の
目標追尾システムは、レーダ・センサー系と光学センサ
ー系とからの情報をもとに1つの測的プロセッサにより
運動している物体の追尾を継続維持しつつ、別の測的プ
ロセッサにより運動している物体の変進運動を予測し、
過去のデータ軌跡から曲率を刻々と算出するものであ
る。これにより、物体の現在の運動ベクトルの直線延長
線上にその物体の未来位置を求めると共に、運動の曲率
を考慮した未来位置を求めることができる。
目標追尾システムは、レーダ・センサー系と光学センサ
ー系とからの情報をもとに1つの測的プロセッサにより
運動している物体の追尾を継続維持しつつ、別の測的プ
ロセッサにより運動している物体の変進運動を予測し、
過去のデータ軌跡から曲率を刻々と算出するものであ
る。これにより、物体の現在の運動ベクトルの直線延長
線上にその物体の未来位置を求めると共に、運動の曲率
を考慮した未来位置を求めることができる。
【0016】しかし、この技術によっては、運動してい
る物体の捕捉に一旦失敗すると、その後に物体の捕捉が
可能となったとしても、連続した軌跡を求めることがで
きない。即ち、常に運動している物体が捕捉できない限
り、追尾を継続することができない。
る物体の捕捉に一旦失敗すると、その後に物体の捕捉が
可能となったとしても、連続した軌跡を求めることがで
きない。即ち、常に運動している物体が捕捉できない限
り、追尾を継続することができない。
【0017】これに対し、特開平1−182776に記
載の速度予測回路は、運動する物体の位置の検出が良好
なときに追尾により得られた位置情報および速度情報を
順次格納する。そして、一定時間物体の位置を検出でき
なくなっても、現在の物体の位置を検出できたときに、
過去の位置情報を原点とし、そのときの速度方向を0度
とする直行座標に現在の位置情報を座標変換する。即
ち、この技術ではしばらくの間物体の捕捉ができなくな
っても、その後に捕捉することができれば、追尾を継続
することができる。
載の速度予測回路は、運動する物体の位置の検出が良好
なときに追尾により得られた位置情報および速度情報を
順次格納する。そして、一定時間物体の位置を検出でき
なくなっても、現在の物体の位置を検出できたときに、
過去の位置情報を原点とし、そのときの速度方向を0度
とする直行座標に現在の位置情報を座標変換する。即
ち、この技術ではしばらくの間物体の捕捉ができなくな
っても、その後に捕捉することができれば、追尾を継続
することができる。
【0018】しかし、この技術では物体の運動を等速円
運動など単純な線形運動と仮定する必要があるが、物体
が仮定した線形運動とは全く異なった動きをしている場
合には追尾を継続することができない。また、過去の位
置情報と現在の位置情報とから物体の速度ベクトルを予
測するものであるが、現在の位置情報を得るのに具体的
な手段はなく、現在の位置情報が得られない場合や偶然
に得られたとしても同一の物体であると判別できない場
合は、追尾を継続することができない。
運動など単純な線形運動と仮定する必要があるが、物体
が仮定した線形運動とは全く異なった動きをしている場
合には追尾を継続することができない。また、過去の位
置情報と現在の位置情報とから物体の速度ベクトルを予
測するものであるが、現在の位置情報を得るのに具体的
な手段はなく、現在の位置情報が得られない場合や偶然
に得られたとしても同一の物体であると判別できない場
合は、追尾を継続することができない。
【0019】これに対し、衛星は所定の軌道上を所定の
周期で周回していることを前提としているので、予測追
尾にそれほど複雑な計算を必要としない。また、衛星を
捕捉できなくなる可能性も小さい。従って、上記2つの
関連技術は、衛星の追尾に適用してもコストがかかる割
には大きな効果が得られず、ほとんど実用性がない。
周期で周回していることを前提としているので、予測追
尾にそれほど複雑な計算を必要としない。また、衛星を
捕捉できなくなる可能性も小さい。従って、上記2つの
関連技術は、衛星の追尾に適用してもコストがかかる割
には大きな効果が得られず、ほとんど実用性がない。
【0020】
【発明が解決しようとする課題】本発明は、上記従来技
術の欠点を解消するためになされたものであり、追尾精
度を向上することのできる衛星追尾方式を提供すること
を目的とする。
術の欠点を解消するためになされたものであり、追尾精
度を向上することのできる衛星追尾方式を提供すること
を目的とする。
【0021】
【課題を解決するための手段】上記目的を達成するた
め、本発明の衛星追尾方式は、衛星から送られる信号の
受信および衛星の位置の捕捉を行うアンテナ手段と、こ
のアンテナ手段で捕捉した衛星の位置を衛星軌跡情報と
して蓄積する第1の記憶手段と、この第1の記憶手段に
蓄積された衛星軌跡情報を用いて所定時間毎に衛星の軌
道関数を求める分析手段と、この分析手段で求めた所定
時間毎の衛星の軌道関数の係数を複数時間分記憶する第
2の記憶手段と、前記第2の記憶手段に記憶された複数
時間分の軌道関数の係数を最適化して衛星の予測追尾に
使用する軌道関数の係数を決定する評価手段と、この評
価手段で決定した係数を有する軌道関数を用いて衛星の
方向を計算する計算手段と、この計算手段で計算された
衛星の方向に前記アンテナ手段を向けるように制御し、
また、一定時間毎に前記アンテナ手段が衛星の位置の捕
捉を行うように制御するアンテナ制御部とを備えたもの
とした。
め、本発明の衛星追尾方式は、衛星から送られる信号の
受信および衛星の位置の捕捉を行うアンテナ手段と、こ
のアンテナ手段で捕捉した衛星の位置を衛星軌跡情報と
して蓄積する第1の記憶手段と、この第1の記憶手段に
蓄積された衛星軌跡情報を用いて所定時間毎に衛星の軌
道関数を求める分析手段と、この分析手段で求めた所定
時間毎の衛星の軌道関数の係数を複数時間分記憶する第
2の記憶手段と、前記第2の記憶手段に記憶された複数
時間分の軌道関数の係数を最適化して衛星の予測追尾に
使用する軌道関数の係数を決定する評価手段と、この評
価手段で決定した係数を有する軌道関数を用いて衛星の
方向を計算する計算手段と、この計算手段で計算された
衛星の方向に前記アンテナ手段を向けるように制御し、
また、一定時間毎に前記アンテナ手段が衛星の位置の捕
捉を行うように制御するアンテナ制御部とを備えたもの
とした。
【0022】ここで、前記分析手段が衛星の軌道関数を
求めるには所定時間(少なくとも衛星の公転周期の1周
期分)の衛星軌跡情報が前記第1の記憶手段に蓄積され
ていなければならない。このため、前記第1の記憶手段
に蓄積された衛星軌跡情報が所定時間分以上あるかどう
かを判定する判定手段を設け、この判定手段の判定結果
に応じて次のように処理を分けて行うことができる。
求めるには所定時間(少なくとも衛星の公転周期の1周
期分)の衛星軌跡情報が前記第1の記憶手段に蓄積され
ていなければならない。このため、前記第1の記憶手段
に蓄積された衛星軌跡情報が所定時間分以上あるかどう
かを判定する判定手段を設け、この判定手段の判定結果
に応じて次のように処理を分けて行うことができる。
【0023】(1)衛星の軌跡情報が所定時間分以上で
あると判定したときに、前記分析手段は衛星の軌道関数
を求める。
あると判定したときに、前記分析手段は衛星の軌道関数
を求める。
【0024】(2)衛星の軌跡情報が所定時間分未満で
あると判定したときに、前記アンテナ制御部は衛星から
の受信信号により前記アンテナ手段が衛星の位置の捕捉
を行うように制御する。
あると判定したときに、前記アンテナ制御部は衛星から
の受信信号により前記アンテナ手段が衛星の位置の捕捉
を行うように制御する。
【0025】本発明の衛星追尾方式によれば、前記分析
手段で求めた衛星の軌道関数の係数を評価手段で最適化
し、最適化された係数を有する軌道関数を使用して予測
追尾を行うことにより、衛星の予測追尾の誤差を収束さ
せ、追尾の精度を向上させることができる。
手段で求めた衛星の軌道関数の係数を評価手段で最適化
し、最適化された係数を有する軌道関数を使用して予測
追尾を行うことにより、衛星の予測追尾の誤差を収束さ
せ、追尾の精度を向上させることができる。
【0026】また、前記評価手段による衛星の軌道関数
の最適化にも種々の方法が考えられるが、通常は複数時
間分の衛星の軌道関数の各係数の平均値をとることによ
る。こうすることにより、最適化された係数を有する軌
道関数は、本来の静止軌道のものへと収束していき、衛
星の予測追尾の誤差が少なくなる。このように衛星の軌
道関数の各係数の平均値をとることによって最適化を行
うには、例えば、次のような2つの方法によることがで
きる。
の最適化にも種々の方法が考えられるが、通常は複数時
間分の衛星の軌道関数の各係数の平均値をとることによ
る。こうすることにより、最適化された係数を有する軌
道関数は、本来の静止軌道のものへと収束していき、衛
星の予測追尾の誤差が少なくなる。このように衛星の軌
道関数の各係数の平均値をとることによって最適化を行
うには、例えば、次のような2つの方法によることがで
きる。
【0027】(1)前記評価手段は、前記第2の記憶手
段に記憶された複数時間分の衛星の軌道関数の係数を平
均して、衛星の軌道関数の最適化を行う。
段に記憶された複数時間分の衛星の軌道関数の係数を平
均して、衛星の軌道関数の最適化を行う。
【0028】この場合、衛星の軌道関数の係数の数が前
記第2の記憶手段の容量を越えた場合には、前記分析手
段で新たに求められた軌道関数の係数を、最も古くに前
記第2の記憶手段に記憶された軌道関数の係数に上書き
していけばよい。
記第2の記憶手段の容量を越えた場合には、前記分析手
段で新たに求められた軌道関数の係数を、最も古くに前
記第2の記憶手段に記憶された軌道関数の係数に上書き
していけばよい。
【0029】(2)前記第2の記憶手段は、前記評価手
段で最適化された複数時間分の軌道関数の係数と前記分
析手段がそれまでに衛星の軌道関数を求めた回数を記憶
するものであり、前記評価手段は、前記分析手段が新た
に軌道関数を求める毎に、 {(記憶された最適化係数)×(記憶された分析回数)
+(新たに求めた軌道関数の係数)}÷(記憶された分
析回数+1) の計算を行うことにより軌道関数の最適化を行うと共
に、前記第2の記憶手段の記憶内容を新たに最適化され
た係数に更新し、分析回数に1を加える。
段で最適化された複数時間分の軌道関数の係数と前記分
析手段がそれまでに衛星の軌道関数を求めた回数を記憶
するものであり、前記評価手段は、前記分析手段が新た
に軌道関数を求める毎に、 {(記憶された最適化係数)×(記憶された分析回数)
+(新たに求めた軌道関数の係数)}÷(記憶された分
析回数+1) の計算を行うことにより軌道関数の最適化を行うと共
に、前記第2の記憶手段の記憶内容を新たに最適化され
た係数に更新し、分析回数に1を加える。
【0030】この場合、前記第2の記憶手段の記憶容量
は僅かで済む。さらに、前記分析手段で軌道関数を求め
た回数に関係なく、すべて同じ式、しかも簡単な計算で
軌道関数の係数を最適化することができる。もっとも、
衛星は時々軌道修正を行うため、一定期間で係数を最適
化した軌道関数も更新しなければならない。そこで、上
記の式で計算を行うのは一定回数のみに限り、以後は、 {(記憶された最適化係数)×(記憶された分析回数−
1)+(新たに求めた軌道関数の係数)}÷(記憶され
た分析回数) の計算で最適化を行えばよい。
は僅かで済む。さらに、前記分析手段で軌道関数を求め
た回数に関係なく、すべて同じ式、しかも簡単な計算で
軌道関数の係数を最適化することができる。もっとも、
衛星は時々軌道修正を行うため、一定期間で係数を最適
化した軌道関数も更新しなければならない。そこで、上
記の式で計算を行うのは一定回数のみに限り、以後は、 {(記憶された最適化係数)×(記憶された分析回数−
1)+(新たに求めた軌道関数の係数)}÷(記憶され
た分析回数) の計算で最適化を行えばよい。
【0031】
【発明の実施の形態】以下、本発明の一実施形態につい
て図面を参照して説明する。ここで説明する衛星は、静
止軌道に乗せられた静止衛星の例である。
て図面を参照して説明する。ここで説明する衛星は、静
止軌道に乗せられた静止衛星の例である。
【0032】図1は、本発明の一実施形態の衛星追尾方
式の機能ブロック図である。この衛星追尾方式は、衛星
軌跡情報記憶部11と、回帰分析処理部12と、回帰係
数記憶部13と、回帰係数評価部14と、衛星位置・指
令角計算処理部15と、アンテナ制御部16と、アンテ
ナ部17とから構成されている。
式の機能ブロック図である。この衛星追尾方式は、衛星
軌跡情報記憶部11と、回帰分析処理部12と、回帰係
数記憶部13と、回帰係数評価部14と、衛星位置・指
令角計算処理部15と、アンテナ制御部16と、アンテ
ナ部17とから構成されている。
【0033】衛星軌跡情報記憶部11は、ステップトラ
ック方式などによりアンテナ部17が捕捉した衛星の位
置を順次記憶するものである。回帰分析処理部12は、
衛星軌跡情報記憶部11に記憶された衛星軌跡情報から
最小二乗法による回帰分析により衛星の軌道関数を求
め、その軌道関数の係数(以下、「回帰係数」という)
を回帰係数評価部14へと送るものである。
ック方式などによりアンテナ部17が捕捉した衛星の位
置を順次記憶するものである。回帰分析処理部12は、
衛星軌跡情報記憶部11に記憶された衛星軌跡情報から
最小二乗法による回帰分析により衛星の軌道関数を求
め、その軌道関数の係数(以下、「回帰係数」という)
を回帰係数評価部14へと送るものである。
【0034】回帰係数記憶部13は、回帰分析処理部1
2で求めた回帰係数を順次記憶するものである。回帰係
数評価部14は、回帰係数記憶部13にすでに記憶され
ている過去の回帰係数および回帰分析処理部12で求め
た現在の回帰係数についての最適化を行って、適切な回
帰係数を決定するものである。
2で求めた回帰係数を順次記憶するものである。回帰係
数評価部14は、回帰係数記憶部13にすでに記憶され
ている過去の回帰係数および回帰分析処理部12で求め
た現在の回帰係数についての最適化を行って、適切な回
帰係数を決定するものである。
【0035】衛星位置・指令角計算処理部15は、回帰
係数評価部14で最適化された回帰係数を用いて衛星の
軌道関数を得て、この軌道関数から衛星の予測位置およ
びアンテナ部17への指令角を求めるものである。アン
テナ制御部16は、衛星位置・指令角計算処理部15で
求めた指令角にアンテナ部17を向けるために、また衛
星軌跡情報を得るために所定時間毎に衛星を捕捉するス
テップトラックを行うために、アンテナ部17の向きを
変えるように制御するものである。アンテナ部17は、
衛星から送られる信号を受信し、衛星の位置の捕捉を行
うものである。
係数評価部14で最適化された回帰係数を用いて衛星の
軌道関数を得て、この軌道関数から衛星の予測位置およ
びアンテナ部17への指令角を求めるものである。アン
テナ制御部16は、衛星位置・指令角計算処理部15で
求めた指令角にアンテナ部17を向けるために、また衛
星軌跡情報を得るために所定時間毎に衛星を捕捉するス
テップトラックを行うために、アンテナ部17の向きを
変えるように制御するものである。アンテナ部17は、
衛星から送られる信号を受信し、衛星の位置の捕捉を行
うものである。
【0036】次に、本発明の一実施形態の動作について
説明する。図2は、本発明の一実施形態の動作を示すフ
ローチャートである。
説明する。図2は、本発明の一実施形態の動作を示すフ
ローチャートである。
【0037】まず、不図示の判定手段が過去の衛星軌跡
情報が衛星軌跡情報記憶部11に少なくとも1周期分記
憶されているかどうかを判定する(S11)。この判定
は、記憶された衛星軌跡情報の数をカウントすることに
より行う。衛星軌跡情報が1周期分まで記憶されていな
ければ、アンテナ制御部16はステップトラック方式に
よりアンテナ部17を駆動して衛星の位置の捕捉を行い
(S18)、その位置を衛星軌跡情報として衛星軌跡情
報記憶部11に記憶する(M11)。
情報が衛星軌跡情報記憶部11に少なくとも1周期分記
憶されているかどうかを判定する(S11)。この判定
は、記憶された衛星軌跡情報の数をカウントすることに
より行う。衛星軌跡情報が1周期分まで記憶されていな
ければ、アンテナ制御部16はステップトラック方式に
よりアンテナ部17を駆動して衛星の位置の捕捉を行い
(S18)、その位置を衛星軌跡情報として衛星軌跡情
報記憶部11に記憶する(M11)。
【0038】一方、衛星軌跡情報が1周期分以上記憶さ
れているときは、衛星軌跡情報記憶部11から一定期間
分の衛星軌跡情報を読み出す(S12)。そして、回帰
分析処理部12は最小二乗法による線形回帰分析により
衛星の軌道関数を求める(S13)。
れているときは、衛星軌跡情報記憶部11から一定期間
分の衛星軌跡情報を読み出す(S12)。そして、回帰
分析処理部12は最小二乗法による線形回帰分析により
衛星の軌道関数を求める(S13)。
【0039】ここで、回帰分析処理部12が最小二乗法
による線形回帰分析で求める衛星の軌道関数には様々な
ものが使用できるが、その一例として数1をあげること
ができる。
による線形回帰分析で求める衛星の軌道関数には様々な
ものが使用できるが、その一例として数1をあげること
ができる。
【0040】
【数1】 y=A(1+Et)sin(t+D)+Bt+C
【0041】式中、yは衛星の位置、sin(t)は2
4時間を周期とする正弦関数、Aは正弦関数sin
(t)の振幅、Dは正弦関数sin(t)の位相、Eは
正弦関数sin(t)の振幅速度偏差、Bは速度偏差、
Cは位相偏差をそれぞれ示す。
4時間を周期とする正弦関数、Aは正弦関数sin
(t)の振幅、Dは正弦関数sin(t)の位相、Eは
正弦関数sin(t)の振幅速度偏差、Bは速度偏差、
Cは位相偏差をそれぞれ示す。
【0042】ここで求めた衛星の軌道関数について、そ
の回帰係数A,D,E,B,Cを回帰係数記憶部13に
設けられたテーブルに図3に示すような形で記憶する
(M12)。
の回帰係数A,D,E,B,Cを回帰係数記憶部13に
設けられたテーブルに図3に示すような形で記憶する
(M12)。
【0043】次に、回帰係数評価部14は、回帰係数記
憶部13に記憶された回帰係数から回帰係数の最適化を
行って、軌道関数の係数を決定する(S14)。ここで
は、回帰係数の最適化は回帰係数記憶部13に記憶され
た各々の回帰係数の平均値を取ることにより行う。
憶部13に記憶された回帰係数から回帰係数の最適化を
行って、軌道関数の係数を決定する(S14)。ここで
は、回帰係数の最適化は回帰係数記憶部13に記憶され
た各々の回帰係数の平均値を取ることにより行う。
【0044】ここで、時間t=t0〜tnについての回帰
係数が図3に示すように得られたとする。衛星軌跡情報
はn+1個あるから、例えば振幅Aを最適化したA’
は、数2の式で求められる。
係数が図3に示すように得られたとする。衛星軌跡情報
はn+1個あるから、例えば振幅Aを最適化したA’
は、数2の式で求められる。
【数2】A’=(A0+A1+A2+ …… +An-1+A
n)÷(n+1) 他の係数についても同様にして求められる。
n)÷(n+1) 他の係数についても同様にして求められる。
【0045】次に、新しい衛星軌跡情報を記憶する時刻
(例えば、1時間毎)となっているかどうかを判定する
(S15)。記憶時刻となっていなければ、衛星位置・
指令角計算処理部15は、ステップS14で求めた最適
化を行った回帰係数を有する軌道関数を使用して現在の
衛星位置を計算し、アンテナ部17に与える指令角を計
算する(S16)。ここで求めた指令角をアンテナ制御
部16に入力し、アンテナ制御部16からの指示により
アンテナ部17を駆動する(S17)。
(例えば、1時間毎)となっているかどうかを判定する
(S15)。記憶時刻となっていなければ、衛星位置・
指令角計算処理部15は、ステップS14で求めた最適
化を行った回帰係数を有する軌道関数を使用して現在の
衛星位置を計算し、アンテナ部17に与える指令角を計
算する(S16)。ここで求めた指令角をアンテナ制御
部16に入力し、アンテナ制御部16からの指示により
アンテナ部17を駆動する(S17)。
【0046】ステップS15の処理で新しい衛星軌跡情
報の記憶時刻となっていると判断したときは、ステップ
トラック方式によりアンテナ部17を駆動して衛星の位
置の捕捉を行い(S18)、その位置を新たな衛星軌跡
情報として衛星軌跡情報記憶部11に追加して記憶する
(M11)。
報の記憶時刻となっていると判断したときは、ステップ
トラック方式によりアンテナ部17を駆動して衛星の位
置の捕捉を行い(S18)、その位置を新たな衛星軌跡
情報として衛星軌跡情報記憶部11に追加して記憶する
(M11)。
【0047】以上の動作を繰り返し、常に衛星軌跡情報
から最小二乗法による回帰分析で求めた回帰係数を平均
化しながら使用することにより、衛星の予測追尾を行
う。
から最小二乗法による回帰分析で求めた回帰係数を平均
化しながら使用することにより、衛星の予測追尾を行
う。
【0048】本実施形態の衛星追尾方式では最小二乗法
による線形回帰分析で求めた衛星の軌道関数の係数を平
均化し、誤差を収束させて最適化することにより、衛星
の追尾の精度を向上させることができる。即ち、衛星の
軌道の誤差は南北および高低のいずれの場合にもほぼ同
じように生じるので、回帰係数の平均をとることによ
り、衛星の軌道関数は本来の軌道のものへと収束してい
き、予測追尾の誤差を小さくすることができる。
による線形回帰分析で求めた衛星の軌道関数の係数を平
均化し、誤差を収束させて最適化することにより、衛星
の追尾の精度を向上させることができる。即ち、衛星の
軌道の誤差は南北および高低のいずれの場合にもほぼ同
じように生じるので、回帰係数の平均をとることによ
り、衛星の軌道関数は本来の軌道のものへと収束してい
き、予測追尾の誤差を小さくすることができる。
【0049】なお、上記実施形態の衛星追尾方式におい
ては、過去の回帰係数を回帰係数記憶部13に記憶して
おき、回帰係数評価部14は過去の回帰係数を平均化す
ることにより回帰係数の最適化を行うように説明した。
しかし、回帰分析処理部12で求めた軌道関数の係数の
数が回帰係数記憶部13の容量を越えることがあり得
る。このような場合は、回帰分析処理部12で求めた最
新の軌道関数の回帰係数を回帰係数記憶部13に最も古
くに記憶された回帰係数に上書きして次々と更新してい
くようにすればよい。衛星軌跡情報の数が衛星軌跡情報
記憶部11の容量が越えた場合も同様にできる。
ては、過去の回帰係数を回帰係数記憶部13に記憶して
おき、回帰係数評価部14は過去の回帰係数を平均化す
ることにより回帰係数の最適化を行うように説明した。
しかし、回帰分析処理部12で求めた軌道関数の係数の
数が回帰係数記憶部13の容量を越えることがあり得
る。このような場合は、回帰分析処理部12で求めた最
新の軌道関数の回帰係数を回帰係数記憶部13に最も古
くに記憶された回帰係数に上書きして次々と更新してい
くようにすればよい。衛星軌跡情報の数が衛星軌跡情報
記憶部11の容量が越えた場合も同様にできる。
【0050】最適化の方法として他に、回帰係数記憶部
13に回帰分析処理部12で求めた回帰係数を複数記憶
するのではなく、回帰分析処理部12で軌道関数を求め
た回数と、回帰係数評価部14で求めた回帰係数の平均
値とを記憶し、回帰分析処理部で新たな軌道関数が求め
られる度に、 {(記憶平均値)×(記憶回数)+(新回帰係数)}÷
(記憶回数+1) の計算をして新たな平均値を求め、新たに回帰係数を最
適化して、衛星位置・指令角計算処理部15に送ると共
に、回帰係数記憶部13の内容を更新するようにしても
よい。
13に回帰分析処理部12で求めた回帰係数を複数記憶
するのではなく、回帰分析処理部12で軌道関数を求め
た回数と、回帰係数評価部14で求めた回帰係数の平均
値とを記憶し、回帰分析処理部で新たな軌道関数が求め
られる度に、 {(記憶平均値)×(記憶回数)+(新回帰係数)}÷
(記憶回数+1) の計算をして新たな平均値を求め、新たに回帰係数を最
適化して、衛星位置・指令角計算処理部15に送ると共
に、回帰係数記憶部13の内容を更新するようにしても
よい。
【0051】こうすることによって、回帰係数記憶部1
3の記憶容量が僅かで済む。さらに、求めた回帰係数の
数に関係なく、すべて同じ式、しかも簡単な計算で最適
化された回帰係数を求めることができる。
3の記憶容量が僅かで済む。さらに、求めた回帰係数の
数に関係なく、すべて同じ式、しかも簡単な計算で最適
化された回帰係数を求めることができる。
【0052】また、本発明で使用する線形回帰分析の式
は、上述したものに限られるものではなく、様々なもの
を使用することができる。衛星軌跡情報の蓄積された量
に従って、式を変更して使用することもできる。もっと
も、衛星の軌道関数を求める方法は最小二乗法による線
形回帰分析に限られるものではない。衛星の軌道関数の
最適化の方法も上述した方法だけに限られるものではな
い。
は、上述したものに限られるものではなく、様々なもの
を使用することができる。衛星軌跡情報の蓄積された量
に従って、式を変更して使用することもできる。もっと
も、衛星の軌道関数を求める方法は最小二乗法による線
形回帰分析に限られるものではない。衛星の軌道関数の
最適化の方法も上述した方法だけに限られるものではな
い。
【0053】なお、上記実施形態においては、静止軌道
上に乗せられた静止衛星の追尾を例として説明を行った
が、静止軌道以外の他の軌道に乗せられた人工衛星の追
尾にも本発明を適用することは可能である。
上に乗せられた静止衛星の追尾を例として説明を行った
が、静止軌道以外の他の軌道に乗せられた人工衛星の追
尾にも本発明を適用することは可能である。
【0054】
【発明の効果】以上説明したように、本発明の衛星追尾
方式によれば、衛星の軌道関数の係数の誤差を収束させ
て係数の最適化を行うことにより、衛星の予測追尾にお
ける誤差を収束させ、追尾精度を向上させることができ
る。
方式によれば、衛星の軌道関数の係数の誤差を収束させ
て係数の最適化を行うことにより、衛星の予測追尾にお
ける誤差を収束させ、追尾精度を向上させることができ
る。
【0055】さらに、請求項4または5の衛星追尾方式
によれば、衛星の軌道関数は本来の静止軌道のものへと
収束させることができ、追尾精度が非常に向上する。
によれば、衛星の軌道関数は本来の静止軌道のものへと
収束させることができ、追尾精度が非常に向上する。
【図1】本発明の一実施形態の衛星追尾方式の機能ブロ
ック図である。
ック図である。
【図2】本発明の一実施形態の衛星追尾方式のフローチ
ャートである。
ャートである。
【図3】本発明の一実施形態の回帰係数テーブルを示す
図である。
図である。
【図4】従来例の衛星追尾方式の機能ブロック図であ
る。
る。
【図5】従来例の衛星追尾方式のフローチャートであ
る。
る。
11 衛星軌跡情報記憶部 12 回帰分析処理部 13 回帰係数記憶部 14 回帰係数評価部 15 衛星位置・指令角計算処理部 16 アンテナ制御部 17 アンテナ部
Claims (5)
- 【請求項1】 衛星から送られる信号の受信および衛星
の位置の捕捉を行うアンテナ手段と、このアンテナ手段
で捕捉した衛星の位置を衛星軌跡情報として蓄積する第
1の記憶手段と、この第1の記憶手段に蓄積された衛星
軌跡情報を用いて所定時間毎に衛星の軌道関数を求める
分析手段と、この分析手段で求めた所定時間毎の衛星の
軌道関数の係数を複数時間分記憶する第2の記憶手段
と、前記第2の記憶手段に記憶された複数時間分の軌道
関数の係数を最適化して衛星の予測追尾に使用する軌道
関数の係数を決定する評価手段と、この評価手段で決定
した係数を有する軌道関数を用いて衛星の方向を計算す
る計算手段と、この計算手段で計算された衛星の方向に
前記アンテナ手段を向けるように制御し、また、一定時
間毎に前記アンテナ手段が衛星の位置の捕捉を行うよう
に制御するアンテナ制御部とを備えた衛星追尾方式。 - 【請求項2】 前記第1の記憶手段に蓄積された衛星軌
跡情報の情報量が所定時間分以上であるかを判定する判
定手段を備え、この判定手段で衛星の軌跡情報が所定時
間分以上であると判定したときに、前記分析手段は衛星
の軌道関数を求めることを特徴とする請求項1に記載の
衛星追尾方式。 - 【請求項3】 前記判定手段で衛星の軌跡情報が所定時
間分未満であると判定したときに、前記アンテナ制御部
は衛星からの受信信号により前記アンテナ手段が衛星の
位置の捕捉を行うように制御することを特徴とする請求
項2に記載の衛星追尾方式。 - 【請求項4】 前記評価手段は、前記第2の記憶手段に
記憶された複数時間分の衛星の軌道関数の係数を平均し
て、衛星の軌道関数の最適化を行うことを特徴とする請
求項1ないし3のいずれかに記載の衛星追尾方式。 - 【請求項5】 前記第2の記憶手段は、前記評価手段で
最適化された複数時間分の軌道関数の係数と前記分析手
段がそれまでに衛星の軌道関数を求めた回数を記憶する
ものであり、前記評価手段は、前記分析手段が新たに軌
道関数を求める毎に、 {(記憶された最適化係数)×(記憶された分析回数)
+(新たに求めた軌道関数の係数)}÷(記憶された分
析回数+1) の計算を行うことにより軌道関数の最適化を行うと共
に、前記第2の記憶手段の記憶内容を新たに最適化され
た係数に更新し、分析回数に1を加えることを特徴とす
る請求項1ないし3のいずれかに記載の衛星追尾方式。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP16921696A JPH1020010A (ja) | 1996-06-28 | 1996-06-28 | 衛星追尾方式 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP16921696A JPH1020010A (ja) | 1996-06-28 | 1996-06-28 | 衛星追尾方式 |
Publications (1)
Publication Number | Publication Date |
---|---|
JPH1020010A true JPH1020010A (ja) | 1998-01-23 |
Family
ID=15882375
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP16921696A Pending JPH1020010A (ja) | 1996-06-28 | 1996-06-28 | 衛星追尾方式 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JPH1020010A (ja) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH11326477A (ja) * | 1998-05-12 | 1999-11-26 | Mitsubishi Electric Corp | アンテナ制御装置及びアンテナ制御方法 |
ES2170725A1 (es) * | 2000-06-30 | 2002-08-01 | Mitsubishi Electric Corp | Metodo para adquirir un satelite. |
KR100432312B1 (ko) * | 2001-12-06 | 2004-05-22 | 한국전자통신연구원 | 위성의 궤도 예측에 의한 자동화된 영상촬영 계획 수립장치 및 그 방법 |
KR100662535B1 (ko) | 2004-09-24 | 2006-12-28 | 세이코 엡슨 가부시키가이샤 | 측위 시스템, 측위장치, 측위정보 제공장치, 측위방법,측위 시스템의 제어 프로그램을 기록한 컴퓨터읽기 가능한기록매체 |
JP2009220622A (ja) * | 2008-03-13 | 2009-10-01 | Fujitsu Ltd | 中高度人工衛星捕捉方法及び装置 |
JP2010173033A (ja) * | 2009-01-30 | 2010-08-12 | Sodick Co Ltd | 工作機械における変位補正式の生成装置および工作機械における変位補正式の生成方法 |
JP2018059844A (ja) * | 2016-10-06 | 2018-04-12 | 富士通株式会社 | 観測計画立案プログラム、観測計画立案装置、及び観測計画立案方法 |
-
1996
- 1996-06-28 JP JP16921696A patent/JPH1020010A/ja active Pending
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH11326477A (ja) * | 1998-05-12 | 1999-11-26 | Mitsubishi Electric Corp | アンテナ制御装置及びアンテナ制御方法 |
ES2170725A1 (es) * | 2000-06-30 | 2002-08-01 | Mitsubishi Electric Corp | Metodo para adquirir un satelite. |
KR100432312B1 (ko) * | 2001-12-06 | 2004-05-22 | 한국전자통신연구원 | 위성의 궤도 예측에 의한 자동화된 영상촬영 계획 수립장치 및 그 방법 |
KR100662535B1 (ko) | 2004-09-24 | 2006-12-28 | 세이코 엡슨 가부시키가이샤 | 측위 시스템, 측위장치, 측위정보 제공장치, 측위방법,측위 시스템의 제어 프로그램을 기록한 컴퓨터읽기 가능한기록매체 |
JP2009220622A (ja) * | 2008-03-13 | 2009-10-01 | Fujitsu Ltd | 中高度人工衛星捕捉方法及び装置 |
JP2010173033A (ja) * | 2009-01-30 | 2010-08-12 | Sodick Co Ltd | 工作機械における変位補正式の生成装置および工作機械における変位補正式の生成方法 |
JP2018059844A (ja) * | 2016-10-06 | 2018-04-12 | 富士通株式会社 | 観測計画立案プログラム、観測計画立案装置、及び観測計画立案方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US5862495A (en) | Real time position correction to ground generated spacecraft ephemeris | |
US7579984B2 (en) | Ultra-tightly coupled GPS and inertial navigation system for agile platforms | |
US8838374B2 (en) | Automatic correction of past position errors for location and inventory tracking | |
US5422829A (en) | Closed-loop control for scanning application | |
US7242345B2 (en) | Method for controlling a radar antenna | |
CN109639341B (zh) | 无星历自主接入低轨广播卫星的方法 | |
JPH1020010A (ja) | 衛星追尾方式 | |
JPH04326079A (ja) | Gps受信機 | |
JPH1120800A (ja) | 衛星捕捉・追尾装置 | |
JP7364857B2 (ja) | 位置検知システム、センサ端末および位置検知方法 | |
CN116576868A (zh) | 一种多传感器融合精确定位及自主导航方法 | |
JP4322829B2 (ja) | サイクルスリップ検出装置及びサイクルスリップ検出方法 | |
JP3427325B2 (ja) | Gps受信装置 | |
EP3631515A1 (en) | Method for estimating a position of a mobile device using gnss signals | |
JPH06281717A (ja) | 衛星追尾方式 | |
JP4400330B2 (ja) | 位置検出装置及び位置検出方法 | |
JPH1114729A (ja) | 追尾型アンテナ装置 | |
US8571506B2 (en) | Systems and methods for sharing an oscillator between receivers | |
JP2001237629A (ja) | 衛星追尾アンテナ駆動制御方法、および衛星追尾アンテナ駆動制御装置 | |
JPH08271622A (ja) | 目標追尾装置 | |
JPH0464074A (ja) | 追尾型アンテナ装置 | |
JPS6210903A (ja) | ステツプトラツク装置 | |
JP2950248B2 (ja) | アンテナ追尾制御方法とその装置 | |
US20240159912A1 (en) | Method and apparatus for determining a frequency related parameter of a frequency source using predictive control | |
JP3492403B2 (ja) | 自動レーダ追尾装置 |