(57)【要約】【課題】  離隔配置による時刻誤差の影響を軽減し、高精度な位置を出力する。【解決手段】  実施形態のレーダシステムにおいて、送信装置は、パルス圧縮された送信パルスを通信期間とレンジ期間に分割し、通信期間のパルスに送信に関する通信情報を重畳して送信する。受信装置は、前記送信装置から送信される直接波または目標からの反射波を受信し、その受信信号を前記レンジ期間の参照信号で圧縮し、信号対雑音電力比の高い反射点を用いて前記通信期間を抽出して相関処理することにより前記通信情報を抽出し、信号対雑音電力比の低い他の反射点についても、前記レンジ期間のみまたは前記通信情報を用いて、前記通信期間及び前記レンジ期間のパルス全体で相関処理することによって目標の距離を出力する。【選択図】図１
2019095391.tif 000002
【特許請求の範囲】
【請求項１】
パルス圧縮用の送信パルスを通信期間とレンジ期間に分割し、通信期間のパルスに送信に関する通信情報を重畳して送信する送信装置と、
前記送信装置から送信される直接波または目標からの反射波を受信し、その受信信号を前記レンジ期間の参照信号で圧縮し、信号対雑音電力比（ＳＮ）の高い反射点を用いて前記通信期間を抽出して相関処理することにより前記通信情報を抽出し、信号対雑音電力比の低い他の反射点についても、前記レンジ期間のみまたは前記通信情報を用いて、前記通信期間及び前記レンジ期間のパルス全体で相関処理することによって目標の距離を出力する受信装置と
を具備するレーダシステム。
【請求項２】
パルスを繰り返すレーダのパルス列を通信期間とドップラ期間とレンジ期間に時分割し、通信期間のパルスに送信に関する通信情報を重畳して送信する送信装置と、
前記送信装置から送信される直接波または目標からの反射波を受信し、その受信信号を前記ドップラ期間のＦＦＴ結果のうち、信号対雑音電力比（ＳＮ）の高い反射点のドップラを用いて前記レンジ期間の参照信号を補正し、補正した参照信号により前記レンジ期間を相関処理してレンジを抽出し、通信期間の受信期間を抽出して前記ドップラで補正した参照信号により相関処理することによって前記通信情報を抽出し、信号対雑音電力比の低い他の反射点についても、前記通信期間は前記通信情報と前記ドップラを用いて参照信号を補正し、前記レンジ期間は前記ドップラを用いて参照信号を補正し、前記通信期間、前記ドップラ期間、前記レンジ期間のうち、少なくとも前記レンジ期間を含むパルス列で用いて相関処理して目標の距離を出力する受信装置と
を具備するレーダシステム。
【請求項３】
パルスを繰り返すレーダのパルス列を通信期間とドップラ期間とレンジ期間の３通り、または通信期間及びドップラ期間とレンジ期間の２通りのパルス繰り返し周波数のパルス列とし、通信期間のパルスに送信に関する通信情報を重畳して、前記通信期間、前記ドップラ期間、前記レンジ期間のパルス列を混合して送信する送信装置と、
前記送信装置から送信される直接波または目標からの反射波を受信し、その受信信号を、ドップラ期間のＦＦＴ結果のうち、ＳＮ比の高い反射点のドップラを用いて前記レンジ期間の参照信号を補正し、補正した参照信号により前記レンジ期間を相関処理してレンジを抽出し、前記通信期間の受信信号を抽出して前記ドップラで補正した参照信号により相関処理することによって前記通信情報を抽出し、信号対雑音電力比の低い他の反射点についても、前記レンジ期間は前記ドップラを用いて参照信号を補正して相関処理することで目標の距離を出力する受信装置と
を具備するレーダシステム。
【請求項４】
パルスを繰り返すレーダのパルス列をドップラ期間と通信期間及びレンジ期間の２通りのパルス繰り返し周波数のパルス列とし、通信期間のパルスに送信に関する通信情報を重畳して、前記通信期間、前記ドップラ期間、前記レンジ期間のパルス列を混合して送信する送信装置と、
前記送信装置から送信される直接波または目標からの反射波を受信し、その受信信号を、ドップラ期間のＦＦＴ結果のうち、信号対雑音電力比（ＳＮ）の高い反射点のドップラを用いて前記レンジ期間の参照信号を補正し、補正した参照信号により前記レンジ期間を相関処理してレンジを抽出し、前記通信期間の受信信号を抽出して前記ドップラで補正した参照信号により相関処理することによって前記通信情報を抽出し、信号対雑音電力比の低い他の反射点についても、前記通信期間は前記通信情報と前記ドップラを用いて第１の参照信号を補正し、前記レンジ期間は前記ドップラを用いて第２の参照信号を補正して、第１の参照信号と第２の参照信号を合わせた参照信号を用いて相関処理することで目標の距離を出力する受信装置と
を具備するレーダシステム。
【請求項５】
送信装置が、パルス圧縮用の送信パルスを通信期間とレンジ期間に分割し、通信期間のパルスに送信に関する通信情報を重畳して送信し、
受信装置が、前記送信装置から送信される直接波または目標からの反射波を受信し、その受信信号を前記レンジ期間の参照信号で圧縮し、信号対雑音電力比（ＳＮ）の高い反射点を用いて前記通信期間を抽出して相関処理することにより前記通信情報を抽出し、信号対雑音電力比の低い他の反射点についても、前記レンジ期間のみまたは前記通信情報を用いて、前記通信期間及び前記レンジ期間のパルス全体で相関処理することによって目標の距離を出力する
レーダシステムのレーダ信号処理方法。
【請求項６】
送信装置が、パルスを繰り返すレーダのパルス列を通信期間とドップラ期間とレンジ期間に時分割し、通信期間のパルスに送信に関する通信情報を重畳して送信し、
受信装置が、前記送信装置から送信される直接波または目標からの反射波を受信し、その受信信号を前記ドップラ期間のＦＦＴ結果のうち、信号対雑音電力比（ＳＮ）の高い反射点のドップラを用いて前記レンジ期間の参照信号を補正し、補正した参照信号により前記レンジ期間を相関処理してレンジを抽出し、通信期間の受信期間を抽出して前記ドップラで補正した参照信号により相関処理することによって前記通信情報を抽出し、信号対雑音電力比の低い他の反射点についても、前記通信期間は前記通信情報と前記ドップラを用いて参照信号を補正し、前記レンジ期間は前記ドップラを用いて参照信号を補正し、前記通信期間、前記ドップラ期間、前記レンジ期間のうち、少なくとも前記レンジ期間を含むパルス列で用いて相関処理して目標の距離を出力する
レーダシステムのレーダ信号処理方法。
【請求項７】
送信装置が、パルスを繰り返すレーダのパルス列を通信期間とドップラ期間とレンジ期間の３通り、または通信期間及びドップラ期間とレンジ期間の２通りのパルス繰り返し周波数のパルス列とし、通信期間のパルスに送信に関する通信情報を重畳して、前記通信期間、前記ドップラ期間、前記レンジ期間のパルス列を混合して送信し、
受信装置が、前記送信装置から送信される直接波または目標からの反射波を受信し、その受信信号を、ドップラ期間のＦＦＴ結果のうち、信号対雑音電力比（ＳＮ）の高い反射点のドップラを用いて前記レンジ期間の参照信号を補正し、補正した参照信号により前記レンジ期間を相関処理してレンジを抽出し、前記通信期間の受信信号を抽出して前記ドップラで補正した参照信号により相関処理することによって前記通信情報を抽出し、信号対雑音電力比の低い他の反射点についても、前記レンジ期間は前記ドップラを用いて参照信号を補正して相関処理することで目標の距離を出力する
レーダシステムのレーダ信号処理方法。
【請求項８】
送信装置が、パルスを繰り返すレーダのパルス列をドップラ期間と通信期間及びレンジ期間の２通りのパルス繰り返し周波数のパルス列とし、通信期間のパルスに送信に関する通信情報を重畳して、前記通信期間、前記ドップラ期間、前記レンジ期間のパルス列を混合して送信し、
受信装置が、前記送信装置から送信される直接波または目標からの反射波を受信し、その受信信号を、ドップラ期間のＦＦＴ結果のうち、信号対雑音電力比（ＳＮ）の高い反射点のドップラを用いて前記レンジ期間の参照信号を補正し、補正した参照信号により前記レンジ期間を相関処理してレンジを抽出し、前記通信期間の受信信号を抽出して前記ドップラで補正した参照信号により相関処理することによって前記通信情報を抽出し、信号対雑音電力比の低い他の反射点についても、前記通信期間は前記通信情報と前記ドップラを用いて第１の参照信号を補正し、前記レンジ期間は前記ドップラを用いて第２の参照信号を補正して、第１の参照信号と第２の参照信号を合わせた参照信号を用いて相関処理することで目標の距離を出力する
レーダシステムの信号処理方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本実施形態は、マルチスタティック方式のレーダシステム及びそのレーダ信号処理方法に関する。
【背景技術】
【０００２】
近時、レーダシステムにあっては、複数の送受信レーダ装置または複数の送信装置及び受信レーダ装置を互いに離間して配置し、各レーダ装置の観測結果により目標の位置を検出するマルチスタティック方式が開発されている。ただし、この種のレーダシステムでは、離隔したレーダ装置間の時刻同期が不十分な場合や距離精度が不十分な場合には、観測位置の誤差が大きくなる。また、送信装置が見通し外にある場合には、直接波を受信することができず、同期したマルチスタティック動作ができなくなってしまう。
【先行技術文献】
【非特許文献】
【０００３】
【非特許文献１】パルス圧縮、大内、‘リモートセンシングのための合成開口レーダの基礎’、東京電機大学出版局、pp.131-149(2003)
【非特許文献２】位相モノパルス（位相比較モノパルス）方式、吉田、‘改訂レーダ技術’、電子情報通信学会、pp.262-264（1996）
【非特許文献３】ＣＦＡＲ（Constant False Alarm Rate）、吉田、‘改訂レーダ技術’、電子情報通信学会、pp.87-89（1996）
【非特許文献４】ＢＰＳＫ（Binary Phase Shift Keying）、ＱＰＳＫ（Quadrature Phase Shift Keying）、西村、‘ディジタル信号処理による通信システム設計、ＣＱ出版社、pp.222-226(2006)
【非特許文献５】テーラー分布、吉田、‘改訂レーダ技術’、電子情報通信学会、pp.134-135(1996)
【発明の概要】
【発明が解決しようとする課題】
【０００４】
以上述べたように、従来のマルチスタティック方式によるレーダシステムでは、レーダ装置間の時刻同期ずれや中心周波数のずれ等の影響で、ビート周波数を精度よく観測できず、目標の位置を算出する上で距離精度が不十分であるという課題があった。
【０００５】
本実施形態は上記課題に鑑みなされたもので、レーダ装置間の時刻同期ずれや中心周波数のずれ等の影響を軽減してビート周波数の観測精度を向上させめことができ、目標の位置を算出する上で十分な距離精度が得られるレーダシステム及びそのレーダ信号処理方法を提供することを目的とする。
【課題を解決するための手段】
【０００６】
上記の課題を解決するために、本実施形態に係るレーダシステムは、送信装置と受信装置を備える。送信装置は、パルス圧縮用の送信パルスを通信期間とレンジ期間に分割し、通信期間のパルスに送信に関する通信情報を重畳して送信する。受信装置は、前記送信装置から送信される直接波または目標からの反射波を受信し、その受信信号を前記レンジ期間の参照信号で圧縮し、信号対雑音電力比（ＳＮ）の高い反射点を用いて前記通信期間を抽出して相関処理することにより前記通信情報を抽出し、信号対雑音電力比の低い他の反射点についても、前記レンジ期間のみまたは前記通信情報を用いて、前記通信期間及び前記レンジ期間のパルス全体で相関処理することによって目標の距離を出力する。
【図面の簡単な説明】
【０００７】
【図１】第１の実施形態に係るレーダシステムの全体構成を示すブロック図。
【図２】第１の実施形態に係るレーダシステムの送信装置の構成を示すブロック図。
【図３】第１の実施形態に係るレーダシステムの受信装置の構成を示すブロック図。
【図４】第１の実施形態において、目標に対して送信ビーム及び受信ビームが形成される様子を示す図。
【図５】第１の実施形態に係るレーダシステムにおいて、通信情報の送受信による目標位置同定の距離観測及びＡＺ，ＥＬ観測を説明するための座標系を示す概念図。
【図６】第１の実施形態に係るレーダシステムにおいて、レンジによるフィッティングを用いて目標の位置同定を行う様子を３次元座標系で示す概念図。
【図７】第１の実施形態に係るレーダシステムの送信波形を示す波形図。
【図８】第１の実施形態に係るレーダシステムにおいて、送信パルス内の通信変調信号を生成する様子を示す波形図。
【図９】第１の実施形態に係るレーダシステムにおいて、受信信号から通信期間を抽出して通信信号を復調する様子を示す波形図。
【図１０】第１の実施形態に係るレーダシステムにおいて、受信信号のレンジ期間Ｐ２列の相関処理から通信期間を抽出して通信信号を復調する様子を示す波形図。
【図１１】第１の実施形態に係るレーダシステムにおいて、受信信号のレンジ期間Ｐ２＋通信期間Ｐ３の相関処理から最大バンク信号の距離を抽出して通信信号を復調する様子を示す波形図。
【図１２】第２の実施形態に係るレーダシステムの送信装置の構成を示すブロック図。
【図１３】第２の実施形態に係るレーダシステムの受信装置の構成を示すブロック図。
【図１４】第２の実施形態に係るレーダシステムにおいて、ドップラ期間、レンジ期間、通信期間それぞれのパルス列の時分割合成による送信波形を示す波形図。
【図１５】第２の実施形態に係るレーダシステムにおいて、ドップラ期間、レンジ期間、通信期間それぞれのパルス列を時分割合成した送信信号に、通信情報の変調信号でパルス間を変調する様子を示す波形図。
【図１６】第２の実施形態に係るレーダシステムにおいて、受信信号からドップラを抽出する様子を示す波形図。
【図１７】第２の実施形態に係るレーダシステムにおいて、受信信号から通信期間を抽出して通信情報を復調する様子を示す波形図。
【図１８】第２の実施形態に係るレーダシステムにおいて、受信信号から距離を抽出する様子を示す波形図。
【図１９】第３の実施形態に係るレーダシステムの送信装置の構成を示すブロック図。
【図２０】第３の実施形態に係るレーダシステムの受信装置の構成を示すブロック図。
【図２１】第３の実施形態に係るレーダシステムにおいて、ドップラ期間、レンジ期間、通信期間それぞれのパルス列を合成して送信信号を生成する様子を示す波形図。
【図２２】第３の実施形態に係るレーダシステムにおいて、パルス列を合成した送信信号に通信情報の変調信号でパルス間を変調する様子を示す波形図。
【図２３】第３の実施形態に係るレーダシステムにおいて、受信信号からＰＲＦが異なるパルス列を分離する様子を示す波形図。
【図２４】第３の実施形態に係るレーダシステムにおいて、受信信号から通信期間を抽出して通信情報を復調する様子を示す波形図。
【図２５】第３の実施形態に係るレーダシステムにおいて、受信信号からドップラを抽出する様子を示す波形図。
【図２６】第３の実施形態に係るレーダシステムにおいて、受信信号から距離を抽出する様子を示す波形図。
【発明を実施するための形態】
【０００８】
以下、実施形態について、図面を参照して説明する。尚、各実施形態の説明において、同一部分には同一符号を付して示し、重複する説明を省略する。
【０００９】
（第１の実施形態）
図１乃至図１１を参照して、第１の実施形態に係るレーダシステムを説明する。
【００１０】
図１は本実施形態に係るレーダシステムの全体構成を示し、図２はその送信装置の構成を示し、図３はその受信装置の構成を示している。
【００１１】
図１に示すレーダシステムは、互いに同構成のＮ個の送信装置Ｔ１〜ＴＮ（図１ではＮ＝３）と受信装置Ｒを備えるマルチスタティック方式を採用している。送信装置Ｔ１〜ＴＮは、送信位置、送信周波数、送信タイミング等の通信情報をレーダ送信波に重畳し、目標Ｔに向けて送出する。一方、受信装置Ｒは、目標Ｔで反射されたレーダ送信波を受信して、目標Ｔの距離、角度、速度を観測するとともに、受信信号から通信情報を抽出し、その通信情報に基づいて目標の観測結果を補正する機能を有する。
【００１２】
上記送信装置Ｔｉ（ｉは１〜Ｎのいずれか）は、図２に示すように、送信アンテナ装置Ｔ１１、パルス変調器Ｔ１２、周波数変換器Ｔ１３、変調器Ｔ１４、制御器Ｔ１５、チャープ信号生成器Ｔ１６、ＩＦＦＴ（Inverse Fast Fourier Transform）処理器Ｔ１７及び通信変調信号生成器Ｔ１８を備える。
【００１３】
上記送信アンテナ装置Ｔ１１は繰り返し供給される特定周波数の送信信号を指定方向に送出する。一方、通信変調信号生成器Ｔ１８は、送信位置、送信周波数、送信タイミング、送信ビーム指向方向等を含む通信情報をビット（０，１）で表し、ＢＰＳＫ（非特許文献４参照）等により変調することで、通信変調信号を生成する。この通信変調信号は、ＩＦＦＴ処理器Ｔ１７で逆フーリエ変換され、チャープ信号生成器Ｔ１６で生成されるチャープ信号と共に制御器Ｔ１５に送られる。制御器Ｔ１５は、送信パルスを通信期間とレンジ期間に分割し、Ｎヒット（Ｎ＞０）送信する。レンジ期間では、チャープ信号と通信変調信号について、パルス内制御によりパルス内のタイミングを制御し、変調器Ｔ１４により送信パルスを生成し、周波数変換器Ｔ１３で送信周波数に変換し、パルス変調器Ｔ１２でパルスゲートによるパルス変調を施し、送信アンテナ装置Ｔ１１で高出力増幅して送信する。
【００１４】
受信装置Ｒ１は、図３に示すように、受信アンテナ装置Ｒ１１、周波数変換器Ｒ１２、ＡＤ（Analog-Digital）変換器Ｒ１３、Ｐ２抽出器Ｒ１４、相関処理パルス圧縮器Ｒ１５、ＦＦＴ（Slow-time）処理器Ｒ１６、検出器Ｒ１７、Ｐ３範囲選定器Ｒ１８、Ｐ３抽出器Ｒ１９、相関処理器（Slow-time）Ｒ２０、復調器Ｒ２１、Ｐ２＋Ｐ３選定器Ｒ２２、通信情報参照信号補正器Ｒ２３、相関処理器Ｒ２４、ＦＦＴ（Slow-time）処理器Ｒ２５、検出器Ｒ２６、レンジ抽出器Ｒ２７、出力処理器Ｒ２８を備える。
【００１５】
すなわち、受信アンテナ装置Ｒ１１で受信した信号は、周波数変換器Ｒ１２でベースバンドに周波数変換され、ＡＤ変換器Ｒ１３でデジタル信号に変換される。そして、Ｐ２抽出器Ｒ２４により、レンジ期間のパルスＰ２の変調信号（チャープ等）が抽出され、相関処理器Ｒ１５で相関処理（パルス圧縮）が施され、ＦＦＴ処理器Ｒ１６でパルス間（slow-time軸）のＦＦＴが行われ、検出器Ｒ１７にて、最大振幅のバンク信号を用いてＣＦＡＲ（非特許文献４）等による検出が行われる。ここで、複数の検出点がある場合には、ＳＮ（信号対雑音電力比）の大きな信号を抽出することができるので、Ｐ３範囲選定器Ｒ１８において、その時間軸を基準に通信期間Ｐ３の時間範囲を選定し、Ｐ３抽出器Ｒ１９で通信期間Ｐ３を抽出する。
【００１６】
次に、ＦＦＴ処理器Ｒ２０で、slow-time軸のＦＦＴ処理を施して高ＳＮとした後、復調器Ｒ２１で復調して通信情報を得る。通信情報の復調により、送信源の送信時刻、送信周波数、送信位置、通信期間の参照信号を取得できる。
【００１７】
また、Ｐ２＋Ｐ３選定器Ｒ２２において、ＡＤ変換出力から、レンジ期間Ｐ２と通信期間Ｐ３を選定し、通信情報参照信号補正器Ｒ２３で通信期間Ｐ３の参照信号をレンジ期間Ｐ２に合わせ、送信処理器Ｒ２４でＰ２＋Ｐ３の相関処理を行い、ＦＦＴ処理器Ｒ２５でslow-time軸のＦＦＴを行い、検出器Ｒ２６でＣＦＡＲ等により目標を検出し、レンジ抽出器Ｒ２７でレンジを抽出し、出力処理器Ｒ２８で所定の形式に変換して出力する。
【００１８】
本実施形態に係るレーダシステムでは、送信装置Ｔｉと受信装置Ｒとの間が離隔しているため、送信位置、送信周波数、送信タイミング、送信ビーム指向方向等の情報を共有する必要がある。特に、送信位置が移動する場合には、リアルタイムに情報を共有する必要がある。送信装置と受信装置との間で見通しがある場合には、通信回線を確保することも可能であるが、見通し外の場合には、通信回線の確保は困難になる。この対策として、レーダ送信波に通信情報を重畳し、目標からの反射信号として復調すれば、通信が可能となる。これは、見通し内の場合でも、レーダの直接波を用いて通信できれば、通信回線を別に備える必要がなくなるので有用である。本実施形態では、パルス内で通信情報を重畳し、測距を行う手法について述べる。なお、測角処理については、ΣとΔビームによるモノパルス処理等の一般的な手法（非特許文献２）を用いることができるので、記述を割愛する。
【００１９】
上記構成によるレーダシステムにおいて、図４にマルチスタティックシステムの目標Ｔに対する送信ビーム及び受信ビームの様子を示し、図５に通信情報の送受信による目標位置同定の距離観測及びＡＺ，ＥＬ観測を説明するための座標系を示し、図６にレンジによるフィッティングを用いて目標の位置同定を行う様子を３次元座標系で示し、図７に送信波形を示し、図８に送信パルス内の通信変調信号を生成する様子を示し、図９に受信信号から通信期間を抽出して通信信号を復調する様子を示し、図１０に受信信号のレンジ期間Ｐ２列の相関処理から通信期間を抽出して通信信号を復調する様子を示し、図１０ａに受信信号のレンジ期間Ｐ２＋通信期間Ｐ３の相関処理から最大バンク信号の距離を抽出して通信信号を復調する様子を示す。なお、図４では、送信ビームをペンシルビームで表現しているが、送信ビームは無指向性を含むファンビームでもよい。
【００２０】
本実施形態では、受信装置Ｒで、目標Ｔの距離とＡＺ及びＥＬの角度を観測すれば、図５に示すように、目標Ｔの位置を同定できる。または、図６に示すように、複数の送信〜受信の距離を観測すれば、距離によるフィッティングによって楕円球面の交点等から目標位置を同定できる。そこで、本実施形態は、送信装置及び受信装置間の時刻同期を高精度に調整する手法を提供する。
【００２１】
送信装置Ｔｉにおいて、送信信号は、図７に示すように、送信パルスが通信期間Ｐ３とレンジ期間Ｐ２に分割され、Ｎヒット（Ｎ＞０）送信とする。通信情報としては、送信位置、送信周波数、送信タイミング、送信ビーム指向方向等をビット（０，１）で表し、図８に示すように、ＢＰＳＫ（非特許文献４）等により変換した信号を生成し（Ｔ１８）、ＩＦＦＴして（Ｔ１７）、パルス内の通信変調信号を生成する。レンジ期間については、チャープ信号を生成し（Ｔ１６）、通信変調信号とともに、パルス内制御によりパルス内のタイミングを制御し（Ｔ１５）、通信期間Ｐ３とレンジ期間Ｐ２を含む送信パルスを生成し（Ｔ１４）、送信周波数に周波数変換して（Ｔ１３）、パルスゲートによるパルス変調を施して（Ｔ１２）、高出力増幅し、目標に向けて送出する（Ｔ１１）。
【００２２】
受信装置Ｒにおいて、受信アンテナ装置Ｒ１１で受信した信号は、ベースバンド信号に周波数変換され（Ｒ１２）、ＡＤ変換される（Ｒ１３）。ＡＤ変換された受信信号は。まず、図９に示すように、レンジ期間Ｐ２が抽出され（Ｒ１４）、変調信号（チャープ等）による相関処理（パルス圧縮）が施され（Ｒ１５）、図１０に示すようにパルス間（slow-time軸）のＦＦＴが行われ（Ｒ１６）、最大振幅のバンク信号を用いてＣＦＡＲ（非特許文献４）等による検出が行われる（Ｒ１７）。ここで、複数の検出点がある場合には、ＳＮ（信号対雑音電力比）の大きな信号を抽出できるので、その時間軸を基準に通信期間Ｐ３の時間範囲を選定し（Ｒ１８）、通信期間Ｐ３を抽出する（Ｒ１９）。次に、slow-time軸で相関処理して高ＳＮとした後（Ｒ２０）、復調して通信情報を得る（Ｒ２１）。
【００２３】
ここで、通信情報を復調することにより、送信源の送信時刻を取得できれば、送信と受信の時刻同期を行うことができる。また、送信周波数がわかれば、ローカル信号の同調を行うことができ、正しいドップラを得ることができる。また、送信位置がわかれば、図５や図６のマルチスタティックによる目標の位置同定ができる。また、通信情報を復調できれば、通信期間の参照信号がわかるので、受信信号からレンジ期間Ｐ２と通信期間Ｐ３の選定を行い（Ｒ２２）、通信期間Ｐ３の参照信号補正を行って、図１１に示すように、レンジ期間Ｐ２を合わせた参照信号を生成し（Ｒ２３）、Ｐ２＋Ｐ３の相関処理を行い（Ｒ２４）、slow-time軸のＦＦＴを行い（Ｒ２５）、ＣＦＡＲ等により目標を検出し（Ｒ２６）、レンジを抽出し（Ｒ２７）、出力する（Ｒ２８）。
【００２４】
以下に、変調信号について定式化して説明する。まず、レンジ期間Ｐ２の変調信号は、チャープ変調とすると、次式で表現できる。  
【００２５】
【数１】
2019095391.tif 000003
【００２６】
本実施形態では、チャープ変調で述べるが、符号変調方式等他の変調方式でもよい。
【００２７】
次に通信期間Ｐ３の変調波形を考えると、例えばＢＰＳＫ（非特許文献４）の場合は次式となる。  
【００２８】
【数２】
2019095391.tif 000004
【００２９】
【数３】
2019095391.tif 000005
【００３０】
これを用いて、パルス幅内の通信期間Ｐ３の変調を行う。
【００３１】
次に、図８に示すように、レンジ期間Ｐ２は（１）式、通信期間Ｐ３は（３）式の変調信号で変調を用いて変調する。通信期間Ｐ３の変調は、ほかの手法でもよい。
【００３２】
次に、受信信号のうちレンジ期間の信号を圧縮する手法について定式化して説明する。これは、入力信号とパルス圧縮用信号の相関処理であり、これを周波数領域で行う場合について定式化すると次の通りである。
【００３３】
まず、入力信号は、次式となる。  
【００３４】
【数４】
2019095391.tif 000006
【００３５】
次に、参照信号は次の通りである。  
【００３６】
【数５】
2019095391.tif 000007
【００３７】
周波数軸の信号で乗算すると、次式となる。  
【００３８】
【数６】
2019095391.tif 000008
【００３９】
パルス圧縮出力は逆ＦＦＴにより算出できる。  
【００４０】
【数７】
2019095391.tif 000009
【００４１】
以上により、図９に示すように、レンジ抽出ができる。レンジ抽出ができると、時間軸で通信期間Ｐ３の抽出ができるため、次に通信期間Ｐ３の復調を行う。このためには、通信期間Ｐ３の信号scom（ｔ）のＦＦＴを行う。
【００４２】
【数８】
2019095391.tif 000010
【００４３】
通信期間Ｐ３を高精度に抽出し、その期間のみによる復調処理ができるため、通信復調の精度が高まることになる。この信号の位相情報を抽出することにより、ビットを抽出して復調できる。
【００４４】
なお、図９では、簡単のために１ＰＲＩ分の処理について述べたが、Ｎヒットの複数ＰＲＩの場合には、図１０に示すように、Ｐ２列の相関処理（Ｒ１５）とslow-time軸のＦＦＴ（Ｒ１６）を行い、ドップラ軸の最大バンク信号に対して、距離を抽出（Ｒ１７）すれば、高いＳＮで通信期間Ｐ３を抽出でき、復調精度も高まる。
【００４５】
次に、通信復調に用いた目標以外の目標のレンジングを行うことを検討する。この場合は、小目標も含まれるため、レンジングのための相関処理時には、通信期間Ｐ３＋レンジ期間Ｐ２の相関処理を行う方がＳＮが高くなる。このため、参照信号を連結する。
【００４６】
通信期間Ｐ３の参照信号は、復調信号が既知となるため、（３）式のScomと同じである。これと、（５）式のレンジ期間の参照信号Srefを連結すればよい。
【００４７】
【数９】
2019095391.tif 000011
【００４８】
このScom_rを用いて相関処理を行う。
【００４９】
まず、入力信号は、次式となる。  
【００５０】
【数１０】
2019095391.tif 000012
【００５１】
次に、参照信号は次の通りである。  
【００５２】
【数１１】
2019095391.tif 000013
【００５３】
周波数軸の信号で乗算すると、次式となる。  
【００５４】
【数１２】
2019095391.tif 000014
【００５５】
パルス圧縮出力は逆ＦＦＴにより算出できる。  
【００５６】
【数１３】
2019095391.tif 000015
【００５７】
以上により、信号対雑音電力比の低い反射点についても、通信期間Ｐ３＋レンジ期間Ｐ２全体のパルスを用いた相関処理を行い、Ｎヒットの複数ＰＲＩがある場合には、図１１に示すようにslow-time軸のＦＦＴ（Ｒ２５）を行い、ＣＦＡＲ等により検出（Ｒ２６）した後、高いＳＮでレンジ抽出（Ｒ２７）できる。本実施形態は、通信期間Ｐ３のパルスの受信期間を、レンジ期間Ｐ２のパルスの相関処理結果を用いて抽出し、高いＳＮで通信復調して送信装置Ｔと受信装置Ｒを同期同調するともに、通信期間Ｐ３を含めた相関処理により、目標Ｔの距離を高いＳＮで抽出することに特徴を有する。
【００５８】
なお、通信期間Ｐ３を含めずレンジ期間Ｐ２のみの相関処理で、目標Ｔが検出できれば、レンジ期間Ｐ２のみの相関処理によるレンジを出力してもよい。
【００５９】
以上により、図５に示すように、送信〜受信の距離と、（１３）式で検出した信号の測角処理より、受信装置からみたAZ、EL角を観測できるのので、目標の位置を同定することができる。
【００６０】
また、図６に示すように送信装置を３台とすると、送信１〜受信、送信２〜受信、送信３〜受信までの各々の距離として、Ｒ１、Ｒ２、Ｒ３を得ることができる。この距離を用いて、図６に示すように、目標Ｔの位置（ｘ、ｙ、ｚ）を算出する。この手法としては、Ｒ１の球面とＲ２及びＲ３の楕円球面の交点となる。その中で、受信装置Ｒにより観測したＡＺ角、ＥＬ角方向の３次元の位置を中心に、所定の範囲内を目標存在領域として、その中の交点を算出する。解で算出できない場合は、目標存在領域内の点を（ｘ、ｙ、ｚ）の格子点に分割し、各々の点で次式の値が最小となる点（ｘ、ｙ、ｚ）を算出する。
【００６１】
【数１４】
2019095391.tif 000016
【００６２】
なお、受信装置に送受信機能がある場合には、送信装置は２台により、同様の方法での交線の中点を目標の３次元の観測位置として出力することができる。
【００６３】
以上のように、第１の実施形態に係るレーダシステムでは、送信装置はパルス圧縮用の送信パルスτ（連続波含む）を通信期間Ｐ３とレンジ期間Ｐ２に分割し、通信期間のパルスに送信に関する通信情報を重畳して送信し、受信装置は送信装置から送信される直接波または目標からの反射波を受信し、その受信信号をレンジ期間の参照信号で圧縮して、ＳＮの高い反射点を用いて、通信期間を抽出して相関処理することにより、通信情報（送信位置、送信時刻、送信周波数、送信ビーム指向方向等）を抽出し、信号対雑音電力比の低い他の反射点についても、レンジ期間のみか、または、通信情報を用いて、通信期間＋レンジ期間のパルス全体で相関処理して、目標の距離を出力するようにしている。
【００６４】
上記の構成によれば、レンジ期間Ｐ２のパルスで通信期間Ｐ３を抽出し、復調して送信装置の位置、送信タイミング、送信周波数、送信ビーム指向方向等を通信情報から抽出することにより、送信装置と受信装置を同期同調し、補正した参照信号を用いて、通信期間Ｐ３、レンジ期間Ｐ２の長いパルス列で相関処理することで、高いＳＮで目標の距離を出力することができる。
【００６５】
（第２の実施形態）
第１の実施形態では、パルスを通信期間とレンジ期間に分割する手法について述べた。しかしながら、ＨＰＲＦ（High Pulse Repetition Frequency）のような短パルスでは、パルスを分割することが困難である。そこで、第２の実施形態では、パルス列を分割する手法に適用する例について述べる。
【００６６】
本実施形態に係るレーダシステムの送信装置Ｔｉの構成を図１２に、受信装置Ｒの構成を図１３に示す。
【００６７】
図１２に示す送信装置Ｔｉにおいて、図２と異なる点は、チャープ信号生成器Ｔ１６に代わって基準信号生成器Ｔ２０を用い、制御器Ｔ１５に代わって時分割処理を行う制御器Ｔ１９を用いたことにある。
【００６８】
また、図１３に示す受信装置Ｒにおいて、図３と異なる点は、Ｐ２抽出器Ｒ１４及び相関処理器（パルス圧縮）Ｒ１５に代わってＰ１選定器Ｒ２９を用い、検出器Ｒ１７の後段にドップラ抽出器Ｒ３０を追加し、Ｐ３抽出器Ｒ１９に代わってドップラ参照信号補正器Ｒ３１を用い、Ｐ２＋Ｐ３選定器Ｒ２２に代わってＰ２選定器Ｒ３２、ＰＲＩ列並べ替え器Ｒ３３、ドップラ参照信号補正器Ｒ３４を用いたことにある。
【００６９】
図１２に示す送信装置Ｔｉにおいて、送信パルス列は、図１４に示すように、ドップラ期間Ｐ１、レンジ期間Ｐ２、通信期間Ｐ３に分割される。このＰ１、Ｐ２、Ｐ３のそれぞれのパルス列は、ＰＲＦ（Pulse Repetition Frequency）が同じでも異なっていてもよい。パルス送信変調は、図１５に示すように、通信情報で符号変調したＢＰＳＫ（非特許文献４）をＦＦＴして、サンプリングした変調信号でパルス間を変調する。
【００７０】
受信は、図１６に示すように、ドップラ期間Ｐ１のＰＲＦ１のＦＦＴにより、ＳＮの高い目標で信号検出する。ドップラ抽出できれば、レンジ期間Ｐ２の参照信号を補正して、レンジ圧縮し、レンジ抽出を行う。レンジ抽出ができれば、通信期間Ｐ３を時間軸で抽出できるため、図１７に示すように通信期間のパルス列Ｐ３のＩＦＦＴにより、復調し、通信情報を得る。通信情報として、送信位置、送信周波数、送信タイミング、送信ビーム指向方向等のマルチスタティックレーダとしての使い方は、第１の実施形態と同様である。
【００７１】
復調信号（通信情報）がわかれば、第１の実施形態と同様に通信期間Ｐ３の補正信号がわかる。レンジ抽出のためには、補正した参照信号と入力信号の相関処理を行えばよい。相関処理としては、slow-time軸×fast-time軸の全セルを用いた相関処理手法もあるが、ＳＮ（信号対雑音電力比）向上を図るために、レンジセル毎のslow-time軸の相関処理を行う。この場合、Ｐ１列、Ｐ２列、Ｐ３列のＰＲＦが同じ場合には、図１８に示すように、方式（１）ではＰ２列、方式（２）ではＰ１＋Ｐ２列、または方式（３）ではＰ１＋Ｐ２＋Ｐ３列による相関処理がある。まず、パルス列の相関処理の一般形を定式化して説明する。
【００７２】
まず、入力信号は、次式となる。  
【００７３】
【数１５】
2019095391.tif 000017
【００７４】
次に、参照信号は次の通りである。
【００７５】
【数１６】
2019095391.tif 000018
【００７６】
周波数軸の信号で乗算すると、次式となる。
【００７７】
【数１７】
2019095391.tif 000019
【００７８】
相関処理は逆ＦＦＴにより算出できる。  
【００７９】
【数１８】
2019095391.tif 000020
【００８０】
以上により、fast-time軸のレンジセル（レンジセル数Cellfast）毎でslow-time軸（セル数Cellslow）の相関処理出力が得られる。全体時間は、Cellfast×Cellslowのセル数で、１セルの時間はfast-time軸のサンプリング時間の時間である。この全体時間の中で図１８に示すように、レンジを抽出することができる。
【００８１】
次に、各方式（１）〜（３）は、（１６）式の参照信号sref0（tslow）が異なるのみである。
【００８２】
方式（１）のＰ２列のみの場合は、レンジ期間の変調信号をドップラにより補正した参照信号を用いる。
【００８３】
【数１９】
2019095391.tif 000021
【００８４】
方式（２）のＰ１＋Ｐ２列の場合は、ドップラ期間＋レンジ期間の補正を行う。
【００８５】
【数２０】
2019095391.tif 000022
【００８６】
方式（３）のＰ１＋Ｐ２＋Ｐ３列の場合は、方式（２）にＰ１列を付加する。
【００８７】
【数２１】
2019095391.tif 000023
【００８８】
これらの参照信号を用いて、全体の相関処理を行う。これにより、信号対雑音電力比の低い反射点についても、高いＳＮで小目標を検出でき、距離と速度を出力できる。なお、Ｐ１列、Ｐ２列及びＰ３列のＰＲＦにおいて、異なるＰＲＦがある場合には、少なくともレンジ期間のＰ２列を含む同じＰＲＦの信号列を用いて、レンジセル毎のslow-time軸における相関処理を行えばよい。
【００８９】
以上述べたように、本実施形態では、通信期間のパルスの受信期間を、ドップラ期間とレンジ期間のパルスの相関処理結果を用いて抽出し、高いＳＮで通信復調して送信装置と受信装置を同期同調するともに、信号対雑音電力比の低い他の反射点についても、通信期間、ドップラ期間、レンジ期間を含めた相関処理により、目標の距離を高いＳＮで抽出することに特徴を有する。
【００９０】
すなわち、第２の実施形態に係るレーダシステムでは、送信装置は、パルスを繰り返すレーダのパルス列を通信期間Ｐ３とドップラ期間Ｐ１とレンジ期間Ｐ２に時分割し、通信期間のパルスに送信に関する通信情報を重畳して送信し、受信装置は送信装置から送信される直接波または目標からの反射波を受信し、その受信信号を、ドップラ期間Ｐ１のＦＦＴ結果のうち、ＳＮ比の高い反射点を用いてレンジ期間の参照信号を補正し、補正した参照信号によりレンジ期間Ｐ２を相関処理してレンジを抽出し、通信期間Ｐ３の受信信号を抽出してドップラで補正した参照信号により相関処理することによって通信情報（送信位置、送信時刻、送信周波数、送信ビーム指向方向等）を抽出する。信号対雑音電力比の低い他の反射点についても、通信期間Ｐ３は、通信情報とドップラを用いて参照信号を補正し、レンジ期間Ｐ２はドップラを用いて参照信号を補正し、通信期間Ｐ３、ドップラ期間Ｐ１、レンジ期間Ｐ２のうち、少なくともレンジ期間Ｐ２を含むパルス列で用いて相関処理して、目標の距離を出力する。
【００９１】
上記の構成によれば、ドップラ期間Ｐ１とレンジ期間Ｐ２のパルス列で通信期間Ｐ３を抽出し、復調して送信装置の位置、送信タイミング、送信周波数を通信情報から抽出することにより、送信装置と受信装置を同期同調し、補正した参照信号を用いて、通信期間Ｐ３、ドップラ期間Ｐ１、レンジ期間Ｐ２の長いパルス列で相関処理することで、高いＳＮで目標の距離を出力できる。
【００９２】
（第３の実施形態）
第２の実施形態では、パルス列をドップラ期間、通信期間及びレンジ期間に分割する手法について述べた。この場合、時分割であるため、パルス列が長くなり、観測時間が長くなる。そこで、第３の実施形態では、これを避けるために、Ｐ１、Ｐ２、Ｐ３列を異なるＰＲＦとして、同時送信する手法について述べる。
【００９３】
本実施形態に係るレーダシステムの送信装置Ｔｉの構成を図１９に、受信装置Ｒの構成を図２０に示す。
【００９４】
図１９に示す送信装置Ｔｉにおいて、図１２と異なる点は、時分割処理を行う制御器Ｔ１９に代わって統合処理を行う制御器Ｔ２１を用いたことにある。
【００９５】
また、図２０に示す受信装置Ｒにおいて、図１３と異なる点は、通信情報参照信号補正器Ｒ２３を削除したことにある。
【００９６】
図１９に示す送信装置Ｔｉにおいて、送信パルス列は、図２１に示すように、通信期間Ｐ３、ドップラ期間Ｐ１、レンジ期間Ｐ２があり、これを同時に送信するように合成する（Ｔ２１）。送信変調は、図２２に示すように、通信情報で符号変調したＢＰＳＫ（非特許文献４）をＦＦＴして、サンプリングした変調信号でパルス間を変調する。
【００９７】
図２０に示す受信装置Ｒでは、図２３に示すように、各信号列Ｐ１、Ｐ２、Ｐ３のＰＲＦが異なる場合は、ＰＲＦを用いて信号を分離できる。また、通信期間Ｐ３とドップラ期間Ｐ１のＰＲＦが同じ場合は、図２４に示すように、Ｐ１＋Ｐ３とＰ２の２通りに分離してもよい。また、通信期間Ｐ３とレンジ期間Ｐ２のＰＲＦが同じ場合は、Ｐ１とＰ２＋Ｐ３の２通りに分離してもよい。
【００９８】
分離した信号列を用いて、図２５に示すように、ドップラ期間Ｐ１のＰＲＦ１のＦＦＴにより、ＳＮの高い目標で信号検出する。Ｐ１＋Ｐ３の場合は、通信期間は通信変調信号によりパルス間で位相が異なるため加算されないが、ドップラ期間を長くすれば、高いＳＮでドップラを抽出できる。ドップラ抽出ができれば、レンジ期間の参照信号を補正して、レンジ圧縮し、レンジ抽出を行う、レンジ抽出ができれば、時間軸で、通信期間Ｐ３を時間軸で抽出できるため、通信期間のパルス列Ｐ３のＩＦＦＴにより、復調し、通信情報を得る。通信情報として、送信位置、送信周波数、送信タイミング等のマルチスタティックレーダとしての使い方は、第１の実施形態と同様である。
【００９９】
次に、レンジ抽出のためには、補正した参照信号と入力信号の相関処理を行えばよい。この場合、第２の実施形態と同様に、レンジセル毎にslow-time軸で相関処理を行うためには、同じＰＲＦ列が必要である。このため、図２６に示すように、Ｐ２列のみを用いて相関処理を行う。なお、通信期間＋レンジ期間を同じＰＲＦとしたＰ２＋Ｐ３の場合は、信号対雑音電力比の低い他の反射点についても、通信情報を復調した信号とドップラ期間のドップラ信号を用いて補正した通信期間Ｐ３の参照信号（第１の参照信号）と、ドップラ信号により補正したレンジ期間Ｐ２の参照信号（第２の参照信号）とを合わせて全体の参照信号として、Ｐ２＋Ｐ３の相関処理を行えばよい。
【０１００】
これらの参照信号のドップラによる補正を含めた生成方法は、各パルス列の時間軸が図２３の上図であり、時分割パルス列である第２の実施形態（図１４）と異なる。このため、図２３に合わせたパルス列の時間軸にする必要がある以外は、第２の実施形態と同様であるため、ここではその説明を割愛する。これらの参照信号を用いて、全体の相関処理を行う。これにより、高いＳＮで小目標を検出でき、距離と速度を出力できる。
【０１０１】
以上述べたように、本実施形態では、通信期間のパルスの受信期間を、ドップラ期間とレンジ期間のパルスの相関処理結果を用いて抽出し、高いＳＮで通信復調して送信装置と受信装置を同期同調するともに、信号対雑音電力比の低い他の反射点についても、通信期間、ドップラ期間、レンジ期間を含めた相関処理により、目標の距離を高いＳＮで抽出することに特徴を有する。
【０１０２】
すなわち、本実施形態に係るレーダシステムでは、送信装置は、パルスを繰り返すレーダのパルス列を通信期間Ｐ３とドップラ期間Ｐ１とレンジ期間Ｐ２の３通り（Ｐ３，Ｐ１，Ｐ２）、または通信期間Ｐ３及びドップラ期間Ｐ１とレンジ期間Ｐ２の２通り（Ｐ３＋Ｐ１，Ｐ２）、またはドップラ期間Ｐ１と通信期間Ｐ３及びレンジ期間Ｐ２の２通り（Ｐ１，Ｐ３＋Ｐ２）のＰＲＦのパルス列とし、通信期間のパルスに送信に関する通信情報を重畳して、通信期間Ｐ３、ドップラ期間Ｐ１、レンジ期間Ｐ２のパルス列を混合して送信し、受信装置は送信装置から送信される直接波または目標からの反射波を受信し、その受信信号を、ドップラ期間Ｐ１のＦＦＴ結果のうち、ＳＮの高い反射点のドップラを用いてレンジ期間Ｐ２の参照信号を補正し、補正した参照信号によりレンジ期間Ｐ２を相関処理してレンジを抽出し、通信期間Ｐ３の受信信号を抽出して、ドップラで補正した参照信号で相関処理することによって通信情報（送信位置、送信時刻、送信周波数、送信ビーム指向方向等）を抽出し、信号対雑音電力比の低い他の反射点についても、レンジ期間Ｐ２において参照信号をドップラで補正したパルス列を用いて相関処理して、目標の距離を出力する。
【０１０３】
このように、本実施形態によれば、３通りのＰＲＦのパルス列を混合して同時に送信することで、時分割の場合に比べて送信時間を短縮化するとともに、ＬＰＩ（Low Probability of intercept）性も高まり、ドップラ期間Ｐ１とレンジ期間Ｐ２のパルス列で通信期間を抽出し、復調して送信装置の位置、送信タイミング、送信周波数、送信ビーム指向方向等を通信情報から抽出することにより、送信装置と受信装置を同期同調することができる。
【０１０４】
なお、本実施形態は上記実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。また、上記実施形態に開示されている複数の構成要素の適宜な組み合わせにより、種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。更に、異なる実施形態にわたる構成要素を適宜組み合わせてもよい。
【符号の説明】
【０１０５】
Ｔ１〜ＴＮ…送信装置、Ｒ…受信装置、Ｔ…目標、
Ｔ１１…送信アンテナ装置、Ｔ１２…パルス変調器、Ｔ１３…周波数変換器、Ｔ１４…変調器、Ｔ１５…制御器、Ｔ１６…チャープ信号生成器、Ｔ１７…ＩＦＦＴ（Inverse Fast Fourier Transform）処理器、Ｔ１８…通信変調信号生成器、Ｔ１９…制御器、Ｔ２０…基準信号生成器、Ｔ２１…制御器、
Ｒ１１…受信アンテナ装置、Ｒ１２…周波数変換器、Ｒ１３…ＡＤ（Analog-Digital）変換器、Ｒ１４…Ｐ２抽出器、Ｒ１５…相関処理パルス圧縮器、Ｒ１６…ＦＦＴ（Slow-time）処理器、Ｒ１７…検出器、Ｒ１８…Ｐ３範囲選定器、Ｒ１９…Ｐ３抽出器、Ｒ２０…相関処理器（Slow-time）、Ｒ２１…復調器、Ｒ２２…Ｐ２＋Ｐ３選定器、Ｒ２３…通信情報参照信号補正器、Ｒ２４…相関処理器、Ｒ２５…ＦＦＴ（Slow-time）処理器、Ｒ２６…検出器、Ｒ２７…レンジ抽出器、Ｒ２８…出力処理器、Ｒ２９…Ｐ１選定器、Ｒ３０…ドップラ抽出器、Ｒ３１…ドップラ参照信号補正器、Ｒ３２…Ｐ２選定器、Ｒ３３…ＰＲＩ列並べ替え器、Ｒ３４…ドップラ参照信号補正器。
【図１】
2019095391.tif 000024
【図２】
2019095391.tif 000025
【図３】
2019095391.tif 000026
【図４】
2019095391.tif 000027
【図５】
2019095391.tif 000028
【図６】
2019095391.tif 000029
【図７】
2019095391.tif 000030
【図８】
2019095391.tif 000031
【図９】
2019095391.tif 000032
【図１０】
2019095391.tif 000033
【図１１】
2019095391.tif 000034
【図１２】
2019095391.tif 000035
【図１３】
2019095391.tif 000036
【図１４】
2019095391.tif 000037
【図１５】
2019095391.tif 000038
【図１６】
2019095391.tif 000039
【図１７】
2019095391.tif 000040
【図１８】
2019095391.tif 000041
【図１９】
2019095391.tif 000042
【図２０】
2019095391.tif 000043
【図２１】
2019095391.tif 000044
【図２２】
2019095391.tif 000045
【図２３】
2019095391.tif 000046
【図２４】
2019095391.tif 000047
【図２５】
2019095391.tif 000048
【図２６】
2019095391.tif 000049