MIMOレーダは、例えば、時分割、周波数分割又は符号分割を用いて多重した信号(レーダ送信波)を複数の送信アンテナ(又は送信アレーアンテナと呼ぶ)から送信し、周辺物体において反射された信号(レーダ反射波)を複数の受信アンテナ(又は受信アレーアンテナと呼ぶ)を用いて受信し、それぞれの受信信号から、多重された送信信号を分離して受信する。このような処理により、MIMOレーダは、送信アンテナ数と受信アンテナ数との積で示される伝搬路応答を取り出すことができ、これらの受信信号を仮想受信アレーとしてアレー信号処理を行う。
また、MIMOレーダでは、送受信アレーアンテナにおける素子間隔を適切に配置することにより、仮想的にアンテナ開口を拡大し、角度分解能の向上を図ることができる。
例えば、特許文献1には、MIMOレーダの多重送信方法として、送信アンテナ毎に送信時間をずらして信号を送信する時分割多重送信を用いたMIMOレーダ(以下、「時分割多重MIMOレーダ」と呼ぶ)が開示されている。時分割多重送信は、周波数多重送信又は符号多重送信と比較し、簡易な構成で実現できる。また、時分割多重送信は、送信時間の間隔を十分に広げることにより、送信信号間の直交性を良好に保つことができる。時分割多重MIMOレーダは、送信アンテナを所定の周期で逐次的に切り替えながら、送信信号の一例である送信パルスを出力する。時分割多重MIMOレーダは、送信パルスが物体で反射された信号を複数の受信アンテナで受信し、受信信号と送信パルスとの相関処理後に、例えば、空間的なFFT(Fast Fourier Transforma)処理(反射波の到来方向推定処理)を行う。
時分割多重MIMOレーダは、送信信号(例えば送信パルス又はレーダ送信波)を送信する送信アンテナを、所定の周期で逐次的に切り替えていく。したがって、時分割多重送信は、周波数分割送信又は符号分割送信と比較し、全ての送信アンテナから送信信号を送信し終えるまでに要する時間が長くなり得る。このため、例えば、特許文献2のように、各送信アンテナから送信信号を送信し、それらの受信位相変化からドップラ周波数(つまり、ターゲットの相対速度)の検出を行う場合、ドップラ周波数を検出するためにフーリエ周波数解析を適用するにあたり、受信位相変化の観測の時間間隔(例えば、サンプリング間隔)が長くなる。よって、折り返しなしでドップラ周波数を検出できるドップラ周波数範囲(つまり、検出できるターゲットの相対速度範囲)が低減する。
また、折り返しなしでドップラ周波数を検出できるドップラ周波数範囲(換言すると、相対速度範囲)を超えるターゲットからの反射波信号が想定される場合、レーダ装置は、反射波信号が折り返し成分か否かを特定できず、ドップラ周波数(換言すると、ターゲットの相対速度)の曖昧性(不確定性、Ambiguity)が生じる。
例えば、レーダ装置が、Nt個の送信アンテナを所定の周期Trで逐次的に切り替えながら送信信号(送信パルス)を送信する場合、全ての送信アンテナから送信信号を送信し終えるまでにTr×Ntの送信時間が必要となる。このような時分割多重送信をNc回繰り返して、ドップラ周波数の検出のためにフーリエ周波数解析を適用すると、折り返しなしでドップラ周波数を検出できるドップラ周波数範囲は、サンプリング定理より、±1/(2Tr×Nt)となる。したがって、折り返しなしでドップラ周波数を検出できるドップラ周波数範囲は、送信アンテナ数Ntが増大するほど低減し、より低速な相対速度でもドップラ周波数の曖昧性が生じやすくなる。
時分割多重MIMOレーダには上述したようなドップラ周波数の曖昧性が生じる恐れがあるため、以下では、一例として、複数の送信アンテナから送信信号を同時に多重して送信する方法に着目する。
複数の送信アンテナから送信信号を同時に多重して送信する方法として、例えば、受信側においてドップラ周波数軸上で複数の送信信号を分離できるように信号を送信する方法(以下、ドップラ多重送信と呼ぶ)がある(例えば、非特許文献3を参照)。
ドップラ多重送信において、送信側では、例えば、基準となる送信アンテナから送信される送信信号に対して、基準となる送信アンテナと異なる送信アンテナから送信される送信信号に、受信信号のドップラ周波数帯域幅よりも大きなドップラシフト量が与えられ、複数の送信アンテナから送信信号が同時に送信される。ドップラ多重送信において、受信側では、ドップラ周波数軸上でフィルタリングすることにより、各送信アンテナから送信された送信信号が分離して受信される。
ドップラ多重送信では、複数の送信アンテナから送信信号を同時に送信することにより、時分割多重送信と比較して、ドップラ周波数(又は、相対速度)の検出のためにフーリエ周波数解析を適用する際の受信位相変化を観測する時間間隔を短縮できる。しかし、ドップラ多重送信では、ドップラ周波数軸上でフィルタリングすることにより各送信アンテナの送信信号を分離するため、送信信号あたりの実効的なドップラ周波数帯域幅が制限されてしまう。
例えば、ドップラ多重送信において、レーダ装置が、Nt個の送信アンテナから周期Trで送信信号を送信する場合について説明する。このようなドップラ多重送信をNc回繰り返して、ドップラ周波数(又は、相対速度)の検出のためにフーリエ周波数解析を適用すると、折り返しなしでドップラ周波数を検出できるドップラ周波数範囲は、サンプリング定理より±1/(2×Tr)となる。つまり、ドップラ多重送信において折り返しなしでドップラ周波数を検出できるドップラ周波数範囲は、時分割多重送信の場合(例えば、±1/(2Tr×Nt))と比較してNt倍に拡大される。
ただし、ドップラ多重送信では、上述したように、ドップラ周波数軸上でフィルタリングすることによって送信信号が分離される。そのため、送信信号あたりの実効的なドップラ周波数帯域幅は、1/(Tr×Nt)に制限されるので、時分割多重送信を行った場合と同様なドップラ周波数範囲となる。また、ドップラ多重送信において、送信信号あたりの実効的なドップラ周波数範囲を超えたドップラ周波数帯域では、当該送信信号と異なる他の送信信号のドップラ周波数帯域の信号と混在するため、送信信号を正しく分離できない可能性がある。
そこで、本開示に係る一実施例では、ドップラ多重送信において、折り返しが発生しない(換言すると、曖昧性が生じない)ドップラ周波数の範囲を拡大させる方法について説明する。これにより、本開示に係る一実施例のレーダ装置は、より広いドップラ周波数範囲において、物標を精度良く検知できる。
以下、本開示の一実施例に係る実施の形態について、図面を参照して詳細に説明する。なお、実施の形態において、同一の構成要素には同一の符号を付し、その説明は重複するので省略する。
以下では、レーダ装置において、送信ブランチにおいて、複数の送信アンテナから同時に多重された異なる送信信号を送出し、受信ブランチにおいて、各送信信号を分離して受信処理を行う構成(換言すると、MIMOレーダ構成)について説明する。
また、以下では、一例として、チャープ(chirp)パルスのような周波数変調したパルス波を用いたレーダ方式(例えば、チャープパルス送信(fast chirp modulation)とも呼ぶ)の構成について説明する。ただし、変調方式は、周波数変調に限定されない。例えば、本開示の一実施例は、パルス列を位相変調又は振幅変調して送信するパルス圧縮レーダを用いたレーダ方式についても適用可能である。
[レーダ装置の構成]
図1は、本実施の形態に係るレーダ装置10の構成を示すブロック図である。
レーダ装置10は、レーダ送信部(送信ブランチ)100と、レーダ受信部(受信ブランチ)200と、を有する。
レーダ送信部100は、レーダ信号(レーダ送信信号)を生成し、複数の送信アンテナ105-1~105-Ntによって構成される送信アレーアンテナを用いて、レーダ送信信号を所定の送信周期にて送信する。
レーダ受信部200は、物標(ターゲット。図示せず)により反射したレーダ送信信号である反射波信号を、複数の受信アンテナ202-1~202-Naを含む受信アレーアンテナを用いて受信する。レーダ受信部200は、各受信アンテナ202において受信した反射波信号を信号処理し、例えば、物標の有無検出又は反射波信号の到来方向の推定を行う。
なお、物標はレーダ装置10が検出する対象の物体であり、例えば、車両(4輪及び2輪を含む)、人、ブロック又は縁石などを含む。
[レーダ送信部100の構成]
レーダ送信部100は、レーダ送信信号生成部101と、ドップラシフト部104-1~104-Ntと、送信アンテナ105-1~105-Ntと、を有する。すなわち、レーダ送信部100は、Nt個の送信アンテナ105を有し、各送信アンテナ105は、それぞれ個別のドップラシフト部104に接続されている。
レーダ送信信号生成部101は、レーダ送信信号を生成する。レーダ送信信号生成部101は、例えば、変調信号発生部102及びVCO(Voltage Controlled Oscillator:電圧制御発信器)103を有する。以下、レーダ送信信号生成部101における各構成部について説明する。
変調信号発生部102は、例えば、図2に示すように、のこぎり歯形状の変調信号を周期的に発生させる。ここで、レーダ送信周期をTrとする。
VCO103は、変調信号発生部102から出力されるレーダ送信信号に基づいて、周波数変調信号(以下、例えば、周波数チャープ信号又はチャープ信号と呼ぶ)をドップラシフト部104-1~104-Nt、及び、レーダ受信部200(後述するミキサ部204)へ出力する。
ドップラシフト部104は、VCO103から入力されるチャープ信号に対して、ドップラシフト量DOPnを付与するために、位相回転φnを付与し、ドップラシフト後の信号を送信アンテナ105に出力する。ここで、n=1,…,Ntである。なお、ドップラシフト部104におけるドップラシフト量DOPn(換言すると、位相回転φn)を付与する方法の一例については後述する。
ドップラシフト部104-1~104-Ntの出力信号は、所定の送信電力に増幅され各送信アンテナ105から空間に放射される。
[レーダ受信部200の構成]
図1において、レーダ受信部200は、Na個の受信アンテナ202を備え、アレーアンテナを構成する。また、レーダ受信部200は、Na個のアンテナ系統処理部201-1~201-Naと、CFAR(Constant False Alarm Rate)部210と、ドップラ多重分離部211と、方向推定部212と、を有する。
各受信アンテナ202は、物標(ターゲット)に反射したレーダ送信信号である反射波信号を受信し、受信した反射波信号を、対応するアンテナ系統処理部201へ受信信号として出力する。
各アンテナ系統処理部201は、受信無線部203と、信号処理部206とを有する。
受信無線部203は、ミキサ部204と、LPF(low pass filter)205と、を有する。受信無線部203は、ミキサ部204において、受信した反射波信号に対して、送信信号であるチャープ信号をミキシングし、LPF205を通過させる。これにより、反射波信号の遅延時間に応じた周波数となるビート信号が取り出される。例えば、図2に示すように、送信信号(送信周波数変調波)の周波数と、受信信号(受信周波数変調波)の周波数との差分周波数がビート周波数として得られる。
各アンテナ系統処理部201-z(ただし、z=1~Naの何れか)の信号処理部206は、AD変換部207と、ビート周波数解析部208と、ドップラ解析部209と、を有する。
LPF205から出力された信号(例えば、ビート信号)は、信号処理部206において、AD変換部207によって、離散的にサンプリングされた離散サンプルデータに変換される。
ビート周波数解析部208は、送信周期Tr毎に、所定の時間範囲(レンジゲート)において得られたNdata個の離散サンプルデータをFFT処理する。これにより、信号処理部206では、反射波信号(レーダ反射波)の遅延時間に応じたビート周波数にピークが現れる周波数スペクトラムが出力される。なお、FFT処理の際、ビート周波数解析部208は、例えば、Han窓又はHamming窓等の窓関数係数を乗算してもよい。窓関数係数を用いることにより、ビート周波数ピーク周辺に発生するサイドローブを抑圧できる。
ここで、第m番目のチャープパルス送信によって得られる第z番目の信号処理部206におけるビート周波数解析部208から出力されるビート周波数応答をRFTz(fb, m)で表す。ここで、fbはビート周波数インデックスを表し、FFTのインデックス(ビン番号)に対応する。例えば、fb=0,…,Ndata/2であり、z=0,…,Naであり、m=1,…,NCである。なお、以下では、NC回のチャープパルス送信を送信フレーム単位と呼ぶ。ビート周波数インデックスfbが小さいほど、反射波信号の遅延時間が小さい(換言すると、物標との距離が近い)ビート周波数を示す。
また、ビート周波数インデックスf
bは、次式を用いて距離情報R(f
b)に変換できる。そのため、以下では、ビート周波数インデックスf
bを「距離インデックスf
b」と呼ぶ。
ここで、Bwは、チャープ信号におけるレンジゲート内での周波数変調帯域幅を表し、C0は光速度を表す。
ドップラ解析部209は、ビート周波数解析部208から出力される、NC回のチャープパルス送信によって得られるビート周波数応答RFTz(fb, 1)、RFTz(fb, 2)、…、RFTz(fb, NC)を用いて、距離インデックスfb毎にドップラ解析を行う。
例えば、Ncが2のべき乗値である場合、ドップラ解析においてFFT処理を適用できる。この場合、FFTサイズはNcであり、サンプリング定理から導出される折り返しが発生しない最大ドップラ周波数は±1/(2Tr)である。また、ドップラ周波数インデックスfsのドップラ周波数間隔は1/(Nc×Tr)であり、ドップラ周波数インデックスfsの範囲はfs = -Nc/2, …, 0, …, Nc/2-1である。
以下では、一例として、Ncが2のべき乗値である場合について説明する。なお、Ncが2のべき乗でない場合には、例えば、ゼロ埋めしたデータを含めることで2のべき乗個のデータサイズとしてFFT処理が可能である。また、ドップラ解析部209は、FFT処理の際に、Han窓又はHamming窓などの窓関数係数を乗算してもよい。窓関数を適用することでビート周波数ピーク周辺に発生するサイドローブを抑圧できる。
例えば、第z番目の信号処理部206のドップラ解析部209の出力VFT
z(f
b, f
s)は、次式に示す。なお、jは虚数単位であり、z=1~Naである。
以上、信号処理部206の各構成部における処理について説明した。
図1において、CFAR部210は、第1~第Na番目の信号処理部206のドップラ解析部209からの出力を用いて、CFAR処理(換言すると、適応的な閾値判定)を行い、ピーク信号を与える距離インデックスfb_cfar及びドップラ周波数インデックスfs_cfarを抽出する。
CFAR部210は、例えば、次式のように、第1~第Na番目の信号処理部206のドップラ解析部209の出力VFT
1(f
b, f
s)、VFT
2(f
b, f
s)、…、VFT
Na(f
b, f
s)を電力加算し、距離軸とドップラ周波数軸(相対速度に相当)とからなる2次元のCFAR処理、又は、1次元のCFAR処理を組み合わせたCFAR処理を行う。2次元のCFAR処理又は1次元のCFAR処理を組み合わせたCFAR処理については、例えば、非特許文献2に開示された処理が適用されてよい。
CFAR部210は、適応的に閾値を設定し、閾値よりも大きい受信電力となる距離インデックスfb_cfar、ドップラ周波数インデックスfs_cfar、及び、受信電力情報PowerFT(fb_cfar, fs_cfar)をドップラ多重分離部211に出力する。
ドップラ多重分離部211は、CFAR部210から入力される情報(例えば、距離インデックスfb_cfar、ドップラ周波数インデックスfs_cfar、及び、受信電力情報PowerFT(fb_cfar, fs_cfar))に基づいて、各ドップラ解析部209からの出力を用いて、ドップラ多重送信された信号(以下、ドップラ多重信号と呼ぶ)から、各送信アンテナ105から送信される送信信号(換言すると、当該送信信号に対する反射波信号)を分離する。ドップラ多重分離部211は、例えば、分離した信号に関する情報を、方向推定部212に出力する。分離した信号に関する情報には、例えば、分離した信号に対応する距離インデックスfb_cfar、及び、ドップラ周波数インデックス(以下、分離インデックス情報と呼ぶこともある)(fdemul_Tx#1, fdemul_Tx#2, …, fdemul_Tx#Nt)が含まれてよい。また、ドップラ多重分離部211は、各ドップラ解析部209からの出力を方向推定部212に出力する。
以下、ドップラ多重分離部211の動作例について、ドップラシフト部104の動作とともに説明する。
[ドップラシフト量の設定方法]
まず、ドップラシフト部104において付与されるドップラシフト量の設定方法の一例について説明する。
ドップラシフト部104-1~104-Ntは、各々に入力されるチャープ信号に対して異なるドップラシフト量DOPnを付与する。本開示の一実施例では、ドップラシフト部104-1~104-Nt間(換言すると、送信アンテナ105-1~105-Nt間)において、ドップラシフト量DOPnの間隔(ドップラシフト間隔)は、等間隔ではなく、少なくとも一つのドップラ間隔が異なる。
換言すれば、ドップラシフト量DOPnは、サンプリング定理を満たすドップラ周波数範囲(-1/(2Tr) ~1/(2Tr))を等間隔に分割するのではなく、少なくとも一つの間隔が異なるように分割する。ここで、送信周期Tr毎の位相回転が-πからπの範囲であればサンプリング定理を満たす。従って、ドップラシフト量DOPnは、-πからπの範囲、すなわち2πの位相範囲内を等間隔ではなく、少なくとも一つの間隔が異なるように分割する位相回転φn(m)を用いる。
例えば、Nt=2の場合、φ1(m)=π/2×m、φ2(m)=-π/2×mと設定すると、|φ1(m)-φ2(m)|=πとなり、2πの位相範囲を等間隔で分割することになる。本開示の一実施例は、ドップラシフト量として、このような2πの位相範囲内を等間隔に分割するような位相回転を用いない。本開示の一実施例は、ドップラシフト量DOP1、DOP2として|φ1(m)-φ2(m)|≠πとなる位相回転φ1(m)、φ2(m)を用いる。また、Nt≧2の場合、本開示の一実施例は、ドップラシフト量DOPnとして、|φn(m)-φadjacent(n)(m)|≠2π/Ntとなる位相回転が含まれる。ここで、nは1からNtの範囲の整数値である。また、adjacent(n)は、φn(m)に隣接する位相回転のインデックスを表し、φn(m)との位相回転の差分(φn(m)-φn1(m))が、2πのモジュロ演算を用いて最小となるインデックスn1を表す。
例えば、第n番目のドップラシフト部104は、入力された第m番目のチャープ信号に対して、互いに異なるドップラシフト量DOPnとなる位相回転φn(m)を付与して出力する。これにより、複数の送信アンテナ105から送信される送信信号には、それぞれ異なるドップラシフト量が付与される。すなわち、一実施例では、ドップラ多重数NDM=Ntである。ここで、m=1,…,NCであり、n=1,…,Ntである。
また、ドップラ解析部209において、サンプリング定理から導出される、折り返しが発生しないドップラ周波数fdの範囲は-1/(2Tr) ≦ fd <1/(2Tr)である。
このことから、仮に、Nt個の送信アンテナ105から送信される送信信号に対して、ドップラシフト間隔が等間隔1/(Nt×Tr)となる位相回転φ
n(m)は、次式で表される。
ここで、φ0は初期位相であり、Δφ0は基準ドップラシフト位相である。また、round(x)は実数値xに対して四捨五入した整数値を出力するラウンド関数である。なお、round(NC/Nt)の項は、位相回転量を、ドップラ解析部209におけるドップラ周波数間隔の整数倍とする目的で導入されている。
仮に、例えば、式(4)に示す位相回転φn(m)を用いる場合、第m番目のチャープ信号に対して付与される送信信号間の位相回転の間隔は、全て等しくなり、2πround(NC/Nt)/NCとなる。
一例として、式(4)において、Nt=2、Δφ0=0、φ0=0、NCを偶数として位相回転φn(m)が付与される場合、ドップラシフト量は、DOP1=0、DOP2=1/(2Tr)となる。
換言すると、複数の送信アンテナ105から送信される送信信号に対して付与されるドップラシフト量の各間隔は、レーダ装置10(レーダ受信部200)においてドップラ周波数の範囲(例えば、折り返しが発生しないドップラ周波数範囲))において等間隔に設定される。例えば、Nt=2個の送信アンテナ105から送信される送信信号に対して付与されるドップラシフト量の間隔は、折り返しが発生しないドップラ周波数範囲(例えば、-1/(2Tr) ≦ fd <1/(2Tr))を送信アンテナ105の数(例えば、Nt=2)で分割した間隔(上記例では1/(2Tr))に設定される。
図3は、仮に、Nt=2個の送信アンテナ105(以下、Tx#1及びTx#2と呼ぶ)から送信される送信信号に対して、DOP1=0、DOP2=1/(2Tr)のドップラシフト量を用いた場合に、ドップラ解析部209でのドップラ解析により得られるドップラピークの一例を示す。
図3に示すように、測定する1つのターゲットのドップラ周波数(target doppler)fd_TargetDopplerに対して、Nt個(図3ではNt=2)のドップラピークが発生する。
以下、一例として、図3において、測定するターゲットのドップラ周波数fd_TargetDoppler = ‐1/(4Tr)の場合及びfd_TargetDoppler = 1/(4Tr)の場合の送信アンテナTx#1から送信された送信信号に対する反射波信号を受信した際に発生するドップラピークと、送信アンテナTx#2からの送信信号に対する反射波信号を受信した際に発生するドップラピークと、の位置関係を比較する。
<ターゲットのドップラ周波数fd_TargetDoppler=‐1/(4Tr)の場合>
fd_TargetDoppler=‐1/(4Tr)の場合、図3に示すように、送信アンテナTx#1からの送信信号に対する反射波信号を受信した際に発生するドップラピーク(P1)と、送信アンテナTx#2からの送信信号に対する反射波信号を受信した際に発生するドップラピーク(P2)との位置関係となる。ドップラピークP1とドップラピークP2との間のドップラ間隔は1/(2Tr)である。
<ターゲットのドップラ周波数fd_TargetDoppler=1/(4Tr)の場合>
fd_TargetDoppler=1/(4Tr)の場合、図3に示すように、送信アンテナTx#2からの送信信号に対する反射波信号を受信した際に発生するドップラピークは、折り返した信号のピーク(P2A)としてFFT出力される。そのため、fd_TargetDoppler=1/(4Tr)の場合、図3に示すように、送信アンテナTx#1からの送信信号に対する反射波信号を受信した際に発生するドップラピーク(P1)と、上記折り返した信号のドップラピーク(P2A)との位置関係となる。ドップラピーク(P1)とドップラピーク(P2A)との間のドップラ間隔は1/(2Tr)である。
このように、fd_TargetDoppler=‐1/(4Tr)の場合、及び、fd_TargetDoppler=1/(4Tr)の場合において、送信アンテナTx#1に対応するドップラピーク(P1)と、送信アンテナTx#2に対応するドップラピーク(P2又はP2A)との間のドップラ間隔は双方とも1/(2Tr)となる。このため、fd_TargetDoppler=‐1/(4Tr)及び1/(4Tr)では、Tx#1及びTx#2に対応するドップラピークの位置関係の区別がつかなくなり、曖昧性が生じる。従って、図3に示す例では、曖昧性が生じないターゲットのドップラ周波数範囲は、例えば、-1/(4Tr) ≦ fd_TargetDoppler < 1/(4Tr)となる。
これに対して、本開示の一実施例に係るドップラシフト部104では、上述したように、送信アンテナ105から送信される送信信号に対して付与されるドップラシフト量DOPn(又は位相回転φn(m))の間隔は少なくとも1つ異なる。
また、例えば、ドップラシフト部104は、Nt個の送信アンテナ105から送信される送信信号に付与されるドップラシフト量の間隔を可能な限り離しつつ、位相回転φn(m)の少なくとも1つの間隔が異なるようにドップラシフトDOPnを付与する。これにより、ドップラ多重の分離性能を向上できる。
例えば、第n番目のドップラシフト部104は、入力された第m番目のチャープ信号に対して、互いに異なるドップラシフト量DOP
nとなる、次式のような位相回転φ
n(m)を付与する。
ここで、Aは1又は‐1の正負の極性を与える係数である。また、δは1以上の正数である。なお、round(NC/(Nt+δ))の項は、位相回転量を、ドップラ解析部209におけるドップラ周波数間隔の整数倍とする目的で導入されている。
一例として、式(5)において、Nt=2、Δφ0=0、φ0=0、A=1、δ=1、NCを3の倍数として位相回転φn(m)が付与される場合、ドップラシフト量は、DOP1=0、DOP2=1/(3Tr)となる。
図4は、Nt=2個の送信アンテナ105(以下、Tx#1及びTx#2と呼ぶ)から送信される送信信号に対して、DOP1=0、DOP2=1/(3Tr)のドップラシフト量を用いた場合に、ドップラ解析部209でのドップラ解析により得られるドップラピークの一例を示す。
図4に示すように、測定する1つのターゲットのドップラ周波数(target doppler)fd_TargetDopplerに対して、Nt個(図4ではNt=2)のドップラピークが発生する。
以下、一例として、図4において、測定するターゲットのドップラ周波数fd_TargetDoppler = ‐1/(4Tr)の場合及びfd_TargetDoppler = 1/(4Tr)の場合の送信アンテナTx#1から送信された送信信号に対する反射波信号を受信した際に発生するドップラピークと、送信アンテナTx#2から送信された送信信号に対する反射波信号を受信した際に発生するドップラピークと、の位置関係を比較する。
<ターゲットのドップラ周波数fd_TargetDoppler=‐1/(4Tr)の場合>
fd_TargetDoppler=‐1/(4Tr)の場合、図4に示すように、送信アンテナTx#1からの送信信号に対する反射波信号を受信した際に発生するドップラピーク(P1)と、送信アンテナTx#2からの送信信号に対する反射波信号を受信した際に発生するドップラピーク(P2)との位置関係となる。ドップラピークP1とドップラピークP2との間のドップラ間隔は1/(3Tr)である。
<ターゲットのドップラ周波数fd_TargetDoppler = 1/(4Tr)の場合>
fd_TargetDoppler = 1/(4Tr)の場合、送信アンテナTx#2からの送信信号に対する反射波信号を受信した際に発生するドップラピークは、折り返した信号のピーク(P2A)としてFFT出力される。そのため、fd_TargetDoppler = 1/(4Tr)の場合、送信アンテナTx#1からの送信信号に対する反射波信号を受信した際に発生するドップラピーク(P1)と、上記折り返した信号のドップラピーク(P2A)との位置関係となる。ドップラピーク(P1)とピーク(P2A)との間のドップラ間隔は2/(3Tr)である。
図4に示すように、ターゲットのドップラ周波数fd_TargetDoppler = ‐1/(4Tr)の場合、及び、fd_TargetDoppler = 1/(4Tr)の場合において、送信アンテナTx#1に対応するドップラピーク(P1)と、送信アンテナTx#2に対応するドップラピーク(P2又はP2A)との位置関係は互いに異なる。
このように、複数の送信アンテナ105から送信される送信信号に対して付与されるドップラシフト量の各間隔は、ドップラ解析の対象となるドップラ周波数の範囲(例えば、折り返しが発生しないドップラ周波数範囲)において不等間隔に設定される。例えば、Nt=2個の送信アンテナ105から送信される送信信号に対して付与されるドップラシフト量の間隔は、折り返しが発生しないドップラ周波数範囲(例えば、-1/(2Tr) ≦ fd <1/(2Tr))を、送信アンテナ105の数(例えば、Nt=2)にδ=1を加算した値で分割した間隔(上記例では1/(3Tr))に設定される。
これにより、例えば、図4に示すように、折り返しが無い場合(例えば、ドップラピーク(P1)及びドップラピーク(P2))のドップラ間隔(1/(3Tr))と、折り返しが有る場合(例えば、ドップラピーク(P1)及びドップラピーク(P2A))のドップラ間隔(2/(3Tr))とが異なる。
よって、図4に示す例では、ドップラ多重分離部211は、ターゲットのドップラ周波数fd_TargetDoppler=‐1/(4Tr)の場合(換言すると、折り返し無しの場合)、及び、fd_TargetDoppler=1/(4Tr)の場合(換言すると、折り返し有りの場合)とを区別できる。
例えば、想定するターゲットドップラ周波数が-1/(2Tr) ≦ fd_TargetDoppler < 1/(2Tr)の場合、ドップラ多重分離部211は、ターゲットのドップラ周波数fd_TargetDoppler = -1/(4Tr)の場合に、折り返し信号を含まないと判別できる。よって、例えば、図4に示すfd_TargetDoppler = -1/(4Tr)の場合、ドップラ多重分離部211は、折り返し信号を含まず、周波数が最も小さいドップラピークから、それぞれ送信アンテナTx#1、Tx#2からの送信信号に対する反射波信号であると判別できる。
また、例えば、想定するターゲットドップラ周波数が-1/(2Tr) ≦ fd_TargetDoppler < 1/(2Tr)の場合、ドップラ多重分離部211は、ターゲットのドップラ周波数fd_TargetDoppler = 1/(4Tr)の場合に、折り返したドップラピーク(例えば、P2A)が含まれると判別でき、ドップラ周波数fd_TargetDoppler= 1/(4Tr)であると判定できる。よって、例えば、図4に示すfd_TargetDoppler= 1/(4Tr)の場合、折り返し信号(P2A)が含まれるので、ドップラ多重分離部211は、ドップラピークの間隔が2/(3Tr)となるドップラピークのうち高い方のドップラピークが送信アンテナTx#1に対応する反射波信号であり、低い方のドップラピークが送信アンテナTx#2に対応する反射波信号であると判別できる。
次に、他の例として、図4において、測定するターゲットのドップラ周波数fd_TargetDoppler = ‐1/(2Tr)の場合及びfd_TargetDoppler = 1/(2Tr)の場合の送信アンテナTx#1から送信された送信信号に対する反射波信号を受信した際に発生するドップラピークと、送信アンテナTx#2から送信された送信信号に対する反射波信号を受信した際に発生するドップラピークと、の位置関係を比較する。
<ターゲットのドップラ周波数fd_TargetDoppler=‐1/(2Tr)の場合>
fd_TargetDoppler=‐1/(2Tr)の場合、図4に示すように、送信アンテナTx#1からの送信信号に対する反射波信号を受信した際に発生するドップラピーク(P1)と、送信アンテナTx#2からの送信信号に対する反射波信号を受信した際に発生するドップラピーク(P2)との位置関係となる。ドップラピーク(P1)とドップラピーク(P2)との間のドップラ間隔は1/(3Tr)である。
<ターゲットのドップラ周波数fd_TargetDoppler=1/(2Tr)の場合>
fd_TargetDoppler=1/(2Tr)の場合、図4に示すように、送信アンテナTx#2からの送信信号に対する反射波信号を受信した際に発生するドップラピークは、折り返した信号のドップラピーク(P2A)としてFFT出力される。そのため、送信アンテナTx#1からの送信信号に対する反射波信号を受信した際に発生するドップラピーク(P1)と、上記折り返した信号のドップラピーク(P2A)との位置関係となる。ドップラピーク(P1)とドップラピーク(P2A)との間のドップラ間隔は1/(3Tr)である。
このように、ターゲットのドップラ周波数fd_TargetDoppler=‐1/(2Tr)の場合、及び、fd_TargetDoppler=1/(2Tr)の場合において、送信アンテナTx#1に対応するドップラピーク(P1)と、送信アンテナTx#2に対応するドップラピーク(P2又はP2A)との間のドップラ間隔は双方とも1/(3Tr)となる。このため、fd_TargetDoppler=‐1/(2Tr)の場合及びfd_TargetDoppler=1/(2Tr)では、Tx#1及びTx#2に対応するドップラピークの位置関係の区別がつかなくなり、曖昧性が生じる。従って、図4に示す例では、曖昧性が生じないターゲットのドップラ周波数範囲は、例えば、-1/(2Tr) ≦ fd_TargetDoppler < 1/(2Tr)となる。
したがって、本実施の形態によれば、曖昧性が生じないターゲットのドップラ周波数範囲を、時分割多重、又はドップラシフト量を等間隔にする場合のドップラ多重(例えば、図3を参照)と比較して、Nt倍(例えば、図4では2倍)に拡大できる。
次に、ドップラ多重分離部211における各送信アンテナ105に対応する信号の分離方法の一例について説明する。
一例として、Nt=2の場合のドップラ多重分離部211の動作について説明する。
以下では、ドップラシフト部104において、一例として、式(5)に示す位相回転φn(m)が付与される場合について説明する。なお、以下では、一例として、Δφ0=0、φ0=0、δ=1、NCを3の倍数とする。A=1の場合、各送信アンテナ105に対するドップラシフト量はDOP1=0、DOP2=1/(3Tr)であり、A=-1の場合、各送信アンテナ105に対するドップラシフト量はDOP1=0、DOP2=-1/(3Tr)である。
この場合、ドップラ多重分離部211は、CFAR部210から入力される閾値よりも大きい受信電力となるピーク(距離インデックスfb_cfar及びドップラ周波数インデックスfs_cfar)を用いて、ドップラ多重信号を分離する。
例えば、ドップラ多重分離部211は、距離インデックスfb_cfarが同一の複数のドップラ周波数インデックスfs_cfarに対して、送信アンテナTx#1~Tx#Ntから送信される送信信号の何れに対応する反射波信号であるかを判定する。ドップラ多重分離部211は、判定した送信アンテナTx#1~Tx#Nt毎の反射波信号を分離して出力する。
以下では、距離インデックスfb_cfarが同一の複数のドップラ周波数インデックスfs_cfarがNs個ある場合の動作について説明する。例えば、fs_cfar∈{fd#1,fd#2…,fd#Ns}とする。
ドップラ多重分離部211は、例えば、距離インデックスfb_cfarが同一の複数のドップラ周波数インデックスfs_cfar ∈{fd#1,fd#2…,fd#Ns}に対して、ドップラインデックス間隔を算出する。
ここで、送信アンテナTx#1及びTx#2からそれぞれ送信される送信信号に付与されるドップラシフト量DOP
1、DOP
2によって、1つのターゲットドップラ周波数f
d_TargetDopplerに対して、Nt=2個のドップラピークが発生する。このドップラピーク間のドップラ間隔に相当するドップラインデックス間隔は、次式に示す送信アンテナTx#1に対する位相回転φ
1(m)と送信アンテナTx#2に対する位相回転φ
2(m)との差分から、round(N
c/(Nt+1))となる。また、折り返し信号を含む場合、ドップラピーク間のドップラ間隔に相当するドップラインデックス間隔は、N
c-round(N
c/(Nt+1))となる。
そして、ドップラ多重分離部211は、折り返し信号を含まない場合のドップラシフト量の間隔に相当するドップラインデックス間隔round(Nc/(Nt+1))と一致するドップラ周波数インデックス、又は、折り返し信号を含む場合のドップラシフト量の間隔に相当するドップラインデックス間隔(Nc-round(Nc/(Nt+1)))と一致するドップラ周波数インデックスを探索する。
ドップラ多重分離部211は、上述した探索の結果に基づいて、以下の処理を行う。
(1)折り返し信号を含まない場合のドップラシフト量の間隔に相当するインデックス間隔round(Nc/(Nt+1))と一致するドップラ周波数インデックスがある場合、ドップラ多重分離部211は、それらのドップラ周波数インデックスのペア(例えば、fd#p,fd#qと表す)を、ドップラ多重信号の分離インデックス情報(fdemul_Tx#1, fdemul_Tx#2)として出力する。
ここで、送信アンテナTx#1及びTx#2に対するドップラシフト量がDOP1<DOP2の関係の場合、ドップラ多重分離部211は、fd#p,fd#qのうち大きい方をTx#2に対応するドップラ周波数インデックスfdemul_Tx#2と判定し、低い方をTx#1に対応するドップラ周波数インデックスfdemul_Tx#1と判定する。一方、送信アンテナTx#1及びTx#2に対するドップラシフト量がDOP1>DOP2の関係の場合、ドップラ多重分離部211は、fd#p,fd#qのうち大きい方をTx#1に対応するドップラ周波数インデックスfdemul_Tx#1と判定し、低い方をTx#2に対応するドップラ周波数インデックスfdemul_Tx#2と判定する。
(2)折り返し信号を含む場合のドップラシフト量の間隔に相当するインデックス間隔Nc- round(Nc/(Nt+1))と一致するドップラ周波数インデックスがある場合、ドップラ多重分離部211は、それらのドップラ周波数インデックスのペア(例えば、fd#p,fd#q)を、ドップラ多重信号の分離インデックス情報(fdemul_Tx#1, fdemul_Tx#2)として出力する。
ここで、送信アンテナTx#1及びTx#2に対するドップラシフト量がDOP1<DOP2の関係の場合、ドップラ多重分離部211は、fd#p,fd#qのうち大きい方をTx#1に対応するドップラ周波数インデックスfdemul_Tx#1と判定し、低い方をTx#2に対応するドップラ周波数インデックスfdemul_Tx#2と判定する。一方、送信アンテナTx#1及びTx#2に対するドップラシフト量がDOP1>DOP2の関係の場合、ドップラ多重分離部211は、fd#p,fd#qのうち大きい方をTx#2に対応するドップラ周波数インデックスfdemul_Tx#2と判定し、低い方をTx#1に対応するドップラ周波数インデックスfdemul_Tx#1と判定する。
(3)折り返し信号を含まない場合のドップラシフト量の間隔に相当するインデックス間隔round(Nc/(Nt+1))と一致するドップラ周波数インデックス、及び、折り返し信号を含む場合のドップラシフト量の間隔に相当するインデックス間隔Nc-round(Nc/(Nt+1))と一致するドップラ周波数インデックスが無い場合、ドップラ多重分離部211は、発生したドップラピークをノイズ成分と判定する。この場合、ドップラ多重分離部211は、ドップラ多重信号の分離インデックス情報(fdemul_Tx#1, fdemul_Tx#2)を出力しなくてよい。
(4)折り返し信号を含まない場合のドップラシフト量の間隔に相当するインデックス間隔round(Nc/(Nt+1))と一致し、かつ、折り返し信号を含む場合のドップラシフト量の間隔に相当するインデックス間隔Nc-round(Nc/(Nt+1))と一致するドップラ周波数インデックスがある場合、ドップラ多重分離部211は、例えば、以下のような重複除去処理を行う。
例えば、折り返し信号を含まない場合のドップラシフト量の間隔に相当するインデックス間隔round(Nc/(Nt+1))と一致するドップラ周波数インデックスのペアを(fd#p,fd#q1)と表す。また、折り返し信号を含む場合のドップラシフト量の間隔に相当するインデックス間隔Nc-round(Nc/(Nt+1))と一致するドップラ周波数インデックスのペアを(fd#p,fd#q2)と表す。
ドップラ多重分離部211は、例えば、ドップラ周波数インデックスのペア(fd#p,fd#q1)の電力差分|PowerFT(fb_cfar, fd#q1)-PowerFT(fb_cfar, fd#p)|、及び、ドップラ周波数インデックスのペア(fd#p,fd#q2)の電力差分|PowerFT(fb_cfar, fd#q2)-PowerFT(fb_cfar, fd#p)|を算出する。ドップラ多重分離部211は、それらの電力差分間の電力(換言すると差)が所定電力閾値TPLよりも大きい場合、ドップラ周波数インデックスのペア間の電力差分が小さい方のペアを採用する。
例えば、次式を満たす場合、ドップラ多重分離部211は、ドップラ周波数インデックスのペア(fd#p,fd#q2)を採用し、上述した(2)の処理を行う。
|PowerFT(fb_cfar, fd#q1)-PowerFT(fb_cfar, fd#p)|
-|PowerFT(fb_cfar, fd#q2)-PowerFT(fb_cfar, fd#p)|>TPL (7)
また、例えば、次式を満たす場合、ドップラ多重分離部211は、ドップラ周波数インデックスのペア(fd#p,fd#q1)を採用し、上述した(1)の処理を行う。
|PowerFT(fb_cfar, fd#q2)-PowerFT(fb_cfar, fd#p)|
-|PowerFT(fb_cfar, fd#q1)-PowerFT(fb_cfar, fd#p)|>TPL (8)
式(7)及び式(8)を満たさない場合、ドップラ多重分離部211は、何れのドップラ周波数インデックスのペアも採用せずに、上述した(3)の処理を行う。
以上のようにして、ドップラ多重分離部211は、ドップラ多重信号を分離できる。
以上、ドップラ多重分離部211の動作例について説明した。
図1において、方向推定部212は、ドップラ多重分離部211から入力される情報(例えば、距離インデックスfb_cfar、及び、分離インデックス情報(fdemul_Tx#1, fdemul_Tx#2, …, fdemul_Tx#Nt))に基づいて、ターゲットの方向推定処理を行う。
例えば、方向推定部212は、ドップラ多重分離部211の出力から、距離インデックスfb_cfar、及び、分離インデックス情報(fdemul_Tx#1, fdemul_Tx#2, …, fdemul_Tx#Nt)に対応する出力を抽出し、次式に示すような仮想受信アレー相関ベクトルh(fb_cfar, fdemul_Tx#1, fdemul_Tx#2, …, fdemul_Tx#Nt)を生成し、方向推定処理を行う。
仮想受信アレー相関ベクトルh(f
b_cfar, f
demul_Tx#1, f
demul_Tx#2, …, f
demul_Tx#Nt)は、送信アンテナ数Ntと受信アンテナ数Naとの積であるNt×Na個の要素を含む。仮想受信アレー相関ベクトルh(f
b_cfar, f
demul_Tx#1, f
demul_Tx#2, …, f
demul_Tx#Nt)は、ターゲットからの反射波信号に対して各受信アンテナ202間の位相差に基づく方向推定を行う処理に用いる。ここで、z=1,…,Naである。
式(9)において、hcal[b]は、送信アレーアンテナ間及び受信アレーアンテナ間の位相偏差及び振幅偏差を補正するアレー補正値である。b=1,…,(Nt×Na)である。
方向推定部212は、例えば、方向推定評価関数値PH(θ, fb_cfar, fdemul_Tx#1, fdemul_Tx#2, …, fdemul_Tx#Nt)における方位方向θを所定の角度範囲内で可変として空間プロファイルを算出する。方向推定部212は、算出した空間プロファイルの極大ピークを大きい順に所定数抽出し、極大ピークの方位方向を到来方向推定値(例えば、測位出力)として出力する。
なお、方向推定評価関数値PH(θBEAM_cfar, fb_cfar, fdemul_Tx#1, fdemul_Tx#2, …, fdemul_Tx#Nt)は、到来方向推定アルゴリズムによって各種の方法がある。例えば、非特許文献3に開示されているアレーアンテナを用いた推定方法を用いてもよい。
例えば、Nt×Na個の仮想受信アレーが等間隔d
Hで直線状に配置される場合、ビームフォーマ法は次式のように表すことができる。他にも、Capon, MUSICといった手法も同様に適用可能である。
ここで、式(10)において、上付き添え字Hはエルミート転置演算子である。また、a(θu)は、方位方向θuの到来波に対する仮想受信アレーの方向ベクトルを示す。
また、方位方向θuは到来方向推定を行う方位範囲内を所定の方位間隔β1で変化させたベクトルである。例えば、θuは以下のように設定される。
θu=θmin + uβ1、u=0,…, NU
NU=floor[(θmax-θmin)/β1]+1
ここでfloor(x)は、実数xを超えない最大の整数値を返す関数である。
なお、ドップラ周波数情報は相対速度成分に変換して出力されてもよい。ドップラ周波数インデックスf
sを相対速度成分v
d(f
s)に変換するには、次式を用いて変換することができる。ここで、λは送信無線部(図示せず)から出力されるRF信号のキャリア周波数の波長である。また、Δ
fは、ドップラ解析部209におけるFFT処理でのドップラ周波数間隔である。例えば、本実施の形態では、Δ
f=1/(N
CTr)である。
以上のように、本実施の形態では、レーダ装置10は、送信信号を送信する複数の送信アンテナ105と、複数の送信アンテナ105毎に送信信号に異なるドップラシフト量を付与するドップラシフト部104を備える。また、レーダ装置10において、複数の送信アンテナ105から送信される送信信号に対して付与されるドップラシフト量の各間隔は、レーダ装置10においてドップラ周波数の範囲において不等間隔に設定される。
これにより、レーダ装置10では、折り返しが有る場合と無い場合とで、各送信信号に対応するドップラピークの間隔を異ならせることができる。換言すると、レーダ装置10は、ドップラピークの折り返しの有無を判定できる。これにより、レーダ装置10では、折り返しが発生する場合のターゲットドップラ周波数(target doppler)と、折り返しが発生しない場合のターゲットドップラ周波数と、を区別して、ドップラ多重信号を分離できる。よって、レーダ装置10では、ドップラ多重信号を分離可能なドップラ周波数範囲(又は相対速度の最大値)を拡大させることができる。
以上のように、本実施の形態によれば、曖昧性が生じないドップラ周波数範囲(又は相対速度の最大値)を拡大させることができる。これにより、レーダ装置10は、より広いドップラ周波数範囲において、物標(例えば、到来方向)を精度良く検知することができる。
(バリエーション1)
上記実施の形態では、Nt=2の場合のドップラ多重の動作例について説明した。しかし、送信アンテナ数Ntは2個に限らず、3個以上でもよい。
バリエーション1では、他の例として、Nt=3の場合のレーダ装置10の動作について説明する。
以下では、ドップラシフト部104において、一例として、式(5)に示す位相回転φn(m)が付与される場合について説明する。なお、以下では、一例として、Δφ0=0、φ0=0、δ=1、NCを偶数とする。例えば、A=1の場合、各送信アンテナ105に対するドップラシフト量はDOP1=0、DOP2=1/(4Tr)、DOP3=1/(2Tr)である。また、例えば、A=-1の場合、各送信アンテナ105に対するドップラシフト量はDOP1=0、DOP2=-1/(4Tr)、DOP3=-1/(2Tr)である。
このようなドップラシフト量を用いる場合、例えば、図5に示すように、測定する1つのターゲットドップラ周波数fd_TargetDopplerに対して、Nt個(図5では3つ)のドップラピークが発生する。なお、図5は、横軸にターゲットドップラ周波数を示し、縦軸にドップラ解析部209(FFT)の出力を示した場合のNt=3のドップラピークの変化を示した図である。
<ターゲットドップラ周波数が0≦ fd_TargetDoppler <1/(2Tr)の場合>
図5に示すように、送信アンテナTx#1からの送信信号に対する反射波信号を受信した際に発生するドップラピーク(実線)と、送信アンテナTx#3からの送信信号に対する反射波信号を受信した際に発生するドップラピーク(破線)との間のドップラ間隔は1/(2Tr)である。
また、この場合、Tx#3に関して折り返し信号を含む。そのため、ドップラ多重分離部211は、ドップラピークの間隔が1/(2Tr)となるドップラピークのうち高い方のドップラピークが送信アンテナTx#1に対応する反射波信号であり、低い方のドップラピークが送信アンテナTx#3に対応する反射波信号であり、残りのドップラピークが送信アンテナTx#2からの反射波信号であると判別できる。
<ターゲットドップラ周波数が-1/(2Tr)≦ fd_TargetDoppler <0の場合>
図5に示すように、送信アンテナTx#1からの送信信号に対する反射波信号を受信した際に発生するドップラピーク(実線)と、送信アンテナTx#2からの送信信号に対する反射波信号を受信した際に発生するドップラピーク(点線)との間のドップラ間隔は1/(4Tr)である。また、送信アンテナTx#2からの送信信号に対する反射波信号を受信した際に発生するドップラピーク(点線)と、送信アンテナTx#3からの送信信号に対する反射波信号を受信した際に発生するドップラピーク(破線)との間のドップラ間隔は1/(4Tr)である。
また、この場合、何れの送信アンテナTx#1, Tx#2, Tx#3についても折り返し信号を含まない。そのため、ドップラ多重分離部211は、周波数が低いドップラピークから、それぞれ、送信アンテナTx#1、Tx#2、Tx#3からの送信信号に対する反射波信号であると判別できる。
このように、複数の送信アンテナ105から送信される送信信号に対して付与されるドップラシフト量の各間隔は、ドップラ周波数範囲(図5の例では、例えば、-1/(2Tr) ≦ fd <1/(2Tr))において不等間隔に設定される。例えば、Nt=3個の送信アンテナから送信される送信信号に対して付与されるドップラシフト量の各間隔は、折り返しが発生しないドップラ周波数範囲(例えば、-1/(2Tr) ≦ fd <1/(2Tr))を、送信アンテナの数(例えば、Nt=3)にδ=1を加算した値で分割した間隔(上記例では1/(4Tr))に設定される。
これにより、例えば、図5に示すように、折り返しが無い場合のドップラ間隔(1/(4Tr))と、折り返しが有る場合のドップラ間隔(1/(4Tr)及び1/(2Tr))とが異なる。
よって、図5に示す例では、ドップラ多重分離部211は、ターゲットのドップラ周波数が-1/(2Tr)≦ fd_TargetDoppler <0の場合(換言すると、折り返し無しの場合)、及び、ターゲットのドップラ周波数が0≦ fd_TargetDoppler <1/(2Tr)の場合(換言すると、折り返し有りの場合)とを区別できる。
したがって、図5に示す例では、曖昧性が生じないターゲットのドップラ周波数範囲は、例えば、-1/(2Tr) ≦ fd_TargetDoppler<1/(2Tr)となる。
よって、曖昧性が生じないターゲットのドップラ周波数範囲を、時分割多重、又はドップラシフト量を等間隔にする場合(図5では1/(3Tr)の場合)のドップラ多重と比較して、Nt倍(例えば、図5では3倍)に拡大できる。
次に、ドップラ多重分離部211における各送信アンテナ105に対応する信号の分離方法の一例について説明する。
ドップラ多重分離部211は、CFAR部210から入力される閾値よりも大きい受信電力となるピーク(距離インデックスfb_cfar及びドップラ周波数インデックスfs_cfar)を用いて、ドップラ多重信号の分離を行う。
例えば、ドップラ多重分離部211は、距離インデックスfb_cfarが同一の複数のドップラ周波数インデックスfs_cfarに対して、送信アンテナTx#1~Tx#Ntから送信される送信信号の何れに対応する反射波信号であるかを判定する。ドップラ多重分離部211は、判定した送信アンテナTx#1~Tx#Nt毎の反射波信号を分離して出力する。
ドップラ多重分離部211は、例えば、距離インデックスfb_cfarが同一の複数のドップラ周波数インデックスfs_cfar ∈{fd#1,fd#2…,fd#Ns}に対して、ドップラインデックス間隔を算出する。
そして、ドップラ多重分離部211は、3つのドップラ周波数インデックスを小さい順に見た場合の2つのドップラインデックス間隔が、折り返し信号を含まない場合のドップラシフト量の間隔に相当するインデックス間隔round(Nc/(Nt+1))及びround(Nc/(Nt+1))と一致するドップラ周波数インデックスの組み合わせを探索する。又は、ドップラ多重分離部211は、3つのドップラ周波数インデックスを小さい順に見た場合の2つのドップラインデックス間隔が、折り返し信号を含む場合のドップラシフト量の間隔に相当するインデックス間隔round(Nc/(Nt+1))及びNc-round(Nc/(Nt+1))、又は、Nc-round(Nc/(Nt+1))及びround(Nc/(Nt+1))と一致するドップラ周波数インデックスの組み合わせを探索する。
ドップラ多重分離部211は、上述した探索の結果に基づいて、以下の処理を行う。
(1)折り返し信号を含まない場合のドップラシフト量の間隔に相当するインデックス間隔round(Nc/(Nt+1))及びround(Nc/(Nt+1))と一致するドップラ周波数インデックスの組み合わせがある場合、ドップラ多重分離部211は、それらのドップラ周波数インデックスの組(例えば、fd#p1,fd#p2,fd#p3と表す)を、ドップラ多重信号の分離インデックス情報(fdemul_Tx#1, fdemul_Tx#2, fdemul_Tx#3)として出力する。
ここで、送信アンテナTx#1~Tx#3に対するドップラシフト量がDOP1<DOP2<DOP3の関係の場合、ドップラ多重分離部211は、fd#p1,fd#p2,fd#p3のうち大きい方から、Tx#3、Tx#2、Tx#1にそれぞれ対応するドップラ周波数インデックスfdemul_Tx#3, fdemul_Tx#2, fdemul_Tx#1と判定する。また、送信アンテナTx#1~Tx#3に対するドップラシフト量がDOP1>DOP2>DOP3の関係の場合、ドップラ多重分離部211は、fd#p1,fd#p2,fd#p3のうち大きい方から、Tx#1、Tx#2、Tx#3にそれぞれ対応するドップラ周波数インデックスfdemul_Tx#1, fdemul_Tx#2, fdemul_Tx#3と判定する。
(2)折り返し信号を含む場合のドップラシフト量の間隔に相当するインデックス間隔Nc-round(Nc/(Nt+1))及びround(Nc/(Nt+1))と一致するドップラ周波数インデックスの組み合わせがある場合、ドップラ多重分離部211は、それらのドップラ周波数インデックスの組(例えば、fd#q1,fd#q2,fd#q3と表す)を、ドップラ多重信号の分離インデックス情報(fdemul_Tx#1, fdemul_Tx#2, fdemul_Tx#3)として出力する。
ここで、送信アンテナTx#1~Tx#3に対するドップラシフト量がDOP1<DOP2<DOP3の関係の場合、ドップラ多重分離部211は、fd#q1,fd#q2,fd#q3のうち大きい方からTx#2、Tx#1、Tx#3にそれぞれ対応するドップラ周波数インデックスfdemul_Tx#2, fdemul_Tx#1, fdemul_Tx#3と判定する。また、送信アンテナTx#1~Tx#3に対するドップラシフト量がDOP1>DOP2>DOP3の関係の場合、ドップラ多重分離部211は、fd#q1,fd#q2,fd#q3のうち大きい方からTx#2、Tx#3、Tx#1にそれぞれ対応するfdemul_Tx#2, fdemul_Tx#3, fdemul_Tx#1と判定する。
(3)折り返し信号を含む場合のドップラシフト量の間隔に相当するインデックス間隔round(Nc/(Nt+1))及びNc-round(Nc/(Nt+1))と一致するドップラ周波数インデックスの組み合わせがある場合、ドップラ多重分離部211は、それらのドップラ周波数インデックスの組(例えば、fd#u1,fd#u2,fd#u3と表す)を、ドップラ多重信号の分離インデックス情報(fdemul_Tx#1, fdemul_Tx#2, fdemul_Tx#3)として出力する。
ここで、送信アンテナTx#1~Tx#3に対するドップラシフト量がDOP1<DOP2<DOP3の関係の場合、ドップラ多重分離部211は、fd#u1,fd#u2,fd#u3のうち大きい方からTx#1、Tx#3、Tx#2にそれぞれ対応するfdemul_Tx#3, fdemul_Tx#2, fdemul_Tx#1と判定する。また、送信アンテナTx#1~Tx#3に対するドップラシフト量がDOP1>DOP2>DOP3の関係の場合、ドップラ多重分離部211は、fd#u1,fd#u2,fd#u3のうち大きい方からTx#3、Tx#1、Tx#2にそれぞれ対応するfdemul_Tx#3, fdemul_Tx#1, fdemul_Tx#2と判定する。
(4)ドップラ多重分離部211は、上記の(1)、(2)及び(3)の何れにも該当しないドップラ周波数インデックスに対応するドップラピークをノイズ成分と判定する。この場合、ドップラ多重分離部211は、ドップラ多重信号の分離インデックス情報(fdemul_Tx#1, fdemul_Tx#2, fdemul_Tx#3)を出力しなくてよい。
(5)ドップラ多重分離部211は、上記の(1)、(2)及び(3)に重複して該当するドップラ周波数インデックスを含む場合、例えば、以下のような重複除去処理を行う。
例えば、上記の(1)及び(2)に重複して該当するドップラ周波数インデックスを含むドップラ周波数インデックスの組がそれぞれ(fd#p1,fd#p2,fd#p3)、(fd#q1,fd#q2,fd#q3)である場合、ドップラ多重分離部211は、それぞれの組の中でのドップラ周波数インデックスの受信電力{例えば、PowerFT(fb_cfar, fd#p1)、PowerFT(fb_cfar, fd#p2)、PowerFT(fb_cfar, fd#p3)}、{PowerFT(fb_cfar, fd#q1)、PowerFT(fb_cfar, fd#q2)、PowerFT(fb_cfar, fd#q3)}を比較して、最小となる受信電力をそれぞれの組から抽出する。そして、ドップラ多重分離部211は、例えば、各組の最小電力となる電力間の電力差が所定電力閾値TPLよりも大きいドップラ周波数インデックスの組を採用する。
例えば、次式を満たす場合、ドップラ多重分離部211は、ドップラ周波数インデックス(fd#p1,fd#p2,fd#p3)の組を採用し、上述した(1)の処理を行う。
Min({PowerFT(fb_cfar, fd#p1), PowerFT(fb_cfar, fd#p2), PowerFT(fb_cfar, fd#p3) })
- Min({PowerFT(fb_cfar, fd#q1), PowerFT(fb_cfar, fd#q2), PowerFT(fb_cfar, fd#q3) })
> TPL (13)
また、例えば、次式を満たす場合、ドップラ多重分離部211は、ドップラ周波数インデックス(fd#q1,fd#q2,fd#q3)の組を採用し、上述した(2)の処理を行う。
Min({PowerFT(fb_cfar, fd#q1), PowerFT(fb_cfar, fd#q2), PowerFT(fb_cfar, fd#q3) })
- Min({PowerFT(fb_cfar, fd#p1), PowerFT(fb_cfar, fd#p2), PowerFT(fb_cfar, fd#p3) })
> TPL (14)
ドップラ多重分離部211は、式(13)又は式(14)を満たさない場合、何れのドップラ周波数インデックスの組も採用せずに、上述した(4)の処理を行う。また、ドップラ多重分離部211は、(1)及び(2)以外の重複の組あわせに対しても同様な重複判定処理を行う。
以上のようにして、ドップラ多重分離部211は、ドップラ多重信号を分離できる。
(バリエーション2)
上記実施の形態では、送信信号に付与されるドップラシフト量に対応する位相回転の一例として、式(5)に示す位相回転φn(m)を用いる場合について説明した。しかし、位相回転は、式(5)に示す位相回転φn(m)に限定されない。
他の例として、第n番目のドップラシフト部104は、入力された第m番目のチャープ信号(送信信号)に対して、式(5)を用いる場合と異なるドップラシフト量DOP
nとなる次式の位相回転φ
n(m)を付与してよい。
ここで、dpnは位相回転をドップラ周波数範囲において不等間隔とする成分である。例えば、dp1、dp2、…、dpNtは、-round(NC/Nt)/2< dpn < round(NC/Nt)/2の範囲内の値であり、全てが同一の値ではなく、少なくとも一つは異なる値の成分を含む。なお、round(NC/Nt)の項は、位相回転量を、ドップラ解析部209におけるドップラ周波数間隔の整数倍とする目的で導入されている。
一例として、式(15)において、Nt=2、Δφ0=0、φ0=0、A=1、dp1=0、dp2=π/5、NCを偶数とした場合の位相回転φn(m)が付与される場合、ドップラシフト量は、DOP1=0、DOP2=1/(2Tr)+1/(10Tr)=6/(10Tr)となる。
図6は、横軸にターゲットドップラ周波数を示し、縦軸にドップラ解析部209(FFT)の出力を示した場合のNt=2、DOP1=0、DOP2=6/(10Tr)のドップラピークの変化を示した図である。
<ターゲットドップラ周波数が-1/(10Tr)≦ fd_TargetDoppler <1/(2Tr)の場合>
図6に示すように、送信アンテナTx#1からの送信信号に対する反射波信号を受信した際に発生するドップラピーク(実線)と、送信アンテナTx#2からの送信信号に対する反射波信号を受信した際に発生するドップラピーク(点線)との間のドップラ間隔は4/(10Tr)である。
また、この場合、Tx#2に関して折り返し信号を含む。そのため、ドップラ多重分離部211は、ドップラピークの間隔が4/(10Tr)となるドップラピークのうち高い方のドップラピークが送信アンテナTx#1に対応する反射波信号であり、低い方のドップラピークが送信アンテナTx#2に対応する反射波信号であると判別できる。
<ターゲットドップラ周波数が-1/(2Tr)≦ fd_TargetDoppler < -1/(10Tr)の場合>
図6に示すように、送信アンテナTx#1からの送信信号に対する反射波信号を受信した際に発生するドップラピーク(実線)と、送信アンテナTx#2からの送信信号に対する反射波信号を受信した際に発生するドップラピーク(点線)との間のドップラ間隔は6/(10Tr)である。
また、この場合、何れの送信アンテナTx#1、Tx#2についても折り返し信号を含まない。そのため、ドップラ多重分離部211は、例えば、周波数が低いドップラピークから、それぞれ送信アンテナTx#1、Tx#2からの送信信号に対する反射波信号であると判別できる。
このように、バリエーション2では、送信アンテナ105のドップラシフト量の各間隔は、ドップラ周波数範囲(図6では、例えば、-1/(2Tr) ≦ fd <1/(2Tr))を複数の送信アンテナ105の数(例えば、Nt=2)で分割した間隔にDOP2=6/(10Tr)のオフセットを加えた間隔に設定される。
これにより、例えば、図6に示すように、折り返しが無い場合のドップラ間隔(6/(10Tr))と、折り返しが有る場合のドップラ間隔(4/(10Tr))とが異なる。
したがって、図6に示す例では、曖昧性が生じないターゲットのドップラ周波数範囲は、例えば、-1/(2Tr) ≦ fd_TargetDoppler<1/(2Tr)となる。
よって、バリエーション2によれば、曖昧性が生じないターゲットのドップラ周波数範囲を、時分割多重又はドップラ多重と比較して、Nt倍(例えば、図6では2倍)に拡大できる。
(バリエーション3)
ドップラ多重において、複数のターゲットのドップラピークの受信レベルがほぼ等しく、ドップラピークの間隔がドップラシフト量の間隔に一致する場合、ドップラ多重分離部211において分離判定ができなくなる可能性がある。
しかしながら、複数のターゲット間においてドップラ周波数が異なる場合、ターゲットとレーダ装置10との間の相対移動速度は異なる。そのため、レーダ装置10においてレーダ観測を続けて行うことにより、或るレーダ装置の測位出力において複数のターゲットのドップラピークの受信レベルがほぼ等しく、ドップラピークの間隔がドップラシフト量の間隔に一致する場合でも、続くレーダ装置の測位出力では、複数のターゲット間の距離が異なって測定される可能性が高い。よって、続くレーダ装置の測位出力において、複数のターゲットが分離された出力を得ることができると考えられる。
バリエーション3では、レーダ装置10の測位出力において、複数のターゲットをより確実に分離するために、例えば、レーダ観測毎にドップラシフト量が可変に設定される場合について説明する。なお、レーダ観測の単位は、例えば、送信フレーム単位でもよく、他の単位でもよい。
例えば、バリエーション3では、ドップラシフト量DOPnに対応する位相回転φn(m)として式(5)を用いてよい。
レーダ装置10は、レーダ観測毎に式(5)におけるδの値を可変に設定することにより、各送信アンテナ105に対するドップラシフト量の間隔を可変に設定できる。例えば、レーダ観測毎にδは1、2、1、2と周期的に可変されてもよい。
また、ドップラシフト量DOPnに対応する位相回転φn(m)として式(15)を用いてよい。例えば、レーダ装置10は、位相回転を不等間隔とする成分dp1、dp2、…、dpNtをレーダ観測毎に異なる値に設定することにより、各送信アンテナ105に対するドップラシフト量の間隔を可変に設定できる。
バリエーション3によれば、1つのターゲットに対する複数の送信アンテナ105に対応するドップラピークの間隔がレーダ観測毎に異なるので、複数のターゲットを分離しやすくなる。
(バリエーション4)
バリエーション4では、レーダ装置の送信アンテナがサブアレー構成である場合について説明する。
送信アンテナのうちのいくつかを組み合わせてサブアレーとして用いることにより、送信指向性ビームパターンのビーム幅を狭めて、送信指向性利得を向上できる。これにより、検知可能な角度範囲は狭まるが、検知可能な距離範囲を増加できる。また、指向性ビームを生成するビームウェイト係数を可変にすることにより、ビーム方向を可変制御できる。
図7は、バリエーション4に係るレーダ送信部100aの構成例を示すブロック図である。なお、図7において、図1に示すレーダ送信部100と同様の動作を行う構成には同一の符号を付し、その説明を省略する。
また、バリエーション4に係るレーダ受信部は、図1に示すレーダ受信部200と基本構成が共通するので、図1を援用して説明する。
図7において、NDMはドップラ多重数を示す。
図7では、各ドップラシフト部104の出力に対して、NSA個の送信アンテナ105を用いたサブアレーが構成される。したがって、送信アンテナ105の数Nt=NSA×NDMとなる。なお、送信アンテナ105のサブアレー構成は、図7に示す例に限定されない。例えば、各ドップラシフト部104の出力に対するサブアレーに含まれる送信アンテナ数は、ドップラシフト104間で同数でなくてもよい。ここで、NSAは1以上の整数である。なお、NSA=1の場合は、図1と同様になる。なお、ドップラシフト部104は、例えば、サブアレー構成の送信アンテナ105(例えば、NSA個の送信アンテナ105)から送信されるレーダ送信信号に、同一のドップラシフト量を付与する。
図7において、ビームウェイト生成部106は、サブアレーを用いて送信ビームの主ビーム方向を所定方向に向けるビームウェイトを生成する。例えば、N
SA個の送信アンテナを用いたサブアレーが、素子間隔d
SAで直線配置される場合の送信ビーム方向をθ
TxBFと表す。この場合、ビームウェイト生成部106は、例えば、次式のようなビームウェイトW
Tx(Index_TxSubArray, θ
TxBF)を生成する。
ここで、Index_TxSubArrayは、サブアレーの要素インデックスを示し、Index_TxSubArray =1,…, NSAである。また、λはレーダ送信信号の波長を示し、dSAはサブアレーアンテナ間隔を示す。
例えば、第ndm番目のビームウェイト乗算部107は、第ndm番目のドップラシフト部104からの出力に対して、ビームウェイト生成部106から入力されるビームウェイト係数WTx(Index_TxSubArray, θTxBF)を乗算する。ビームウェイトWTx(Index_TxSubArray, θTxBF)が乗算された送信信号は、第{NSA×(ndm-1)+ Index_TxSubArray}番目の送信アンテナ105から送信される。ここで、Index_TxSubArray =1,…, NSAであり、ndm=1,…,NDMである。
以上の動作により、レーダ送信部100aは、ドップラシフト部104からの出力に対して、サブアレーを用いて所定方向に送信指向性ビームを向ける送信が可能となる。これにより、所定方向の送信指向性利得を向上でき、検知可能な距離範囲を拡大できる。
また、レーダ送信部100aは、送信指向性ビームを生成するビームウェイト係数を可変に設定することにより、ビーム方向を可変制御できる。
なお、バリエーション4において説明したサブアレー送信を行う構成は、他のバリエーション又は実施の形態においても同様に適用できる。
(バリエーション5)
バリエーション5では、例えば、同一周波数帯又は一部の周波数帯が重複する複数のレーダ装置からの干渉の影響を低減する方法について説明する。
図8は、バリエーション5に係るレーダ装置10bの構成例を示すブロック図である。なお、図8において、図1と同一の構成には同一の符号を付し、その説明を省略する。例えば、図8に示すレーダ装置10bは、図1に示すレーダ装置10に対して、レーダ送信部100bにおいてランダム符号生成部108及びランダム符号乗算部109を追加し、レーダ受信部200bにおいてランダム符号乗算部213を追加した構成である。
図8において、ランダム符号生成部108は、例えば、擬似ランダムな符号系列RCode={RC(1),RC(2),…,RC(NLRC)}を生成する。例えば、擬似ランダムな符号には、PN(pseudo random noise)符号、M系列符号、又は、Gold符号を用いてよい。また、ランダム符号生成部108は、擬似ランダムな符号系列の符号要素{1、-1}に対し、例えば、{π、-π}の位相回転を与える信号を生成する。
擬似ランダムな符号系列の符号長NLRCは、Nc以下である。また、ランダム符号生成部108は、送信周期m毎に擬似ランダムな符号系列の符号要素インデックスをRC_INDEX(m) =mのように可変して、擬似ランダムな符号系列RCodeのランダム符号要素RC(RC_INDEX(m))をランダム符号乗算部109及び213に出力する。
レーダ送信部100bのランダム符号乗算部109は、送信周期mのチャープ信号cp(t)に対して、ランダム符号生成部108から入力されるランダム符号要素RC(RC_INDEX)を乗算する。ランダム符号乗算部109は、RC(RC_INDEX(m))×cp(t)で表される信号を各ドップラシフト部104に出力する。
レーダ受信部200bのランダム符号乗算部213は、送信周期mにおけるビート周波数解析部208の出力信号RFTz(fb, m)に対して、ランダム符号生成部108から入力されるランダム符号要素RC(RC_INDEX)を乗算する。ランダム符号乗算部213は、RC(RC_INDEX (m))×RFTz(fb、m)で表される信号をドップラ解析部209に出力する。ここで、z=1、…、Naである。
以上のような動作により、同一周波数帯又は一部の周波数帯が重複する複数のレーダ装置からの干渉が存在する場合でも、レーダ装置10bにおいて、干渉信号は、ドップラ解析部209に入力される前にランダム符号乗算部213によって擬似ランダムな信号に変換できる。これにより、ドップラ解析部209の出力において、干渉波の信号電力をドップラ周波数領域に拡散する効果が得られる。例えば、擬似ランダムな符号系列の乗算によって、干渉波のピーク電力を1/Nc程度に低減できる。したがって、後段のCFAR部210において、干渉波のピークが誤って検出さる確率を大幅に低減できる。
(バリエーション6)
例えば、ドップラシフト量DOPnとして式(5)に示す位相回転を用いる場合、ドップラ周波数範囲をドップラ多重数NDMより多い数(NDM+δ)で等分割した間隔(ΔFD=round(NC/(NDM+δ))に対して、ドップラシフト量の間隔には、ΔFDの間隔、及び、(δ+1)ΔFDの間隔が用いられる。
そのため、ドップラ多重された各信号は、ドップラ解析部209(例えば、図1を参照)の出力では、ドップラ周波数領域においてΔFDの間隔で折り返したように検出される。
このような性質を用いると、例えば、CFAR部210及びドップラ多重分離部211の動作を以下のように簡易化できる。
[CFAR部210の動作]
CFAR部210は、例えば、CFAR処理の対象のドップラ周波数範囲のうち、レーダ送信信号にそれぞれ付与されるドップラシフト量の各間隔に対応する範囲(例えば、ΔFD)毎の反射波信号の受信電力を加算した電力加算値に対して閾値を用いて、ドップラピークを検出する。
例えば、CFAR部210は、第1~第Na番目の信号処理部206のドップラ解析部209からの出力に対して、次式に示すように、ΔFDの範囲で折り返して加算した電力加算値を算出して、CFAR処理を行う。ここで、f
s_shrink=-N
c,…,-N
c+ΔFD-1である。
これにより、CFAR処理の対象のドップラ周波数範囲を1/(NDM+δ)にすることができ、CFAR処理の演算量削減を図ることができる。
CFAR部210は、適応的に閾値を設定し、閾値よりも大きい受信電力となる距離インデックスfb_cfar、ドップラ周波数インデックスfshrink_cfar、及び、受信電力情報(PowerFT(fb_cfar, fshrink_cfar+ndm×ΔFD)、ndm=1,…,NDM)をドップラ多重分離部211に出力する。
[ドップラ多重分離部211の動作]
ドップラ多重分離部211は、CFAR部210から入力される受信電力情報(PowerFT(fb_cfar, fshrink_cfar+ndm×ΔFD)、ndm=1,…,NDM)を比較して、受信電力の上位NDM個のドップラ周波数インデックスの受信レベルと、上位NDM以外を除いたδ個のドップラ周波数インデックスの受信レベルとの差が大きく異なる場合(例えば、所定閾値以上の場合)、受信レベルが小さいδ個のドップラ周波数インデックスを(δ+1)ΔFDの間隔に含まれるものと判定し、受信電力の上位NDM個のドップラ周波数インデックスをドップラ多重信号の分離インデックス情報(fdemul_Tx#1,…, fdemul_Tx#NDM)として出力する。
換言すると、ドップラ多重分離部211は、ドップラ周波数範囲において検出されたドップラピークのうち、受信電力の高い順にNDM個のドップラピークに対応する受信レベルと、上記NDM個のドップラピーク以外のドップラピーク(例えば、δ個のドップラピーク)に対応する受信レベルとの差が閾値以上の場合、上記NDM個のドップラピークに基づいて、反射波信号からドップラ多重信号をそれぞれ分離する。なお、受信レベルの差は、例えば、NDM個の受信レベル及びδ個の受信レベルのそれぞれの平均値の差でもよい。または、受信レベルの差は、NDM個の受信レベルのうちの最小値と、δ個の受信レベルのうちの最大値との差としてもよい。
なお、上記処理に限らず、例えば、送信アンテナ105と、送信アンテナ105からそれぞれ送信されるレーダ送信信号に付与されるドップラシフト量との関係に基づいて、反射波信号から、ドップラ多重信号をそれぞれ分離してもよい。例えば、ドップラ多重信号の分離インデックス情報は、(δ+1)ΔFDの間隔となるドップラ周波数インデックス情報と受信電力の上位NDM個のドップラ周波数インデックスの相対的な位置関係を用いて決定されてよい。例えば、図5では、NDM=3、δ=1を用いて式(5)に示す位相回転を用いたドップラシフト量が付与されている。従って、ターゲットドップラ周波数はΔFDのドップラ間隔と(δ+1)ΔFDのドップラ間隔とが含まれる。(δ+1)ΔFDのドップラ間隔となるドップラ周波数インデックスは、図5の場合、fdemul_Tx#1とfdemul_Tx#3であることは既知であり、ドップラ多重分離部211は、これを利用して、ドップラ多重信号の分離インデックス情報を決定することができる。すなわち、ドップラ解析部209の出力において、(δ+1)ΔFDのドップラ間隔が0から1/(2T)の範囲にある場合は、(δ+1)ΔFDのドップラ間隔となるドップラ周波数インデックスの高い方がfdemul_Tx#1であり、低い方がfdemul_Tx#3である。また、(δ+1)ΔFDのドップラ間隔が-1/(2T)から0の範囲にある場合は、fdemul_Tx#3のドップラ周波数インデックスが折り返して発生することを考慮して、(δ+1)ΔFDのドップラ間隔となるドップラ周波数インデックスの高い方がfdemul_Tx#3であり、低い方がfdemul_Tx#1となる。受信電力の上位NDM個のドップラ周波数インデックスのうち、残りのドップラ周波数インデックスがfdemul_Tx#2となる。ドップラ多重分離部211は、この結果を用いることで、ドップラシフト量DOPnを判定することができ、ドップラ多重信号を分離できる。
このように、ドップラ多重分離部211では、受信電力情報PowerFT(fb_cfar、fshrink_cfar+ndm×ΔFD)、ndm=1,…,NDMの比較処理によってドップラ多重分離が可能となるため、ドップラ分離処理を削減できる。
(実施の形態2)
本実施の形態では、ドップラ多重送信と、符号分割多重(CDM:Code Division Multiplexing)送信とを併用する場合について説明する。
例えば、実施の形態1(例えば、図1を参照)においてドップラ多重数が多くなると、ドップラ多重分離部211の処理において、折り返しがある場合のドップラシフト量の間隔と、折り返しがない場合のドップラシフト量の間隔とが重複するドップラ周波数インデックスが存在する確率が増加する。したがって、反射物体が多い伝搬環境に依存して、ドップラ多重数には適した範囲があり、上限となるドップラ多重数が存在する。
そこで、本実施の形態では、実施の形態1において説明したドップラ多重を行う構成に、さらに符号多重を併用することにより、送信アンテナ数(例えば、ドップラ多重数)が増加した場合でも、ドップラ領域と符号領域とを用いて多重数を増大可能な構成について説明する。
図9は、本実施の形態に係るレーダ装置10cの構成例を示すブロック図である。なお、図9において、実施の形態1(例えば、図1)と同様の構成には同一の符号を付し、その説明を省略する。例えば、図9に示すレーダ装置10cには、図1に示すレーダ装置10に対して、レーダ送信部100cにおいて直交符号生成部301及び直交符号乗算部302が追加され、レーダ受信部200cにおいて出力切替部401及び符号多重分離部402が追加されている。
以下では、ドップラ多重数をNDMとし、符号多重数をNCMとし、送信アンテナ105の数Nt=NDM×NCMとなるドップラ多重数と符号多重数を用いる場合について説明する。
[レーダ送信部100cの構成例]
レーダ送信部100cにおいて、直交符号生成部301は、直交符号長LocのNCM個の直交符号系列Codencm={OCncm(1), OCncm(2),…, OCncm(Loc)}をそれぞれ生成する。ここで、ncm=1,…, NCMである。
例えば、直交符号生成部301は、レーダ送信周期(Tr)毎に、直交符号系列Code1~CodeNcmの要素を指示する直交符号要素インデックスOC_INDEXを巡回的に可変設定することにより、直交符号系列Code1~CodeNcmの要素OC1(OC_INDEX)~OCNcm(OC_INDEX)を第1~第Ntの直交符号乗算部302に出力する。また、直交符号生成部301は、レーダ送信周期(Tr)毎に、直交符号要素インデックスOC_INDEXを出力切替部401に出力する。
ここで、OC_INDEX=1, 2, …, Locである。例えば、m番目の送信周期において、OC_INDEX =MOD(m-1, Loc)+1である。ここで、MOD(x, y)はモジュロ演算子であり、xをyで割った後の余りを出力する関数である。
また、直交符号生成部301において生成される直交符号系列には、例えば、互いに無相関となる符号が用いられる。例えば、直交符号系列には、Walsh-Hadamard-符号が用いられてよい。
一例として、NCM=2の場合、Walsh-Hadamard-符号の直交符号長Loc=2であり、直交符号生成部301は、OC1={1,1}、OC2={1,-1}となる直交符号系列を生成する。
また、他の例として、NCM =4の場合、直交符号長Loc =4であり、直交符号生成部301は、OC1={1,1, 1, 1}、OC2={1,-1, 1, -1}, OC3={1,1, -1, -1}、OC4={1,-1, -1, 1}となる直交符号系列を生成する。
なお、直交符号系列を構成する要素は実数に限らず、複素数が含まれてもよく、次式のような位相回転を用いた直交符号でもよい。
式(19)において、例えば、Nt=3の場合、直交符号長Loc =Ntであり、直交符号生成部301は、OC1={1,1,1}、OC2={1, exp(j2π/3),exp(j4π/3)}, OC3={1, exp(-j2π/3),exp(-j4π/3)}となる直交符号系列を生成する。
また、他の例として、Nt=4の場合、直交符号長Loc =Ntであり、直交符号生成部301は、OC1={1,1,1, 1}、OC2={1, j,-1 ,-j},OC3={1, -1, 1, -1}, OC4={1, -j,-1, j}となる直交符号系列を生成する。
図9に示すレーダ送信部100cは、例えば、ドップラ多重数をNDMとした場合、NDM個のドップラシフト部104-1~104-NDMを備える。また、レーダ送信部100cは、ドップラシフト部104と同数のNDM個の直交符号乗算部302を備える。
各ドップラシフト部104は、レーダ送信信号生成部101から入力されるチャープ信号に対して、所定のドップラシフト量DOPndmを付与するために、所定の位相回転φndmを付与して、位相回転を付与したチャープ信号を、対応する直交符号乗算部302に出力する。ここで、ndm=1,…, NDMである。
各直交符号乗算部302は、符号多重数NCMに相当する数の乗算器を備える。直交符号乗算部302は、ドップラシフト部104の出力に対して、NCM個の直交符号系列Code1、Code2、…、CodeNcmをそれぞれ乗算し、NCM個の信号を送信アンテナ105に出力する。
以上のようなドップラシフト部104及び直交符号乗算部302の動作により、Nt個の送信アンテナ105のうち、第n番目の送信アンテナ105からは、レーダ送信信号生成部101の出力に対して、第floor[(n-1)/NCM]+1番目のドップラシフト部104によるドップラシフトDOPfloor[(n-1)/ NCM]+1が付与され、さらに、第floor[(n-1)/ NCM]+1番目の直交符号乗算部302における第mod(n-1, NCM)+1番目の直交符号Codemod(n-1, NCM)+1が乗算された信号が出力される。
例えば、Nt=6の送信アンテナ105、ドップラ多重数をNDM=3、符号多重数をNCM=2とした場合について説明する。この場合、3個(=NDM)のドップラシフト部104は、チャープ信号に対して、それぞれドップラシフト量DOP1、DOP2、DOP3を付与する。また、3個(=NDM)の直交符号乗算部302は、ドップラシフト部104の出力に対して、2(=NCM)個の直交符号系列Code1、及び、Code2を乗算する。
この場合、例えば、第1番目の送信アンテナ105からは、送信周期Tr毎に以下の信号が出力される。
ここで、cp(t)は送信周期Tr毎のチャープ信号を表す。また、ドップラシフト部104における位相回転φ
ndm(m)を付与する際の乗算値を次式に示すΛ
ndm(m)と表記している。
同様に、第2番目の送信アンテナ105からは、送信周期Tr毎に以下の信号が出力される。
同様に、第3番目の送信アンテナ105からは、送信周期Tr毎に以下の信号が出力される。
同様に、第4番目の送信アンテナ105からは、送信周期Tr毎に以下の信号が出力される。
同様に、第5番目の送信アンテナ105からは、送信周期Tr毎に以下の信号が出力される。
同様に、第6番目の送信アンテナ105からは、送信周期Tr毎に以下の信号が出力される。
また、レーダ送信部100cは、チャープパルス送信回数が直交符号長Locの整数倍(Ncode倍)となるように信号送信する。例えば、NC=LOC×Ncodeとする。
なお、レーダ装置10cにおけるレーダ送信部の構成は、図9に示す構成に限定されない。例えば、図10に示すレーダ送信部100dのように、図9に示すドップラシフト部104の位相回転と直交符号乗算部302における符号乗算とを同時に行う構成でもよい。なお、図10に示すレーダ受信部200dは、図9に示すレーダ受信部200cの構成と同様である。
例えば、図10において、ドップラシフト及び直交符号生成部303は、送信周期Tr毎のドップラシフト及び直交符号を行う乗算係数を生成する。例えば、ドップラシフト及び直交符号生成部303は、Nt個の送信アンテナ105のうち、第n番目の送信アンテナに接続されている乗算部304に対して、第floor[(n-1)/ NCM]+1番目のドップラシフトDOPfloor[(n-1)/ NCM]+1を付与するための位相回転と、第mod(n-1, NCM)+1番目の直交符号Codemod(n-1, NCM)+1を乗算した乗算係数を出力する。
乗算部304は、レーダ送信信号生成部101の出力信号(チャープ信号)に対して、ドップラシフ及び直交符号生成部303から入力される乗算係数を乗算する。
[レーダ受信部200cの構成例]
次に、図9に示すレーダ受信部200cの構成例について説明する。
第z番目の信号処理部206cにおいて、出力切替部401は、直交符号生成部301から入力される直交符号要素インデックスOC_INDEXに基づいて、送信周期Tr毎のビート周波数解析部208の出力を、Loc個のドップラ解析部209-1~209-LOCのうち、OC_INDEX番目のドップラ解析部209に選択的に切り替えて出力する。すなわち、出力切替部401は、第m番目の送信周期Trにおいて、OC_INDEX番目のドップラ解析部209を選択する。
第z番目の信号処理部206cはLoc個のドップラ解析部209を有する。
第z番目の信号処理部206cにおける第nol番目のドップラ解析部209には、出力切替部401により、Loc回の送信周期毎(LOC×Tr)にデータが入力される。そのため、第nol番目のドップラ解析部209は、NC回の送信周期のうち、Ncode回の送信周期のデータを用いてドップラ解析を行う。ここで、nol=1,…, LOCである。
また、ドップラ解析部209は、Ncodeが2のべき乗値である場合、次式に示すようなFFT(高速フーリエ変換)処理を適用することができる。
ここで、FFTサイズはNcodeであり、サンプリング定理から導出される折り返しが発生しない最大ドップラ周波数は、±1/(2Loc×Tr)である。また、ドップラ周波数インデックスfsのドップラ周波数間隔は1/(Ncode× Loc×Tr)であり、ドップラ周波数インデックスfsの範囲はfs=-Ncode/2,…,0,…, Ncode/2-1である。
なお、Ncodeが2のべき乗でない場合、例えば、ゼロ埋めしたデータを含めることで2のべき乗個のFFTサイズとしてFFT処理が可能である。また、ドップラ解析部209は、FFT処理の際に、Han窓又はHamming窓などの窓関数係数を乗算してもよく、窓関数を適用することでビート周波数ピーク周辺に発生するサイドローブを抑圧できる。
符号多重分離部402は、直交符号を用いて多重送信された信号を分離する。
例えば、符号多重分離部402は、次式のように送信時に用いた直交符号要素OC
ncmを複素共役(*で表す)して、直交符号要素インデックスOC_INDEX毎のドップラ解析結果に乗算して加算する。これにより、直交符号Code
ncmを用いて符号多重送信された信号を分離した信号が得られる。ここで、ncm=1,…,N
CMである。
CFAR部210cは、符号多重分離部402の出力を用いて、CFAR処理(換言すると、適応的な閾値判定)を行い、ピーク信号を与える距離インデックスfb_cfar及びドップラ周波数インデックスfs_cfarを抽出する。
CFAR部210cは、例えば、次式のように、符号多重分離部402の出力を電力加算し、距離軸及びドップラ周波数軸(相対速度に相当)とからなる2次元のCFAR処理、又は、1次元のCFAR処理を組み合わせたCFAR処理を行う。2次元のCFAR処理又は1次元のCFAR処理を組み合わせたCFAR処理については、例えば、非特許文献2に開示されている処理が適用されてよい。
CFAR部210cは、適応的に閾値を設定し、閾値よりも大きい受信電力となる距離インデックスfb_cfar、ドップラ周波数インデックスfs_cfar、及び、受信電力情報PowerFT(fb_cfar, fs_cfar)をドップラ多重分離部211cに出力する。
なお、図9においては、CFAR部210cは、符号多重分離部402の出力を用いた構成を示したが、これに限定されない。別な構成として、CFAR部210cは、ドップラ解析部209の出力を用いてCFAR処理を行ってもよい。この場合、CFAR部210cは、例えば、次式のようにドップラ解析部209の出力を電力加算し、距離軸及びドップラ周波数軸(相対速度に相当)とからなる2次元のCFAR処理、又は、1次元のCFAR処理を組み合わせたCFAR処理を行ってよい。2次元のCFAR処理又は1次元のCFAR処理を組み合わせたCFAR処理については、例えば、非特許文献2に開示されている処理が適用されてよい。
また、CFAR部210cがドップラ解析部209の出力を用いてCFAR処理を行う場合、符号多重分離部402は、CFAR部210cで指示された閾値よりも大きい受信電力となる距離インデックスfb_cfar、ドップラ周波数インデックスfs_cfar、及び、受信電力情報PowerFT(fb_cfar, fs_cfar)を用いて符号多重分離の動作を行ってもよく、これにより、CFAR部210cで指示された閾値よりも大きい受信電力となる距離インデックスfb_cfar、ドップラ周波数インデックスfs_cfarに限定して、符号多重分離の動作を行うことができるため、符号多重分離部402の処理量を削減することができる。
ドップラ多重分離部211cは、CFAR部210cから入力される情報(例えば、距離インデックスfb_cfar、ドップラ周波数インデックスfs_cfar、及び、受信電力情報PowerFT(fb_cfar, fs_cfar))に基づいて、各符号多重分離部402からの出力を用いて、各送信アンテナ105から送信される送信信号を分離する。
以下、ドップラ多重分離部211cの動作について、ドップラシフト部104の動作とともに説明する。
第1~第NDM番目のドップラシフト部104は、入力されたチャープ信号に対して、異なるドップラシフト量DOP1, DOP2, …, DOPNDMをそれぞれ付与する。ここで、ドップラシフト量DOP1, DOP2, …, DOPNDMの各間隔(ドップラシフト間隔)は、実施の形態1と同様、例えば、折り返しが発生しないドップラ周波数範囲を等間隔に分割した値ではなく、不等間隔に分割した間隔である(例えば、少なくとも一つのドップラ間隔が異なる)。例えば、ドップラシフト量DOPndmの各間隔は、ドップラ周波数範囲(例えば、-1/(2Loc×Tr) ≦ fd <1/(2Loc×Tr))を複数の送信アンテナ105の数Ntを符号多重数NCMで除算した値(換言すると、ドップラ多重数NDM)に、1以上(例えば、δ)を加算した整数値で分割した間隔に設定されてよい。
なお、実施の形態1では、ドップラ多重数を送信アンテナ数Ntに等しい場合(すなわち、Nt=NDM)について説明した。これに対して、本実施の形態では、ドップラ多重に対して符号多重を併用するので、ドップラ多重数NDMは、送信アンテナ数Ntより少ない多重数となる(例えば、Nt>NDM)。よって、例えば、ドップラシフト量DOPndmの各間隔は、折り返しが発生しないドップラ周波数範囲(例えば、-1/(2Loc×Tr) ≦ fd <1/(2Loc×Tr))を、送信アンテナ105の数Nt以下で分割した間隔に設定されてよい。
したがって、本実施の形態では、ドップラシフト量DOPndmに対して、実施の形態1で用いた式(5)又は式(15)の中のNtをNDMで置き換えた式を用いる。また、直交符号系列を乗算する送信周期(LOC×Tr)では同一の位相回転とするために、直交符号長Locの送信周期間(LOC×Tr)では同じ位相回転φndm(m)が繰り返して出力される。
すなわち、第ndm番目のドップラシフト部104は、入力された第m番目のチャープ信号に対して、異なるドップラシフト量DOP
ndmとなる、次式のような位相回転φ
ndm(m)を付与する。
ここで、Aは1又は‐1の正負の極性を与える係数である。また、δは1以上の正数である。また、φ0は初期位相であり、Δφ0は基準ドップラシフト位相である。なお、round(x)は実数値xに対し、四捨五入した整数値を出力するラウンド関数である。また、floor[x]は、実数x以下で最も近い整数を出力する演算子である。また、round(Ncode/(NDM+1))の項は、位相回転量を、ドップラ解析部209におけるドップラ周波数間隔の整数倍とする目的で導入している。
上述したように、実施の形態1ではドップラ多重数NDMが送信アンテナ数Ntに等しい場合について説明したのに対して、本実施の形態では、ドップラ多重数NDMは、送信アンテナ数Ntより少ない。ドップラ多重分離部211cでは、実施の形態1に係るドップラ多重分離部211(例えば、図1を参照)において用いたパラメータNtをNDMで置き換える。
また、実施の形態1では、ドップラ解析部209(例えば、図1を参照)におけるFFTサイズがNCであるのに対して、本実施の形態では、Ncodeである。そこで、ドップラ多重分離部211cでは、実施の形態1に係るドップラ多重分離部211において用いたパラメータNCをNcodeで置き換える。
また、実施の形態1では、ドップラ解析部209におけるFFTのサンプリング周期がTrであるのに対して、本実施の形態では、LOC×Trである。そこで、ドップラ多重分離部211cでは、実施の形態1に係るドップラ多重分離部211において用いたパラメータTrをLOC×Trで置き換える。
一例として、NDM=2、Δφ0=0、φ0=0、δ=1、Ncodeを3の倍数として位相回転φndm(m)(例えば、式(31))が付与される場合、ドップラシフト量は、A=1の場合、DOP1=0、DOP2=1/(3LOC×Tr)となり、A=-1の場合、DOP1=0、DOP2=-1/(3LOC×Tr)となる。
ドップラ多重分離部211cは、CFAR部210cから入力される閾値よりも大きい受信電力となるピーク(距離インデックスfb_cfar、ドップラ周波数インデックスfs_cfar)を用いて、ドップラ多重信号を分離する。
例えば、ドップラ多重分離部211cは、距離インデックスfb_cfarが同一の複数のドップラ周波数インデックスfs_cfarに対して、ドップラ多重送信信号#1~#NDMの何れの反射波信号であるかを判定する。ドップラ多重分離部211cは、判定したドップラ多重送信信号毎の反射波信号を分離して出力する。
以下では、距離インデックスfb_cfarが同一の複数のドップラ周波数インデックスfs_cfarがNs個ある場合の動作について説明する。例えば、fs_cfar ∈{fd#1,fd#2…,fd#Ns}とする。
ドップラ多重分離部211cは、例えば、距離インデックスfb_cfarが同一の複数のドップラ周波数インデックスfs_cfar ∈{fd#1,fd#2…,fd#Ns}に対して、ドップラインデックス間隔を算出する。
ここで、ドップラシフト量DOP
ndmによって、1つのターゲットドップラ周波数f
d_TargetDopplerに対して、ドップラ解析部でのドップラ解析により得られるドップラスペクトラムには、N
DM個(ただし、N
DM<Nt)のドップラピークが発生する。このドップラピーク間のドップラ間隔に相当するドップラインデックス間隔は、次式に示す、位相回転φ
1(m)と位相回転φ
2(m)との差分から、round(Ncode/(N
DM+1))となる。また、折り返し信号を含む場合、ドップラピーク間のドップラ間隔に相当するドップラインデックス間隔は、N
c-round(Ncode/(N
DM+1))となる。
そして、ドップラ多重分離部211cは、折り返し信号を含まない場合のドップラシフト量の間隔に相当するインデックス間隔round(Ncode/(NDM+1))と一致するドップラ周波数インデックス、又は、折り返し信号を含む場合のドップラシフト量の間隔に相当するインデックス間隔Nc-round(Ncode/(NDM+1))と一致するドップラ周波数インデックスを探索する。
ドップラ多重分離部211cは、上述した探索の結果に基づいて、以下の処理を行う。
(1)折り返し信号を含まない場合のドップラシフト量の間隔に相当するインデックス間隔round(Ncode/(NDM+1))と一致するドップラ周波数インデックスがある場合、ドップラ多重分離部211cは、それらのドップラ周波数インデックスのペア(例えば、fd#p,fd#qと表す)を、ドップラ多重信号の分離インデックス情報(fdemul_DS#1, fdemul_ DS #2)として出力する。
ここで、ドップラシフト量がDOP1<DOP2の関係の場合、ドップラ多重分離部211cは、fd#p,fd#qのうち大きい方を第2番目のドップラシフト部104の出力(DS#2)と判定し、低い方を第1番目のドップラシフト部104の出力(DS#1)と判定する。一方、ドップラシフト量がDOP1>DOP2の関係の場合、ドップラ多重分離部211cは、fd#p,fd#qのうち大きい方を第1番目のドップラシフト部104の出力(DS#1)と判定し、低い方を第2番目のドップラシフト部104の出力(DS#2)と判定する。
(2)折り返し信号を含む場合のドップラシフト量の間隔に相当するインデックス間隔Nc-round(Ncode/(NDM+1))と一致するドップラ周波数インデックスがある場合、ドップラ多重分離部211cは、それらのドップラ周波数インデックスのペア(例えば、fd#p,fd#q)を、ドップラ多重信号の分離インデックス情報(fdemul_DS#1, fdemul_DS#2)として出力する。
ここで、ドップラシフト量がDOP1<DOP2の関係の場合、ドップラ多重分離部211cは、fd#p,fd#qのうち大きい方を第1番目のドップラシフト部104の出力(DS#1)と判定し、低い方を第2番目のドップラシフト部104の出力(DS#2)と判定する。一方、ドップラシフト量がDOP1>DOP2の関係の場合、ドップラ多重分離部211cは、fd#p,fd#qのうち大きい方を第2番目のドップラシフト部104の出力(DS#2)と判定し、低い方を第1番目のドップラシフト部104の出力(DS#1)と判定する。
(3)折り返し信号を含まない場合のドップラシフト量の間隔に相当するインデックス間隔round(Ncode/(NDM+1))と一致するドップラ周波数インデックス、及び、折り返し信号を含む場合のドップラシフト量の間隔に相当するインデックス間隔Nc-round(Ncode/(NDM+1))と一致するドップラ周波数インデックスが無い場合、ドップラ多重分離部211cは、発生したドップラピークをノイズ成分と判定する。この場合、ドップラ多重分離部211cは、ドップラ多重信号の分離インデックス情報(fdemul_Tx#1, fdemul_Tx#2)を出力しなくてよい。
(4)折り返し信号を含まない場合のドップラシフト量の間隔に相当するインデックス間隔round(Ncode/(NDM+1))と一致し、かつ、折り返し信号を含む場合のドップラシフト量の間隔に相当するインデックス間隔Nc-round(Ncode/(NDM+1))と一致するドップラ周波数インデックスがある場合、ドップラ多重分離部211cは、例えば、以下のような重複除去処理を行う。
例えば、折り返し信号を含まない場合のドップラシフト量の間隔に相当するインデックス間隔round(Ncode/(NDM+1))と一致するドップラ周波数インデックスのペアを(fd#p,fd#q1)と表す。また、折り返し信号を含む場合のドップラシフト量の間隔に相当するインデックス間隔Nc-round(Ncode/(NDM+1))と一致するドップラ周波数インデックスのペアを(fd#p,fd#q2)と表す。
ドップラ多重分離部211cは、例えば、ドップラ周波数インデックスのペア(fd#p,fd#q1)の電力差分|PowerFT(fb_cfar, fd#q1)-PowerFT(fb_cfar, fd#p)|、及び、ドップラ周波数インデックスのペア(fd#p,fd#q2)の電力差分|PowerFT(fb_cfar, fd#q2)-PowerFT(fb_cfar, fd#p)|を算出する。ドップラ多重分離部211cは、それらの電力差分間の電力(換言すると差)が所定電力閾値TPLよりも大きい場合、ドップラ周波数インデックスのペア間の電力差分が小さい方のペアを採用する。
例えば、次式を満たす場合、ドップラ多重分離部211cは、ドップラ周波数インデックスのペア(fd#p,fd#q2)を採用し、上述した(2)の処理を行う。
|PowerFT(fb_cfar, fd#q1)-PowerFT(fb_cfar, fd#p)|
-|PowerFT(fb_cfar, fd#q2)-PowerFT(fb_cfar, fd#p)|>TPL (33)
また、例えば、次式を満たす場合、ドップラ多重分離部211cは、ドップラ周波数インデックスのペア(fd#p,fd#q1)を採用し、上述した(1)の処理を行う。
|PowerFT(fb_cfar, fd#q2)-PowerFT(fb_cfar, fd#p)|
-|PowerFT(fb_cfar, fd#q1)-PowerFT(fb_cfar, fd#p)|>TPL (34)
式(33)及び式(34)を満たさない場合、ドップラ多重分離部211cは、何れのドップラ周波数インデックスのペアも採用せずに、上述した(3)の処理を行う。
以上のようにして、ドップラ多重分離部211cは、ドップラ多重信号を分離できる。
なお、本実施の形態において、式(31)の代わりに、次式の位相回転φ
ndm(m)を用いてもよい。
ここで、dpndmは位相回転をドップラ周波数範囲において不等間隔とする成分である。例えば、dp1、dp2、…、dpDMは、-round(NC/NDM)/2<dpn<round(NC/NDM)/2の範囲の値であり、全てが同一の値ではなく、少なくとも一つは異なる値の成分を含む。なお、round(NC/NDM)の項は、位相回転量を、ドップラ解析部209におけるドップラ周波数間隔の整数倍とする目的で導入されている。
以上、ドップラ多重分離部211cの動作例について説明した。
図9において、方向推定部212cは、ドップラ多重分離部211cから入力される情報(例えば、距離インデックスfb_cfar及び分離インデックス情報(fdemul_DS#1, fdemul_DS#2, …,fdemul_DS#NDM)に基づいて、ターゲットの方向推定処理を行う。
例えば、方向推定部212cは、各符号多重分離部402の出力から、距離インデックスfb_cfar、及び、分離インデックス情報(fdemul_DS#1, fdemul_DS#2, …,fdemul_DS#NDM)に対応する出力を抽出し、次式に示すような仮想受信アレー相関ベクトルh(fb_cfar, fdemul_DSx#1, fdemul_DSx#2, …, fdemul_DS#NDM)を生成し、方向推定処理を行う。
仮想受信アレー相関ベクトルh(f
b_cfar, f
demul_DSx#1, f
demul_DSx#2, …, f
demul_DS#NDM)は、送信アンテナ数Ntと受信アンテナ数Naとの積であるNt×Na個の要素を含む。仮想受信アレー相関ベクトルh(f
b_cfar, f
demul_DSx#1, f
demul_DSx#2, …, f
demul_DS#NDM)は、ターゲットからの反射波信号に対して各受信アンテナ202間の位相差に基づく方向推定を行う処理に用いる。ここで、z=1,…,Naである。なお、方向推定方法は、例えば、実施の形態1と同様の方法が適用されてよい。
式(36)において、hcal[b]は、送信アレーアンテナ間及び受信アレーアンテナ間の位相偏差及び振幅偏差を補正するアレー補正値である。b=1,…, (Nt×Na)である。
以上のように、本実施の形態では、ドップラ多重と符号多重とを併用する構成により、実施の形態1と同様の効果に加え、同時に多重送信する信号数を増大することができ、送信アンテナ数が増大したMIMOアレー構成への適応が可能になる。
なお、以上の説明において、ドップラ多重数をNDMとし、符号多重数をNCMとし、送信アンテナ105の数Nt=NDM×NCMとなるドップラ多重数と符号多重数を用いたが、これに限定されない。例えば、NDM個のドップラ多重信号に対し、符号多重数を同一にせずに異なる符号多重数を用いてもよい。例えば、直交符号生成部301で、直交符号長LocのNCM個の直交符号系列Codencmを生成し、各直交符号乗算部302は、符号多重数NCM以下の乗算器を備えてよい。直交符号乗算部302は、ドップラシフト部104の出力に対して、NCM個の直交符号系列Code1、Code2、…、CodeNcmのうちNCM以下の直交符号系列をそれぞれ乗算し、NCM個以下の信号を送信アンテナ105に出力する構成を用いてもよい。
例えば、Nt=5の送信アンテナ105、ドップラ多重数をNDM=3、符号多重数NCMを2以下とした場合について説明する。この場合、3個(=NDM)のドップラシフト部104は、チャープ信号に対して、それぞれドップラシフト量DOP1、DOP2、DOP3を付与する。また、3個(=NDM)の直交符号乗算部302は、ドップラシフト部104-1およびドップラシフト部104-2の出力に対して、2(=NCM)個の直交符号系列Code1、及び、Code2を乗算し、ドップラシフト部104-3の出力に対して、1(≦NCM)個の直交符号系列Code1を乗算する構成を用いる。換言すると、複数の送信アンテナ105から送信されるレーダ送信信号間で、レーダ送信信号に適用される符号多重数NCMが異なる。この場合のレーダ受信部200cの処理は、ドップラシフト部104-3の出力に対して直交符号系列Code2を乗算して送信した信号の符号多重分離が不要となることを除き、これまで説明した処理(Nt=6の送信アンテナ105、ドップラ多重数をNDM=3、符号多重数をすべてNCM=2とした場合の処理)と同様な処理で、Nt=5の送信アンテナからの送信信号を分離することができる。このように、NDM個のドップラ多重信号に対し、符号多重数を同一にせずに異なる符号多重数を用いることで、ドップラ多重数NDMを超えた送信アンテナ数(換言すれば同時多重送信数)の適用範囲を増やすことができる。例えば、ドップラ多重数をNDM=3、符号多重数NCMを2以下とした場合、送信アンテナ数Nt(換言すれば同時多重送信数)は4、5、6の範囲で用いることが可能である。より一般的に記載すれば、送信アンテナ数Nt(換言すれば同時多重送信数)は、NDM+1≦Nt≦NDM×NCMの範囲で適用が可能である。
また、直交符号乗算部302は、複数のドップラシフト部104の出力のうち、少なくとも1つのドップラシフト部104の出力に対して、NCM個の直交符号系列Code1、Code2、…、CodeNcmのうち1つの直交符号系列を乗算し、送信アンテナ105に出力する構成を用いてもよい。このように複数のドップラシフト部104の出力のうち、少なくとも1つのドップラシフト部104の出力に対して、符号多重を用いない信号を送信アンテナから出力する構成を用いることで、レーダ受信部200cにおいて、ドップラ解析部209の出力にドップラ折り返し信号が含まれるか否かを検出できる。すなわち、ドップラ解析部209がサンプリング定理から導出される折り返しが発生しない最大ドップラ周波数は、±1/(2Loc×Tr)であるが、このように複数のドップラシフト部104の出力のうち、少なくとも1つのドップラシフト部104の出力に対して、符号多重を用いない信号を送信アンテナから出力する構成を用いることで、ドップラ解析部209がサンプリング定理から導出される折り返しが発生しない最大ドップラ周波数を、±1/(2×Tr)とすることができ、曖昧性なく検出できるドップラ周波数範囲を拡大する効果も得ることができる。
なお、ドップラ多重と符号多重とを併用する場合に、実施の形態1のバリエーション5と同様、擬似ランダムな符号系列が送信信号に乗算されてもよい。擬似ランダムな符号系列の符号長NLRCは、Ncode以下として、符号多重周期毎にランダム符号要素インデックスをRC_INDEX(m)=floor[(m-1)/NLOC]+1のように可変して、擬似ランダムな符号系列RCodeのランダム符号要素RC(RC_INDEX(m))を出力してもよい。
(実施の形態3)
本実施の形態では、ドップラ多重送信と、時分割多重(TDM:Time Division Multiplexing)送信とを併用する場合について説明する。
例えば、実施の形態1(例えば、図1を参照)においてドップラ多重数が多くなると、ドップラ多重分離部211の処理において、折り返しがある場合のドップラシフト量の間隔と、折り返しがない場合のドップラシフト量の間隔とが重複するドップラ周波数インデックスが存在する確率が増加する。したがって、反射物体が多い伝搬環境に依存して、ドップラ多重数には適した範囲があり、上限となるドップラ多重数が存在する。
そこで、本実施の形態では、実施の形態1において説明したドップラ多重を行う構成に、さらに時分割多重を併用することにより、送信アンテナ数(例えば、ドップラ多重数)が増加した場合でも、ドップラ領域と時間領域とを用いて多重数を増大可能な構成について説明する。
図11は、本実施の形態に係るレーダ装置10eの構成例を示すブロック図である。なお、図11において、実施の形態1(例えば、図1)と同様の構成には同一の符号を付し、その説明を省略する。例えば、図11に示すレーダ装置10eには、図1に示すレーダ装置10に対して、レーダ送信部100eにおいて送信切替制御部501及び送信切替部502が追加され、レーダ受信部200eにおいて出力切替部601が追加されている。
以下では、ドップラ多重数をNDMとし、時分割多重数をNTMとし、送信アンテナ105の数Nt=NDM×NTMとなるドップラ多重数と時分割多重数を用いる場合について説明する。
[レーダ送信部100eの構成例]
送信切替制御部501は、レーダ送信周期(Tr)毎に、時間多重で用いる、送信アンテナ105の切り替えを指示する時分割多重インデックスTM_INDEXを生成し、時分割多重インデックスTM_INDEXを送信切替部502及び出力切替部601に出力する。
ここで、TM_INDEX=1、2、…、NTMである。例えば、m番目の送信周期において、TM_INDEX =MOD(m-1, NTM)+1である。ここで、MOD(x、y)はモジュロ演算子であり、xをyで割った後の余りを出力する関数である。
図11に示すレーダ送信部100eは、例えば、ドップラ多重数をNDMとした場合、NDM個のドップラシフト部104-1~104-NDMを備える。また、レーダ送信部100eは、ドップラシフト部104と同数のNDM個の送信切替部502を備える。
各ドップラシフト部104は、レーダ送信信号生成部101から入力されるチャープ信号に対して、所定のドップラシフト量DOPndmを付与するために、所定の位相回転φndmを付与して、位相回転を付与したチャープ信号を、対応する送信切替部502に出力する。ここで、ndm=1,…, NDMである。
第ndm番目の送信切替部502は、時分割多重インデックスTM_INDEXの指示に従って、第ndm番目のドップラシフト部104の出力を、第{(ndm-1)×NTM+TM_INDEX}番目の送信アンテナ105に切り替えて出力する。
以上のようなドップラシフト部104及び送信切替部502の動作により、Nt個の送信アンテナ105のうち、第n番目の送信アンテナ105からは、レーダ送信信号生成部101の出力に対して、第floor[(n-1)/NTM]+1番目のドップラシフト部104によるドップラシフトDOPfloor[(n-1)/ NTM]+1が付与された信号が、第floor[(n-1)/ NTM]+1番目の送信切替部502により時分割多重インデックスTM_INDEXがmod(n-1, NTM)+1となるときに出力される。
例えば、Nt=6の送信アンテナ105、ドップラ多重数をNDM=3、時分割多重数をNTM=2とした場合について説明する。この場合、3(=NDM)個のドップラシフト部104は、チャープ信号に対して、それぞれドップラシフト量DOP1、DOP2、DOP3を付与する。また、3(=NDM)個の送信切替部502のそれぞれの時分割多重インデックスTM_INDEXは2(=NTM)個の要素からなる。
この場合、例えば、第1番目の送信アンテナ105からは、送信周期Tr毎に以下の信号が出力される。
ここでcp(t)は送信周期Tr毎のチャープ信号を表す。また、ドップラシフト部104における位相回転φ
ndm(m)を付与する際の乗算値を次式に示すΛ
ndm(m)と表記し、送信信号がない場合はゼロとしている。
同様に、例えば、第2番目の送信アンテナ105からは、送信周期Tr毎に以下の信号が出力される。
同様に、例えば、第3番目の送信アンテナ105からは、送信周期Tr毎に以下の信号が出力される。
同様に、例えば、第4番目の送信アンテナ105からは、送信周期Tr毎に以下の信号が出力される。
同様に、例えば、第5番目の送信アンテナ105からは、送信周期Tr毎に以下の信号が出力される。
同様に、例えば、第6番目の送信アンテナ105からは、送信周期Tr毎に以下の信号が出力される。
また、レーダ送信部100eは、チャープパルス送信回数がNTMの整数倍(Ncode倍)となるように信号送信する。例えば、NC=NTM×Ncodeとする。
[レーダ受信部200eの構成例]
次に、図11に示すレーダ受信部200eの構成例について説明する。
第z番目の信号処理部206eにおいて、出力切替部601は、送信切替制御部501から入力される時分割多重インデックスTM_INDEXに基づいて、送信周期Tr毎のビート周波数解析部208の出力を、NTM個のドップラ解析部209-1~209-NTMのうち、TM_INDEX番目のドップラ解析部209に選択的に切り替えて出力する。すなわち、出力切替部601は、第m番目の送信周期Trにおいて、TM_INDEX番目のドップラ解析部209を選択する。
第z番目の信号処理部206eはNTM個のドップラ解析部209を有する。
第z番目の信号処理部206eにおける第ntm番目のドップラ解析部209には、出力切替部601により、NTM回の送信周期毎(NTM×Tr)にデータが入力される。そのため、第ntm番目のドップラ解析部209は、NC回の送信周期のうち、Ncode回の送信周期のデータを用いてドップラ解析を行う。ここで、ntm=1,…, NTMである。
また、ドップラ解析部209は、Ncodeが2のべき乗値である場合、次式に示すようなFFT(高速フーリエ変換)処理を適用することができる。
ここで、FFTサイズはNcodeであり、サンプリング定理から導出される折り返しが発生しない最大ドップラ周波数は、±1/(2 NTM×Tr)である。また、ドップラ周波数インデックスfsのドップラ周波数間隔は1/(Ncode× NTM×Tr)であり、ドップラ周波数インデックスfsの範囲はfs=-Ncode/2,..,0,..., Ncode/2-1である。
なお、Ncodeが2のべき乗でない場合、例えば、ゼロ埋めしたデータを含めることで2のべき乗個のFFTサイズとしてFFT処理が可能である。また、FFT処理の際に、Han窓又はHamming窓などの窓関数係数を乗算してもよく、窓関数を適用することでビート周波数ピーク周辺に発生するサイドローブを抑圧できる。
CFAR部210eは、全ての信号処理部206eの第1~第NTM番目のドップラ解析部209の出力を用いて、CFAR処理(換言すると、適応的な閾値判定)を行い、ピーク信号を与える距離インデックスfb_cfar及びドップラ周波数インデックスfs_cfarを抽出する。
CFAR部210eは、例えば、次式のように、ドップラ解析部209の出力を電力加算し、距離軸及びドップラ周波数軸(相対速度に相当)とからなる2次元のCFAR処理、又は、1次元のCFAR処理を組み合わせたCFAR処理を行う。2次元のCFAR処理又は1次元のCFAR処理を組み合わせたCFAR処理については、例えば、非特許文献2に開示されている処理が適用されてよい。
CFAR部210eは、適応的に閾値を設定し、閾値よりも大きい受信電力となる距離インデックスfb_cfar、ドップラ周波数インデックスfs_cfar、及び、受信電力情報をPowerFT(fb_cfar, fs_cfar)をドップラ多重分離部211eに出力する。
ドップラ多重分離部211eは、CFAR部210eから入力される情報(例えば、距離インデックスfb_cfar、ドップラ周波数インデックスfs_cfar、及び、受信電力情報PowerFT(fb_cfar, fs_cfar))に基づいて、各ドップラ解析部209からの出力を用いて、各送信アンテナ105から送信される送信信号を分離する。
以下、ドップラ多重分離部211eの動作について、ドップラシフト部104の動作とともに説明する。
第1~第NDMのドップラシフト部104は、入力されたチャープ信号に対して、異なるドップラシフト量DOP1,DOP2,…,DOPNDMをそれぞれ付与する。ここで、ドップラシフト量DOP1,DOP2,…,DOPNDMの各間隔(ドップラシフト間隔)は、実施の形態1と同様、例えば、折り返しが発生しないドップラ周波数範囲を等間隔に分割した値ではなく、不等間隔に分割した間隔である(例えば、少なくとも一つのドップラ間隔が異なる)。例えば、ドップラシフト量DOPndmの各間隔は、ドップラ周波数範囲(例えば、-1/(2NTM×Tr) ≦ fd <1/(2NTM×Tr))を複数の送信アンテナ105の数Ntを時分割多重数NTMで除算した値(換言すると、ドップラ多重数NDM)に、1以上(例えば、δ)を加算した整数値で分割した間隔に設定されてよい。
なお、実施の形態1では、ドップラ多重数を送信アンテナ数Ntに等しい場合(すなわち、Nt=NDM)について説明した。これに対して、本実施の形態では、ドップラ多重に対して時分割多重を併用するので、ドップラ多重数NDMは、送信アンテナ数Ntより少ない多重数となる(例えば、Nt>NDM)。
したがって、本実施の形態では、ドップラシフト量DOPndmに対して、実施の形態1で用いた式(5)又は式(15)の中のNtをNDMで置き換えた式を用いる。また、時分割多重を行う送信周期(NTM×Tr)では同一の位相回転とするために、時分割多重を行う送信周期間(NTM×Tr)では同じ位相回転φndm(m)が繰り返して出力される。
すなわち、第ndm番目のドップラシフト部104は、入力された第m番目のチャープ信号に対して、異なるドップラシフト量DOP
ndmとなる、次式のような位相回転φ
ndm(m)を付与する。
ここで、Aは1又は‐1の正負の極性を与える係数である。また、δは1以上の正数である。また、φ0は初期位相であり、Δφ0は基準ドップラシフト位相である。なお、round(x)は実数値xに対し、四捨五入した整数値を出力するラウンド関数である。また、floor[x]は、実数x以下で最も近い整数を出力する演算子である。また、round(Ncode/(NDM+1))の項は、位相回転量を、ドップラ解析部209におけるドップラ周波数間隔の整数倍とする目的で導入されている。
または、本実施の形態において、式(46)の代わりに、次式の位相回転φ
ndm(m)を用いてもよい。
ここで、dpndmは位相回転を不等間隔とする成分である。例えば、dp1、dp2、 …、dpDM、-round(NC/NDM)/2<dpn<round(NC/NDM)/2の範囲の値であり、全てが同一の値ではでなく、少なくとも一つは異なる値の成分を含む。なお、round(NC/NDM)の項は、位相回転量を、ドップラ解析部209におけるドップラ周波数間隔の整数倍とする目的で導入されている。
なお、本実施の形態におけるドップラ多重分離部211eにおける動作は、実施の形態2においてドップラ多重と符号多重とを併用した場合のドップラ多重分離部211c(例えば、図9を参照)の動作におけるLOCをNTMで置き換えた場合の動作と同一であるため、その動作の説明を省略する。
以上のようにして、ドップラ多重分離部211eは、ドップラ多重信号を分離できる。
以上、ドップラ多重分離部211eの動作例について説明した。
図11において、方向推定部212eは、ドップラ多重分離部211eから入力される情報(例えば、距離インデックスfb_cfar及び分離インデックス情報(fdemul_DS#1, fdemul_DS#2, …,fdemul_DS#NDM)に基づいて、ターゲットの方向推定処理を行う。
例えば、方向推定部212eは、各ドップラ解析部209の出力から、距離インデックスfb_cfar、及び、分離インデックス情報(fdemul_DS#1, fdemul_DS#2, …,fdemul_DS#NDM)に対応する出力を抽出し、次式に示すような仮想受信アレー相関ベクトルh(fb_cfar, fdemul_DSx#1, fdemul_DSx#2, …, fdemul_DS#NDM)を生成し、方向推定処理を行う。
仮想受信アレー相関ベクトルh(f
b_cfar, f
demul_DSx#1, f
demul_DSx#2, …, f
demul_DS#NDM)は、送信アンテナ数Ntと受信アンテナ数Naとの積であるNt×Na個の要素を含む。仮想受信アレー相関ベクトルh(f
b_cfar, f
demul_DSx#1, f
demul_DSx#2, …, f
demul_DS#NDM)は、ターゲットからの反射波信号に対して各受信アンテナ202間の位相差に基づく方向推定を行う処理に用いる。ここで、z=1,…,Naである。なお、方向推定方法は、例えば、実施の形態1と同様の方法が適用されてよい。
式(48)において、h
cal[b]は、送信アレーアンテナ間及び受信アレーアンテナ間の位相偏差及び振幅偏差を補正するアレー補正値である。b=1,…, (Nt×Na)である。また、Txc
ntm(f
s)は、送信アンテナを時分割で切り替えることで、ドップラ周波数インデックスf
sによって異なる位相回転が発生するため、位相回転を補正し基準送信アンテナの位相に一致させる送信位相補正係数である。例えば、第1番目の時分割多重インデックス(ntm=1)を基準送信アンテナとした場合、次式となる。ここで、ntm=1,…, N
TMである。
以上のように、本実施の形態では、ドップラ多重と時分割多重とを併用する構成により、実施の形態1と同様の効果に加え、同時に多重送信する信号数を増大することができ、送信アンテナ数が増大したMIMOアレー構成への適応が可能になる。
なお、以上の説明において、ドップラ多重数をNDMとし、時分割多重数をNTMとし、送信アンテナ105の数Nt=NDM×NTMとなるドップラ多重数と時分割多重数を用いたが、これに限定されない。例えば、NDM個のドップラ多重信号に対し、時分割多重数を同一にせずに時分割多重数NTM以下としてもよい。
例えば、Nt=5の送信アンテナ105、ドップラ多重数をNDM=3、時分割多重数をNTM=2とした場合について説明する。この場合、3(=NDM)個のドップラシフト部104は、チャープ信号に対して、それぞれドップラシフト量DOP1、DOP2、DOP3を付与する。また、3(=NDM)個の送信切替部502は、ドップラシフト部104-1およびドップラシフト部104-2の出力に対して、2(=NTM)個の送信アンテナを切り替えて出力し、ドップラシフト部104-3の出力に対して、1(≦NTM)個の送信アンテナから出力する構成を用いる。換言すると、複数の送信アンテナ105から送信されるレーダ送信信号間で、レーダ送信信号に適用される時分割多重数NTMが異なる。この場合のレーダ受信部200eの処理は、これまで説明した処理(Nt=6の送信アンテナ105、ドップラ多重数をNDM=3、時分割多重数をすべてNTM=2とした場合の処理)と同様な処理で、Nt=5の送信アンテナからの送信信号を分離することができる。このように、NDM個のドップラ多重信号に対し、時分割多重数を同一にせずに時分割多重数NTM以下を用いることで、ドップラ多重数NDMを超えた送信アンテナ数(換言すれば多重送信数)の適用範囲を増やすことができる。例えば、ドップラ多重数をNDM=3、時分割多重数NTMを2以下とした場合、送信アンテナ数Nt(換言すれば多重送信数)は4、5、6の範囲で用いることが可能である。より一般的に記載すれば、送信アンテナ数Nt(換言すれば多重送信数)は、NDM+1≦Nt≦NDM×NTMの範囲で適用が可能である。
また、複数のドップラシフト部104の出力のうち、少なくとも1つのドップラシフト部104の出力に対して、送信切替部502を用いずに、送信アンテナ105に出力する構成を用いてもよい。このように複数のドップラシフト部104の出力のうち、少なくとも1つのドップラシフト部104の出力に対して、時分割多重を用いない信号を送信アンテナから出力する構成を用いることで、レーダ受信部200eにおいて、ドップラ解析部209の出力にドップラ折り返し信号が含まれるかどうかを検出できる。すなわち、ドップラ解析部209がサンプリング定理から導出される折り返しが発生しない最大ドップラ周波数は、±1/(2Loc×Tr)であるが、このように複数のドップラシフト部104の出力のうち、少なくとも1つのドップラシフト部104の出力に対して、時分割多重を用いない信号を送信アンテナから出力する構成を用いることで、ドップラ解析部209がサンプリング定理から導出される折り返しが発生しない最大ドップラ周波数を、±1/(2×Tr)とすることができ、曖昧性なく検出できるドップラ周波数範囲を拡大する効果も得ることができる。
なお、ドップラ多重と時分割多重とを併用する場合に、実施の形態1のバリエーション5と同様、擬似ランダムな符号系列が送信信号に乗算されてもよい。擬似ランダムな符号系列の符号長NLRCは、Ncode以下として、時分割周期毎にランダム符号要素インデックスをRC_INDEX(m)=floor[(m-1)/NTM]+1のように可変して、擬似ランダムな符号系列RCodeのランダム符号要素RC(RC_INDEX(m))を出力してもよい。
以上、本開示に係る一実施例について説明した。
[他の実施の形態]
(バリエーション7)
バリエーション7では、例えば、レーダ装置は、送信周期毎にドップラシフト量の各間隔を可変に設定し、送信アンテナに対するドップラ多重の割り当てを変える。
なお、バリエーション7に係るレーダ装置は、図1に示すレーダ装置10と基本構成が共通するので、図1を援用して説明する。例えば、バリエーション7では、図1に示すレーダ装置10において、ドップラシフト部104、ドップラ解析部209、CFAR部210及びドップラ多重分離部211の動作が実施の形態1と異なる。
例えば、ドップラ多重において、複数のターゲットのドップラピークの受信レベルがほぼ等しく、ドップラピークの間隔がドップラシフト量の間隔に一致する場合、ドップラ多重分離部211において分離判定ができなくなる可能性がある。
例えば、バリエーション3では、レーダ装置10の測位出力において、複数のターゲットをより確実に分離するために、レーダ観測毎にドップラシフト量が可変に設定される場合について説明した。
バリエーション7では、レーダ装置10の測位出力において、複数のターゲットをより確実に分離するために、送信周期毎にドップラシフト量の各間隔が可変に設定される場合について説明する。バリエーション7によれば、1つのターゲットに対して、複数の送信アンテナ105に対応するドップラピークの間隔が送信周期毎に異なるので、レーダ装置10は、複数のターゲットを1回のレーダ観測で分離しやすくなる。
以下、バリエーション7に係るドップラシフト部104において付与されるドップラシフト量の設定方法の一例について説明する。
ドップラシフト部104-1~104-Ntは、各々に入力されるチャープ信号に対して異なるドップラシフト量DOPnを付与する。ここで、n=1、…、Ntである。
また、ドップラシフト部104-1~104-Ntは、送信周期Tr毎にドップラシフト量DOPnを可変に設定する。例えば、ドップラシフト部104-1~104-Ntは、奇数番目の送信周期Tr毎のドップラシフト量DOPn
odd、及び、偶数番目の送信周期Tr毎のドップラシフト量DOPn
evenをそれぞれ設定する。
例えば、第n番目のドップラシフト部104は、次式に従って、入力された第m番目のチャープ信号に対して、奇数番目の送信周期Tr毎にドップラシフト量DOP
n
oddに対応する位相回転量φ
n(m)を付与し、偶数番目の送信周期Tr毎にドップラシフト量DOP
n
evenに対応する位相回転量φ
n(m)を付与する。
ここで、δodd及びδevenは1以上の正数であり、互いに異なる値に設定される。δodd及びδevenの設定により、奇数番目の送信周期Tr毎のドップラシフト量DOPn
oddと、偶数番目の送信周期Tr毎のドップラシフト量DOPn
evenとが異なる設定になる。換言すると、ドップラシフト量の各間隔が送信周期Tr毎に可変に設定される。
なお、位相回転量φnは、式(50)に示す値に限定されず、ドップラシフト量DOPn
odd及びドップラシフト量DOPn
evenの間隔が異なるような位相回転量であればよい。
また、ドップラシフト部104がレーダ送信信号(例えば、チャープ信号)に対して位相回転量を付与する際に、位相回転誤差が含まれる場合、ドップラ領域にスプリアスが発生する。ここで、例えば、スプリアスレベルがドップラピークレベルと比較して-20dB程度以下であれば、レーダ装置10におけるレーダ検出性能に顕著な劣化影響を与えない。そのため、位相回転時の位相回転誤差として、スプリアスレベルがドップラピークと比較して-20dB程度以下の範囲内(例えば、5°~10°程度の範囲)の位相回転誤差が含まれてもよい。なお、他の実施の形態(又はバリエーション)においても同様に、スプリアスレベルがドップラピークと比較して-20dB程度以下の範囲内(例えば、5°~10°程度の範囲)の位相回転誤差を含んでもよい。
図1において、ドップラ解析部209は、ビート周波数解析部208から出力される、NC回のチャープパルス送信によって得られるビート周波数応答RFTz(fb, 1)、RFTz(fb, 2)、…、RFTz(fb, NC)を用いて、距離インデックスfb毎にドップラ解析を行う。
バリエーション7では、レーダ送信信号(例えば、チャープ信号)に対して、奇数番目の送信周期Tr毎のドップラシフト量と、偶数番目の送信周期Tr毎のドップラシフト量とが異なる設定の位相回転φnが付与されている。そのため、ドップラ解析部209は、例えば、奇数番目の送信周期Tr毎のビート周波数応答を用いて、距離インデックスfb毎にドップラ解析を行う。同様に、ドップラ解析部209は、例えば、偶数番目の送信周期Tr毎のビート周波数応答を用いて、距離インデックスfb毎にドップラ解析を行う。
例えば、Ncが2のべき乗値である場合、ドップラ解析においてFFT処理を適用できる。この場合、FFTサイズはNc/2であり、ドップラ解析部209は、奇数番目又は偶数番目の送信周期Tr毎(換言すると、2Tr毎)に得られたデータに基づいてFFT処理を行う。このため、サンプリング定理から導出される折り返しが発生しない最大ドップラ周波数は±1/(4Tr)である。また、ドップラ周波数インデックスfsのドップラ周波数間隔は1/(Nc×Tr)であり、ドップラ周波数インデックスfsの範囲はfs = -Nc/4, …, 0, …, Nc/4-1である。
以下では、一例として、Ncが2のべき乗値である場合について説明する。なお、Ncが2のべき乗でない場合には、例えば、ゼロ埋めしたデータを含めることで2のべき乗個のデータサイズとしてFFT処理が可能である。また、ドップラ解析部209は、FFT処理の際に、Han窓又はHamming窓などの窓関数係数を乗算してもよい。窓関数を適用することでビート周波数ピーク周辺に発生するサイドローブを抑圧できる。
例えば、第z番目の信号処理部206における、奇数番目の送信周期Tr毎のビート周波数応答に対するドップラ解析部209の出力VFT
z
odd(f
b, f
s)、及び、偶数番目の送信周期Tr毎のビート周波数応答に対するドップラ解析部209の出力VFT
z
even(f
b, f
s)は、次式に示される。なお、jは虚数単位であり、z=1~Naである。
CFAR部210は、第1~第Na番目の信号処理部206のドップラ解析部209からの出力を用いて、CFAR処理(換言すると、適応的な閾値判定)を行い、ピーク信号を与える距離インデックスfb_cfar及びドップラ周波数インデックスfs_cfarを抽出する。
バリエーション7に係るCFAR部210は、例えば、奇数番目の送信周期Tr毎のビート周波数応答に対するドップラ解析部209の出力VFTz
odd(fb, fs)に対してCFAR処理を行うことにより、適応的に閾値を設定し、閾値よりも大きい受信電力となる距離インデックスfb_cfar
odd、ドップラ周波数インデックスfs_cfar
odd、及び、受信電力情報PowerFTodd(fb_cfar
odd, fs_cfar
odd)をドップラ多重分離部211に出力する。
同様に、バリエーション7に係るCFAR部210は、例えば、偶数番目の送信周期Tr毎のビート周波数応答に対するドップラ解析部209の出力VFTz
even(fb, fs)に対しCFAR処理を行うことにより、適応的に閾値を設定し、閾値よりも大きい受信電力となる距離インデックスfb_cfar
even、ドップラ周波数インデックスfs_cfar
even、及び、受信電力情報PowerFTeven(fb_cfar
even, fs_cfar
even)をドップラ多重分離部211に出力する。
ドップラ多重分離部211は、CFAR部210から入力される情報(例えば、奇数番目の送信周期Tr毎のビート周波数応答に対する距離インデックスfb_cfar
odd、ドップラ周波数インデックスfs_cfar
odd、及び、受信電力情報PowerFTodd(fb_cfar
odd, fs_cfar
odd)、及び、偶数番目の送信周期Tr毎のビート周波数応答に対する距離インデックスfb_cfar
even、ドップラ周波数インデックスfs_cfar
even、及び、受信電力情報PowerFTeven(fb_cfar
even, fs_cfar
even))に基づいて、各ドップラ解析部209からの出力を用いて、ドップラ多重送信された信号(以下、ドップラ多重信号と呼ぶ)から、各送信アンテナ105から送信される送信信号(換言すると、当該送信信号に対する反射波信号)を分離する。
ドップラ多重分離部211は、例えば、分離した信号に関する情報を、方向推定部212に出力する。分離した信号に関する情報には、例えば、分離した信号に対応する距離インデックスfb_cfar、及び、ドップラ周波数インデックス(以下、分離インデックス情報と呼ぶこともある)(fdemul_Tx#1, fdemul_Tx#2, …, fdemul_Tx#Nt)が含まれてよい。また、ドップラ多重分離部211は、各ドップラ解析部209からの出力を方向推定部212に出力する。
一例として、式(50)において、Nt=3、Δφ
0=0、φ
0=0、A=1、δ
odd=1、δ
even=2、N
Cを4の倍数とする場合、レーダ送信信号に対して次式のような位相回転量φ
n(m)が付与される。
また、ドップラ解析部209において、式(51)のFFT処理を行う場合、ドップラシフト量は、DOP1
odd=0、DOP1
even=0、DOP2
odd=1/(8Tr)、DOP2
even=1/(10Tr) 、DOP3
odd=1/(4Tr)、DOP3
even=1/(5Tr)となる。
このようなドップラシフト量を用いる場合、例えば、図12に示すように、測定する1つのターゲットドップラ周波数fd_TargetDopplerに対して、Nt個(図12では3つ)のドップラピークが発生する。なお、図12は、横軸にターゲットドップラ周波数を示し、縦軸にドップラ解析部209(FFT)の出力を示した場合のNt=3のドップラピークの変化を示した図である。
図12(a)は、奇数番目の送信周期Tr毎のビート周波数応答に対するドップラ解析部209の出力の一例を示し、図12(b)は、偶数番目の送信周期Tr毎のビート周波数応答に対するドップラ解析部209の出力の一例を示す。
図12(a)及び図12(b)では、測定する1つのターゲットドップラ周波数fd_TargetDopplerに対して、Nt個(図12では3つ)のドップラピークが発生するが、ドップラピークの間隔が異なる。例えば、図12(a)では、ドップラピークの間隔は、1/(8Tr)又は1/(4Tr)である。一方、例えば、図12(b)では、ドップラピークの間隔は、1/(10Tr)又は3/(10Tr)である。
したがって、同じ距離インデックスfbにおいて、2つのターゲットが存在する場合に、2つのターゲットのドップラ周波数の差が、例えば、奇数番目の送信周期Tr毎のドップラシフト量の間隔に一致する場合でも、偶数番目の送信周期Tr毎のドップラシフト量の間隔とは一致しないため、ドップラ多重分離部211は、2つのターゲットに対応する信号を分離して検出できる。
同様に、同じ距離インデックスfbにおいて、2つのターゲットが存在する場合に、2つのターゲットのドップラ周波数の差が、例えば、偶数番目の送信周期Tr毎のドップラシフト量の間隔に一致する場合でも、奇数番目の送信周期Tr毎のドップラシフト量の間隔とは一致しないため、ドップラ多重分離部211は、2つのターゲットに対応する信号を分離して検出できる。
よって、レーダ装置10では、1回のレーダ観測内において複数のターゲットを分離しやすくなる。
例えば、同じ距離インデックスfbにおいて、図13に示すように、ターゲット#1のドップラ周波数が0であり、ターゲット#2のドップラ周波数が1/(8Tr)である場合について説明する。
この場合、例えば、図13(a)に示すように、ターゲット#1、#2のドップラ周波数の差1/(8Tr)(換言すると、間隔)は、奇数番目の送信周期Tr毎のドップラシフト量の間隔(例えば、1/(8Tr))に一致する。このため、例えば、図13(a)に示すように、奇数番目の送信周期Tr毎のビート周波数応答に対するドップラ解析部209の出力では、ターゲット#1、#2のドップラピークが重なるため、ドップラ多重分離部211は、ターゲット#1、#2の信号を分離しにくくなる。
一方、例えば、図13(b)に示すように、ターゲット#1、#2のドップラ周波数の差1/(8Tr)(換言すると、間隔)は、偶数番目の送信周期Tr毎のドップラシフト量の間隔(例えば、1/(10Tr))とは一致しない。このため、例えば、図13(b)に示すように、偶数番目の送信周期Tr毎のビート周波数応答に対するドップラ解析部209の出力では、ターゲット#1、#2のドップラピークが重ならないため、ドップラ多重分離部211は、ターゲット#1、#2の信号を分離し、検出しやすくなる。
このように、レーダ装置10は、ドップラシフト量の間隔が異なる送信周期Trの何れかにおいて複数のターゲットに対応する信号を分離できる可能性が高くなる。これにより、レーダ装置10では、1回のレーダ観測内において複数のターゲットが分離しやすくなる。
以上のように、バリエーション7では、レーダ装置10は、送信周期Tr毎にドップラシフト量の間隔を可変に設定する。これにより、1つのターゲットに対する複数の送信アンテナ105に対応するドップラピークの間隔が送信周期毎に異なるので、レーダ装置10は、複数のターゲットを1回のレーダ観測で分離しやすくなる。
(バリエーション8)
バリエーション8では、例えば、レーダ装置は、送信周期毎にドップラシフト量を可変に設定し、送信アンテナに対するドップラ多重の割り当てを変える。
なお、バリエーション8に係るレーダ装置は、図1に示すレーダ装置10と基本構成が共通するので、図1を援用して説明する。例えば、バリエーション8では、図1に示すレーダ装置10において、ドップラシフト部104、ドップラ解析部209、CFAR部210及びドップラ多重分離部211の動作が実施の形態1と異なる。なお、バリエーション8に係るドップラ解析部209、CFAR部210及びドップラ多重分離部211の動作は、バリエーション7と同様の動作であるので、ここでは説明を省略する。
バリエーション8では、レーダ装置10の測位出力において、送信周期毎にドップラシフト量が可変に設定される場合について説明する。バリエーション8によれば、1つのターゲットに対して、複数の送信アンテナ105に対応するドップラピークの位置が送信周期毎に異なるので、レーダ装置10は、ドップラ領域で有色な干渉成分が存在する場合でも、1回のレーダ観測内においてターゲットを分離しやすくなる。
以下、バリエーション8に係るドップラシフト部104において付与されるドップラシフト量の設定方法の一例について説明する。
ドップラシフト部104-1~104-Ntは、各々に入力されるチャープ信号に対して異なるドップラシフト量DOPnを付与する。ここで、n=1、…、Ntである。
また、ドップラシフト部104-1~104-Ntは、送信周期Tr毎にドップラシフト量DOPnを可変に設定する。例えば、ドップラシフト部104-1~104-Ntは、奇数番目の送信周期Tr毎のドップラシフト量DOPn
odd、及び、偶数番目の送信周期Tr毎のドップラシフト量DOPn
evenをそれぞれ設定する。
例えば、第n番目のドップラシフト部104は、次式に従って、入力された第m番目のチャープ信号に対して、奇数番目の送信周期Tr毎にドップラシフト量DOP
n
oddに対応する位相回転量φ
n(m)を付与し、偶数番目の送信周期Tr毎にドップラシフト量DOP
n
evenに対応する位相回転量φ
n(m)を付与する。
ここで、δは1以上の正数である。式(55)のような位相回転φnが付与される。δodd及びδevenの設定により、奇数番目の送信周期Tr毎のドップラシフト量DOPn
oddと、偶数番目の送信周期Tr毎のドップラシフト量DOPn
evenとが異なる設定になる。換言すると、ドップラシフト量が送信周期Tr毎に可変に設定される。これにより、送信アンテナ105に対するドップラ多重の割り当てが送信周期Tr毎に可変に設定される。
なお、位相回転量φnは、式(55)に示す値に限定されず、ドップラシフト量DOPn
odd及びドップラシフト量DOPn
evenの位置(換言すれば、割り当て)が異なるように設定される位相回転量であればよい。
一例として、式(55)において、Nt=3、Δφ
0=0、φ
0=0、A=1、δ=1、N
Cを4の倍数とする場合、レーダ送信信号に対して次式のような位相回転φ
n(m)が付与される。
またドップラ解析部209において、式(51)のFFT処理を行う場合、ドップラシフト量は、DOP1
odd=0、DOP1
even=1/(8Tr)、DOP2
odd=1/(8Tr)、DOP2
even=1/(4Tr) 、DOP3
odd=1/(4Tr)、DOP3
even=-1/(8Tr)となる。
このようなドップラシフト量を用いる場合、例えば、図14に示すように、測定する1つのターゲットドップラ周波数fd_TargetDopplerに対して、Nt個(図14では3つ)のドップラピークが発生する。なお、図14は、横軸にターゲットドップラ周波数を示し、縦軸にドップラ解析部209(FFT)の出力を示した場合のNt=3のドップラピークの変化を示した図である。
図14(a)は、奇数番目の送信周期Tr毎のビート周波数応答に対するドップラ解析部209の出力の一例を示し、図14(b)は、偶数番目の送信周期Tr毎のビート周波数応答に対するドップラ解析部209の出力の一例を示す。
図14(a)及び図14(b)では、測定する1つのターゲットドップラ周波数fd_TargetDopplerに対して、Nt個(図14では3つ)のドップラピークが発生するが、ドップラピークの位置が異なる。例えば、図14(a)に示すドップラ解析部209の出力と、図14(b)に示すドップラ解析部209の出力とは、ドップラ領域において1/8Trずれている。
したがって、同じ距離インデックスfbにおいて、ドップラ領域で有色な干渉成分が存在する場合(換言すると、干渉成分が限定されたドップラ領域で発生する場合)に、例えば、奇数番目の送信周期及び偶数番目の送信周期の何れか一方の送信周期において、干渉成分が存在するドップラ領域にドップラピークが発生する場合でも、他方の送信周期において、干渉成分が存在するドップラ領域と異なるドップラ領域にドップラピークが発生する可能性を高くすることができる。よって、ドップラ多重分離部211は、1回のレーダ観測内において干渉影響を受けずに分離して検出しやすくなる。
例えば、同じ距離インデックスfbにおいて、図15に示すように、ドップラ領域で有色な干渉成分がドップラ周波数-1/(16Tr)~1/(16Tr)の範囲に存在する場合について説明する。また、図15では、一例として、ターゲット#1のドップラ周波数が0である場合について説明する。
この場合、例えば、図15(a)に示すように、奇数番目の送信周期Tr毎のビート周波数応答に対するドップラ解析部209の出力では、ターゲット#1のドップラピークの一部が、有色な干渉成分と重なるため、ドップラ多重分離部211は、ターゲット#1の信号を分離しにくくなる。
一方、例えば、図15(b)に示すように、偶数番目の送信周期Tr毎のビート周波数応答に対するドップラ解析部209の出力では、ターゲット#1のドップラピークは有色な干渉成分と重ならないため、ドップラ多重分離部211は、ターゲット#1の信号を分離しやすくなる。
このように、レーダ装置10は、ドップラシフト量(換言すると、ドップラ周波数範囲における位置)が異なる送信周期Trの何れかにおいて複数のターゲットに対応する信号を分離できる可能性が高くなる。これにより、レーダ装置10では、1回のレーダ観測内においてドップラ領域に有色な干渉成分が存在する場合でもターゲットを分離しやすくなる。
以上のように、バリエーション8では、レーダ装置10は、送信周期Tr毎にドップラシフト量を可変に設定する。これにより、1つのターゲットに対する複数の送信アンテナ105に対応するドップラピークの位置が送信周期毎に異なるので、レーダ装置10は、ドップラ領域で有色な干渉成分が存在する場合でも、1回のレーダ観測内においてターゲットを分離しやすくなる。
以上、バリエーション8について説明した。なお、バリエーション7及びバリエーション8を組み合わせてもよい。すなわち、1つのターゲットに対する複数の送信アンテナ105に対応するドップラピークの間隔及び位置を送信周期Tr毎に異ならせるように、ドップラシフト量(換言すると、位相回転量)が設定されてもよい。
本開示の一実施例に係るレーダ装置において、レーダ送信部及びレーダ受信部は、物理的に離れた場所に個別に配置されてもよい。また、本開示の一実施例に係るレーダ受信部において、方向推定部と、他の構成部とは、物理的に離れた場所に個別に配置されてもよい。
また、本開示の一実施例において用いた数値(例えば、送信アンテナ数Nt、受信アンテナ数Na、ドップラ多重数NDM、位相回転に関する値(δ、φ0、δ、Δφ0、dpn等))は一例であり、それらの値に限定されない。
本開示の一実施例に係るレーダ装置は、図示しないが、例えば、CPU(Central Processing Unit)、制御プログラムを格納したROM(Read Only Memory)等の記憶媒体、およびRAM(Random Access Memory)等の作業用メモリを有する。この場合、上記した各部の機能は、CPUが制御プログラムを実行することにより実現される。但し、レーダ装置のハードウェア構成は、かかる例に限定されない。例えば、レーダ装置の各機能部は、集積回路であるIC(Integrated Circuit)として実現されてもよい。各機能部は、個別に1チップ化されてもよいし、その一部または全部を含むように1チップ化されてもよい。
以上、図面を参照しながら各種の実施形態について説明したが、本開示はかかる例に限定されないことは言うまでもない。当業者であれば、特許請求の範囲に記載された範疇内において、各種の変更例または修正例に想到し得ることは明らかであり、それらについても当然に本開示の技術的範囲に属するものと了解される。また、開示の趣旨を逸脱しない範囲において、上記実施形態における各構成要素を任意に組み合わせてもよい。
上記各実施形態では、本開示はハードウェアを用いて構成する例にとって説明したが、本開示はハードウェアとの連携においてソフトウェアでも実現することも可能である。
また、上記各実施形態の説明に用いた各機能ブロックは、典型的には集積回路であるLSIとして実現される。集積回路は、上記実施の形態の説明に用いた各機能ブロックを制御し、入力端子と出力端子を備えてもよい。これらは個別に1チップ化されてもよいし、一部または全てを含むように1チップ化されてもよい。ここでは、LSIとしたが、集積度の違いにより、IC、システムLSI、スーパーLSI、ウルトラLSIと呼称されることもある。
また、集積回路化の手法はLSIに限るものではなく、専用回路または汎用プロセッサを用いて実現してもよい。LSI製造後に、プログラムすることが可能なFPGA(Field Programmable Gate Array)、LSI内部の回路セルの接続又は設定を再構成可能なリコンフィギュラブル プロセッサ(Reconfigurable Processor)を利用してもよい。
さらには、半導体技術の進歩又は派生する別技術により、LSIに置き換わる集積回路化の技術が登場すれば、当然、その技術を用いて機能ブロックを集積化してもよい。バイオ技術の適用等が可能性としてありえる。
<本開示のまとめ>
本開示の一実施例に係るレーダ装置は、送信信号を送信する複数の送信アンテナと、前記複数の送信アンテナから送信される前記送信信号にドップラシフト量を付与する回路と、を具備し、前記ドップラシフト量の各間隔は、ドップラ解析の対象となるドップラ周波数範囲を不等間隔に分割した間隔に設定される。
本開示の一実施例において、前記ドップラシフト量の各間隔は、前記ドップラ周波数範囲を前記複数の送信アンテナの数に1以上の整数を加算した値で分割した間隔である。
本開示の一実施例において、前記ドップラシフト量の各間隔は、前記ドップラ周波数範囲を前記複数の送信アンテナの数で分割した間隔にオフセットを加えた間隔である。
本開示の一実施例において、前記ドップラシフト量は、前記送信信号が送信されるフレーム毎に可変に設定される。
本開示の一実施例において、前記ドップラシフト量は、前記送信信号が送信される送信周期毎に可変に設定される。
本開示の一実施例において、前記ドップラシフト量の各間隔は、前記送信信号が送信される送信周期毎に可変に設定される。
本開示の一実施例において、前記回路は、前記送信信号に、疑似ランダムな符号系列を乗算する。
本開示の一実施例において、前記送信アンテナは、サブアレー構成である。
本開示の一実施例において、前記回路は、前記サブアレー構成の前記送信アンテナから送信される前記送信信号に同一のドップラシフト量を付与する。
本開示の一実施例において、前記回路は、更に、時分割送信及び符号分割送信の少なくとも1つを適用して、前記送信信号を送信する。
本開示の一実施例において、前記ドップラシフト量の各間隔は、前記ドップラ周波数範囲を前記複数の送信アンテナの数以下の値で分割した間隔である。
本開示の一実施例において、前記回路は、更に、符号分割送信を適用して、前記送信信号を送信し、前記ドップラシフト量の各間隔は、前記ドップラ周波数範囲を前記複数の送信アンテナの数を符号多重数で除算した値に、1以上を加算した整数値で分割した間隔である。
本開示の一実施例において、前記回路は、更に、符号分割送信を適用して、前記送信信号を送信し、前記複数の送信アンテナから送信される前記送信信号間で、前記送信信号に適用される符号分割多重数が異なる。
本開示の一実施例において、前記回路は、更に、時分割送信を適用して、前記送信信号を送信し、前記ドップラシフト量の各間隔は、前記ドップラ周波数範囲を前記複数の送信アンテナの数を時分割数で除算した値に、1以上を加算した整数値で分割した間隔である。
本開示の一実施例において、前記回路は、更に、時分割送信を適用して、前記送信信号を送信し、前記複数の送信アンテナから送信される前記送信信号間で、前記送信信号に適用される時分割多重数が異なる。
本開示の一実施例において、レーダ装置は、前記送信信号がターゲットに反射した反射波信号を受信する複数の受信アンテナと、前記ドップラ周波数範囲のうち、前記ドップラシフト量の各間隔に対応する範囲毎の前記反射波信号の受信電力を加算した電力加算値に対して閾値を用いて、前記反射波信号のピークを検出する受信回路、を更に具備する。
本開示の一実施例において、前記ドップラシフト量の各間隔は、前記ドップラ周波数範囲をドップラ多重数より多い数で分割した間隔であり、前記受信回路は、前記ドップラ周波数範囲において、前記検出されたピークのうち、前記受信電力の高い順に、前記ドップラ多重数に相当する数の第1のピークの受信レベルと、前記第1のピーク以外の第2のピークの受信レベルとの差が閾値以上の場合、前記第1のピークに基づいて、前記反射波信号から前記送信信号をそれぞれ分離する。
本開示の一実施例において、前記受信回路は、前記送信アンテナと、前記送信アンテナからそれぞれ送信される前記送信信号に付与される前記ドップラシフト量との関係に基づいて、前記反射波信号から、前記送信信号をそれぞれ分離する。