JP6706056B2 - Signal processing device, radar device, underwater detection device, signal processing method, and program - Google Patents
Signal processing device, radar device, underwater detection device, signal processing method, and program Download PDFInfo
- Publication number
- JP6706056B2 JP6706056B2 JP2015245005A JP2015245005A JP6706056B2 JP 6706056 B2 JP6706056 B2 JP 6706056B2 JP 2015245005 A JP2015245005 A JP 2015245005A JP 2015245005 A JP2015245005 A JP 2015245005A JP 6706056 B2 JP6706056 B2 JP 6706056B2
- Authority
- JP
- Japan
- Prior art keywords
- signal processing
- spectrum
- processing device
- correlation matrix
- calculation unit
- 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.)
- Active
Links
- 238000012545 processing Methods 0.000 title claims description 279
- 238000001514 detection method Methods 0.000 title claims description 22
- 238000003672 processing method Methods 0.000 title description 6
- 238000001228 spectrum Methods 0.000 claims description 235
- 238000004364 calculation method Methods 0.000 claims description 148
- 239000011159 matrix material Substances 0.000 claims description 104
- 239000013598 vector Substances 0.000 claims description 103
- 238000000034 method Methods 0.000 claims description 55
- 238000012935 Averaging Methods 0.000 claims description 16
- 238000004148 unit process Methods 0.000 claims description 2
- 230000005540 biological transmission Effects 0.000 description 33
- 238000010586 diagram Methods 0.000 description 17
- 238000000354 decomposition reaction Methods 0.000 description 8
- 230000006870 function Effects 0.000 description 5
- 238000009434 installation Methods 0.000 description 5
- 230000003044 adaptive effect Effects 0.000 description 4
- 238000012217 deletion Methods 0.000 description 4
- 230000037430 deletion Effects 0.000 description 4
- 230000007261 regionalization Effects 0.000 description 3
- 238000004891 communication Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000008054 signal transmission Effects 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 241000251468 Actinopterygii Species 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000012447 hatching Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000017105 transposition Effects 0.000 description 1
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Description
本発明は、信号処理装置、レーダ装置、水中探知装置、信号処理方法、および、プログラムに関する。 The present invention relates to a signal processing device, a radar device, an underwater detection device, a signal processing method, and a program.
たとえば、レーダ装置は、レーダアンテナで受信した到来波を処理することで、到来波が到来する方向（以下、到来方向という）、および到来波の強度を算出する。到来方向をより精度よく算出するための方法として、超解像法が知られている。（たとえば、特許文献１〜５参照）。超解像法として、Ｃａｐｏｎ法、ＭＵＳＩＣ（MUltiple SIgnal Classification）法等が知られている。
For example, the radar device processes the incoming wave received by the radar antenna to calculate the direction in which the incoming wave arrives (hereinafter referred to as the incoming direction) and the intensity of the incoming wave. A super-resolution method is known as a method for more accurately calculating the arrival direction. (For example, refer to
特許文献１に記載の合成開口レーダ装置は、目標画像のクロスレンジ方向各成分を受信し、その後、固有値解析、および評価関数を得る処理を行う。
The synthetic aperture radar device described in
特許文献２に記載の方向検出装置は、アレーアンテナを有している。アレーアンテナのアンテナ素子毎に、受信回路が設けられている。これらの受信回路は、共分散行列作成部と接続されている。共分散行列作成部は、固有値分解処理部に接続されている。固有値分解処理部は、二次元ＭＵＳＩＣスペクトラム描画部に接続されている。
The direction detection device described in
特許文献３に記載の受信装置は、アレーアンテナを有している。アレーアンテナは、複数のセンサ（受信素子）を有している。また、センサ毎に、帯域制限ろ波器、中間周波数変換器、混合器、Ａ／Ｄ変換器、およびデジタルフィルタが設けられている。各デジタルフィルタからの出力は、１つの到来方位演算部へ与えられる。到来方位演算部は、ＭＵＳＩＣ法を用いて、入射信号の到来方位を算出する。
The receiving device described in
特許文献４，５に記載の受信装置は、アレーアンテナ装置を用いて、電波到来角を算出する。この算出過程において、ＭＵＳＩＣ法が用いられる。
The receiving devices described in
超解像法を用いて、到来波の到来方向を推定する場合、特許文献２〜５に記載されているように、アレーアンテナを用いることが前提となっている。アレーアンテナは、静止した状態で、多数の方位からの到来波を受信することができる。
When estimating the arrival direction of an incoming wave using the super-resolution method, it is premised that an array antenna is used as described in
一方で、たとえば、レーダ装置のアンテナとして、鉛直方向に延びる回転軸線回りを自転するアンテナ（以下、回転アンテナともいう。）が用いられることがある。回転アンテナは、回転軸線回りを回転しながら、パルス状電波であるビームの送信と、この送信によって得られる受信信号の受信と、を繰り返す。回転アンテナから送信されるビームの幅（ビーム幅）が狭いほど、ビームの指向性を高くできる。そして、ビーム幅が狭いほど、物標検出の分解能を高くできる。 On the other hand, for example, an antenna that rotates about a rotation axis extending in the vertical direction (hereinafter, also referred to as a rotating antenna) may be used as an antenna of a radar device. The rotating antenna repeats the transmission of the beam of the pulsed radio wave and the reception of the reception signal obtained by this transmission while rotating around the rotation axis. The narrower the width of the beam transmitted from the rotating antenna (beam width), the higher the directivity of the beam. Then, the narrower the beam width, the higher the resolution of the target detection.
そして、回転アンテナの全長を長くすることで、上記ビーム幅を、より狭くでき、その結果、レーダ装置における物標検出の分解能は、より高くなる。しかしながら、全長の大きな回転アンテナの製造コストは、高い。また、全長の大きな回転アンテナは、小型船舶等、設置スペースに制限がある場所への設置に適さない。このため、回転アンテナの全長に拘わらず、高い分解能を実現できる構成が求められている。 By increasing the total length of the rotary antenna, the beam width can be made narrower, and as a result, the resolution of target detection in the radar device becomes higher. However, the manufacturing cost of the rotating antenna having a large total length is high. Further, the rotating antenna having a large total length is not suitable for installation in a place where the installation space is limited, such as a small boat. Therefore, there is a demand for a configuration that can realize high resolution regardless of the total length of the rotating antenna.
同様の課題は、ソナー装置等、到来波を受信する他の装置においても存在する。 Similar problems also exist in other devices that receive incoming waves, such as sonar devices.
そこで、本発明は、機械的に変位する受信装置を用いて得られた受信信号の処理に際して、より高い分解能を実現できる、信号処理装置、レーダ装置、水中探知装置、信号処理方法、および、プログラムを提供することを目的とする。 Therefore, the present invention provides a signal processing device, a radar device, an underwater detection device, a signal processing method, and a program capable of realizing higher resolution when processing a reception signal obtained using a mechanically displaced receiving device. The purpose is to provide.
（１）上記課題を解決するために、この発明のある局面に係わる信号処理装置は、相関行列算出部と、スペクトラム算出部と、を備えている。前記相関行列算出部は、所定の副走査方向に沿って変位し且つ前記副走査方向と交差する主走査方向からの信号を順次受信する受信装置を用いて得られた受信信号に基づいて、前記副走査方向に沿って得られる複数の前記受信信号の相関行列を算出する。前記スペクトラム算出部は、前記相関行列に基づいて超解像処理を行うことで、前記副走査方向に沿って得られる前記受信信号の信号強度のスペクトラムを算出する。 (1) In order to solve the above problems, a signal processing device according to an aspect of the present invention includes a correlation matrix calculation unit and a spectrum calculation unit. The correlation matrix calculation unit is based on a reception signal obtained by using a reception device that sequentially receives signals from a main scanning direction that is displaced along a predetermined sub-scanning direction and intersects the sub-scanning direction, A correlation matrix of the plurality of received signals obtained along the sub-scanning direction is calculated. The spectrum calculation unit calculates the spectrum of the signal intensity of the received signal obtained along the sub-scanning direction by performing super-resolution processing based on the correlation matrix.
（２）ある場合において、前記信号処理装置は、前記受信装置のステアリングベクトルを算出する、ステアリングベクトル算出部をさらに備える。前記スペクトラム算出部は、前記相関行列、および前記ステアリングベクトルに基づいて、前記スペクトラムを算出する。 (2) In some cases, the signal processing device further includes a steering vector calculation unit that calculates a steering vector of the reception device. The spectrum calculation unit calculates the spectrum based on the correlation matrix and the steering vector.
（３）ある場合において、前記信号処理装置は、前記相関行列の固有ベクトルを算出する固有ベクトル算出部をさらに備える。前記スペクトラム算出部は、前記固有ベクトル、および、前記ステアリングベクトルに基づいて、前記スペクトラムを計算する。 (3) In some cases, the signal processing device further includes an eigenvector calculation unit that calculates an eigenvector of the correlation matrix. The spectrum calculation unit calculates the spectrum based on the eigenvector and the steering vector.
（４）ある場合において、前記ステアリングベクトル算出部は、前記受信装置のアンテナパターンに基づいて前記ステアリングベクトルを算出する。 (4) In some cases, the steering vector calculation unit calculates the steering vector based on the antenna pattern of the reception device.
（５）ある場合において、前記アンテナパターンは、前記受信装置のメインローブのパターンに相当する。 (5) In some cases, the antenna pattern corresponds to a main lobe pattern of the receiving device.
（６）ある場合において、前記信号処理装置は、前記副走査方向に沿って複数の領域を設定する領域設定部をさらに備える。前記相関行列算出部は、前記領域毎に前記相関行列を算出する。前記スペクトラム算出部は、前記領域毎に、前記スペクトラムを算出する。 (6) In one case, the signal processing device further includes a region setting unit that sets a plurality of regions along the sub-scanning direction. The correlation matrix calculation unit calculates the correlation matrix for each area. The spectrum calculation unit calculates the spectrum for each area.
（７）ある場合において、前記信号処理装置は、前記相関行列の固有ベクトルを算出する固有ベクトル算出部と、前記受信装置のステアリングベクトルを算出する、ステアリングベクトル算出部と、をさらに備える。前記固有ベクトル算出部および前記ステアリングベクトル算出部は、それぞれ、前記領域毎に前記固有ベクトルおよび前記ステアリングベクトルを算出する。前記スペクトラム算出部は、前記領域毎に、前記固有ベクトルと前記ステアリングベクトルとを用いて前記スペクトラムを算出する。 (7) In one case, the signal processing device further includes an eigenvector calculation unit that calculates an eigenvector of the correlation matrix, and a steering vector calculation unit that calculates a steering vector of the reception device. The eigenvector calculation unit and the steering vector calculation unit respectively calculate the eigenvector and the steering vector for each region. The spectrum calculation unit calculates the spectrum for each of the regions using the eigenvector and the steering vector.
（８）ある場合において、前記信号処理装置は、相関行列平均化処理部をさらに備える。前記相関行列平均化処理部は、前記領域毎に前記相関行列を平均化することで平均相関行列を算出する。前記固有ベクトル算出部は、前記平均相関行列を固有値分解することで、前記固有ベクトルを算出する。前記ステアリングベクトル算出部は、各前記領域において平均化された平均ステアリングベクトルを、前記ステアリングベクトルとして算出する。 (8) In some cases, the signal processing device further includes a correlation matrix averaging processing unit. The correlation matrix averaging processing unit calculates an average correlation matrix by averaging the correlation matrix for each of the regions. The eigenvector calculation unit calculates the eigenvector by performing eigenvalue decomposition on the average correlation matrix. The steering vector calculation unit calculates an average steering vector averaged in each of the regions as the steering vector.
（９）ある場合において、前記スペクトラム算出部は、前記領域毎に、予め定められた数のピークが現れる手法を選択し、この手法を用いて前記スペクトラムを算出する。 (9) In some cases, the spectrum calculation unit selects a method in which a predetermined number of peaks appear in each of the regions, and calculates the spectrum using this method.
（１０）ある場合において、前記手法は、前記予め定められた数が１の場合には超解像法であり、前記予め定められた数が複数の場合にはＭＵＳＩＣ法である。 (10) In some cases, the method is a super-resolution method when the predetermined number is 1, and a MUSIC method when the predetermined number is plural.
（１１）ある場合において、前記信号処理装置は、フィルタ処理部を更に備える。前記フィルタ処理部は、前記領域毎に算出された複数の前記スペクトラムのなかから、一部の前記スペクトラムを選択する。 (11) In one case, the signal processing device further includes a filter processing unit. The filter processing unit selects a part of the spectrum from the plurality of spectra calculated for each area.
（１２）ある場合において、前記信号処理装置は、ビームパターン加工部をさらに備える。前記ビームパターン加工部は、前記スペクトラムの包絡線の形状を加工する。 (12) In some cases, the signal processing device further includes a beam pattern processing unit. The beam pattern processing unit processes the shape of the envelope of the spectrum.
（１３）ある場合において、前記ビームパターン加工部は、前記包絡線の形状を、予め定められたパターン形状に置き換え、且つ、当該パターン形状の前記信号強度のピークを、前記スペクトラムの前記信号強度のピークと同じ値に設定する。 (13) In some cases, the beam pattern processing unit replaces the shape of the envelope with a predetermined pattern shape, and sets the peak of the signal strength of the pattern shape to the signal strength of the spectrum. Set to the same value as the peak.
（１４）ある場合において、前記パターン形状における前記ピークを中心とする所定範囲の幅は、前記スペクトラムの包絡線における前記ピークを中心とする所定範囲の幅よりも大きい。 (14) In some cases, the width of the predetermined range centered on the peak in the pattern shape is larger than the width of the predetermined range centered on the peak in the envelope of the spectrum.
（１５）上記課題を解決するために、この発明のある局面に係わるレーダ装置は、送受信装置と、前述の信号処理装置と、を備える。前記送受信装置は、所定の方位方向を副走査方向として回転し且つ前記方位方向と交差する距離方向に沿って電磁波を順次送受信するアンテナ部を有している。前記信号処理装置の前記相関行列算出部は、前記アンテナ部を用いて得られた受信信号に基づいて、前記相関行列を算出する。 (15) In order to solve the above problems, a radar device according to an aspect of the present invention includes a transmission/reception device and the signal processing device described above. The transmission/reception device has an antenna unit that rotates in a predetermined azimuth direction as a sub-scanning direction and sequentially transmits/receives electromagnetic waves along a distance direction intersecting the azimuth direction. The correlation matrix calculation unit of the signal processing device calculates the correlation matrix based on a received signal obtained using the antenna unit.
（１６）上記課題を解決するために、この発明のある局面に係わる水中探知装置は、送受信装置と、前述の信号処理装置と、を備える。前記送受信装置は、所定の方位方向を副走査方向として回転し且つ前記方位方向と交差する距離方向に沿って超音波を順次送受信する超音波振動子を有している。前記信号処理装置の前記相関行列算出部は、前記超音波振動子を用いて得られた受信信号に基づいて、前記相関行列を算出する。 (16) In order to solve the above problems, an underwater detection apparatus according to an aspect of the present invention includes a transmission/reception apparatus and the above-described signal processing apparatus. The transmission/reception device includes an ultrasonic transducer that rotates in a predetermined azimuth direction as a sub-scanning direction and sequentially transmits/receives ultrasonic waves along a distance direction intersecting the azimuth direction. The correlation matrix calculation unit of the signal processing device calculates the correlation matrix based on a received signal obtained using the ultrasonic transducer.
（１７）上記課題を解決するために、この発明のある局面に係わる信号処理方法は、相関行列算出ステップと、スペクトラム算出ステップと、を含んでいる。前記相関行列算出ステップは、所定の副走査方向に沿って変位し且つ前記副走査方向と交差する主走査方向からの信号を順次受信する受信装置を用いて得られた受信信号に基づいて、前記副走査方向に沿って得られる複数の前記受信信号の相関行列を算出する。前記スペクトラム算出ステップは、前記相関行列に基づいて超解像処理を行うことで、前記副走査方向に沿って得られる前記受信信号の信号強度のスペクトラムを算出する。 (17) In order to solve the above problems, a signal processing method according to an aspect of the present invention includes a correlation matrix calculating step and a spectrum calculating step. The correlation matrix calculation step is based on a reception signal obtained by using a reception device that is displaced along a predetermined sub-scanning direction and sequentially receives signals from a main scanning direction intersecting the sub-scanning direction, A correlation matrix of the plurality of received signals obtained along the sub-scanning direction is calculated. In the spectrum calculation step, a super-resolution process is performed based on the correlation matrix to calculate a spectrum of signal intensity of the received signal obtained along the sub-scanning direction.
（１８）上記課題を解決するために、この発明のある局面に係わるプログラムは、相関行列算出ステップと、スペクトラム算出ステップと、をコンピュータに実行させるプログラムである。前記相関行列算出ステップは、所定の副走査方向に沿って変位し且つ前記副走査方向と交差する主走査方向からの信号を順次受信する受信装置を用いて得られた受信信号に基づいて、前記副走査方向に沿って得られる複数の前記受信信号の相関行列を算出する。前記スペクトラム算出ステップは、前記相関行列に基づいて超解像処理を行うことで、前記副走査方向に沿って得られる前記受信信号の信号強度のスペクトラムを算出する。 (18) In order to solve the above problems, a program according to an aspect of the present invention is a program for causing a computer to execute a correlation matrix calculating step and a spectrum calculating step. The correlation matrix calculation step is based on a reception signal obtained by using a reception device that is displaced along a predetermined sub-scanning direction and sequentially receives signals from a main scanning direction intersecting the sub-scanning direction, A correlation matrix of the plurality of received signals obtained along the sub-scanning direction is calculated. In the spectrum calculation step, a super-resolution process is performed based on the correlation matrix to calculate a spectrum of signal intensity of the received signal obtained along the sub-scanning direction.
本発明によると、機械的に変位する受信装置を用いて得られた受信信号の処理に際して、より高い分解能を実現できる。 According to the present invention, a higher resolution can be realized when processing a received signal obtained using a mechanically displaced receiving device.
以下、本発明を実施するための形態について、図面を参照しつつ説明する。本発明は、分解能を高めるための信号処理装置（到来波方向推定装置）として広く適用することができる。なお、以下では、図中同一または相当部分には、同一符号を付してその説明は繰り返さない。 Hereinafter, embodiments for carrying out the present invention will be described with reference to the drawings. INDUSTRIAL APPLICABILITY The present invention can be widely applied as a signal processing device (arrival wave direction estimation device) for increasing resolution. In the following, the same or corresponding parts in the drawings will be denoted by the same reference numerals and description thereof will not be repeated.
［第１実施形態］
［レーダ装置の概略構成］
図１は、本発明の第１実施形態にかかる信号処理装置３を備えるレーダ装置１の概略構成を示すブロック図である。レーダ装置１は、漁船などの船舶に搭載される。以下、レーダ装置１が搭載される船舶を、自船という。レーダ装置１は、自船の周囲の物標を探知する。
[First Embodiment]
[Schematic configuration of radar device]
FIG. 1 is a block diagram showing a schematic configuration of a
レーダ装置１は、送受信装置２と、信号処理装置３と、操作・表示装置４と、を有している。
The
送受信装置２は、指向性を有するパルス状の電波を送信信号として送信する。また、送受信装置２は、この送信信号に対するエコー信号を受信する。送受信装置２は、後述するように、単一のアンテナを有している。送受信装置２は、アンテナで受信された受信信号Ｓ１１をデジタルの受信データｘ（φｉ，Ｒｊ）に変換し、この受信データｘ（φｉ，Ｒｊ）を信号処理装置３へ出力する。
The transmission/
信号処理装置３は、受信データｘ（φｉ，Ｒｊ）（受信信号Ｓ１１）の方位φと信号強度との関係を、角度スペクトラムＰＳＳＡ（φ）として算出するように構成されている。このような構成により、レーダ装置１は、当該レーダ装置１の周囲に存在する物標を探知する。信号処理装置３は、受信データｘ（φｉ，Ｒｊ）の位相と強度との関係を示す画像データＧを、操作・表示装置４へ出力する。
The
操作・表示装置４は、操作部５と、表示部６とを有している。操作部５は、レーダ装置１のオペレータによって操作され得る。操作部５は、種々の入力キー等を備えており、電磁波の送受信、および映像表示等に必要な、種々の設定値等を入力できるように構成されている。操作部５の設定に応じて、角度スペクトラム算出に必要な定数などが設定される構成であってもよい。表示部６は、信号処理装置３から出力された画像データＧに応じた映像を表示する。
The operation/
［レーダ装置の詳細な構成］
図２は、レーダ装置１の詳細な構成を示すブロック図である。図３は、レーダ装置１が自船１０に設置された状態を示す、模式的な平面図である。図３は、自船１０の周囲に１つの物標１１が存在している状態を例示している。
[Detailed configuration of radar device]
FIG. 2 is a block diagram showing a detailed configuration of the
図２および図３を参照して、本実施形態では、平面視において自船１０を中心とする時計回り方向を、方位方向Ｄ１という。また、平面視において、自船１０を中心とする径方向を、距離方向Ｄ２という。
With reference to FIGS. 2 and 3, in the present embodiment, the clockwise direction centered on the
送受信装置２は、送信部２１と、サーキュレータ２２と、回転アンテナ部（受信装置）２３と、受信部２４と、を有している。
The transmitter/
送信部２１は、マグネトロン等を有している。より具体的には、送信部２１は、マイクロ波を発振する電子素子としてたとえばＤ／Ａコンバータ、周波数変換器、電力増幅器を有している。送信部２１は、定期的にパルス状のレーダ送信信号を生成し、このレーダ送信信号をサーキュレータ２２へ出力する。
The
送信部２１からレーダ送信信号が出力されるタイミングを示すデータは、送信部２１から信号処理装置３へ出力される。これにより、信号処理装置３は、回転アンテナ部２３から出力された電磁波が反射して生じたエコー信号が回転アンテナ部２３で受信されるまでの時間を、算出できる。その結果、信号処理装置３は、物標１１と自船１０との間の距離を算出することができる。
Data indicating the timing at which the radar transmission signal is output from the
サーキュレータ２２は、送信部２１から出力されたレーダ送信信号を、回転アンテナ部２３へ出力するように構成されている。また、サーキュレータ２２は、回転アンテナ部２３で受信された受信信号Ｓ１１を、受信部２４へ出力する。
The
回転アンテナ部２３は、たとえば、自船１０のマスト（図示せず）に設置されている。回転アンテナ部２３は、たとえば、スロットアレイアンテナである。
The
回転アンテナ部２３は、全体として、水平方向に細長い形状（長尺形状）に形成されており、自船１０の上下方向に延びる回転軸線回りを自転するように構成されている。回転アンテナ部２３は、当該回転アンテナ部２３の方位を示す方位信号を、信号処理装置３に出力する。
The
回転アンテナ部２３は、方位方向Ｄ１に変位し、且つ、方位方向Ｄ１と交差する距離方向Ｄ２において電磁波（信号）を順次送受信するように構成されている。なお、方位方向Ｄ１は、本発明の「副走査方向」の一例であり、距離方向Ｄ２は、本発明の「主走査方向」の一例である。
The
前述したように、回転アンテナ部２３は、スロットアレイアンテナなどであり、高い指向性を有する送信信号を送信する。また、回転アンテナ部２３は、受信信号Ｓ１１を受信するように構成されている。受信信号Ｓ１１は、送信信号の反射信号としてのエコー信号と、ノイズ信号と、を含んでいる。
As described above, the rotating
図３は、回転アンテナ部２３におけるビームパターンＰ１の一例を示している。ビームパターンＰ１は、メインローブＰ１１と、サイドローブＰ１２と、を有している。回転アンテナ部２３は、実質的にメインローブＰ１１で信号を受信し、この信号を受信信号Ｓ１１として出力する。
FIG. 3 shows an example of the beam pattern P1 in the
回転アンテナ部２３は、所定の方位（たとえば、方位方向Ｄ１における０．１度）毎に、電磁波の送信と受信とを繰り返す。本実施形態では、レーダ装置１が電磁波を送信してから次の電磁波を送信するまでの動作を「スイープ」という。送受信装置２は、１スイープ毎に、受信信号Ｓ１１を受信部２４へ出力する。本実施形態では、自船１０を中心とする方位方向Ｄ１に沿う、電磁波（パルス信号）の送信角度間隔Δφは、０．１度である。この角度間隔Δφは、回転アンテナ部２３の回転速度と、レーダ送信信号のパルス信号の繰返周波数に依存する。
The
受信部２４は、回転アンテナ部２３から受信信号Ｓ１１を取り込み、増幅した後に中間周波数に変換する処理などを行うことで、受信信号Ｓ１１を受信データｘ（φｉ，Ｒｊ）に変換する。本実施形態では、受信部２４は、１つのみ設けられている。この点、複数の受信素子が備えられ且つ受信素子毎に受信部が設けられるアレイアンテナ装置とは、異なる構成である。本実施形態では、受信部２４の数は、回転アンテナ部２３の特性に拘らず、１つである。
The
受信部２４は、ＬＮＡ（ローノイズアンプ）２５と、ミキサ２６と、とＡＭＰ（リニアアンプ）２７と、Ａ／Ｄコンバータ２８と、局部発振器２９とを有している。
The receiving
ＬＮＡ２５は、レーダ受信信号Ｓ１１を増幅するように構成されている。ＬＮＡ２５で増幅された受信信号Ｓ１１は、ミキサ２６へ出力される。ミキサ２６は、局部発振器２９から出力されたローカル信号と、受信信号Ｓ１１とをミキシングすることで、レーダ受信信号Ｓ１１の周波数を、中間周波数に変換する。周波数変換された受信信号Ｓ１１’は、ＡＭＰ２７へ出力される。
The
ＡＭＰ２７は、中間周波数に変換されたレーダ受信信号Ｓ１１’を対数増幅し、包絡線信号を生成する。Ａ／Ｄコンバータ２８は、ＡＭＰ２７が出力したレーダ受信信号Ｓ１１’’をＩＱ検波等でサンプリングする。これにより、Ａ／Ｄコンバータ２８は、レーダ受信信号Ｓ１１’’を、複素デジタル信号に変換する。即ち、Ａ／Ｄコンバータ２８は、レーダ受信信号Ｓ１１’’を、デジタルの受信データｘ（φｉ，Ｒｊ）に変換する。受信データｘ（φｉ，Ｒｊ）は、信号強度（振幅値）を示す。また、φｉ，Ｒｊは、それぞれ、方位方向Ｄ１と、距離方向Ｄ２に沿った位置（経過時間）とを示す。
The
図４は、受信データｘ（φｉ，Ｒｊ）について説明するための模式的な平面図であり、各受信データｘ（φｉ，Ｒｊ）の方位および経過時間について、自船１０を中心に視覚的に示す図である。経過時間は、たとえば、回転アンテナ部２３が電磁波を送信してから受信信号Ｓ１１を受信するまでの時間をいう。
Figure 4 is a schematic plan view for describing received data x (φ i, R j) , the orientation and the elapsed time of each received data x (φ i, R j) , around the
図４を参照して、図中の白抜きの丸字は、各受信データｘ（φｉ，Ｒｊ）の到来方位および経過時間を示している。同一方位上に存在する各受信データｘ（φｉ，Ｒｊ）は、自船１０から延びる直線上に並んでいる。即ち、同一スイープにおける各受信データｘ（φｉ，Ｒｊ）は、自船１０から時系列に沿って直線状に並んでいる。なお、自船１０に近い位置ほど、回転アンテナ部２３に対する受信データｘ（φｉ，Ｒｊ）の受信開始からの経過時間が短いことを示している。
Referring to FIG. 4, white circles in the figure indicate the arrival azimuth and the elapsed time of each reception data x(φ i , R j ). The respective reception data x(φ i , R j ) existing in the same direction are arranged on a straight line extending from the
各受信データｘ（φｉ，Ｒｊ）の経過時間に対して、短い経過時間から順に経過時間番号ｊ（ｊ＝１，２，３，…，Ｂ−１，Ｂ）が付されている。Ｂの値は、受信部２４における受信データｘ（φｉ，Ｒｊ）のサンプリング周波数等に依存して決まる。各スイープにおいて、受信部２４によるレーダ受信信号Ｓ１１のサンプリングは、一定の時間間隔で行われており、各受信データｘ（φｉ，Ｒｊ）の時間間隔は、一定である。
Elapsed time numbers j (j=1, 2, 3,..., B-1, B) are sequentially added to the elapsed time of each received data x(φ i , R j ) from the shortest elapsed time. The value of B is determined depending on the sampling frequency or the like of the reception data x(φ i , R j ) in the
受信データｘ（φｉ，Ｒｊ）には、自船１０回りの単位方位毎に、時計周りの順で方位番号ｉ（ｉ＝１，２，３，…Ｃ−１，Ｃ）が付されている。上記受信データｘ（φｉ，Ｒｊ）は、図５に示すように、行列として示すことができる。図５は、受信データｘ（φｉ，Ｒｊ）を行列状に配置した状態を示す図である。図５において、横方向は、方位方向Ｄ１を示し、縦方向は距離方向Ｄ２を示している。
The received data x (φ i , R j ) is assigned a bearing number i (i=1, 2, 3,... C-1, C) in clockwise order for each unit bearing around the
図２、図３および図５を参照して、送受信装置２は、一回のスイープによって、方位方向Ｄ１が同じである複数の受信信号Ｓ１１を受信する。送受信装置２は、複数回のスイープによって、方位方向Ｄ１の位置が異なる複数の受信信号Ｓ１１を取得する。
With reference to FIGS. 2, 3 and 5, transmitting/receiving
信号処理装置３は、送受信装置２から出力された受信データｘ（φｉ，Ｒｊ）を用いて、超解像処理を行う。即ち、信号処理装置３は、回転アンテナ部２３の解像度を超える解像度で、物標（物標１１など）と物標以外とを識別するように構成されている。
The
［信号処理装置の構成］
信号処理装置３は、方位領域設定部３１と、相関行列算出部３２と、空間平均処理部３３と、固有ベクトル算出部３４と、角度スペクトラム算出部３５と、アンテナパターン記憶部３６と、平均ステアリングベクトル算出部３７と、表示用処理部３８と、を有している。
[Configuration of signal processing device]
The
本実施形態では、信号処理装置３は、自船１０を中心とする方位方向Ｄ１において区切られた所定の領域Ｆを単位として、超解像処理（高分解能処理）を行う。より具体的には、信号処理装置３は、所定領域Ｆにおいて、経過時間番号ｊが同じである複数の受信データｘ（φｉ，Ｒｊ）を用いて、超解像処理を行う。
In the present embodiment, the
次に、信号処理装置３における、受信データｘ（φｉ，Ｒｊ）の領域化の一例を説明する。図６は、受信データｘ（φｉ，Ｒｊ）の領域化の一例を説明するための図である。図２および図６を参照して、信号処理装置３は、方位方向Ｄ１に関して、受信データｘ（φｉ，Ｒｊ）を複数の領域Ｆ（Ｆ＝１，２，…Ｕ）に分けた状態で、高分解能処理を行う。各領域Ｆは、複数の連続するスイープによって得られた受信データｘ（φｉ，Ｒｊ）を含んでいる。
Next, an example of regionalization of the reception data x(φ i , R j ) in the
図６では、領域化の例が示されている。この例では、隣り合う領域Ｆ（たとえば、Ｆ１，Ｆ２）が共通の受信データｘ（φｉ，Ｒｊ）を有しておらず、且つ、隣り合う領域Ｆ（たとえば、Ｆ１，Ｆ２）の間には何れの領域Ｆにも属していない受信データｘ（φｉ，Ｒｊ）が存在していない。本実施形態では、方位方向Ｄ１に関して一定の間隔毎に、領域Ｆが分けられている。本実施形態では、信号処理装置３は、各領域Ｆ１，Ｆ２，・・・，ＦＵに基づいて、処理を行う。
FIG. 6 shows an example of regionalization. In this example, the adjacent regions F (for example, F1 and F2) do not have the common reception data x (φ i , R j ) and between the adjacent regions F (for example, F1 and F2). Does not include received data x(φ i , R j ) that does not belong to any area F. In this embodiment, the region F is divided at regular intervals in the azimuth direction D1. In the present embodiment, the
［方位領域設定部の構成］
方位領域設定部３１は、方位方向Ｄ１において複数の領域Ｆ（たとえばＦ１，Ｆ２，Ｆ３）を設定するように構成されている。具体的には、方位領域設定部３１は、経過時間番号ｊが同じである受信データｘｉを、受信部２４から読み出す。
[Configuration of azimuth region setting section]
The azimuth
図７は、受信データｘｉの一例を示すグラフである。図７では、受信データｘｉで特定される信号強度が示されている。図７では、一例として、方位番号ｉ＝１〜３ｎ（ｎは自然数）までの各受信データ受信データｘｉが示されている。なお、以下では、１つの経過時間番号ｊ（距離）における受信データｘｉについて説明し、経過時間番号ｊの記載は省略する。受信データｘ１，ｘ２，・・・ｘ３ｎにおいては、受信データｘ３ｎ／２がピークの振幅を示している。 FIG. 7 is a graph showing an example of the reception data x i . In FIG. 7, the signal strength specified by the reception data x i is shown. In FIG. 7, as an example, the respective received data received data x i up to azimuth numbers i=1 to 3 n (n is a natural number) are shown. In the following, the received data x i at one elapsed time number j (distance) will be described, and the description of the elapsed time number j will be omitted. In the reception data x 1 , x 2 ,... X 3n , the reception data x 3n/2 shows the peak amplitude.
図２および図７を参照して、本実施形態では、方位領域設定部３１は、たとえば、方位番号ｉ＝１〜３ｎまでの各受信データｘｉを受信部２４から読み出す。方位領域設定部３１は、これらの受信データｘｉを、複数の領域Ｆ（本実施形態では、Ｆ１，Ｆ２，Ｆ３）に分ける。
With reference to FIG. 2 and FIG. 7, in the present embodiment, the azimuth
領域Ｆ１は、受信データｘ１〜ｘｎを含んでいる。領域１のベクトルデータｘ｛１｝は、下記式で表される。
ｘ｛１｝＝［ｘ１，ｘ２，ｘ３，・・・，ｘｎ］Ｔ
なお、Ｔは、転置を表している。
The area F1 includes received data x 1 to x n . The vector data x {1} of the
x {1} = [x 1 , x 2 , x 3 ,..., X n ] T
In addition, T represents transposition.
領域Ｆ２は、受信データｘｎ＋１〜ｘ２ｎを含んでいる。領域２のベクトルデータｘ｛２｝は、下記式で表される。
ｘ｛２｝＝［ｘｎ＋１，ｘｎ＋２，ｘｎ＋３，・・・，ｘ２ｎ］Ｔ
The area F2 includes received data x n+1 to x 2n . The vector data x {2} of the
x {2} =[ xn+1 , xn+2 , xn+3 ,..., x2n ] T
領域Ｆ３は、受信データｘ２ｎ＋１〜ｘ３ｎを含んでいる。領域３のベクトルデータｘ｛３｝は、下記式で表される。
ｘ｛３｝＝［ｘ２ｎ＋１，ｘ２ｎ＋２，ｘ２ｎ＋３，・・・，ｘ３ｎ］Ｔ
The area F3 includes received data x 2n+1 to x 3n . The vector data x {3} of the
x {3} = [x 2n+1 , x 2n+2 , x 2n+3 ,..., x 3n ] T
方位領域設定部３１は、ベクトルデータｘ（ｘ｛１｝，ｘ｛２｝，ｘ｛３｝）を、相関行列算出部３２へ出力する。
The azimuth
［相関行列算出部の構成］
相関行列算出部３２は、受信信号Ｓ１１に基づいて、方位方向Ｄ１に沿って複数の受信信号Ｓ１１の相関行列Ｒｘｘを算出する。相関行列Ｒｘｘは、１つの回転アンテナ部２３による複数のスイープで得られた受信データｘｉに対して、適応ビームフォーミング法を適用するために算出される。この場合、複数回のスイープで得られた受信データｘｉは、複数の仮想アンテナによって一度に得られた受信データｘｉとして捉えることができる。相関行列算出部３２は、領域Ｆ１，Ｆ２，Ｆ３毎に、相関行列（分散行列）Ｒｘｘ（Ｒ｛１｝ ｘｘ，Ｒ｛２｝ ｘｘ，Ｒ｛３｝ ｘｘ）を算出する。
[Configuration of correlation matrix calculation unit]
The correlation
ベクトルデータｘ｛１｝の相関行列Ｒ｛１｝ ｘｘは、下記（１）式で表される。
なお、Ｅ［・］は、期待値（アンサンブル平均）を求める操作を表し、Ｈは複素共役転置を表す。
The correlation matrix R {1} xx of the vector data x {1} is represented by the following equation (1).
Note that E[·] represents an operation for obtaining an expected value (average ensemble), and H represents a complex conjugate transpose.
ベクトルデータｘ｛２｝の相関行列Ｒ｛２｝ ｘｘは、下記（２）式で表される。
ベクトルデータｘ｛３｝の相関行列Ｒ｛３｝ ｘｘは、下記式（３）式で表される。
相関行列算出部３２は、各相関行列Ｒ｛１｝ ｘｘ，Ｒ｛２｝ ｘｘ，Ｒ｛３｝ ｘｘのデータを、空間平均処理部３３へ出力する。なお、本実施形態では、ｘ・ｘ＊を算出する形態を例に説明している。しかしながら、この通りでなくてもよい。たとえば、ｘ・ｘ＊に代えてｘ・ｘを用いてもよい。
The correlation
［空間平均処理部の構成］
空間平均処理部３３は、本発明の「相関行列平均化処理部」の一例である。空間平均処理部３３は、相関行列Ｒ｛１｝ ｘｘ，Ｒ｛２｝ ｘｘ，Ｒ｛３｝ ｘｘを個別に（領域Ｆ１，Ｆ２，…毎に）平均化することで、空間平均相関行列Ｒ(＿)ｘｘ（Ｒ(＿)｛１｝ ｘｘ，Ｒ(＿)｛２｝ ｘｘ，Ｒ(＿)｛３｝ ｘｘ）を算出する。空間平均処理部３３は、適応ビームフォーミング法におけるサブアレイの手法を用いることで、空間平均相関行列Ｒ(＿)ｘｘを算出する。
[Configuration of spatial averaging processing unit]
The spatial
空間平均処理部３３は、算出した空間平均相関行列Ｒ(＿)ｘｘのデータを、固有ベクトル算出部３４へ出力する。
The spatial
［固有ベクトル算出部の構成］
固有ベクトル算出部３４は、領域Ｆ１，Ｆ２，Ｆ３毎に、固有値分解を行う。具体的には、固有ベクトル算出部３４は、空間平均相関行列Ｒ(＿)｛１｝ ｘｘ，Ｒ(＿)｛２｝ ｘｘ，Ｒ(＿)｛３｝ ｘｘのそれぞれについて、固有値分解を行う。固有ベクトル算出部３４は、たとえばべき乗法、Jacobi法、ＱＲ法などを用いて、固有値分解を行う。
[Configuration of eigenvector calculation unit]
The
固有ベクトル算出部３４は、空間平均相関行列Ｒ(＿)｛１｝ ｘｘの固有ベクトルｑ｛１｝ １，ｑ｛１｝ ２，…，ｑ｛１｝ Ｍと、固有値λ｛１｝ １，λ｛１｝ ２，…，λ｛１｝ Ｍとを算出する。なお、Ｍは変数である。
The
上記と同様にして、固有ベクトル算出部３４は、空間平均相関行列Ｒ(＿)｛２｝ ｘｘの固有ベクトルｑ｛２｝ １，…，ｑ｛２｝ Ｍと、固有値λ｛２｝ １，…，λ｛２｝ Ｍとを算出する。また、固有ベクトル算出部３４は、空間平均相関行列Ｒ(＿)｛３｝ ｘｘの固有ベクトルｑ｛３｝ １，…，ｑ｛３｝ Ｍと、固有値λ｛３｝ １，…，λ｛３｝ Ｍとを算出する。固有ベクトル算出部３４は、算出した固有ベクトルｑ（ｑ｛１｝ １〜ｑ｛１｝ Ｍ，ｑ｛２｝ １〜ｑ｛２｝ Ｍ，ｑ｛３｝ １〜ｑ｛３｝ Ｍ）のデータを、角度スペクトラム算出部３５へ出力する。
In the same manner as described above, the
角度スペクトラム算出部３５は、領域Ｆ１，Ｆ２，Ｆ３毎に、角度スペクトラムＰＳＳＡ（φ）（ＰＳＳＡ ｛１｝（φ），ＰＳＳＡ ｛２｝（φ），ＰＳＳＡ ｛３｝（φ））を算出する。角度スペクトラムＰＳＳＡ（φ）は、方位方向Ｄ１の位置と受信データｘｉ（受信信号Ｓ１１）の信号強度（振幅）との関係を示す。角度スペクトラム算出部３５は、相関行列Ｒｘｘと、平均ステアリングベクトル算出部３７から与えられる平均ステアリングベクトルＡ(＿)（φ）（Ａ(＿)｛１｝（φ），Ａ(＿)｛２｝（φ），Ａ(＿)｛３｝（φ）]と、に基づいて超解像処理を行うことで、角度スペクトラムＰＳＳＡ（φ）を算出する。
The angle
より具体的には、角度スペクトラム算出部３５は、固有ベクトル算出部３４から与えられる固有ベクトルｑと、平均ステアリングベクトル算出部３７から与えられる平均ステアリングベクトルＡ(＿)（φ）と、を用いて、角度スペクトラムＰＳＳＡ（φ）を算出する。本実施形態では、平均ステアリングベクトル算出部３７は、アンテナパターン記憶部３６に格納されているアンテナパターンＡＰを基に、平均ステアリングベクトルＡ(＿)（φ）を算出する。
More specifically, the angle
図８は、回転アンテナ部２３のアンテナパターンＡＰを示すグラフである。図２、図７および図８を参照して、本実施形態では、アンテナパターンＡＰは、回転アンテナ部２３のメインローブＰ１１のパターンに相当する。アンテナパターン記憶部３６は、１種類のアンテナパターンＡＰを記憶している。
FIG. 8 is a graph showing the antenna pattern AP of the
このアンテナパターンＡＰは、たとえば、信号処理装置３の製造時にアンテナパターン記憶部３６に記憶されている。なお、アンテナパターンＡＰは、信号処理装置３の流通時にアンテナパターン記憶部３６に記憶されてもよい。また、アンテナパターンＡＰは、信号処理装置３が自船１０に設置された後にアンテナパターン記憶部３６に記憶されてもよい。
The antenna pattern AP is stored in the antenna
アンテナパターンＡＰは、自船１０回りの角度θと、信号の送受信強度との関係として示されている。角度θは、アンテナパターンの正面角度であり、平面視において回転アンテナ部２３からの信号送信方向の角度をゼロとして規定されている。アンテナパターンＡＰ＝ｈｎ ２（φ）＝ｈ２（θｎ，φ）である。アンテナパターンＡＰは、十分に強い指向性を有するように設定されている。たとえば、アンテナパターンＡＰのビーム幅は、−３ｄＢにおいて２度程度に設定されている。アンテナパターン記憶部３６は、アンテナパターンＡＰのデータを、平均ステアリングベクトル算出部３７へ出力する。
The antenna pattern AP is shown as a relationship between the angle θ around the
［平均ステアリングベクトル算出部の構成］
平均ステアリングベクトル算出部３７は、本発明の「ステアリングベクトル算出部」の一例である。平均ステアリングベクトル算出部３７は、回転アンテナ部２３のステアリングベクトルとしての平均ステアリングベクトルＡ(＿)（Ａ(＿)｛１｝（φ），Ａ(＿)｛２｝（φ），Ａ(＿)｛３｝（φ））を算出する。本実施形態では、平均ステアリングベクトル算出部３７は、アンテナパターンＡＰ＝ｈｎ ２（φ）に基づいて、平均ステアリングベクトルＡ(＿)（φ）を算出する。
[Configuration of average steering vector calculation unit]
The average steering
具体的には、平均ステアリングベクトル算出部３７は、領域Ｆ１，Ｆ２，Ｆ３に対応する平均ステアリングベクトルＡ(＿)｛１｝（φ），Ａ(＿)｛２｝（φ），Ａ(＿)｛３｝（φ）を算出する。
Specifically, the average steering
平均ステアリングベクトル算出部３７は、平均ステアリングベクトルＡ(＿)（Ａ(＿)｛１｝（φ），Ａ(＿)｛２｝（φ），Ａ(＿)｛３｝（φ））のデータを、角度スペクトラム算出部３５へ出力する。
The average steering
［角度スペクトラム算出部の構成］
角度スペクトラム算出部３５は、相関行列Ｒｘｘに基づく超解像処理を行うことで、領域Ｆ１，Ｆ２，Ｆ３について、個別に角度スペクトラムＰＳＳＡ（φ）を計算する。具体的には、角度スペクトラム算出部３５は、領域Ｆ１の角度スペクトラムＰ｛１｝ ＳＳＡ（φ）を、対応する固有ベクトルｑ｛１｝ １，…，ｑ｛１｝ Ｍおよび平均ステアリングベクトルＡ(＿)｛１｝（φ）を用いて、算出する。同様に、角度スペクトラム算出部３５は、領域Ｆ２，Ｆ３の角度スペクトラムＰ｛２｝ ＳＳＡ（φ），Ｐ｛３｝ ＳＳＡ（φ）を、対応する固有ベクトルｑ｛２｝ １，…，ｑ｛２｝ Ｍ；ｑ｛３｝ １，…，ｑ｛３｝ Ｍおよび平均ステアリングベクトルＡ(＿)｛２｝（φ），Ａ(＿)｛３｝（φ）を用いて、算出する。
[Configuration of angle spectrum calculation unit]
The angle
本実施形態では、スペクトラム算出部３５は、領域Ｆ１，Ｆ２，Ｆ３毎に、予め定められた数のピークが現れる手法を選択し、この手法を用いて角度スペクトラムＰ｛１｝ ＳＳＡ（φ），角度スペクトラムＰ｛２｝ ＳＳＡ（φ），Ｐ｛３｝ ＳＳＡ（φ）を算出する。この手法は、上記の予め定められた数が１の場合には後で詳述するＳＳＡ法であり、上記の予め定められた数が複数の場合（すなわち、２以上の自然数の場合）にはＭＵＳＩＣ法である。
In the present embodiment, the
まず、角度スペクトラム算出部３５による、角度スペクトラムＰ｛１｝ ＳＳＡ（φ）の算出について説明する。図９は、角度スペクトラムＰ｛１｝ ＳＳＡ（φ）の処理について説明するためのグラフである。図９のグラフの横軸は、方位φである。また、図９のグラフの縦軸は、信号強度を示している。図９は、ステアリングベクトルＡ｛１｝（φ）の各要素ｈ１ ２（φ）〜ｈｎ ２（φ）のうち、ｈ１ ２（φ），ｈｎ ２（φ）を図示している。
First, the calculation of the angle spectrum P {1} SSA (φ) by the angle
図２および図９を参照して、角度スペクトラム算出部３５は、ピークの数を１として設定した場合、ＭＵＳＩＣ法に基づく次式（４）を用いて、角度スペクトラムＰ｛１｝ ＳＳＡ（φ）を算出する。この場合、スペクトラムを算出するためのＳＳＡ (Single Signal Arrangement) 法と称する方法が用いられることで、角度スペクトラム算出部３５は、角度スペクトラムＰ｛１｝ ＳＳＡ（φ）を算出する。
なお、Ａ(＿)｛１｝Ｈ（φ）は、平均ステアリングベクトルＡ(＿)｛１｝（φ）の複素共役転置である。 Note that A(_) {1}H (φ) is a complex conjugate transpose of the average steering vector A(_) {1} (φ).
式（４）は、固有ベクトルｑ｛１｝ ２〜ｑ｛１｝ Ｍの値に拘わらず、固有ベクトルｑ｛１｝ ２〜ｑ｛１｝ Ｍの全てを、ノイズ領域のベクトルとして扱っている。即ち、式（４）では、領域Ｆ１に含まれる到来波の波源は、予め定められた数としての１のみであるとみなして、角度スペクトラムＰ｛１｝ ＳＳＡ（φ）を算出している。その結果、角度スペクトラムＰ｛１｝ ＳＳＡ（φ）に含まれるピークの数は、１つである。本願発明者は、ＭＵＳＩＣ法をベースに上述のＳＳＡ法を編み出した。なお、前述したように、上記の予め定められた数が２以上の自然数である場合、上記の式（４）に代えて、すなわち、ＳＳＡ (Single Signal Arrangement) 法と称する方法に代えて、ＭＵＳＩＣ（MUltiple SIgnal Classification）法で用いられる式を用いて、角度スペクトラムＰ｛１｝ ＳＳＡ（φ）を算出する。 Equation (4), regardless of the value of the eigenvector q {1} 2 ~q {1 } M, all eigenvectors q {1} 2 ~q {1 } M, are treated as a vector of the noise region. That is, in the equation (4), the angle spectrum P {1} SSA (φ) is calculated by regarding that the number of arriving waves included in the region F1 is only one, which is a predetermined number. As a result, the number of peaks included in the angle spectrum P {1} SSA (φ) is one. The present inventor has devised the above-mentioned SSA method based on the MUSIC method. As described above, when the predetermined number is a natural number of 2 or more, the MUSIC is replaced with the method called the SSA (Single Signal Arrangement) method instead of the expression (4). The angle spectrum P {1} SSA (φ) is calculated using the formula used in the (MUltiple SIgnal Classification) method.
角度スペクトラム算出部３５が上記式（４）を計算することで得られた角度スペクトラムＰ｛１｝ ＳＳＡ（φ）は、図９に示されている。角度スペクトラムＰ｛１｝ ＳＳＡ（φ）の波形は、山形形状である。なお、角度スペクトラムＰ｛１｝ ＳＳＡ（φ）の波形の一端ｅｎｄ１１は、ステアリングベクトルＡ｛１｝（φ）における１番目の要素ｈ１ ２（φ）の波形の一点に重なる。また、角度スペクトラムＰ｛１｝ ＳＳＡ（φ）の波形の他端ｅｎｄ１２は、ステアリングベクトルＡ｛１｝（φ）における最後の要素ｈｎ ２（φ）の波形の一点に重なる。また、角度スペクトラムＰ｛１｝ ＳＳＡ（φ）のピークは、φ＝約０．１°に存在しており、φ＝約０．１度に物標１１が存在していることを示している。
The angle spectrum P {1} SSA (φ) obtained by the angle
角度スペクトラム算出部３５は、前述したのと同様の構成で、領域Ｆ２に対応する角度スペクトラムＰ｛２｝ ＳＳＡ（φ）と、領域Ｆ３に対応する角度スペクトラムＰ｛３｝ ＳＳＡ（φ）と、を算出する。
The angle
次に、上記予め定められた数が１である場合の、角度スペクトラム算出部３５による、角度スペクトラムＰ｛２｝ ＳＳＡ（φ）の算出について説明する。図１０は、角度スペクトラムＰ｛２｝ ＳＳＡ（φ）の処理について説明するためのグラフである。図１０は、ステアリングベクトルＡ｛２｝（φ）の各要素ｈｎ＋１ ２（φ）〜ｈ２ｎ ２（φ）のうち、要素ｈｎ＋１ ２（φ），ｈ２ｎ ２（φ）を図示している。
Next, calculation of the angle spectrum P {2} SSA (φ) by the angle
図２および図１０を参照して、角度スペクトラム算出部３５は、領域Ｆ２に関する角度スペクトラムＰ｛２｝ ＳＳＡ（φ）を、式（４）を用いて算出する。ただし、この場合、式（４）における｛１｝は、｛２｝に置き換えられる。
With reference to FIGS. 2 and 10, the angle
角度スペクトラム算出部３５が上記式（４）を計算することで、得られた角度スペクトラムＰ｛２｝ ＳＳＡ（φ）は、図１０に示されている。角度スペクトラムＰ｛２｝ ＳＳＡ（φ）の波形は、山形形状である。なお、角度スペクトラムＰ｛２｝ ＳＳＡ（φ）の波形の一端ｅｎｄ２１は、ステアリングベクトルＡ｛２｝（φ）における１番目の要素ｈ２ ｎ＋１（φ）の波形の一点に重なる。また、角度スペクトラムＰ｛２｝ ＳＳＡ（φ）の波形の他端ｅｎｄ２２は、ステアリングベクトルＡ｛２｝（φ）における最後の要素ｈ２ｎ ２（φ）の波形の一点に重なる。また、角度スペクトラムＰ｛２｝ ＳＳＡ（φ）のピークは、φ＝約０．１°に存在しており、φ＝約０．１度に物標１１が存在していることを示している。
FIG. 10 shows the angle spectrum P {2} SSA (φ) obtained by the
次に、上記予め定められた数が１である場合の、角度スペクトラム算出部３５による、角度スペクトラムＰ｛３｝ ＳＳＡ（φ）の算出について説明する。図１１は、角度スペクトラムＰ｛３｝ ＳＳＡ（φ）の処理について説明するためのグラフである。図１１は、ステアリングベクトルＡ｛３｝（φ）の各要素ｈ２ｎ＋１ ２（φ）〜ｈ３ｎ ２（φ）のうち、要素ｈ２ｎ＋１ ２（φ），ｈ３ｎ ２（φ）を図示している。
Next, the calculation of the angle spectrum P {3} SSA (φ) by the angle
図２および図１１を参照して、角度スペクトラム算出部３５は、領域Ｆ３に関する角度スペクトラムＰ｛３｝ ＳＳＡ（φ）を、式（４）を用いて算出する。ただし、この場合、式（４）における｛１｝は、｛３｝に置き換えられる。
With reference to FIG. 2 and FIG. 11, the angle
角度スペクトラム算出部３５が上記式（４）を計算することで、得られた角度スペクトラムＰ｛３｝ ＳＳＡ（φ）は、図１１に示されている。角度スペクトラムＰ｛３｝ ＳＳＡ（φ）の波形は、山形形状である。なお、角度スペクトラムＰ｛３｝ ＳＳＡ（φ）の波形の一端ｅｎｄ３１は、ステアリングベクトルＡ｛３｝（φ）における１番目の要素ｈ２ｎ＋１ ２（φ）の波形の一点に重なる。また、角度スペクトラムＰ｛３｝ ＳＳＡ（φ）の波形の他端ｅｎｄ３２は、ステアリングベクトルＡ｛３｝（φ）における最後の要素ｈ３ｎ ２（φ）の波形の一点に重なる。また、角度スペクトラムＰ｛３｝ ＳＳＡ（φ）のピークは、φ＝約０．１°に存在しており、φ＝約０．１度に物標１１が存在していることを示している。
An angle spectrum P {3} SSA (φ) obtained by the angle
なお、角度スペクトラムＰ｛１｝ ＳＳＡ（φ），Ｐ｛２｝ ＳＳＡ（φ），Ｐ｛３｝ ＳＳＡ（φ）を１つのグラフで表示すると、図１２のようになる。図１２は、角度スペクトラム算出部３５で算出された角度スペクトラムＰ｛１｝ ＳＳＡ（φ），Ｐ｛２｝ ＳＳＡ（φ），Ｐ｛３｝ ＳＳＡ（φ）を説明するための図である。
Note that the angle spectra P {1} SSA (φ), P {2} SSA (φ), and P {3} SSA (φ) are displayed in one graph as shown in FIG. 12. FIG. 12 is a diagram for explaining the angle spectra P {1} SSA (φ), P {2} SSA (φ), and P {3} SSA (φ) calculated by the angle
図１２では、ステアリングベクトルＡ｛１｝（φ）〜Ａ｛３｝（φ）のそれぞれの要素ｈ１ ２（φ）〜ｈ３ｎ ２（φ）のうち、ｈ１ ２（φ），ｈ３ｎ／２ ２（φ），ｈ３ｎ ２（φ）が例示されている。
In Figure 12, steering vectors A {1} (φ) ~A {3} (φ) of each element h 1 2 (φ) ~h 3n 2 (φ),
図２および図１２を参照して、角度スペクトラムＰ｛１｝ ＳＳＡ（φ），Ｐ｛２｝ ＳＳＡ（φ），Ｐ｛３｝ ＳＳＡ（φ），は、何れも、１つのピークが存在していることを示している。前述したように、本実施形態では、１つの物標１１を探知する場合を説明しているので、信号処理装置３は、物標１１を精度よく探知できていることとなる。角度スペクトラム算出部３５は、算出した各角度スペクトラムＰＳＳＡ（φ）（Ｐ｛１｝ ＳＳＡ（φ），Ｐ｛２｝ ＳＳＡ（φ），Ｐ｛３｝ ＳＳＡ（φ））のデータを、表示用処理部３８へ出力する。
With reference to FIG. 2 and FIG. 12, each of the angular spectra P {1} SSA (φ), P {2} SSA (φ), and P {3} SSA (φ) has one peak. It indicates that As described above, in the present embodiment, the case where one
［表示用処理部の構成］
表示用処理部３８は、各角度スペクトラムＰ｛１｝ ＳＳＡ（φ），Ｐ｛２｝ ＳＳＡ（φ），Ｐ｛３｝ ＳＳＡ（φ）のデータを基に、表示部６に表示される画像データＧを生成するように構成されている。
[Configuration of display processing unit]
The
表示用処理部３８は、フィルタ処理部４１を有している。
The
［フィルタ処理部の構成］
フィルタ処理部４１は、各角度スペクトラムＰ｛１｝ ＳＳＡ（φ），Ｐ｛２｝ ＳＳＡ（φ），Ｐ｛３｝ ＳＳＡ（φ）のなかから、一部の角度スペクトラムを選択するように構成されている。
[Configuration of filter processing unit]
The
本実施形態では、方位φ＝約０．１度の位置に、各角度スペクトラムＰ｛１｝ ＳＳＡ（φ），Ｐ｛２｝ ＳＳＡ（φ），Ｐ｛３｝ ＳＳＡ（φ）の信号強度のピークＰｅａｋ｛１｝，Ｐｅａｋ｛２｝，Ｐｅａｋ｛３｝が存在している。なお、以下では、ピークＰｅａｋ｛１｝，Ｐｅａｋ｛２｝，Ｐｅａｋ｛３｝などのピークを総称していう場合には、単に、ピークＰｅａｋという。たとえば、方位方向Ｄ１における所定角度範囲内に２つ以上のピークが存在している場合、フィルタ処理部４１は、これらのピークＰｅａｋのうちの何れか１つ、たとえば、Ｐｅａｋ｛３｝を、信号強度の変化の急峻さなどに基づいて選択する。そして、フィルタ処理部４１は、選択したピークＰｅａｋを有する角度スペクトラムＰＳＳＡ（φ）のデータを、操作・表示装置４の表示部６へ出力する。
In the present embodiment, the signal strength of each angular spectrum P {1} SSA (φ), P {2} SSA (φ), P {3} SSA (φ) is located at the position of azimuth φ=about 0.1 degree. The peaks Peak {1} , Peak {2} , and Peak {3} are present. In the following description, peaks such as Peak {1} , Peak {2} , and Peak {3} are collectively referred to simply as Peak Peak. For example, when there are two or more peaks within the predetermined angle range in the azimuth direction D1, the
表示部６は、たとえば、ＰＰＩ表示装置である。表示部６は、信号処理装置３から出力された画像データＧに基づき、表示画面に画像を表示する。角度スペクトラムＰ｛３｝ ＳＳＡ（φ）の画像データＧは、方位方向Ｄ１に沿う画像１ライン分の濃淡データになる（角度スペクトラムＰ｛３｝ ＳＳＡ（φ）の信号強度が、その色と対応する）。自船１０の周囲の全周に亘って、目的の距離までの角度スペクトラムＰ｛３｝ ＳＳＡ（φ）を取得すれば、角度−距離の２次元濃淡画像を形成できる。たとえば、信号処理装置３が、角度−距離の２次元濃淡画像のデータを座標変換することで、水平距離−垂直距離の２次元画像データを生成することができる。これにより、図１３（ａ）に示されているような画像が、表示部６に表示される。
The
図１３（ａ）は、表示部６に表示される画像の一例を示している。表示部６に表示された映像では、物標１１の物標像ｅ１１は、方位方向Ｄ１における幅が狭くされており、方位方向Ｄ１における分解能が高い。すなわち、信号処理装置３は、物標１１を探知する分解能が十分に高い。なお、図１３（ｂ）は、信号処理装置３に代えて、公知のマルチビームフォーマ法によって角度スペクトラムを出力する信号処理装置を用いた場合の、表示部に表示される画像の一例を示している。図１３（ａ）および図１３（ｂ）は、ハッチングの斜線の間隔が短いほど、角度スペクトラムの信号強度が大きいことを示している。図１３（ａ）および図１３（ｂ）から明らかなように、公知の信号処理装置を用いて得られた物標像ｅ１１’は、方位方向Ｄ１にぼやけており、方位方向Ｄ１に必要以上に広い。この場合、方位方向Ｄ１における分解能は、低い。これに対して、信号処理装置３を用いて得られた物標像ｅ１１の輪郭は、より明確である。
FIG. 13A shows an example of an image displayed on the
［信号処理装置における処理の流れ］
次に、信号処理装置３における処理の流れの一例について、図１４を参照しつつ説明する。図１４は、信号処理装置３における処理の流れの一例を説明するためのフローチャートである。なお、フローチャートを参照して説明する場合は、フローチャート以外の図も適宜参照する。
[Process Flow in Signal Processing Device]
Next, an example of the processing flow in the
図１４を参照して、信号処理装置３の方位領域設定部３１が、受信部２４から受信データｘｉを読み出す（ステップＳ１０１）。次に、方位領域設定部３１は、所定の経過時間番号ｊについて、方位方向Ｄ１における領域Ｆ（本実施形態では、Ｆ１，Ｆ２，Ｆ３）を設定する（ステップＳ１０２）。本実施形態では、方位領域設定部３１は、領域Ｆ１のベクトルデータｘ｛１｝と、領域Ｆ２のベクトルデータｘ｛２｝と、領域Ｆ３のベクトルデータｘ｛３｝と、を設定する。
Referring to FIG. 14, azimuth
次に、相関行列算出部３２が、相関行列を算出する（ステップＳ１０３）。より具体的には、相関行列算出部３２は、ベクトルデータｘ｛１｝，ｘ｛２｝，ｘ｛３｝を基に、相関行列Ｒ｛１｝ ｘｘ，Ｒ｛２｝ ｘｘ，Ｒ｛３｝ ｘｘを算出する。
Next, the correlation
次に、空間平均処理部３３が、相関行列Ｒ｛１｝ ｘｘ，Ｒ｛２｝ ｘｘ，Ｒ｛３｝ ｘｘのそれぞれについて、空間平均相関行列Ｒ(＿)｛１｝ ｘｘ，Ｒ(＿)｛２｝ ｘｘ，Ｒ(＿)｛３｝ ｘｘを算出する（ステップ１０４）。
Next, the spatial
次に、固有ベクトル算出部３４が、空間平均相関行列Ｒ(＿)｛１｝ ｘｘ，Ｒ(＿)｛２｝ ｘｘ，Ｒ(＿)｛３｝ ｘｘのそれぞれについて、固有値分解を行う。これにより、固有ベクトル算出部３４は、固有ベクトルｑ｛１｝〜ｑ｛１｝ Ｍ，ｑ｛２｝ １〜ｑ｛２｝ Ｍ，ｑ｛３｝ １〜ｑ｛３｝ Ｍを算出する（ステップＳ１０５）。
Next, the
また、信号処理装置３において、平均ステアリングベクトル算出部３７が、平均ステアリングベクトルＡ(＿)｛１｝（φ），Ａ(＿)｛２｝（φ），Ａ(＿)｛３｝（φ）を算出する（ステップＳ１０６）。
Further, in the
次に、角度スペクトラム算出部３５は、領域Ｆ１，Ｆ２，Ｆ３について、個別に角度スペクトラムＰ｛１｝ ＳＳＡ（φ），Ｐ｛２｝ ＳＳＡ（φ），Ｐ｛３｝ ＳＳＡ（φ）を算出する（ステップＳ１０７）。
Next, the angle
次に、表示用処理部３８のフィルタ処理部４１が、表示用フィルタ処理を行う（ステップＳ１０８）。具体的には、フィルタ処理部４１は、角度スペクトラムＰ｛１｝ ＳＳＡ（φ），Ｐ｛２｝ ＳＳＡ（φ），Ｐ｛３｝ ＳＳＡ（φ）のなかから、１つの角度スペクトラムを適宜選択する。本実施形態では、たとえば、前述したように、フィルタ処理部４１は、角度スペクトラムＰ｛３｝ ＳＳＡ（φ）を、表示用の角度スペクトラムとして選択する。フィルタ処理部４１は、この角度スペクトラムＰ｛３｝ ＳＳＡ（φ）に基づく画像データＧを生成し、この画像データＧを表示部６へ出力する。
Next, the
信号処理装置３は、経過時間番号ｊ毎に、上記の処理を繰り返すことで、経過時間番号毎の画像データを生成する。操作・表示装置４の表示部６は、これらの画像データを基に、自船１０を中心とする、二次元のエコー画像を表示する。
The
［プログラム］
本実施形態にかかるプログラムは、コンピュータに、信号処理装置３の処理を実行させるプログラムであればよい。このプログラムをコンピュータにインストールし、実行することによって、本実施形態における信号処理装置３と、信号処理方法と、を実現することができる。この場合、コンピュータのＣＰＵ（Central Processing Unit）は、方位領域設定部３１、相関行列算出部３２、空間平均処理部３３、固有ベクトル算出部３４、角度スペクトラム算出部３５、平均ステアリングベクトル算出部３７、および、表示用処理部３８として機能し、処理を行う。なお、ＣＰＵ以外に、ＤＳＰ、ＳＰＵ、ＡＳＩＣ、ＦＰＧＡ、ＰＡＬなどが用いられてもよい。また、コンピュータのＲＯＭ（Read Only Memory）は、アンテナパターン記憶部３６として機能する。なお、信号処理装置３は、このようにソフトウェアとハードウェアとの協働によって実現されてもよいし、ハードウェアによって実現されてもよい。また、本実施形態にかかるプログラムは、ＤＶＤ（Digital Versatile Disc）等の記録媒体に記録された状態で流通されてもよいし、有線又は無線を用いた通信回線によって流通されてもよい。
[program]
The program according to the present embodiment may be any program that causes a computer to execute the processing of the
以上説明したように、信号処理装置３によると、相関行列算出部３２は、方位方向Ｄ１に沿って変位し且つ距離方向Ｄ２からの信号を順次受信する回転アンテナ部２３で得られた受信信号Ｓ１１に基づいて、相関行列Ｒｘｘ（Ｒ｛１｝ ｘｘ，Ｒ｛２｝ ｘｘ，Ｒ｛３｝ ｘｘ）を算出する。また、角度スペクトラム算出部３５は、相関行列Ｒｘｘに基づいて超解像処理を行うことで、角度スペクトラムＰＳＳＡ（φ）（Ｐ｛１｝ ＳＳＡ（φ），Ｐ｛２｝ ＳＳＡ（φ），Ｐ｛３｝ ＳＳＡ（φ））を算出する。このような構成により、信号処理装置３は、機械的に変位する回転アンテナ部２３を用いて得られた受信信号Ｓ１１の処理に際して、より高い分解能を実現できる。
As described above, according to the
より具体的には、信号処理装置３は、１つの回転アンテナ部２３による複数回のスイープによって得られた受信信号Ｓ１１（受信データｘｉ）を、アレイアンテナによって得られた受信信号と同様に扱うことができる。これにより、信号処理装置３は、アレイアンテナによって得られた受信信号に対する超分解能処理と同様の処理を、１つの回転アンテナ部２３によって得られた受信信号Ｓ１１（受信データｘｉ）に対して行うことができる。したがって、信号処理装置３は、より高い分解能を有する角度スペクトラムＰＳＳＡ（φ）を算出できる。また、本実施形態では、信号処理装置３での処理によって、高分解能化を達成する構成であり、回転アンテナ部２３の全長を長くすることで送信ビーム幅を狭くして高分解能化を達成する構成ではない。このため、回転アンテナ部２３として、全長の比較的短いアンテナを用いることができる。これにより、小型であることから安価な回転アンテナを、回転アンテナ部２３として用いることができる。さらに、全長の短い回転アンテナ部２３が採用されることで、小型船舶などの設置スペースに制限がある場所への回転アンテナ部２３の設置を実現できる。
More specifically, the
また、信号処理装置３によると、角度スペクトラム算出部３５は、相関行列Ｒｘｘ、およびステアリングベクトルとしての平均ステアリングベクトルＡ(＿)（φ）（Ａ(＿)｛１｝（φ），Ａ(＿)｛２｝（φ），Ａ(＿)｛３｝（φ）]に基づいて、角度スペクトラムＰＳＳＡ（φ）を算出する。より具体的には、本実施形態では、角度スペクトラム算出部３５は、固有ベクトルｑ、および、平均ステアリングベクトルＡ(＿)（φ）に基づいて、角度スペクトラムＰＳＳＡ（φ）を算出する。このような構成により、信号処理装置３は、アレイアンテナで得られた受信信号への処理と同様の処理で、高分解能処理を実現できる。よって、信号処理装置３は、受信信号Ｓ１１に含まれるノイズを、より確実に除去できる。
Further, according to the
また、信号処理装置３によると、平均ステアリングベクトル算出部３７は、回転アンテナ部２３のアンテナパターンＡＰに基づいて、平均ステアリングベクトルＡ(＿)（φ）を算出する。このような構成により、適応ビームフォーミング法を用いる場合のステアリングベクトルに相当するステアリングベクトルを、平均ステアリングベクトル算出部３７で算出することができる。
Further, according to the
また、信号処理装置３によると、アンテナパターンＡＰは、回転アンテナ部２３のメインローブＰ１１のパターンに相当する。このような構成により、信号処理装置３は、適切なステアリングベクトル（本実施形態では、平均ステアリングベクトルＡ(＿)（φ））を算出できる。
Further, according to the
また、信号処理装置３によると、角度スペクトラム算出部３５は、方位方向Ｄ１に沿った領域Ｆ毎に、角度スペクトラムＰＳＳＡ（φ）を算出する。より具体的には、信号処理装置３において、角度スペクトラム算出部３５は、領域Ｆ毎に、対応する固有ベクトルｑと平均ステアリングベクトルＡ(＿)（φ）とを用いて、各角度スペクトラムＰＳＳＡ ｛１｝（φ），ＰＳＳＡ ｛２｝（φ），ＰＳＳＡ ｛３｝（φ）を算出する。このような構成により、角度スペクトラム算出部３５は、方位方向Ｄ１に沿った範囲が狭く設定された領域Ｆ毎に角度スペクトラムＰＳＳＡ（φ）を算出すればよい。よって、角度スペクトラム算出部３５の単位時間あたり計算負荷を、より少なくできる。
Further, according to the
また、信号処理装置３によると空間平均処理部３３は、領域Ｆ毎に相関行列Ｒｘｘを平均化することで、空間平均相関行列Ｒ(＿)ｘｘを算出する。また、平均ステアリングベクトル算出部３７は、各領域Ｆにおいて平均化された平均ステアリングベクトルＡ(＿)（φ）を算出する。このような構成により、信号処理装置３は、相関行列Ｒｘｘに関するランクが回復された状態で、角度スペクトラムＰＳＳＡ（φ）を算出することができる。これにより、信号処理装置３は、より正確に、受信信号Ｓ１１で特定される到来波の到来方向を検出することができる。
Further, according to the
また、信号処理装置３によると、角度スペクトラム算出部３５は、領域Ｆ毎に、予め定められた数（たとえば、１）のピークが現れるように角度スペクトラムＰＳＳＡ（φ）を算出する。このような構成により、角度スペクトラム算出部３５の計算負荷を、より少なくできる。
Moreover, according to the
より具体的には、信号処理装置３によると、角度スペクトラム算出部３５は、相関行列Ｒｘｘにおける固有ベクトルｑの数に拘わらず、１つのピークが現れるとみなして、角度スペクトラムＰＳＳＡ（φ）を算出する。このような構成により、ヌル走査の手法を用いた、より分解能が高く且つ計算負荷の小さい、信号処理装置３を実現できる。
More specifically, according to the
また、信号処理装置３によると、フィルタ処理部４１は、領域Ｆ毎に算出された複数の角度スペクトラムＰＳＳＡ ｛１｝（φ），ＰＳＳＡ ｛２｝（φ），ＰＳＳＡ ｛３｝（φ）のなかから、一部の角度スペクトラムＰＳＳＡ ｛３｝（φ）を選択する。このような構成により、信号処理装置３は、表示部６での表示に適した角度スペクトラムＰＳＳＡ ｛３｝（φ）を選択できる。
Further, according to the
［第２実施形態］
図１５は、本発明の第２実施形態にかかるレーダ装置１Ａの信号処理装置３Ａの構成を示すブロック図である。本実施形態では、信号処理装置３Ａの表示用処理部３８Ａの構成が、第１実施形態における表示用処理部３８の構成と異なっている。この点以外について、信号処理装置３Ａは、信号処理装置３と同様の構成を有している。
[Second Embodiment]
FIG. 15 is a block diagram showing the configuration of the signal processing device 3A of the
表示用処理部３８Ａは、ビームパターン加工部３９を有している。ビームパターン加工部３９は、角度スペクトラムＰＳＳＡ（φ）で特定される包絡線の形状を加工するように構成されている。
The
図１６（ａ）は、自船１０および自船１０の周囲に存在する物標１３，１４を示す、模式的な平面図である。第２実施形態では、図１６（ａ）に示すように、自船１０から所定の距離に２つの物標１３，１４が存在している状態での、レーダ装置１Ａを例に説明する。図１６（ｂ）は、送受信装置２の回転アンテナ部２３の回転に伴って得られた受信データｘＡ（φｉ，Ｒｊ）の包絡線を示している。
FIG. 16A is a schematic plan view showing the
信号処理装置３Ａの方位領域設定部３１は、方位方向Ｄ１に沿った方位領域を設定する。これにより、第２実施形態では、たとえば、２つの領域Ｆ１Ａ，Ｆ２Ａが設定される。次に、相関行列算出部３２は、領域Ｆ１Ａのベクトルデータｘ｛１Ａ｝の相関行列Ｒ｛１Ａ｝ ｘｘと、領域Ｆ２Ａのベクトルデータｘ｛２Ａ｝の相関行列Ｒ｛２Ａ｝ ｘｘと、を算出する。
The azimuth
次に、空間平均処理部３３は、領域Ｆ１Ａにおける空間平均相関行列Ｒ(＿)｛１Ａ｝ ｘｘと、領域Ｆ２Ａにおける空間平均相関行列Ｒ(＿)｛２Ａ｝ ｘｘと、を算出する。固有ベクトル算出部３４は、空間平均相関行列Ｒ(＿)｛１Ａ｝ ｘｘ，Ｒ(＿)｛２Ａ｝ ｘｘのそれぞれについて、固有値分解を行う。これにより、固有ベクトル算出部３４は、固有ベクトルｑ｛１Ａ｝ １〜ｑ｛１Ａ｝ Ｍと、固有ベクトルｑ｛２Ａ｝ １〜ｑ｛２Ａ｝ Ｍとを算出する。
Next, the spatial
一方、平均ステアリングベクトル算出部３７は、アンテナパターン記憶部３６に格納されているアンテナパターンＡＰを基に、平均ステアリングベクトルＡ(＿)｛１Ａ｝（φ），Ａ(＿)｛２Ａ｝（φ）を算出する。
On the other hand, the average steering
角度スペクトラム算出部３５は、上記固有ベクトルｑ｛１Ａ｝ １〜ｑ｛１Ａ｝ Ｍと、平均ステアリングベクトルＡ(＿)｛１Ａ｝（φ）と、を式（４）に適用することで、角度スペクトラムＰ｛１Ａ｝ ＳＳＡ（φ）を算出する。ただし、この場合、式（４）における｛１｝は、｛１Ａ｝に置き換えられる。同様に、角度スペクトラム算出部３５は、上記固有ベクトルｑ｛２Ａ｝ １〜ｑ｛２Ａ｝-Ｍと、平均ステアリングベクトルＡ(＿)｛２Ａ｝（φ）と、を式（４）に適用することで、角度スペクトラムＰ｛２Ａ｝ ＳＳＡ（φ）を算出する。ただし、この場合、式（４）における｛１｝は、｛２Ａ｝に置き換えられる。
The angle
角度スペクトラム算出部３５は、角度スペクトラムＰ｛１Ａ｝ ＳＳＡ（φ），Ｐ｛２Ａ｝ ＳＳＡ（φ）のデータを、表示用処理部３８のビームパターン加工部３９へ出力する。
The angle
なお、角度スペクトラムＰ｛１Ａ｝ ＳＳＡ（φ）は、図１７に示されている。角度スペクトラムＰ｛１Ａ｝ ＳＳＡ（φ）の包絡線は、ピークＰｅａｋ｛１Ａ｝と、２つの極小値Ｍｉｎ｛１Ａ１｝，Ｍｉｎ｛１Ａ２｝と、を有している。以下では、角度スペクトラムＰＳＳＡ（φ）の包絡線のことを、単に角度スペクトラムという場合がある。方位方向Ｄ１において、極小値Ｍｉｎ｛１Ａ１｝，Ｍｉｎ｛１Ａ２｝が存在している方位の間に、ピークＰｅａｋ｛１Ａ｝の方位が存在している。このため、方位方向Ｄ１において、極小値Ｍｉｎ｛１Ａ１｝，Ｍｉｎ｛１Ａ２｝で挟まれた領域ΔＦ１Ａは、物標１３を特定する角度スペクトラムとして重要である。一方、方位方向Ｄ１において、領域ΔＦ１Ａ以外の領域は、角度スペクトラムとして重要ではない。
The angular spectrum P {1A} SSA (φ) is shown in FIG. The envelope of the angular spectrum P {1A} SSA (φ) has a peak Peak {1A} and two local minimum values Min {1A1} and Min {1A2} . Hereinafter, the envelope of the angle spectrum P SSA (φ) may be simply referred to as the angle spectrum. In the azimuth direction D1, the azimuth of the peak Peak {1A} exists between the azimuths having the minimum values Min {1A1} and Min {1A2} . Therefore, in the azimuth direction D1, the region ΔF1A sandwiched between the minimum values Min {1A1} and Min {1A2} is important as an angular spectrum for specifying the
そこで、ビームパターン加工部３９は、角度スペクトラムＰ｛１Ａ｝ ＳＳＡ（φ）について、領域ΔＦ１Ａ以外の成分を、削除する。この削除処理が行われた後の角度スペクトラムＰ｛１Ａ’｝ ＳＳＡ（φ）は、図１８に示されている。このように、角度スペクトラムＰ｛１Ａ’｝ ＳＳＡ（φ）の形状は、角度スペクトラムＰ｛１Ａ｝ ＳＳＡ（φ）の形状の一部を省略した形状となっている。
Therefore, the beam
また、角度スペクトラムＰ｛２Ａ｝ ＳＳＡ（φ）は、図１７に示されている。図１７に示されているように、角度スペクトラムＰ｛２Ａ｝ ＳＳＡ（φ）の曲線は、ピークＰｅａｋ｛２Ａ｝と、２つの極小値Ｍｉｎ｛２Ａ１｝，Ｍｉｎ｛２Ａ２｝と、を有している。角度スペクトラムＰ｛２Ａ｝ ＳＳＡ（φ）においても、方位方向Ｄ１において、極小値Ｍｉｎ｛２Ａ１｝，Ｍｉｎ｛２Ａ２｝で挟まれた領域ΔＦ２Ａは、ピークＰｅａｋ｛２Ａ｝の方位を含んでおり、物標１４を特定する角度スペクトラムとして重要である。一方、方位方向Ｄ１において、領域ΔＦ２Ａ以外の領域は、角度スペクトラムとして重要ではない。
The angle spectrum P {2A} SSA (φ) is shown in FIG. As shown in FIG. 17, the curve of the angular spectrum P {2A} SSA (φ) has a peak Peak {2A} and two local minimum values Min {2A1} , Min {2A2}. There is. Also in the angular spectrum P {2A} SSA (φ), in the azimuth direction D1, the region ΔF2A sandwiched between the minimum values Min {2A1} and Min {2A2} includes the azimuth of the peak Peak {2A} , and It is important as an angular spectrum for identifying the
そこで、表示用処理部３８Ａのビームパターン加工部３９は、角度スペクトラムＰ｛２Ａ｝ ＳＳＡ（φ）について領域ΔＦ２Ａ以外の成分を、削除する。この削除処理が行われた後の角度スペクトラムＰ２｛２Ａ’｝ ＳＳＡ（φ）は、図１９に示されている。図１９は、角度スペクトラムＰ｛１Ａ’｝ ＳＳＡ（φ），Ｐ｛２Ａ’｝ ＳＳＡ（φ）を示すグラフである。このように、角度スペクトラムＰ２｛２Ａ’｝ ＳＳＡ（φ）の形状は、角度スペクトラムＰ｛２Ａ｝ ＳＳＡ（φ）の形状の一部を省略した形状となっている。
Therefore, the beam
図２０は、角度スペクトラムＰＳＳＡ（φ）Ａ（Ｐ｛１Ａ’’｝ ＳＳＡ（φ），Ｐ｛２Ａ’’｝ ＳＳＡ（φ））を示すグラフである。図１９および図２０を参照して、更に、本実施形態では、角度スペクトラムＰ２｛１Ａ’’｝ ＳＳＡ（φ），Ｐ２｛２Ａ’’｝ ＳＳＡ（φ）のそれぞれの一部は、方位方向Ｄ１における位置が重なっている。この場合、ビームパターン加工部３９は、角度スペクトラムＰ２｛１Ａ’’｝ ＳＳＡ（φ），Ｐ２｛２Ａ’’｝ ＳＳＡ（φ）における互いの重なり領域ＲＡを削除する。これにより、ビームパターン加工部３９は、適宜角度スペクトラムＰ｛１Ａ｝ ＳＳＡ（φ），Ｐ｛２Ａ｝ ＳＳＡ（φ）の包絡線の形状を加工することで、図２０に示すような角度スペクトラムＰＳＳＡ（φ）Ａを算出する。角度スペクトラムＰＳＳＡ（φ）Ａは、表示対象となる物標１３，１４のエコー像の輪郭をより明確にするように構成されたスペクトラムである。
FIG. 20 is a graph showing the angular spectrum P SSA (φ)A(P {1A″} SSA (φ), P {2A″} SSA (φ)). With reference to FIGS. 19 and 20, further, in the present embodiment, a part of each of the angular spectra P2 {1A″} SSA (φ) and P2 {2A″} SSA (φ) is the azimuth direction D1. The positions in are overlapping. In this case, the beam
次に、表示用処理部３８Ａは、物標１３，１４のそれぞれに対応するピークＰｅａｋ｛１Ａ｝，ピークＰｅａｋ｛２Ａ｝を示す画像データＧＡを生成し、この画像データＧＡを表示部６へ出力する。
Next, the
図２１（ａ）は、信号処理装置３Ａからの映像信号ＧＡによって表示部６に表示される画像の一例を示す図である。図２１（ａ）を参照して、表示部６には、物標１３に対応する物標像ｅ１３と、物標１４に対応する物標像ｅ１４と、が示されている。図２１（ａ）から明らかなように、物標像ｅ１３の中心と、物標像ｅ１４の中心とは、明確に区別されている。よって、オペレータは、表示部６を通じて、隣接する物標１３，１４を、明確に識別できる。
FIG. 21A is a diagram showing an example of an image displayed on the
なお、図１９（ｂ）は、受信データｘＡｉが超解像処理されることなく表示部６へ出力された場合の、表示部６に表示される画像の一例を示す図である。この場合、物標像ｅ１３’と物標像ｅ１４’とは、区別不能な程度に融合している。このため、オペレータは、表示部を通じて、隣接する物標を識別することが困難である。
Note that FIG. 19B is a diagram showing an example of an image displayed on the
以上説明したように、本発明の第２実施形態にかかる信号処理装置３Ａによると、ビームパターン加工部３９は、角度スペクトラムＰＳＳＡ ｛１Ａ｝（φ），ＰＳＳＡ ｛２Ａ｝（φ）で特定される包絡線の形状を加工するように構成されている。このような構成により、信号処理装置３Ａは、隣接する２つの物標１３，１４についての物標像ｅ１３，ｅ１４を、より明確に区別できる態様で表示部６に表示させることができる。
As described above, according to the signal processing device 3A according to the second embodiment of the present invention, the beam
［第３実施形態］
図２２は、本発明の第３実施形態における信号処理装置３Ｂの構成を示すブロック図である。信号処理装置３Ｂが信号処理装置３Ａと異なっているのは、角度スペクトラムＰ｛１Ａ｝ ＳＳＡ（φ），Ｐ｛２Ａ｝ ＳＳＡ（φ）が、表示部６での表示のために、より一層適した波形に変換される点にある。
[Third Embodiment]
FIG. 22 is a block diagram showing the configuration of the
より具体的には、図２０および図２１（ａ）を参照して、信号処理装置３Ａが用いられた場合、表示部６において、隣接する物標像ｅ１３，ｅ１４は、明確に識別可能に表示される。しかしながら、角度スペクトラムＰ｛１Ａ’’｝ ＳＳＡ（φ），Ｐ｛２Ａ’’｝ ＳＳＡ（φ）の波形において、ピークＰｅａｋ｛１Ａ｝，ピークＰｅａｋ｛２Ａ｝への立ち上がり部分は、方位方向Ｄ１に沿った範囲が狭い。その結果、物標像ｅ１３，ｅ１４は、それぞれ、表示部６において、方位方向Ｄ１に沿った幅が極めて狭くなる。一方で、物標１３，１４の実際の形状は、方位方向Ｄ１に広がりを有している。このため、オペレータに対する違和感の少ない表示を表示部６で行うためには、物標像ｅ１３，１４の形状についても、物標１３，１４の形状に対応して方位方向Ｄ１に広がりを有していることが好ましい。このような効果を発揮するための構成が、信号処理装置３Ｂに設けられている。
More specifically, referring to FIGS. 20 and 21(a), when the signal processing device 3A is used, the adjacent target images e13 and e14 are displayed clearly and identifiable on the
図２０を参照して、具体的には、信号処理装置３Ｂは、信号処理装置３Ａの表示用処理部３８Ａに代えて、表示用処理部３８Ｂを有している。表示用処理部３８Ｂは、ビームパターン加工部３９Ｂを有している。また、信号処理装置３Ｂは、表示用パターン記憶部４２を有している。
Referring to FIG. 20, specifically, the
ビームパターン加工部３９Ｂは、角度スペクトラムＰＳＳＡ（φ）の形状を、予め定められたパターン形状に置き換える。また、ビームパターン加工部３９Ｂは、置き換えたパターン形状の信号強度のピークを、角度スペクトラムＰＳＳＡ（φ）の信号強度のピークＰｅａｋと同じに設定する。
The beam
たとえば、第２実施形態の場合と同様に、自船１０の周囲に物標１３，１４が存在している場合、信号処理装置３Ｂは、受信データｘＡｉを読み出す。そして、信号処理装置３Ｂは、信号処理装置３Ａと同様に、角度スペクトラムＰ｛１Ａ’｝ ＳＳＡ（φ），Ｐ｛２Ａ’｝ ＳＳＡ（φ）（図１９参照）を算出する。
For example, as in the case of the second embodiment, when the
次に、フィルタ処理部４１Ｂのビームパターン加工部３９Ｂは、各角度スペクトラムＰ｛１Ａ’｝ ＳＳＡ（φ），Ｐ｛２Ａ’｝ ＳＳＡ（φ）のピークＰｅａｋ｛１Ａ｝，Ｐｅａｋ｛２Ａ｝の値および方位φ１Ａ，φ２Ａを検出する。なお、図２３は、Ｐｅａｋ｛１Ａ｝，Ｐｅａｋ｛２Ａ｝の大きさと方位φ１Ａ，φ２Ａとの関係を示すグラフである。
Next, the beam
ビームパターン加工部３９Ｂは、表示用パターン記憶部４２に記憶されているパターンＰＴのデータを読み出す。パターンＰＴのデータは、１又は複数の波形のパターンのデータである。
The beam
本実施形態では、表示用パターン記憶部４２は、３つのパターンＰＴ（ＰＴ１，ＰＴ２，ＰＴ３）のデータを記憶している。
In the present embodiment, the display
パターンＰＴ１は、回転アンテナ部２３に関するアンテナパターンＡＰと同様のパターンである。パターンＰＴ２は、矩形波のパターンである。パターンＰＴ３は、三角波のパターンである。
The pattern PT1 is a pattern similar to the antenna pattern AP for the
パターンＰＴ１，ＰＴ２，ＰＴ３のビーム幅は、それぞれ、方位方向Ｄ１に所定の長さを有している。これにより、表示部６で方位方向Ｄ１に十分な長さのエコー像を表示させることが可能である。より具体的には、パターンＰＴ１，ＰＴ２，ＰＴ３のピークを中心とする所定範囲のビーム幅は、角度スペクトラムＰＳＳＡ ｛１Ａ’｝（φ）のピークＰｅａｋ｛１Ａ｝を中心とする上記所定範囲のビーム幅よりも大きい。同様に、パターンＰＴ１，ＰＴ２，ＰＴ３のピークを中心とする所定範囲のビーム幅は、角度スペクトラムＰＳＳＡ（φ）｛２Ａ’｝のピークＰｅａｋ｛２Ａ｝を中心とする上記所定範囲のビーム幅よりも大きい。
The beam widths of the patterns PT1, PT2, and PT3 each have a predetermined length in the azimuth direction D1. As a result, it is possible to display an echo image having a sufficient length in the azimuth direction D1 on the
パターンＰＴのデータは、たとえば、信号処理装置３Ｂの製造時に表示用パターン記憶部４２に記憶される。なお、パターンＰＴのデータは、信号処理装置３Ｂの流通時に表示用パターン記憶部４２に記憶されてもよい。また、パターンＰＴのデータは、信号処理装置３Ｂが自船１０に設置された後に表示用パターン記憶部４２に記憶されてもよい。
The data of the pattern PT is stored in the display
ビームパターン加工部３９Ｂは、パターンＰＴ１，ＰＴ２，ＰＴ３の何れかを選択する。この場合、ビームパターン加工部３９Ｂは、たとえば、予めオペレータによって設定されたパターンＰＴを選択する。なお、ビームパターン加工部３９Ｂは、角度スペクトラムＰ｛１Ａ’｝ ＳＳＡ（φ），Ｐ｛２Ａ’｝ ＳＳＡ（φ）の波形に応じて、適宜パターンＰＴを設定してもよい。本実施形態では、ビームパターン加工部３９Ｂは、パターンＰＴ１を設定する。
The beam
次に、図２４に示すように、ビームパターン加工部３９Ｂは、パターンＰＴ１のピークの大きさおよび方位を、ピークＰｅａｋ｛１Ａ｝の大きさおよび方位φ１Ａに設定する。同様に、ビームパターン加工部３９Ｂは、もう一つのパターンＰＴ１のピークの大きさおよび方位を、ピークＰｅａｋ｛２Ａ｝の大きさおよび方位φ２Ａに設定する。
Next, as shown in FIG. 24, the beam
これにより、ビームパターン加工部３９Ｂは、物標１３，１４の数と同じ数（２つ）のパターンＰＴ１，ＰＴ１を基にしたパターンのデータを生成する。次に、ビームパターン加工部３９Ｂは、このパターンにおいて方位方向Ｄ１の位置が重なる領域を削除する。その結果、ビームパターン加工部３９Ｂは、図２５に示すような、角度スペクトラムＰＳＳＡ（φ）Ｂを算出する。この角度スペクトラムＰＳＳＡ（φ）Ｂのデータに基づく画像データＧＢが表示部６に与えられた場合、表示部６は、図２６に示すような画像を表示する。
Accordingly, the beam
図２６は、信号処理装置３Ｂからの画像データＧＢによって表示部６に表示される画像の一例を示す図である。図２１（ａ）および図２６を参照して、表示部６の表示画面には、物標１３に対応する物標像ｅ１３Ｂと、物標１４に対応する物標像ｅ１４Ｂと、が示されている。図２６から明らかなように、物標像ｅ１３Ｂ，ｅ１４Ｂは、方位方向Ｄ１の幅が適度に確保されている。このため、信号処理装置３Ａの処理に基づいて得られる物標像ｅ１３，１４と比較して、物標像ｅ１３Ｂ，ｅ１４Ｂは、より自然な態様で、表示部６に表示される。
FIG. 26 is a diagram showing an example of an image displayed on the
なお、たとえば、物標１３，１４のなかで物標１３のみが、自船１０の周囲に存在している場合、角度スペクトラム角度スペクトラムＰＳＳＡ（φ）Ｂの波形は、図２７に示されるパターンＰＴ１のように、方位方向Ｄ１に対称となる。
Note that, for example, when only the
以上説明したように、本発明の第３実施形態にかかる信号処理装置３Ｂによると、ビームパターン加工部３９Ｂは、角度スペクトラムＰＳＳＡ（φ）の包絡線の形状を、予め定められたパターンＰＴの形状に置き換える。さらに、ビームパターン加工部３９Ｂは、パターンＰＴの形状の信号強度のピークを、角度スペクトラムＰＳＳＡ（φ）のピークＰｅａｋと同じに設定する。このような構成により、ビームパターン加工部３９は、角度スペクトラムＰＳＳＡ（φ）のピークＰｅａｋ｛１Ａ｝，Ｐｅａｋ｛２Ａ｝の値を変更することなく、物標１３，１４の物標像ｅ１３，ｅ１４を、より自然な形状となるようにすることができる。
As described above, according to the
また、信号処理装置３Ｂによると、パターンＰＴ（ＰＴ１，ＰＴ２，ＰＴ３）のピークを中心とする所定範囲のビーム幅は、角度スペクトラムＰＳＳＡ（φ）の包絡線におけるピークＰｅａｋ｛１Ａ｝，Ｐｅａｋ｛２Ａ｝を中心とする所定範囲の幅よりも大きく設定されている。角度スペクトラム算出部３５で算出された角度スペクトラムＰＳＳＡ ｛１Ａ｝（φ），ＰＳＳＡ ｛２Ａ｝（φ）の形状は、ピークＰｅａｋ｛１Ａ｝，Ｐｅａｋ｛２Ａ｝の周辺で信号強度が急峻に変化する形状である。その結果、角度スペクトラムＰＳＳＡ（φ）で特定される物標像の形状は、方位方向Ｄ１に広がりの少ない形状となり、オペレータに違和感を与えてしまう。そこで、ビームパターン加工部３９Ｂが、ピークＰｅａｋ｛１Ａ｝，Ｐｅａｋ｛２Ａ｝の周辺でのビーム幅を広くなるような処理を施す。これにより、角度スペクトラムＰＳＳＡ（φ）Ｂで特定される物標像ｅ１３Ｂ，ｅ１４Ｂの形状は、方位方向Ｄ１に広がりのある形状となり、表示部６で表示されたときに、オペレータに違和感を与えることを抑制できる。
Further, according to the
また、信号処理装置３Ｂによると、ビームパターン加工部３９Ｂは、オペレータの指示などに基づいて、複数のパターンＰＴのなかから、１つのパターンを選択することができる。その結果、画像データで表示される画像が、表示部６での表示に適した画像となるように、最適なパターンＰＴを選択できる。
Further, according to the
［第４実施形態］
図２８は、本発明の第４実施形態にかかる水中探知装置１００の構成を示すブロック図である。図２８を参照して、水中探知装置１００は、超音波探知装置であり、たとえば、自船に備えられる。水中探知装置１００は、水中の物標を探知するために用いられる。本実施形態において、水中探知装置１００は、サーチライトソナーであり、自船を中心とする所定の単位方位毎に、超音波の送受信を繰り返す。
[Fourth Embodiment]
FIG. 28 is a block diagram showing the configuration of the
水中探知装置１００は、送受信装置２Ｈと、信号処理装置３Ｈと、操作・表示装置４Ｈと、を備えている。
The
送受信装置２Ｈは、超音波を送受信する機能を有している。送受信装置２Ｈは、自船の船底に配置された取り付けられた受信装置としての超音波振動子１０１を、１つ有している。超音波振動子１０１は、鉛直軸線回りの方位方向を副走査方向として、自転するように構成されている。超音波振動子１０１は、自転した状態で、方位方向と交差する距離方向に向かって超音波を水中に送信するとともに、水中からのエコーを受信する。送受信装置２Ｈは、このエコー信号に基づく受信信号を電気信号に変換することで、受信データｘＨ（φｉ，Ｒｊ）信号を生成する。
The
このような構成により、送受信装置２は、自船回りの各方位において、距離方向に沿う超音波を順次送受信する。この受信データｘＨ（φｉ，Ｒｊ）は、単位方位毎に生成される。この受信データｘＨ（φｉ，Ｒｊ）は、信号処理装置３Ｈへ出力される。信号処理装置３Ｈは、受信データｘＨ（φｉ，Ｒｊ）を基に、画像データＧＨを生成する。信号処理装置３ＧＨは、この画像データＧＨを、操作・表示装置４へ出力する。操作・表示装置４は、画像データに基づく画像を表示する。
With such a configuration, the transmission/
このような構成によると、受信データｘＨ（φｉ，Ｒｊ）の処理に際し、より高い分解能を実現できる。 With such a configuration, a higher resolution can be realized when processing the received data xH(φ i , R j ).
なお、本実施形態では、超音波振動子の数が１である形態を例に説明した。しかしながら、この通りでなくてもよい。たとえば、複数の超音波振動子と、これらの超音波振動子が受信した信号の位相を一致させる整相器と、を設けてもよい。 It should be noted that the present embodiment has been described by taking the form in which the number of ultrasonic transducers is one as an example. However, this need not be the case. For example, a plurality of ultrasonic transducers and a phase adjuster for matching the phases of the signals received by these ultrasonic transducers may be provided.
なお、第４実施形態では、超音波探知装置の一例として、水中探知装置１００を例に説明した。しかしながら、この形態に限定されない。本発明は、超音波探知装置としての魚群探知機に適用されてもよい。
In the fourth embodiment, the
以上、本発明の実施形態について複数説明したけれども、本発明は上述の実施の形態に限らず、特許請求の範囲に記載した限りにおいて様々な変更が可能である。たとえば、次のように変更して実施してもよい。 Although a plurality of embodiments of the present invention have been described above, the present invention is not limited to the above-described embodiments, and various modifications can be made as long as they are set forth in the claims. For example, the following modifications may be performed.
（１）上述の各実施形態では、信号処理装置に種々の構成が備えられる形態を例に説明した。しかしながら、この通りでなくてもよい。信号処理装置は、少なくとも、相関行列算出部と、角度スペクトラム算出部とを備えていればよく、他の構成は、設けられていなくてもよい。 (1) In each of the above-described embodiments, the mode in which the signal processing device is provided with various configurations has been described as an example. However, this need not be the case. The signal processing device only needs to include at least a correlation matrix calculation unit and an angle spectrum calculation unit, and other configurations may not be provided.
（２）また、上述の各実施形態では、角度スペクトラム算出部が、ＭＵＳＩＣ法をベースとした評価関数を用いて角度スペクトラムを算出する構成を説明した。しかしながら、この通りでなくてもよい。角度スペクトラム算出部は、たとえば、Ｃａｐｏｎ法、Ｐｉｓａｒｅｎｋｏ法、最小ノルム法など他の適応ビームフォーミング法を用いて、角度スペクトラムを算出してもよい。 (2) Further, in each of the above-described embodiments, the configuration has been described in which the angle spectrum calculation unit calculates the angle spectrum using the evaluation function based on the MUSIC method. However, this need not be the case. The angle spectrum calculation unit may calculate the angle spectrum using another adaptive beamforming method such as the Capon method, the Pisarenko method, or the minimum norm method.
（３）また、上述の各実施形態では、信号処理装置が、方位方向において複数の領域を設定する形態を例に説明した。しかしながら、この通りでなくてもよい。たとえば、信号処理装置は、方位方向の３６０度に設定された１つの領域について、角度スペクトラムを算出してもよい。 (3) Further, in each of the above-described embodiments, the signal processing device has been described as an example in which a plurality of regions are set in the azimuth direction. However, this need not be the case. For example, the signal processing device may calculate the angle spectrum for one area set to 360 degrees in the azimuth direction.
（４）また、上述の各実施形態のうち、レーダ装置では、回転アンテナ部が、方位方向において１回に１つの受信信号を出力する形態を例に説明した。この限りにおいて、回転アンテナ部は、アレイアンテナであってもよいし、アレイアンテナでなくてもよい。また、上述の実施形態のうち、水中探知装置では、超音波振動子が方位方向において１回に１つの受信信号を出力する形態を例に説明した。この限りにおいて、超音波振動子は、アレイであってもよいし、アレイでなくてもよい。 (4) Further, among the above-described respective embodiments, in the radar device, the rotation antenna unit outputs one received signal at a time in the azimuth direction, as an example. As long as this is the case, the rotating antenna unit may or may not be an array antenna. Further, among the above-described embodiments, the underwater detection device has been described by taking the form in which the ultrasonic transducer outputs one reception signal at a time in the azimuth direction. As long as this is the case, the ultrasonic transducer may or may not be an array.
（５）また、上述の各実施形態では、信号処理装置が、レーダ装置または水中探知装置に適用された形態を例に説明した。しかしながら、この通りでなくてもよい。たとえば、信号処理装置は、超音波などを利用する医療機器、通信機器、ネットワーク機器、または、ＧＮＳＳ（Global Navigation Satellite System）システムの利用者受信機など、他の機器に備えられていてもよい。 (5) Moreover, in each of the above-described embodiments, the signal processing device is described as an example in which the signal processing device is applied to a radar device or an underwater detection device. However, this need not be the case. For example, the signal processing device may be provided in another device such as a medical device utilizing ultrasonic waves, a communication device, a network device, or a user receiver of a GNSS (Global Navigation Satellite System) system.
本発明は、信号処理装置、レーダ装置、水中探知装置、信号処理方法、および、プログラムとして、広く適用することができる。 The present invention can be widely applied as a signal processing device, a radar device, an underwater detection device, a signal processing method, and a program.
２ 送受信装置
３ 信号処理装置
２３ 回転アンテナ部（受信装置）
３１ 方位領域設定部（領域設定部）
３２ 相関行列算出部
３３ 空間平均処理部（相関行列平均化処理部）
３４ 固有ベクトル算出部
３５ 角度スペクトラム算出部（スペクトラム算出部）
３７ 平均ステアリングベクトル算出部（ステアリングベクトル算出部）
３９ ビームパターン加工部
４１ フィルタ処理部
１０１ 超音波振動子（受信装置）
Ａ（φ） ステアリングベクトル
Ａ(＿)（φ） 平均ステアリングベクトル
Ｄ１ 方位方向（副走査方向）
Ｄ２ 距離方向（主走査方向）
Ｆ 領域
Ｐｅａｋ ピーク
ＰＴ パターン
ＰＳＳＡ（φ） 角度スペクトラム（スペクトラム）
ｑ 固有ベクトル
Ｒ(＿)ｘｘ 空間平均相関行列（平均相関行列）
Ｒｘｘ 相関行列
Ｓ１１ 受信信号
2 transmitter/
31 Direction Area Setting Section (Area Setting Section)
32
34
37 Average Steering Vector Calculation Unit (Steering Vector Calculation Unit)
39 Beam
A(φ) Steering vector A(_)(φ) Average steering vector D1 Azimuth direction (sub-scanning direction)
D2 Distance direction (main scanning direction)
F region Peak Peak PT pattern P SSA (φ) Angle spectrum
q eigenvector R(_) xx spatial mean correlation matrix (mean correlation matrix)
R xx correlation matrix S11 received signal
Claims (16)
前記相関行列に基づいて超解像処理を行うことで、前記副走査方向に沿って得られる前記受信信号の信号強度のスペクトラムを算出する、スペクトラム算出部と、
前記受信装置のステアリングベクトルを算出する、ステアリングベクトル算出部と、
前記相関行列の固有ベクトルを算出する固有ベクトル算出部と、
を備え、
前記スペクトラム算出部は、前記固有ベクトル、および、前記ステアリングベクトルに基づいて、前記スペクトラムを算出することを特徴とする、信号処理装置。 Obtained along the sub-scanning direction based on a reception signal obtained by using a receiving device that is displaced along a predetermined sub-scanning direction and sequentially receives signals from the main scanning direction intersecting the sub-scanning direction. A correlation matrix calculating unit for calculating a correlation matrix of the plurality of received signals,
By performing super-resolution processing based on the correlation matrix, a spectrum of the signal intensity of the received signal obtained along the sub-scanning direction is calculated, a spectrum calculation unit ,
A steering vector calculation unit for calculating a steering vector of the receiving device;
An eigenvector calculation unit that calculates an eigenvector of the correlation matrix,
Bei to give a,
The signal processing device, wherein the spectrum calculation unit calculates the spectrum based on the eigenvector and the steering vector .
前記ステアリングベクトル算出部は、前記受信装置のアンテナパターンに基づいて前記ステアリングベクトルを算出することを特徴とする、信号処理装置。 The signal processing device according to claim 1 , wherein
The signal processing device, wherein the steering vector calculation unit calculates the steering vector based on an antenna pattern of the reception device.
前記アンテナパターンは、前記受信装置のメインローブのパターンに相当することを特徴とする、信号処理装置。 The signal processing device according to claim 2 , wherein
The signal processing device, wherein the antenna pattern corresponds to a main lobe pattern of the receiving device.
前記相関行列に基づいて超解像処理を行うことで、前記副走査方向に沿って得られる前記受信信号の信号強度のスペクトラムを算出する、スペクトラム算出部と、
を備え、
前記副走査方向に沿って複数の領域を設定する領域設定部をさらに備え、
前記相関行列算出部は、前記領域毎に前記相関行列を算出し、
前記スペクトラム算出部は、前記領域毎に、前記スペクトラムを算出することを特徴とする、信号処理装置。 Obtained along the sub-scanning direction based on a reception signal obtained by using a receiving device that is displaced along a predetermined sub-scanning direction and sequentially receives signals from the main scanning direction intersecting the sub-scanning direction. A correlation matrix calculating unit for calculating a correlation matrix of the plurality of received signals,
By performing super-resolution processing based on the correlation matrix, a spectrum of the signal intensity of the received signal obtained along the sub-scanning direction is calculated, a spectrum calculation unit,
Equipped with
Further comprising an area setting unit that sets a plurality of areas along the front Symbol sub-scanning direction,
The correlation matrix calculation unit calculates the correlation matrix for each of the regions,
The signal processing device, wherein the spectrum calculation unit calculates the spectrum for each of the regions.
前記相関行列の固有ベクトルを算出する固有ベクトル算出部と、
前記受信装置のステアリングベクトルを算出する、ステアリングベクトル算出部と、
をさらに備え、
前記固有ベクトル算出部および前記ステアリングベクトル算出部は、それぞれ、前記領域毎に前記固有ベクトルおよび前記ステアリングベクトルを算出し、
前記スペクトラム算出部は、前記領域毎に、前記固有ベクトルと前記ステアリングベクトルとを用いて前記スペクトラムを算出することを特徴とする、信号処理装置。 The signal processing device according to claim 4 , wherein
An eigenvector calculation unit that calculates an eigenvector of the correlation matrix,
A steering vector calculation unit for calculating a steering vector of the receiving device;
Further equipped with,
The eigenvector calculation unit and the steering vector calculation unit respectively calculate the eigenvector and the steering vector for each of the regions,
The signal processing device, wherein the spectrum calculation unit calculates the spectrum for each of the regions using the eigenvector and the steering vector.
相関行列平均化処理部をさらに備え、
前記相関行列平均化処理部は、前記領域毎に前記相関行列を平均化することで平均相関行列を算出し、
前記固有ベクトル算出部は、前記平均相関行列を固有値分解することで、前記固有ベクトルを算出し、
前記ステアリングベクトル算出部は、各前記領域において平均化された平均ステアリングベクトルを、前記ステアリングベクトルとして算出することを特徴とする、信号処理装置。 The signal processing device according to claim 5 , wherein
Further comprising a correlation matrix averaging processing unit,
The correlation matrix averaging processing unit calculates an average correlation matrix by averaging the correlation matrix for each of the regions,
The eigenvector calculation unit calculates the eigenvector by decomposing the average correlation matrix into eigenvalues,
The signal processing device, wherein the steering vector calculation unit calculates an average steering vector averaged in each of the regions as the steering vector.
前記スペクトラム算出部は、前記領域毎に、予め定められた数のピークが現れる手法を選択し、この手法を用いて前記スペクトラムを算出することを特徴とする、信号処理装置。 The signal processing device according to any one of claims 4 to 6 , wherein:
The signal processing device, wherein the spectrum calculation unit selects a method in which a predetermined number of peaks appear in each of the regions and calculates the spectrum using this method.
前記手法は、前記予め定められた数が１の場合にはＳＳＡ（Ｓｉｎｇｌｅ Ｓｉｇｎａｌ Ａｒｒａｎｇｅｍｅｎｔ）法であり、前記予め定められた数が複数の場合にはＭＵＳＩＣ法であることを特徴とする、信号処理装置。 The signal processing device according to claim 7 , wherein
The signal processing is characterized in that the method is an SSA (Single Signal Arrangement) method when the predetermined number is 1 and a MUSIC method when the predetermined number is plural. apparatus.
フィルタ処理部を更に備え、
前記フィルタ処理部は、前記領域毎に算出された複数の前記スペクトラムのなかから、一部の前記スペクトラムを選択することを特徴とする、信号処理装置。 The signal processing device according to any one of claims 4 to 8 , wherein:
Further equipped with a filter processing unit,
The signal processing device, wherein the filter processing unit selects a part of the spectrum from the plurality of spectra calculated for each area.
ビームパターン加工部をさらに備え、
前記ビームパターン加工部は、前記スペクトラムの包絡線の形状を加工することを特徴とする、信号処理装置。 The signal processing device according to any one of claims 1 to 9 , wherein
Further equipped with a beam pattern processing unit,
The signal processing device, wherein the beam pattern processing unit processes the shape of the envelope of the spectrum.
前記ビームパターン加工部は、前記包絡線の形状を、予め定められたパターン形状に置き換え、且つ、当該パターン形状の前記信号強度のピークを、前記スペクトラムの前記信号強度のピークと同じ値に設定することを特徴とする、信号処理装置。 The signal processing device according to claim 10 , wherein
The beam pattern processing unit replaces the shape of the envelope with a predetermined pattern shape, and sets the peak of the signal intensity of the pattern shape to the same value as the peak of the signal intensity of the spectrum. A signal processing device, characterized in that
前記パターン形状における前記ピークを中心とする所定範囲の幅は、前記スペクトラムの包絡線における前記ピークを中心とする所定範囲の幅よりも大きいことを特徴とする、信号処理装置。 The signal processing device according to claim 11 , wherein
A signal processing device, wherein a width of a predetermined range centered on the peak in the pattern shape is larger than a width of a predetermined range centered on the peak in the envelope of the spectrum.
請求項１乃至請求項１２の何れか１項に記載の信号処理装置と、を備え、
前記相関行列算出部は、前記アンテナ部を用いて得られた受信信号に基づいて、前記相関行列を算出することを特徴とする、レーダ装置。 A transmitting/receiving device, which is provided with an antenna unit that rotates in a predetermined azimuth direction as a sub-scanning direction and sequentially transmits and receives electromagnetic waves along a distance direction that intersects the azimuth direction,
A signal processing device according to any one of claims 1 to 12 ,
The radar device, wherein the correlation matrix calculation unit calculates the correlation matrix based on a received signal obtained using the antenna unit.
請求項１乃至請求項１２の何れか１項に記載の信号処理装置と、を備え、
前記相関行列算出部は、前記超音波振動子を用いて得られた受信信号に基づいて、前記相関行列を算出することを特徴とする、水中探知装置。 A transmitting and receiving device provided with an ultrasonic transducer that rotates in a predetermined azimuth direction as a sub-scanning direction and sequentially transmits and receives ultrasonic waves along a distance direction that intersects with the azimuth direction,
A signal processing device according to any one of claims 1 to 12 ,
The underwater detection apparatus, wherein the correlation matrix calculation unit calculates the correlation matrix based on a received signal obtained using the ultrasonic transducer.
前記相関行列に基づいて超解像処理を行うことで、前記副走査方向に沿って得られる前記受信信号の信号強度のスペクトラムを算出する、スペクトラム算出ステップと、
前記受信装置のステアリングベクトルを算出する、ステアリングベクトル算出ステップと、
前記相関行列の固有ベクトルを算出する固有ベクトル算出ステップと、
を含み、
前記スペクトラム算出ステップにおいては、前記固有ベクトル、および、前記ステアリングベクトルに基づいて、前記スペクトラムを算出することを特徴とする、信号処理方法。 Obtained along the sub-scanning direction based on a reception signal obtained by using a receiving device that is displaced along a predetermined sub-scanning direction and sequentially receives signals from the main scanning direction intersecting the sub-scanning direction. Calculating a correlation matrix of the plurality of received signals that are
By performing super-resolution processing based on the correlation matrix, calculating the spectrum of the signal intensity of the received signal obtained along the sub-scanning direction, a spectrum calculation step ,
Calculating a steering vector of the receiving device, a steering vector calculation step,
An eigenvector calculating step of calculating an eigenvector of the correlation matrix,
Only including,
In the spectrum calculation step, the spectrum is calculated based on the eigenvector and the steering vector .
前記相関行列に基づいて超解像処理を行うことで、前記副走査方向に沿って得られる前記受信信号の信号強度のスペクトラムを算出する、スペクトラム算出ステップと、
前記受信装置のステアリングベクトルを算出する、ステアリングベクトル算出ステップと、
前記相関行列の固有ベクトルを算出する固有ベクトル算出ステップと、
をコンピュータに実行させ、
前記スペクトラム算出ステップにおいては、前記固有ベクトル、および、前記ステアリングベクトルに基づいて、前記スペクトラムを算出することを特徴とする、プログラム。 Obtained along the sub-scanning direction based on a reception signal obtained by using a receiving device that is displaced along a predetermined sub-scanning direction and sequentially receives signals from the main scanning direction intersecting the sub-scanning direction. Calculating a correlation matrix of the plurality of received signals that are
By performing super-resolution processing based on the correlation matrix, calculating the spectrum of the signal intensity of the received signal obtained along the sub-scanning direction, a spectrum calculation step ,
Calculating a steering vector of the receiving device, a steering vector calculation step,
An eigenvector calculating step of calculating an eigenvector of the correlation matrix,
To run on your computer ,
A program, wherein in the spectrum calculating step, the spectrum is calculated based on the eigenvector and the steering vector .
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2015245005A JP6706056B2 (en) | 2015-12-16 | 2015-12-16 | Signal processing device, radar device, underwater detection device, signal processing method, and program |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2015245005A JP6706056B2 (en) | 2015-12-16 | 2015-12-16 | Signal processing device, radar device, underwater detection device, signal processing method, and program |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2017110992A JP2017110992A (en) | 2017-06-22 |
| JP6706056B2 true JP6706056B2 (en) | 2020-06-03 |
Family
ID=59081519
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2015245005A Active JP6706056B2 (en) | 2015-12-16 | 2015-12-16 | Signal processing device, radar device, underwater detection device, signal processing method, and program |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP6706056B2 (en) |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2159593B1 (en) * | 2008-08-26 | 2012-05-02 | Nuance Communications, Inc. | Method and device for locating a sound source |
| JP6031268B2 (en) * | 2012-06-21 | 2016-11-24 | 古野電気株式会社 | Radar apparatus, target detection method, and target detection program |
| JP2014115160A (en) * | 2012-12-07 | 2014-06-26 | Furuno Electric Co Ltd | Detection device, detection method, and detection program |
| JP6179973B2 (en) * | 2013-02-15 | 2017-08-16 | 古野電気株式会社 | Signal processing device, underwater detection device, signal processing method, and program |
| WO2015075992A1 (en) * | 2013-11-22 | 2015-05-28 | 古野電気株式会社 | Detection device, underwater detection device, and radar device |
-
2015
- 2015-12-16 JP JP2015245005A patent/JP6706056B2/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| JP2017110992A (en) | 2017-06-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11340342B2 (en) | Automotive radar using 3D printed luneburg lens | |
| KR100749560B1 (en) | Radar apparatus | |
| JP2013015522A (en) | Broad band beam formation device, broad band beam steering device and corresponding method | |
| KR20190084593A (en) | Method and apparatus for estimating direction of arrival using generation of virtual received signals | |
| JP2019152488A (en) | Antenna device and radar device | |
| JP2012108075A (en) | Radar device and target detection method | |
| Blomberg et al. | Adaptive beamforming applied to a cylindrical sonar array using an interpolated array transformation | |
| JP6712313B2 (en) | Signal processing device and radar device | |
| JP6939981B2 (en) | Object detection device and object detection method | |
| EP2317335B1 (en) | Improved beamforming method for analysing signals received by a transducer arrray, and relative detection system | |
| JP5116590B2 (en) | Wave number estimation device | |
| JP6179973B2 (en) | Signal processing device, underwater detection device, signal processing method, and program | |
| EP2963443B1 (en) | Device, method and program for processing signal, underwater detecting device and radar apparatus | |
| JP6706056B2 (en) | Signal processing device, radar device, underwater detection device, signal processing method, and program | |
| JP4481078B2 (en) | Radar equipment | |
| JP7302293B2 (en) | Direction estimation device, direction estimation method, and program | |
| JP2009250952A (en) | Radar system | |
| WO2023003017A1 (en) | Water level detection device | |
| JP2007333404A (en) | Impulse radar | |
| JP6703829B2 (en) | Signal processing device, radar device, underwater detection device, signal processing method, and program | |
| JP4275425B2 (en) | Radar equipment | |
| JP6321472B2 (en) | Signal processing device, underwater detection device, radar device, signal processing method, and signal processing program | |
| JP2005189107A (en) | Radar system | |
| JP4143007B2 (en) | Radar equipment | |
| WO2020160546A1 (en) | Systems and methods for synthetic aperture focusing in ultrasound imaging |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20181206 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20191112 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20191224 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20200218 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20200512 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20200515 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 6706056 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 |