JP2017053721A - レーダシステム、信号処理システム、信号処理方法及びプログラム - Google Patents
レーダシステム、信号処理システム、信号処理方法及びプログラム Download PDFInfo
- Publication number
- JP2017053721A JP2017053721A JP2015177821A JP2015177821A JP2017053721A JP 2017053721 A JP2017053721 A JP 2017053721A JP 2015177821 A JP2015177821 A JP 2015177821A JP 2015177821 A JP2015177821 A JP 2015177821A JP 2017053721 A JP2017053721 A JP 2017053721A
- Authority
- JP
- Japan
- Prior art keywords
- signal
- detection target
- reception
- radar
- receiving station
- 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.)
- Granted
Links
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
【課題】検出対象までの距離の検出精度を高めることができるレーダシステム、信号処理システム、信号処理方法及びプログラムを提供することである。【解決手段】レーダシステムは、受信局と、補正部と、推定部とを持つ。受信局は、送信局から送信された電波が検出対象において反射されて、前記反射された電波の周波数成分を含む受信信号を受信する。補正部は、前記送信局が送信した電波の周波数と、前記受信信号に含まれる前記反射された電波の周波数との差として生じたドップラ遷移を打ち消す方向に、前記受信信号の周波数成分を補正する。推定部は、移動する前記検出対象から前記受信局までの距離を、前記補正部により補正された前記受信信号の周波数成分に基づいて推定する。【選択図】図４
Description
本発明の実施形態は、レーダシステム、信号処理システム、信号処理方法及びプログラムに関する。
従来、レーダシステムにおいて、送信時にレーダ波に対して変調処理を施し、受信時に復調処理を行うことでレンジング処理の分解能向上を図るパルスレーダが知られている。ただし、検出対象が移動する場合、変調されたレーダ波は、検出対象の相対速度に対応するドップラ遷移の影響を受ける。ドップラ遷移の影響を受けたレーダ波を受信レーダが復調処理をした際に、復調後の信号にドップラ遷移の影響が表れて、レンジング処理により得られた位置と実際の検出対象の位置との間に誤差が生じる。
また、送信波の反射信号を複数の受信レーダでそれぞれ受信するマルチスタティックレーダが知られている。マルチスタティックレーダにおける各受信レーダは、同一の検出対象において反射されたレーダ波であっても、検出対象の相対速度が互いに異なることにより、ドップラ遷移の影響度合いが異なるレーダ波をそれぞれ受信する。そのため、各受信レーダにおいて共通の復調処理を実施した場合、レンジング処理により得られた検出対象の位置と実際の検出対象の位置との間に、受信レーダごとに互いに異なる誤差が生じることにより、検出対象までの距離を所望の検出精度で得ることは容易ではない。
また、送信波の反射信号を複数の受信レーダでそれぞれ受信するマルチスタティックレーダが知られている。マルチスタティックレーダにおける各受信レーダは、同一の検出対象において反射されたレーダ波であっても、検出対象の相対速度が互いに異なることにより、ドップラ遷移の影響度合いが異なるレーダ波をそれぞれ受信する。そのため、各受信レーダにおいて共通の復調処理を実施した場合、レンジング処理により得られた検出対象の位置と実際の検出対象の位置との間に、受信レーダごとに互いに異なる誤差が生じることにより、検出対象までの距離を所望の検出精度で得ることは容易ではない。
M. I. Skolnik, "Radar Handbook," Third edition, pp. 23.1-23.13, McGraw-Hill, New York,
N. J. Wills, "Bistatic Radar," Second edition, Silver Spring, MD: Technology Service Corp., 1995. Corrected and republished by Raleigh, NC: SciTech Publishing, Inc., 2005
本発明が解決しようとする課題は、検出対象までの距離の検出精度を高めることができるレーダシステム、信号処理システム、信号処理方法及びプログラムを提供することである。
実施形態のレーダシステムは、受信局と、補正部と、推定部とを持つ。受信局は、送信局から送信された電波が検出対象において反射されて、前記反射された電波の周波数成分を含む受信信号を受信する。補正部は、前記送信局が送信した電波の周波数と、前記受信信号に含まれる前記反射された電波の周波数との差として生じたドップラ遷移を打ち消す方向に、前記受信信号の周波数成分を補正する。推定部は、移動する前記検出対象から前記受信局までの距離を、前記補正部により補正された前記受信信号の周波数成分に基づいて推定する。
以下、実施形態のレーダシステムを、図面を参照して説明する。以下、前述の構成と同じ構成には同じ符号を附す。
（第１の実施形態）
（レーダシステムの概要）
図１は、第１の実施形態に係るレーダシステム１の概略を示す図である。レーダシステム１は、例えば、送信レーダ（以下、「ＴＸ」という。）１０と、Ｎ局の受信レーダ（以下、「ＲＸ」という）２０−１、２０−２、…、２０−Ｎと、演算処理部５０とを含む。以下、いずれの受信レーダであるかを区別しないときは、単にＲＸ２０という。レーダシステム１は、ＴＸ１０により送信された電波が検出対象３０によって反射された反射波を、ＴＸ１０とは異なる位置に配置された複数のＲＸ２０で受信することで、検出対象３０の位置を特定するマルチスタティックレーダである。この場合、Ｎは、例えば３以上の整数である。なお、レーダシステム１は、ＲＸ２０を２つ備えるバイスタティックレーダであってもよい。
（レーダシステムの概要）
図１は、第１の実施形態に係るレーダシステム１の概略を示す図である。レーダシステム１は、例えば、送信レーダ（以下、「ＴＸ」という。）１０と、Ｎ局の受信レーダ（以下、「ＲＸ」という）２０−１、２０−２、…、２０−Ｎと、演算処理部５０とを含む。以下、いずれの受信レーダであるかを区別しないときは、単にＲＸ２０という。レーダシステム１は、ＴＸ１０により送信された電波が検出対象３０によって反射された反射波を、ＴＸ１０とは異なる位置に配置された複数のＲＸ２０で受信することで、検出対象３０の位置を特定するマルチスタティックレーダである。この場合、Ｎは、例えば３以上の整数である。なお、レーダシステム１は、ＲＸ２０を２つ備えるバイスタティックレーダであってもよい。
ＲＸ２０のそれぞれは、互いに距離が隔てて配置される。ＲＸ２０同士の距離は、ＲＸ２０のそれぞれからの検出対象３０に対する測角方向が互いに異なる程度に決定される。また、少なくとも複数のＲＸ２０のそれぞれは、ＴＸ１０に対する距離が隔てて配置される。ＴＸ１０とＲＸ２０との距離は、上記の複数のＲＸ２０のそれぞれからの検出対象３０に対する測角方向が、ＴＸ１０からの検出対象３０の方向と異なる程度に決定される。
ＴＸ１０は、例えば、変調された信号を送信波として出力する。ＴＸ１０が送信した送信波は、検出対象３０において反射され、その反射波がＮ個のＲＸ２０によりそれぞれ受信される。
（マルチスタティックレーダの方式）
図２と図３とを参照して、実施形態におけるマルチスタティックレーダの方式について説明する。図２と図３は、実施形態において採用し得るマルチスタティックレーダを例示した図である。マルチスタティックレーダは、検出対象を検出する処理の違いにより、次の２つの方式に大別される。
図２と図３とを参照して、実施形態におけるマルチスタティックレーダの方式について説明する。図２と図３は、実施形態において採用し得るマルチスタティックレーダを例示した図である。マルチスタティックレーダは、検出対象を検出する処理の違いにより、次の２つの方式に大別される。
図２に示す第１の方式は、各受信レーダ（例えばＲＸ２０）で検出対象３０の検出処理は行わず、得られたＲａｗデータのレベルの情報を中央装置（例えば演算処理部５０）に伝送し、中央装置側で検出対象３０の検出の有無を判定する方式である。このような方式を適用したものは、Centralized型と呼ばれている。
図３に示す第２の方式は、各受信レーダ（例えば後述のＲＸ２０Ａ）で検出対象３０の検出処理を行い、検出結果であるＰｌｏｔレベルの情報を中央装置（例えば後述の演算処理部５０Ａ）に伝送し、中央装置が各受信レーダの検出結果を基にして改めて総合的な検出判定行う方式である。このような方式を適用したものは、Decentralized型と呼ばれている。
第１の実施形態のレーダシステム１は、上記の第１の方式のCentralized型に分類される。ＲＸ２０は、それぞれ受信した反射波に基づいた信号処理を行った後のデータを演算処理部５０に提供する。演算処理部５０は、ＲＸ２０のそれぞれから提供されたデータを解析して、検出対象３０を検出するとともに、検出対象３０の位置を推定する。
（レーダシステムの構成）
図４を参照して、本実施形態におけるレーダシステム１の構成について説明する。同図は、本実施形態におけるレーダシステム１の構成図である。
レーダシステム１（信号処理システム）は、例えば、複数のＲＸ２０と、演算処理部５０とを備える。レーダシステム１は、さらに、ＴＸ１０を備えていてもよい。
図４を参照して、本実施形態におけるレーダシステム１の構成について説明する。同図は、本実施形態におけるレーダシステム１の構成図である。
レーダシステム１（信号処理システム）は、例えば、複数のＲＸ２０と、演算処理部５０とを備える。レーダシステム１は、さらに、ＴＸ１０を備えていてもよい。
［ＴＸ］
ＴＸ１０は、例えば、所定の規則に従って変調した変調波をパルス波（電波）にして、アンテナ１１から送出する。変調波の変調方式は、周波数変調であってもよく、予め定められた符号則に従ったパルス変調であってもよい。なお、変調方式の規則は、各ＲＸ２０と共有される。
ＴＸ１０は、例えば、所定の規則に従って変調した変調波をパルス波（電波）にして、アンテナ１１から送出する。変調波の変調方式は、周波数変調であってもよく、予め定められた符号則に従ったパルス変調であってもよい。なお、変調方式の規則は、各ＲＸ２０と共有される。
上記のパルス波は、独立したパルスとして送出されてもよく、或いは、所定の期間に断続的に複数送出され、予め定められた送出タイミングに従って繰り返されてもよい。以下の説明では、所定の期間に断続的に複数送出されるパルス波の場合を例示して説明する。
［ＲＸ］
ＲＸ２０は、ＴＸ１０から送信された変調波（電波）が検出対象３０において反射されて、反射された電波をアンテナ２１によって受信する。ＲＸ２０は、アンテナ２１によって受信した電波から受信信号（信号ＲＳ１）を取り出し、ＴＸ１０が送信した変調波（電波）の周波数と、信号ＲＳ１に含まれる変調波の周波数との差として生じたドップラ遷移を打ち消す方向に、信号ＲＳ１の周波数成分を補正する。
ＲＸ２０は、ＴＸ１０から送信された変調波（電波）が検出対象３０において反射されて、反射された電波をアンテナ２１によって受信する。ＲＸ２０は、アンテナ２１によって受信した電波から受信信号（信号ＲＳ１）を取り出し、ＴＸ１０が送信した変調波（電波）の周波数と、信号ＲＳ１に含まれる変調波の周波数との差として生じたドップラ遷移を打ち消す方向に、信号ＲＳ１の周波数成分を補正する。
ＲＸ２０は、例えば、信号変換部２２と、復調処理部２３と、パルス積分処理部２６と、通信部２７とを備える。なお、図４ではＲＸ２０を代表してＲＸ２０−１の構成についてのみ記載し、その他のＲＸの構成については省略している。
信号変換部２２は、信号ＲＳ１に対して周波数変換とアナログ・ディジタル変換とを実施し、信号ＲＳ１に含まれる所望の信号成分を含む信号ＲＳ２をディジタル信号として出力する。信号変換部２２は、さらに、ディジタル信号に対して所定の重み係数を乗じることにより、所望の受信ビームを形成してもよい。
信号変換部２２は、信号ＲＳ１に対して周波数変換とアナログ・ディジタル変換とを実施し、信号ＲＳ１に含まれる所望の信号成分を含む信号ＲＳ２をディジタル信号として出力する。信号変換部２２は、さらに、ディジタル信号に対して所定の重み係数を乗じることにより、所望の受信ビームを形成してもよい。
復調処理部２３（補正部）は、信号変換部２２によって変換された信号ＲＳ２の復調処理を実施する。復調処理部２３は、検出対象３０の移動にともなって検出対象３０における反射波に生じるドップラ遷移であって、ＴＸ１０が送信した変調波の周波数と、信号ＲＳ１に含まれる変調波の周波数との差として生じたドップラ遷移を打ち消す方向に、信号ＲＳ１の周波数成分を補正する。
復調処理部２３は、例えば、復調係数生成部２４と変換部２５とを備える。
復調係数生成部２４は、例えば、復調処理の系統数と同じ個数の復調係数を生成する。復調処理の系統とは、上記のドップラ遷移を打ち消すことについて複数の結果を得るために、予め周波数成分の補正量の傾向を変えて用意した構成のことである。復調係数生成部２４は、ドップラ遷移を打ち消すように調整された復調係数を、復調処理の系統に対応させて生成する。
変換部２５は、復調係数生成部２４によって生成された復調係数に基づいて、信号ＲＳ２の周波数成分を系統ごとに変換する。
変換部２５は、復調処理を実施する際に、復調係数生成部２４によって生成された復調係数に基づいてパルス圧縮処理を実施する。上記のように、ドップラ遷移を打ち消すように調整された復調係数に基づいて復調処理を実施することにより、ドップラ遷移の影響を低減した復調処理を実施する。
復調係数生成部２４は、例えば、復調処理の系統数と同じ個数の復調係数を生成する。復調処理の系統とは、上記のドップラ遷移を打ち消すことについて複数の結果を得るために、予め周波数成分の補正量の傾向を変えて用意した構成のことである。復調係数生成部２４は、ドップラ遷移を打ち消すように調整された復調係数を、復調処理の系統に対応させて生成する。
変換部２５は、復調係数生成部２４によって生成された復調係数に基づいて、信号ＲＳ２の周波数成分を系統ごとに変換する。
変換部２５は、復調処理を実施する際に、復調係数生成部２４によって生成された復調係数に基づいてパルス圧縮処理を実施する。上記のように、ドップラ遷移を打ち消すように調整された復調係数に基づいて復調処理を実施することにより、ドップラ遷移の影響を低減した復調処理を実施する。
なお、復調処理部２３は、複数の系統に分けて、信号変換部２２によって変換された信号ＲＳ２の復調処理を実施する。それぞれの系統に対応させて、復調係数生成部２４と変換部２５とを組にして複数設けておき系統ごとの処理を並列処理してもよく、１対の復調係数生成部２４と変換部２５とを系統ごとに時分割処理をしてもよい。
パルス積分処理部２６は、復調係数に基づいて復調された後の信号ＲＳ３を、復調処理の系統ごとにパルス積分する。
通信部２７は、パルス積分処理部２６によるパルス積分の結果(信号ＲＳ４)を演算処理部５０宛に通知する。
ＲＸ２０−２からＲＸ２０−Ｎは、例えば、ＲＸ２０−１と同様の構成を備える。それらの説明はＲＸ２０−１の説明を参照する。
演算処理部５０は、ＲＸ２０のそれぞれから通知されたパルス積分の結果（信号ＲＳ４）に基づいて、検出対象３０の位置を推定する。演算処理部５０は、例えば、通信部５１と、選択部５２と、インコヒーレント積分処理部５３と、閾値判定処理部５４と、推定処理部５５と、出力処理部５６とを備える。
通信部５１は、ＲＸ２０からのデータを受信する。
選択部５２は、複数のＲＸ２０によってそれぞれ受信され、復調処理部２３によりそれぞれ補正された後の信号ＲＳ４のうちから、受信信号の信号強度に基づいて受信信号をＲＸ２０ごとに選択し、選択された結果を信号ＲＳ５とする。
インコヒーレント積分処理部５３（積分処理部）は、ＲＸ２０ごとに選択された結果の信号ＲＳ５に対してインコヒーレント積分処理を実施して、インコヒーレント積分処理の結果を信号ＲＳ６とする。
閾値判定処理部５４は、予め定められた閾値を判定条件にして、信号ＲＳ６の判定を実施する。当該判定により、検出対象３０として判定し得る確率が比較的高いものを検出対象３０の候補として判定する。
推定処理部５５は、インコヒーレント積分処理部５３による演算処理（積分処理）された後の信号ＲＳ６に基づいて、判定処理により選択された検出対象３０の位置を推定する。
出力処理部５６は、推定された検出対象３０の位置を出力する。
通信部５１は、ＲＸ２０からのデータを受信する。
選択部５２は、複数のＲＸ２０によってそれぞれ受信され、復調処理部２３によりそれぞれ補正された後の信号ＲＳ４のうちから、受信信号の信号強度に基づいて受信信号をＲＸ２０ごとに選択し、選択された結果を信号ＲＳ５とする。
インコヒーレント積分処理部５３（積分処理部）は、ＲＸ２０ごとに選択された結果の信号ＲＳ５に対してインコヒーレント積分処理を実施して、インコヒーレント積分処理の結果を信号ＲＳ６とする。
閾値判定処理部５４は、予め定められた閾値を判定条件にして、信号ＲＳ６の判定を実施する。当該判定により、検出対象３０として判定し得る確率が比較的高いものを検出対象３０の候補として判定する。
推定処理部５５は、インコヒーレント積分処理部５３による演算処理（積分処理）された後の信号ＲＳ６に基づいて、判定処理により選択された検出対象３０の位置を推定する。
出力処理部５６は、推定された検出対象３０の位置を出力する。
（レーダシステムにおける推定処理について）
以下、レーダシステム１の推定処理について説明する。
以下、レーダシステム１の推定処理について説明する。
「レーダシステム１における推定処理の原理」
上記のとおりレーダシステム１は、１局の送信レーダ（ＴＸ１０）と複数の受信レーダ（ＲＸ２０）を含み、ＴＸ１０とＲＸ２０とを異なる位置に設けたCentralized型のマルチスタティックレーダである。
上記のとおりレーダシステム１は、１局の送信レーダ（ＴＸ１０）と複数の受信レーダ（ＲＸ２０）を含み、ＴＸ１０とＲＸ２０とを異なる位置に設けたCentralized型のマルチスタティックレーダである。
本実施形態のレーダシステム１では、ＲＸ２０の測角性能に影響されない方法で位置の推定処理を実施する。本実施形態の位置の推定処理では、ＲＸ２０で受信される直接波の受信時刻ｔ１と、当該ＲＸ２０で反射波が受信された時刻をｔ２とし、受信時刻の時間差ΔＴｒｔを利用する。推定処理に用いるＲＸ２０の局数を複数とし、例えば３局以上にする。以下の説明では、ＲＸ２０の局数を３局とした場合の推定処理について説明する。各ＲＸ２０と検出対象３０までの距離を式（１）として定式化する。
式（１）において、ＲＸ２０−１（第１の受信局）と検出対象３０の距離をＲＲ１で示し、ＲＸ２０−２（第２の受信局）と検出対象３０の距離をＲＲ２で示し、ＲＸ２０−３（第３の受信局）と検出対象３０の距離をＲＲ３で示す。また、（ｘ１,ｙ１,ｚ１）がＲＸ２０−１の位置、（ｘ２,ｙ２,ｚ２）がＲＸ２０−２の位置、（ｘ３,ｙ３,ｚ３）がＲＸ２０−３の位置をそれぞれ示し、（ｘＴ,ｙＴ,ｚＴ）が検出対象３０の位置を示す。ｃは、伝搬速度を示す。
上記の式（１）の連立方程式を解くことで、検出対象３０の位置（ｘＴ,ｙＴ,ｚＴ）が求まる。上記の式（１）が非線形連立方程式であり、ＲＸ２０による測距誤差の影響により、同連立方程式を解けない場合がある。そのような場合には、ニュートン法などの反復による逐次計算法を適用して解を求めてもよい。
この方法では、求めるべき未知数を、初期値（０）と修正値（Δ）の和を用いて、式（２）のように仮定する。
上記の式（１）の連立方程式を解くことで、検出対象３０の位置（ｘＴ,ｙＴ,ｚＴ）が求まる。上記の式（１）が非線形連立方程式であり、ＲＸ２０による測距誤差の影響により、同連立方程式を解けない場合がある。そのような場合には、ニュートン法などの反復による逐次計算法を適用して解を求めてもよい。
この方法では、求めるべき未知数を、初期値（０）と修正値（Δ）の和を用いて、式（２）のように仮定する。
次に、上記の式（２）の修正値（Δ）の項を式（３）のように線形化して、ΔＲｉとする。
上記の式（３）において、ｉは、各ＲＸ２０を識別する識別情報であり、１以上の整数とする。ここで、ΔＲｉは、「擬似距離Ｒｉ」と「初期値から算出される推定距離」との差分として、次の式（４）のように表される。
上記の式（４）の偏微分項は、次の式（５）になる。
上記の式（５）を上から順にαｉ、βｉ、γｉとおくと、ΔＲｉは、次の式（６）のように表すことができる。
次の式（７）のようにΔＲ、Ａ、Δｘをおくことにより、上記の式（６）を式（８）と表すことができる。
上記の式（８）をΔｘについて整理して式（９）に示す。
上記の式（９）において、Δｘが十分に小さい値であれば、初期値と修正値とを加算した（ｘＴ,ｙＴ,ｚＴ）を解にする。一方、Δｘが十分に小さい値でなければ、（Δｘ,Δｙ,Δｚ）を（ｘ０,ｙ０,ｚ０）に代入して、Δｘの計算を繰り返し実施する。
次に、ＲＸ２０が４局以上にする場合について、ΔＲを次の式（１０）のように定義する。
上記の式（１０）において、εは、各ＲＸ２０における誤差を要素に含む行列（ベクトル）であり、ＲＸ２０をＮ局とする場合、式（１１）のように表すことができる。
上記の式（１１）において、「（・）Ｔ」は、転置行列を示す。ここで、最適化処理の評価手法として、この各要素のεｉの２乗和を最小にするように解を決定する最小２乗法を適用する。式（１０）を式（１２）のように変形する。
各要素のεｉの２乗和を式（１３）に示し、同式を式（１４）のように変換する。
上記のｆが変数Δｘで偏微分可能であり、ｆを最小にする変数Δｘで極値をとるとすれば、次の式（１５）と式（１６）に示す結果が導かれる。
上記の式（１６）の両辺の転置をとり、（ＡＴＢ)＝ＢＴＡＴという変換則を用いて、更に、（ＡＴＡ）が対称行列であることから、上記の式（１６）を式（１７）のように変換できる。
式（１７）をΔｘについて解くことで、式（１８）に示す修正値を導出する。
上記の式（１８）により導出したΔｘに基づいた推定処理を実施する。これにより、各ＲＸ２０による測角性能の影響を受けることなく、ＴＸ１０における送信パルスの送信時刻と検出対象３０における反射信号の受信時刻の差分ΔＴｒｔ１からΔＴｒｔＮの観測値のみを用いて検出対象３０の位置の推定が行える。上記の位置検出の方法によれば、推定誤差に影響する要因を各ＲＸ２０が有する測距誤差の影響のみに抑え込むことが可能となる。
図５と図６とを参照して、レーダシステムの検出性能について説明する。
（所要ＳＮＲに対する検出確率）
図５は、レーダシステムの検出性能として所要ＳＮＲ（Signal-to-Noise Ratio）に対する検知確率の推移を示す図である。同図には、ＲＸ２０の局数を２局から４局にした場合のレーダシステム１の検出性能を試算した結果を検知確率（Detection Probability）として示す。レーダシステム１の検出性能として、レーダシステム１における１局のＲＸ２０で得られる受信信号電力Ｓを、モノスタティックレーダの受信レーダにより得られる受信信号電力Ｓに基づいて規定する。
モノスタティックレーダの受信レーダにより得られる受信信号電力Ｓの値を、ＰＴＧＴＧＲとする。ＰＴは送信電力を示し、ＧＴは送信時のアンテナの利得を示し、ＧＲは受信時のアンテナの利得を示す。
上記のレーダシステム１における１局のＲＸ２０で得られる受信信号電力Ｓの値を、ＰＴＧＴ（ＮｘＧＲ）としたものである。なお、Ｎは、ＲＸ２０の局数を示す。なお、同図に、比較例としてモノスタティックレーダの検出特性を示す。同図に示す結果から、レーダシステム１のＲＸ２０の局数を増やすことにより、検出性能が改善できることがわかる。
図５は、レーダシステムの検出性能として所要ＳＮＲ（Signal-to-Noise Ratio）に対する検知確率の推移を示す図である。同図には、ＲＸ２０の局数を２局から４局にした場合のレーダシステム１の検出性能を試算した結果を検知確率（Detection Probability）として示す。レーダシステム１の検出性能として、レーダシステム１における１局のＲＸ２０で得られる受信信号電力Ｓを、モノスタティックレーダの受信レーダにより得られる受信信号電力Ｓに基づいて規定する。
モノスタティックレーダの受信レーダにより得られる受信信号電力Ｓの値を、ＰＴＧＴＧＲとする。ＰＴは送信電力を示し、ＧＴは送信時のアンテナの利得を示し、ＧＲは受信時のアンテナの利得を示す。
上記のレーダシステム１における１局のＲＸ２０で得られる受信信号電力Ｓの値を、ＰＴＧＴ（ＮｘＧＲ）としたものである。なお、Ｎは、ＲＸ２０の局数を示す。なお、同図に、比較例としてモノスタティックレーダの検出特性を示す。同図に示す結果から、レーダシステム１のＲＸ２０の局数を増やすことにより、検出性能が改善できることがわかる。
（受信レーダの局数に対する所要ＳＮＲの推移）
図６は、レーダシステムの検出性能として受信レーダの局数に対する所要ＳＮＲの推移を示す図である。同図には、ＲＸ２０の局数を１局から１０局にした場合のレーダシステム１の検出性能を試算した結果を所要ＳＮＲとして示す。同図に示す結果から、レーダシステム１のＲＸ２０の局数を増やすことにより、ある検出性能に必要とされる１局のＲＸ２０の所要ＳＮＲを低減することが可能となることが確認でき、検出性能が改善できることがわかる。
図６は、レーダシステムの検出性能として受信レーダの局数に対する所要ＳＮＲの推移を示す図である。同図には、ＲＸ２０の局数を１局から１０局にした場合のレーダシステム１の検出性能を試算した結果を所要ＳＮＲとして示す。同図に示す結果から、レーダシステム１のＲＸ２０の局数を増やすことにより、ある検出性能に必要とされる１局のＲＸ２０の所要ＳＮＲを低減することが可能となることが確認でき、検出性能が改善できることがわかる。
上記のとおり、マルチスタティックレーダとして構成したレーダシステム１は、ＲＸ２０の局数を増やすことで、位置の推定性能と検出性能を向上させることができる。
（分解性能）
レーダシステム１が、前述した推定性能や検出性能を得るためには、複数のＲＸ２０が検出対象３０からの反射信号を示すデータ（Ｒａｗデータ）を、当該ＲＸ２０から検出対象が存在する真の位置までの距離としてレンジング処理する必要がある。
レーダシステム１は、送信するレーダ波に対して変調処理を施し、受信した信号の復調処理を行うことでレンジング処理の分解能向上を図る。ただし、検出対象３０が移動する場合、変調したレーダ波の反射波には、移動する検出対象３０の相対速度に対応するドップラ遷移が生じる。ドップラ遷移の影響を加味せずに復調処理をした場合には、レンジング処理により得られた位置と実際の位置に誤差が生じてしまう。そこで、本実施形態のレーダシステム１は、ドップラ遷移の影響を低減するように復調処理を実施する。
レーダシステム１が、前述した推定性能や検出性能を得るためには、複数のＲＸ２０が検出対象３０からの反射信号を示すデータ（Ｒａｗデータ）を、当該ＲＸ２０から検出対象が存在する真の位置までの距離としてレンジング処理する必要がある。
レーダシステム１は、送信するレーダ波に対して変調処理を施し、受信した信号の復調処理を行うことでレンジング処理の分解能向上を図る。ただし、検出対象３０が移動する場合、変調したレーダ波の反射波には、移動する検出対象３０の相対速度に対応するドップラ遷移が生じる。ドップラ遷移の影響を加味せずに復調処理をした場合には、レンジング処理により得られた位置と実際の位置に誤差が生じてしまう。そこで、本実施形態のレーダシステム１は、ドップラ遷移の影響を低減するように復調処理を実施する。
さらに、複数の受信レーダを持つマルチスタティックレーダでは、検出対象３０に対する各受信レーダの視線角度が受信レーダごとにそれぞれ異なる。
例えば、図７を参照して各受信レーダにおけるドップラ遷移について説明する。同図は、複数の受信レーダが配置された領域を検出対象３０が移動する状況を模式化した説明図である。同図には、ＲＸ２０−１とＲＸ２０−２とＲＸ２０−３と、検出対象３０が示されている。検出対象３０は、ＲＸ２０−１とＲＸ２０−２とＲＸ２０−３のそれぞれから検出可能な範囲にあるとする。
ＲＸ２０−１と検出対象３０とを結ぶ方向に軸ＬＴ１を定め、ＲＸ２０−２と検出対象３０とを結ぶ方向に軸ＬＴ２を定め、ＲＸ２０−３と検出対象３０とを結ぶ方向に軸ＬＴ３を定める。
検出対象３０は、例えば、ＲＸ２０−１とＲＸ２０−２とＲＸ２０−３とに対して相対速度ＶＴで矢印の方向に移動する。相対速度ＶＴの軸ＬＴ１と軸ＬＴ２と軸ＬＴ３の各方向の成分は、それぞれ相対速度ＶＴ１と相対速度ＶＴ２と相対速度ＶＴ３になる。図示するように、相対速度ＶＴ１と相対速度ＶＴ２と相対速度ＶＴ３は互いに異なる大きさになる。上記のとおり各ＲＸ２０は、検出対象３０に対する相対速度が異なることにより、異なるドップラ遷移を受けたレーダ波を受信する。
このように、マルチスタティックレーダにおいて推定性能と検出性能に影響を与える「誤差」の要因となるドップラ遷移は、ＲＸ２０ごとに異なるものとなる。
例えば、図７を参照して各受信レーダにおけるドップラ遷移について説明する。同図は、複数の受信レーダが配置された領域を検出対象３０が移動する状況を模式化した説明図である。同図には、ＲＸ２０−１とＲＸ２０−２とＲＸ２０−３と、検出対象３０が示されている。検出対象３０は、ＲＸ２０−１とＲＸ２０−２とＲＸ２０−３のそれぞれから検出可能な範囲にあるとする。
ＲＸ２０−１と検出対象３０とを結ぶ方向に軸ＬＴ１を定め、ＲＸ２０−２と検出対象３０とを結ぶ方向に軸ＬＴ２を定め、ＲＸ２０−３と検出対象３０とを結ぶ方向に軸ＬＴ３を定める。
検出対象３０は、例えば、ＲＸ２０−１とＲＸ２０−２とＲＸ２０−３とに対して相対速度ＶＴで矢印の方向に移動する。相対速度ＶＴの軸ＬＴ１と軸ＬＴ２と軸ＬＴ３の各方向の成分は、それぞれ相対速度ＶＴ１と相対速度ＶＴ２と相対速度ＶＴ３になる。図示するように、相対速度ＶＴ１と相対速度ＶＴ２と相対速度ＶＴ３は互いに異なる大きさになる。上記のとおり各ＲＸ２０は、検出対象３０に対する相対速度が異なることにより、異なるドップラ遷移を受けたレーダ波を受信する。
このように、マルチスタティックレーダにおいて推定性能と検出性能に影響を与える「誤差」の要因となるドップラ遷移は、ＲＸ２０ごとに異なるものとなる。
そこで、本実施形態のレーダシステム１では、以下に示す処理を実施することにより分解性能を向上させる。
図８と図９を参照して、レーダシステム１により実行される処理について説明する。図８は、レーダシステムにより実行される処理の流れを示すシーケンス図である。図９は、ＲＸ２０における処理を模式化して示す図である。
以下に示すステップＳ１０からステップＳ２３までの処理は、所定の周期で繰り返し実施されるが、図８には１回分を明示して、繰り返し実施される処理の記載を省略している。一方、図９においては、上記の繰り返しがＺ回繰り返されるものとして示し、繰り返し実施される途中の回の記載を省略している。
［ＴＸ１０における処理］
（ステップＳ１０）：ＴＸ１０は、予め定められた変調方式により変調したレーダ波（変調波）を送出する。ＴＸ１０から送信されたレーダ波が検出対象３０において反射される。
（ステップＳ１０）：ＴＸ１０は、予め定められた変調方式により変調したレーダ波（変調波）を送出する。ＴＸ１０から送信されたレーダ波が検出対象３０において反射される。
［ＲＸ２０における処理］
ＲＸ２０のそれぞれは、以下に示すステップＳ２１からステップＳ２５の処理を実施する。代表してＲＸ２０−１の処理について説明する。
（ステップＳ２１）：ＲＸ２０−１は、アンテナ２１により反射波を受信する。
ＲＸ２０のそれぞれは、以下に示すステップＳ２１からステップＳ２５の処理を実施する。代表してＲＸ２０−１の処理について説明する。
（ステップＳ２１）：ＲＸ２０−１は、アンテナ２１により反射波を受信する。
（ステップＳ２２）：信号変換部２２は、例えば、アンテナ２１により受信した反射波に対して、ＩＦ変換、および、アナログ・ディジタル変換により、ディジタルデータへ変換する。
上記のディジタルデータは、検出対象３０の相対速度に対応するドップラ遷移ｆｄを受けた信号として次の式（１９）のように表される。
上記のディジタルデータは、検出対象３０の相対速度に対応するドップラ遷移ｆｄを受けた信号として次の式（１９）のように表される。
ｘ(t)= exp (+j (ω0+ωd)t) ・・・（１９）
上記の式（１９）において、角速度ω０は、ＴＸ１０のレーダ波の中心周波数ｆ０と次の式（２０）に示す関係にある。
ω０＝２πｆ０ ・・・（２０）
また、角速度ωｄは、ドップラ遷移周波数ｆｄと次の式（２１）に示す関係にある。
ωｄ＝２πｆｄ ・・・（２１）
（ステップＳ２３）：上記のディジタルデータに変換された信号ＲＳ２に対して、ＲＸ２０は、以下に示す復調処理をＬ個の系統についてそれぞれ実施する。
復調処理部２３は、ＴＸ１０が送信したレーダ波の周波数と、ＲＸ２０の受信信号に含まれる変調波の周波数との差として生じたドップラ遷移を打ち消す方向に、受信信号の周波数成分を補正する復調処理を実施して、補正程度を異ならせた複数の補正結果を出力する。
復調処理部２３は、ＴＸ１０が送信したレーダ波の周波数と、ＲＸ２０の受信信号に含まれる変調波の周波数との差として生じたドップラ遷移を打ち消す方向に、受信信号の周波数成分を補正する復調処理を実施して、補正程度を異ならせた複数の補正結果を出力する。
例えば、復調係数生成部２４は、上記のドップラ遷移を打ち消すように調整された復調係数として、復調処理の系統数と同数のＬ個の復調係数を生成する。Ｌ個の復調係数は、目標相対速度により生じるドップラ遷移を打ち消すよう、ドップラ補正係数を含んだ係数とする。ドップラ補正の係数を次の式（２１）に示す。
Ｗｄ１(ｉ)= exp (−ｊ(ω０+ωｄ１)(ＴＳ×(ｉ−１)))
Ｗｄ２(ｉ)= exp (−ｊ(ω０+ωｄ２)(ＴＳ×(ｉ−１)))
・・・
ＷｄＬ(ｉ)= exp (−ｊ(ω０+ωｄＬ)(ＴＳ×(ｉ−１))) ・・・（２２）
Ｗｄ２(ｉ)= exp (−ｊ(ω０+ωｄ２)(ＴＳ×(ｉ−１)))
・・・
ＷｄＬ(ｉ)= exp (−ｊ(ω０+ωｄＬ)(ＴＳ×(ｉ−１))) ・・・（２２）
上記の式（２２）において、ＴＳは受信レーダにおけるサンプリング周期であり、ｉは１からＰまで変化する変数であり、Ｐは復調係数の係数長である。また、検出対象３０の移動により生じるドップラ遷移の範囲を予測して、予測した範囲から選択したドップラ遷移ｆｄｌに対応する角周波数をωｄｌで示す。角周波数をωｄｌは、次の式（２３）に示す関係を有する。
ωｄｌ＝２πｆｄｌ ・・・（２３）
なお、上記の式（２２）において、ｌは１からＬまでの整数である。
また、ＴＸ１０によるパルス繰り返し周波数をｆＰＲＦとした場合、ｆｄｌを式（２４）に示すように定義してもよい。
また、ＴＸ１０によるパルス繰り返し周波数をｆＰＲＦとした場合、ｆｄｌを式（２４）に示すように定義してもよい。
ｆｄｌ＝（ｆＰＲＦ/Ｌ）×(ｌ−１)） ・・・（２４）
なお、上記の式（２４）からもわかるように、１／ｆＰＲＦを超えるパルスの遅延は検出することができない。
ＴＸ１０から送出されたレーダ波に対応するパルス圧縮係数をＨｉとすると、Ｌ個の復調係数は以下の式（２５）により求まる。
ＴＸ１０から送出されたレーダ波に対応するパルス圧縮係数をＨｉとすると、Ｌ個の復調係数は以下の式（２５）により求まる。
Ｈωｌ＝Ｈｉ×Ｗｄｌ ・・・（２５）
上記の式（２５）において、復調係数Ｈωｌとパルス圧縮係数Ｈｉは、ともに行列であり、Ｗｄｌは、Ｗｄ１(ｉ)からＷｄＬ(ｉ)を要素とするベクトルとする。
変換部２５は、上記の式（２５）による復調係数Ｈωｌを用いて、受信ディジタルデータを復調する。
変換部２５は、上記の式（２５）による復調係数Ｈωｌを用いて、受信ディジタルデータを復調する。
なお、適した補正量で周波数遷移の補正がなされた系統の信号レベルは、他の系統の信号レベルより大きくなる。
上記のステップＳ１０からステップＳ２３までの処理は、ＴＸ１０からパルス波が送出されるごとに一連の処理として実施される。ＴＸ１０から複数のパルス波が送出される場合には、上記の一連の処理をパルス数に応じて繰り返す。
ここで図９を参照して、ＴＸ１０からＺ個のパルス波が送出された場合について説明する。
上記のとおり、ＴＸ１０からＺ個のパルス波の送出に応じて、ディジタルデータに変換された信号ＲＳ２−１、ＲＳ２−２、・・・、ＲＳ２−Ｚが、復調処理部２３に供給される。上記の信号ＲＳ２−ｋの「ｋ」は、順に送出されたパルス波の送出回数を示す。復調処理部２３は、復調処理をＬ個の系統の復調処理をそれぞれ実施して、系統ごとに復調した結果を出力する。
例えば、第１回目のパルス波の送出に対応して、信号ＲＳ２−１が復調処理部２３に供給され、復調処理部２３は、復調処理をＬ個の系統の復調処理をそれぞれ実施して、第１の系統に信号ＲＳ３＿１−１を、第２の系統に信号ＲＳ３＿２−１を、第Ｌの系統に信号ＲＳ３＿Ｌ−１を、系統ごと出力する。信号ＲＳ３＿ｌ−ｋの「ｋ」は、順に送出されたパルス波の送出回数を示し、「ｌ」は、系統を識別する識別情報の値を示す。
復調処理部２３に供給される信号ＲＳ２−ｋは、パルス圧縮処理がなされる前の信号であり、復調処理部２３の機能により復調処理されるとともにパルス圧縮処理がなされた信号ＲＳ３＿ｌ−ｋに変換される。図中の三角形は、受信信号に基づいて生成されたパルス信号の振幅の時系列的変化を示す。同図の場合、信号ＲＳ３＿１−１と信号ＲＳ３＿２−１と信号ＲＳ３＿Ｌ−１のうちで、信号ＲＳ３＿２−１の振幅が他の系統の信号より大きい。第１回目から第Ｚ回目まで通してみても、第２の系統の信号ＲＳ３＿２−ｋの振幅が他の系統の信号より大きい。上記のような傾向がみられる場合には、第２の系統の補正量が、検出対象３０の移動により生じたドップラ遷移量に適した補正がなされたと判断できる。
ここで図９を参照して、ＴＸ１０からＺ個のパルス波が送出された場合について説明する。
上記のとおり、ＴＸ１０からＺ個のパルス波の送出に応じて、ディジタルデータに変換された信号ＲＳ２−１、ＲＳ２−２、・・・、ＲＳ２−Ｚが、復調処理部２３に供給される。上記の信号ＲＳ２−ｋの「ｋ」は、順に送出されたパルス波の送出回数を示す。復調処理部２３は、復調処理をＬ個の系統の復調処理をそれぞれ実施して、系統ごとに復調した結果を出力する。
例えば、第１回目のパルス波の送出に対応して、信号ＲＳ２−１が復調処理部２３に供給され、復調処理部２３は、復調処理をＬ個の系統の復調処理をそれぞれ実施して、第１の系統に信号ＲＳ３＿１−１を、第２の系統に信号ＲＳ３＿２−１を、第Ｌの系統に信号ＲＳ３＿Ｌ−１を、系統ごと出力する。信号ＲＳ３＿ｌ−ｋの「ｋ」は、順に送出されたパルス波の送出回数を示し、「ｌ」は、系統を識別する識別情報の値を示す。
復調処理部２３に供給される信号ＲＳ２−ｋは、パルス圧縮処理がなされる前の信号であり、復調処理部２３の機能により復調処理されるとともにパルス圧縮処理がなされた信号ＲＳ３＿ｌ−ｋに変換される。図中の三角形は、受信信号に基づいて生成されたパルス信号の振幅の時系列的変化を示す。同図の場合、信号ＲＳ３＿１−１と信号ＲＳ３＿２−１と信号ＲＳ３＿Ｌ−１のうちで、信号ＲＳ３＿２−１の振幅が他の系統の信号より大きい。第１回目から第Ｚ回目まで通してみても、第２の系統の信号ＲＳ３＿２−ｋの振幅が他の系統の信号より大きい。上記のような傾向がみられる場合には、第２の系統の補正量が、検出対象３０の移動により生じたドップラ遷移量に適した補正がなされたと判断できる。
（ステップＳ２４）：ステップＳ１０からステップＳ２３までの処理を、ＴＸ１０が送出したパルス数に応じた回数分繰り返した後、パルス積分処理部２６は、復調処理の系統ごとに、信号ＲＳ３における所定数の受信パルスをパルス積分する。所定数の受信パルスとは、ＴＸ１０が周期的に繰り返し送出したパルスの数に応じた数の受信パルスのことである。例えば、ＴＸ１０が周期的に繰り返し送出したパルスの数に、パルス積分の対象にするパルスの数を一致させる。なお、パルス積分を実施することにより、単一パルスの信号を処理する場合に比べてＳＮＲを改善することができる。
（ステップＳ２５）：通信部２７は、パルス積分処理部２６によるパルス積分の結果(信号ＲＳ４)を演算処理部５０宛に通知する。
ＲＸ２０−２からＲＸ２０−Ｎは、ＲＸ２０−１と同様に、ステップＳ２１からステップＳ２５までの処理をそれぞれ実施する。
上記のステップＳ２１からステップＳ２５までの処理により、ＲＸ２０のそれぞれから通知されたパルス積分の結果（信号ＲＳ４）には、何れかの系統のパルス積分の結果にドップラ遷移の影響を最も低減できた系統が含まれる。
上記のステップＳ２１からステップＳ２５までの処理により、ＲＸ２０のそれぞれから通知されたパルス積分の結果（信号ＲＳ４）には、何れかの系統のパルス積分の結果にドップラ遷移の影響を最も低減できた系統が含まれる。
［演算処理部５０における処理］
（ステップＳ５１）：演算処理部５０において、通信部５１は、ＲＸ２０からそれぞれ通知されたデータをそれぞれ受信する。
（ステップＳ５１）：演算処理部５０において、通信部５１は、ＲＸ２０からそれぞれ通知されたデータをそれぞれ受信する。
（ステップＳ５２）：選択部５２は、ＲＸ２０の復調処理部２３によりそれぞれ補正された後の信号ＲＳ４のうちの何れかの受信信号を、受信信号の信号強度に基づいてＲＸ２０ごとに選択し、選択された結果を信号ＲＳ５とする。上記の選択処理では、信号レベルが最も大きくなった系統の補正量が適していると評価する評価基準に基づいて、適切な補正がなされた系統を信号レベルの大きさに基づいて選択し、その系統の信号を選択する。例えば、図９の場合では、第２の系統が適した補正処理を実施した例を示した。そのような場合には、選択部５２は、第２の系統の信号を選択することになる。
（ステップＳ５３）：インコヒーレント積分処理部５３（積分処理部）は、ステップＳ５２においてＲＸ２０ごとに選択された系統の信号ＲＳ５に対してのインコヒーレント積分処理を実施して、その結果を信号ＲＳ６とする。インコヒーレント積分では、位相成分を用いることなく信号強度に基づいた積分を実施する。例えば、インコヒーレント積分の対象単位は、レンジ方向に定めた所定の距離と、ＲＸ２０の測角方向に定めた所定の角度とにより囲まれる領域を単位とする。複数のＲＸ２０の信号をインコヒーレント積分するに当たり、それぞれのＲＸ２０のインコヒーレント積分の結果を共通の座標系のセルにマッピングする。上記のインコヒーレント積分処理では、ＲＸ２０ごとに選択された系統の信号ＲＳ５の信号強度を、全てのＲＸ２０についてインコヒーレント積分して、その結果を信号ＲＳ６にする。
（ステップＳ５４）：閾値判定処理部５４は、予め定められた閾値を判定条件にして、信号ＲＳ６を閾値判定する。当該判定により、当該閾値を超える信号強度の地点を検出対象３０として判定し得る確率が比較的高いセルを検出対象３０の候補として判定する。
（ステップＳ５５）：推定処理部５５は、インコヒーレント積分処理（積分処理）された後の信号ＲＳ６から選択されたセルの位置情報から、検出対象３０の位置を推定する。
（ステップＳ５６）：出力処理部５６は、推定された検出対象３０の位置を出力する。
レーダシステム１は、上記の手順に従って、検出対象３０の相対速度に対する適正な補正量により補正された後の信号を、補正後の信号強度から特定する。レーダシステム１は、補正後の信号強度が最大値を示した場合は、その補正量が適していると判定し、適した補正量で補正した信号に基づいて、検出対象３０の相対速度によるドップラ遷移の影響を低減させたレンジング処理を実施する。
（レーダシステムの処理による効果について）
レーダシステム１は、１局の送信レーダ（ＴＸ１０）と複数の受信レーダ（ＲＸ２０）を含み、ＴＸ１０とＲＸ２０とを異なる位置に設けたマルチスタティックレーダである。
レーダシステム１は、１局の送信レーダ（ＴＸ１０）と複数の受信レーダ（ＲＸ２０）を含み、ＴＸ１０とＲＸ２０とを異なる位置に設けたマルチスタティックレーダである。
比較例として、送信レーダと受信レーダとを異なる位置に設ける方式のレーダシステムの一種であるバイスタティックレーダを挙げて、レーダシステム１による効果を説明する。バイスタティックレーダは、送信レーダ１局と受信レーダ１局とを組にして構成され、マルチスタティックレーダに比べると受信レーダの局数が異なる。バイスタティックレーダに関する詳細は、Nicholas, J Willis, Bistatic Radar, Artech House 1991.を参照する。
マルチスタティックレーダは、バイスタティックレーダと対比すると受信レーダの局数が異なるものであるが、バイスタティックレーダと同様の推定処理を適用しても十分な性能を得ることができない。
ここで、バイスタティックレーダによる推定原理を、マルチスタティックレーダに適用した場合の問題について説明する。
前述の図１を参照して、受信レーダの局数をＮとした場合、マルチスタティックレーダによる推定処理の座標系について説明する。
ＴＸ１０とＲＸ２０のそれぞれの間の距離をベースラインと呼び、Ｌｋで示す。但し、ｋは、ＲＸ２０の識別情報を示し、例えば１以上Ｎ以下の整数である。
ＲＸ２０−ｋと検出対象３０との距離を、ＲＲｋで示す。ＲＲｋとして、例えば、ＲＲ１とＲＲ２とＲＲＮとが図示されている。
ＴＸ１０と検出対象３０までの距離をＲＴで示す。
ＲＸ２０−ｋにおける検出対象３０からの反射波が到来する方向を到来方向θＲｋで示す。到来方向θｋとして、例えば、θＲ１とθＲ２とθＲＮとが図示されている。
ＲＸ２０−ｋにおいて、ベースラインＬｋを基軸とし、検出対象３０からの反射波の到来方向θＲｋまでの角度をαｋで示す。例えば、ＲＸ２０−１において、ベースラインＬ１を基軸とし、検出対象３０からの反射波の到来方向θＲ１までの角度をα１で示す。
ＲＸ２０で受信される直接波の時刻をｔ１、当該ＲＸ２０で反射波が受信された時刻をｔ２とし、受信時刻の時間差ΔＴｒｔを式（２６）とおくことにより、当該受信時刻に対する推定距離ＲＲｋは、式（２７）で計算される。
前述の図１を参照して、受信レーダの局数をＮとした場合、マルチスタティックレーダによる推定処理の座標系について説明する。
ＴＸ１０とＲＸ２０のそれぞれの間の距離をベースラインと呼び、Ｌｋで示す。但し、ｋは、ＲＸ２０の識別情報を示し、例えば１以上Ｎ以下の整数である。
ＲＸ２０−ｋと検出対象３０との距離を、ＲＲｋで示す。ＲＲｋとして、例えば、ＲＲ１とＲＲ２とＲＲＮとが図示されている。
ＴＸ１０と検出対象３０までの距離をＲＴで示す。
ＲＸ２０−ｋにおける検出対象３０からの反射波が到来する方向を到来方向θＲｋで示す。到来方向θｋとして、例えば、θＲ１とθＲ２とθＲＮとが図示されている。
ＲＸ２０−ｋにおいて、ベースラインＬｋを基軸とし、検出対象３０からの反射波の到来方向θＲｋまでの角度をαｋで示す。例えば、ＲＸ２０−１において、ベースラインＬ１を基軸とし、検出対象３０からの反射波の到来方向θＲ１までの角度をα１で示す。
ＲＸ２０で受信される直接波の時刻をｔ１、当該ＲＸ２０で反射波が受信された時刻をｔ２とし、受信時刻の時間差ΔＴｒｔを式（２６）とおくことにより、当該受信時刻に対する推定距離ＲＲｋは、式（２７）で計算される。
ΔＴｒｔ＝ｔ２−ｔ１ ・・・（２６）
なお、ＲＸ２０でそれぞれ推定された複数の目標位置に対して複数点の重心を算出し、当該重心位置を統合推定位置として検出する方式を適用する。統合推定位置と、実際の目標位置との差分を推定誤差として評価する。
図１０を参照して、比較例のレーダシステムにおける誤差要因について説明する。同図は、誤差要因をモデル化して示す説明図である。
上記の式（２７）からもわかるように、各ＲＸ２０の受信ビームには、観測時間から導かれる測距成分の誤差（±√（σＲ））と、角度成分の誤差（±√（σａｚ））と、高度成分の誤差（±√（σｅｌ））とが含まれる。上記の誤差が、検出対象３０の位置の検出性能に影響する。
図１０を参照して、比較例のレーダシステムにおける誤差要因について説明する。同図は、誤差要因をモデル化して示す説明図である。
上記の式（２７）からもわかるように、各ＲＸ２０の受信ビームには、観測時間から導かれる測距成分の誤差（±√（σＲ））と、角度成分の誤差（±√（σａｚ））と、高度成分の誤差（±√（σｅｌ））とが含まれる。上記の誤差が、検出対象３０の位置の検出性能に影響する。
一方、本実施形態のレーダシステム１の場合について、図１１を参照して説明する。同図は、本実施形態のレーダシステムにおける測距成分の誤差について示す説明図である。同図に示すように、各ＲＸ２０の距離成分の誤差（±√（σＲ））に制限される。
マルチスタティックレーダに、バイスタティックレーダの推定原理に基づいた定式と同様の定式を当てはめた場合の推定精度は、送信レーダと受信レーダの配置によって異なるものとなる。そこで、送信レーダと受信レーダと検出対象との位置をモデル化し、当該モデルにおけるマルチスタティックレーダによる推定精度を定量的に評価する。上記のモデルでは、ＴＸ１０から検出対象３０までの距離（ＲＴ）と、検出対象３０からＲＸ２０までの距離（ＲＲｋ）とを一律に３０ｋｍとする。上記のモデルを、半径が３０ｋｍの円を用いて説明する。円の中心に検出対象３０が配置され、円周上にＴＸ１０とＲＸ２０のそれぞれが配置される。検出対象３０の位置からＴＸ１０の方向とＲＸ２０の方向とがなす角をバイスタティックアングルという。バイスタティックレーダにおけるバイスタティックアングルを、マルチスタティックレーダのマルチスタティックアングルとみなす。このモデルを用いて、マルチスタティックアングルに対する推定誤差推移を検討する。
検討の条件として定めたパラメータ諸元を図１２に示す。同図は、検討の条件として定めたパラメータ諸元について示す説明図である。同図に示す通り、パラメータ諸元として、例えば、送信電力「ＰＴ（ｄＢ）」と、送信アンテナ利得「ＧＴ（ｄＢ）」と、受信アンテナ利得「ＧＴ−１０ｌｏｇＮ（ｄＢ）」と、ＳＮＲ「（１２．８−１０ｌｏｇＮ）（ｄＢ）」と、３ｄＢビーム幅「（１．５ｘＮ）（ｄｅｇ）」と、レンジ分解能「３７．５（ｍ）」とを含む。パラメータ諸元として、さらに、モノパルス係数について、「方位」、「迎角」、「距離」の各項目を含む。図示する数値は、一例を示すものであり、適宜変更してもよい。
図１３を参照して、レーダシステムのマルチスタティックアングルによる推定誤差の推移について説明する。同図は、レーダシステムのマルチスタティックアングルによる推定誤差の推移を示す説明図である。なお、検討に当たり、ＴＸ１０の送信電力やＲＸ２０の受信アンテナ利得等のレーダリソースは各種形態で同一とした。
同図に本実施形態のレーダシステム１の場合の推定誤差をＭＳＥで示す。また、同図に比較例として、バイスタティックレーダの場合の推定誤差をＢＥで示し、マルチスタティックレーダの場合の推定誤差をＭＢＥ（Ｍ）で示す。比較例として示すマルチスタティックレーダは、バイスタティックレーダと同様の推定原理を用いたものである。なお、「（Ｍ）」は、受信レーダの局数を示す。
同図に本実施形態のレーダシステム１の場合の推定誤差をＭＳＥで示す。また、同図に比較例として、バイスタティックレーダの場合の推定誤差をＢＥで示し、マルチスタティックレーダの場合の推定誤差をＭＢＥ（Ｍ）で示す。比較例として示すマルチスタティックレーダは、バイスタティックレーダと同様の推定原理を用いたものである。なお、「（Ｍ）」は、受信レーダの局数を示す。
同図に示すように、バイスタティックレーダの推定誤差ＢＥと、バイスタティックレーダの推定原理に基づくマルチスタティックレーダの推定誤差ＭＢＥ２から推定誤差ＭＢＥ４は、マルチスタティックアングル（バイスタティックアングルアングル）が大きくなるほど、また、マルチスタティックレーダの受信レータの局数が多くなるほど推定誤差が大きくなる。マルチスタティックレーダは、マルチスタティックアングルを有することが特徴である。ただし、上記の式（２７）に示した定式を用いた場合に生じる上記の傾向は、マルチスタティックアングル（バイスタティックアングルアングル）を大きくするほど推定誤差が大きくなるということであり、望ましいものではない。このように、比較例として示したバイスタティックレーダと、バイスタティックレーダの推定原理に基づくマルチスタティックレーダの双方に満足な結果が得られない。
一方、本実施形態のレーダシステム１の指定誤差ＭＳＥは、マルチスタティックアングルが大きくなるにしたがって減少する。このような傾向は、マルチスタティックレーダにとって適したものであり、指定誤差ＭＳＥは、マルチスタティックアングルが５度を超える領域で、比較例の推定誤差に対して充分に小さな値とる。ＴＸ１０とＲＸ２０を分離して配置するマルチスタティックレーダでは、マルチスタティックアングルが５度以下になるようにＴＸ１０とＲＸ２０を配置することが必要とされる場合は極めて少ない。マルチスタティックアングルが５度以下の領域で、指定誤差ＭＳＥが比較例より大きな値を示すものの、上記のＴＸ１０とＲＸ２０の配置の条件をふまえれば、実際の構成で問題になることはない。
以上に説明した、第１の実施形態によれば、レーダシステム１は、ＲＸ２０と、復調処理部２３（補正部）と、演算処理部５０（推定部）とを持つ。ＲＸ２０は、ＴＸ１０から送信された電波が検出対象３０において反射されて、前記反射された電波の周波数成分を含む受信信号を受信する。復調処理部２３は、ＴＸ１０が送信した電波の周波数と、前記受信信号に含まれる前記反射された電波の周波数との差として生じたドップラ遷移を打ち消す方向に、前記受信信号の周波数成分を補正する。演算処理部５０は、移動する検出対象３０から前記ＲＸ２０までの距離を、復調処理部２３により補正された前記受信信号の周波数成分に基づいて推定することによって、検出対象３０までの距離の検出精度を高めることができる。
（第１の実施形態の変形例）
以下、第１の実施形態の変形例について説明する。第１の実施形態の変形例において、パルスの遅延量が大きくなる比較的速い速度で移動する検出対象３０を検出する構成について説明する。より具体的には、第１の実施形態の場合にドップラ遷移折り返しにより制限されていた範囲を、本変形例では検出対象範囲にする。以下、この点を中心に説明する。
以下、第１の実施形態の変形例について説明する。第１の実施形態の変形例において、パルスの遅延量が大きくなる比較的速い速度で移動する検出対象３０を検出する構成について説明する。より具体的には、第１の実施形態の場合にドップラ遷移折り返しにより制限されていた範囲を、本変形例では検出対象範囲にする。以下、この点を中心に説明する。
ＴＸ１０が送出するパルスの送出周期より、ドップラ遷移による伝搬遅延が大きくなると、ドップラ遷移折り返しが発生する。そこで、復調処理部２３における処理の系統数、すなわち、変調係数の個数を、ドップラ遷移折り返しを考慮した個数にすることで、移動速度の比較的速い検出対象３０の検出が可能になる。
レーダシステム１では、復調係数の個数を、ドップラ遷移折り返しの折り返し数を想定し、想定する折り返し数をＤとする拡張により、（ＬｘＤ）個にする。具体的には、折り返し数をＤにした場合、ドップラ補正係数は、式（２９）によって示される条件のもとで式（２８）として表すことができる。
レーダシステム１では、復調係数の個数を、ドップラ遷移折り返しの折り返し数を想定し、想定する折り返し数をＤとする拡張により、（ＬｘＤ）個にする。具体的には、折り返し数をＤにした場合、ドップラ補正係数は、式（２９）によって示される条件のもとで式（２８）として表すことができる。
Ｗｄｌ(ｉ)= exp (−ｊ(ω０+ωｄｌ)(ＴＳ×(ｉ−１)))・・・（２８）
ｆｄｌ＝（ｆＰＲＦ/Ｌ）×ｄ(ｌ−１)） ・・・（２９）
なお、上記の式（２８）と式（２９）において、ＴＸ１０によるパルス繰り返し周波数をｆＰＲＦとした場合、ｄは１からＤまでの整数であり、ｌは１から（Ｌ×Ｄ）までの整数である。
なお、レーダシステム１Ａは、上記のように拡張したことにより、１／ｆＰＲＦを超えるパルスの遅延も検出することができる。
上記の拡張により、復調処理部２３から選択部５２までの間で扱う系統数が増加するが、上記に示した以外の基本的な処理は、第１の実施形態と同様に実施することができる。
なお、レーダシステム１Ａは、上記のように拡張したことにより、１／ｆＰＲＦを超えるパルスの遅延も検出することができる。
上記の拡張により、復調処理部２３から選択部５２までの間で扱う系統数が増加するが、上記に示した以外の基本的な処理は、第１の実施形態と同様に実施することができる。
以上に説明した、第１の実施形態の変形例によれば、第１の実施形態と同様の効果を奏するものに加え、相対移動速度が比較的速く、第１の実施形態で検出可能な相対速度範囲を超える検出対象３０についても検出対象にすることができる。この結果、第１の実施形態の変形例によれば、検出対象３０の移動状況に応じて、より柔軟な処理を実現することができる。
（第２の実施形態）
以下、第２の実施形態について説明する。第２の実施形態において、Decentralized型のマルチスタティックレーダへの適用について説明する。以下、この点を中心に説明する。
本実施形態のレーダシステム１Ａは、上記の図３に示した第２の方式のDecentralized型に分類される。
ＲＸ２０Ａは、それぞれ受信した反射波に基づいた信号処理を行った後に、検出対象３０の検出を実施して、検出対象３０の検出結果を示すデータをＰｌｏｔデータとして演算処理部５０に提供する。演算処理部５０は、ＲＸ２０Ａのそれぞれから提供されたＰｌｏｔデータを解析して、検出対象３０の位置を推定する。
以下、第２の実施形態について説明する。第２の実施形態において、Decentralized型のマルチスタティックレーダへの適用について説明する。以下、この点を中心に説明する。
本実施形態のレーダシステム１Ａは、上記の図３に示した第２の方式のDecentralized型に分類される。
ＲＸ２０Ａは、それぞれ受信した反射波に基づいた信号処理を行った後に、検出対象３０の検出を実施して、検出対象３０の検出結果を示すデータをＰｌｏｔデータとして演算処理部５０に提供する。演算処理部５０は、ＲＸ２０Ａのそれぞれから提供されたＰｌｏｔデータを解析して、検出対象３０の位置を推定する。
（レーダシステムの構成）
図１４は、第２の実施形態に係るレーダシステム１Ａの構成図である。レーダシステム１Ａは、例えば、複数のＲＸ２０Ａと、演算処理部５０Ａとを備える。レーダシステム１Ａは、さらに、ＴＸ１０を備えていてもよい。
図１４は、第２の実施形態に係るレーダシステム１Ａの構成図である。レーダシステム１Ａは、例えば、複数のＲＸ２０Ａと、演算処理部５０Ａとを備える。レーダシステム１Ａは、さらに、ＴＸ１０を備えていてもよい。
[ＲＸ]
ＲＸ２０Ａは、ＴＸ１０から送信された変調波（電波）が検出対象３０において反射されて、反射された電波の周波数成分を含む受信信号をアンテナ２１によって受信する。ＲＸ２０Ａは、アンテナ２１によって受信した受信信号（信号ＲＳ１）に基づいて、ＴＸ１０が送信した変調波の周波数と、信号ＲＳ１に含まれる変調波の周波数との差として生じたドップラ遷移を打ち消す方向に、信号ＲＳ１の周波数成分を補正する。ＲＸ２０Ａは、例えば、信号変換部２２と、復調処理部２３と、パルス積分処理部２６と、閾値判定処理部２８と、通信部２７Ａとを備える。
ＲＸ２０Ａは、ＴＸ１０から送信された変調波（電波）が検出対象３０において反射されて、反射された電波の周波数成分を含む受信信号をアンテナ２１によって受信する。ＲＸ２０Ａは、アンテナ２１によって受信した受信信号（信号ＲＳ１）に基づいて、ＴＸ１０が送信した変調波の周波数と、信号ＲＳ１に含まれる変調波の周波数との差として生じたドップラ遷移を打ち消す方向に、信号ＲＳ１の周波数成分を補正する。ＲＸ２０Ａは、例えば、信号変換部２２と、復調処理部２３と、パルス積分処理部２６と、閾値判定処理部２８と、通信部２７Ａとを備える。
閾値判定処理部２８は、予め定められた閾値を判定条件にして、パルス積分処理部２６の演算処理の結果について判定する。当該判定により、検出対象３０として判定し得る確率が比較的高いものが検出対象３０の候補として判定される。閾値判定処理部２８は、判定の結果を、Ｐｌｏｔデータとして出力する。
通信部２７Ａは、閾値判定処理部２８による判定の結果（Ｐｌｏｔデータ）を演算処理部５０Ａ宛に通知する。
演算処理部５０Ａは、ＲＸ２０Ａのそれぞれから通知されたＰｌｏｔデータに基づいて、検出対象３０の位置を推定する。演算処理部５０は、例えば、通信部５１Ａと、推定処理部５５Ａと、出力処理部５６とを備える。
通信部５１Ａは、ＲＸ２０Ａからのデータを受信する。
推定処理部５５Ａは、ＲＸ２０Ａのそれぞれから通知されたＰｌｏｔデータに基づいて、検出対象３０の位置を推定する。
出力処理部５６は、推定された検出対象３０の位置を出力する。
通信部５１Ａは、ＲＸ２０Ａからのデータを受信する。
推定処理部５５Ａは、ＲＸ２０Ａのそれぞれから通知されたＰｌｏｔデータに基づいて、検出対象３０の位置を推定する。
出力処理部５６は、推定された検出対象３０の位置を出力する。
図１５を参照して、レーダシステム１Ａの処理について説明する。同図は、レーダシステム１Ａの処理を示すシーケンス図である。
同図に示すステップＳ１０からステップＳ２４までの処理は、前述の図８を参照し、レーダシステム１をレーダシステム１Ａに、ＲＸ２０をＲＸ２０Ａに、演算処理部５０を演算処理部５０Ａに読み替える。
［ＲＸ２０における処理］
（ステップＳ２６）：ステップＳ２４の処理を終えた後、閾値判定処理部２８は、予め定められた閾値を判定条件にして、パルス積分処理部２６の演算処理の結果について判定する。閾値判定処理部２８は、上記の判定において、検出対象３０として判定し得る確率が比較的高いものを検出対象３０の候補として判定し、判定の結果を、Ｐｌｏｔデータとして出力する。
例えば、閾値判定処理部２８は、パルス積分処理部２６の演算処理の結果が示す信号の信号成分が最も大きな値を示す系統を選択する。閾値判定処理部２８は、選択した系統の信号から、検出対象３０として判定し得る確率が比較的高いものを検出対象３０の候補として判定する。
（ステップＳ２６）：ステップＳ２４の処理を終えた後、閾値判定処理部２８は、予め定められた閾値を判定条件にして、パルス積分処理部２６の演算処理の結果について判定する。閾値判定処理部２８は、上記の判定において、検出対象３０として判定し得る確率が比較的高いものを検出対象３０の候補として判定し、判定の結果を、Ｐｌｏｔデータとして出力する。
例えば、閾値判定処理部２８は、パルス積分処理部２６の演算処理の結果が示す信号の信号成分が最も大きな値を示す系統を選択する。閾値判定処理部２８は、選択した系統の信号から、検出対象３０として判定し得る確率が比較的高いものを検出対象３０の候補として判定する。
（ステップＳ２７）：通信部２７Ａは、閾値判定処理部２８による判定の結果（Ｐｌｏｔデータ）を演算処理部５０Ａ宛に通知する。
ＲＸ２０Ａ−２からＲＸ２０Ａ−Ｎは、ＲＸ２０Ａ−１と同様に、ステップＳ２１からステップＳ２７までの処理をそれぞれ実施する。
上記のステップＳ２１からステップＳ２７までの処理により、ＲＸ２０Ａのそれぞれから通知されたＰｌｏｔデータには、ドップラ遷移の影響を最も低減できた系統のＰｌｏｔデータ含まれる。
上記のステップＳ２１からステップＳ２７までの処理により、ＲＸ２０Ａのそれぞれから通知されたＰｌｏｔデータには、ドップラ遷移の影響を最も低減できた系統のＰｌｏｔデータ含まれる。
［演算処理部５０Ａにおける処理］
（ステップＳ５１Ａ）：演算処理部５０Ａにおいて、通信部５１Ａは、ＲＸ２０Ａからそれぞれ通知されたデータをそれぞれ受信する。
（ステップＳ５１Ａ）：演算処理部５０Ａにおいて、通信部５１Ａは、ＲＸ２０Ａからそれぞれ通知されたデータをそれぞれ受信する。
（ステップＳ５５Ａ）：推定処理部５５Ａは、ＲＸ２０Ａのそれぞれから通知されたＰｌｏｔデータに基づいて、検出対象３０の位置を推定する。
（ステップＳ５６）：出力処理部５６は、推定された検出対象３０の位置を出力する。
レーダシステム１Ａは、上記の手順に従って、検出対象３０の相対速度に対する適正な補正量により補正された後の信号を、補正後の信号強度から特定する。レーダシステム１Ａは、補正後の信号強度が最大値を示した場合は、その補正量が適していると判定し、適した補正量で補正した信号に基づいて、検出対象３０の相対速度によるドップラ遷移の影響を低減させたレンジング処理を実施する。
以上に説明した、第２の実施形態によれば、レーダシステム１Ａは、ＲＸ２０Ａと、復調処理部２３（補正部）と、演算処理部５０Ａ（推定部）とを持つ。ＲＸ２０Ａは、ＴＸ１０から送信された電波が検出対象３０において反射されて、前記反射された電波の周波数成分を含む受信信号を受信する。復調処理部２３は、ＴＸ１０が送信した電波の周波数と、前記受信信号に含まれる前記反射された電波の周波数との差として生じたドップラ遷移を打ち消す方向に、前記受信信号の周波数成分を補正する。演算処理部５０Ａは、移動する検出対象３０から前記ＲＸ２０Ａまでの距離を、復調処理部２３により補正された前記受信信号の周波数成分に基づいて推定することによって、検出対象３０までの距離の検出精度を高めることができる。
第２の実施形態によれば、上記の第１の実施形態と同様の効果を奏するものに加え、Decentralized型のマルチスタティックレーダにおいて、検出対象３０までの距離の検出精度を高めることができることにより、検出対象３０の移動状況に応じて、より柔軟な処理を実現することができる。
第２の実施形態によれば、上記の第１の実施形態と同様の効果を奏するものに加え、Decentralized型のマルチスタティックレーダにおいて、検出対象３０までの距離の検出精度を高めることができることにより、検出対象３０の移動状況に応じて、より柔軟な処理を実現することができる。
（第３の実施形態）
以下、第３の実施形態について説明する。第３の実施形態において、レーダシステム１Ｂへの適用を例示して、本実施形態の信号処理システムについて説明する。以下、この点を中心に説明する。
以下、第３の実施形態について説明する。第３の実施形態において、レーダシステム１Ｂへの適用を例示して、本実施形態の信号処理システムについて説明する。以下、この点を中心に説明する。
（信号処理システムの構成）
図１６は、第３の実施形態に係る信号処理システム６０の構成図である。
同図に示された信号処理システム６０は、例えばレーダシステム１Ｂの一部として構成される。このようなレーダシステム１Ｂは、上記の図１に示した第１の方式のCentralized型に分類される。
レーダシステム１Ｂは、例えば、複数のＲＸ２０Ｂと、信号処理システム６０とを含む。レーダシステム１Ｂは、さらに、ＴＸ１０を備えていてもよい。
図１６は、第３の実施形態に係る信号処理システム６０の構成図である。
同図に示された信号処理システム６０は、例えばレーダシステム１Ｂの一部として構成される。このようなレーダシステム１Ｂは、上記の図１に示した第１の方式のCentralized型に分類される。
レーダシステム１Ｂは、例えば、複数のＲＸ２０Ｂと、信号処理システム６０とを含む。レーダシステム１Ｂは、さらに、ＴＸ１０を備えていてもよい。
［ＲＸ］
ＲＸ２０Ｂは、ＴＸ１０から送信された変調波（電波）が検出対象３０において反射されて、反射された電波を含む受信信号をアンテナ２１によって受信する。ＲＸ２０Ｂは、それぞれ受信した反射波に基づいた信号変換処理を行い、ディジタルデータを信号処理システム６０に提供する。
ＲＸ２０Ｂは、ＴＸ１０から送信された変調波（電波）が検出対象３０において反射されて、反射された電波を含む受信信号をアンテナ２１によって受信する。ＲＸ２０Ｂは、それぞれ受信した反射波に基づいた信号変換処理を行い、ディジタルデータを信号処理システム６０に提供する。
ＲＸ２０Ｂ−１（ＲＸ２０Ｂ）は、例えば、信号変換部２２Ｂを備える。
信号変換部２２Ｂは、アンテナ２１によって受信した受信信号（信号ＲＳ１）に基づいて、信号ＲＳ１の周波数変換とアナログ・ディジタル変換とを実施する。信号変換部２２は、信号ＲＳ１に含まれた所望の信号成分を含む信号ＲＳ２をディジタル信号として出力する。
信号変換部２２Ｂは、アンテナ２１によって受信した受信信号（信号ＲＳ１）に基づいて、信号ＲＳ１の周波数変換とアナログ・ディジタル変換とを実施する。信号変換部２２は、信号ＲＳ１に含まれた所望の信号成分を含む信号ＲＳ２をディジタル信号として出力する。
ＲＸ２０Ｂ−２からＲＸ２０Ｂ−Ｎは、例えば、ＲＸ２０Ｂ−１と同様の構成を備える。それらの説明はＲＸ２０Ｂ−１の説明を参照する。
［信号処理システム］
信号処理システム６０は、例えば、複数の補正部２９と、演算処理部５０とを備える。
補正部２９は、ＲＸ２０Ｂに対応させて設けられ、ＲＸ２０Ｂにより生成されたディジタル信号を補正する。
補正部２９−１（補正部２９）は、例えば、復調処理部２３と、パルス積分処理部２６と、通信部２７とを備える。
補正部２９−１（補正部２９）は、ＲＸ２０Ｂ−１の信号変換部２２Ｂにより生成されたディジタル信号に基づいて、ＴＸ１０が送信した変調波（電波）の周波数と、信号ＲＳ１に含まれる変調波の周波数との差として生じたドップラ遷移を打ち消す方向に、信号ＲＳ１の周波数成分を、復調処理部２３により補正する。
補正部２９−１は、パルス積分処理部２６によるパルス積分の結果(信号ＲＳ４)を、通信部２７を介して演算処理部５０宛に通知する。
信号処理システム６０は、例えば、複数の補正部２９と、演算処理部５０とを備える。
補正部２９は、ＲＸ２０Ｂに対応させて設けられ、ＲＸ２０Ｂにより生成されたディジタル信号を補正する。
補正部２９−１（補正部２９）は、例えば、復調処理部２３と、パルス積分処理部２６と、通信部２７とを備える。
補正部２９−１（補正部２９）は、ＲＸ２０Ｂ−１の信号変換部２２Ｂにより生成されたディジタル信号に基づいて、ＴＸ１０が送信した変調波（電波）の周波数と、信号ＲＳ１に含まれる変調波の周波数との差として生じたドップラ遷移を打ち消す方向に、信号ＲＳ１の周波数成分を、復調処理部２３により補正する。
補正部２９−１は、パルス積分処理部２６によるパルス積分の結果(信号ＲＳ４)を、通信部２７を介して演算処理部５０宛に通知する。
補正部２９−２から補正部２９−Ｎは、例えば、補正部２９−１と同様の構成を備える。それらの説明は補正部２９−１の説明を参照する。
演算処理部５０は、ＲＸ２０のそれぞれから通知されたパルス積分の結果（信号ＲＳ４）に基づいて、検出対象３０の位置を推定する。
以上のように構成された信号処理システム６０は、第１の実施形態のレーダシステム１と同様の処理を実施する。
以上に説明した、第３の実施形態によれば、信号処理システム６０は、補正部２９と、演算処理部５０（推定部）とを持つ。補正部２９は、ＴＸ１０から送信された電波の周波数と、検出対象において反射されて受信局により受信された受信信号に含まれる前記反射された電波の周波数との差として生じたドップラ遷移を打ち消す方向に、前記受信信号の周波数成分を補正する。演算処理部５０は、移動する検出対象３０から前記ＲＸ２０Ｂまでの距離を、補正部２９により補正された前記受信信号の周波数成分に基づいて推定することによって、検出対象３０までの距離の検出精度を高めることができる。
第３の実施形態によれば、上記の第１の実施形態と同様の効果を奏するものに加え、新御具処理システム６０は、ＲＸ２０Ｂの信号変換部２２により生成されたディジタル信号に対する信号処理を実施することにより、検出対象３０の移動状況に応じて、より柔軟な処理を実現することができる。
第３の実施形態によれば、上記の第１の実施形態と同様の効果を奏するものに加え、新御具処理システム６０は、ＲＸ２０Ｂの信号変換部２２により生成されたディジタル信号に対する信号処理を実施することにより、検出対象３０の移動状況に応じて、より柔軟な処理を実現することができる。
以上説明した少なくともひとつの実施形態によれば、レーダシステムは、受信局と、補正部と、推定部とを持つ。受信局は、送信局から送信された電波が検出対象において反射されて、前記反射された電波の周波数成分を含む受信信号を受信する。補正部は、前記送信局が送信した電波の周波数と、前記受信信号に含まれる前記反射された電波の周波数との差として生じたドップラ遷移を打ち消す方向に、前記受信信号の周波数成分を補正する。推定部は、移動する前記検出対象から前記受信局までの距離を、前記補正部により補正された前記受信信号の周波数成分に基づいて推定することにより、検出対象までの距離の検出精度を高めることができる。
本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。
１、１Ａ…レーダシステム、１０…送信レーダ（ＴＸ、送信局）、２０、２０Ａ…受信レーダ（ＲＸ、受信局）、２２…信号変換部、２３…復調処理部（補正部）、２６…パルス積分処理部、２７…通信部、２９…補正部、５０、５０Ａ…演算処理部、５１、５１Ａ…通信部、５２…選択部、５３…インコヒーレント積分処理部、５４…閾値判定処理部、５５、５５Ａ…推定処理部、５６…出力処理部、６０…信号処理システム
Claims (10)
- 送信局から送信された電波が検出対象において反射されて、前記反射された電波の周波数成分を含む受信信号を受信する受信局と、
前記送信局が送信した電波の周波数と、前記受信信号に含まれる前記反射された電波の周波数との差として生じたドップラ遷移を打ち消す方向に、前記受信信号の周波数成分を補正する補正部と、
移動する前記検出対象から前記受信局までの距離を、前記補正部により補正された前記受信信号の周波数成分に基づいて推定する推定部と
を備えるレーダシステム。 - 前記推定部は、
前記受信局に対する相対移動速度が異なる前記検出対象から前記受信局までの距離を推定する
請求項１記載のレーダシステム。 - 前記補正部は、
前記ドップラ遷移を打ち消すように調整された復調係数を生成する復調係数生成部と、
前記復調係数に基づいて前記受信信号の周波数成分を変換するとともに、前記復調係数に従ったパルス圧縮処理を実施する変換部と、
を備える
請求項１または請求項２記載のレーダシステム。 - 複数の前記受信局を備え、
前記推定部は、
前記複数の受信局によってそれぞれ受信され、前記補正部によりそれぞれ補正された前記受信信号の周波数成分に基づいて、少なくとも前記検出対象に対する第１の受信局の測角方向と第２の前記受信局の測角方向とが異なる前記検出対象から前記第１の受信局までの距離と前記第２の受信局までの距離とを推定する
請求項１から請求項３の何れか１項に記載のレーダシステム。 - 複数の前記受信局を備え、
前記推定部は、
前記複数の受信局によってそれぞれ受信され、前記補正部によりそれぞれ補正された前記受信信号の周波数成分に基づいて、少なくとも前記検出対象の第１の前記受信局に対する相対移動速度と前記検出対象の第２の前記受信局に対する相対移動速度とが異なる前記検出対象から前記第１の受信局までの距離と前記第２の受信局までの距離とを推定する
請求項１から請求項３の何れか１項に記載のレーダシステム。 - 前記複数の受信局として、前記第１の受信局と前記第２の受信局と第３の前記受信局とを備え、
前記推定部は、
前記第１の受信局と前記第２の受信局と第３の前記受信局によってそれぞれ受信され、前記補正部によりそれぞれ補正された前記受信信号の周波数成分に基づいて、前記検出対象から前記第１の受信局までの距離と前記第２の受信局までの距離と前記第３の受信局までの距離とを推定する
請求項４又は請求項５記載のレーダシステム。 - 前記補正部は、補正程度を異ならせた複数の補正結果を出力し、
前記推定部は、
前記複数の受信局によってそれぞれ受信され、前記補正部によりそれぞれ補正された前記受信信号のうちの何れかの受信信号を、前記受信信号の信号強度に基づいて前記受信局ごとに選択する選択部と、
前記受信局ごとに前記選択された信号をインコヒーレント積分する積分処理部と、
前記インコヒーレント積分された後の信号から前記検出対象の位置を推定する推定部と、
を備える請求項４から請求項６の何れか１項に記載のレーダシステム。 - 送信局から送信された電波の周波数と、検出対象において反射されて受信局により受信された受信信号に含まれる前記反射された電波の周波数との差として生じたドップラ遷移を打ち消す方向に、前記受信信号の周波数成分を補正する補正部と、
移動する前記検出対象から前記受信局までの距離を、前記補正部により補正された前記受信信号の周波数成分に基づいて推定する推定部と
を備える信号処理システム。 - 送信局から送信された電波の周波数と、検出対象において反射されて受信局により受信された受信信号に含まれる前記反射された電波の周波数との差として生じたドップラ遷移を打ち消す方向に、前記受信信号の周波数成分を補正するステップと、
移動する前記検出対象から前記受信局までの距離を、前記補正された前記受信信号の周波数成分に基づいて推定するステップと
を含む信号処理方法。 - 送信局から送信された電波の周波数と、検出対象において反射されて受信局により受信された受信信号に含まれる前記反射された電波の周波数との差として生じたドップラ遷移を打ち消す方向に、前記受信信号の周波数成分を補正するステップと、
移動する前記検出対象から前記受信局までの距離を、前記補正された前記受信信号の周波数成分に基づいて推定するステップと
をコンピュータに実行させるためのプログラム。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2015177821A JP6716223B2 (ja) | 2015-09-09 | 2015-09-09 | レーダシステム、および信号処理システム |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2015177821A JP6716223B2 (ja) | 2015-09-09 | 2015-09-09 | レーダシステム、および信号処理システム |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2017053721A true JP2017053721A (ja) | 2017-03-16 |
| JP6716223B2 JP6716223B2 (ja) | 2020-07-01 |
Family
ID=58317753
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2015177821A Active JP6716223B2 (ja) | 2015-09-09 | 2015-09-09 | レーダシステム、および信号処理システム |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP6716223B2 (ja) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPWO2018038128A1 (ja) * | 2016-08-26 | 2019-06-24 | 日本電気株式会社 | 移動目標探知システム及び移動目標探知方法 |
| JP2020027047A (ja) * | 2018-08-13 | 2020-02-20 | 株式会社東芝 | レーダシステム及びそのレーダ信号処理方法 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0566268A (ja) * | 1991-09-06 | 1993-03-19 | Mitsubishi Electric Corp | デイジタル・パルス圧縮装置 |
| US5410314A (en) * | 1993-11-30 | 1995-04-25 | University Corporation For Atmospheric Research | Bistatic multiple-doppler radar network |
| JPH0875846A (ja) * | 1994-09-02 | 1996-03-22 | Mitsubishi Heavy Ind Ltd | レーダ装置 |
| JP2651054B2 (ja) * | 1990-03-07 | 1997-09-10 | エイチイー・ホールディングス・インコーポレーテッド・ディービーエー・ヒューズ・エレクトロニクス | ポリスタティック相関レーダ |
| JP2003227870A (ja) * | 2002-02-05 | 2003-08-15 | Mitsubishi Electric Corp | レーダシステム及びターゲット検出方法 |
| JP2012255667A (ja) * | 2011-06-07 | 2012-12-27 | Nippon Soken Inc | 物体検出装置 |
-
2015
- 2015-09-09 JP JP2015177821A patent/JP6716223B2/ja active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2651054B2 (ja) * | 1990-03-07 | 1997-09-10 | エイチイー・ホールディングス・インコーポレーテッド・ディービーエー・ヒューズ・エレクトロニクス | ポリスタティック相関レーダ |
| JPH0566268A (ja) * | 1991-09-06 | 1993-03-19 | Mitsubishi Electric Corp | デイジタル・パルス圧縮装置 |
| US5410314A (en) * | 1993-11-30 | 1995-04-25 | University Corporation For Atmospheric Research | Bistatic multiple-doppler radar network |
| JPH0875846A (ja) * | 1994-09-02 | 1996-03-22 | Mitsubishi Heavy Ind Ltd | レーダ装置 |
| JP2003227870A (ja) * | 2002-02-05 | 2003-08-15 | Mitsubishi Electric Corp | レーダシステム及びターゲット検出方法 |
| JP2012255667A (ja) * | 2011-06-07 | 2012-12-27 | Nippon Soken Inc | 物体検出装置 |
Non-Patent Citations (1)
| Title |
|---|
| 東海林 慶一 他２名: "マルチスタティックレーダによる標定処理精度に関する検討 The Localization Accuracy for the Multistati", 電子情報通信学会２０１５年通信ソサイエティ大会講演論文集１ PROCEEDINGS OF THE 2015 IEICE COMMUNICAT, vol. B-2-19, JPN6019012411, 25 August 2015 (2015-08-25), JP, pages 177, ISSN: 0004013225 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPWO2018038128A1 (ja) * | 2016-08-26 | 2019-06-24 | 日本電気株式会社 | 移動目標探知システム及び移動目標探知方法 |
| US11125870B2 (en) | 2016-08-26 | 2021-09-21 | Nec Corporation | Moving-target detection system and moving-target detection method |
| JP2020027047A (ja) * | 2018-08-13 | 2020-02-20 | 株式会社東芝 | レーダシステム及びそのレーダ信号処理方法 |
| JP7166101B2 (ja) | 2018-08-13 | 2022-11-07 | 株式会社東芝 | レーダシステム及びそのレーダ信号処理方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| JP6716223B2 (ja) | 2020-07-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109407093B (zh) | 解决宽孔径雷达到达角模糊度的多普勒测量 | |
| CN107831480B (zh) | 弹载雷达和差通道稳健自适应杂波抑制方法 | |
| US10175348B2 (en) | Use of range-rate measurements in a fusion tracking system via projections | |
| DK3123197T3 (en) | METHODS AND APPARATUS FOR DETERMINING THE ARRANGEMENT (AOA) OF A RADAR WARNING RECEIVER | |
| US20150247924A1 (en) | Radar apparatus | |
| Goh et al. | A weighted measurement fusion Kalman filter implementation for UAV navigation | |
| EP3494404B1 (en) | System and method for detecting heading and velocity of a target object | |
| US10613197B2 (en) | Antenna specification estimation device and radar device | |
| CN102156279A (zh) | 基于mimo的双基地雷达地面动目标检测方法 | |
| CN111948619B (zh) | 一种多类型外辐射源照射下目标被动协同探测方法及系统 | |
| EP3324205A1 (en) | Decentralised radar system | |
| CN110488293B (zh) | 一种非均匀空间构型的分布式sar系统 | |
| JP6352837B2 (ja) | レーダシステム及びそのレーダ信号処理方法 | |
| US20180203108A1 (en) | Iterative approach to achieve angular ambiguity resolution | |
| CN108776342A (zh) | 一种高速平台sar慢速动目标检测与速度估计方法 | |
| CN115128594A (zh) | 用于雷达跟踪中的速度估计的部分学习模型 | |
| Zhang et al. | Weak target detection within the nonhomogeneous ionospheric clutter background of HFSWR based on STAP | |
| Yang et al. | Joint position and velocity estimation of a moving target in multistatic radar by bistatic range, TDOA, and Doppler shifts | |
| CN110133646B (zh) | 基于nlcs成像的双基前视sar的多通道两脉冲杂波对消方法 | |
| JP6716223B2 (ja) | レーダシステム、および信号処理システム | |
| Page et al. | Detection and tracking of moving vehicles with Gotcha radar systems | |
| CN109521418B (zh) | 基于干涉场的地基雷达测角方法 | |
| Zhu et al. | Robust moving targets detection and velocity estimation using multi-channel and multi-look SAR images | |
| KR102158740B1 (ko) | 레이다 DoA 추정 시스템 및 방법 | |
| JP2001272464A (ja) | レーダ装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A711 | Notification of change in applicant |
Free format text: JAPANESE INTERMEDIATE CODE: A711 Effective date: 20170912 Free format text: JAPANESE INTERMEDIATE CODE: A712 Effective date: 20170912 |
|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20180521 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20190327 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20190409 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20190606 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20191112 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20191129 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20200512 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20200610 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 6716223 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |