目標測位装置
【課題】狭帯な受信機帯域で低速信号処理且つ高分解能を有し、送信局と受信局間の時刻の同期が必要なく、小口径のアンテナを備える小型の受信局からなり、運用が柔軟で低コストの目標測位装置を提供する。
【解決手段】目標測位装置は、送信局はFMCW電波を放射し、3つ以上の受信局は、FMCW電波と同様なローカル信号で直接波、または目標反射波をミキシングして直接ビート信号または目標反射ビート信号を求め、直接ビート信号の出力が最大になる直接波周波数を算出するとともに目標反射ビート信号の出力が最大になる目標反射波周波数を算出し、算出された直接波周波数と目標反射波周波数とに基づいて送信局と受信局との距離和と相対速度和とを推定し、目標位置・速度推定器は、距離和と相対速度和とから目標の3次元位置座標と3次元相対速度成分とを推定する。
【解決手段】目標測位装置は、送信局はFMCW電波を放射し、3つ以上の受信局は、FMCW電波と同様なローカル信号で直接波、または目標反射波をミキシングして直接ビート信号または目標反射ビート信号を求め、直接ビート信号の出力が最大になる直接波周波数を算出するとともに目標反射ビート信号の出力が最大になる目標反射波周波数を算出し、算出された直接波周波数と目標反射波周波数とに基づいて送信局と受信局との距離和と相対速度和とを推定し、目標位置・速度推定器は、距離和と相対速度和とから目標の3次元位置座標と3次元相対速度成分とを推定する。
【発明の詳細な説明】
【技術分野】
【0001】
この発明は、送信アンテナと受信アンテナとが異なる地点に設置され、目標の物体の位置を測位する目標測位装置に関する。
【背景技術】
【0002】
バイスタティックレーダは、送信アンテナと受信アンテナとを異なる地点に設置し、目標の物体の位置を測位する。そして、目標の物体との距離を計測するために送信局と受信局間の時刻が同期されていなければならないので、ロラン(Long Range Navigation)電波を用いる時刻同期用のパルスを送信する、送信パルス周期、時刻、位相を示す送信パルス情報を生成し送受信する方法が採用されている(例えば、特許文献1、特許文献2参照)。
また、送信局と受信局間の時刻同期を行わず、受信局を複数を備え、直接波と目標反射波との相関処理により目標反射パルスの時間差を求め、複数の時刻差から目標の位置を特定する(例えば、特許文献3参照)。
【0003】
【特許文献1】特開平7−140124号公報
【特許文献2】特開2003−156557号公報
【特許文献3】特開平6−148318号公報
【発明の開示】
【発明が解決しようとする課題】
【0004】
しかし、送信局と受信局間の時刻を同期するためには、同期パルサや送信パルス情報を必要としており、高い距離精度を得るためには、広帯域の受信系と高速の信号処理が必要である。また、受信系において角度を高精度に計測するために開口径の大きい大規模な高価なアンテナが必要となるという問題がある。
また、受信局にて直接波パルスと目標反射パルスの時刻差を高精度に求めるためには、広帯域受信機と高速の相関処理が必要であり装置規模が大きくなるという問題がある。また、ドップラ周波数を利用した高精度な相対速度の計測ができないという問題がある。
【0005】
この発明の目的は、狭帯な受信機帯域で低速信号処理且つ高分解能を有し、送信局と受信局間の時刻の同期が必要なく、小口径のアンテナを備える小型の受信局からなり、運用が柔軟で低コストの目標測位装置を提供することである。
【課題を解決するための手段】
【0006】
この発明に係わる目標測位装置は、電波を放射する送信局、直接にまたは目標物体で反射してから到達する電波を受信する3つ以上の受信局および受信する電波に基づいて目標の位置と速度とを推定する目標位置・速度推定器を備える目標測位装置において、上記送信局は、送信FMCW信号を発生する送信FMCW信号発生器と、電波を上記送信FMCW信号により変調する送信信号変換器と、を備え、上記受信局は、上記FMCW信号と同様な掃引傾斜および掃引時間のローカル信号を発生する受信FMCW信号発生器と、上記送信局から受信される電波と上記ローカル信号とをミキシングしてビート信号を求める受信信号変換器と、上記ビート信号の出力が最大になるビート周波数と時刻誤差を算出する周波数分析手段と、上記ビート周波数と上記時刻誤差とに基づいて上記送信局と上記受信局との距離和と相対速度和とを推定する距離和速度和推定手段と、を備え、上記目標位置・速度推定器は、上記受信局から入力される上記距離和と上記時刻誤差とに基づいて上記目標の3次元位置座標を推定し、上記3次元位置座標と上記相対速度和とに基づいて上記目標の3次元相対速度成分を推定する。
【発明の効果】
【0007】
この発明に係わる目標測位装置の効果は、FMCW信号により変調された電波を用い、受信局の時刻を送信局の時刻に同期する必要がなく、同期パルサや送信パルス情報を必要としないので、狭帯域の受信系と低速の信号処理であっても高精度の位置分解能を実現でき、小型で移動可能な受信局となり、運用の柔軟性を確保できる。また、FMCW信号により変調された電波を用いているので、ドップラ現象を利用して高精度な相対速度の計測ができる。
【発明を実施するための最良の形態】
【0008】
実施の形態1.
図1は、この発明の実施の形態1に係わる目標測位装置の構成図である。図2は、送信局、受信局および目標のXY座標面上の位置を示す図である。図3は、送信局および受信局でのFMCW信号の周波数の変化の様子を示す図である。
この発明の実施の形態1に係わる目標測位装置1は、図1に示すように、1つの送信局2、3つの受信局3a、3b、3cおよび目標位置・速度推定器4を備える。各受信局3a〜3cにおいて求められる距離和および相対速度和は一般的な通信手段を用いて目標位置・速度推定器4に送られる。なお、実施の形態1に係わる目標測位装置1は、3つの受信局3を備えるが、4つ以上の受信局3を備えてもよい。
【0009】
送信局2、受信局3a〜3cおよび目標20は、図2に示すように、配置されている。送信局2の座標をXYZ座標系の原点に置くと受信局3a〜3cの位置ベクトルはXnになる。nは受信局3a〜3cを識別する番号であり、3つの受信局3a〜3cが備えられているので、n=1、2、3である。この位置ベクトルXnは、既知である。目標20の位置ベクトルXは、この目標測位装置1により求める未知の値である。
また、図2におけるR0は、送信局2と目標20間の距離、Rnは、目標20と識別番号nの受信局3間の距離、R0nは、送信局2と識別番号nの受信局3間の距離である。
また、v0は、送信局2と目標20との相対速度、vnは、目標20と識別番号nの受信局3との相対速度、v0nは、送信局2と識別番号nの受信局3との相対速度である。
【0010】
送信局2は、図1に示すように、送信Frequency Modulated Continuous Wave信号(以下、FMCW信号と略記する。)を発生する送信FMCW信号発生器5、送信FMCW信号発生器5から入力される送信FMCW信号に従って電波を変調する送信信号変換器6、電波を空中に放射する送信アンテナ7を備える。
送信FMCW信号発生器5は、送信信号変換器6において行われる周波数変調処理の周波数を指定する送信FMCW信号を発生する。送信FMCW信号は、図3(a)に示すように、中心周波数f、掃引帯域幅B、掃引時間TCであり、掃引傾斜が掃引時間TC毎にμ1=B/TCとμ2=−B/TCとに切り替わる。
【0011】
送信信号変換器6は、送信FMCW信号により指定される周波数の連続波に変調して送信アンテナ7に出力する。そして、放射される電波Txiは、式(1)で表される。iは、掃引傾斜μ1または掃引傾斜μ2を識別する番号であり、1または2である。なお、説明を理解しやすくするためにエンベロープ1の一定振幅とした。
【0012】
【数1】
【0013】
各受信局3a〜3cは、図1に示すように、空中を伝搬して到達する電波を受信する受信アンテナ8、受信FMCW信号に基づいてローカル信号を発生する受信FMCW信号発生器9、受信する電波にローカル信号をミキシングし、ローパスフィルターにより周波数の和信号を除去し、ビート信号を生成する受信信号変換器10、ビート信号を処理して距離和と相対速度和とを推定するFMCW信号処理器11を備える。
【0014】
受信アンテナ8には、送信アンテナ7から空中に放射された電波が直接または目標20の物体で反射してから到達する。以下の説明において目標20の物体で反射してから到達する電波を目標反射波、直接に到達する電波を直接波と称す。
目標反射波Rx1,nは、式(2)で表される。τ0は、時刻t=0での送信局2と目標20間の電波の伝搬時間であり、送信局2と目標20間の距離R0を用いると式(3)で表される。τnは、時刻t=0での目標20と識別番号nの受信局3間の電波の伝搬時間であり、目標20と識別番号nの受信局3間の距離Rnを用いると式(4)で表される。v0は、時刻t=0での送信局2と目標20との相対速度、vnは、時刻t=0での目標20と識別番号nの受信局3との相対速度である。cは光速である。
【0015】
【数2】
【0016】
また、直接波Rx0i,nは、式(5)で表される。なお、直接波を目標反射波と区別する方法は,時分割して目標20がない時に計測したり、アンテナビームスキャンによる利得差を計測したりするなどにより区別する。なお、同時に観測されたとした場合も、送信局2と受信局3間の相対速度がある程度既知(たとえばともに静止しているなど)とすることで、目標反射波と直接波を分離可能である。以下の説明では、別個に受信したものとして説明する。
【0017】
【数3】
【0018】
受信FMCW信号発生器9は、受信した電波を受信信号変換器10においてミキシングするためのローカル信号を生成する受信FMCW信号を生成し、その受信FMCW信号を用いてローカル信号を生成して受信信号変換器10に出力する。受信FMCW信号は、図3(b)に示すように、送信FMCW信号と同様に中心周波数f、掃引帯域幅B、掃引時間TCであり、掃引傾斜が掃引時間TC毎にμ1=B/TCとμ2=−B/TCとに切り替わる。但し、送信局2と受信局3との間には時刻ずれδtnがある。
受信FMCW信号発生器9が生成するローカル信号Li,nは、式(6)で表される。
【0019】
【数4】
【0020】
受信信号変換器10は、目標反射波Rxi,nをローカル信号Li,nでミキシンングし、ローパスフィルターを通過して周波数の和信号が除去され、目標反射ビート信号Bi,nが得られる。目標反射ビート信号Bi,nは、式(7)で表される。なお、光速c比べて速度v0、vnは十分小さいと仮定している。
また、直接波Rx0i,nをローカル信号Li,nでミキシンングし、ローパスフィルターを通過して周波数の和信号が除去され、直接ビート信号B0i,nが得られる。直接ビート信号B0i,nは、式(8)で表される。
【0021】
【数5】
【0022】
FMCW信号処理器11は、図1に示すように、周波数分析手段13としての目標反射波周波数分析手段14、直接波周波数分析手段15および時刻誤差推定手段16と距離和速度和推定手段17とを有する。
目標反射波周波数分析手段14は、目標反射ビート信号Bi,nをフーリエ変換し、出力がピークになる目標反射ビート周波数fi,n(バー)を、式(9)から求める。
【0023】
【数6】
【0024】
直接波周波数分析手段15は、直接ビート信号B0i,nをフーリエ変換し、出力がピークになる直接ビート周波数f0i,n(バー)を、式(10)から求める。
【0025】
【数7】
【0026】
時刻誤差推定手段16は、検出された目標反射ビート周波数fi,n(バー)と直接ビート周波数f0i,n(バー)を用いて、式(11)から送信局2と各受信局3a〜3bの間の時刻誤差δtn(ハット)を求める。
【0027】
【数8】
【0028】
距離和速度和推定手段17は、送信局2と各受信局3a〜3bの間の時刻誤差δtn(ハット)を用いて、式(12)から送信局2と目標20および目標20と各受信局3a〜3c間の距離和ξn(ハット)を求める。
式(12)において、周波数μ1、μ2は設定値であり、目標反射ビート周波数f1,n(バー)、目標反射ビート周波数f2,n(バー)、時刻誤差δtn(ハット)はフーリエ変換により求められている。そして、送信局2と各受信局3間の距離R0nは既知であるので、時間遅延τ0nは既知となり、式(12)から距離和ξn(ハット)が求められる。
【0029】
【数9】
【0030】
また、距離和速度和推定手段17は、時刻誤差δtn(ハット)、距離和ξn(ハット)および目標反射ビート周波数f1,n(バー)を用いて、各受信局3a〜3cでの送信局2と目標20および目標20と各受信局3a〜3c間の相対速度和Vn(ハット)を式(13)から求める。
【0031】
【数10】
【0032】
目標位置・速度推定器4は、以下の手順により目標20の3次元座標と3次元相対速度成分を推定する。
既知の送信局2の位置ベクトル、既知の各受信局3の位置ベクトルおよび推定対象の未知の目標20の位置ベクトルは、それぞれX0、Xn、Xであるので、式(12)の距離和ξnを用いることにより、式(14)の関係式が成り立つ。
そして、目標20の位置ベクトルXの推定値X(ハット)は、3つの受信局3a〜3cでの距離和ξn(ハット)を用いて、共役勾配法、準Newton法、Levenberg−Marquardt法などの一般的な非線形手法を用いて求められる。
また、目標20の相対速度ベクトルをVとすると、式(13)の線形方程式は式(15)に書き直せる。
そして、目標20の相対速度ベクトルVの推定値V(ハット)は、3つの受信局3a〜3cでの相対速度和Vn(ハット)と推定値X(ハット)から、逆行列や一般逆行列を用いて求められる。
【0033】
【数11】
【0034】
このような目標測位装置1は、電波としてFMCW電波を用い、3つ以上の受信局を備えており、受信局3の時刻を送信局2の時刻に同期する必要がなく、同期パルサや送信パルス情報を必要としないので、狭帯域の受信系と低速の信号処理系であっても高精度の位置分解能を実現でき、小型で移動可能な受信局3となり、運用の柔軟性を確保できる。また、FMCW電波を用いているので、ドップラ現象を利用して高精度な相対速度の推定ができる。
【0035】
なお、送信局2が送信機能とともに受信機能を有するときには、3つの受信局3の1つを送信局2に負わせることができるので、受信局3は2つでよい。
また、目標20の位置座標および相対速度の推定を水平面座標上の2次元推定を行うときには、受信局3を2つとしてもよい。
また、受信局3が2つのときに、式(13)で求めた2つの距離和ξnを用い、目標位置・速度推定器4では、回転楕円体の交点として目標20の位置を測位することも可能である。すなわち、目標20は送信局2と受信局3とを焦点とした距離和ξnの2つの回転楕円体の交点に存在するが、水平面座標を正とするなどの制約条件を付加することで、2つの回転楕円体の交点の中で目標20の位置として妥当な解を選択することにより目標20の位置座標を推定できる。
【0036】
また、FMCWをパルス化したFMパルス(FMICW(FM Interrupted CW)とも呼ばれる。)とし、距離ゲート化により、更なる狭帯域、直接波と目標反射波の分離、多目標の分離、近距離クラッタの回避などができる。
【0037】
実施の形態2.
図4は、この発明の実施の形態2に係わる目標測位装置の構成図である。図5は、送信局、受信局および目標のXY座標面上の位置を示す図である。
この発明の実施の形態2に係わる目標測位装置1Bは、実施の形態1に係わる目標測位装置1と送信局2Ba〜2Bcと受信局3Bの局数が異なっており、それ以外は同様であるので、同様な部分に同じ符号を付記して説明は省略する。
実施の形態2に係わる目標測位装置1Bは、図4に示すように、3つの送信局2Ba〜2Bcと1つの受信局3Bとを備える。なお、実施の形態2に係わる目標測位装置1Bは、3つの送信局2Ba〜2Bcを備えるが、4つ以上の送信局を備えても同様な効果が得られる。
【0038】
送信局2Ba〜2Bc、受信局3Bおよび目標20は、図5に示すように、配置されている。受信局3Bの座標をXYZ座標系の原点に置くと送信局2Ba〜2Bcの位置ベクトルはXnになる。nは送信局2Ba〜2Bcを識別する番号であり、3つの送信局2Ba〜2Bcが備えられているので、n=1、2、3である。この位置ベクトルXnは、既知である。目標20の位置ベクトルXは、この目標測位装置1Bにより求める未知の値である。
【0039】
各送信局2Ba〜2Bcからは、3つの目標反射波を分離するために、送信を時分割している。なお、送信局2Ba〜2Bc毎に異なる中心周波数のFMCW信号の電波を放射し、受信局3Bの受信信号変換器10でそれぞれのFMCW信号のローカル信号に対しミキシングしてもよい。
【0040】
FMCW信号処理器11は、図4に示すように、周波数分析手段13Bとしての目標反射波周波数分析手段14、直接波周波数分析手段15および時刻誤差推定手段16と距離和速度和推定手段17とを有する。
各送信局2Ba〜2Bcと受信局3BとのFMCW信号の時刻誤差に起因する位相誤差δtnとすると、目標反射ビート信号Bi,nは、式(7)と同様に表される。
そして、目標反射波周波数分析手段14は、実施の形態1と同様にして、目標反射ビート信号Bi,nをフーリエ変換し、出力がピークになる目標反射ビート周波数fi,n(バー)を式(9)から求める。
【0041】
また、直接波周波数分析手段15は、実施の形態1と同様にして、直接ビート信号B0i,nをフーリエ変換し、出力がピークになる直接ビート周波数f0i,n(バー)を式(10)から求める。
【0042】
このようにして求められた目標反射ビート周波数fi,n(バー)、直接ビート周波数f0i,n(バー)を用いて、実施の形態1と同様に、目標20の位置ベクトルの推定値X(ハット)、相対速度ベクトルの推定値V(ハット)を推定することができる。
【0043】
このような目標測位装置1Bは、送信局2を3つ以上備えれば、実施の形態1と同様に、受信局3の時刻を送信局2の時刻に同期する必要がなく、同期パルサや送信パルス情報を必要としないので、狭帯域の受信系と低速の信号処理系であっても高精度の位置分解能を実現でき、小型で移動可能な受信局3となり、運用の柔軟性を確保できる。
なお、送信局数が2であっても、2つの距離和ξnから回転楕円体の交点を求めることで目標20の位置座標を推定可能である。
【0044】
実施の形態3.
図6は、この発明の実施の形態3に係わる目標測位装置の構成図である。図7は、送信局、受信局および目標のXY座標面上の位置を示す図である。
この発明の実施の形態3に係わる目標測位装置1Cは、図6に示すように、実施の形態1に係わる目標測位装置1の3つの受信局3a〜3cの受信信号変換器10に共通の受信FMCW信号発生器9からローカル信号が供給され、それに伴ってFMCW信号処理器11Cが異なっており、それ以外は同様であるので、同様な部分に同じ符号を付記して説明は省略する。
実施の形態3に係わる周波数分析手段13Cは、実施の形態1に係わる周波数分析手段13から直接波周波数分析手段15および時刻誤差推定手段16が削除されている。
【0045】
実施の形態3に係わる目標測位装置1Cでは、図7に示すように、目標反射波だけを受信する。
各受信局3Ca〜3Ccには、1つの受信FMCW信号発生器9から受信FMCW信号がケーブル長の違いによる遅延などを補正して位相が一定値となるよう調整されながらケーブル18で接続されている。
【0046】
この実施の形態3に係わる受信FMCW信号は、各受信局3Ca〜3Cc間で位相が揃っているので、送信局2と各受信局3Ca〜3Ccの時刻誤差δtnは等しく、時刻誤差δtとなる。従って、各受信局3Ca〜3Ccでの受信信号変換器10でミキシングされた後の目標反射ビート信号Bi,nは、式(16)で表される。
【0047】
目標反射波周波数分析手段14Cは、目標反射ビート信号Bi,nをフーリエ変換し、出力がピークになる目標反射ビート周波数fi,n(バー)を式(17)から求める。
距離和速度和推定手段17Cは、検出した目標反射ビート周波数fi,n(バー)を用いて式(18)から送信局2と目標20間と目標20と各受信局3Ca〜3Cc間の時刻誤差δtを含む距離和ξn(ハット)を求める。
また、距離和速度和推定手段17Cは、式(19)から相対速度和Vn(ハット)を求める。
【0048】
【数12】
【0049】
目標位置・速度推定器4Cは、以下の手順により目標20の3次元座標と3次元相対速度成分を推定する。
既知の送信局2の位置ベクトル、既知の各受信局3の位置ベクトルおよび推定対象の未知の目標20の位置ベクトルを、それぞれX0、Xn、Xとすると、式(18)の距離和ξnを用いることにより、時刻誤差δtを消去可能な式(20)〜式(22)が成り立つ。
そして、目標20の位置ベクトルXの推定値X(ハット)は、3つの受信局3Ca〜3Ccでの距離和ξn(ハット)を用いて、共役勾配法、準Newton法、Levenberg−Marquardt法などの一般的な非線形手法を用いて求められる。
【0050】
【数13】
【0051】
また、目標20の相対速度ベクトルをVとすると、式(19)の線形方程式は式(23)に書き直せる。
そして、目標20の相対速度ベクトルVの推定値V(ハット)は、3つの受信局3での相対速度和Vn(ハット)と推定値X(ハット)を用いて、逆行列や一般逆行列を用いて求められる。
【0052】
【数14】
【0053】
このような目標測位装置1Cは、各受信局3Ca〜3Ccへの受信FMCW信号間に時刻誤差の差異が発生しないので、直接波を計測する必要がなくなる。
なお、複数の受信FMCW信号を備え、それらが時刻同期している構成でも同様である。
【0054】
実施の形態4.
図8は、この発明の実施の形態4に係わる目標測位装置の構成図である。図9は、送信局、受信局および目標のXY座標面上の位置を示す図である。
この発明の実施の形態4に係わる目標測位装置1Dは、実施の形態2に係わる目標測位装置1Bと送信FMCW信号が同期されて供給され、それにともなってFMCW信号処理器11Dが異なっており、それ以外は同様であるので、同様な部分に同じ符号を付記して説明は省略する。
実施の形態4に係わる目標測位装置1Dは、3つの送信局2Da〜2Dcと1つの受信局3Dとを備える。なお、実施の形態4に係わる目標測位装置1Dは、3つの送信局2Da〜2Dcを備えるが、4つ以上の送信局を備えても同様な効果が得られる。
各送信局2Da〜2Dcには、1つの送信FMCW信号発生器5から送信FMCW信号がケーブル長の違いによる遅延などを補正して位相が一定値となるよう調整されながらケーブル18で接続されている。また、各送信局2Da〜2Dcからは、目標反射波を分離するために、送信を時分割している。
【0055】
この実施の形態4に係わる送信FMCW信号は、各送信局2Da〜2Dc間で位相が揃っているので、各送信局2Da〜2Dcと受信局3Dの時刻誤差δtnは等しく、時刻誤差δtとなる。従って、受信局3Dの受信信号変換器10でミキシングされた後の目標反射ビート信号Bi,nは、式(16)で表される。
【0056】
実施の形態4に係わる周波数分析手段13Dは、実施の形態2に係わる周波数分析手段13Bから直接波周波数分析手段15および時刻誤差推定手段16が削除されている。
目標反射波周波数分析手段14Dは、目標反射ビート信号Bi,nをフーリエ変換し、出力がピークになる目標反射ビート周波数fi,n(バー)を式(17)から求める。
【0057】
距離和速度和推定手段17Dは、検出した目標反射ビート周波数fi,n(バー)を用いて式(18)から各送信局2と目標20間と目標20と受信局3間の時刻誤差δtを含む距離和ξn(ハット)を求める。
また、距離和速度和推定手段17Dは、式(19)から相対速度和Vn(ハット)を求める。
【0058】
目標位置・速度推定器4Dは、以下の手順により目標20の3次元座標と3次元相対速度成分を推定する。
既知の受信局3Dの位置ベクトル、既知の各送信局2Da〜2Dcの位置ベクトルおよび推定対象の未知の目標20の位置ベクトルを、それぞれX0、Xn、Xとすると、式(18)の距離和ξnを用いることにより、時刻誤差δtが消去可能な式(20)〜式(22)が成り立つ。
そして、目標20の位置ベクトルXの推定値X(ハット)は、3つの送信局2Da〜2Dcからの距離和ξn(ハット)を用いて、共役勾配法、準Newton法、Levenberg−Marquardt法などの一般的な非線形手法を用いて求められる。
【0059】
また、目標20の相対速度ベクトルをVとすると、式(19)の線形方程式は式(23)に書き直せる。
そして、目標20の相対速度ベクトルVの推定値V(ハット)は、3つの送信局2Da〜2Dcでの相対速度和Vn(ハット)と推定値X(ハット)を用いて、逆行列や一般逆行列を用いて求められる。
【0060】
このような目標測位装置1Dは、各送信局2Da〜2Dcへの送信FMCW信号間に時刻誤差の差異が発生しないので、直接波を計測する必要がなくなる。
なお、複数の送信FMCW信号を備え、それらが時刻同期している構成でも同様である。
【0061】
実施の形態5.
図10は、この発明の実施の形態5に係わる目標測位装置の構成図である。図11は、目標である送信局および受信局のXY座標面上の位置を示す図である。
この発明の実施の形態5に係わる目標測位装置1Eは、1つの送信局2、3つの受信局3Ea〜3Ecおよび目標位置・速度推定器4Eを備え、送信局2の位置と相対速度とを推定する。各受信局3Ea〜3Ecにて求められる距離和および相対速度和は通信手段を用いて目標位置・速度推定器4Eに送られる。なお、実施の形態5に係わる目標測位装置1Eは、3つの受信局3Ea〜3Ecを備えるが、4つ以上の受信局を備えても同様な効果が得られる。
【0062】
また、実施の形態5に係わる目標測位装置1Eでは、3つの受信局3Ea〜3Ecの受信信号変換器10に共通の受信FMCW信号発生器9から受信FMCW信号が供給されている。各受信局3Ea〜3Ecには、1つの受信FMCW信号発生器9から受信FMCW信号がケーブル長の違いによる遅延などを補正して位相が一定値となるよう調整されながらケーブル18で接続されている。この実施の形態5に係わる受信FMCW信号は、各受信局3Ea〜3Ec間で位相が揃っているので、送信局2と各受信局3Ea〜3Ecの時刻誤差δtnは等しく、時刻誤差δtとなる。
【0063】
送信局2は、実施の形態1に係わる送信局2と同様であるので、説明は省略する。そして、放射される電波Txiは、式(24)で表される。
【0064】
【数15】
【0065】
受信局3Ea〜3Ecは、それぞれ、空中を伝搬して到達する電波を受信する受信アンテナ8、受信する電波に受信FMCW信号をミキシングし、ローパスフィルターにより周波数の和信号を除去し、ビート信号を生成する受信信号変換器10、ビート信号を処理して距離和と相対速度和とを推定するFMCW信号処理器11Eを備える。また、受信局3Eaは、受信FMCW信号を発生する受信FMCW信号発生器9を備え、受信信号変換器10にはケーブル18が接続されている。
【0066】
受信アンテナ8には、送信アンテナ7から空中に放射された電波が直接に到達する。以下の説明において、直接に到達する電波を直接波と称す。直接波Rx0i,nは、式(25)で表される。なお、τ0nは、送信局と受信局間を電波が伝搬する伝搬時間であり、送信局と受信局間の距離R0nと式(26)の関係式が成り立つ。
【0067】
【数16】
【0068】
受信FMCW信号発生器9は、受信した直接波を受信信号変換器10においてミキシングするためのローカル信号を生成する受信FMCW信号を生成し、その受信FMCW信号を用いてローカル信号を生成して受信信号変換器10に出力する。
受信FMCW信号発生器9が生成するローカル信号Li,nは、式(27)で表される。
【0069】
【数17】
【0070】
受信信号変換器10は、直接波Rx0i,nをローカル信号Li,nでミキシンングし、ローパスフィルターを通過して周波数の和信号が除去され、直接ビート信号B0i,nが得られる。直接ビート信号B0i,nは、式(28)で表される。
【0071】
【数18】
【0072】
FMCW信号処理器11Eは、周波数分析手段13Eとしての直接波周波数分析手段15および距離和速度和推定手段17Eを有する。
直接波周波数分析手段15Eは、直接ビート信号B0i,nをフーリエ変換し、出力がピークになる直接ビート周波数f0i,n(バー)を式(29)から求める。
【0073】
【数19】
【0074】
距離和速度和推定手段17Eは、検出した直接ビート周波数f0i,n(バー)を用いて、式(30)から送信局2と目標20間と目標20と各受信局3間の時刻誤差δtを含む距離和ξn(ハット)を求める。
また、距離和速度和推定手段17Eは、式(31)から相対速度和Vn(ハット)を求める。
【0075】
【数20】
【0076】
目標位置・速度推定器4Eは、以下の手順により目標である送信局2の3次元座標と3次元相対速度成分とを推定する。
推定対象の未知の送信局2Eの位置ベクトルおよび既知の各受信局3の位置ベクトルを、それぞれX、Xnとすると、式(30)の距離和ξnを用いることにより、時刻誤差δtが消去される式(32)〜式(34)の関係式が成り立つ。
そして、送信局2の位置ベクトルXの推定値X(ハット)は、3局の受信局3Ea〜3Ecでの距離和ξn(ハット)を用いて、共役勾配法、準Newton法、Levenberg−Marquardt法などの一般的な非線形手法を用いて求められる。
【0077】
【数21】
【0078】
また、送信局2の相対速度ベクトルをVとすると、相対速度ベクトルVは、式(35)の関係式が成り立つ。
そして、送信局2の相対速度ベクトルVの推定値V(ハット)は、3つの受信局3Ea〜3Ecでの相対速度和Vn(ハット)と推定値X(ハット)を用いて、逆行列や一般逆行列を用いて求められる。
【0079】
【数22】
【0080】
このような目標測位装置1Eは、推定する対象の目標が送信局2である運用形態においても、実施の形態1と同様に、受信局3Ea〜3Ecの時刻を送信局2の時刻に同期する必要がなく、同期パルサや送信パルス情報を必要としないので、狭帯域の受信系と低速の信号処理系であっても高精度の位置分解能を実現でき、小型で移動可能な受信局3Ea〜3Ecとなり、運用の柔軟性を確保できる。
【0081】
実施の形態6.
図12は、この発明の実施の形態6に係わる目標測位装置の構成図である。図13は、目標である受信局および送信局のXY座標面上の位置を示す図である。
この発明の実施の形態6に係わる目標測位装置1Fは、3つの送信局2Fa〜2Fc、1つの受信局3Fおよび目標位置・速度推定器4Fを備え、受信局3Fの位置と相対速度とを推定する。受信局3Fにて求められる距離和および相対速度和は通信手段を用いて目標位置・速度推定器4Fに送られる。なお、実施の形態6に係わる目標測位装置1Fは、3つの送信局2Fa〜2Fcを備えるが、4つ以上の送信局を備えても同様な効果が得られる。
【0082】
また、実施の形態6に係わる目標測位装置1Fでは、3つの送信局2Fa〜2Fcの送信信号変換器6に共通の送信FMCW信号発生器5から送信FMCW信号が供給されている。各送信局2Fa〜2Fcには、1つの送信FMCW信号発生器5から送信FMCW信号がケーブル長の違いによる遅延などを補正して位相が一定値となるよう調整されながらケーブル18で接続されている。この実施の形態6に係わる送信FMCW信号は、各送信局2Fa〜2Fc間で位相が揃っているので、送信局2Fa〜2Fcと受信局3Fの時刻誤差δtnは等しく、時刻誤差δtとなる。
【0083】
受信FMCW信号発生器9は、受信した直接波を受信信号変換器10においてミキシングするためのローカル信号を生成する受信FMCW信号を生成し、その受信FMCW信号を用いてローカル信号を生成して受信信号変換器10に出力する。
受信FMCW信号発生器9が生成するローカル信号Li,nは、式(27)で表される。
【0084】
受信信号変換器10は、直接波Rx0i,nをローカル信号Li,nでミキシンングし、ローパスフィルターを通過して周波数の和信号が除去され、直接ビート信号B0i,nが得られる。直接ビート信号B0i,nは、式(28)で表される。
【0085】
直接波周波数分析手段15Fは、直接ビート信号B0i,nをフーリエ変換し、出力がピークになる直接ビート周波数f0i,n(ハット)を式(29)から求める。
【0086】
距離和速度和推定手段17Fは、検出した直接ビート周波数f0i,n(ハット)を用いて式(30)から各送信局2Fa〜2Fcと目標20間と目標20と受信局3間の距離の和ξn(ハット)を求める。
また、距離和速度和推定手段17Fは、式(31)から相対速度和Vn(ハット)を求める。
【0087】
目標位置・速度推定器4Fは、以下の手順により受信局3Fの3次元座標と3次元相対速度成分を推定する。
推定対象の未知の受信局3Fの位置ベクトルおよび既知の各送信局2Fa〜2Fcの位置ベクトルを、それぞれX、Xnとすると、式(30)の距離和ξnを用いることにより、時刻誤差δtが消去されて、式(32)〜式(34)の関係式が成り立つ。
そして、受信局3Fの位置ベクトルXの推定値X(ハット)は、3つの送信局2Fa〜2Fcとの距離和ξn(ハット)を用いて、共役勾配法、準Newton法、Levenberg−Marquardt法などの一般的な非線形手法を用いて求められる。
【0088】
また、受信局3Fの相対速度ベクトルをVとすると、相対速度ベクトルVは、式(35)の関係式が成り立つ。
そして、受信局3Fの相対速度ベクトルVの推定値V(ハット)は、3つの送信局2Fa〜2Fcでの相対速度和Vn(ハット)と推定値X(ハット)を用いて、逆行列や一般逆行列を用いて求められる。
【0089】
このような目標測位装置1Fは、位置および速度を推定する目標が受信局3Fであるような運用形態でも、送信局2Fa〜2Fcの時刻を受信局3Fの時刻に同期する必要がなく、同期パルサや送信パルス情報を必要としないので、狭帯域の受信系と低速の信号処理系であっても高精度の位置分解能を実現でき、小型で移動可能な受信局3Fとなり、運用の柔軟性を確保できる。
【0090】
実施の形態7.
図14は、この発明の実施の形態7に係わる目標測位装置の構成図である。図15は、送信局および受信局でのFMCW信号の周波数の変化の様子を示す図である。
この発明の実施の形態7に係わる目標測位装置1Gは、実施の形態1に係わる目標測位装置1とFMCW信号が異なるとともに受信局3Ga〜3Gcで目標反射波と直接波とを分離せずに受信することが異なっており、それ以外は同様であるので、同様な部分に同じ符号を付記して説明は省略する。また、送信局2G、受信局3Ga〜3Gcおよび目標20の位置や相対速度は、図2と同様である。
この実施の形態7に係わる目標測位装置1Gは、受信局3Ga〜3Gcが送信局2Gからの直接波と目標反射波とを分離できずに受信するときの目標測位に適用される。但し、実施の形態7に係わる目標測位装置1Gは、送信局と受信局間の相対速度τ0nが既知のときに適用することができる。
【0091】
送信FMCW信号発生器5Gは、図15(a)に示すように、掃引傾斜μ1、掃引傾斜μ2=−μ1、掃引傾斜零の3つの掃引時間TCが連続する周波数パターンの送信FMCW信号を発生する。
受信FMCW信号発生器9Gは、図15(b)に示すように、送信FMCW信号を時刻ずれδtnだけシフトした周波数パターンの受信FMCW信号を発生し、受信FMCW信号に基づいてローカル信号を生成し、受信信号変換器10に送信する。
受信信号変換器10は、受信した電波の信号をローカル信号でミキシングし、ローパスフィルターにより周波数の和信号を除去し、ビート信号を求める。
周波数分析手段19は、ビート信号をフーリエ変換し、出力がピークになる2つのビート周波数を検出する。この2つのビート周波数は、式(36)と式(37)で表すことができる。しかし、このビート周波数のどちらが直接波または目標反射波に係わるものであるか区別することができない。
【0092】
【数23】
【0093】
そこで、距離和速度和推定手段17Gは、式(38)、式(39)のようにビート周波数の差の絶対値を求めることにより、送信局2Gと目標20の間の伝搬時間τ0と目標20と各受信局3Ga〜3Gcの間の伝搬時間τnとの和(τ0+τn)と、送信局2Gと目標20との相対速度v0と目標20と各受信局3Ga〜3Gcとの相対速度vnとの和(v0+vn)を求めることができる。なお、この運用形態では送信局2Gと受信局3Ga〜3Gc間の相対速度τ0nが既知であり、送信局2Gと受信局3Ga〜3Gcの位置も既知である。
【0094】
【数24】
【0095】
式(38)、式(39)の観測値が得られることで、例えば、In−phaseチャネル(実数のビート信号)のみを計測可能なモノスタティックレーダにおいて、文献(三本著、「同相信号だけを検波するFMCWレーダの目標距離・速度計測方法」、電子情報通信学会、論文誌B、Vol.J82−B、No.12、p.2355−2363、1999年12月)に記載の距離と相対速度のアンビギュイティ(1つの目標時)を取り除く手法により距離和ξn=c(τ0+τn)と相対速度和Vn(ハット)を求めることができる。
そして、実施の形態1と同様に、距離和ξn=c(τ0+τn)と相対速度和Vn(ハット)から目標20の位置座標と相対速度を推定できる。
【0096】
このような目標測位装置1Gは、一部に周波数掃引が行われないパターンをFMCW信号として採用しているので、ビート信号の差の絶対値から距離和と相対速度和を求めることができ、直接波と目標反射波の受信が分離できないときにも目標の位置座標と相対速度を推定することができる。
【0097】
実施の形態8.
図16は、この発明の実施の形態8に係わる目標測位装置の構成図である。
この発明の実施の形態8に係わる目標測位装置1Hは、実施の形態2に係わる目標測位装置1BとFMCW信号が異なるとともに受信局3Hで目標反射波と直接波とを分離せずに受信することが異なっており、それ以外は同様であるので、同様な部分に同じ符号を付記して説明は省略する。また、送信局2Ha〜2Hc、受信局3Hおよび目標20の位置や相対速度は、図4と同様である。
この実施の形態8に係わる目標測位装置1Hは、受信局3Hが送信局2Ha〜2Hcからの直接波と目標反射波とを分離できずに受信するときの目標測位に適用される。但し、実施の形態8に係わる目標測位装置1Hは、送信局2Ha〜2Hcと受信局3H間の相対速度τ0nが既知のときに適用することができる。
【0098】
送信FMCW信号発生器5Hは、図15(a)に示すように、掃引傾斜μ1、掃引傾斜μ2=−μ1、掃引傾斜零の3つの掃引時間TCが連続する周波数パターンの送信FMCW信号を発生する。
受信FMCW信号発生器9Hは、図15(b)に示すように、送信FMCW信号を時刻ずれδtnだけシフトした周波数パターンの受信FMCW信号を発生し、受信FMCW信号に基づいてローカル信号を生成し、受信信号変換器10に送信する。
受信信号変換器10は、受信した電波の信号をローカル信号でミキシングし、ローパスフィルターにより周波数の和信号を除去し、ビート信号を求める。
周波数分析手段19は、ビート信号をフーリエ変換し、出力がピークになる2つのビート周波数を検出する。この2つのビート周波数は、式(36)と式(37)で表すことができる。しかし、このビート周波数のどちらが直接波または目標反射波に係わるものであるか区別することができない。
【0099】
そこで、距離和速度和推定手段17Hは、式(38)、式(39)のようにビート周波数の差の絶対値を求めることにより、送信局2Ha〜2Hcと目標20の間の伝搬時間τnと目標20と受信局3Hの間の伝搬時間τ0との和(τ0+τn)と、送信局2Ha〜2Hcと目標20との相対速度vnと目標20と受信局3Hとの相対速度v0との和(v0+vn)を求めることができる。なお、この運用形態では送信局2Ha〜2Hcと受信局3H間の相対速度τ0nが既知であり、送信局2Ha〜2Hcと受信局3Hの位置も既知である。
【0100】
式(38)、式(39)の観測値が得られることで、例えば、In−phaseチャネル(実数のビート信号)のみを計測可能なモノスタティックレーダにおいて、文献(三本著、「同相信号だけを検波するFMCWレーダの目標距離・速度計測方法」、電子情報通信学会、論文誌B、Vol.J82−B、No.12、p.2355−2363、1999年12月)に記載の距離と相対速度のアンビギュイティ(1つの目標時)を取り除く手法により距離和ξn=c(τ0+τn)と相対速度和Vn(ハット)を求めることができる。
そして、実施の形態1と同様に、距離和ξn=c(τ0+τn)と相対速度和Vn(ハット)から目標20の位置と相対速度を推定できる。
【0101】
このような目標測位装置1Hは、一部に周波数掃引が行われないパターンをFMCW信号として採用しているので、ビート信号の差の絶対値から距離和と相対速度和を求めることができ、直接波と目標反射波の受信が分離できないときにも目標の位置と相対速度を推定することができる。
【0102】
なお、実施の形態1乃至8における距離和速度和推定手段を受信局で角度計測を行うことが可能とした場合の目標位置推定法と組み合わせることもできる。
【図面の簡単な説明】
【0103】
【図1】この発明の実施の形態1に係わる目標測位装置の構成図である。
【図2】送信局、受信局および目標のXY座標面上の位置を示す図である。
【図3】送信局および受信局でのFMCW信号の周波数の変化の様子を示す図である。
【図4】この発明の実施の形態2に係わる目標測位装置の構成図である。
【図5】送信局、受信局および目標のXY座標面上の位置を示す図である。
【図6】この発明の実施の形態3に係わる目標測位装置の構成図である。
【図7】送信局、受信局および目標のXY座標面上の位置を示す図である。
【図8】この発明の実施の形態4に係わる目標測位装置の構成図である。
【図9】送信局、受信局および目標のXY座標面上の位置を示す図である。
【図10】この発明の実施の形態5に係わる目標測位装置の構成図である。
【図11】目標である送信局および受信局のXY座標面上の位置を示す図である。
【図12】この発明の実施の形態6に係わる目標測位装置の構成図である。
【図13】目標である受信局および送信局のXY座標面上の位置を示す図である。
【図14】この発明の実施の形態7に係わる目標測位装置の構成図である。
【図15】送信局および受信局でのFMCW信号の周波数の変化の様子を示す図である。
【図16】この発明の実施の形態8に係わる目標測位装置の構成図である。
【符号の説明】
【0104】
1 目標測位装置、2 送信局、3 受信局、4 目標位置・速度推定器、5 送信FMCW信号発生器、6 送信信号変換器、7 送信アンテナ、8 受信アンテナ、9 受信FMCW信号発生器、10 受信信号変換器、11 FMCW信号処理器、13、19 周波数分析手段、14 目標反射波周波数分析手段、15 直接波周波数分析手段、16 時刻誤差推定手段、17 距離和速度和推定手段、18 ケーブル、20 目標。
【技術分野】
【0001】
この発明は、送信アンテナと受信アンテナとが異なる地点に設置され、目標の物体の位置を測位する目標測位装置に関する。
【背景技術】
【0002】
バイスタティックレーダは、送信アンテナと受信アンテナとを異なる地点に設置し、目標の物体の位置を測位する。そして、目標の物体との距離を計測するために送信局と受信局間の時刻が同期されていなければならないので、ロラン(Long Range Navigation)電波を用いる時刻同期用のパルスを送信する、送信パルス周期、時刻、位相を示す送信パルス情報を生成し送受信する方法が採用されている(例えば、特許文献1、特許文献2参照)。
また、送信局と受信局間の時刻同期を行わず、受信局を複数を備え、直接波と目標反射波との相関処理により目標反射パルスの時間差を求め、複数の時刻差から目標の位置を特定する(例えば、特許文献3参照)。
【0003】
【特許文献1】特開平7−140124号公報
【特許文献2】特開2003−156557号公報
【特許文献3】特開平6−148318号公報
【発明の開示】
【発明が解決しようとする課題】
【0004】
しかし、送信局と受信局間の時刻を同期するためには、同期パルサや送信パルス情報を必要としており、高い距離精度を得るためには、広帯域の受信系と高速の信号処理が必要である。また、受信系において角度を高精度に計測するために開口径の大きい大規模な高価なアンテナが必要となるという問題がある。
また、受信局にて直接波パルスと目標反射パルスの時刻差を高精度に求めるためには、広帯域受信機と高速の相関処理が必要であり装置規模が大きくなるという問題がある。また、ドップラ周波数を利用した高精度な相対速度の計測ができないという問題がある。
【0005】
この発明の目的は、狭帯な受信機帯域で低速信号処理且つ高分解能を有し、送信局と受信局間の時刻の同期が必要なく、小口径のアンテナを備える小型の受信局からなり、運用が柔軟で低コストの目標測位装置を提供することである。
【課題を解決するための手段】
【0006】
この発明に係わる目標測位装置は、電波を放射する送信局、直接にまたは目標物体で反射してから到達する電波を受信する3つ以上の受信局および受信する電波に基づいて目標の位置と速度とを推定する目標位置・速度推定器を備える目標測位装置において、上記送信局は、送信FMCW信号を発生する送信FMCW信号発生器と、電波を上記送信FMCW信号により変調する送信信号変換器と、を備え、上記受信局は、上記FMCW信号と同様な掃引傾斜および掃引時間のローカル信号を発生する受信FMCW信号発生器と、上記送信局から受信される電波と上記ローカル信号とをミキシングしてビート信号を求める受信信号変換器と、上記ビート信号の出力が最大になるビート周波数と時刻誤差を算出する周波数分析手段と、上記ビート周波数と上記時刻誤差とに基づいて上記送信局と上記受信局との距離和と相対速度和とを推定する距離和速度和推定手段と、を備え、上記目標位置・速度推定器は、上記受信局から入力される上記距離和と上記時刻誤差とに基づいて上記目標の3次元位置座標を推定し、上記3次元位置座標と上記相対速度和とに基づいて上記目標の3次元相対速度成分を推定する。
【発明の効果】
【0007】
この発明に係わる目標測位装置の効果は、FMCW信号により変調された電波を用い、受信局の時刻を送信局の時刻に同期する必要がなく、同期パルサや送信パルス情報を必要としないので、狭帯域の受信系と低速の信号処理であっても高精度の位置分解能を実現でき、小型で移動可能な受信局となり、運用の柔軟性を確保できる。また、FMCW信号により変調された電波を用いているので、ドップラ現象を利用して高精度な相対速度の計測ができる。
【発明を実施するための最良の形態】
【0008】
実施の形態1.
図1は、この発明の実施の形態1に係わる目標測位装置の構成図である。図2は、送信局、受信局および目標のXY座標面上の位置を示す図である。図3は、送信局および受信局でのFMCW信号の周波数の変化の様子を示す図である。
この発明の実施の形態1に係わる目標測位装置1は、図1に示すように、1つの送信局2、3つの受信局3a、3b、3cおよび目標位置・速度推定器4を備える。各受信局3a〜3cにおいて求められる距離和および相対速度和は一般的な通信手段を用いて目標位置・速度推定器4に送られる。なお、実施の形態1に係わる目標測位装置1は、3つの受信局3を備えるが、4つ以上の受信局3を備えてもよい。
【0009】
送信局2、受信局3a〜3cおよび目標20は、図2に示すように、配置されている。送信局2の座標をXYZ座標系の原点に置くと受信局3a〜3cの位置ベクトルはXnになる。nは受信局3a〜3cを識別する番号であり、3つの受信局3a〜3cが備えられているので、n=1、2、3である。この位置ベクトルXnは、既知である。目標20の位置ベクトルXは、この目標測位装置1により求める未知の値である。
また、図2におけるR0は、送信局2と目標20間の距離、Rnは、目標20と識別番号nの受信局3間の距離、R0nは、送信局2と識別番号nの受信局3間の距離である。
また、v0は、送信局2と目標20との相対速度、vnは、目標20と識別番号nの受信局3との相対速度、v0nは、送信局2と識別番号nの受信局3との相対速度である。
【0010】
送信局2は、図1に示すように、送信Frequency Modulated Continuous Wave信号(以下、FMCW信号と略記する。)を発生する送信FMCW信号発生器5、送信FMCW信号発生器5から入力される送信FMCW信号に従って電波を変調する送信信号変換器6、電波を空中に放射する送信アンテナ7を備える。
送信FMCW信号発生器5は、送信信号変換器6において行われる周波数変調処理の周波数を指定する送信FMCW信号を発生する。送信FMCW信号は、図3(a)に示すように、中心周波数f、掃引帯域幅B、掃引時間TCであり、掃引傾斜が掃引時間TC毎にμ1=B/TCとμ2=−B/TCとに切り替わる。
【0011】
送信信号変換器6は、送信FMCW信号により指定される周波数の連続波に変調して送信アンテナ7に出力する。そして、放射される電波Txiは、式(1)で表される。iは、掃引傾斜μ1または掃引傾斜μ2を識別する番号であり、1または2である。なお、説明を理解しやすくするためにエンベロープ1の一定振幅とした。
【0012】
【数1】
【0013】
各受信局3a〜3cは、図1に示すように、空中を伝搬して到達する電波を受信する受信アンテナ8、受信FMCW信号に基づいてローカル信号を発生する受信FMCW信号発生器9、受信する電波にローカル信号をミキシングし、ローパスフィルターにより周波数の和信号を除去し、ビート信号を生成する受信信号変換器10、ビート信号を処理して距離和と相対速度和とを推定するFMCW信号処理器11を備える。
【0014】
受信アンテナ8には、送信アンテナ7から空中に放射された電波が直接または目標20の物体で反射してから到達する。以下の説明において目標20の物体で反射してから到達する電波を目標反射波、直接に到達する電波を直接波と称す。
目標反射波Rx1,nは、式(2)で表される。τ0は、時刻t=0での送信局2と目標20間の電波の伝搬時間であり、送信局2と目標20間の距離R0を用いると式(3)で表される。τnは、時刻t=0での目標20と識別番号nの受信局3間の電波の伝搬時間であり、目標20と識別番号nの受信局3間の距離Rnを用いると式(4)で表される。v0は、時刻t=0での送信局2と目標20との相対速度、vnは、時刻t=0での目標20と識別番号nの受信局3との相対速度である。cは光速である。
【0015】
【数2】
【0016】
また、直接波Rx0i,nは、式(5)で表される。なお、直接波を目標反射波と区別する方法は,時分割して目標20がない時に計測したり、アンテナビームスキャンによる利得差を計測したりするなどにより区別する。なお、同時に観測されたとした場合も、送信局2と受信局3間の相対速度がある程度既知(たとえばともに静止しているなど)とすることで、目標反射波と直接波を分離可能である。以下の説明では、別個に受信したものとして説明する。
【0017】
【数3】
【0018】
受信FMCW信号発生器9は、受信した電波を受信信号変換器10においてミキシングするためのローカル信号を生成する受信FMCW信号を生成し、その受信FMCW信号を用いてローカル信号を生成して受信信号変換器10に出力する。受信FMCW信号は、図3(b)に示すように、送信FMCW信号と同様に中心周波数f、掃引帯域幅B、掃引時間TCであり、掃引傾斜が掃引時間TC毎にμ1=B/TCとμ2=−B/TCとに切り替わる。但し、送信局2と受信局3との間には時刻ずれδtnがある。
受信FMCW信号発生器9が生成するローカル信号Li,nは、式(6)で表される。
【0019】
【数4】
【0020】
受信信号変換器10は、目標反射波Rxi,nをローカル信号Li,nでミキシンングし、ローパスフィルターを通過して周波数の和信号が除去され、目標反射ビート信号Bi,nが得られる。目標反射ビート信号Bi,nは、式(7)で表される。なお、光速c比べて速度v0、vnは十分小さいと仮定している。
また、直接波Rx0i,nをローカル信号Li,nでミキシンングし、ローパスフィルターを通過して周波数の和信号が除去され、直接ビート信号B0i,nが得られる。直接ビート信号B0i,nは、式(8)で表される。
【0021】
【数5】
【0022】
FMCW信号処理器11は、図1に示すように、周波数分析手段13としての目標反射波周波数分析手段14、直接波周波数分析手段15および時刻誤差推定手段16と距離和速度和推定手段17とを有する。
目標反射波周波数分析手段14は、目標反射ビート信号Bi,nをフーリエ変換し、出力がピークになる目標反射ビート周波数fi,n(バー)を、式(9)から求める。
【0023】
【数6】
【0024】
直接波周波数分析手段15は、直接ビート信号B0i,nをフーリエ変換し、出力がピークになる直接ビート周波数f0i,n(バー)を、式(10)から求める。
【0025】
【数7】
【0026】
時刻誤差推定手段16は、検出された目標反射ビート周波数fi,n(バー)と直接ビート周波数f0i,n(バー)を用いて、式(11)から送信局2と各受信局3a〜3bの間の時刻誤差δtn(ハット)を求める。
【0027】
【数8】
【0028】
距離和速度和推定手段17は、送信局2と各受信局3a〜3bの間の時刻誤差δtn(ハット)を用いて、式(12)から送信局2と目標20および目標20と各受信局3a〜3c間の距離和ξn(ハット)を求める。
式(12)において、周波数μ1、μ2は設定値であり、目標反射ビート周波数f1,n(バー)、目標反射ビート周波数f2,n(バー)、時刻誤差δtn(ハット)はフーリエ変換により求められている。そして、送信局2と各受信局3間の距離R0nは既知であるので、時間遅延τ0nは既知となり、式(12)から距離和ξn(ハット)が求められる。
【0029】
【数9】
【0030】
また、距離和速度和推定手段17は、時刻誤差δtn(ハット)、距離和ξn(ハット)および目標反射ビート周波数f1,n(バー)を用いて、各受信局3a〜3cでの送信局2と目標20および目標20と各受信局3a〜3c間の相対速度和Vn(ハット)を式(13)から求める。
【0031】
【数10】
【0032】
目標位置・速度推定器4は、以下の手順により目標20の3次元座標と3次元相対速度成分を推定する。
既知の送信局2の位置ベクトル、既知の各受信局3の位置ベクトルおよび推定対象の未知の目標20の位置ベクトルは、それぞれX0、Xn、Xであるので、式(12)の距離和ξnを用いることにより、式(14)の関係式が成り立つ。
そして、目標20の位置ベクトルXの推定値X(ハット)は、3つの受信局3a〜3cでの距離和ξn(ハット)を用いて、共役勾配法、準Newton法、Levenberg−Marquardt法などの一般的な非線形手法を用いて求められる。
また、目標20の相対速度ベクトルをVとすると、式(13)の線形方程式は式(15)に書き直せる。
そして、目標20の相対速度ベクトルVの推定値V(ハット)は、3つの受信局3a〜3cでの相対速度和Vn(ハット)と推定値X(ハット)から、逆行列や一般逆行列を用いて求められる。
【0033】
【数11】
【0034】
このような目標測位装置1は、電波としてFMCW電波を用い、3つ以上の受信局を備えており、受信局3の時刻を送信局2の時刻に同期する必要がなく、同期パルサや送信パルス情報を必要としないので、狭帯域の受信系と低速の信号処理系であっても高精度の位置分解能を実現でき、小型で移動可能な受信局3となり、運用の柔軟性を確保できる。また、FMCW電波を用いているので、ドップラ現象を利用して高精度な相対速度の推定ができる。
【0035】
なお、送信局2が送信機能とともに受信機能を有するときには、3つの受信局3の1つを送信局2に負わせることができるので、受信局3は2つでよい。
また、目標20の位置座標および相対速度の推定を水平面座標上の2次元推定を行うときには、受信局3を2つとしてもよい。
また、受信局3が2つのときに、式(13)で求めた2つの距離和ξnを用い、目標位置・速度推定器4では、回転楕円体の交点として目標20の位置を測位することも可能である。すなわち、目標20は送信局2と受信局3とを焦点とした距離和ξnの2つの回転楕円体の交点に存在するが、水平面座標を正とするなどの制約条件を付加することで、2つの回転楕円体の交点の中で目標20の位置として妥当な解を選択することにより目標20の位置座標を推定できる。
【0036】
また、FMCWをパルス化したFMパルス(FMICW(FM Interrupted CW)とも呼ばれる。)とし、距離ゲート化により、更なる狭帯域、直接波と目標反射波の分離、多目標の分離、近距離クラッタの回避などができる。
【0037】
実施の形態2.
図4は、この発明の実施の形態2に係わる目標測位装置の構成図である。図5は、送信局、受信局および目標のXY座標面上の位置を示す図である。
この発明の実施の形態2に係わる目標測位装置1Bは、実施の形態1に係わる目標測位装置1と送信局2Ba〜2Bcと受信局3Bの局数が異なっており、それ以外は同様であるので、同様な部分に同じ符号を付記して説明は省略する。
実施の形態2に係わる目標測位装置1Bは、図4に示すように、3つの送信局2Ba〜2Bcと1つの受信局3Bとを備える。なお、実施の形態2に係わる目標測位装置1Bは、3つの送信局2Ba〜2Bcを備えるが、4つ以上の送信局を備えても同様な効果が得られる。
【0038】
送信局2Ba〜2Bc、受信局3Bおよび目標20は、図5に示すように、配置されている。受信局3Bの座標をXYZ座標系の原点に置くと送信局2Ba〜2Bcの位置ベクトルはXnになる。nは送信局2Ba〜2Bcを識別する番号であり、3つの送信局2Ba〜2Bcが備えられているので、n=1、2、3である。この位置ベクトルXnは、既知である。目標20の位置ベクトルXは、この目標測位装置1Bにより求める未知の値である。
【0039】
各送信局2Ba〜2Bcからは、3つの目標反射波を分離するために、送信を時分割している。なお、送信局2Ba〜2Bc毎に異なる中心周波数のFMCW信号の電波を放射し、受信局3Bの受信信号変換器10でそれぞれのFMCW信号のローカル信号に対しミキシングしてもよい。
【0040】
FMCW信号処理器11は、図4に示すように、周波数分析手段13Bとしての目標反射波周波数分析手段14、直接波周波数分析手段15および時刻誤差推定手段16と距離和速度和推定手段17とを有する。
各送信局2Ba〜2Bcと受信局3BとのFMCW信号の時刻誤差に起因する位相誤差δtnとすると、目標反射ビート信号Bi,nは、式(7)と同様に表される。
そして、目標反射波周波数分析手段14は、実施の形態1と同様にして、目標反射ビート信号Bi,nをフーリエ変換し、出力がピークになる目標反射ビート周波数fi,n(バー)を式(9)から求める。
【0041】
また、直接波周波数分析手段15は、実施の形態1と同様にして、直接ビート信号B0i,nをフーリエ変換し、出力がピークになる直接ビート周波数f0i,n(バー)を式(10)から求める。
【0042】
このようにして求められた目標反射ビート周波数fi,n(バー)、直接ビート周波数f0i,n(バー)を用いて、実施の形態1と同様に、目標20の位置ベクトルの推定値X(ハット)、相対速度ベクトルの推定値V(ハット)を推定することができる。
【0043】
このような目標測位装置1Bは、送信局2を3つ以上備えれば、実施の形態1と同様に、受信局3の時刻を送信局2の時刻に同期する必要がなく、同期パルサや送信パルス情報を必要としないので、狭帯域の受信系と低速の信号処理系であっても高精度の位置分解能を実現でき、小型で移動可能な受信局3となり、運用の柔軟性を確保できる。
なお、送信局数が2であっても、2つの距離和ξnから回転楕円体の交点を求めることで目標20の位置座標を推定可能である。
【0044】
実施の形態3.
図6は、この発明の実施の形態3に係わる目標測位装置の構成図である。図7は、送信局、受信局および目標のXY座標面上の位置を示す図である。
この発明の実施の形態3に係わる目標測位装置1Cは、図6に示すように、実施の形態1に係わる目標測位装置1の3つの受信局3a〜3cの受信信号変換器10に共通の受信FMCW信号発生器9からローカル信号が供給され、それに伴ってFMCW信号処理器11Cが異なっており、それ以外は同様であるので、同様な部分に同じ符号を付記して説明は省略する。
実施の形態3に係わる周波数分析手段13Cは、実施の形態1に係わる周波数分析手段13から直接波周波数分析手段15および時刻誤差推定手段16が削除されている。
【0045】
実施の形態3に係わる目標測位装置1Cでは、図7に示すように、目標反射波だけを受信する。
各受信局3Ca〜3Ccには、1つの受信FMCW信号発生器9から受信FMCW信号がケーブル長の違いによる遅延などを補正して位相が一定値となるよう調整されながらケーブル18で接続されている。
【0046】
この実施の形態3に係わる受信FMCW信号は、各受信局3Ca〜3Cc間で位相が揃っているので、送信局2と各受信局3Ca〜3Ccの時刻誤差δtnは等しく、時刻誤差δtとなる。従って、各受信局3Ca〜3Ccでの受信信号変換器10でミキシングされた後の目標反射ビート信号Bi,nは、式(16)で表される。
【0047】
目標反射波周波数分析手段14Cは、目標反射ビート信号Bi,nをフーリエ変換し、出力がピークになる目標反射ビート周波数fi,n(バー)を式(17)から求める。
距離和速度和推定手段17Cは、検出した目標反射ビート周波数fi,n(バー)を用いて式(18)から送信局2と目標20間と目標20と各受信局3Ca〜3Cc間の時刻誤差δtを含む距離和ξn(ハット)を求める。
また、距離和速度和推定手段17Cは、式(19)から相対速度和Vn(ハット)を求める。
【0048】
【数12】
【0049】
目標位置・速度推定器4Cは、以下の手順により目標20の3次元座標と3次元相対速度成分を推定する。
既知の送信局2の位置ベクトル、既知の各受信局3の位置ベクトルおよび推定対象の未知の目標20の位置ベクトルを、それぞれX0、Xn、Xとすると、式(18)の距離和ξnを用いることにより、時刻誤差δtを消去可能な式(20)〜式(22)が成り立つ。
そして、目標20の位置ベクトルXの推定値X(ハット)は、3つの受信局3Ca〜3Ccでの距離和ξn(ハット)を用いて、共役勾配法、準Newton法、Levenberg−Marquardt法などの一般的な非線形手法を用いて求められる。
【0050】
【数13】
【0051】
また、目標20の相対速度ベクトルをVとすると、式(19)の線形方程式は式(23)に書き直せる。
そして、目標20の相対速度ベクトルVの推定値V(ハット)は、3つの受信局3での相対速度和Vn(ハット)と推定値X(ハット)を用いて、逆行列や一般逆行列を用いて求められる。
【0052】
【数14】
【0053】
このような目標測位装置1Cは、各受信局3Ca〜3Ccへの受信FMCW信号間に時刻誤差の差異が発生しないので、直接波を計測する必要がなくなる。
なお、複数の受信FMCW信号を備え、それらが時刻同期している構成でも同様である。
【0054】
実施の形態4.
図8は、この発明の実施の形態4に係わる目標測位装置の構成図である。図9は、送信局、受信局および目標のXY座標面上の位置を示す図である。
この発明の実施の形態4に係わる目標測位装置1Dは、実施の形態2に係わる目標測位装置1Bと送信FMCW信号が同期されて供給され、それにともなってFMCW信号処理器11Dが異なっており、それ以外は同様であるので、同様な部分に同じ符号を付記して説明は省略する。
実施の形態4に係わる目標測位装置1Dは、3つの送信局2Da〜2Dcと1つの受信局3Dとを備える。なお、実施の形態4に係わる目標測位装置1Dは、3つの送信局2Da〜2Dcを備えるが、4つ以上の送信局を備えても同様な効果が得られる。
各送信局2Da〜2Dcには、1つの送信FMCW信号発生器5から送信FMCW信号がケーブル長の違いによる遅延などを補正して位相が一定値となるよう調整されながらケーブル18で接続されている。また、各送信局2Da〜2Dcからは、目標反射波を分離するために、送信を時分割している。
【0055】
この実施の形態4に係わる送信FMCW信号は、各送信局2Da〜2Dc間で位相が揃っているので、各送信局2Da〜2Dcと受信局3Dの時刻誤差δtnは等しく、時刻誤差δtとなる。従って、受信局3Dの受信信号変換器10でミキシングされた後の目標反射ビート信号Bi,nは、式(16)で表される。
【0056】
実施の形態4に係わる周波数分析手段13Dは、実施の形態2に係わる周波数分析手段13Bから直接波周波数分析手段15および時刻誤差推定手段16が削除されている。
目標反射波周波数分析手段14Dは、目標反射ビート信号Bi,nをフーリエ変換し、出力がピークになる目標反射ビート周波数fi,n(バー)を式(17)から求める。
【0057】
距離和速度和推定手段17Dは、検出した目標反射ビート周波数fi,n(バー)を用いて式(18)から各送信局2と目標20間と目標20と受信局3間の時刻誤差δtを含む距離和ξn(ハット)を求める。
また、距離和速度和推定手段17Dは、式(19)から相対速度和Vn(ハット)を求める。
【0058】
目標位置・速度推定器4Dは、以下の手順により目標20の3次元座標と3次元相対速度成分を推定する。
既知の受信局3Dの位置ベクトル、既知の各送信局2Da〜2Dcの位置ベクトルおよび推定対象の未知の目標20の位置ベクトルを、それぞれX0、Xn、Xとすると、式(18)の距離和ξnを用いることにより、時刻誤差δtが消去可能な式(20)〜式(22)が成り立つ。
そして、目標20の位置ベクトルXの推定値X(ハット)は、3つの送信局2Da〜2Dcからの距離和ξn(ハット)を用いて、共役勾配法、準Newton法、Levenberg−Marquardt法などの一般的な非線形手法を用いて求められる。
【0059】
また、目標20の相対速度ベクトルをVとすると、式(19)の線形方程式は式(23)に書き直せる。
そして、目標20の相対速度ベクトルVの推定値V(ハット)は、3つの送信局2Da〜2Dcでの相対速度和Vn(ハット)と推定値X(ハット)を用いて、逆行列や一般逆行列を用いて求められる。
【0060】
このような目標測位装置1Dは、各送信局2Da〜2Dcへの送信FMCW信号間に時刻誤差の差異が発生しないので、直接波を計測する必要がなくなる。
なお、複数の送信FMCW信号を備え、それらが時刻同期している構成でも同様である。
【0061】
実施の形態5.
図10は、この発明の実施の形態5に係わる目標測位装置の構成図である。図11は、目標である送信局および受信局のXY座標面上の位置を示す図である。
この発明の実施の形態5に係わる目標測位装置1Eは、1つの送信局2、3つの受信局3Ea〜3Ecおよび目標位置・速度推定器4Eを備え、送信局2の位置と相対速度とを推定する。各受信局3Ea〜3Ecにて求められる距離和および相対速度和は通信手段を用いて目標位置・速度推定器4Eに送られる。なお、実施の形態5に係わる目標測位装置1Eは、3つの受信局3Ea〜3Ecを備えるが、4つ以上の受信局を備えても同様な効果が得られる。
【0062】
また、実施の形態5に係わる目標測位装置1Eでは、3つの受信局3Ea〜3Ecの受信信号変換器10に共通の受信FMCW信号発生器9から受信FMCW信号が供給されている。各受信局3Ea〜3Ecには、1つの受信FMCW信号発生器9から受信FMCW信号がケーブル長の違いによる遅延などを補正して位相が一定値となるよう調整されながらケーブル18で接続されている。この実施の形態5に係わる受信FMCW信号は、各受信局3Ea〜3Ec間で位相が揃っているので、送信局2と各受信局3Ea〜3Ecの時刻誤差δtnは等しく、時刻誤差δtとなる。
【0063】
送信局2は、実施の形態1に係わる送信局2と同様であるので、説明は省略する。そして、放射される電波Txiは、式(24)で表される。
【0064】
【数15】
【0065】
受信局3Ea〜3Ecは、それぞれ、空中を伝搬して到達する電波を受信する受信アンテナ8、受信する電波に受信FMCW信号をミキシングし、ローパスフィルターにより周波数の和信号を除去し、ビート信号を生成する受信信号変換器10、ビート信号を処理して距離和と相対速度和とを推定するFMCW信号処理器11Eを備える。また、受信局3Eaは、受信FMCW信号を発生する受信FMCW信号発生器9を備え、受信信号変換器10にはケーブル18が接続されている。
【0066】
受信アンテナ8には、送信アンテナ7から空中に放射された電波が直接に到達する。以下の説明において、直接に到達する電波を直接波と称す。直接波Rx0i,nは、式(25)で表される。なお、τ0nは、送信局と受信局間を電波が伝搬する伝搬時間であり、送信局と受信局間の距離R0nと式(26)の関係式が成り立つ。
【0067】
【数16】
【0068】
受信FMCW信号発生器9は、受信した直接波を受信信号変換器10においてミキシングするためのローカル信号を生成する受信FMCW信号を生成し、その受信FMCW信号を用いてローカル信号を生成して受信信号変換器10に出力する。
受信FMCW信号発生器9が生成するローカル信号Li,nは、式(27)で表される。
【0069】
【数17】
【0070】
受信信号変換器10は、直接波Rx0i,nをローカル信号Li,nでミキシンングし、ローパスフィルターを通過して周波数の和信号が除去され、直接ビート信号B0i,nが得られる。直接ビート信号B0i,nは、式(28)で表される。
【0071】
【数18】
【0072】
FMCW信号処理器11Eは、周波数分析手段13Eとしての直接波周波数分析手段15および距離和速度和推定手段17Eを有する。
直接波周波数分析手段15Eは、直接ビート信号B0i,nをフーリエ変換し、出力がピークになる直接ビート周波数f0i,n(バー)を式(29)から求める。
【0073】
【数19】
【0074】
距離和速度和推定手段17Eは、検出した直接ビート周波数f0i,n(バー)を用いて、式(30)から送信局2と目標20間と目標20と各受信局3間の時刻誤差δtを含む距離和ξn(ハット)を求める。
また、距離和速度和推定手段17Eは、式(31)から相対速度和Vn(ハット)を求める。
【0075】
【数20】
【0076】
目標位置・速度推定器4Eは、以下の手順により目標である送信局2の3次元座標と3次元相対速度成分とを推定する。
推定対象の未知の送信局2Eの位置ベクトルおよび既知の各受信局3の位置ベクトルを、それぞれX、Xnとすると、式(30)の距離和ξnを用いることにより、時刻誤差δtが消去される式(32)〜式(34)の関係式が成り立つ。
そして、送信局2の位置ベクトルXの推定値X(ハット)は、3局の受信局3Ea〜3Ecでの距離和ξn(ハット)を用いて、共役勾配法、準Newton法、Levenberg−Marquardt法などの一般的な非線形手法を用いて求められる。
【0077】
【数21】
【0078】
また、送信局2の相対速度ベクトルをVとすると、相対速度ベクトルVは、式(35)の関係式が成り立つ。
そして、送信局2の相対速度ベクトルVの推定値V(ハット)は、3つの受信局3Ea〜3Ecでの相対速度和Vn(ハット)と推定値X(ハット)を用いて、逆行列や一般逆行列を用いて求められる。
【0079】
【数22】
【0080】
このような目標測位装置1Eは、推定する対象の目標が送信局2である運用形態においても、実施の形態1と同様に、受信局3Ea〜3Ecの時刻を送信局2の時刻に同期する必要がなく、同期パルサや送信パルス情報を必要としないので、狭帯域の受信系と低速の信号処理系であっても高精度の位置分解能を実現でき、小型で移動可能な受信局3Ea〜3Ecとなり、運用の柔軟性を確保できる。
【0081】
実施の形態6.
図12は、この発明の実施の形態6に係わる目標測位装置の構成図である。図13は、目標である受信局および送信局のXY座標面上の位置を示す図である。
この発明の実施の形態6に係わる目標測位装置1Fは、3つの送信局2Fa〜2Fc、1つの受信局3Fおよび目標位置・速度推定器4Fを備え、受信局3Fの位置と相対速度とを推定する。受信局3Fにて求められる距離和および相対速度和は通信手段を用いて目標位置・速度推定器4Fに送られる。なお、実施の形態6に係わる目標測位装置1Fは、3つの送信局2Fa〜2Fcを備えるが、4つ以上の送信局を備えても同様な効果が得られる。
【0082】
また、実施の形態6に係わる目標測位装置1Fでは、3つの送信局2Fa〜2Fcの送信信号変換器6に共通の送信FMCW信号発生器5から送信FMCW信号が供給されている。各送信局2Fa〜2Fcには、1つの送信FMCW信号発生器5から送信FMCW信号がケーブル長の違いによる遅延などを補正して位相が一定値となるよう調整されながらケーブル18で接続されている。この実施の形態6に係わる送信FMCW信号は、各送信局2Fa〜2Fc間で位相が揃っているので、送信局2Fa〜2Fcと受信局3Fの時刻誤差δtnは等しく、時刻誤差δtとなる。
【0083】
受信FMCW信号発生器9は、受信した直接波を受信信号変換器10においてミキシングするためのローカル信号を生成する受信FMCW信号を生成し、その受信FMCW信号を用いてローカル信号を生成して受信信号変換器10に出力する。
受信FMCW信号発生器9が生成するローカル信号Li,nは、式(27)で表される。
【0084】
受信信号変換器10は、直接波Rx0i,nをローカル信号Li,nでミキシンングし、ローパスフィルターを通過して周波数の和信号が除去され、直接ビート信号B0i,nが得られる。直接ビート信号B0i,nは、式(28)で表される。
【0085】
直接波周波数分析手段15Fは、直接ビート信号B0i,nをフーリエ変換し、出力がピークになる直接ビート周波数f0i,n(ハット)を式(29)から求める。
【0086】
距離和速度和推定手段17Fは、検出した直接ビート周波数f0i,n(ハット)を用いて式(30)から各送信局2Fa〜2Fcと目標20間と目標20と受信局3間の距離の和ξn(ハット)を求める。
また、距離和速度和推定手段17Fは、式(31)から相対速度和Vn(ハット)を求める。
【0087】
目標位置・速度推定器4Fは、以下の手順により受信局3Fの3次元座標と3次元相対速度成分を推定する。
推定対象の未知の受信局3Fの位置ベクトルおよび既知の各送信局2Fa〜2Fcの位置ベクトルを、それぞれX、Xnとすると、式(30)の距離和ξnを用いることにより、時刻誤差δtが消去されて、式(32)〜式(34)の関係式が成り立つ。
そして、受信局3Fの位置ベクトルXの推定値X(ハット)は、3つの送信局2Fa〜2Fcとの距離和ξn(ハット)を用いて、共役勾配法、準Newton法、Levenberg−Marquardt法などの一般的な非線形手法を用いて求められる。
【0088】
また、受信局3Fの相対速度ベクトルをVとすると、相対速度ベクトルVは、式(35)の関係式が成り立つ。
そして、受信局3Fの相対速度ベクトルVの推定値V(ハット)は、3つの送信局2Fa〜2Fcでの相対速度和Vn(ハット)と推定値X(ハット)を用いて、逆行列や一般逆行列を用いて求められる。
【0089】
このような目標測位装置1Fは、位置および速度を推定する目標が受信局3Fであるような運用形態でも、送信局2Fa〜2Fcの時刻を受信局3Fの時刻に同期する必要がなく、同期パルサや送信パルス情報を必要としないので、狭帯域の受信系と低速の信号処理系であっても高精度の位置分解能を実現でき、小型で移動可能な受信局3Fとなり、運用の柔軟性を確保できる。
【0090】
実施の形態7.
図14は、この発明の実施の形態7に係わる目標測位装置の構成図である。図15は、送信局および受信局でのFMCW信号の周波数の変化の様子を示す図である。
この発明の実施の形態7に係わる目標測位装置1Gは、実施の形態1に係わる目標測位装置1とFMCW信号が異なるとともに受信局3Ga〜3Gcで目標反射波と直接波とを分離せずに受信することが異なっており、それ以外は同様であるので、同様な部分に同じ符号を付記して説明は省略する。また、送信局2G、受信局3Ga〜3Gcおよび目標20の位置や相対速度は、図2と同様である。
この実施の形態7に係わる目標測位装置1Gは、受信局3Ga〜3Gcが送信局2Gからの直接波と目標反射波とを分離できずに受信するときの目標測位に適用される。但し、実施の形態7に係わる目標測位装置1Gは、送信局と受信局間の相対速度τ0nが既知のときに適用することができる。
【0091】
送信FMCW信号発生器5Gは、図15(a)に示すように、掃引傾斜μ1、掃引傾斜μ2=−μ1、掃引傾斜零の3つの掃引時間TCが連続する周波数パターンの送信FMCW信号を発生する。
受信FMCW信号発生器9Gは、図15(b)に示すように、送信FMCW信号を時刻ずれδtnだけシフトした周波数パターンの受信FMCW信号を発生し、受信FMCW信号に基づいてローカル信号を生成し、受信信号変換器10に送信する。
受信信号変換器10は、受信した電波の信号をローカル信号でミキシングし、ローパスフィルターにより周波数の和信号を除去し、ビート信号を求める。
周波数分析手段19は、ビート信号をフーリエ変換し、出力がピークになる2つのビート周波数を検出する。この2つのビート周波数は、式(36)と式(37)で表すことができる。しかし、このビート周波数のどちらが直接波または目標反射波に係わるものであるか区別することができない。
【0092】
【数23】
【0093】
そこで、距離和速度和推定手段17Gは、式(38)、式(39)のようにビート周波数の差の絶対値を求めることにより、送信局2Gと目標20の間の伝搬時間τ0と目標20と各受信局3Ga〜3Gcの間の伝搬時間τnとの和(τ0+τn)と、送信局2Gと目標20との相対速度v0と目標20と各受信局3Ga〜3Gcとの相対速度vnとの和(v0+vn)を求めることができる。なお、この運用形態では送信局2Gと受信局3Ga〜3Gc間の相対速度τ0nが既知であり、送信局2Gと受信局3Ga〜3Gcの位置も既知である。
【0094】
【数24】
【0095】
式(38)、式(39)の観測値が得られることで、例えば、In−phaseチャネル(実数のビート信号)のみを計測可能なモノスタティックレーダにおいて、文献(三本著、「同相信号だけを検波するFMCWレーダの目標距離・速度計測方法」、電子情報通信学会、論文誌B、Vol.J82−B、No.12、p.2355−2363、1999年12月)に記載の距離と相対速度のアンビギュイティ(1つの目標時)を取り除く手法により距離和ξn=c(τ0+τn)と相対速度和Vn(ハット)を求めることができる。
そして、実施の形態1と同様に、距離和ξn=c(τ0+τn)と相対速度和Vn(ハット)から目標20の位置座標と相対速度を推定できる。
【0096】
このような目標測位装置1Gは、一部に周波数掃引が行われないパターンをFMCW信号として採用しているので、ビート信号の差の絶対値から距離和と相対速度和を求めることができ、直接波と目標反射波の受信が分離できないときにも目標の位置座標と相対速度を推定することができる。
【0097】
実施の形態8.
図16は、この発明の実施の形態8に係わる目標測位装置の構成図である。
この発明の実施の形態8に係わる目標測位装置1Hは、実施の形態2に係わる目標測位装置1BとFMCW信号が異なるとともに受信局3Hで目標反射波と直接波とを分離せずに受信することが異なっており、それ以外は同様であるので、同様な部分に同じ符号を付記して説明は省略する。また、送信局2Ha〜2Hc、受信局3Hおよび目標20の位置や相対速度は、図4と同様である。
この実施の形態8に係わる目標測位装置1Hは、受信局3Hが送信局2Ha〜2Hcからの直接波と目標反射波とを分離できずに受信するときの目標測位に適用される。但し、実施の形態8に係わる目標測位装置1Hは、送信局2Ha〜2Hcと受信局3H間の相対速度τ0nが既知のときに適用することができる。
【0098】
送信FMCW信号発生器5Hは、図15(a)に示すように、掃引傾斜μ1、掃引傾斜μ2=−μ1、掃引傾斜零の3つの掃引時間TCが連続する周波数パターンの送信FMCW信号を発生する。
受信FMCW信号発生器9Hは、図15(b)に示すように、送信FMCW信号を時刻ずれδtnだけシフトした周波数パターンの受信FMCW信号を発生し、受信FMCW信号に基づいてローカル信号を生成し、受信信号変換器10に送信する。
受信信号変換器10は、受信した電波の信号をローカル信号でミキシングし、ローパスフィルターにより周波数の和信号を除去し、ビート信号を求める。
周波数分析手段19は、ビート信号をフーリエ変換し、出力がピークになる2つのビート周波数を検出する。この2つのビート周波数は、式(36)と式(37)で表すことができる。しかし、このビート周波数のどちらが直接波または目標反射波に係わるものであるか区別することができない。
【0099】
そこで、距離和速度和推定手段17Hは、式(38)、式(39)のようにビート周波数の差の絶対値を求めることにより、送信局2Ha〜2Hcと目標20の間の伝搬時間τnと目標20と受信局3Hの間の伝搬時間τ0との和(τ0+τn)と、送信局2Ha〜2Hcと目標20との相対速度vnと目標20と受信局3Hとの相対速度v0との和(v0+vn)を求めることができる。なお、この運用形態では送信局2Ha〜2Hcと受信局3H間の相対速度τ0nが既知であり、送信局2Ha〜2Hcと受信局3Hの位置も既知である。
【0100】
式(38)、式(39)の観測値が得られることで、例えば、In−phaseチャネル(実数のビート信号)のみを計測可能なモノスタティックレーダにおいて、文献(三本著、「同相信号だけを検波するFMCWレーダの目標距離・速度計測方法」、電子情報通信学会、論文誌B、Vol.J82−B、No.12、p.2355−2363、1999年12月)に記載の距離と相対速度のアンビギュイティ(1つの目標時)を取り除く手法により距離和ξn=c(τ0+τn)と相対速度和Vn(ハット)を求めることができる。
そして、実施の形態1と同様に、距離和ξn=c(τ0+τn)と相対速度和Vn(ハット)から目標20の位置と相対速度を推定できる。
【0101】
このような目標測位装置1Hは、一部に周波数掃引が行われないパターンをFMCW信号として採用しているので、ビート信号の差の絶対値から距離和と相対速度和を求めることができ、直接波と目標反射波の受信が分離できないときにも目標の位置と相対速度を推定することができる。
【0102】
なお、実施の形態1乃至8における距離和速度和推定手段を受信局で角度計測を行うことが可能とした場合の目標位置推定法と組み合わせることもできる。
【図面の簡単な説明】
【0103】
【図1】この発明の実施の形態1に係わる目標測位装置の構成図である。
【図2】送信局、受信局および目標のXY座標面上の位置を示す図である。
【図3】送信局および受信局でのFMCW信号の周波数の変化の様子を示す図である。
【図4】この発明の実施の形態2に係わる目標測位装置の構成図である。
【図5】送信局、受信局および目標のXY座標面上の位置を示す図である。
【図6】この発明の実施の形態3に係わる目標測位装置の構成図である。
【図7】送信局、受信局および目標のXY座標面上の位置を示す図である。
【図8】この発明の実施の形態4に係わる目標測位装置の構成図である。
【図9】送信局、受信局および目標のXY座標面上の位置を示す図である。
【図10】この発明の実施の形態5に係わる目標測位装置の構成図である。
【図11】目標である送信局および受信局のXY座標面上の位置を示す図である。
【図12】この発明の実施の形態6に係わる目標測位装置の構成図である。
【図13】目標である受信局および送信局のXY座標面上の位置を示す図である。
【図14】この発明の実施の形態7に係わる目標測位装置の構成図である。
【図15】送信局および受信局でのFMCW信号の周波数の変化の様子を示す図である。
【図16】この発明の実施の形態8に係わる目標測位装置の構成図である。
【符号の説明】
【0104】
1 目標測位装置、2 送信局、3 受信局、4 目標位置・速度推定器、5 送信FMCW信号発生器、6 送信信号変換器、7 送信アンテナ、8 受信アンテナ、9 受信FMCW信号発生器、10 受信信号変換器、11 FMCW信号処理器、13、19 周波数分析手段、14 目標反射波周波数分析手段、15 直接波周波数分析手段、16 時刻誤差推定手段、17 距離和速度和推定手段、18 ケーブル、20 目標。
【特許請求の範囲】
【請求項1】
電波を放射する送信局、直接にまたは目標物体で反射してから到達する電波を受信する3つ以上の受信局および受信する電波に基づいて目標の位置と速度とを推定する目標位置・速度推定器を備える目標測位装置において、
上記送信局は、
送信FMCW信号を発生する送信FMCW信号発生器と、
電波を上記送信FMCW信号により変調する送信信号変換器と、
を備え、
上記受信局は、
上記FMCW信号と同様な掃引傾斜および掃引時間のローカル信号を発生する受信FMCW信号発生器と、
上記送信局から受信される電波と上記ローカル信号とをミキシングしてビート信号を求める受信信号変換器と、
上記ビート信号の出力が最大になるビート周波数と時刻誤差を算出する周波数分析手段と、
上記ビート周波数と上記時刻誤差とに基づいて上記送信局と上記受信局との距離和と相対速度和とを推定する距離和速度和推定手段と、
を備え、
上記目標位置・速度推定器は、上記受信局から入力される上記距離和と上記時刻誤差とに基づいて上記目標の3次元位置座標を推定し、上記3次元位置座標と上記相対速度和とに基づいて上記目標の3次元相対速度成分を推定することを特徴とする目標測位装置。
【請求項2】
上記受信信号変換器は、上記送信局から直接受信される電波、または目標で反射してから受信される電波と上記ローカル信号とをミキシングして直接ビート信号または目標反射ビート信号を求め、
上記周波数分析手段は、上記直接ビート信号の出力が最大になる直接ビート周波数を算出する直接波周波数分析手段と、上記目標反射ビート信号の出力が最大になる目標反射ビート周波数を算出する目標反射波周波数分析手段と、上記直接ビート周波数に基づいて時刻誤差を推定する時刻誤差推定手段とを備え、
上記距離和速度和推定手段は、上記目標反射ビート周波数と上記時刻誤差とに基づいて上記送信局と上記受信局との距離和と相対速度和とを推定することを特徴とする請求項1に記載する目標測位装置。
【請求項3】
各上記受信局の上記受信信号変換器に、他の上記受信局の上記受信信号変換器に供給される上記受信FMCW信号と同期する受信FMCW信号が供給され、上記目標で反射してから受信される電波と上記ローカル信号とをミキシングして目標反射ビート信号を求め、
上記周波数分析手段は、上記目標反射ビート信号の出力が最大になる目標反射ビート周波数を算出する目標反射波周波数分析手段を備え、
上記距離和速度和推定手段は、上記目標反射波周波数分析手段で算出される目標ビート周波数に基づいて上記送信局と上記受信局との時刻誤差を含む距離和と相対速度和とを推定し、
上記目標位置・速度推定器は、上記受信局から入力される上記時刻誤差を含む距離和から時刻誤差を消去することにより上記目標の3次元位置座標を推定し、上記3次元位置座標と上記相対速度和とに基づいて上記目標の3次元相対速度を推定することを特徴とすることを特徴とする請求項1に記載する目標測位装置。
【請求項4】
各上記受信局の上記受信信号変換器に、他の上記受信局の上記受信信号変換器に供給される上記受信FMCW信号と同期する受信FMCW信号が供給され、上記送信局から直接に受信される電波と上記ローカル信号とをミキシングして直接ビート信号を求め、
上記周波数分析手段は、上記直接ビート信号の出力が最大になる直接ビート周波数を算出する直接波周波数分析手段を備え、
上記距離和速度和推定手段は、上記直接波周波数分析手段で算出される直接ビート周波数に基づいて上記送信局と上記受信局との時刻誤差を含む距離和と相対速度和とを推定し、
上記目標位置・速度推定器は、上記受信局から入力される上記時刻誤差を含む距離和から時刻誤差を消去することにより上記送信局の3次元位置座標を推定し、上記3次元位置座標と上記相対速度和とに基づいて上記送信局の3次元相対速度を推定することを特徴とすることを特徴とする請求項1に記載する目標測位装置。
【請求項5】
電波を放射する3つ以上の送信局、直接にまたは目標で反射してから到達する電波を受信する受信局および受信する電波に基づいて目標の位置座標と相対速度とを推定する目標位置・速度推定器を備える目標測位装置において、
上記送信局は、
送信FMCW信号を発生する送信FMCW信号発生器と、
電波を上記送信FMCW信号により変調する送信信号変換器と、
を備え、
上記受信局は、
上記FMCW信号と同様な掃引傾斜および掃引時間のローカル信号を発生する受信FMCW信号発生器と、
上記送信局から受信される電波と上記ローカル信号とをミキシングしてビート信号を求める受信信号変換器と、
上記ビート信号の出力が最大になるビート周波数と時刻誤差を算出する周波数分析手段と、
上記ビート周波数と上記時刻誤差とに基づいて上記送信局と上記受信局との距離和と相対速度和とを推定する距離和速度和推定手段と、
を備え、
上記目標位置・速度推定器は、上記受信局から入力される上記距離和と上記時刻誤差とに基づいて上記目標の3次元位置座標を推定し、上記3次元位置座標と上記相対速度和とに基づいて上記目標の3次元相対速度成分を推定することを特徴とする目標測位装置。
【請求項6】
上記受信信号変換器は、各上記送信局から直接受信される電波、または目標で反射してから受信される電波と上記ローカル信号とをミキシングして直接ビート信号または目標反射ビート信号を求め、
上記周波数分析手段は、上記直接ビート信号の出力が最大になる直接ビート周波数を算出する直接波周波数分析手段と、上記目標反射ビート信号の出力が最大になる目標反射ビート周波数を算出する目標反射波周波数分析手段と、上記直接ビート周波数に基づいて時刻誤差を推定する時刻誤差推定手段とを備え、
上記距離和速度和推定手段は、上記目標反射ビート周波数と上記時刻誤差とに基づいて上記送信局と上記受信局との距離和と相対速度和とを推定することを特徴とする請求項5に記載する目標測位装置。
【請求項7】
各上記送信局の上記送信信号変換器に、他の上記送信局の上記送信信号変換器に供給される上記送信FMCW信号と同期する送信FMCW信号が供給され、
上記受信信号変換器は、上記目標で反射してから受信される電波と上記ローカル信号とをミキシングして目標反射ビート信号を求め、
上記周波数分析手段は、上記目標反射ビート信号の出力が最大になる目標反射ビート周波数を算出する目標反射波周波数分析手段を備え、
上記距離和速度和推定手段は、上記目標反射波周波数分析手段で算出される目標ビート周波数に基づいて上記送信局と上記受信局との時刻誤差を含む距離和と相対速度和とを推定し、
上記目標位置・速度推定器は、上記受信局から入力される上記時刻誤差を含む距離和から時刻誤差を消去することにより上記目標の3次元位置座標を推定し、上記3次元位置座標と上記相対速度和とに基づいて上記目標の3次元相対速度を推定することを特徴とすることを特徴とする請求項5に記載する目標測位装置。
【請求項8】
各上記送信局の上記送信信号変換器に、他の上記送信局の上記送信信号変換器に供給される上記送信FMCW信号と同期する送信FMCW信号が供給され、
上記受信信号変換器は、各上記送信局から直接に受信される電波と上記ローカル信号とをミキシングして直接ビート信号を求め、
上記周波数分析手段は、上記直接ビート信号の出力が最大になる直接ビート周波数を算出する直接波周波数分析手段を備え、
上記距離和速度和推定手段は、上記直接波周波数分析手段で算出される直接ビート周波数に基づいて上記送信局と上記受信局との時刻誤差を含む距離和と相対速度和とを推定し、
上記目標位置・速度推定器は、上記受信局から入力される上記時刻誤差を含む距離和から時刻誤差を消去することにより上記受信局の3次元位置座標を推定し、上記3次元位置座標と上記相対速度和とに基づいて上記受信局の3次元相対速度を推定することを特徴とすることを特徴とする請求項5に記載する目標測位装置。
【請求項9】
上記送信FMCW信号および上記受信FMCW信号は、周波数が一定に維持される時間を含み、
上記送信局と上記受信局との間の距離が既知であり、
上記周波数分析手段は、上記ビート信号から2つのビート周波数を算出し、
上記距離和速度和推定手段は、2つの上記ビート周波数の差の絶対値に基づいて上記送信局と上記受信局との距離和と相対速度和とを推定することを特徴とする請求項1または5に記載する目標測位装置。
【請求項1】
電波を放射する送信局、直接にまたは目標物体で反射してから到達する電波を受信する3つ以上の受信局および受信する電波に基づいて目標の位置と速度とを推定する目標位置・速度推定器を備える目標測位装置において、
上記送信局は、
送信FMCW信号を発生する送信FMCW信号発生器と、
電波を上記送信FMCW信号により変調する送信信号変換器と、
を備え、
上記受信局は、
上記FMCW信号と同様な掃引傾斜および掃引時間のローカル信号を発生する受信FMCW信号発生器と、
上記送信局から受信される電波と上記ローカル信号とをミキシングしてビート信号を求める受信信号変換器と、
上記ビート信号の出力が最大になるビート周波数と時刻誤差を算出する周波数分析手段と、
上記ビート周波数と上記時刻誤差とに基づいて上記送信局と上記受信局との距離和と相対速度和とを推定する距離和速度和推定手段と、
を備え、
上記目標位置・速度推定器は、上記受信局から入力される上記距離和と上記時刻誤差とに基づいて上記目標の3次元位置座標を推定し、上記3次元位置座標と上記相対速度和とに基づいて上記目標の3次元相対速度成分を推定することを特徴とする目標測位装置。
【請求項2】
上記受信信号変換器は、上記送信局から直接受信される電波、または目標で反射してから受信される電波と上記ローカル信号とをミキシングして直接ビート信号または目標反射ビート信号を求め、
上記周波数分析手段は、上記直接ビート信号の出力が最大になる直接ビート周波数を算出する直接波周波数分析手段と、上記目標反射ビート信号の出力が最大になる目標反射ビート周波数を算出する目標反射波周波数分析手段と、上記直接ビート周波数に基づいて時刻誤差を推定する時刻誤差推定手段とを備え、
上記距離和速度和推定手段は、上記目標反射ビート周波数と上記時刻誤差とに基づいて上記送信局と上記受信局との距離和と相対速度和とを推定することを特徴とする請求項1に記載する目標測位装置。
【請求項3】
各上記受信局の上記受信信号変換器に、他の上記受信局の上記受信信号変換器に供給される上記受信FMCW信号と同期する受信FMCW信号が供給され、上記目標で反射してから受信される電波と上記ローカル信号とをミキシングして目標反射ビート信号を求め、
上記周波数分析手段は、上記目標反射ビート信号の出力が最大になる目標反射ビート周波数を算出する目標反射波周波数分析手段を備え、
上記距離和速度和推定手段は、上記目標反射波周波数分析手段で算出される目標ビート周波数に基づいて上記送信局と上記受信局との時刻誤差を含む距離和と相対速度和とを推定し、
上記目標位置・速度推定器は、上記受信局から入力される上記時刻誤差を含む距離和から時刻誤差を消去することにより上記目標の3次元位置座標を推定し、上記3次元位置座標と上記相対速度和とに基づいて上記目標の3次元相対速度を推定することを特徴とすることを特徴とする請求項1に記載する目標測位装置。
【請求項4】
各上記受信局の上記受信信号変換器に、他の上記受信局の上記受信信号変換器に供給される上記受信FMCW信号と同期する受信FMCW信号が供給され、上記送信局から直接に受信される電波と上記ローカル信号とをミキシングして直接ビート信号を求め、
上記周波数分析手段は、上記直接ビート信号の出力が最大になる直接ビート周波数を算出する直接波周波数分析手段を備え、
上記距離和速度和推定手段は、上記直接波周波数分析手段で算出される直接ビート周波数に基づいて上記送信局と上記受信局との時刻誤差を含む距離和と相対速度和とを推定し、
上記目標位置・速度推定器は、上記受信局から入力される上記時刻誤差を含む距離和から時刻誤差を消去することにより上記送信局の3次元位置座標を推定し、上記3次元位置座標と上記相対速度和とに基づいて上記送信局の3次元相対速度を推定することを特徴とすることを特徴とする請求項1に記載する目標測位装置。
【請求項5】
電波を放射する3つ以上の送信局、直接にまたは目標で反射してから到達する電波を受信する受信局および受信する電波に基づいて目標の位置座標と相対速度とを推定する目標位置・速度推定器を備える目標測位装置において、
上記送信局は、
送信FMCW信号を発生する送信FMCW信号発生器と、
電波を上記送信FMCW信号により変調する送信信号変換器と、
を備え、
上記受信局は、
上記FMCW信号と同様な掃引傾斜および掃引時間のローカル信号を発生する受信FMCW信号発生器と、
上記送信局から受信される電波と上記ローカル信号とをミキシングしてビート信号を求める受信信号変換器と、
上記ビート信号の出力が最大になるビート周波数と時刻誤差を算出する周波数分析手段と、
上記ビート周波数と上記時刻誤差とに基づいて上記送信局と上記受信局との距離和と相対速度和とを推定する距離和速度和推定手段と、
を備え、
上記目標位置・速度推定器は、上記受信局から入力される上記距離和と上記時刻誤差とに基づいて上記目標の3次元位置座標を推定し、上記3次元位置座標と上記相対速度和とに基づいて上記目標の3次元相対速度成分を推定することを特徴とする目標測位装置。
【請求項6】
上記受信信号変換器は、各上記送信局から直接受信される電波、または目標で反射してから受信される電波と上記ローカル信号とをミキシングして直接ビート信号または目標反射ビート信号を求め、
上記周波数分析手段は、上記直接ビート信号の出力が最大になる直接ビート周波数を算出する直接波周波数分析手段と、上記目標反射ビート信号の出力が最大になる目標反射ビート周波数を算出する目標反射波周波数分析手段と、上記直接ビート周波数に基づいて時刻誤差を推定する時刻誤差推定手段とを備え、
上記距離和速度和推定手段は、上記目標反射ビート周波数と上記時刻誤差とに基づいて上記送信局と上記受信局との距離和と相対速度和とを推定することを特徴とする請求項5に記載する目標測位装置。
【請求項7】
各上記送信局の上記送信信号変換器に、他の上記送信局の上記送信信号変換器に供給される上記送信FMCW信号と同期する送信FMCW信号が供給され、
上記受信信号変換器は、上記目標で反射してから受信される電波と上記ローカル信号とをミキシングして目標反射ビート信号を求め、
上記周波数分析手段は、上記目標反射ビート信号の出力が最大になる目標反射ビート周波数を算出する目標反射波周波数分析手段を備え、
上記距離和速度和推定手段は、上記目標反射波周波数分析手段で算出される目標ビート周波数に基づいて上記送信局と上記受信局との時刻誤差を含む距離和と相対速度和とを推定し、
上記目標位置・速度推定器は、上記受信局から入力される上記時刻誤差を含む距離和から時刻誤差を消去することにより上記目標の3次元位置座標を推定し、上記3次元位置座標と上記相対速度和とに基づいて上記目標の3次元相対速度を推定することを特徴とすることを特徴とする請求項5に記載する目標測位装置。
【請求項8】
各上記送信局の上記送信信号変換器に、他の上記送信局の上記送信信号変換器に供給される上記送信FMCW信号と同期する送信FMCW信号が供給され、
上記受信信号変換器は、各上記送信局から直接に受信される電波と上記ローカル信号とをミキシングして直接ビート信号を求め、
上記周波数分析手段は、上記直接ビート信号の出力が最大になる直接ビート周波数を算出する直接波周波数分析手段を備え、
上記距離和速度和推定手段は、上記直接波周波数分析手段で算出される直接ビート周波数に基づいて上記送信局と上記受信局との時刻誤差を含む距離和と相対速度和とを推定し、
上記目標位置・速度推定器は、上記受信局から入力される上記時刻誤差を含む距離和から時刻誤差を消去することにより上記受信局の3次元位置座標を推定し、上記3次元位置座標と上記相対速度和とに基づいて上記受信局の3次元相対速度を推定することを特徴とすることを特徴とする請求項5に記載する目標測位装置。
【請求項9】
上記送信FMCW信号および上記受信FMCW信号は、周波数が一定に維持される時間を含み、
上記送信局と上記受信局との間の距離が既知であり、
上記周波数分析手段は、上記ビート信号から2つのビート周波数を算出し、
上記距離和速度和推定手段は、2つの上記ビート周波数の差の絶対値に基づいて上記送信局と上記受信局との距離和と相対速度和とを推定することを特徴とする請求項1または5に記載する目標測位装置。
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
【図14】
【図15】
【図16】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
【図14】
【図15】
【図16】
【公開番号】特開2007−192575(P2007−192575A)
【公開日】平成19年8月2日(2007.8.2)
【国際特許分類】
【出願番号】特願2006−8801(P2006−8801)
【出願日】平成18年1月17日(2006.1.17)
【出願人】(000006013)三菱電機株式会社 (33,312)
【Fターム(参考)】
【公開日】平成19年8月2日(2007.8.2)
【国際特許分類】
【出願日】平成18年1月17日(2006.1.17)
【出願人】(000006013)三菱電機株式会社 (33,312)
【Fターム(参考)】
[ Back to top ]