JPWO2019030893A1 - レーダ信号処理装置及びレーダ信号処理方法 - Google Patents
レーダ信号処理装置及びレーダ信号処理方法 Download PDFInfo
- Publication number
- JPWO2019030893A1 JPWO2019030893A1 JP2019531353A JP2019531353A JPWO2019030893A1 JP WO2019030893 A1 JPWO2019030893 A1 JP WO2019030893A1 JP 2019531353 A JP2019531353 A JP 2019531353A JP 2019531353 A JP2019531353 A JP 2019531353A JP WO2019030893 A1 JPWO2019030893 A1 JP WO2019030893A1
- Authority
- JP
- Japan
- Prior art keywords
- unit
- resolution
- cell
- target
- signal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/50—Systems of measurement based on relative movement of target
- G01S13/58—Velocity or trajectory determination systems; Sense-of-movement determination systems
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/66—Radar-tracking systems; Analogous systems
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/87—Combinations of radar systems, e.g. primary radar and secondary radar
Abstract
Description
しかし、レーダ信号処理装置は、例えば、2つのレーダの受信信号を振幅加算する場合、2つのレーダの観測領域が同一であるとは限らないため、一方のレーダの分解能セルと他方のレーダの分解能セルとの対応関係が分からない。
このため、レーダ信号処理装置は、同一の目標に係る信号の振幅加算を行う際、一方のレーダの分解能セルと他方のレーダの分解能セルとの対応付けを行う必要がある。
まず、このレーダ信号処理装置は、一方のレーダにおける複数の分解能セルのうち、任意の分解能セルを検定セルに設定する(処理1)。
次に、このレーダ信号処理装置は、目標の運動に関する速度の仮説を設定し、設定した仮説と、2つのレーダの設置位置とを用いて、他方のレーダにおける複数の分解能セルの中で、検定セルに対応する分解能セルを特定する(処理2)。
次に、このレーダ信号処理装置は、特定した分解能セルを中心とする一定領域内の複数の分解能セルを積分セル候補に設定する(処理3)。
次に、このレーダ信号処理装置は、検定セルに含まれている信号と積分セルに含まれている信号とを振幅加算する(処理5)。
このレーダ信号処理装置は、一方のレーダにおける複数の分解能セルのそれぞれを検定セルに設定して、各々の検定セルに含まれている信号と積分セルに含まれている信号とを振幅加算するまで、上記の(処理1)〜(処理5)を繰り返し実施する。
図1は、この発明の実施の形態1によるレーダ信号処理装置を示す構成図である。
図2は、この発明の実施の形態1によるレーダ信号処理装置を示すハードウェア構成図である。
図1及び図2において、信号処理部1は、第1の信号処理部1a及び第2の信号処理部1bを備えており、例えば、図2に示す信号処理回路21で実現される。
第1の信号処理部1aは、第1のレーダの受信信号をアナログ信号からデジタル信号に変換するアナログデジタル変換器を備えている。第1の信号処理部1aに入力される第1のレーダの受信信号が、デジタル信号である場合には、第1の信号処理部1aがアナログデジタル変換器を備えている必要はない。
第1の信号処理部1aは、第1のレーダのデジタルの受信信号から第1のレンジドップラマップを生成し、第1のレンジドップラマップを尤度算出部6に出力する処理を実施する。
第2の信号処理部1bは、第2のレーダのデジタルの受信信号から第2のレンジドップラマップを生成し、第2のレンジドップラマップを尤度算出部6に出力する処理を実施する。
仮説生成部2は、目標の3次元運動に関する仮説を生成する処理を実施する。
目標の3次元運動は、直交座標での目標の運動であり、例えば、目標までの距離、目標の方位角、目標の高度、目標の速度及び目標の進行方向で定義される。
運動予測部3は、仮説生成部2により生成された仮説を用いて、第1のレーダにより信号が受信された時刻である第1の受信時刻での目標の位置及び速度をそれぞれ予測し、予測した目標の位置及び速度を中心セル算出部5に出力する処理を実施する。
また、運動予測部3は、仮説生成部2により生成された仮説を用いて、第2のレーダにより信号が受信された時刻である第2の受信時刻での目標の位置及び速度をそれぞれ予測し、予測した目標の位置及び速度を中心セル算出部5に出力する処理を実施する。
データベース部4は、第1のレーダが存在している位置の座標及び第2のレーダが存在している位置の座標をそれぞれ格納している。
中心セル算出部5は、第1のレーダが存在している位置の座標を用いて、第1の信号処理部1aから出力された第1のレンジドップラマップにおける複数の分解能セルの中で、運動予測部3により予測された第1の受信時刻での目標の位置及び速度に対応する分解能セルを第1の中心セルとして算出する処理を実施する。
また、中心セル算出部5は、第2のレーダが存在している位置の座標を用いて、第2の信号処理部1bから出力された第2のレンジドップラマップにおける複数の分解能セルの中で、運動予測部3により予測された第2の受信時刻での目標の位置及び速度に対応する分解能セルを第2の中心セルとして算出する処理を実施する。
中心セル算出部5は、算出した第1及び第2の中心セルを尤度算出部6に出力する処理を実施する。
尤度算出部6は、中心セル算出部5から出力された第1の中心セルを含む第1の領域内の各々の分解能セルと、中心セル算出部5から出力された第2の中心セルを含む第2の領域内の各々の分解能セルとの間の尤度をそれぞれ算出する処理を実施する。
組み合わせ選択部8は、尤度算出部6によりそれぞれ算出された尤度に基づいて、第1の領域内の各々の分解能セルと第2の領域内の各々の分解能セルとの組み合わせの中から、少なくとも1つ以上の組み合わせを選択する処理を実施する。
振幅加算部9は、組み合わせ選択部8により選択された少なくとも1つ以上の組み合わせに係る分解能セルに含まれている信号同士を加算する処理を実施する。
ここで、記録回路24は、例えば、RAM(Random Access Memory)、ROM(Read Only Memory)、フラッシュメモリ、EPROM(Erasable Programmable Read Only Memory)、EEPROM(Electrically Erasable Programmable Read Only Memory)などの不揮発性又は揮発性の半導体メモリ、磁気ディスク、フレキシブルディスク、光ディスク、コンパクトディスク、ミニディスク、あるいは、DVD(Digital Versatile Disc)が該当する。
また、信号処理回路21、仮説生成回路22、運動予測回路23、中心セル算出回路25、尤度算出回路26及び信号加算回路27は、例えば、単一回路、複合回路、プログラム化したプロセッサ、並列プログラム化したプロセッサ、ASIC(Application Specific Integrated Circuit)、FPGA(Field−Programmable Gate Array)、または、これらを組み合わせたものが該当する。
ソフトウェア又はファームウェアはプログラムとして、コンピュータのメモリに格納される。コンピュータは、プログラムを実行するハードウェアを意味し、例えば、CPU(Central Processing Unit)、中央処理装置、処理装置、演算装置、マイクロプロセッサ、マイクロコンピュータ、プロセッサ、あるいは、DSP(Digital Signal Processor)が該当する。
レーダ信号処理装置がソフトウェア又はファームウェアなどで実現される場合、データベース部4をコンピュータのメモリ41上に構成するとともに、信号処理部1、仮説生成部2、運動予測部3、中心セル算出部5、尤度算出部6及び信号加算部7の処理手順をコンピュータに実行させるためのプログラムをメモリ41に格納し、コンピュータのプロセッサ42がメモリ41に格納されているプログラムを実行するようにすればよい。
図4は、レーダ信号処理装置がソフトウェア又はファームウェアなどで実現される場合の処理手順であるレーダ信号処理方法を示すフローチャートである。
第1の信号処理部1aは、第1のレーダの受信信号をアナログ信号からデジタル信号に変換する。
第1の信号処理部1aは、第1のレーダのデジタルの受信信号に対する公知の信号処理として、ビーム形成処理、不要波抑圧処理及び信号積分処理などを実施することで、第1のレンジドップラマップを生成する(図4のステップST1)。
第1の信号処理部1aは、生成した第1のレンジドップラマップを尤度算出部6に出力する。
第2の信号処理部1bは、第2のレーダのデジタルの受信信号に対する公知の信号処理として、ビーム形成処理、不要波抑圧処理及び信号積分処理などを実施することで、第2のレンジドップラマップを生成する(図4のステップST1)。
第2の信号処理部1bは、生成した第2のレンジドップラマップを尤度算出部6に出力する。
図5は、第1の信号処理部1a及び第2の信号処理部1bにより生成される第1及び第2のレンジドップラマップの一例を示す説明図である。
第1及び第2のレンジドップラマップは、図5に示すように、目標のドップラ速度と目標までの距離との関係を示す2次元マップである。
図5では、第1及び第2のレンジドップラマップのそれぞれが110個(=10×11個)の分解能セルを備える例を示しており、各々の分解能セルは、ドップラ分解能及び距離分解能で定義される大きさを有している。
目標の3次元運動は、直交座標での目標の運動であり、例えば、目標までの距離、目標の方位角、目標の高度、目標の速度及び目標の進行方向で定義される。
以下の式(1)は、仮説生成部2により生成される時刻tにおける仮説hi(t)の一例を示している。iは、仮説を識別する仮説番号である。ただし、この実施の形態1では、説明の簡単化のため、仮説生成部2により生成される仮説hi(t)の数は1つであるとする。このため、この実施の形態1では、i=1である。
式(1)において、riは、或る基準位置を原点とする北基準直交座標系における目標までの距離、Aiは、目標の方位角、hiは、目標の高度である。
また、Viは、目標の速度、φiは、目標の北を基準とする場合の水平面における進行方向、εiは、目標の進行方向における鉛直成分である。
また、運動予測部3は、仮説生成部2により生成された仮説hi(t)を用いて、第2のレーダにより信号が受信された時刻である第2の受信時刻T2での目標の位置及び速度をそれぞれ予測し、予測した目標の位置及び速度を中心セル算出部5に出力する(図4のステップST3)。
まず、運動予測部3は、仮説生成部2により生成された仮説hi(t)を用いて、以下の式(2)〜(7)に示すように、北基準直交座標系における時刻tでの目標の位置ベクトル及び速度ベクトルをそれぞれ算出する。
xi(t)は、北基準直交座標系における目標のx軸方向の位置ベクトル、yi(t)は、北基準直交座標系における目標のy軸方向の位置ベクトル、zi(t)は、北基準直交座標系における目標のz軸方向の位置ベクトルである。
Vxi(t)は、北基準直交座標系における目標のx軸方向の速度ベクトル、Vyi(t)は、北基準直交座標系における目標のy軸方向の速度ベクトル、Vzi(t)は、北基準直交座標系における目標のz軸方向の速度ベクトルである。
xi(t+ΔT1)は、第1の受信時刻T1での目標のx軸方向の位置ベクトル、yi(t+ΔT1)は、第1の受信時刻T1での目標のy軸方向の位置ベクトル、zi(t+ΔT1)は、第1の受信時刻T1での目標のz軸方向の位置ベクトルである。
Vxi(t+ΔT1)は、第1の受信時刻T1での目標のx軸方向の速度ベクトル、Vyi(t+ΔT1)は、第1の受信時刻T1での目標のy軸方向の速度ベクトル、Vzi(t+ΔT1)は、第1の受信時刻T1での目標のz軸方向の速度ベクトルである。
xi(t+ΔT2)は、第2の受信時刻T2での目標のx軸方向の位置ベクトル、yi(t+ΔT2)は、第2の受信時刻T2での目標のy軸方向の位置ベクトル、zi(t+ΔT2)は、第2の受信時刻T2での目標のz軸方向の位置ベクトルである。
Vxi(t+ΔT2)は、第2の受信時刻T2での目標のx軸方向の速度ベクトル、Vyi(t+ΔT2)は、第2の受信時刻T2での目標のy軸方向の速度ベクトル、Vzi(t+ΔT2)は、第2の受信時刻T2での目標のz軸方向の速度ベクトルである。
次に、中心セル算出部5は、以下の式(20)〜(22)に示すように、第1のレーダが存在している位置の座標を用いて、式(8)〜(10)に示す第1の受信時刻T1での目標の位置の座標を、第1のレーダが存在している位置からの相対座標に変換する。
式(20)〜(22)において、x1は、第1のレーダが存在しているx軸方向の位置の座標、y1は、第1のレーダが存在しているy軸方向の位置の座標、z1は、第1のレーダが存在しているz軸方向の位置の座標である。
xi,1(t+ΔT1)は、第1のレーダが存在している位置からの目標のx軸方向の相対座標、yi,1(t+ΔT1)は、第1のレーダが存在している位置からの目標のy軸方向の相対座標、zi,1(t+ΔT1)は、第1のレーダが存在している位置からの目標のz軸方向の相対座標である。
式(23)〜(25)において、x2は、第2のレーダが存在しているx軸方向の位置の座標、y2は、第2のレーダが存在しているy軸方向の位置の座標、z2は、第2のレーダが存在しているz軸方向の位置の座標である。
xi,2(t+ΔT2)は、第2のレーダが存在している位置からの目標のx軸方向の相対座標、yi,2(t+ΔT2)は、第2のレーダが存在している位置からの目標のy軸方向の相対座標、zi,2(t+ΔT2)は、第2のレーダが存在している位置からの目標のz軸方向の相対座標である。
具体的には、中心セル算出部5は、以下の式(26)に示すように、式(20)〜(22)に示す相対座標を用いて、第1の中心セルc1における距離Ri,1(t+ΔT1)を算出する。
また、中心セル算出部5は、以下の式(27)に示すように、式(11)〜(13)に示す第1の受信時刻T1での目標の速度ベクトル及び式(20)〜(22)に示す相対座標を用いて、第1の中心セルc1におけるドップラ速度VRi,1(t+ΔT1)を算出する。
式(27)において、PFRは、第1のレーダにおけるパルス繰り返し周波数である。
mod(a,b)は、aをbで除算したときの余りを算出する関数である。
中心セル算出部5は、算出した第1の中心セルc1(R,VR)を尤度算出部6に出力する。
具体的には、中心セル算出部5は、以下の式(28)に示すように、式(23)〜(25)に示す相対座標を用いて、第2の中心セルc2おける距離Ri,2(t+ΔT2)を算出する。
また、中心セル算出部5は、以下の式(29)に示すように、式(17)〜(19)に示す第2の受信時刻T2での目標の速度ベクトル及び式(23)〜(25)に示す相対座標を用いて、第2の中心セルc2におけるドップラ速度VRi,2(t+ΔT2)を算出する。
中心セル算出部5は、算出した第2の中心セルc2(R,VR)を尤度算出部6に出力する。
図6は、尤度算出部6により設定される第1の領域及び第2の領域を示す説明図である。
図6の例では、第1の領域は、第1の中心セルc1(R,VR)が中心に配置されている15個(=5×3個)の分解能セルを含む領域である。
また、第2の領域は、第2の中心セルc2(R,VR)が中心に配置されている9個(=3×3個)の分解能セルを含む領域である。
第1の領域のサイズは、第1のレンジドップラマップのサイズよりも小さければよく、どのようなサイズでもよい。また、第2の領域のサイズは、第2のレンジドップラマップのサイズよりも小さければよく、どのようなサイズでもよい。
以下、第1の領域に含まれている各々の分解能セルをc1,(a,b)、第2の領域に含まれている各々の分解能セルをc2,(d,e)で表す。
aは、第1の領域のVR方向の分解能セルを示す変数、bは、第1の領域のR方向の分解能セルを示す変数である。
dは、第2の領域のVR方向の分解能セルを示す変数、eは、第2の領域のR方向の分解能セルを示す変数である。
図6の例では、第1の領域が15個の分解能セルc1,(a,b)を含み、第2の領域が9個の分解能セルc2,(d,e)を含んでいるので、尤度算出部6は、全部で135個(=15×9個)の尤度I(a,b)−(d,e)を算出する。
z1は、分解能セルc1,(a,b)が有する振幅値、z2は、分解能セルc2,(d,e)が有する振幅値、zバーは、第1の領域内の複数の分解能セルc1,(a,b)及び第2の領域内の複数の分解能セルc2,(d,e)における平均的な振幅値の期待値である。明細書中では、電子出願の都合上、文字の上“−”の記号を付することができないため、zバーのように表記している。
σは、分解能セルc1,(a,b)及び分解能セルc2,(d,e)に含まれている信号のバラツキを表すパラメータであり、事前に設定される。I0は、事前に設定された初期値である。
方位角Aの座標は、以下の式(35)に示す方位角Ai,1(t+ΔT1)を方位角の分解能で除算することで得られ、仰角Eの座標は、以下の式(36)に示す仰角Ei,1(t+ΔT1)を仰角の分解能で除算することで得られる。
方位角Aの座標は、以下の式(37)に示す方位角Ai,2(t+ΔT2)を方位角の分解能で除算することで得られ、仰角Eの座標は、以下の式(38)に示す仰角Ei,2(t+ΔT2)を仰角の分解能で除算することで得られる。
組み合わせ選択部8は、尤度I(a,b)−(d,e)の比較結果に基づいて、第1の領域内の各々の分解能セルc1,(a,b)と第2の領域内の各々の分解能セルc2,(d,e)との組み合わせの中から、少なくとも1つ以上の組み合わせを選択する(図4のステップST6)。
例えば、組み合わせ選択部8は、第1の領域内の各々の分解能セルc1,(a,b)と第2の領域内の各々の分解能セルc2,(d,e)との組み合わせの中から、尤度I(a,b)−(d,e)が最大の組み合わせを選択する。
あるいは、組み合わせ選択部8は、第1の領域内の各々の分解能セルc1,(a,b)と第2の領域内の各々の分解能セルc2,(d,e)との組み合わせの中から、尤度I(a,b)−(d,e)が高い上位N個の組み合わせを選択する。N個は、事前に設定される値である。
あるいは、組み合わせ選択部8は、第1の領域内の各々の分解能セルc1,(a,b)と第2の領域内の各々の分解能セルc2,(d,e)との組み合わせの中から、尤度I(a,b)−(d,e)が事前に設定された閾値以上の組み合わせを選択する。
振幅加算部9は、加算した信号を、目標の検出処理等を行う後段の処理部に出力する。
組み合わせ選択部8により選択された組み合わせが、例えば、分解能セルc1,(3,2)と分解能セルc2,(3,2)との組み合わせであれば、分解能セルc1,(3,2)に含まれている信号と分解能セルc2,(3,2)に含まれている信号とを加算する。このとき、信号同士の加算は、加算対象の信号の振幅値だけを加算するようにしてもよい。
また、組み合わせ選択部8により選択された組み合わせが、例えば、分解能セルc1,(3,2)と分解能セルc2,(3,2)との組み合わせと、分解能セルc1,(4,3)と分解能セルc2,(3,3)との組み合わせであれば、分解能セルc1,(3,2)に含まれている信号と、分解能セルc2,(3,2)に含まれている信号と、分解能セルc1,(4,3)に含まれている信号と、分解能セルc2,(3,3)に含まれている信号とを加算する。
これにより、特許文献1に開示されているレーダ信号処理装置と比べて、目標に係る信号同士を加算するまでに要する処理の負荷を大幅に削減することができる。
また、第2の信号処理部1bが、第2のレーダの受信信号として、1つの受信ビームの受信信号を取得する例を示しているが、第2のレーダの受信信号として、複数の受信ビームの受信信号を取得するようにしてもよい。
図7の例では、第2の信号処理部1bによって、受信ビームb−1についての第2のレンジドップラマップと、受信ビームbについての第2のレンジドップラマップと、受信ビームb+1についての第2のレンジドップラマップとが生成されている。
この場合、尤度算出部6は、中心セル算出部5から出力された第1の中心セルc1(R,VR)を含む第1の領域を設定する。
また、尤度算出部6は、図7に示すように、受信ビームb−1の第2のレンジドップラマップにおいて、中心セル算出部5から出力された第2の中心セルc2(R,VR)を含む第2の領域を設定する。
また、尤度算出部6は、図7に示すように、受信ビームbの第2のレンジドップラマップにおいて、第2の中心セルc2(R,VR)を含む第2の領域を設定し、受信ビームb+1の第2のレンジドップラマップにおいて、第2の中心セルc2(R,VR)を含む第2の領域を設定する。
また、尤度算出部6は、第1の領域内の各々の分解能セルc1,(a,b)と、受信ビームbにおける第2の領域内の各々の分解能セルc2,(d,e)との間の尤度I(a,b)−(d,e)をそれぞれ算出する。
また、尤度算出部6は、第1の領域内の各々の分解能セルc1,(a,b)と、受信ビームb+1における第2の領域内の各々の分解能セルc2,(d,e)との間の尤度I(a,b)−(d,e)をそれぞれ算出する。
図7の例では、第1の領域が15個の分解能セルc1,(a,b)を含み、3つの第2の領域が15個の分解能セルc2,(d,e)を含んでいるので、尤度算出部6は、全部で675個(=15×(15×3)個)の尤度I(a,b)−(d,e)を算出する。
例えば、組み合わせ選択部8は、第1の領域内の各々の分解能セルc1,(a,b)と、3つの第2の領域内の各々の分解能セルc2,(d,e)との組み合わせの中から、尤度I(a,b)−(d,e)が最大の組み合わせを選択する。
あるいは、組み合わせ選択部8は、第1の領域内の各々の分解能セルc1,(a,b)と、3つの第2の領域内の各々の分解能セルc2,(d,e)との組み合わせの中から、尤度I(a,b)−(d,e)が高い上位N個の組み合わせを選択する。N個は、事前に設定される値である。
あるいは、組み合わせ選択部8は、第1の領域内の各々の分解能セルc1,(a,b)と、3つの第2の領域内の各々の分解能セルc2,(d,e)との組み合わせの中から、尤度I(a,b)−(d,e)が事前に設定された閾値以上の組み合わせを選択する。
振幅加算部9は、組み合わせ選択部8により選択された少なくとも1つ以上の組み合わせに係る分解能セルに含まれている信号同士を加算する。
この実施の形態1では、レーダ信号処理装置が、第1及び第2のレーダの受信信号のそれぞれから第1及び第2のレンジドップラマップを生成する信号処理部1を備えている例を示している。
この実施の形態2では、N(Nは3以上の整数)個のレーダの受信信号のそれぞれからN個のレンジドップラマップを生成する信号処理部1を備えているレーダ信号処理装置について説明する。
第Nの信号処理部1Nは、第Nのレーダの受信信号をアナログ信号からデジタル信号に変換するアナログデジタル変換器を備えている。第Nの信号処理部1Nに入力される第Nのレーダの受信信号が、デジタル信号である場合には、第Nの信号処理部1Nがアナログデジタル変換器を備えている必要はない。
第Nの信号処理部1Nは、第Nのレーダのデジタルの受信信号から第Nのレンジドップラマップを生成し、第Nのレンジドップラマップを尤度算出部6に出力する処理を実施する。
図8では、第2の信号処理部1bと第Nの信号処理部1Nとの間の信号処理部を省略しているが、第3の信号処理部,・・・,第(N−1)の信号処理部などを実装している例が考えられる。
第3の信号処理部,・・・,第(N−1)の信号処理部の処理内容は、第Nの信号処理部1Nの処理内容と同様であるため説明を省略する。
第Nの信号処理部1Nは、第Nのレーダの受信信号をアナログ信号からデジタル信号に変換する。
第Nの信号処理部1Nは、第1の信号処理部1a及び第2の信号処理部1bと同様の処理を実施することで、第Nのレーダのデジタルの受信信号から第Nのレンジドップラマップを生成する。
第Nの信号処理部1Nは、生成した第Nのレンジドップラマップを尤度算出部6に出力する。
また、運動予測部3は、同様の方法で、第Nのレーダにより信号が受信された時刻である第Nの受信時刻TNでの目標の位置及び速度をそれぞれ予測する。
中心セル算出部5は、上記実施の形態1と同様の方法で、第1のレーダが存在している位置の座標を用いて、第1のレンジドップラマップにおける複数の分解能セルの中で、運動予測部3により予測された第1の受信時刻T1での目標の位置及び速度に対応する分解能セルを第1の中心セルc1(R,VR)として算出する。
中心セル算出部5は、上記実施の形態1と同様の方法で、第2のレーダが存在している位置の座標を用いて、第2のレンジドップラマップにおける複数の分解能セルの中で、運動予測部3により予測された第2の受信時刻T2での目標の位置及び速度に対応する分解能セルを第2の中心セルc2(R,VR)として算出する。
中心セル算出部5は、同様の方法で、第Nのレーダが存在している位置の座標を用いて、第Nのレンジドップラマップにおける複数の分解能セルの中で、運動予測部3により予測された第Nの受信時刻TNでの目標の位置及び速度に対応する分解能セルを第Nの中心セルcN(R,VR)として算出する。
中心セル算出部5は、第1の中心セルc1(R,VR)〜第Nの中心セルcN(R,VR)を尤度算出部6に出力する。
また、尤度算出部6は、中心セル算出部5から出力された第Nの中心セルcN(R,VR)を含む第Nの領域を設定する。
この実施の形態2では、第Nの領域に含まれている各々の分解能セルをcN,(f,g)で表す。
また、尤度算出部6は、同様の方法で、第1の中心セルc1(R,VR)を含む第1の領域内の各々の分解能セルc1,(a,b)と、中心セル算出部5から出力された第Nの中心セルcN(R,VR)を含む第Nの領域内の各々の分解能セルcN,(f,g)との間の尤度I(a,b)−(f,g)をそれぞれ算出する。
組み合わせ選択部8は、尤度I(a,b)−(d,e),I(a,b)−(f,g)の比較結果に基づいて、第1の領域内の各々の分解能セルc1,(a,b)と第2の領域内の各々の分解能セルc2,(d,e)との組み合わせ及び第1の領域内の各々の分解能セルc1,(a,b)と第Nの領域内の各々の分解能セルcN,(f,g)との組み合わせの中から、少なくとも1つ以上の組み合わせを選択する。
例えば、組み合わせ選択部8は、第1の領域内の各々の分解能セルc1,(a,b)と第2の領域内の各々の分解能セルc2,(d,e)との組み合わせ及び第1の領域内の各々の分解能セルc1,(a,b)と第Nの領域内の各々の分解能セルcN,(f,g)との組み合わせの中から、尤度が最大の組み合わせを選択する。
あるいは、組み合わせ選択部8は、第1の領域内の各々の分解能セルc1,(a,b)と第2の領域内の各々の分解能セルc2,(d,e)との組み合わせ及び第1の領域内の各々の分解能セルc1,(a,b)と第Nの領域内の各々の分解能セルcN,(f,g)との組み合わせの中から、尤度が高い上位N個の組み合わせを選択する。
あるいは、組み合わせ選択部8は、第1の領域内の各々の分解能セルc1,(a,b)と第2の領域内の各々の分解能セルc2,(d,e)との組み合わせ及び第1の領域内の各々の分解能セルc1,(a,b)と第Nの領域内の各々の分解能セルcN,(f,g)との組み合わせの中から、尤度が事前に設定された閾値以上の組み合わせを選択する。
信号加算部7の振幅加算部9は、組み合わせ選択部8により選択された少なくとも1つ以上の組み合わせに係る分解能セルに含まれている信号同士を加算する。
上記実施の形態1では、仮説生成部2が1つの仮説を生成する例を示している。
この実施の形態3では、仮説生成部2が複数の仮説をする例を説明する。
図10は、この発明の実施の形態3によるレーダ信号処理装置を示すハードウェア構成図である。
図9及び図10において、図1及び図2と同一符号は同一または相当部分を示すので説明を省略する。
仮説生成部11は、例えば、図10に示す仮説生成回路31で実現される。
仮説生成部11は、目標の3次元運動に関する複数の仮説を生成する処理を実施する。
仮説生成部11は、例えば、目標までの距離、目標の方位角、目標の高度、目標の速度及び目標の進行方向のそれぞれを乱数で変えながら、距離、方位角、高度、速度及び進行方向に基づいて、目標の3次元運動に関する仮説を生成する。
振幅加算部13は、仮説生成部11により生成される仮説単位に、組み合わせ選択部8により選択された少なくとも1つ以上の組み合わせに係る分解能セルに含まれている信号同士を加算する処理を実施する。
また、振幅加算部13は、仮説単位で加算した信号を互いに比較し、加算した信号の比較結果に基づいて、仮説単位のいずれかの加算信号を選択する処理を実施する。
ここで、信号処理回路21、仮説生成回路31、運動予測回路23、中心セル算出回路25、尤度算出回路26及び信号加算回路32は、例えば、単一回路、複合回路、プログラム化したプロセッサ、並列プログラム化したプロセッサ、ASIC、FPGA、または、これらを組み合わせたものが該当する。
レーダ信号処理装置がソフトウェア又はファームウェアなどで実現される場合、データベース部4を図3に示すコンピュータのメモリ41上に構成するとともに、信号処理部1、仮説生成部11、運動予測部3、中心セル算出部5、尤度算出部6及び信号加算部12の処理手順をコンピュータに実行させるためのプログラムをメモリ41に格納し、コンピュータのプロセッサ42がメモリ41に格納されているプログラムを実行するようにすればよい。
図11は、レーダ信号処理装置がソフトウェア又はファームウェアなどで実現される場合の処理手順であるレーダ信号処理方法を示すフローチャートである。
第1の信号処理部1aは、上記実施の形態1と同様に、第1のレーダの受信信号をアナログ信号からデジタル信号に変換する。
第1の信号処理部1aは、上記実施の形態1と同様に、第1のレーダのデジタルの受信信号に対する公知の信号処理を実施することで、第1のレンジドップラマップを生成する(図11のステップST11)。
第1の信号処理部1aは、生成した第1のレンジドップラマップを尤度算出部6に出力する。
第2の信号処理部1bは、上記実施の形態1と同様に、第2のレーダのデジタルの受信信号に対する公知の信号処理を実施することで、第2のレンジドップラマップを生成する(図11のステップST11)。
第2の信号処理部1bは、生成した第2のレンジドップラマップを尤度算出部6に出力する。
仮説生成部11は、仮説hi(t)を生成する際、例えば、目標までの距離r、目標の方位角A、目標の高度h、目標の速度V、目標の進行方向φ及び進行方向の鉛直成分εのそれぞれの平均値と上下限範囲を取得する。
仮説生成部11は、それぞれの平均値を中心にして、それぞれの上下限範囲内で、距離r、方位角A、高度h、速度V、進行方向φ及び鉛直成分εのそれぞれの値を乱数によって決定することで、仮説hi(t)を生成する。
乱数として、例えば、正規乱数または一様乱数を用いることができる。
また、運動予測部3は、上記実施の形態1と同様に、仮説生成部11により生成された仮説hi(t)を用いて、第2のレーダにより信号が受信された時刻である第2の受信時刻T2での目標の位置及び速度をそれぞれ予測し、予測した目標の位置及び速度を中心セル算出部5に出力する(図11のステップST13)。
中心セル算出部5は、上記実施の形態1に、第2のレンジドップラマップにおける複数の分解能セルの中で、運動予測部3により予測された第2の受信時刻T2での目標の位置及び速度に対応する分解能セルを第2の中心セルc2(R,VR)として算出する(図11のステップST14)。
組み合わせ選択部8は、上記実施の形態1と同様に、尤度I(a,b)−(d,e)の比較結果に基づいて、第1の領域内の各々の分解能セルc1,(a,b)と第2の領域内の各々の分解能セルc2,(d,e)との組み合わせの中から、少なくとも1つ以上の組み合わせを選択する(図11のステップST16)。
信号加算部7の振幅加算部9は、上記実施の形態1と同様に、組み合わせ選択部8により選択された少なくとも1つ以上の組み合わせに係る分解能セルに含まれている信号同士を加算する(図11のステップST17)。
仮説hi(t)の生成は、前回決定した距離r、方位角A、高度h、速度V、進行方向φ及び鉛直成分εのそれぞれの値を乱数で振りなおすことで生成する。
以下、ステップST13〜ST18の処理が実施される。既に生成した仮説hi(t)の個数が規定数Jに到達するまで、ステップST12〜ST18の処理が繰り返し実施される。
振幅加算部13は、J個の加算信号の比較結果に基づいて、J個の加算信号の中から、最大の加算信号を選択する(図11のステップST19)。
あるいは、振幅加算部13は、J個の加算信号の中から、大きい上位N個の加算信号を選択する。
あるいは、振幅加算部13は、J個の加算信号の中から、事前に設定された閾値以上の加算信号を選択する。
振幅加算部13は、選択した少なくとも1つ以上の加算信号を、目標の検出処理等を行う後段の処理部に出力する。
また、上記実施の形態1よりも、目標以外の信号を加算してしまう確率を下げて、後段の処理部での目標の検出精度を高めることができる。
上記実施の形態1では、組み合わせ選択部8が、第1の領域内の各々の分解能セルc1,(a,b)と第2の領域内の各々の分解能セルc2,(d,e)との組み合わせの中から、尤度I(a,b)−(d,e)が最大の組み合わせ、あるいは、尤度I(a,b)−(d,e)が高い上位N個の組み合わせ等を選択する例を示している。
しかし、振幅加算部9が、組み合わせ選択部8により選択された1つ以上の組み合わせに係る分解能セルに含まれている信号の中には、信号対雑音電力比が小さい信号が含まれている場合もある。
この実施の形態4におけるレーダ信号処理装置の構成図は、上記実施の形態1における図1の構成図と同様である。
式(39)において、vsは、組み合わせ選択部8により選択された組み合わせに係る分解能セルに含まれている信号同士を加算した加算信号である。
vnは、雑音電力であり、例えば、第1のレンジドップラマップの雑音領域から計算される振幅の平均値である。
図12は、第1のレンジドップラマップの雑音領域を示す説明図である。
図12の例では、図中、左上の領域と右上の領域を雑音領域としている。
振幅加算部9は、加算信号の信号対雑音電力比SNRが閾値SNRth以上であれば、加算信号を有効として、後段の処理部に出力する。
Claims (6)
- 第1及び第2のレーダの受信信号のそれぞれから第1及び第2のレンジドップラマップを生成する信号処理部と、
目標の3次元運動に関する仮説を生成する仮説生成部と、
前記仮説を用いて、前記第1及び第2のレーダによる信号の受信時刻である第1及び第2の受信時刻での前記目標の位置及び速度をそれぞれ予測する運動予測部と、
前記第1及び第2のレンジドップラマップにおける複数の分解能セルの中で、前記運動予測部により予測された第1及び第2の受信時刻での目標の位置及び速度に対応する分解能セルのそれぞれを第1及び第2の中心セルとして算出する中心セル算出部と、
前記第1の中心セルを含む第1の領域内の各々の分解能セルと、前記第2の中心セルを含む第2の領域内の各々の分解能セルとの間の尤度をそれぞれ算出する尤度算出部と、
前記尤度算出部によりそれぞれ算出された尤度に基づいて、前記第1の領域内の各々の分解能セルと前記第2の領域内の各々の分解能セルとの組み合わせの中から、少なくとも1つ以上の組み合わせを選択し、選択した組み合わせに係る分解能セルに含まれている信号同士を加算する信号加算部と
を備えたレーダ信号処理装置。 - 前記信号加算部は、前記第1の領域内の各々の分解能セルと前記第2の領域内の各々の分解能セルとの組み合わせの中から、前記尤度算出部により算出された尤度が高い上位N個の組み合わせを選択し、選択した上位N個の組み合わせに係る分解能セルに含まれている信号同士を加算することを特徴とする請求項1記載のレーダ信号処理装置。
- 前記仮説生成部は、目標の3次元運動に関する複数の仮説をそれぞれ生成し、
前記運動予測部は、前記仮説生成部により生成された各々の仮説を用いて、前記目標の位置及び速度をそれぞれ予測し、
前記中心セル算出部は、前記仮説生成部により生成される仮説単位に、前記第1及び第2の中心セルを算出し、
前記尤度算出部は、前記仮説生成部により生成される仮説単位に、前記第1の中心セルを含む第1の領域内の各々の分解能セルと、前記第2の中心セルを含む第2の領域内の各々の分解能セルとの間の尤度をそれぞれ算出し、
前記信号加算部は、前記仮説生成部により生成される仮説単位に、選択した組み合わせに係る分解能セルに含まれている信号同士を加算し、仮説単位で加算した信号を互いに比較し、前記加算した信号の比較結果に基づいて、仮説単位のいずれかの加算信号を選択することを特徴とする請求項1記載のレーダ信号処理装置。 - 前記仮説生成部は、目標までの距離、目標の方位角、目標の高度、目標の速度及び目標の進行方向のそれぞれを乱数で変えながら、前記距離、前記方位角、前記高度、前記速度及び前記進行方向に基づいて、目標の3次元運動に関する仮説を生成することを特徴とする請求項3記載のレーダ信号処理装置。
- 前記信号加算部は、選択した組み合わせに係る分解能セルに含まれている信号同士を加算して、加算した信号の信号対雑音電力比を算出し、前記信号対雑音電力比が閾値未満であれば、前記加算した信号を無効とし、前記信号対雑音電力比が閾値以上であれば、前記加算した信号を有効とすることを特徴とする請求項1記載のレーダ信号処理装置。
- 信号処理部が、第1及び第2のレーダの受信信号のそれぞれから第1及び第2のレンジドップラマップを生成し、
仮説生成部が、目標の3次元運動に関する仮説を生成し、
運動予測部が、前記仮説を用いて、前記第1及び第2のレーダによる信号の受信時刻である第1及び第2の受信時刻での前記目標の位置及び速度をそれぞれ予測し、
中心セル算出部が、前記第1及び第2のレンジドップラマップにおける複数の分解能セルの中で、前記運動予測部により予測された第1及び第2の受信時刻での目標の位置及び速度に対応する分解能セルのそれぞれを第1及び第2の中心セルとして算出し、
尤度算出部が、前記第1の中心セルを含む第1の領域内の各々の分解能セルと、前記第2の中心セルを含む第2の領域内の各々の分解能セルとの間の尤度をそれぞれ算出し、
信号加算部が、前記尤度算出部によりそれぞれ算出された尤度に基づいて、前記第1の領域内の各々の分解能セルと前記第2の領域内の各々の分解能セルとの組み合わせの中から、少なくとも1つ以上の組み合わせを選択し、選択した組み合わせに係る分解能セルに含まれている信号同士を加算する
レーダ信号処理方法。
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/JP2017/029084 WO2019030893A1 (ja) | 2017-08-10 | 2017-08-10 | レーダ信号処理装置及びレーダ信号処理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP6587784B2 JP6587784B2 (ja) | 2019-10-09 |
JPWO2019030893A1 true JPWO2019030893A1 (ja) | 2019-11-07 |
Family
ID=65271941
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2019531353A Active JP6587784B2 (ja) | 2017-08-10 | 2017-08-10 | レーダ信号処理装置及びレーダ信号処理方法 |
Country Status (2)
Country | Link |
---|---|
JP (1) | JP6587784B2 (ja) |
WO (1) | WO2019030893A1 (ja) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2021185350A (ja) * | 2020-05-25 | 2021-12-09 | ソニーセミコンダクタソリューションズ株式会社 | 信号処理装置、および信号処理方法、並びにプログラム |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006003348A (ja) * | 2004-05-20 | 2006-01-05 | Mitsubishi Electric Corp | レーダ装置 |
JP2012194044A (ja) * | 2011-03-16 | 2012-10-11 | Mitsubishi Electric Corp | レーダ装置 |
JP2013130473A (ja) * | 2011-12-21 | 2013-07-04 | Mitsubishi Electric Corp | レーダ信号処理装置およびレーダ装置 |
WO2015173891A1 (ja) * | 2014-05-13 | 2015-11-19 | 三菱電機株式会社 | レーダ装置 |
JP2016017929A (ja) * | 2014-07-11 | 2016-02-01 | 三菱電機株式会社 | レーダ装置 |
JP2016017928A (ja) * | 2014-07-11 | 2016-02-01 | 三菱電機株式会社 | レーダ装置 |
JP2016183953A (ja) * | 2015-03-25 | 2016-10-20 | パナソニック株式会社 | 物体検出装置および物体検出方法 |
WO2016199296A1 (ja) * | 2015-06-12 | 2016-12-15 | 三菱電機株式会社 | レーダ装置 |
US20170102458A1 (en) * | 2014-06-09 | 2017-04-13 | Src, Inc. | Multiplatform GMTI Radar With Adaptive Clutter Suppression |
JP2017125802A (ja) * | 2016-01-15 | 2017-07-20 | パナソニック株式会社 | 物体検出装置および物体検出方法 |
-
2017
- 2017-08-10 JP JP2019531353A patent/JP6587784B2/ja active Active
- 2017-08-10 WO PCT/JP2017/029084 patent/WO2019030893A1/ja active Application Filing
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006003348A (ja) * | 2004-05-20 | 2006-01-05 | Mitsubishi Electric Corp | レーダ装置 |
JP2012194044A (ja) * | 2011-03-16 | 2012-10-11 | Mitsubishi Electric Corp | レーダ装置 |
JP2013130473A (ja) * | 2011-12-21 | 2013-07-04 | Mitsubishi Electric Corp | レーダ信号処理装置およびレーダ装置 |
WO2015173891A1 (ja) * | 2014-05-13 | 2015-11-19 | 三菱電機株式会社 | レーダ装置 |
US20170102458A1 (en) * | 2014-06-09 | 2017-04-13 | Src, Inc. | Multiplatform GMTI Radar With Adaptive Clutter Suppression |
JP2016017929A (ja) * | 2014-07-11 | 2016-02-01 | 三菱電機株式会社 | レーダ装置 |
JP2016017928A (ja) * | 2014-07-11 | 2016-02-01 | 三菱電機株式会社 | レーダ装置 |
JP2016183953A (ja) * | 2015-03-25 | 2016-10-20 | パナソニック株式会社 | 物体検出装置および物体検出方法 |
WO2016199296A1 (ja) * | 2015-06-12 | 2016-12-15 | 三菱電機株式会社 | レーダ装置 |
JP2017125802A (ja) * | 2016-01-15 | 2017-07-20 | パナソニック株式会社 | 物体検出装置および物体検出方法 |
Also Published As
Publication number | Publication date |
---|---|
JP6587784B2 (ja) | 2019-10-09 |
WO2019030893A1 (ja) | 2019-02-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10175348B2 (en) | Use of range-rate measurements in a fusion tracking system via projections | |
Folster et al. | Data association and tracking for automotive radar networks | |
EP1610152A1 (en) | Tracking of a moving object for a self-defence system | |
CN111386476B (zh) | 确定车辆雷达系统中的对象运动和加速度矢量 | |
JPH07225275A (ja) | 車両用距離データ処理装置 | |
US20170146638A1 (en) | Radar device | |
JP6384018B2 (ja) | 車載用レーダ装置 | |
JP2018205174A (ja) | レーダ装置及びそのレーダ信号処理方法 | |
CA2092106C (en) | Three-dimensional maximum a posteriori (map) tracking | |
WO2021219023A1 (zh) | 定位方法、装置、电子设备及存储介质 | |
AU2010326314B2 (en) | System and method for discriminating targets at the water surface from targets below the water surface | |
JP2018205175A (ja) | レーダ装置及びそのレーダ信号処理方法 | |
Anderson et al. | Track association for over-the-horizon radar with a statistical ionospheric model | |
KR101303417B1 (ko) | 정보 처리 장치, 정보 처리 방법 및 기록 매체 | |
JP6587784B2 (ja) | レーダ信号処理装置及びレーダ信号処理方法 | |
JP2004251660A (ja) | 目標追尾装置 | |
JP6215543B2 (ja) | レーダ装置 | |
US20220317276A1 (en) | Radar signal processing device, radar system, and signal processing method | |
KR101844131B1 (ko) | 하모니 탐색 알고리즘 기반 tdoa 위치탐지 시스템에서의 레이더 펄스 모호성 해결 방법 및 장치 | |
Choi et al. | Doppler coherent focusing doa method for efficient radar map generation | |
JP5704983B2 (ja) | レーダ装置 | |
CN113670253B (zh) | 空间目标姿态反演方法、装置、计算设备及存储介质 | |
RU2196341C1 (ru) | Способ определения параметров движения маневрирующего объекта | |
Belyaev et al. | Object detection in an urban environment using 77GHz radar | |
JPWO2019220574A1 (ja) | 合成開口レーダ信号解析装置、合成開口レーダ信号解析方法および合成開口レーダ信号解析プログラム |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20190611 |
|
A871 | Explanation of circumstances concerning accelerated examination |
Free format text: JAPANESE INTERMEDIATE CODE: A871 Effective date: 20190611 |
|
A975 | Report on accelerated examination |
Free format text: JAPANESE INTERMEDIATE CODE: A971005 Effective date: 20190703 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20190813 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20190910 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 6587784 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |