(57)【特許請求の範囲】
【請求項１】
衛星が受信する受信観測信号であって基準信号が含まれる受信観測信号を取得し、前記基準信号を用いて前記受信観測信号に対して遅延検波を実行することにより、前記衛星の移動により生じる前記受信観測信号のドップラー量を位相差で表す第１ドップラー量を算出し、算出した前記第１ドップラー量と前記基準信号とに基づいて相関を実行することにより前記第１ドップラー量の確からしさを第１相関値として算出する遅延検波部と、
前記遅延検波部から取得した前記第１ドップラー量と前記基準信号を用いて前記受信観測信号に対して同期検波を実行することにより前記受信観測信号のドップラー量を周波数で表す第２ドップラー量を算出し、算出した前記第２ドップラー量と前記基準信号とに基づいて相関を実行することにより前記第２ドップラー量の確からしさを第２相関値として算出する同期検波部と、
前記第１相関値と前記第２相関値とに基づいて、前記第１ドップラー量と前記第２ドップラー量とのいずれを前記受信観測信号の信号ドップラー量とするかを判定する検波出力判定部と
を備えるドップラー推定装置。
【請求項２】
前記同期検波部は、
前記遅延検波部により算出された前記第１ドップラー量を第１ドップラー周波数として周波数に変換し、前記第１ドップラー周波数に基づいて第１レプリカ信号を生成し、生成した第１レプリカ信号を用いて前記同期検波を実行することにより前記第２ドップラー量を算出する請求項１に記載のドップラー推定装置。
【請求項３】
前記ドップラー推定装置は、
前記信号ドップラー量を判定するための相関閾値を記憶する相関閾値記憶部を備え、
前記検波出力判定部は、
前記第１相関値と前記第２相関値との各々を前記相関閾値と比較すると共に、前記第１相関値と前記第２相関値とを比較し、前記相関閾値より大きく、かつ、他方の相関値よりも大きい相関値を特定し、特定した相関値に対応するドップラー量を前記信号ドップラー量とする請求項１または２に記載のドップラー推定装置。
【請求項４】
前記基準信号は、前記受信観測信号を検出するためのトレーニングシーケンスである請求項１から３のいずれか１項に記載のドップラー推定装置。
【請求項５】
前記遅延検波部は、前記トレーニングシーケンスを用いて前記第１相関値を算出し、
前記同期検波部は、前記トレーニングシーケンスを用いて前記第２相関値を算出する請求項４に記載のドップラー推定装置。
【請求項６】
前記ドップラー推定装置は、
前記検波出力判定部により判定された前記信号ドップラー量に基づいて、前記受信観測信号を補正する信号補正部を備える請求項１から５のいずれか１項に記載のドップラー推定装置。
【請求項７】
衛星が受信する受信観測信号であって基準信号が含まれる受信観測信号を取得し、前記基準信号を用いて前記受信観測信号に対して遅延検波を実行することにより、前記衛星の移動により生じる前記受信観測信号のドップラー量を位相差で表す第１ドップラー量を算出し、算出した前記第１ドップラー量と前記基準信号とに基づいて相関を実行することにより前記第１ドップラー量の確からしさを第１相関値として算出する遅延検波処理と、
前記遅延検波処理から取得した前記第１ドップラー量と前記基準信号を用いて前記受信観測信号に対して同期検波を実行することにより前記受信観測信号のドップラー量を周波数で表す第２ドップラー量を算出し、算出した前記第２ドップラー量と前記基準信号とに基づいて相関を実行することにより前記第２ドップラー量の確からしさを第２相関値として算出する同期検波処理と、
前記第１相関値と前記第２相関値とに基づいて、前記第１ドップラー量と前記第２ドップラー量とのいずれを前記受信観測信号の信号ドップラー量とするかを判定する検波出力判定処理とをコンピュータであるドップラー推定装置に実行させるプログラム。
【請求項８】
前記同期検波処理は、
前記遅延検波処理により算出された前記第１ドップラー量を第１ドップラー周波数として周波数に変換し、変換した前記第１ドップラー周波数に基づいて第１レプリカ信号を生成し、生成した第１レプリカ信号を用いて前記同期検波を実行することにより前記第２ドップラー量を算出する請求項７に記載のプログラム。
【請求項９】
前記ドップラー推定装置は、前記信号ドップラー量を判定するための相関閾値を記憶する相関閾値記憶部を備え、
前記検波出力判定処理は、
前記第１相関値と前記第２相関値との各々を前記相関閾値と比較すると共に、前記第１相関値と前記第２相関値とを比較し、前記相関閾値より大きく、かつ、他方の相関値よりも大きい相関値を特定し、特定した相関値に対応するドップラー量を前記信号ドップラー量とする請求項７または８に記載のプログラム。
【請求項１０】
前記基準信号は、前記受信観測信号を検出するためのトレーニングシーケンスである請求項７から９のいずれか１項に記載のプログラム。
【請求項１１】
前記遅延検波処理は、前記トレーニングシーケンスを用いて前記第１相関値を算出し、
前記同期検波処理は、前記トレーニングシーケンスを用いて前記第２相関値を算出する請求項１０に記載のプログラム。
【請求項１２】
前記プログラムは、さらに、
前記検波出力判定処理により判定された前記信号ドップラー量に基づいて、前記受信観測信号を補正する信号補正処理を備える請求項７から１１のいずれか１項に記載のプログラム。
【請求項１３】
遅延検波部が、衛星が受信する受信観測信号であって基準信号が含まれる受信観測信号を取得し、前記基準信号を用いて前記受信観測信号に対して遅延検波を実行することにより、前記衛星の移動により生じる前記受信観測信号のドップラー量を位相差で表す第１ドップラー量を算出し、算出した前記第１ドップラー量と前記基準信号とに基づいて相関を実行することにより前記第１ドップラー量の確からしさを第１相関値として算出し、
同期検波部が、前記遅延検波部から取得した前記第１ドップラー量と前記基準信号を用いて前記受信観測信号に対して同期検波を実行することにより前記受信観測信号のドップラー量を周波数で表す第２ドップラー量を算出し、算出した前記第２ドップラー量と前記基準信号とに基づいて相関を実行することにより前記第２ドップラー量の確からしさを第２相関値として算出し、
検波出力判定部が、前記第１相関値と前記第２相関値とに基づいて、前記第１ドップラー量と前記第２ドップラー量とのいずれを前記受信観測信号の信号ドップラー量とするかを判定するドップラー推定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ドップラー推定装置、プログラム及びドップラー推定方法に関する。特に、衛星で受信する自動移動体識別システム（ＡＩＳ：Ａｕｔｏｍａｔｉｃ  Ｉｄｅｎｔｉｆｉｃａｔｉｏｎ  Ｓｙｓｔｅｍ）におけるＡＩＳ信号についてドップラーを推定するドップラー推定装置、プログラム及びドップラー推定方法に関する。
【背景技術】
【０００２】
現在、移動体及び港湾施設の保安のため、移動体を自動識別する自動移動体識別システム、すなわちＡＩＳが利用されている。ＡＩＳは、船名、船種、位置、針路、速度、目的地、積載物等のＡＩＳ信号を周辺移動体や陸上局に向けて自動的に送信するシステムである。ＡＩＳでは、海上約５０Ｋｍ半径内におけるＴＤＭＡ（Ｔｉｍｅ  Ｄｉｖｉｓｉｏｎ
Ｍｕｌｔｉｐｌｅ  Ａｃｃｅｓｓ：時分割多元接続）により、移動体相互または移動体と陸上局との通信を行っている。
【０００３】
近年、このＡＩＳを衛星に搭載することによりＡＩＳ信号を衛星で受信し、全地球における移動体の航行情報を得るシステムが開発されている。
しかし、衛星による宇宙空間における受信は、地上に比べて格段に広範囲、具体的には半径約２５００Ｋｍ程度での受信となるため、ＴＤＭＡ管理外に存在する他移動体のＡＩＳ信号が入感し、信号同士の衝突が発生し、信号解読が困難となる。さらに、海上運用では影響が無かったドップラー周波数の影響が発生し、ＡＩＳ信号の解読をより困難にしている。
【０００４】
特許文献１には、船舶の信号を、例えば衛星で観測と同時に航空機でも観測するなどにより、違う状況で同じ信号を同時観測してその信号タイミングまで既知の情報を用いて衝突分離する技術が開示されている（特許文献１参照）。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２０１２−０２９２８５号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
特許文献１では、衛星受信信号からドップラー周波数の影響を効果的に排除することができないという課題がある。
【０００７】
本発明は、衛星受信信号からドップラー周波数の影響を効果的に排除することを目的とする。
【課題を解決するための手段】
【０００８】
本発明に係るドップラー推定装置は、
衛星が受信する受信観測信号であって基準信号が含まれる受信観測信号を取得し、前記基準信号を用いて前記受信観測信号に対して遅延検波を実行することにより、前記衛星の移動により生じる前記受信観測信号のドップラー量を位相差で表す第１ドップラー量を算出し、算出した前記第１ドップラー量と前記基準信号とに基づいて相関を実行することにより前記第１ドップラー量の確からしさを第１相関値として算出する遅延検波部と、
前記遅延検波部から取得した前記第１ドップラー量と前記基準信号を用いて前記受信観測信号に対して同期検波を実行することにより前記受信観測信号のドップラー量を周波数で表す第２ドップラー量を算出し、算出した前記第２ドップラー量と前記基準信号とに基づいて相関を実行することにより前記第２ドップラー量の確からしさを第２相関値として出力する同期検波部と、
前記遅延検波部から出力された前記第１相関値と前記第２相関値とに基づいて、前記第１ドップラー量と前記第２ドップラー量とのいずれを前記受信観測信号の信号ドップラー量とするかを判定する検波出力判定部とを備える。
【０００９】
前記同期検波部は、
前記遅延検波部により算出された前記第１ドップラー量を第１ドップラー周波数として周波数に変換し、前記第１ドップラー周波数に基づいて第１レプリカ信号を生成し、生成した第１レプリカ信号を用いて前記同期検波を実行することにより前記第２ドップラー量を算出する。
【００１０】
前記ドップラー推定装置は、
前記信号ドップラー量を判定するための相関閾値を記憶する相関閾値記憶部を備え、
前記検波出力判定部は、
前記第１相関値と前記第２相関値との各々を前記相関閾値と比較すると共に、前記第１相関値と前記第２相関値とを比較し、前記相関閾値より大きく、かつ、他方の相関値よりも大きい相関値を特定し、特定した相関値に対応するドップラー量を前記信号ドップラー量とする。
【００１１】
前記基準信号は、前記受信観測信号を検出するためのトレーニングシーケンスである。
【００１２】
前記遅延検波部は、前記トレーニングシーケンスを用いて前記第１相関値を算出し、
前記同期検波部は、前記トレーニングシーケンスを用いて前記第２相関値を算出する。
【００１３】
前記ドップラー推定装置は、
前記検波出力判定部により判定された前記信号ドップラー量に基づいて、前記受信観測信号を補正する信号補正部を備える。
【００１４】
本発明に係るプログラムは、
衛星が受信する受信観測信号であって基準信号が含まれる受信観測信号を取得し、前記基準信号を用いて前記受信観測信号に対して遅延検波を実行することにより、前記衛星の移動により生じる前記受信観測信号のドップラー量を位相差で表す第１ドップラー量を算出し、算出した前記第１ドップラー量と前記基準信号とに基づいて相関を実行することにより前記第１ドップラー量の確からしさを第１相関値として算出する遅延検波処理と、
前記遅延検波処理から取得した前記第１ドップラー量と前記基準信号を用いて前記受信観測信号に対して同期検波を実行することにより前記受信観測信号のドップラー量を周波数で表す第２ドップラー量を算出し、算出した前記第２ドップラー量と前記基準信号とに基づいて相関を実行することにより前記第２ドップラー量の確からしさを第２相関値として出力する同期検波処理と、
前記遅延検波処理から出力された前記第１相関値と前記第２相関値とに基づいて、前記第１ドップラー量と前記第２ドップラー量とのいずれを前記受信観測信号の信号ドップラー量とするかを判定する検波出力判定処理とをコンピュータであるドップラー推定装置に実行させる。
【００１５】
前記同期検波処理は、
前記遅延検波処理により算出された前記第１ドップラー量を第１ドップラー周波数として周波数に変換し、変換した前記第１ドップラー周波数に基づいて第１レプリカ信号を生成し、生成した第１レプリカ信号を用いて前記同期検波を実行することにより前記第２ドップラー量を算出する。
【００１６】
前記ドップラー推定装置は、前記信号ドップラー量を判定するための相関閾値を記憶する相関閾値記憶部を備え、
前記検波出力判定処理は、
前記第１相関値と前記第２相関値との各々を前記相関閾値と比較すると共に、前記第１相関値と前記第２相関値とを比較し、前記相関閾値より大きく、かつ、他方の相関値よりも大きい相関値を特定し、特定した相関値に対応するドップラー量を前記信号ドップラー量とする。
【００１７】
前記基準信号は、前記受信観測信号を検出するためのトレーニングシーケンスである。
【００１８】
前記遅延検波処理は、前記トレーニングシーケンスを用いて前記第１相関値を算出し、
前記同期検波処理は、前記トレーニングシーケンスを用いて前記第２相関値を算出する。
【００１９】
前記プログラムは、さらに、
前記検波出力判定処理により判定された前記信号ドップラー量に基づいて、前記受信観測信号を補正する信号補正処理を備える。
【００２０】
本発明に係るドップラー推定方法は、
遅延検波部が、衛星が受信する受信観測信号であって基準信号が含まれる受信観測信号を取得し、前記基準信号を用いて前記受信観測信号に対して遅延検波を実行することにより、前記衛星の移動により生じる前記受信観測信号のドップラー量を位相差で表す第１ドップラー量を算出し、算出した前記第１ドップラー量と前記基準信号とに基づいて相関を実行することにより前記第１ドップラー量の確からしさを第１相関値として算出し、
同期検波部が、前記遅延検波部から取得した前記第１ドップラー量と前記基準信号を用いて前記受信観測信号に対して同期検波を実行することにより前記受信観測信号のドップラー量を周波数で表す第２ドップラー量を算出し、算出した前記第２ドップラー量と前記基準信号とに基づいて相関を実行することにより前記第２ドップラー量の確からしさを第２相関値として出力し、
検波出力判定部が、前記遅延検波部から出力された前記第１相関値と前記第２相関値とに基づいて、前記第１ドップラー量と前記第２ドップラー量とのいずれを前記受信観測信号の信号ドップラー量とするかを判定する。
【発明の効果】
【００２１】
本発明に係るドップラー推定装置では、遅延検波部が、基準信号を用いて受信観測信号に対して遅延検波を実行することにより位相差である第１ドップラー量を算出する。さらに、遅延検波部は、基準信号と第１ドップラー量とに基づいて相関を実行することにより第１ドップラー量の確からしさを第１相関値として算出する。また、同期検波部が、第１ドップラー量を用いて受信観測信号に対して同期検波を実行することにより周波数である第２ドップラー量を算出する。さらに、同期検波部は、第２ドップラー量と基準信号に基づいて相関を実行することにより第２ドップラー量の確からしさを第２相関値として算出する。また、検波出力判定部が、第１相関値と第２相関値とに基づいて、第１ドップラー量と第２ドップラー量とのいずれを受信観測信号の信号ドップラー量とするかを判定する。よって、本発明に係るドップラー推定装置によれば、第１ドップラー量を用いて同期検波部の初期ドップラー推定値を設定することにより同期時間を短縮しつつ、ドップラー推定帯域を十分狭められる。これにより、衛星がＡＩＳ信号を受ける場合に発生するドップラーの同期検波による推定を短時間かつ精度の良いドップラー推定が可能となり、衛星受信信号からドップラー周波数の影響を効率的に排除することができる。
【図面の簡単な説明】
【００２２】
【図１】地上及び海上におけるＡＩＳの概要について説明する図。
【図２】衛星がＡＩＳ信号を受信する場合の信号衝突を説明する図。
【図３】実施の形態１に係る信号検波システムのブロック構成図。
【図４】実施の形態１に係る衛星模擬部の詳細ブロック構成図。
【図５】実施の形態１に係る模擬信号生成部の詳細ブロック構成図。
【図６】実施の形態１に係る信号検波システムのハードウェア構成図。
【図７】実施の形態１に係る信号検波システムの信号検波方法の動作を示すフロー図。
【図８】実施の形態１に係る検波時刻決定処理における検波計画部の動作を示すフロー図。
【図９】実施の形態１に係る模擬信号生成処理における模擬信号生成部の動作を示すフロー図。
【図１０】実施の形態１に係る信号処理における信号処理部の動作を示すフロー図。
【図１１】実施の形態１に係る信号分離処理を示す図。
【図１２】実施の形態１に係る信号相関処理の動作を示すフロー図。
【図１３】実施の形態２に係るドップラー推定装置のブロック構成図。
【図１４】実施の形態２に係るドップラー推定装置におけるドップラー推定方法の動作を示すフロー図。
【図１５】観測計画の具体例（ａ）から（ｅ）についての説明図。
【発明を実施するための形態】
【００２３】
実施の形態１．
図１を用いて、地上及び海上における船舶などの移動体を識別する自動船舶識別システム、すなわちＡＩＳの概要について説明する。
以下の実施の形態では、移動体として船舶を想定している。
図１に示すように、地上及び海上におけるＡＩＳにおいて、移動体６００同士のＡＩＳ信号の送受信は、半径約５０Ｋｍ内のＳＯＴＤＭＡ（Ｓｅｌｆ−Ｏｒｇａｎｉｚｅｄ・Ｔｉｍｅ・Ｄｉｖｉｓｉｏｎ・Ｍｕｌｔｉｐｌｅ・Ａｃｃｅｓｓ）により行われる。半径約５０Ｋｍ内のＳＯＴＤＭＡにより、ＡＩＳ信号の衝突は回避される。また、船舶である移動体６００の速度は高速艇でも約３０ノット程度であり、ドップラー周波数は約１０Ｈｚ程度となり、ほぼ無いに等しい。
【００２４】
図２を用いて、衛星７００がＡＩＳ信号を受信する場合の信号衝突について説明する。
図２に示すように、衛星７００が観測可能な観測エリア７０１は、半径約２５００Ｋｍに及ぶ。
ここで、衛星７００が受信する信号を衛星受信信号４００とし、観測対象の観測移動体６０１から送信されるＡＩＳ信号を観測信号４０１とし、観測対象外の観測外移動体６０２から送信されるＡＩＳ信号を観測外信号４０２とする。
【００２５】
衛星７００には、ＳＯＴＤＭＡで制御できる範囲外のエリアからもＡＩＳ信号が入射される。衛星７００は、観測信号４０１に観測外信号４０２が衝突した信号を衛星受信信号４００として受信する。ＳＯＴＤＭＡ範囲外の信号衝突は無作為であり、このままでは復調あるいは復号は困難である。
また、衛星７００の速度は約８０００ｍ／ｓｅｃであるため、±４ＫＨｚ程度のドップラー周波数が発生する。このドップラー周波数は、データレートである約１０ＫＨｚの幅の範囲内となってしまい、フィルタリングで分離することが困難である。
【００２６】
＊＊＊構成の説明＊＊＊
図３を用いて、本実施の形態に係る信号検波システム５００のブロック構成について説明する。
信号検波システム５００は、衛星情報１１０、衛星模擬部１２０、検波計画部１３０、模擬信号生成部１４０、信号模擬情報１５０、信号処理部１６０を備える。さらに、信号検波システム５００は、解析結果ＤＢ１７０、衛星受信信号記憶部１８１、処理情報１８２、データ加工部１９１、マニュアル入力部１９２を備える。
【００２７】
信号検波システム５００は、衛星７００が受信する衛星受信信号４００から、観測対象の観測移動体６０１により送信される観測信号４０１を検波する。観測信号４０１は観測移動体６０１のＡＩＳ信号である。
【００２８】
衛星模擬部１２０は、観測信号４０１を検波する検波時刻４２３における衛星７００の状態を模擬し、模擬した衛星７００の状態を検波時衛星模擬情報４２０として出力する。
衛星情報１１０は、衛星７００に関する情報を記憶する。
衛星模擬部１２０は、衛星情報１１０を用いて、検波時刻４２３における衛星７００の状態を模擬する。衛星模擬部１２０は、ＡＩＳ受信機を搭載する衛星７００の位置、軌道、及び姿勢を模擬する。すなわち、衛星模擬部１２０は、検波時衛星模擬情報４２０として、検波時刻４２３における衛星７００の位置と速度と姿勢とを算出する。
【００２９】
検波計画部１３０は、衛星模擬部１２０により出力された検波時衛星模擬情報４２０に基づいて、検波時刻４２３において衛星７００が観測する移動体であって観測移動体以外の移動体である観測外移動体６０２の情報を観測外移動体情報４２５として取得する。具体的には、検波計画部１３０は、検波時衛星模擬情報４２０に基づいて衛星７００が観測する観測エリア７０１を算出し、算出した観測エリア７０１に存在する観測外移動体６０２の情報を観測外移動体情報４２５として取得する。すなわち、検波計画部１３０は、移動体の過去の移動体情報に基づいて、次パス観測で衝突することが想定される観測外移動体６０２を同定する。また、検波計画部１３０は、衝突データの状況を推定する。
なお、検波計画部１３０は、過去の移動体情報が蓄積された解析結果ＤＢ１７０から、検波時刻４２３において衛星７００が観測するはずの観測外移動体６０２に関する移動体情報を抽出する。よって、実際の検波時刻４２３に観測する観測移動体６０１の情報が観測外移動体６０２の情報として抽出されることはないものとする。移動体情報とはＡＩＳ情報のことである。また、過去の移動体情報とは、言い換えると既知の移動体情報である。
【００３０】
解析結果ＤＢ１７０には、過去に衛星７００が受信した衛星受信信号４００を信号検波した結果、抽出された観測外信号４０２が示す移動体情報が格納されている。すなわち、解析結果ＤＢ１７０には、衛星７００が過去に受信した観測移動体６０１の移動体情報が蓄積されている。
検波計画部１３０は、検波時衛星模擬情報４２０に基づいて、解析結果ＤＢ１７０から観測エリア７０１に存在する予定の観測外移動体６０２の観測外移動体情報４２５を取得する。
観測外移動体６０２は、複数存在するものとする。
【００３１】
模擬信号生成部１４０は、検波計画部１３０により取得された観測外移動体情報４２５に基づいて、検波時刻４２３に観測外移動体６０２から送信される信号を模擬信号４２６として生成する。
信号模擬情報１５０には、信号を模擬する際に用いられる情報が格納されている。
模擬信号生成部１４０は、観測外移動体情報４２５と信号模擬情報１５０とを用いて、検波時刻４２３に観測外移動体６０２から送信される信号を模擬信号４２６として生成する。
【００３２】
信号処理部１６０は、信号分離部１６１、復調部１６４を備える。復調部１６４は、ドップラー推定部１６２、信号復号部１６３を備える。
【００３３】
信号分離部１６１は、模擬信号生成部１４０により生成された模擬信号４２６を衛星受信信号４００から分離し、模擬信号４２６が分離された衛星受信信号４００を受信観測信号４２７として出力する。
【００３４】
信号検波システム５００は、陸上局に設置される。あるいは、信号検波システム５００は、衛星７００に搭載されていてもよい。
衛星受信信号記憶部１８１は、衛星７００により受信された衛星受信信号４００を、衛星７００から受信し記憶機器に記憶する。
【００３５】
復調部１６４は、信号分離部１６１により出力された受信観測信号４２７を復調する。
ドップラー推定部１６２は、観測移動体６０１におけるドップラー周波数を観測移動体ドップラー周波数４７１として算出する。ドップラー推定部１６２は、算出した観測移動体ドップラー周波数４７１に基づいて、信号分離部１６１から出力された受信観測信号４２７を補正し、補正した受信観測信号４２７を補正受信観測信号４２８として出力する。
ここで、補正受信観測信号４２８は、すなわち、観測移動体６０１から送信されたＡＩＳ信号である観測信号４０１である。
【００３６】
信号復号部１６３は、ドップラー推定部１６２により出力された補正受信観測信号４２８を復号し、復号移動体情報４３１として出力する。信号復号部１６３から出力された復号移動体情報４３１は、解析結果ＤＢ１７０に格納される。復号移動体情報４３１は、観測移動体６０１から送信された観測信号４０１が示すＡＩＳ情報、すなわち観測移動体情報である。
【００３７】
データ加工部１９１は、解析結果ＤＢ１７０に格納された移動体情報を用いて、移動体の運航状況をグローバルに把握するため等の各種のデータ加工を実行する。
マニュアル入力部１９２は、過去の移動体情報を解析結果ＤＢ１７０にマニュアルで入力させる。
【００３８】
図４を用いて、衛星模擬部１２０の詳細構成について説明する。
衛星模擬部１２０は、地球重力モデル部１２１、アンテナ姿勢モデル部１２２、衛星軌道計算部１２３を備える。
また、衛星情報１１０には、衛星テレメータ情報１１１、軌道決定値情報１１２が記憶される。
【００３９】
アンテナ姿勢モデル部１２２は、検波計画部１３０から取得した検波時刻４２３と衛星テレメータ情報１１１とより、検波時刻４２３における衛星７００のアンテナの方向を取得する。アンテナ姿勢モデル部１２２は、アンテナ姿勢モデルをアンテナ姿勢モデルテーブルとして備え、検波時刻４２３とアンテナ姿勢モデルテーブルと衛星テレメータ情報１１１とより、検波時刻４２３における衛星７００のアンテナの方向を取得する。
【００４０】
地球重力モデル部１２１は、衛星７００の軌道伝搬に必要な重力モデルを取得する。地球重力モデル部１２１は、地球重力モデルを地球重力モデルテーブルとして備え、地球重力モデルテーブルを用いて必要な重力モデルを取得する。
【００４１】
衛星軌道計算部１２３は、地球重力モデル、アンテナ姿勢モデル、軌道決定値情報１１２から取得した最新の軌道決定値、及び検波計画部１３０から取得した検波時刻４２３を取得する。衛星軌道計算部１２３は、地球重力モデル、アンテナ姿勢モデル、最新の軌道決定値、及び検波時刻４２３を入力として、衛星７００において観測移動体６０１を観測する検波時刻４２３における衛星７００の位置、速度、及び姿勢を計算する。衛星軌道計算部１２３は、算出した検波時刻４２３における衛星７００の位置、速度、及び姿勢を検波時衛星模擬情報４２０として出力する。
【００４２】
ここで、検波計画部１３０が検波時刻４２３を決定する機能について説明する。
まず、衛星模擬部１２０は、衛星受信信号４００から観測信号４０１を検波する条件を満たす検波期間４２４を取得し、検波期間４２４における衛星７００の状態を検波期間衛星状態４２１として模擬する。衛星模擬部１２０は、検波期間衛星状態４２１を検波計画部１３０に出力する。
【００４３】
次に、検波計画部１３０は、衛星模擬部１２０により模擬された検波期間衛星状態４２１に基づいて、検波期間４２４において衛星７００が観測する移動体であって観測移動体６０１以外の移動体である観測外移動体６０２の情報を検波期間観測外移動体情報４３０として取得する。検波計画部１３０は、取得した検波期間観測外移動体情報４３０に基づいて、検波時刻４２３を決定する。具体的には、検波計画部１３０は、検波期間観測外移動体情報４３０を用いて検波期間４２４における観測外移動体６０２の位置を判定し、検波期間４２４における観測外移動体６０２の位置に基づいて検波時刻４２３を決定する。
【００４４】
図５を用いて、模擬信号生成部１４０の詳細構成について説明する。
模擬信号生成部１４０は、複数信号生成部１４９、信号加算部１４５、ＧＭＳＫ復調部１４６を備える。
複数信号生成部１４９は、移動体信号模擬部１４１、ＧＭＳＫ変調部１４２、ドップラー周波数挿入部１４３、時間遅延補正部１４４を備える。
【００４５】
複数信号生成部１４９は、観測外移動体情報４２５に基づいて、検波時刻４２３における複数の観測外移動体６０２０の各々の情報を複数の移動体送信信号４２９０として生成する。複数の観測外移動体６０２０の各々を観測外移動体６０２と表記する場合がある。複数の移動体送信信号４２９０の各々を移動体送信信号４２９と表記する場合がある。
【００４６】
移動体信号模擬部１４１は、観測外移動体情報４２５に基づいて、検波時刻４２３における複数の観測外移動体６０２の各々の情報を複数の移動体信号１４１１として模擬する。
ＧＭＳＫ変調部１４２は、複数の移動体信号１４１１の各々をガウス最小偏移変調（ＧＭＳＫ：Ｇａｕｓｓｉａｎ・ｍｉｎｉｍｕｍ・ｓｈｉｆｔ・ｋｅｙｉｎｇ）する。
ドップラー周波数挿入部１４３は、移動体信号模擬部１４１により生成された複数の移動体信号１４１１の各々について、観測外移動体情報４２５に基づいてドップラー周波数を算出する。ドップラー周波数挿入部１４３は、算出した複数の移動体信号１４１１の各々のドップラー周波数に基づいて複数の移動体信号１４１１の各々を補正する。
時間遅延補正部１４４は、ドップラー周波数挿入部１４３により補正された複数の移動体信号１４１１の各々について、観測外移動体情報４２５に基づいて時間遅延を算出する。時間遅延補正部１４４は、算出した複数の移動体信号１４１１の各々の時間遅延に基づいて複数の移動体信号１４１１の各々を補正し、補正した複数の移動体信号１４１１の各々を複数の移動体送信信号４２９の各々として出力する。複数の移動体信号１４１１０の各々を移動体信号１４１１と表記する場合がある。
【００４７】
信号加算部１４５は、複数信号生成部１４９により生成された複数の移動体送信信号４２９の各々を加算し、加算した信号を模擬信号４２６として生成する。
ＧＭＳＫ復調部１４６は、信号加算部１４５から出力された信号をＧＭＳＫ復調し、模擬信号４２６として出力する。
【００４８】
なお、上記において説明した移動体情報６１０、観測外移動体情報４２５、検波期間観測外移動体情報４３０は、船舶である移動体のＡＩＳ情報である。
【００４９】
図６を用いて、本実施の形態に係る信号検波システム５００、後述するドップラー推定装置１６２０のハードウェア構成例について説明する。
信号検波システム５００及びドップラー推定装置１６２０は１つの装置から構成される信号検波装置であるものとして説明する。
【００５０】
信号検波システム５００、ドップラー推定装置１６２０はコンピュータである。
信号検波システム５００、ドップラー推定装置１６２０は、プロセッサ９０１、補助記憶装置９０２、メモリ９０３、通信装置９０４、入力インタフェース９０５、ディスプレイインタフェース９０６といったハードウェアを備える。
プロセッサ９０１は、信号線９１０を介して他のハードウェアと接続され、これら他のハードウェアを制御する。
入力インタフェース９０５は、入力装置９０７に接続されている。
ディスプレイインタフェース９０６は、ディスプレイ９０８に接続されている。
【００５１】
プロセッサ９０１は、プロセッシングを行うＩＣ（Ｉｎｔｅｇｒａｔｅｄ  Ｃｉｒｃｕｉｔ）である。
プロセッサ９０１は、例えば、ＣＰＵ（Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）、ＤＳＰ（Ｄｉｇｉｔａｌ  Ｓｉｇｎａｌ  Ｐｒｏｃｅｓｓｏｒ）、ＧＰＵ（Ｇｒａｐｈｉｃｓ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）である。
補助記憶装置９０２は、例えば、ＲＯＭ（Ｒｅａｄ  Ｏｎｌｙ  Ｍｅｍｏｒｙ）、フラッシュメモリ、ＨＤＤ（Ｈａｒｄ  Ｄｉｓｋ  Ｄｒｉｖｅ）である。
メモリ９０３は、例えば、ＲＡＭ（Ｒａｎｄｏｍ  Ａｃｃｅｓｓ  Ｍｅｍｏｒｙ）である。
通信装置９０４は、データを受信するレシーバー９４０１及びデータを送信するトランスミッター９４０２を含む。
通信装置９０４は、例えば、通信チップ又はＮＩＣ（Ｎｅｔｗｏｒｋ  Ｉｎｔｅｒｆａｃｅ  Ｃａｒｄ）である。
入力インタフェース９０５は、入力装置９０７のケーブル９１１が接続されるポートである。
入力インタフェース９０５は、例えば、ＵＳＢ（Ｕｎｉｖｅｒｓａｌ  Ｓｅｒｉａｌ  Ｂｕｓ）端子である。
ディスプレイインタフェース９０６は、ディスプレイ９０８のケーブル９１２が接続されるポートである。
ディスプレイインタフェース９０６は、例えば、ＵＳＢ端子又はＨＤＭＩ（登録商標）（Ｈｉｇｈ  Ｄｅｆｉｎｉｔｉｏｎ  Ｍｕｌｔｉｍｅｄｉａ  Ｉｎｔｅｒｆａｃｅ）端子である。
入力装置９０７は、例えば、マウス、キーボード又はタッチパネルである。  ディスプレイ９０８は、例えば、ＬＣＤ（Ｌｉｑｕｉｄ  Ｃｒｙｓｔａｌ  Ｄｉｓｐｌａｙ）である。
【００５２】
補助記憶装置９０２には、衛星模擬部１２０、検波計画部１３０、模擬信号生成部１４０、信号処理部１６０、遅延検波部１６２１、同期検波部１６２２、検波出力判定部１６２３、信号補正部１６２４等の機能を実現するプログラムが記憶されている。以下、衛星模擬部１２０、検波計画部１３０、模擬信号生成部１４０、信号処理部１６０、遅延検波部１６２１、同期検波部１６２２、検波出力判定部１６２３、信号補正部１６２４等をまとめて「部」と表記する。
このプログラムは、メモリ９０３にロードされ、プロセッサ９０１に読み込まれ、プロセッサ９０１によって実行される。
更に、補助記憶装置９０２には、ＯＳ（Ｏｐｅｒａｔｉｎｇ  Ｓｙｓｔｅｍ）も記憶されている。
そして、ＯＳの少なくとも一部がメモリ９０３にロードされ、プロセッサ９０１はＯＳを実行しながら、「部」の機能を実現するプログラムを実行する。
図６では、１つのプロセッサ９０１が図示されているが、信号検波システム５００、ドップラー推定装置１６２０が複数のプロセッサ９０１を備えていてもよい。
そして、複数のプロセッサ９０１が「部」の機能を実現するプログラムを連携して実行してもよい。
また、「部」の処理の結果を示す情報やデータや信号値や変数値が、メモリ９０３、補助記憶装置９０２、又は、プロセッサ９０１内のレジスタ又はキャッシュメモリにファイルとして記憶される。
【００５３】
「部」を「サーキットリー」で提供してもよい。
また、「部」を「回路」又は「工程」又は「手順」又は「処理」に読み替えてもよい。
「回路」及び「サーキットリー」は、プロセッサ９０１だけでなく、ロジックＩＣ又はＧＡ（Ｇａｔｅ  Ａｒｒａｙ）又はＡＳＩＣ（Ａｐｐｌｉｃａｔｉｏｎ  Ｓｐｅｃｉｆｉｃ  Ｉｎｔｅｇｒａｔｅｄ  Ｃｉｒｃｕｉｔ）又はＦＰＧＡ（Ｆｉｅｌｄ−Ｐｒｏｇｒａｍｍａｂｌｅ  Ｇａｔｅ  Ａｒｒａｙ）といった他の種類の処理回路をも包含する概念である。
【００５４】
また、プログラムプロダクトと称されるものは、「部」として説明している機能を実現するプログラムが記録された記憶媒体、記憶装置などであり、見た目の形式に関わらず、コンピュータ読み取り可能なプログラムをロードしているものである。
【００５５】
＊＊＊動作の説明＊＊＊
図７を用いて、本実施の形態に係る信号検波システム５００の信号検波方法について説明する。
信号検波システム５００は、プロセッサ９０１、上記ハードウェア資源と協働し、以下に示す処理を行うプログラム９９０を実行することにより信号検波方法、信号検波処理を実現する。
信号検波システム５００の信号検波方法は、衛星７００が受信する衛星受信信号４００から観測対象の観測移動体６０１により送信される観測信号４０１を検波する。
【００５６】
＜検波時刻決定処理Ｓ１００＞
Ｓ１００において、衛星模擬部１２０及び検波計画部１３０は、検波時刻決定処理Ｓ１００を実行する。
衛星模擬部１２０は、衛星受信信号４００から観測信号４０１を検波する条件を満たす検波期間４２４を取得し、検波期間４２４における衛星の状態を検波期間衛星状態４２１として模擬する。衛星模擬部１２０は、例えば、衛星受信信号４００から観測信号４０１を検波するのに適した条件を、条件テーブル１１３として備える。検波するのに適した条件とは、例えば、「衛星７００の観測エリア７０１に移動体が２隻存在する」等の条件である。衛星模擬部１２０は、この条件テーブル１１３と衛星テレメータ情報１１１とに基づいて、検波するのに適した検波期間４２４を取得し、検波期間４２４における衛星７００の状態を模擬する。検波期間４２４は、例えば、「開始時刻Ｔ１から９０分間」等の情報である。衛星情報１１０が条件テーブル１１３を備えていてもよい。
そして、検波計画部１３０は、模擬された検波期間衛星状態４２１に基づいて、検波期間において衛星７００が観測する移動体であって観測移動体６０１以外の移動体である観測外移動体６０２の情報を検波期間観測外移動体情報４３０として取得する。検波計画部１３０は、取得した検波期間観測外移動体情報４３０に基づいて、検波時刻４２３を決定する。検波計画部１３０は、検波期間４２４が「開始時刻Ｔ１から９０分間」である場合、検波期間４２４のうち最も検波に適した時刻を検波時刻４２３として決定する。このように、検波時刻決定処理Ｓ１００では、検波期間観測外移動体情報４３０を用いて検波期間４２４における観測外移動体６０２の位置を判定し、検波期間４２４における観測外移動体６０２の位置に基づいて検波時刻４２３を決定する。
【００５７】
図８を用いて、検波時刻決定処理Ｓ１００における検波計画部１３０の動作についてさらに説明する。
Ｓ１２１において、検波計画部１３０は、観測エリア計算処理を実行する。観測エリア計算処理では、検波計画部１３０は衛星模擬部１２０から取得した検波期間衛星状態４２１である衛星７００の位置、速度、及び姿勢に基づき、地表上の観測エリア７０１を計算する。
Ｓ１２２において、検波計画部１３０は、移動体抽出処理を実行する。移動体抽出処理では、検波計画部１３０は、観測エリア計算処理で求めた観測エリア７０１と、解析結果ＤＢ１７０に収められた既知の移動体情報とにより、検波期間４２４のうちのある時刻に観測すれば検出される予定の移動体、すなわち観測外移動体６０２の移動体情報を抽出する。
Ｓ１２３において、アンテナ指向方向計画作成処理を実行する。アンテナ指向方向計画作成処理では、検波計画部１３０は、抽出された移動体情報から、ある時刻におけるアンテナの指向方向を計算し、観測計画を生成する。
観測計画とは信号を観測するにあたって、その要求されるパス、パス内の時刻幅、アンテナを指向させる方向、アンテナの観測モード（静止や回転、場合によってはアンテナビーム幅）を規定した衛星における観測スケジュールを示す情報である。
【００５８】
以下に、図１５の（ａ）から図１５の（ｅ）を用いて、観測計画の具体例について簡単に説明する。
図１５の（ａ）は、観測計画として、衛星に対して特定の方向（例えば衛星直下点方向に対して進行方向にｘ°の方向）を観測する計画を示している。
図１５の（ｂ）は、上記のアンテナ指向方向を初期値として直下点廻りにレートωで回転させる観測計画を示している。
図１５の（ｃ）は、アンテナ動作角度範囲内で既知のＡＩＳ船舶を指向する観測計画を示している。
図１５の（ｄ）は、アンテナ動作角度範囲内で前回のパスで観測したエリアを含む方向を指向する観測計画（観測場所の自動計算によるアンテナ指向方法計画設定）を示している。
図１５の（ｅ）は、既に抽出されている船舶を衛星直下点（ドップラー０）やアンテナ覆域の縁（ドップラー最大）に持つアンテナ指向方向を指向する観測計画を示している。
このような設定を行うことにより、あるパスにおいてある時刻にてどちらの方向にアンテナを向けた観測を行うかを観測計画として作成する。
【００５９】
Ｓ１２４において、情報利用部位指定処理を実行する。情報利用部位指定処理では、検波計画部１３０は、移動体抽出処理による移動体情報の抽出結果に基づき、移動体情報のうち利用が可能なビット範囲、ＡＩＳパラメータなどが含まれる部位指定情報を決定する。ＡＩＳパラメータとは、例えば、ＡＩＳ信号を検出するためのトレーニングシーケンスである。
Ｓ１２５において、衛星軌道判定処理を実行する。衛星軌道判定処理では、検波計画部１３０は、アンテナ指向方向計画作成処理で求められた観測計画、情報利用部位指定処理で求め部位指定情報、解析結果ＤＢ１７０に収められた観測時刻確定判定情報に基づき、検波時刻４２３を決定する。または、検波計画部１３０は、検波時刻４２３の更新を衛星模擬部１２０に要求として出力する。
【００６０】
なお、図８の説明では、検波計画部１３０が、観測エリア計算処理、移動体抽出処理、アンテナ指向方向計画作成処理、情報利用部位指定処理、衛星軌道判定処理を実行するものとした。しかし、上記の各処理を実行する観測エリア計算部、移動体抽出部、アンテナ指向方向計画作成部、情報利用部位指定部、衛星軌道判定部を、検波計画部１３０が備える構成でもよい。
【００６１】
＜衛星模擬処理Ｓ１１０＞
衛星模擬処理Ｓ１１０において、衛星模擬部１２０は、観測信号４０１を検波する検波時刻４２３における衛星７００の状態を模擬し、模擬した衛星７００の状態を検波時衛星模擬情報４２０として出力する衛星模擬処理Ｓ１１０を実行する。衛星模擬部１２０は、衛星模擬処理Ｓ１１０において、検波時衛星模擬情報４２０として、検波時刻４２３における衛星７００の位置と速度と姿勢とを算出する。
【００６２】
図４を用いて、衛星模擬処理Ｓ１１０の動作についてさらに説明する。
アンテナ姿勢モデル部１２２は、検波計画部１３０から取得した検波時刻４２３と衛星テレメータ情報１１１とより、検波時刻４２３における衛星７００のアンテナの方向を取得する。アンテナ姿勢モデル部１２２は、アンテナ姿勢モデルをアンテナ姿勢モデルテーブルとして備え、検波時刻４２３とアンテナ姿勢モデルテーブルと衛星テレメータ情報１１１とより、検波時刻４２３における衛星７００のアンテナの方向を取得する。
【００６３】
地球重力モデル部１２１は、衛星７００の軌道伝搬に必要な重力モデルを取得する。地球重力モデル部１２１は、地球重力モデルを地球重力モデルテーブルとして備え、地球重力モデルテーブルを用いて必要な重力モデルを取得する。
【００６４】
衛星軌道計算部１２３には、地球重力モデル部１２１から地球重力モデル、アンテナ姿勢モデル部１２２からアンテナ姿勢モデル、軌道決定値情報１１２から最新の軌道決定値、検波計画部１３０から検波時刻４２３が入力される。衛星軌道計算部１２３は、地球重力モデル、アンテナ姿勢モデル、軌道決定値、及び検波時刻４２３を入力として、衛星７００において観測移動体６０１を観測する検波時刻４２３における衛星７００の位置、速度、及び姿勢を計算する。衛星軌道計算部１２３は、算出した検波時刻４２３における衛星７００の位置、速度、及び姿勢を検波時衛星模擬情報４２０として出力する。
【００６５】
＜検波計画処理Ｓ１２０＞
検波計画処理Ｓ１２０において、検波計画部１３０は、衛星模擬処理Ｓ１１０により出力された検波時衛星模擬情報４２０に基づいて、検波時刻４２３において衛星７００が観測する移動体であって観測移動体６０１以外の移動体である観測外移動体６０２の情報を観測外移動体情報４２５として取得する検波計画処理Ｓ１２０を実行する。検波計画部１３０は、検波時衛星模擬情報４２０に基づいて衛星７００が観測する観測エリア７０１を算出し、算出した観測エリア７０１に存在する観測外移動体６０２の情報を観測外移動体情報４２５として取得する。
【００６６】
図８を用いて、検波計画処理Ｓ１２０における検波計画部１３０の動作について説明する。
Ｓ１２１の観測エリア計算処理では、検波計画部１３０は衛星模擬部１２０から取得した検波時衛星模擬情報４２０である検波時刻４２３における衛星７００の位置、速度、及び姿勢に基づき、検波時刻４２３の観測エリア７０１を計算する。
Ｓ１２２の移動体抽出処理では、検波計画部１３０は、検波時刻４２３の観測エリア７０１と、解析結果ＤＢ１７０に収められた移動体情報６１０とにより、検波時刻４２３に観測すれば検出される予定の観測外移動体６０２の観測外移動体情報４２５を抽出する。観測外移動体情報４２５は、検波時刻４２３に観測すれば検出される予定の移動体であって観測移動体６０１以外の観測外移動体６０２の移動体情報である。
検波計画部１３０は、検波時刻４２３の観測エリア７０１と、解析結果ＤＢ１７０に収められた移動体情報６１０とにより、検波時刻４２３に観測すれば検出される予定の移動体であって観測移動体６０１以外の観測外移動体６０２の移動体情報を観測外移動体情報４２５として抽出する。
Ｓ１２４の情報利用部位指定処理では、検波計画部１３０は、抽出された観測外移動体情報４２５に基づき、観測外移動体情報４２５のうち利用が可能なビット範囲、ＡＩＳパラメータなどが含まれる観測外移動体情報４２５の部位指定情報を決定する。検波計画部１３０は、観測外移動体情報４２５と部位指定情報とを信号処理部１６０に出力する。
【００６７】
＜模擬信号生成処理Ｓ１３０＞
図５及び図９を用いて、模擬信号生成処理Ｓ１３０における模擬信号生成部１４０の動作について説明する。
模擬信号生成処理Ｓ１３０において、模擬信号生成部１４０は、検波計画処理Ｓ１２０により取得された観測外移動体情報４２５に基づいて、検波時刻４２３に観測外移動体６０２から送信される信号を模擬信号４２６として生成する模擬信号生成処理Ｓ１３０を実行する。観測外移動体６０２は、複数存在する。
【００６８】
模擬信号生成処理Ｓ１３０は、複数信号生成処理Ｓ１３０ａと、信号加算処理Ｓ１３５と、ＧＭＳＫ復調処理Ｓ１３６とを備える。
複数信号生成処理Ｓ１３０ａでは、複数信号生成部１４９は、観測外移動体情報４２５に基づいて、検波時刻４２３における複数の観測外移動体６０２の各々の情報を複数の移動体送信信号４２９として生成する複数信号生成処理Ｓ１３０ａを実行する。
複数信号生成処理Ｓ１３０ａは、移動体信号模擬処理Ｓ１３１、ＧＭＳＫ変調処理Ｓ１３２、ドップラー周波数挿入処理Ｓ１３３、時間遅延補正処理Ｓ１３４を備える。
【００６９】
移動体信号模擬処理Ｓ１３１では、移動体信号模擬部１４１は、観測外移動体情報４２５に基づいて、検波時刻４２３における複数の観測外移動体６０２の各々の情報を複数の移動体信号として模擬する移動体信号模擬処理Ｓ１３１を実行する。移動体信号模擬部１４１は、既知の観測外移動体情報４２５をＡＩＳビットストリームに変換する。
【００７０】
ＧＭＳＫ変調処理Ｓ１３２では、ＧＭＳＫ変調部１４２は、移動体信号模擬処理Ｓ１３１により模擬された複数の移動体信号をＧＭＳＫ変調する。ＧＭＳＫ変調部１４２は、ＡＩＳ信号の規格に基づいて、移動体信号模擬部１４１から出力された複数の移動体信号をＲＦまたはＩＦ信号に変調する。
【００７１】
ドップラー周波数挿入処理Ｓ１３３では、ドップラー周波数挿入部１４３は、移動体信号模擬処理Ｓ１３１により模擬されＧＭＳＫ変調処理Ｓ１３２により変調された複数の移動体信号の各々について観測外移動体情報４２５に基づいてドップラー周波数を算出する。ドップラー周波数挿入部１４３は、算出した複数の移動体信号の各々のドップラー周波数に基づいて複数の移動体信号の各々を補正するドップラー周波数挿入処理Ｓ１３３を実行する。ドップラー周波数挿入部１４３は、衛星模擬部１２０からの検波時衛星模擬情報４２０に含まれる位置及び速度情報と、解析結果ＤＢ１７０からの移動体の位置及び速度情報、すなわち観測外移動体情報４２５とに基づいて、ドップラー周波数を計算する。ドップラー周波数挿入部１４３は、ＲＦまたはＩＦ信号に変調された複数の移動体信号の各々に対し、計算されたドップラー周波数分だけ周波数をシフトする。
【００７２】
時間遅延補正処理Ｓ１３４では、時間遅延補正部１４４は、ドップラー周波数挿入処理Ｓ１３３により補正された複数の移動体信号の各々について、観測外移動体情報４２５に基づいて時間遅延を算出する。時間遅延補正部１４４は、算出した複数の移動体信号の各々の時間遅延に基づいて複数の移動体信号の各々を補正し、補正した複数の移動体信号の各々を複数の移動体送信信号４２９の各々として出力する時間遅延補正処理Ｓ１３４を実行する。時間遅延補正部１４４は、ドップラー周波数挿入部１４３から出力された複数の移動体信号の各々について、以下の時間遅延補正処理Ｓ１３４を実行する。時間遅延補正部１４４は、衛星模擬部１２０からの検波時衛星模擬情報４２０に含まれる位置情報と、解析結果ＤＢ１７０からの移動体の位置情報と、信号模擬情報１５０に含まれる衛星搭載アンテナモデル１５１とに基づいて、生成信号の振幅と時間遅延とを計算し、移動体信号に付加する。
【００７３】
なお、複数信号生成部１４９が、地上で受信した移動体信号である地上受信ＡＩＳ信号から、衛星で受信する移動体信号である衛星受信ＡＩＳ信号を模擬する場合は、ドップラー周波数挿入部１４３には解析結果ＤＢ１７０からＡＩＳ実信号ファイルが入力される。ＡＩＳ実信号ファイルとは、ＡＩＳにより実際に受信された移動体信号である。
【００７４】
信号加算処理Ｓ１３５では、信号加算部１４５は、複数信号生成処理Ｓ１３０ａにより生成された複数の移動体送信信号４２９の各々を加算し、加算した信号を出力する。信号加算部１４５は、複数信号生成処理Ｓ１３０ａが施された想定衝突信号数分の移動体信号を加算し、加算信号として出力する。
【００７５】
ＧＭＳＫ復調処理Ｓ１３６では、ＧＭＳＫ復調部１４６は、信号加算処理Ｓ１３５から出力された加算信号をＧＭＳＫ復調し、模擬信号４２６として出力する。ＧＭＳＫ復調部１４６は、加算信号をＧＭＳＫ復調し、適当なＩＦ帯信号に変換する。なお、信号加算部１４５が、複数信号生成処理Ｓ１３０ａにより生成された複数の移動体送信信号４２９の各々を加算し、加算した信号をＧＭＳＫ復調し、模擬信号４２６として生成してもよい。
【００７６】
図１０を用いて、信号処理部１６０による信号処理Ｓ１４０の動作について説明する。
信号処理Ｓ１４０は、信号分離処理Ｓ１５０と、復調処理Ｓ１６０を備える。
信号分離処理Ｓ１５０には、分離対象となる衛星受信信号４００と、模擬信号生成処理Ｓ１３０により生成された模擬信号４２６とが入力される。
【００７７】
＜信号分離処理Ｓ１５０＞
信号分離処理Ｓ１５０において、信号分離部１６１は、模擬信号生成処理Ｓ１３０により生成された模擬信号４２６を衛星受信信号４００から分離する。信号分離部１６１は、模擬信号４２６が分離された衛星受信信号４００を受信観測信号４２７として出力する信号分離処理Ｓ１５０を実行する。
ここで、信号処理部１６０は、必要があれば衛星受信信号４００を複素信号化する複素信号変換処理を実行する。衛星受信信号４００は複素信号であることを前提にして処理を考えているが、衛星コンフィギュレーションによっては複素信号でなくリアル信号でインターフェースされる場合もある。このためリアル信号がインターフェースされた場合には入力した時点で複素信号への変換を行う。信号を複素信号の形で処理すると言うことは信号振幅と信号位相を別に表現できるため各種の処理（例えばドップラー位相抽出等）が簡便に処理できる。
【００７８】
図１１及び図１２は、信号分離処理Ｓ１５０を説明するための図である。
図１１に示すように、信号分離処理Ｓ１５０では信号相関処理Ｓ１５１が実行される。
信号相関処理Ｓ１５１において、信号分離部１６１は、模擬信号生成部１４０から出力された模擬信号４２６と、分離対象である衛星受信信号４００と、検波計画部１３０から出力される観測外移動体情報４２５の部位指定情報とを取得する。信号分離部１６１は、取得した模擬信号４２６と、衛星受信信号４００と、観測外移動体情報４２５の部位指定情報とを用いて相関処理を実行し、衝突信号タイミングを同定する。
ここで、模擬信号４２６には、検波時刻４２３における衛星７００の位置、速度及び姿勢に基づくアンテナゲイン、既知移動体方向のドップラー情報、電力減衰量が含まれる。また、分離対象である衛星受信信号４００は、衝突データを含む信号である。また、検波計画部１３０から出力される観測外移動体情報４２５の部位指定情報とは、衝突データの状況推定結果である。
【００７９】
図１２を用いて、信号相関処理Ｓ１５１の動作について説明する。
Ｓ１５１１の相関実行処理において、信号分離部１６１は、衝突を含む観測されたＡＩＳ信号である衛星受信信号４００から、既知のデータに基づく観測時点のＡＩＳ衝突予測信号である模擬信号４２６を差引く。模擬信号４２６を差し引いた残りのデータを予測信号とする。信号分離部１６１は、差し引いた残りのデータについて共通データとの相関を取る。Ｓ１５１１における相関において利用可能な共通データとは、以下のようなデータである。検波計画部１３０にて推定されたトレーニングシーケンス、スタートフラグ、エンドフラグ、移動体緯度・経度等相関に利用可能な既知ＡＩＳデータの一部または全部である。すなわち、共通データとは、検波計画部１３０において推定した衝突データの状況結果から求められる相関検出の元となるデータである。
【００８０】
Ｓ１５１２において、信号分離部１６１は、印加時刻シフトにより予測信号の印加タイミングをずらす。
Ｓ１５１３において、信号分離部１６１は、予測信号と共通データとの相関から信号印加時刻を決定すると共に、信号印加時刻に合わせた衝突データ印加タイミングを決定し、衝突データを相殺した真のＡＩＳデータを抽出し、抽出した真のＡＩＳデータを受信観測信号４２７として出力する。
【００８１】
以上のように、信号分離部１６１は、衛星受信信号４００から模擬信号４２６を分離した受信観測信号４２７を出力する。受信観測信号４２７は、衛星７００が受信した衛星受信信号４００から衝突信号を分離した信号であり、衝突がなかったと仮定した場合に、衛星７００が観測移動体６０１から受信する移動体信号に相当する。
【００８２】
＜復調処理Ｓ１６０＞
復調処理Ｓ１６０において、復調部１６４は、信号分離処理Ｓ１５０において信号分離部１６１により出力された受信観測信号４２７を復調する復調処理Ｓ１６０を実行する。
復調処理Ｓ１６０は、ドップラー推定処理Ｓ１６１と、信号復号処理Ｓ１６２とを備える。
【００８３】
ドップラー推定処理Ｓ１６１において、ドップラー推定部１６２は、観測移動体６０１におけるドップラー周波数を観測移動体ドップラー周波数４７１として算出する。ドップラー推定部１６２は、算出した観測移動体ドップラー周波数４７１に基づいて、信号分離処理Ｓ１５０から出力された受信観測信号４２７を補正し、補正した受信観測信号４２７を補正受信観測信号４２８として出力するドップラー推定処理Ｓ１６１を実行する。
ドップラー推定部１６２は、信号分離部１６１から出力された受信観測信号４２７に基づいて、観測移動体ドップラー周波数４７１を推定する。ドップラー推定部１６２は、例えば、同期検波方式を用いて、受信観測信号４２７の観測移動体ドップラー周波数４７１を推定する。
【００８４】
信号復号処理Ｓ１６２において、信号復号部１６３は、ドップラー推定処理Ｓ１６１により出力された補正受信観測信号４２８を復号し、復号移動体情報４３１として出力する信号復号処理Ｓ１６２を実行する。
【００８５】
信号復号処理Ｓ１６２は、シンボルタイミング検知処理Ｓ１６２ａ、ＣＲＣ／ＳＦ／ＥＦ／ビットスタッフィング判定処理Ｓ１６２ｂ、工学値変換処理Ｓ１６２ｃを備える。
シンボルタイミング検知処理Ｓ１６２ａ及びＣＲＣ／ＳＦ／ＥＦ／ビットスタッフィング判定処理Ｓ１６２ｂにおいて、信号復号部１６３は、復号を行うためトレーニングシーケンスを検知する。ビットスタッフィング処理とはＡＩＳの０／１のデータ部にＳＦやＥＦが偶然入り込まないようにＳＦやＥＦと同じビットパターンになったデータを強制的にある手順に基づいて変更するもので、復号において逆の手順で正しいビットに直す処理である。
工学値変換処理Ｓ１６２ｃは、工学値変換を実行することにより、補正受信観測信号４２８を意味ある物理量である復号移動体情報４３１に変換する。工学値変換処理Ｓ１６２ｃは、例えば、ＭＳＧ＿ＩＤ工学値変換を実行する。信号復号部１６３は、復号移動体情報４３１を解析結果ＤＢ１７０に格納する。
信号復号処理Ｓ１６２により出力される復号移動体情報４３１は、検波時刻４２３に観測された観測移動体６０１についての移動体情報、すなわちＡＩＳ情報である。復号移動体情報４３１は、観測移動体６０１から衛星に送信された観測信号４０１が示す移動体情報である。
【００８６】
以上で、信号検波システム５００の信号検波方法についての説明を終わる。
なお、信号検波システム５００は、陸上局に設置された１つの装置であってもよい。あるいは、信号検波システム５００は、複数の装置により構成されたシステムであってもよい。
【００８７】
＊＊＊効果の説明＊＊＊
本実施の形態によれば、移動体や基地局からのＡＩＳ信号を衛星で受信する場合、復調に大きな影響を与える信号同士の衝突を推定し、正常な検波を行うことができる。地上の運用範囲、半径約５０Ｋｍ以外の場所からの衝突信号を被衝突信号至近からと仮定すれば信号強度はほぼ同じと考えられ、ドップラーも同程度と考えられる。この場合、被衝突信号は衝突信号によって万遍なく干渉を受けることになり信号復調／復号はほぼ不可能と考えられる。しかし、本実施の形態によれば、既知のデータに基づいて衝突信号を模擬することができるので、衛星が受信した衛星受信信号から衝突信号を分離することができる。
【００８８】
また、本実施の形態によれば、衝突信号のドップラーを模擬するだけでなく、観測信号のドップラーをも補正することができるので、高い精度で観測信号を検波することができる。
【００８９】
また、本実施の形態によれば、一旦分離した信号は解析結果ＤＢに登録することによって既知信号となり、これを繰り返し利用することによって分離信号、すなわち受信観測信号の数を増加していくことができる。このプロシジャを繰り返すことにより定常状態では地球上のすべての移動体を同定し、且つ新たに衝突した信号に対してより確実な分離が可能となる。
【００９０】
実施の形態２．
本実施の形態では、主に、実施の形態１と異なる点について説明する。
本実施の形態において、実施の形態１で説明した構成と同様の構成には同一の符号を付し、その説明を省略する。
【００９１】
本実施の形態では、信号処理部１６０のドップラー推定部１６２により高精度にドップラーを推定することができるドップラー推定方法について説明する。
本実施の形態では、実施の形態１で説明したドップラー推定部１６２を１つの装置であるドップラー推定装置１６２０として説明する。ドップラー推定部１６２は、ドップラー推定装置１６２０でもよい。以下、ドップラー推定装置１６２０として説明する。
【００９２】
＊＊＊構成の説明＊＊＊
図１３を用いて、ドップラー推定装置１６２０のブロック構成について説明する。
ドップラー推定装置１６２０は、遅延検波部１６２１、同期検波部１６２２、検波出力判定部１６２３、信号補正部１６２４を備える。また、遅延検波部１６２１は、フィルタリング部１６２１１、最適タイミング取得部１６２１２、ドップラー推定補償部１６２１３を備える。また、同期検波部１６２２は、フィルタリング部１６２２１、搬送波捕捉追尾部１６２２２を備える。
【００９３】
遅延検波部１６２１は、衛星７００が受信する受信観測信号４２７であって基準信号４２７１が含まれる受信観測信号４２７を取得する。遅延検波部１６２１は、基準信号４２７１を用いて受信観測信号４２７に対して遅延検波を実行することにより、衛星７００の移動により生じる受信観測信号４２７のドップラー量を位相差で表す第１ドップラー量４８１を算出する。遅延検波部１６２１は、算出した第１ドップラー量４８１と基準信号４２７１とに基づいて相関を実行することにより第１ドップラー量４８１の確からしさを第１相関値４９１として算出する。すなわち、遅延検波部１６２１は、第１ドップラー量４８１を位相として算出する。
ここで、受信観測信号４２７は、ＡＩＳ信号である。受信観測信号４２７には、受信観測信号４２７を検出するための基準信号４２７１が含まれる。
基準信号４２７１は、受信観測信号４２７を検出するためのトレーニングシーケンス４２７１ａである。
【００９４】
同期検波部１６２２は、遅延検波部１６２１から取得した第１ドップラー量４８１と基準信号４２７１とを用いて受信観測信号４２７に対して同期検波を実行することにより受信観測信号４２７のドップラー量を周波数で表す第２ドップラー量４８２を算出する。同期検波部１６２２は、算出した第２ドップラー量４８２と基準信号４２７１とに基づいて相関を実行することにより第２ドップラー量４８２の確からしさを第２相関値４９２として算出する。
【００９５】
遅延検波部１６２１にて抽出される第１ドップラー量４８１は位相量即ちドップラー周波数に検出時間（この場合遅延信号の遅延量）を掛け合わせた形で算出される。
ここで、ドップラー位相（遅延検波検出量）をΘｄｏｐ、ドップラー周波数をＦｄｏｐ、遅延検波遅延時間をΔＴとすると、ΘｄｏｐとＦｄｏｐは以下の関係になる。但しπは円周率である。
（式１）Θｄｏｐ＝２×π×Ｆｄｏｐ×ΔＴ
上記（式１）の関係から、Θｄｏｐが分かるとＦｄｏｐが算出できる。同期検波部１６２２は周波数の変化を検波ループという周波数を少しずつ変化させて受信周波数に合わせる機能（以下チューニングと称する）を持っている。最初の合わせ込みで全く異なる周波数からチューニングを実施すると、合わせ込みすべき周波数に近づいてもそこで終了せず飛び越して異なる周波数まで到達してしまう等の問題が発生する。また合わせ込みの幅（例えば、衛星環境では±４ＫＨｚ）が大きいとそれを追従するために大きな帯域幅が必要となりそれだけ雑音を取り込んでしまい、収束した場合でも誤差が大きく発生してしまう等の問題も発生する。
そこで遅延検波部１６２１で算出したドップラー位相からドップラー周波数を割り出し、それを初期値としてチューニングを行うこととすると、追従のための帯域幅を非常に狭くでき、収束後の精度を高く保つことができる。なお、帯域を狭くしてしまうと収束までの時間がかかってしまい、トレーニングシーケンス内で収束できないという事象が発生する可能性がある。しかし、このような事象もある程度確からしいドップラー周波数を初期値に設定することで収束時間を短縮することができる。
【００９６】
検波出力判定部１６２３は、遅延検波部１６２１から出力された第１相関値４９１と第２相関値４９２とに基づいて、第１ドップラー量４８１と第２ドップラー量４８２とのいずれを受信観測信号４２７の信号ドップラー量４８３とするかを判定する。
【００９７】
同期検波部１６２２は、遅延検波部１６２１により算出された第１ドップラー量４８１を第１ドップラー周波数４８５として周波数に変換する。同期検波部１６２２は、変換した第１ドップラー周波数４８５に基づいて第１レプリカ信号４３３を生成する。同期検波部１６２２は、生成した第１レプリカ信号４３３を用いて同期検波を実行することにより第２ドップラー量４８２を算出する。
【００９８】
ドップラー推定装置１６２０は、信号ドップラー量４８３を判定するための相関閾値４９３を記憶する相関閾値記憶部１６２５を備える。
検波出力判定部１６２３は、第１相関値４９１と第２相関値４９２との各々を相関閾値４９３と比較すると共に、第１相関値４９１と第２相関値４９２とを比較する。検波出力判定部１６２３は、相関閾値４９３より大きく、かつ、他方の相関値よりも大きい相関値を特定し、特定した相関値に対応するドップラー量を信号ドップラー量４８３とする。
【００９９】
遅延検波部１６２１は、トレーニングシーケンス４２７１ａを用いて第１相関値４９１を算出する。また、同期検波部１６２２は、トレーニングシーケンス４２７１ａを用いて第２相関値４９２を算出する。
【０１００】
信号補正部１６２４は、検波出力判定部１６２３により判定された信号ドップラー量４８３に基づいて、受信観測信号４２７を補正する。信号補正部１６２４は、補正した受信観測信号４２７を補正受信観測信号４２８として出力する。信号ドップラー量４８３は、実施の形態１で説明した観測移動体ドップラー周波数４７１に相当する。
【０１０１】
＊＊＊動作の説明＊＊＊
図１４を用いて、本実施の形態に係るドップラー推定装置１６２０におけるドップラー推定方法について説明する。
ドップラー推定方法を実現するドップラー推定処理Ｓ１６２０は、遅延検波処理Ｓ１６２１、同期検波処理Ｓ１６２２、検波出力判定処理Ｓ１６２３、信号補正処理Ｓ１６２４を備える。
また、遅延検波処理Ｓ１６２１は、フィルタリング処理Ｓ１６２１１、最適タイミング取得処理Ｓ１６２１２、ドップラー推定補償処理Ｓ１６２１３を備える。
また、同期検波処理Ｓ１６２２は、フィルタリング処理Ｓ１６２２１、搬送波捕捉追尾処理Ｓ１６２２２を備える。
【０１０２】
遅延検波処理Ｓ１６２１において、遅延検波部１６２１は、衛星７００が受信する受信観測信号４２７であって基準信号４２７１が含まれる受信観測信号４２７を取得する。遅延検波部１６２１は、基準信号４２７１を用いて受信観測信号４２７に対して遅延検波を実行することにより、衛星７００の移動により生じる受信観測信号４２７のドップラー量を位相差で表す第１ドップラー量４８１を算出する。また、遅延検波部１６２１は、算出した第１ドップラー量４８１と基準信号４２７１とに基づいて相関を実行することにより第１ドップラー量４８１の確からしさを第１相関値４９１として算出する。ここで、第１相関値４９１を求めるために基準信号４２７１と相関を取る信号は、第１ドップラー量４８１（位相差）により補正された受信観測信号４２７である。遅延検波部１６２１は、第１ドップラー量４８１を位相差として算出する。
遅延検波部１６２１は、ドップラー推定結果である第１ドップラー量４８１を同期検波部１６２２に出力すると共に、トレーニングシーケンス相関結果である第１相関値４９１を検波出力判定部１６２３に出力する。
【０１０３】
Ｓ１６２１１において、フィルタリング部１６２１１は、ノイズを除去するために受信観測信号４２７をフィルタリングする。
Ｓ１６２１２において、最適タイミング取得部１６２１２は、ＡＩＳ信号である受信観測信号４２７に入っている基準信号４２７１であるトレーニングシーケンス４２７１ａのビットストリームから最適な復調タイミングを取得する。
Ｓ１６２１３において、ドップラー推定補償部１６２１３は、遅延検波によりドップラー周波数が位相に変換されているため位相差を計算することによりドップラー位相を推定する。ドップラー推定補償部１６２１３は、推定したドップラー位相を第１ドップラー量４８１として同期検波部１６２２に受け渡す。
【０１０４】
同期検波処理Ｓ１６２２において、同期検波部１６２２は、遅延検波部１６２１から取得した第１ドップラー量４８１と基準信号４２７１とを用いて受信観測信号４２７に対して同期検波を実行する。同期検波部１６２２は、この同期検波により受信観測信号４２７のドップラー量を周波数で表す第２ドップラー量４８２を算出する。同期検波部１６２２は、算出した第２ドップラー量４８２と基準信号４２７１とに基づいて相関を実行することにより第２ドップラー量４８２の確からしさを第２相関値４９２として算出する。ここで、第２相関値４９２を求めるために基準信号４２７１と相関を取る信号は、第２ドップラー量４８２（周波数）により補正された受信観測信号４２７である。同期検波部１６２２は、遅延検波部１６２１から取得したドップラー推定結果である第１ドップラー量４８１を用いて、同期検波し、搬送波再生を行う。同期検波部１６２２は、搬送波再生時に出力されるトレーニングシーケンス相関結果である第２相関値４９２を検波出力判定部１６２３に出力する。
また、同期検波部１６２２は、遅延検波部１６２１により算出された第１ドップラー量４８１を第１ドップラー周波数４８５として周波数に変換する。同期検波部１６２２は、変換した第１ドップラー周波数４８５に基づいて第１レプリカ信号４３３を生成する。同期検波部１６２２は、生成した第１レプリカ信号４３３を用いて同期検波を実行することにより第２ドップラー量４８２を算出する。
【０１０５】
Ｓ１６２２１において、フィルタリング部１６２２１は、ノイズを除去するために受信観測信号４２７をフィルタリングする。
Ｓ１６２２２において、搬送波捕捉追尾部１６２２２は、遅延検波部１６２１から取得した第１ドップラー量４８１であるドップラー位相を第１ドップラー周波数４８５変換する。搬送波捕捉追尾部１６２２２は、同期検波部１６２２の搬送波再生局部発振器の基準周波数を、変換した第１ドップラー周波数４８５分シフトし、搬送波レプリカを生成する。搬送波捕捉追尾部１６２２２は、生成した搬送波レプリカと受信観測信号４２７との信号同期を行う。搬送波捕捉追尾部１６２２２は、搬送波レプリカと受信観測信号４２７との信号同期の結果を第２相関値４９２として検波出力判定部１６２３に出力する。
【０１０６】
検波出力判定処理Ｓ１６２３において、検波出力判定部１６２３は、遅延検波部１６２１から出力された第１相関値４９１と第２相関値４９２とに基づいて、第１ドップラー量４８１と第２ドップラー量４８２とのいずれを受信観測信号４２７の信号ドップラー量４８３とするかを判定する。検波出力判定処理Ｓ１６２３は、遅延検波部１６２１から出力されたトレーニングシーケンス相関結果である第１相関値４９１と、同期検波部１６２２から出力されたトレーニングシーケンス相関結果である第２相関値４９２とを取得する。検波出力判定部１６２３は、相関閾値記憶部１６２５により記憶された相関閾値４９３を用いて、第１相関値４９１と第２相関値４９２とを判定する。
【０１０７】
検波出力判定部１６２３は、第１相関値４９１と第２相関値４９２との各々を相関閾値４９３と比較すると共に、第１相関値４９１と第２相関値４９２とを比較する。検波出力判定部１６２３は、相関閾値４９３より大きく、かつ、他方の相関値よりも大きい相関値を特定し、特定した相関値に対応するドップラー量を信号ドップラー量４８３とする。なお、第１相関値４９１と第２相関値４９２とのどちらもが相関閾値４９３より小さい場合、すなわち相関値が規定を満たさない場合には次の受信観測信号４２７のフレームを待つ。
なお、信号ドップラー量４８３は、実施の形態１で説明した観測移動体ドップラー周波数４７１に相当する。
【０１０８】
信号補正処理Ｓ１６２４において、信号補正部１６２４は、検波出力判定部１６２３により判定された信号ドップラー量４８３に基づいて、受信観測信号４２７を補正する。信号補正部１６２４は、補正した受信観測信号４２７を補正受信観測信号４２８として出力する。
【０１０９】
＊＊＊効果の説明＊＊＊
ドップラー推定において、同期検波を直接利用した場合、以下のようにドップラーを追従することが困難である。すなわち、９６００ＢＰＳであり周波数に換算すると約２．５ＫＨｚ〜１０ＫＨｚであるＡＩＳ信号にスペクトラムが生成されると、約−４ＫＨｚ〜＋４ＫＨｚのドップラー周波数が混在するためである。約−４ＫＨｚ〜＋４ＫＨｚのドップラー周波数は、衛星を高度約５００Ｋｍの周回軌道上にあると仮定した場合の数値である。
【０１１０】
本実施の形態に示す方式を用いた場合、ドップラー周波数は位相変動の形式で推定され、またＡＩＳ信号の１フレーム間での推定を行うため、その間のドップラーの変動は小さく、定常位相として取り扱うことができ、精度の高いドップラー量の推定が可能である。
【０１１１】
本実施の形態では、ドップラー推定のための検波方式、すなわち遅延検波と、復調のための検波方式、すなわち同期検波とをハイブリッドに構成することにより、遅延検波そのものでは達成できない高性能なＢＥＲを達成できる。ＢＥＲ（Ｂｉｔ  Ｅｒｒｏｒ  Ｒａｔｅ）とは送信ビット総数に対するエラービット数の比である。遅延検波ではＢＥＲ１０−５に対してＥｂ／Ｎｏ約１８〜１９ｄＢであるが、同期検波では同じＢＥＲに対してＥｂ／Ｎｏ約１０ｄＢ程度を実現可能であり８〜９ｄＢ程度の性能改善が可能である。
【０１１２】
また、ドップラー推定のための遅延検波方式と復調のための同期検波方式とをハイブリッドに構成することで、同期検波を直接利用する場合の高速フェージング発生時の壊滅的復調失敗が見込まれる場合には、遅延検波側の出力をそのまま利用した復調が可能となる。よって、復調失敗の発生率を大幅に低減できる。
【０１１３】
信号検波システム５００、ドップラー推定装置１６２０の機能ブロックは、上記実施の形態に記載した機能を実現することができれば、任意である。信号検波システム５００、ドップラー推定装置１６２０の機能ブロックでは、機能ブロックを、他のどのような組み合わせで信号検波システム５００、ドップラー推定装置１６２０を構成しても構わない。
【０１１４】
また、信号検波システム５００、ドップラー推定装置１６２０は、１つの装置でなく、複数の装置から構成されたシステムでもよいし、あるいは１つのシステム、装置に組み込まれた部品でも構わない。
【０１１５】
以上、実施の形態１，２について説明したが、これらの２つの実施の形態のうち、１つの実施の形態を部分的に実施しても構わない。あるいは、これらの２つの実施の形態のうち、複数の部分を組み合わせて実施しても構わない。その他、これらの２つの実施の形態を、全体としてあるいは部分的に、どのように組み合わせて実施しても構わない。
なお、上記の実施の形態は、本質的に好ましい例示であって、本発明、その適用物や用途の範囲を制限することを意図するものではなく、必要に応じて種々の変更が可能である。
【符号の説明】
【０１１６】
１１０  衛星情報、１１１  衛星テレメータ情報、１１２  軌道決定値情報、１１３  条件テーブル、１２０  衛星模擬部、１２１  地球重力モデル部、１２２  アンテナ姿勢モデル部、１２３  衛星軌道計算部、１３０  検波計画部、１４０  模擬信号生成部、１４１  移動体信号模擬部、１４２  ＧＭＳＫ変調部、１４３  ドップラー周波数挿入部、１４４  時間遅延補正部、１４５  信号加算部、１４６  ＧＭＳＫ復調部、１４９  複数信号生成部、１５０  信号模擬情報、１５１  衛星搭載アンテナモデル、１６０  信号処理部、１６１  信号分離部、１６２  ドップラー推定部、１６３  信号復号部、１６４  復調部、１７０  解析結果ＤＢ、１８１  衛星受信信号記憶部、１８２  処理情報、１９１  データ加工部、１９２  マニュアル入力部、４００  衛星受信信号、４０１  観測信号、４０２  観測外信号、４２０  検波時衛星模擬情報、４２１  検波期間衛星状態、４２３  検波時刻、４２４  検波期間、４２５  観測外移動体情報、４２６  模擬信号、４２７  受信観測信号、４２８  補正受信観測信号、４２９  移動体送信信号、４３０  検波期間観測外移動体情報、４３１  復号移動体情報、４３３  第１レプリカ信号、４７１  観測移動体ドップラー周波数、４８１  第１ドップラー量、４８２  第２ドップラー量、４８３  信号ドップラー量、４８５  第１ドップラー周波数、４９１  第１相関値、４９２  第２相関値、４９３  相関閾値、５００  信号検波システム、６００  移動体、６０１  観測移動体、６０２  観測外移動体、６１０  移動体情報、７００  衛星、７０１  観測エリア、９０１  プロセッサ、９０２  補助記憶装置、９０３  メモリ、９０４  通信装置、９０５  入力インタフェース、９０６  ディスプレイインタフェース、９０７  入力装置、９０８  ディスプレイ、９１０  信号線、９１１，９１２  ケーブル、９９０  プログラム、１４１１  移動体信号、１６２０  ドップラー推定装置、１６２１  遅延検波部、１６２２  同期検波部、１６２３  検波出力判定部、１６２４  信号補正部、１６２５  相関閾値記憶部、４２７１  基準信号、４２７１ａ  トレーニングシーケンス、４２９０  複数の移動体送信信号、６０２０  複数の観測外移動体、９４０１  レシーバー、９４０２  トランスミッター、１４１１０  複数の移動体信号、１６２１１  フィルタリング部、１６２１２  最適タイミング取得部、１６２１３  ドップラー推定補償部、１６２２１  フィルタリング部、１６２２２  搬送波捕捉追尾部、Ｓ１００  検波時刻決定処理、Ｓ１１０  衛星模擬処理、Ｓ１２０  検波計画処理、Ｓ１３０  模擬信号生成処理、Ｓ１３０ａ  複数信号生成処理、Ｓ１３１  移動体信号模擬処理、Ｓ１３２  ＧＭＳＫ変調処理、Ｓ１３３  ドップラー周波数挿入処理、Ｓ１３４  時間遅延補正処理、Ｓ１３５  信号加算処理、Ｓ１３６  ＧＭＳＫ復調処理、Ｓ１４０  信号処理、Ｓ１５０  信号分離処理、Ｓ１５１  信号相関処理、Ｓ１６０  復調処理、Ｓ１６１  ドップラー推定処理、Ｓ１６２  信号復号処理、Ｓ１６２ａ  シンボルタイミング検知処理、Ｓ１６２ｂ  ＣＲＣ／ＳＦ／ＥＦ／ビットスタッフィング判定処理、Ｓ１６２ｃ  工学値変換処理、Ｓ１６２０  ドップラー推定処理、Ｓ１６２１  遅延検波処理、Ｓ１６２２  同期検波処理、���１６２３  検波出力判定処理、Ｓ１６２４  信号補正処理、Ｓ１６２１１  フィルタリング処理、Ｓ１６２１２  最適タイミング取得処理、Ｓ１６２１３  ドップラー推定補償処理、Ｓ１６２２１  フィルタリング処理、Ｓ１６２２２  搬送波捕捉追尾処理。
【図１】
0006445329.tif 000002
【図２】
0006445329.tif 000003
【図３】
0006445329.tif 000004
【図４】
0006445329.tif 000005
【図５】
0006445329.tif 000006
【図６】
0006445329.tif 000007
【図７】
0006445329.tif 000008
【図８】
0006445329.tif 000009
【図９】
0006445329.tif 000010
【図１０】
0006445329.tif 000011
【図１１】
0006445329.tif 000012
【図１２】
0006445329.tif 000013
【図１３】
0006445329.tif 000014
【図１４】
0006445329.tif 000015
【図１５】
0006445329.tif 000016