JP2015230284A - Radar apparatus and radar signal processing method of the same - Google Patents

Radar apparatus and radar signal processing method of the same Download PDF

Info

Publication number
JP2015230284A
JP2015230284A JP2014117835A JP2014117835A JP2015230284A JP 2015230284 A JP2015230284 A JP 2015230284A JP 2014117835 A JP2014117835 A JP 2014117835A JP 2014117835 A JP2014117835 A JP 2014117835A JP 2015230284 A JP2015230284 A JP 2015230284A
Authority
JP
Japan
Prior art keywords
range
axis
signal
frequency
doppler
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
Application number
JP2014117835A
Other languages
Japanese (ja)
Other versions
JP6352688B2 (en
Inventor
晋一 竹谷
Shinichi Takeya
晋一 竹谷
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Toshiba Corp
Original Assignee
Toshiba Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Toshiba Corp filed Critical Toshiba Corp
Priority to JP2014117835A priority Critical patent/JP6352688B2/en
Publication of JP2015230284A publication Critical patent/JP2015230284A/en
Application granted granted Critical
Publication of JP6352688B2 publication Critical patent/JP6352688B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

PROBLEM TO BE SOLVED: To provide a pulse compression radar apparatus reducing a loss even at the time of long time integration, suppressing fixed clutter and enhancing target detection capability.SOLUTION: A radar apparatus identifies a substantial range of a target position by inputting a radar pulse reception signal and applying coherent integration processing of a temporal axis to individual range axes of a Σ beam and a Δ beam, focuses on the substantial range of the target position, performs integration processing of coherent integration processing results of the temporal axis in a frequency domain, acquires individual range-Doppler images of the Σ beam and the Δ beam by bringing back to a time domain, detects a target cell by CFAR processing with respect to the range-Doppler image of the Σ beam and performs distance measurement, and performs angle measurement by using a signal of a cell corresponding to the target cell detected by the Σ beam from the range-Doppler image of the Δ beam.

Description

本実施形態は、例えば航空機に搭載して遠距離の小目標を検出するレーダ装置及びそのレーダ信号処理方法に関する。   The present embodiment relates to a radar apparatus that is mounted on, for example, an aircraft and detects a small target at a long distance, and a radar signal processing method thereof.

従来の航空機搭載のレーダ装置では、積分ヒット数の多い長時間積分を行うと、レンジウォーク及び被搭載機の運動による位相の非線形性により積分ロスが生じるという課題があった。また、メインローブが地表面や海面等を向く場合には、固定クラッタを受信してしまうため、目標の検出性能が劣化するという課題があった。   A conventional aircraft-mounted radar device has a problem that integration loss occurs due to the non-linearity of the phase due to the range walk and the motion of the mounted device when long-time integration with a large number of integration hits is performed. Further, when the main lobe faces the ground surface, the sea surface, or the like, the fixed clutter is received, and thus there is a problem that the target detection performance deteriorates.

SAR方式(レンジ圧縮)、大内、“リモートセンシングのための合成開口レーダの基礎”、東京電機大学出版局、pp.131-149(2003)SAR method (range compression), Ouchi, “Basics of Synthetic Aperture Radar for Remote Sensing”, Tokyo Denki University Press, pp.131-149 (2003) SAR方式(Az圧縮)、大内、“リモートセンシングのための合成開口レーダの基礎”、東京電機大学出版局、pp.171-178(2003)SAR method (Az compression), Ouchi, “Basics of Synthetic Aperture Radar for Remote Sensing”, Tokyo Denki University Press, pp.171-178 (2003) CFAR(Constant False Alarm Rate)、吉田、“改定レーダ技術”、電子情報通信学会、pp.87-89(1996)CFAR (Constant False Alarm Rate), Yoshida, “Revised Radar Technology”, IEICE, pp.87-89 (1996) 2次元CFAR、Guy Morris, "AIRBORNE PULSE DOPPLER RADAR 2nd edition", pp.399-404(1996)2D CFAR, Guy Morris, "AIRBORNE PULSE DOPPLER RADAR 2nd edition", pp.399-404 (1996) 位相モノパルス測角、吉田、“改定レーダ技術”、電子情報通信学会、pp.260-264(1996)Phase monopulse angle measurement, Yoshida, “Revised radar technology”, IEICE, pp. 260-264 (1996) MPRF方式、Guy Morris, "AIRBORNE PULSE DOPPLER RADAR 2nd edition", pp.264-270(1996)MPRF, Guy Morris, "AIRBORNE PULSE DOPPLER RADAR 2nd edition", pp.264-270 (1996) HPRF方式、George W.Stimson, "INTRODUCTION TO AIRBORN RADAR", HEGHES, pp231-239 (1983)HPRF system, George W. Stimson, "INTRODUCTION TO AIRBORN RADAR", HEGHES, pp231-239 (1983) FM−CWレーダ、吉田、“改定レーダ技術”、電子情報通信学会、pp.274-275(1996)FM-CW radar, Yoshida, "Revised radar technology", IEICE, pp.274-275 (1996) PGA方式、Charles V.Jakowatz, "Spotlight-Mode Synthetic Aperture Radar: A Signal Processing Approach", Springer, pp.251-256(1996)PGA, Charles V. Jakowatz, "Spotlight-Mode Synthetic Aperture Radar: A Signal Processing Approach", Springer, pp.251-256 (1996)

以上述べたように、従来のレーダ装置では、積分ヒット数の多い長時間積分を行うと、レンジウォーク及び被搭載機の運動による位相の非線形性により積分ロスが生じるという課題があった。また、メインローブが地表面や海面等を向く場合には、固定クラッタを受信してしまうため、目標の検出性能が劣化するという課題があった。   As described above, the conventional radar apparatus has a problem in that when long-time integration with a large number of integration hits is performed, an integration loss occurs due to phase nonlinearity due to the range walk and the motion of the mounted machine. Further, when the main lobe faces the ground surface, the sea surface, or the like, the fixed clutter is received, and thus there is a problem that the target detection performance deteriorates.

本実施形態は上記課題に鑑みなされたもので、長時間積分時にもロスを低減し、固定クラッタを抑圧して目標の検出能力を高めることのできるレーダ装置及びそのレーダ信号処理方法を提供することを目的とする。   The present embodiment has been made in view of the above problems, and provides a radar apparatus and a radar signal processing method thereof that can reduce loss even during long-time integration, suppress fixed clutter, and increase target detection capability. With the goal.

上記の課題を解決するために、本実施形態に係るレーダ装置は、Σビーム画像取得手段と、Δビーム画像取得手段と、測距・測角手段とを具備する。Σビーム画像取得手段は、レーダパルス受信信号を入力してΣビームのレンジ軸に対する時間軸のコヒーレント積分処理を施すことで目標位置の概略範囲を特定し、前記目標位置の概略範囲に焦点を合わせ、前記時間軸のコヒーレント積分処理結果について周波数領域で積分処理して時間領域に戻すことでΣビームのレンジ−ドップラ画像を取得する。Δビーム画像取得手段は、前記レーダパルス受信信号を入力してΔビームのレンジ軸に対する時間軸のコヒーレント積分処理を施すことで目標位置の概略範囲を特定し、前記目標位置の概略範囲に焦点を合わせ、前記時間軸のコヒーレント積分処理結果について周波数領域で積分処理して時間領域に戻すことでΔビームのレンジ−ドップラ画像を取得する。測距・測角手段は、前記Σビーム画像取得手段で取得されたΣビームのレンジ−ドップラ画像に対してCFAR(Constant False Alarm Rate)処理により目標セルを検出して測距し、前記Δビーム画像取得手段で取得されたΔビームのレンジ−ドップラ画像から前記Σビームで検出された目標セルに対応するセルの信号を用いて測角する。   In order to solve the above-described problems, the radar apparatus according to the present embodiment includes a Σ beam image acquisition unit, a Δ beam image acquisition unit, and a distance measurement / angle measurement unit. The Σ-beam image acquisition means inputs a radar pulse reception signal, performs a time-axis coherent integration process with respect to the Σ-beam range axis, identifies an approximate range of the target position, and focuses on the approximate range of the target position. The time-axis coherent integration processing result is integrated in the frequency domain and returned to the time domain to obtain a range-Doppler image of the Σ beam. The Δ beam image acquisition means inputs the radar pulse reception signal, performs a time-axis coherent integration process with respect to the Δ beam range axis, identifies the approximate range of the target position, and focuses on the approximate range of the target position. In addition, a range-Doppler image of Δ beam is acquired by integrating the time-axis coherent integration processing result in the frequency domain and returning it to the time domain. The distance measuring / angle measuring means detects a target cell by CFAR (Constant False Alarm Rate) processing for the Σ beam range-Doppler image acquired by the Σ beam image acquiring means, and measures the distance. The angle is measured using the signal of the cell corresponding to the target cell detected by the Σ beam from the range-Doppler image of the Δ beam acquired by the image acquisition means.

第1の実施形態に係るレーダ装置の構成を示すブロック図。1 is a block diagram showing a configuration of a radar apparatus according to a first embodiment. 図1に示すレーダ装置において、信号処理部の処理の流れを示すフローチャート。2 is a flowchart showing a flow of processing of a signal processing unit in the radar apparatus shown in FIG. 図1に示すレーダ装置において、機体フライト軸と目標の位置関係を示す概念図。The conceptual diagram which shows the positional relationship of a body flight axis | shaft and a target in the radar apparatus shown in FIG. 第2の実施形態に係るレーダ装置の構成を示すブロック図。The block diagram which shows the structure of the radar apparatus which concerns on 2nd Embodiment. 図4に示すレーダ装置において、信号処理部の処理の流れを示すフローチャート。5 is a flowchart showing a processing flow of a signal processing unit in the radar apparatus shown in FIG. 図4に示すレーダ装置において、信号処理部の位相補正を説明するための概念図。The conceptual diagram for demonstrating the phase correction of a signal processing part in the radar apparatus shown in FIG. 図4に示すレーダ装置において、信号処理部の位相補正の処理の流れを示すフローチャート。5 is a flowchart showing a flow of phase correction processing of a signal processing unit in the radar apparatus shown in FIG. 4. 第3の実施形態に係るレーダ装置の構成を示すブロック図。The block diagram which shows the structure of the radar apparatus which concerns on 3rd Embodiment. 図8に示すレーダ装置において、信号処理部の処理の流れを示すフローチャート。9 is a flowchart showing a flow of processing of a signal processing unit in the radar apparatus shown in FIG. レーダ装置におけるLPRFの区分とアンビギュイティとの関係を示す図。The figure which shows the relationship between the division of LPRF in a radar apparatus, and ambiguity. レーダ装置におけるHPRFの区分とアンビギュイティとの関係を示す図。The figure which shows the relationship between the division of HPRF in a radar apparatus, and ambiguity. レーダ装置におけるMPRFの区分とアンビギュイティとの関係を示す図。The figure which shows the relationship between the division of MPRF in a radar apparatus, and ambiguity. 図8に示すレーダ装置において、周波数フィルタの役割を示す図。The figure which shows the role of a frequency filter in the radar apparatus shown in FIG. 図8に示すレーダ装置において、2次元CFARの説明を示す図。The figure which shows description of two-dimensional CFAR in the radar apparatus shown in FIG. 図8に示すレーダ装置において、測距の処理を説明するための図。The figure for demonstrating the ranging process in the radar apparatus shown in FIG. 図8に示すレーダ装置において、測速の処理を説明するための図。The figure for demonstrating the speed measurement process in the radar apparatus shown in FIG. 第4の実施形態に係るレーダ装置の構成を示すブロック図。The block diagram which shows the structure of the radar apparatus which concerns on 4th Embodiment. 図15に示す装置において、信号処理部の処理の流れを示すフローチャート。The flowchart which shows the flow of a process of a signal processing part in the apparatus shown in FIG. 図15に示す装置において、HPRFの測距の説明を示す図。The figure which shows description of the ranging of HPRF in the apparatus shown in FIG.

以下、実施形態について、図面を参照して説明する。   Hereinafter, embodiments will be described with reference to the drawings.

(第1の実施形態)
以下、図1乃至図3を参照して、第1の実施形態に係るレーダ装置を説明する。
(First embodiment)
The radar apparatus according to the first embodiment will be described below with reference to FIGS. 1 to 3.

図1は上記レーダ装置の系統構成を示すブロック図、図2はその具体的な処理の流れを示すフローチャート、図3は上記レーダ装置が搭載される航空機の機体フライト軸と目標との位置関係を示す概念図である。   FIG. 1 is a block diagram showing a system configuration of the radar apparatus, FIG. 2 is a flowchart showing a specific processing flow thereof, and FIG. 3 shows a positional relationship between a fuselage flight axis of an aircraft on which the radar apparatus is mounted and a target. FIG.

図1において、アンテナ1は複数のアンテナ素子を配列して大開口アレイを形成してなるフェーズドアレイアンテナであり、送受信器2の送受信部21から繰り返し供給される特定周波数の送信パルス信号(以下、PRF(Pulse Repetition Frequency)信号)を指定方向に送出してその反射波を受信する。送受信部21はアンテナ1の複数のアンテナ素子でそれぞれ受信された信号をビーム制御部22からの指示に従って位相制御を施し合成することで、任意の方向に受信ビームを形成してPRF受信信号を得る(図2:ステップS11)。受信ビームで得られたPRF受信信号は信号処理器3に送られる。   In FIG. 1, an antenna 1 is a phased array antenna in which a plurality of antenna elements are arranged to form a large aperture array. A PRF (Pulse Repetition Frequency) signal) is transmitted in a designated direction and the reflected wave is received. The transmission / reception unit 21 forms a reception beam in an arbitrary direction to obtain a PRF reception signal by performing phase control and combining the signals respectively received by the plurality of antenna elements of the antenna 1 according to an instruction from the beam control unit 22 (FIG. 2: Step S11). The PRF reception signal obtained by the reception beam is sent to the signal processor 3.

この信号処理器3は、送受信器2からのPRF受信信号をΣビーム用とΔビーム用の2系統に分配する。Σビーム用の系統は、AD(Analog-Digital)変換部311、レンジ軸FFT(Fast Fourier Transformation)部(FFTx)321、レンジ参照信号乗算部331、データ保存部341、クロスレンジ参照信号乗算部351、クロスレンジ軸FFT部(FFTy)361、レンジ軸逆FFT部(逆FFTx)371、CFAR処理部38、測距部39及び測角部3Aを備える。Δビーム用の系統は、AD(Analog-Digital)変換部312、レンジ軸FFT部(FFTx)322、レンジ参照信号乗算部332、データ保存部342、クロスレンジ参照信号乗算部352、クロスレンジ軸FFT部(FFTy)362、レンジ軸逆FFT部(逆FFTx)372、セル抽出部3Bを備える。   The signal processor 3 distributes the PRF received signal from the transceiver 2 into two systems for Σ beam and Δ beam. The Σ beam system includes an AD (Analog-Digital) conversion unit 311, a range axis FFT (Fast Fourier Transformation) unit (FFTx) 321, a range reference signal multiplication unit 331, a data storage unit 341, and a cross range reference signal multiplication unit 351. , A cross range axis FFT unit (FFTy) 361, a range axis inverse FFT unit (inverse FFTx) 371, a CFAR processing unit 38, a distance measuring unit 39, and an angle measuring unit 3A. The Δ beam system includes an AD (Analog-Digital) converter 312, a range axis FFT unit (FFTx) 322, a range reference signal multiplier 332, a data storage unit 342, a cross range reference signal multiplier 352, and a cross range axis FFT. Unit (FFTy) 362, a range axis inverse FFT unit (inverse FFTx) 372, and a cell extraction unit 3B.

上記Σビーム用の系統において、送受信器2から供給されるPRF受信信号はAD変換部311によってディジタル信号に変換される(図2:ステップS12)。その変換結果はレンジ軸FFT部321にてレンジ軸方向にFFT処理され(図2:ステップS13)、レンジ参照信号乗算部331にてレンジ軸方向の線形チャープによるレンジ参照信号が乗算されてレンジ圧縮される(図2:ステップS14)。レンジ圧縮されたPRF受信信号はデータ保存部341に保存され(図2:ステップS15)、長時間のコヒーレント積分に供される。   In the Σ beam system, the PRF reception signal supplied from the transceiver 2 is converted into a digital signal by the AD converter 311 (FIG. 2: step S12). The conversion result is FFT processed in the range axis direction by the range axis FFT unit 321 (FIG. 2: step S13), and the range reference signal multiplication unit 331 multiplies the range reference signal by the linear chirp in the range axis direction to compress the range. (FIG. 2: Step S14). The range-compressed PRF reception signal is stored in the data storage unit 341 (FIG. 2: step S15) and used for long-time coherent integration.

また、上記レンジ圧縮されたPRF受信信号は、処理選択部3C1にて処理選択を受けた場合に、クロスレンジ軸FFT部361にてクロスレンジ軸(PRI軸)方向にFFT処理され(図2:ステップS16)、さらにレンジ軸逆FFT部371にてレンジ軸方向に逆FFT処理される(図2:ステップS17)。このようにしてパルス圧縮された信号に対して、CFAR処理部38によりレンジ−ドップラ軸における目標セルの検出が行われ(図2:ステップS18)、測距部39及び測角部3Aにてそれぞれ検出目標セルの距離、レンジ方向の角度が検出される(図2:ステップS19,S20)。   The PRF received signal subjected to the range compression is subjected to FFT processing in the cross range axis (PRI axis) direction by the cross range axis FFT unit 361 when the processing selection unit 3C1 receives a processing selection (FIG. 2: In step S16), the range axis inverse FFT unit 371 performs inverse FFT processing in the range axis direction (FIG. 2: step S17). The target cell in the range-Doppler axis is detected by the CFAR processing unit 38 for the pulse-compressed signal in this way (FIG. 2: step S18), and the distance measuring unit 39 and the angle measuring unit 3A respectively. The distance of the detection target cell and the angle in the range direction are detected (FIG. 2: steps S19 and S20).

一方、上記Δビーム用の系統において、送受信器2から供給されるPRF受信信号はAD変換部312によってディジタル信号に変換される(図2:ステップS21)。その変換結果はレンジ軸FFT部322にてレンジ軸方向にFFT処理され(図2:ステップS22)、レンジ参照信号乗算部332にてレンジ軸方向の線形チャープによる参照信号が乗算されてレンジ圧縮される(図2:ステップS23)。レンジ圧縮されたPRF受信信号はデータ保存部342に保存され(図2:ステップS24)、長時間のコヒーレント積分に供される。   On the other hand, in the Δ beam system, the PRF reception signal supplied from the transceiver 2 is converted into a digital signal by the AD converter 312 (FIG. 2: step S21). The conversion result is FFT processed in the range axis direction by the range axis FFT unit 322 (FIG. 2: step S22), and the range reference signal multiplication unit 332 multiplies the reference signal by the linear chirp in the range axis direction to perform range compression. (FIG. 2: Step S23). The range-compressed PRF reception signal is stored in the data storage unit 342 (FIG. 2: step S24) and used for long-time coherent integration.

また、上記レンジ圧縮されたPRF受信信号は、処理選択部3C2にて処理選択を受けた場合に、クロスレンジ軸FFT部362にてクロスレンジ軸(PRI軸)方向にFFT処理され(図2:ステップS25)、さらにレンジ軸逆FFT部372にてレンジ軸方向に逆FFT処理される(図2:ステップS26)。このようにしてパルス圧縮された信号に対して、セル抽出部3Bにて、CFAR処理部38からの情報に基づいてΣビーム側で検出したレンジ−ドップラ軸の目標セルと同様のΔビームの目標セルが抽出される(図2:ステップS27)。抽出されたΔビームの目標セルについては、測角部3Aにてクロスレンジ方向の角度が検出される。   Further, the PRF received signal subjected to the range compression is subjected to FFT processing in the cross range axis (PRI axis) direction by the cross range axis FFT unit 362 when the processing selection unit 3C2 receives the processing selection (FIG. 2: In step S25), the range axis inverse FFT unit 372 performs inverse FFT processing in the range axis direction (FIG. 2: step S26). With respect to the pulse-compressed signal in this way, the cell extraction unit 3B detects the Δ beam target similar to the target cell of the range-Doppler axis detected on the Σ beam side based on the information from the CFAR processing unit 38. A cell is extracted (FIG. 2: step S27). For the extracted target cell of the Δ beam, the angle measuring unit 3A detects the angle in the cross range direction.

ここで、上記処理選択部3C1,3C2は、上記測角部3Aにおける測角の検出対象の有無について判断し(図2:ステップS28)、検出対象が得られなかった場合、次のPRI送受信で得られたPRF受信信号について目標セルの検出処理を行い、検出対象が得られた場合には、目標のRCS(レーダ散乱断面積)が小さいことを考慮して、更に高効率な積分効果を得るために、長時間のコヒーレント積分を行う。このために、処理選択部3C1,3C2は、測角部3Aにて目標が検出されたビーム方向に対して、データ保存部341,342に保存していたレンジ圧縮データを読み出して(経過時間によっては再度送受信した信号を用いてもよい)長時間のコヒーレント積分処理に切り替える。   Here, the process selection units 3C1 and 3C2 determine whether or not there is an angle detection target in the angle measurement unit 3A (FIG. 2: step S28), and if no detection target is obtained, the next PRI transmission / reception is performed. A target cell detection process is performed on the obtained PRF received signal, and when a detection target is obtained, an even more efficient integration effect is obtained in consideration of a small target RCS (radar scattering cross section). For this purpose, long-time coherent integration is performed. For this purpose, the process selection units 3C1 and 3C2 read out the range compression data stored in the data storage units 341 and 342 with respect to the beam direction in which the target is detected by the angle measurement unit 3A (depending on the elapsed time). May use a signal transmitted and received again). Switch to a long-time coherent integration process.

具体的には、Σビーム用の系統において、長時間のコヒーレント積分処理に切り替えられると、データ保存部341からレンジ圧縮データが読み出され(図2:ステップS29)、クロスレンジ参照信号乗算部351にてクロスレンジ軸方向の線形チャープによるクロスレンジ参照信号が乗算された後(図2:ステップS30)、上記クロスレンジ軸FFT部361によるクロスレンジ軸(PRI軸)方向のFFT処理(図2:ステップS31)、レンジ軸逆FFT部371によるレンジ軸方向の逆FFT処理(図2:ステップS32)、CFAR処理部38によるレンジ−ドップラ軸における目標セルの検出(図2:ステップS33)、測距部39及び測角部3Aによる検出目標セルの距離、レンジ方向の角度の検出が行われる(図2:ステップS34,S35)。   Specifically, in the Σ beam system, when the long-time coherent integration process is switched, the range compressed data is read from the data storage unit 341 (FIG. 2: step S29), and the cross range reference signal multiplication unit 351 is read out. 2 is multiplied by a cross range reference signal by linear chirp in the cross range axis direction (FIG. 2: step S30), and then the FFT processing in the cross range axis (PRI axis) direction by the cross range axis FFT unit 361 (FIG. 2: Step S31), reverse FFT processing in the range axis direction by the range axis inverse FFT unit 371 (FIG. 2: step S32), target cell detection in the range-Doppler axis by the CFAR processing unit 38 (FIG. 2: step S33), and distance measurement The distance of the detection target cell and the angle in the range direction are detected by the unit 39 and the angle measurement unit 3A (FIG. 2: scan). -Up S34, S35).

一方、Δビーム用の系統においても、長時間のコヒーレント積分処理に切り替えられると、データ保存部342からレンジ圧縮データが読み出され(図2:ステップS36)、クロスレンジ参照信号乗算部352にてクロスレンジ軸方向の線形チャープによるクロスレンジ参照信号が乗算された後(図2:ステップS37)、クロスレンジ軸FFT部362によるクロスレンジ軸(PRI軸)方向のFFT処理(図2:ステップS38)、レンジ軸逆FFT部372によるレンジ軸方向の逆FFT処理(図2:ステップS39)が行われる。続いて、セル抽出部3Bにて、CFAR処理部38からの情報に基づいてΣビーム側で検出したレンジ−ドップラ軸の目標セルと同様の目標セルが抽出される(図2:ステップS40)。抽出されたΔビームの目標セルについては、測角部3Aにてクロスレンジ方向の角度が検出される。このようにして測角部3Aで目標セルのレンジ方向、クロスレンジ方向の測角結果が得られる。   On the other hand, also in the Δ beam system, when switching to the long-time coherent integration process, the range compressed data is read from the data storage unit 342 (FIG. 2: step S36), and the cross range reference signal multiplication unit 352 After the cross range reference signal by the linear chirp in the cross range axis direction is multiplied (FIG. 2: step S37), the cross range axis FFT unit 362 performs FFT processing in the cross range axis (PRI axis) direction (FIG. 2: step S38). Then, the inverse FFT process in the range axis direction (FIG. 2: step S39) is performed by the range axis inverse FFT unit 372. Subsequently, the cell extraction unit 3B extracts a target cell similar to the target cell of the range-Doppler axis detected on the Σ beam side based on the information from the CFAR processing unit 38 (FIG. 2: step S40). For the extracted target cell of the Δ beam, the angle measuring unit 3A detects the angle in the cross range direction. In this way, the angle measurement unit 3A can obtain angle measurement results in the range direction and cross range direction of the target cell.

ここで、図3は、上記構成によるレーダ装置が、飛翔体に搭載された場合に、実開口ビームを形成して目標検出範囲に常に照射するようにビームを向ける場合を示している。この際、メインロ−ブやサイドローブから地表面のクラッタも送受信するので、その座標系も合わせて記述している。合成開口時間(1サイクル)内でPRI(Pulse Repetition Interval)間隔で送信したパルス毎に、PRI内のレンジセル単位でデータを取得する。この取得データを用いて長時間のコヒーレント積分処理を実施する。   Here, FIG. 3 shows a case where the radar apparatus having the above configuration forms a real aperture beam and directs the beam so as to always irradiate the target detection range when mounted on a flying object. At this time, since clutter on the ground surface is also transmitted and received from the main lobe and side lobes, the coordinate system is also described. Data is acquired in units of range cells in the PRI for each pulse transmitted at a PRI (Pulse Repetition Interval) interval within the synthetic aperture time (one cycle). A long-time coherent integration process is performed using the acquired data.

次に、上記の長時間のコヒーレント積分処理について詳述する。
まずレンジ圧縮について述べる(非特許文献1参照)。レンジ圧縮は、入力信号とレンジ圧縮用の参照信号との相関処理であり、これを周波数領域で行う場合について定式化すると次の通りである(図1の321,322,331,332)。

Figure 2015230284
Next, the long-time coherent integration process will be described in detail.
First, range compression will be described (see Non-Patent Document 1). Range compression is a correlation process between an input signal and a reference signal for range compression, and is formulated as follows in the frequency domain (321, 322, 331, 332 in FIG. 1).
Figure 2015230284

(3)式で示すレンジ圧縮信号sを時間軸上の信号にするには逆フーリエ変換すればよい。但し、この後でクロスレンジ圧縮(クロスレンジ参照信号乗算部351,352、非特許文献2参照)を行うために、レンジ圧縮信号sは(ω、u)軸のままとする。   In order to convert the range compressed signal s represented by the expression (3) into a signal on the time axis, inverse Fourier transform may be performed. However, in order to perform cross range compression (cross range reference signal multipliers 351 and 352, non-patent document 2) after this, the range compressed signal s remains on the (ω, u) axis.

次にクロスレンジ圧縮用の参照信号fs0を生成する。

Figure 2015230284
Next, a reference signal fs0 for cross range compression is generated.
Figure 2015230284

尚、検出位置(Xc, Yc, Zc)は、測距部39及び測角部3Aで検出した測距値、測角値より算出した位置である。目標が複数の場合は、検出目標毎に参照信号を生成する。また、オフセットは、検出範囲について、追跡情報等を用いて補正する場合等に用いる。前述のレンジ圧縮信号sとクロスレンジ圧縮用の参照信号fs0を乗算して信号csを得る。

Figure 2015230284
The detection positions (Xc, Yc, Zc) are positions calculated from the distance measurement values and the angle measurement values detected by the distance measurement unit 39 and the angle measurement unit 3A. When there are a plurality of targets, a reference signal is generated for each detection target. The offset is used when the detection range is corrected using tracking information or the like. A signal cs is obtained by multiplying the above-described range compression signal s and the reference signal fs0 for cross range compression.
Figure 2015230284

これを用いて、u軸でFFTして信号fcs(ω, ku)を得る(図1、クロスレンジ軸FFT部361,362)。

Figure 2015230284
Using this, FFT is performed on the u-axis to obtain a signal fcs (ω, ku) (FIG. 1, cross-range axis FFT units 361 and 362).
Figure 2015230284

長時間積分の出力fpは、fcsのω軸に関するレンジ軸逆FFT部371,372により算出できる。

Figure 2015230284
The long-time integration output fp can be calculated by the range axis inverse FFT units 371 and 372 regarding the ω axis of fcs.
Figure 2015230284

fp(t,ku)の出力より、1次元または2次元CFAR処理(非特許文献3、4参照)を用いて検出し、測距部39にて時間をレンジ(距離)に変換して出力する。 From the output of fp (t, ku), it is detected using one-dimensional or two-dimensional CFAR processing (see Non-Patent Documents 3 and 4), and the distance measurement unit 39 converts the time into a range (distance) and outputs it. .

本実施形態は、(Xc, Yc, Zc)に焦点を合わせた処理である。目標が十分遠距離の場合は、焦点を無限遠点に設定してもよく、処理規模を削減することができる。また、位相ずれによる影響を軽減することができる。   In the present embodiment, the processing focuses on (Xc, Yc, Zc). When the target is a sufficiently long distance, the focal point may be set at an infinite point, and the processing scale can be reduced. In addition, the influence of the phase shift can be reduced.

次に測角する手法について述べる。アンテナのΔビーム(ΔAz, ΔEL)の出力について、Σビームと同様にレンジ−ドップラ画像を得る。次にΣビームで検出したセル(レンジ−ドップラ軸のセル)と同様のセルを抽出し、次式の誤差電圧を算出する。

Figure 2015230284
Next, a method for measuring the angle will be described. For the output of the Δ beam (ΔAz, ΔEL) of the antenna, a range-Doppler image is obtained in the same manner as the Σ beam. Next, a cell similar to the cell (range-Doppler axis cell) detected by the Σ beam is extracted, and an error voltage of the following equation is calculated.
Figure 2015230284

ΣとΔのビームパターンにおいて、予め取得した角度特性により、基準となる誤差電圧εrefと角度の関係をテーブル化しておき、εを観測し、テーブルを引用して角度を抽出すれば、観測角度が得られる。Az角、EL角について同様の処理を行うと、両者の軸の測角θAz、θELを実施できる(非特許文献5参照)。   In the beam patterns of Σ and Δ, the relationship between the reference error voltage εref and the angle is tabulated based on the angle characteristics acquired in advance, ε is observed, and if the angle is extracted by quoting the table, the observation angle becomes can get. If the same processing is performed for the Az angle and the EL angle, the measured angles θAz and θEL of both axes can be implemented (see Non-Patent Document 5).

以上のように、第1の実施形態の構成によれば、レンジウォークによるロス及び位相の非線形成を補正して周波数軸で積分するので、ロスを低減することができる。さらにΔビームを用いているので、測角を実施することができる。   As described above, according to the configuration of the first embodiment, the loss and phase nonlinearity due to the range walk are corrected and integrated on the frequency axis, so that the loss can be reduced. Further, since the Δ beam is used, angle measurement can be performed.

(第2の実施形態)
長時間のコヒーレント積分による手法では、その積分の間に、目標の加速度や、被搭載機の加速度による位相変化が生じる場合がある。この対策手法については、合成開口処理のオートフォーカス手法であるPGA(Phase gradient autofocus:非特許文献9)と類似の方法を適用することができる。図4は、第2の実施形態として、この対処手法を適用した場合のレーダ装置の系統構成を示すブロック図、図5はその処理方法を適用した場合の具体的な処理の流れを示すフローチャート、図6は上記対処方法の処理を説明するための概念図である。尚、図4及び図5において、図1及び図2と同一部分には同一符号を付して示す。
(Second Embodiment)
In the method using long-time coherent integration, there may be a phase change due to the target acceleration or the acceleration of the mounted machine during the integration. As this countermeasure method, a method similar to PGA (Phase gradient autofocus: Non-Patent Document 9), which is an autofocus method for synthetic aperture processing, can be applied. FIG. 4 is a block diagram showing a system configuration of a radar apparatus when this coping method is applied as a second embodiment, and FIG. 5 is a flowchart showing a specific processing flow when the processing method is applied. FIG. 6 is a conceptual diagram for explaining the processing of the coping method. 4 and 5, the same parts as those in FIGS. 1 and 2 are denoted by the same reference numerals.

本実施形態において、第1の実施形態と異なる点は、Σ系統のレンジ軸逆FFT処理部371で得られたレンジ−ドップラ(クロスレンジ)画像から補正信号を抽出する補正信号抽出部3Dと、上記補正信号に基づいてクロスレンジ参照信号乗算部351,352のクロスレンジ参照信号を補正するクロスレンジ参照信号補正部3Eを追加した点にある(図5のステップS41,S42)。   In this embodiment, the difference from the first embodiment is that a correction signal extraction unit 3D that extracts a correction signal from a range-Doppler (cross range) image obtained by the range axis inverse FFT processing unit 371 of the Σ system, A cross range reference signal correction unit 3E that corrects the cross range reference signal of the cross range reference signal multipliers 351 and 352 based on the correction signal is added (steps S41 and S42 in FIG. 5).

まず、図6(a)に示すように、レンジ−ドップラ(クロスレンジ)画像の結果より、所定の振幅スレショルドを超えたN個の極大値を抽出する。次に、極大値のドップラ軸に対する位相勾配を除き、飛翔経路及び機体動揺による位相ずれのみを抽出するため、極大値をドップラ軸で0周波数にシフトする。ここで、位相ずれの振動成分を取り除き、安定した補正成分を得るために、図6(b)に示すように、その周囲に窓関数(ウィンドウ)を乗算し、窓関数の外側にゼロ埋めした信号s0を生成し、この信号を逆FFTすることで、図6(c)に示す誤差信号S0を求める。

Figure 2015230284
First, as shown in FIG. 6A, N local maximum values exceeding a predetermined amplitude threshold are extracted from the result of the range-Doppler (cross range) image. Next, the maximum value is shifted to 0 frequency on the Doppler axis in order to extract only the phase shift due to the flight path and body motion, excluding the phase gradient with respect to the Doppler axis of the maximum value. Here, in order to remove the vibration component of the phase shift and obtain a stable correction component, as shown in FIG. 6B, a window function (window) is multiplied around it and zero-filled outside the window function. An error signal S0 shown in FIG. 6C is obtained by generating a signal s0 and performing inverse FFT on the signal.
Figure 2015230284

続いて、図6(d)に示すように、誤差信号S0の逆特性となる補正量Wc(t)を求め、クロスレンジ圧縮信号fs0(ω、u)(レンジの周波数軸、飛行長)を時間軸(PRI軸)に変換したfs0(t)に乗算する。このようにして、図6(e)に示す補正後の参照信号Rcを求める。

Figure 2015230284
Subsequently, as shown in FIG. 6 (d), a correction amount Wc (t) that is an inverse characteristic of the error signal S0 is obtained, and the cross-range compression signal fs0 (ω, u) (range frequency axis, flight length) is obtained. Multiply fs0 (t) converted to the time axis (PRI axis). In this way, the corrected reference signal Rc shown in FIG. 6E is obtained.
Figure 2015230284

この参照信号Fs0を用いて、(5)〜(7)式の処理により、補正後のレンジ−ドップラ画像が得られる。この処理フローを図7に示す。   Using this reference signal Fs0, a corrected range-Doppler image is obtained by the processing of equations (5) to (7). This processing flow is shown in FIG.

すなわち、図7において、SAR処理が指示されると(ステップS51)、クロスレンジ圧縮信号のピーク値を抽出し(ステップS52)、そのピークの前後±Mセルを抽出して(ステップS53)、ゼロシフトを行った後(ステップS54)、ウィンドウ乗算を行ってゼロ埋めを行う(ステップS55)。次に、逆FFTにより周波数軸の信号から時間軸の信号に変換して誤差信号を求め(ステップS56)、その逆特性となる補正値を算出して(ステップS57)、その補正値を時間軸に変換されたクロスレンジ圧縮信号に乗算してSAR処理する(ステップS58)。この時点で検出されたピーク値全てについてSAR処理を行ったか判断し(ステップS59)、未処理のピーク値がある場合にはピーク値を変更してステップS53より処理を繰り返す(ステップS60)。全てのピーク値について処理が完了した場合には、一連の処理を終了する。この一連の処理を必要に応じてM回(M≧1)繰り返すことにより、積分ロスの少ないレンジ−ドップラ画像が得られる。このレンジ−ドップラ画像を用いてCFAR処理することにより、目標セルを確実に検出することができ、より正確な距離を出力することができる。   That is, in FIG. 7, when SAR processing is instructed (step S51), the peak value of the cross-range compressed signal is extracted (step S52), ± M cells before and after the peak are extracted (step S53), and zero shift is performed. (Step S54), window multiplication is performed to perform zero padding (step S55). Next, an error signal is obtained by converting from a frequency axis signal to a time axis signal by inverse FFT (step S56), a correction value having the inverse characteristic is calculated (step S57), and the correction value is converted to the time axis. SAR processing is performed by multiplying the cross-range compressed signal converted into (step S58). It is determined whether the SAR process has been performed for all the peak values detected at this time (step S59). If there is an unprocessed peak value, the peak value is changed and the process is repeated from step S53 (step S60). When the processing is completed for all peak values, a series of processing ends. By repeating this series of processing M times (M ≧ 1) as necessary, a range-Doppler image with little integration loss can be obtained. By performing CFAR processing using this range-Doppler image, the target cell can be reliably detected, and a more accurate distance can be output.

以上のように、第2の実施形態の構成によれば、レンジ−ドップラ画像から抽出した信号を用いて補正信号を生成し、位相補正するので、積分ロスを低減することができる。   As described above, according to the configuration of the second embodiment, the correction signal is generated using the signal extracted from the range-Doppler image and the phase is corrected, so that the integral loss can be reduced.

(第3の実施形態)
図8は、第3の実施形態として、MPRFを適用した場合のレーダ装置の系統構成を示すブロック図、図9はその処理方法を適用した場合の具体的な処理の流れを示すフローチャートである。尚、図8及び図9において、図4及び図5と同一部分には同一符号を付して示す。
(Third embodiment)
FIG. 8 is a block diagram showing a system configuration of a radar apparatus when MPRF is applied as a third embodiment, and FIG. 9 is a flowchart showing a specific processing flow when the processing method is applied. 8 and 9, the same parts as those in FIGS. 4 and 5 are denoted by the same reference numerals.

本実施形態において、第2の実施形態と異なる点は、Σ系統、Δ系統それぞれのクロスレンジ軸FFT処理部361,362で得られた信号をそれぞれ周波数フィルタ3F1,3F2に通すことで、固定クラッタが存在する特定の周波数成分を取り出してクロスレンジ逆FFT処理部371,372に送るようにした点(図9のステップS60,S61,S65,S66)と、Σビーム側のCFAR処理により得られるレンジ−ドップラ軸の目標セルをMPRF測距・測速部3Gに入力して目標セルに対する測距及び測速演算を実行するようにした点にある(図9のステップS64,S69)。   In this embodiment, the difference from the second embodiment is that the signals obtained by the cross-range axis FFT processing units 361 and 362 of the Σ system and Δ system are passed through the frequency filters 3F1 and 3F2, respectively. And a range obtained by CFAR processing on the side of the Σ beam and a point where a specific frequency component in which is present is extracted and sent to the cross-range inverse FFT processing units 371 and 372 (steps S60, S61, S65, and S66 in FIG. 9). The target cell of the Doppler axis is input to the MPRF ranging / speed measuring unit 3G to perform distance measurement and speed measurement calculation on the target cell (steps S64 and S69 in FIG. 9).

ここで、レーダ装置におけるPRFとレンジ−ドップラの関係には、LPRF(Low PRF)、HPRF(High PRF)、MPRF(Middle PRF)の3種類がある。その概要を図10〜図12に示す。LPRFは、図10に示すように、レンジ軸にアンビギュイティが無く、ドップラ周波数軸にアンビギュイティがある。HPRFは、図11に示すように、レンジ軸にアンビギュイティがあり、ドップラ周波数軸にアンビギュイティが無い。MPRFは、図12に示すように、レンジ軸及びドップラ周波数軸の両者にアンビギュイティがある。   Here, there are three types of relationship between PRF and range-Doppler in the radar apparatus: LPRF (Low PRF), HPRF (High PRF), and MPRF (Middle PRF). The outline is shown in FIGS. As shown in FIG. 10, LPRF has no ambiguity on the range axis and ambiguity on the Doppler frequency axis. As shown in FIG. 11, HPRF has ambiguity on the range axis and no ambiguity on the Doppler frequency axis. As shown in FIG. 12, the MPRF has ambiguity on both the range axis and the Doppler frequency axis.

レンジ−ドップラ画像には、固定クラッタによる画像も含まれており、強度が強い場合には、移動目標のみを抽出するために抑圧する必要がある。このためのクラッタ周波数は、図3を参照すれば、メインビームの左右の端において、次式で表現できる。

Figure 2015230284
The range-Doppler image includes an image by a fixed clutter. When the intensity is strong, it is necessary to suppress in order to extract only the moving target. The clutter frequency for this purpose can be expressed by the following equation at the left and right ends of the main beam with reference to FIG.
Figure 2015230284

具体的には、レンジ−ドップラ画像において、図13に示すように固定クラッタが含まれる特定周波数を抑圧する周波数フィルタを用いて、この周波数範囲の信号を抑圧した上で、図14に示す2次元CFARの検出を行う。このCFAR処理は、時間軸、周波数軸に分けて行う1次元CFARを適用してもよいのは言うまでもない。   Specifically, in the range-Doppler image, a signal in this frequency range is suppressed using a frequency filter that suppresses a specific frequency including a fixed clutter as shown in FIG. CFAR detection is performed. Needless to say, this CFAR process may be applied to a one-dimensional CFAR that is divided into a time axis and a frequency axis.

本実施形態は、MPRFを適用した場合のレーダ装置である。すなわち、被搭載機の速度が高く、走査角度が大きい場合には、クラッタ周波数範囲が広がる。クラッタフリー領域を広げるためには、PRFをあげて、距離及び速度ともに折り返しが発生(アンビギュイティ)するMPRFを適用する。MPRFの場合には、周波数プラインド、パルスブラインドを避けて、測距/測速を行うために、N個のPRFにより送受信を行う。   The present embodiment is a radar apparatus when MPRF is applied. That is, when the speed of the mounted machine is high and the scanning angle is large, the clutter frequency range is widened. In order to widen the clutter-free region, PRF is raised and MPRF in which folding (ambiguity) occurs in both distance and speed is applied. In the case of MPRF, transmission / reception is performed by N PRFs in order to perform distance measurement / speed measurement while avoiding frequency-principles and pulse blinds.

MPRFにおける測距/測速は、図15に示すレンジ軸及び図16に示すドップラ軸において、複数の目標A,Bについて複数のPRI1〜PRI3でそれぞれレンジが合致する、あるいはドップラが合致するM/N個の検出があがるレンジ及びドップラ周波数を抽出することにより行う(非特許文献6参照)。ドップラ周波数fdと速度Vの関係は次の通りである。

Figure 2015230284
The range / speed measurement in MPRF is the M / N in which the ranges match in the plurality of PRI1 to PRI3 for the plurality of targets A and B on the range axis shown in FIG. 15 and the Doppler axis shown in FIG. The detection is performed by extracting the range where the detection is performed and the Doppler frequency (see Non-Patent Document 6). The relationship between the Doppler frequency fd and the velocity V is as follows.
Figure 2015230284

以上のように、第3の実施形態の構成によれば、MPRFを適用することにより、クラッタフリー領域を広げることができ、目標を確実に検出し、測距/測速することができる。   As described above, according to the configuration of the third embodiment, by applying MPRF, it is possible to widen the clutter-free region, reliably detect the target, and perform distance measurement / speed measurement.

(第4の実施形態)
図17は、第4の実施形態として、HPRFを適用した場合のレーダ装置の系統構成を示すブロック図、図18はその処理方法を適用した場合の具体的な処理の流れを示すフローチャートである。尚、図17及び図18において、図8及び図9と同一部分には同一符号を付して示す。
(Fourth embodiment)
FIG. 17 is a block diagram showing a system configuration of a radar apparatus when HPRF is applied as the fourth embodiment, and FIG. 18 is a flowchart showing a specific processing flow when the processing method is applied. 17 and 18, the same parts as those in FIGS. 8 and 9 are denoted by the same reference numerals.

本実施形態において、第3の実施形態と異なる点は、Σビーム側のCFAR処理により得られるレンジ−ドップラ軸の目標セルをHPRF測距・測速部3Hに入力してHPRFによって目標セルに対する測距及び測速演算を実行するようにした点にある(図18のステップS72,S75)。   This embodiment is different from the third embodiment in that a target cell of the range-Doppler axis obtained by the CFAR processing on the Σ beam side is input to the HPRF ranging / speed measuring unit 3H, and ranging with respect to the target cell is performed by HPRF. In addition, the speed measurement calculation is performed (steps S72 and S75 in FIG. 18).

すなわち、被搭載機と目標の相対速度が高い場合には、ドップラ周波数軸でクラッタフリー領域を広げるために、HPRFによる送受信を適用すればよい。この場合のメインローブクラッタ周波数は、(11)式と同様に算出することができ、周波数フィルタ3F1,3F2を適用できる。HPRFの場合にも、パルスブラインドを避けるために、必要に応じて複数のPRFにより送受信を行う。複数のPRFは、各PRFの送受信信号の送信パルスが、極力重なり合わないものを選択する。また、検出後、測距/測速が必要な場合は、FMCW(Frequency Modulated continuous Wave)によるレンジングを行う(非特許文献7参照)。   That is, when the relative speed between the mounted machine and the target is high, transmission and reception by HPRF may be applied in order to widen the clutter-free area on the Doppler frequency axis. The main lobe clutter frequency in this case can be calculated in the same manner as the equation (11), and the frequency filters 3F1 and 3F2 can be applied. Also in the case of HPRF, in order to avoid pulse blinds, transmission / reception is performed using a plurality of PRFs as necessary. The plurality of PRFs are selected such that the transmission pulses of the transmission / reception signals of each PRF do not overlap as much as possible. If ranging / speed measurement is required after detection, ranging is performed by FMCW (Frequency Modulated Continuous Wave) (see Non-Patent Document 7).

まず、距離、速度とビート周波数の関係は、次式で与えられる。

Figure 2015230284
First, the relationship between distance, speed and beat frequency is given by the following equation.
Figure 2015230284

観測されるビート周波数を代入して連立方程式をたてて距離、速度を算出すると、次式となる。

Figure 2015230284
Substituting the observed beat frequency to create simultaneous equations and calculate distance and speed, the following equation is obtained.
Figure 2015230284

複数目標の場合には、複数のビ−ト周波数が得られるので、複数スイープ間のビート周波数において、測角値や振幅等が類似のビート周波数同士を選定(ペアリング)して、連立方程式を解く。または、目標数をNとして、N+1個のスイープ信号の目標のビート周波数を用いて、各々のスイープ間で目標信号のペアリングをして、連立方程式を解くことにより、目標の距離と速度を算出できる。   In the case of multiple targets, multiple beat frequencies can be obtained. Therefore, in the beat frequency between multiple sweeps, select the beat frequencies with similar angle values, amplitudes, etc. solve. Alternatively, N is the target number, and the target beat frequency of N + 1 sweep signals is used, the target signal is paired between each sweep, and the simultaneous equations are solved to calculate the target distance and speed. it can.

すなわち、図19(a)に示す変調方式のスイープ1及びスイープ2に対して、図19(b)に示すように、複数の検出バンクが得られ、スイープ1とスイープ2で別々の周波数のバンクから目標を抽出することが可能となる。   That is, as shown in FIG. 19B, a plurality of detection banks are obtained with respect to the sweep 1 and sweep 2 of the modulation scheme shown in FIG. 19A, and banks having different frequencies for the sweep 1 and the sweep 2 are obtained. It is possible to extract a target from

以上のように、第4の実施形態の構成によれば、HPRFを適用したことにより、クラッタフリー領域を広げることができ、これによって目標を検出し、測距/測速することができる。   As described above, according to the configuration of the fourth embodiment, by applying HPRF, it is possible to widen the clutter-free region, thereby detecting a target and performing distance measurement / speed measurement.

尚、第3の実施形態及び第4の実施形態では、それぞれMPRFとHPRFの場合において、周波数フィルタを用いてクラッタを抑圧する手法について述べた。本手法は、LPRFの場合においても、クラッタの広がりが小さい場合には、周波数フィルタによりクラッタを抑圧できるのは言うまでもない。   In the third and fourth embodiments, the method of suppressing clutter using a frequency filter in the case of MPRF and HPRF has been described. Needless to say, this technique can suppress clutter by a frequency filter when the spread of clutter is small even in the case of LPRF.

尚、本実施形態は上記実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。また、上記実施形態に開示されている複数の構成要素の適宜な組み合わせにより、種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。更に、異なる実施形態にわたる構成要素を適宜組み合わせてもよい。   Note that the present embodiment is not limited to the above-described embodiment as it is, and can be embodied by modifying constituent elements without departing from the scope of the invention in the implementation stage. In addition, various inventions can be formed by appropriately combining a plurality of components disclosed in the embodiment. For example, some components may be deleted from all the components shown in the embodiment. Furthermore, constituent elements over different embodiments may be appropriately combined.

1…アンテナ、2…送受信器、21…送受信部、22…ビーム制御部、3…信号処理器、311,312…AD(Analog-Digital)変換部、321,322…レンジ軸FFT部(FFTx)、331,332…レンジ参照信号乗算部、341,342…データ保存部、351,352…クロスレンジ参照信号乗算部、361,362…クロスレンジ軸FFT部(FFTy)、371,372…レンジ軸逆FFT部(逆FFTx)、38…CFAR処理部、39…測距部、3A…測角部、3B…セル抽出部、3C1,3C2…処理選択部、3D…補正信号抽出部、3E…クロスレンジ参照信号補正部、3F1,3F2…周波数フィルタ、3G…MPRF測距・測速部、3H…HPRF測距・測速部   DESCRIPTION OF SYMBOLS 1 ... Antenna, 2 ... Transmitter / receiver, 21 ... Transmitter / receiver, 22 ... Beam control unit, 3 ... Signal processor, 311, 312 ... AD (Analog-Digital) converter, 321, 322 ... Range axis FFT unit (FFTx) 331, 332... Range reference signal multiplier 341 342 Data storage unit 351 352 Cross range reference signal multiplier 361 362 Cross range axis FFT unit (FFTy) 371 372 Reverse range axis FFT unit (inverse FFTx), 38 ... CFAR processing unit, 39 ... ranging unit, 3A ... anging unit, 3B ... cell extracting unit, 3C1, 3C2 ... processing selecting unit, 3D ... correction signal extracting unit, 3E ... cross range Reference signal correction unit, 3F1, 3F2 ... frequency filter, 3G ... MPRF ranging / speed measuring unit, 3H ... HPRF ranging / speed measuring unit

Claims (8)

レーダパルス受信信号を入力してΣビームのレンジ軸に対する時間軸のコヒーレント積分処理を施すことで目標位置の概略範囲を特定し、前記目標位置の概略範囲に焦点を合わせ、前記時間軸のコヒーレント積分処理結果についてFFT(Fast Fourier Transform)処理により周波数領域で積分処理して、逆FFT処理により時間領域に戻すことでΣビームのレンジ−ドップラ画像を取得するΣビーム画像取得手段と、
前記レーダパルス受信信号を入力してΔビームのレンジ軸に対する時間軸のコヒーレント積分処理を施すことで目標位置の概略範囲を特定し、前記目標位置の概略範囲に焦点を合わせ、前記時間軸のコヒーレント積分処理結果について周波数領域で積分処理して時間領域に戻すことでΔビームのレンジ−ドップラ画像を取得するΔビーム画像取得手段と、
前記Σビーム画像取得手段で取得されたΣビームのレンジ−ドップラ画像に対してCFAR(Constant False Alarm Rate)処理により目標セルを検出して測距し、前記Δビーム画像取得手段で取得されたΔビームのレンジ−ドップラ画像から前記Σビームで検出された目標セルに対応するセルの信号を用いて測角する測距・測角手段と
を具備するレーダ装置。
Specify the approximate range of the target position by inputting the radar pulse received signal and applying time-axis coherent integration processing to the range axis of the Σ beam, focus on the approximate range of the target position, and coherent integration of the time axis Σ beam image acquisition means for acquiring the range-Doppler image of the Σ beam by integrating in the frequency domain by FFT (Fast Fourier Transform) processing results and returning to the time domain by inverse FFT processing;
The radar pulse received signal is input and a time axis coherent integration process is performed on the range axis of the Δ beam to identify an approximate range of the target position, focusing on the approximate range of the target position, and coherent on the time axis Δ beam image acquisition means for acquiring the range-Doppler image of Δ beam by integrating in the frequency domain and returning to the time domain for the integration processing result;
The target cell is detected and measured by CFAR (Constant False Alarm Rate) processing for the range-Doppler image of the Σ beam acquired by the Σ beam image acquisition unit, and Δ acquired by the Δ beam image acquisition unit A radar apparatus comprising a range-finding means for measuring an angle using a cell signal corresponding to a target cell detected by the Σ beam from a beam range-Doppler image.
前記Σビーム画像取得手段は、前記Σビームのレンジ−ドップラ画像の結果より、所定の振幅スレショルドを超えた点を通る所定の範囲の周波数軸の信号出力を抽出し、ドップラ周波数を0にシフトして、それ以外は0埋めした信号を逆FFTして時間領域に戻し、その時間領域の信号の共役複素値を補正係数として、前記コヒーレント積分処理する入力信号をパルス繰り返し期間毎に補正し、その補正結果から再度レンジ−ドップラ画像を算出し、これをM回(M≧1)繰り返して最終のレンジ−ドップラ画像を取得する請求項1記載のレーダ装置。   The Σ beam image acquisition means extracts a signal output of a predetermined range of a frequency axis that passes a point exceeding a predetermined amplitude threshold from the result of the range-Doppler image of the Σ beam, and shifts the Doppler frequency to 0. In other cases, the zero-padded signal is inverse FFTed and returned to the time domain, and the input signal to be subjected to the coherent integration process is corrected for each pulse repetition period using the conjugate complex value of the signal in the time domain as a correction coefficient. The radar apparatus according to claim 1, wherein a range-Doppler image is calculated again from the correction result, and this is repeated M times (M ≧ 1) to obtain a final range-Doppler image. 前記Σビーム画像取得手段は、距離−周波数軸でアンビギュイティを持つMPRF(Middle Pulse Repetition Frequency )による送受信信号について、被搭載機の高度、クラッタ距離、ビーム走査角に応じたドップラ周波数を中心として所定の幅をもつドップラ周波数範囲を抑圧したレンジ−ドップラ画像を取得し、
前記測距・測角手段は、前記測距及び測角の処理に、前記CFAR処理によって得られた目標セルについて行う複数のPRFを用いた請求項1または2記載のレーダ装置。
The Σ-beam image acquisition means is centered on the Doppler frequency according to the altitude of the mounted machine, the clutter distance, and the beam scanning angle for the transmission / reception signal by MPRF (Middle Pulse Repetition Frequency) having ambiguity on the distance-frequency axis. Obtain a range-Doppler image in which the Doppler frequency range having a predetermined width is suppressed,
The radar apparatus according to claim 1, wherein the distance measurement / angle measurement means uses a plurality of PRFs performed on a target cell obtained by the CFAR process for the distance measurement and angle measurement processing.
前記Σビーム画像取得手段は、速度軸のアンビギュイティは発生させず、距離軸でアンビギュイティを持つHPRF(High Pulse Repetition Frequency )による送受信信号について、被搭載機の高度、クラッタ距離、ビーム走査角に応じたドップラ周波数を中心として所定の幅をもつドップラ周波数範囲を抑圧したレンジ−ドップラ画像を取得し、
前記測距・測角手段は、前記測距及び測角の処理に、前記CFAR処理によって得られた目標セルについて行う複数のスイープによるFMCW(Frequency Modulated continuous Wave)を用いた請求項1または2記載のレーダ装置。
The Σ beam image acquisition means does not generate velocity axis ambiguity, and transmits and receives signals by HPRF (High Pulse Repetition Frequency) having ambiguity on the distance axis. A range-Doppler image in which a Doppler frequency range having a predetermined width around a Doppler frequency according to the angle is suppressed is acquired,
The distance measurement / angle measurement means uses FMCW (Frequency Modulated Continuous Wave) by a plurality of sweeps performed on the target cell obtained by the CFAR process for the distance measurement and angle measurement processing. Radar equipment.
レーダパルス受信信号を入力してΣビームのレンジ軸に対する時間軸のコヒーレント積分処理を施すことで目標位置の概略範囲を特定し、前記目標位置の概略範囲に焦点を合わせ、前記時間軸のコヒーレント積分処理結果についてFFT(Fast Fourier Transform)処理により周波数領域で積分処理して、逆FFT処理により時間領域に戻すことでΣビームのレンジ−ドップラ画像を取得し、
前記レーダパルス受信信号を入力してΔビームのレンジ軸に対する時間軸のコヒーレント積分処理を施すことで目標位置の概略範囲を特定し、前記目標位置の概略範囲に焦点を合わせ、前記時間軸のコヒーレント積分処理結果について周波数領域で積分処理して時間領域に戻すことでΔビームのレンジ−ドップラ画像を取得し、
前記Σビームのレンジ−ドップラ画像に対してCFAR(Constant False Alarm Rate)処理により目標セルを検出して測距し、前記Δビームのレンジ−ドップラ画像から前記Σビームで検出された目標セルに対応するセルの信号を用いて測角するレーダ装置のレーダ信号処理方法。
Specify the approximate range of the target position by inputting the radar pulse received signal and applying time-axis coherent integration processing to the range axis of the Σ beam, focus on the approximate range of the target position, and coherent integration of the time axis About the processing result, an integration process in the frequency domain is performed by FFT (Fast Fourier Transform) processing, and a range-Doppler image of Σ beam is acquired by returning to the time domain by inverse FFT processing,
The radar pulse received signal is input and a time axis coherent integration process is performed on the range axis of the Δ beam to identify an approximate range of the target position, focusing on the approximate range of the target position, and coherent on the time axis The integration process result is integrated in the frequency domain and returned to the time domain to obtain a Δ beam range-Doppler image,
The target cell is detected and measured by CFAR (Constant False Alarm Rate) processing on the range-Doppler image of the Σ beam and corresponds to the target cell detected by the Σ beam from the range-Doppler image of the Δ beam. A radar signal processing method of a radar apparatus that measures an angle using a signal of a cell to be operated.
前記Σビームのレンジ−ドップラ画像の結果より、所定の振幅スレショルドを超えた点を通る所定の範囲の周波数軸の信号出力を抽出し、ドップラ周波数を0にシフトして、それ以外は0埋めした信号を逆FFTして時間領域に戻し、その時間領域の信号の共役複素値を補正係数として、前記コヒーレント積分処理する入力信号をパルス繰り返し期間毎に補正し、その補正結果から再度レンジ−ドップラ画像を算出し、これをM回(M≧1)繰り返して最終のレンジ−ドップラ画像を取得する請求項5記載のレーダ装置のレーダ信号処理方法。   From the result of the range-Doppler image of the Σ beam, the signal output of the frequency axis in a predetermined range passing through the point exceeding the predetermined amplitude threshold is extracted, the Doppler frequency is shifted to 0, and the others are filled with 0. The signal is inverse FFTed and returned to the time domain, and the input complex signal is corrected for each pulse repetition period using the conjugate complex value of the signal in the time domain as a correction coefficient, and the range-Doppler image is again obtained from the correction result. 6. The radar signal processing method for a radar apparatus according to claim 5, wherein the final range-Doppler image is acquired by repeating the calculation M times (M ≧ 1). 前記Σビームのレンジ−ドップラ画像は、距離−周波数軸でアンビギュイティを持つMPRF(Middle Pulse Repetition Frequency )による送受信信号について、被搭載機の高度、クラッタ距離、ビーム走査角に応じたドップラ周波数を中心として所定の幅をもつドップラ周波数範囲を抑圧して取得し、
前記測距及び測角の処理に、前記CFAR処理によって得られた目標セルについて行う複数のPRFを用いた請求項5または6記載のレーダ装置のレーダ信号処理方法。
The range-Doppler image of the Σ beam shows the Doppler frequency according to the altitude, clutter distance, and beam scanning angle of the mounted machine for the transmission / reception signal by MPRF (Middle Pulse Repetition Frequency) having ambiguity in the distance-frequency axis. Obtained by suppressing the Doppler frequency range with a predetermined width as the center,
The radar signal processing method of a radar apparatus according to claim 5 or 6, wherein a plurality of PRFs used for a target cell obtained by the CFAR processing are used for the ranging and angle measurement processing.
前記Σビームのレンジ−ドップラ画像は、速度軸のアンビギュイティは発生させず、距離軸でアンビギュイティを持つHPRF(High Pulse Repetition Frequency )による送受信信号について、被搭載機の高度、クラッタ距離、ビーム走査角に応じたドップラ周波数を中心として所定の幅をもつドップラ周波数範囲を抑圧して取得し、
前記測距及び測角の処理に、前記CFAR処理によって得られた目標セルについて行う複数のスイープによるFMCW(Frequency Modulated continuous Wave)を用いた請求項5または6記載のレーダ装置のレーダ信号処理方法。
The range-Doppler image of the Σ beam does not generate ambiguity of the velocity axis, and transmits and receives signals by HPRF (High Pulse Repetition Frequency) having ambiguity on the distance axis. Obtained by suppressing a Doppler frequency range having a predetermined width around the Doppler frequency according to the beam scanning angle,
7. The radar signal processing method of a radar apparatus according to claim 5, wherein FMCW (Frequency Modulated Continuous Wave) based on a plurality of sweeps performed on the target cell obtained by the CFAR process is used for the ranging and angle measurement processes.
JP2014117835A 2014-06-06 2014-06-06 Radar apparatus and radar signal processing method thereof Active JP6352688B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2014117835A JP6352688B2 (en) 2014-06-06 2014-06-06 Radar apparatus and radar signal processing method thereof

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2014117835A JP6352688B2 (en) 2014-06-06 2014-06-06 Radar apparatus and radar signal processing method thereof

Publications (2)

Publication Number Publication Date
JP2015230284A true JP2015230284A (en) 2015-12-21
JP6352688B2 JP6352688B2 (en) 2018-07-04

Family

ID=54887112

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2014117835A Active JP6352688B2 (en) 2014-06-06 2014-06-06 Radar apparatus and radar signal processing method thereof

Country Status (1)

Country Link
JP (1) JP6352688B2 (en)

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106597411A (en) * 2016-12-30 2017-04-26 无锡市雷华科技有限公司 Radar signal processing method
CN106646447A (en) * 2017-01-18 2017-05-10 武汉雷博合创电子技术有限公司 Detection method for radar target long-time accumulation based on linear frequency modulation continuous wave
CN107132534A (en) * 2017-06-27 2017-09-05 西安电子科技大学 A kind of optimization method of High-Speed RADAR target frequency domain detection
JP2017161358A (en) * 2016-03-09 2017-09-14 株式会社東芝 Radar device
JP2018205174A (en) * 2017-06-06 2018-12-27 株式会社東芝 Radar device and radar signal processing method thereof
KR20190014917A (en) * 2017-08-04 2019-02-13 국방과학연구소 Tracking prf selection method of mprf mode and aircraft radar using the same
JP2020051752A (en) * 2018-09-21 2020-04-02 株式会社東芝 Radar device and radar device pri control method
CN112462341A (en) * 2020-10-20 2021-03-09 西南石油大学 Small rotor unmanned aerial vehicle target detection method based on multi-pulse accumulation
CN113330324A (en) * 2019-01-31 2021-08-31 三菱电机株式会社 Radar apparatus and signal processing method
CN113341404A (en) * 2020-03-02 2021-09-03 加特兰微电子科技(上海)有限公司 Constant false alarm detection method and device, integrated circuit and radio device
JPWO2022029865A1 (en) * 2020-08-04 2022-02-10
JP7055915B1 (en) 2021-03-23 2022-04-18 株式会社東芝 Target extraction device and target extraction method
CN117572426A (en) * 2024-01-17 2024-02-20 长沙莫之比智能科技有限公司 Vehicle-mounted millimeter wave radar static target detection method

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001194456A (en) * 2000-01-14 2001-07-19 Mitsubishi Electric Corp Composite aperture radar
JP2006343290A (en) * 2005-06-10 2006-12-21 Mitsubishi Electric Corp Imaging radar device
US7391357B1 (en) * 2005-03-25 2008-06-24 Sandia Corporation Correction of motion measurement errors beyond the range resolution of a synthetic aperture radar
JP2008170323A (en) * 2007-01-12 2008-07-24 Mitsubishi Electric Corp Radar device
JP2010145346A (en) * 2008-12-22 2010-07-01 Mitsubishi Electric Corp Pulse compression radar system
JP2013019702A (en) * 2011-07-07 2013-01-31 Mitsubishi Heavy Ind Ltd Radar device and target position calculating method
JP2013148377A (en) * 2012-01-17 2013-08-01 Mitsubishi Electric Corp Signal processing device
JP2013181864A (en) * 2012-03-02 2013-09-12 Mitsubishi Electric Corp Synthetic aperture radar device and signal processor

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001194456A (en) * 2000-01-14 2001-07-19 Mitsubishi Electric Corp Composite aperture radar
US7391357B1 (en) * 2005-03-25 2008-06-24 Sandia Corporation Correction of motion measurement errors beyond the range resolution of a synthetic aperture radar
JP2006343290A (en) * 2005-06-10 2006-12-21 Mitsubishi Electric Corp Imaging radar device
JP2008170323A (en) * 2007-01-12 2008-07-24 Mitsubishi Electric Corp Radar device
JP2010145346A (en) * 2008-12-22 2010-07-01 Mitsubishi Electric Corp Pulse compression radar system
JP2013019702A (en) * 2011-07-07 2013-01-31 Mitsubishi Heavy Ind Ltd Radar device and target position calculating method
JP2013148377A (en) * 2012-01-17 2013-08-01 Mitsubishi Electric Corp Signal processing device
JP2013181864A (en) * 2012-03-02 2013-09-12 Mitsubishi Electric Corp Synthetic aperture radar device and signal processor

Cited By (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2017161358A (en) * 2016-03-09 2017-09-14 株式会社東芝 Radar device
CN106597411B (en) * 2016-12-30 2019-06-04 无锡市雷华科技有限公司 Method for processing radar signals
CN106597411A (en) * 2016-12-30 2017-04-26 无锡市雷华科技有限公司 Radar signal processing method
CN106646447A (en) * 2017-01-18 2017-05-10 武汉雷博合创电子技术有限公司 Detection method for radar target long-time accumulation based on linear frequency modulation continuous wave
CN106646447B (en) * 2017-01-18 2019-03-26 武汉雷博合创电子科技有限公司 Radar target long time integration detection method based on linear frequency modulation continuous wave
JP2018205174A (en) * 2017-06-06 2018-12-27 株式会社東芝 Radar device and radar signal processing method thereof
CN107132534A (en) * 2017-06-27 2017-09-05 西安电子科技大学 A kind of optimization method of High-Speed RADAR target frequency domain detection
KR101993855B1 (en) 2017-08-04 2019-06-27 국방과학연구소 Tracking prf selection method of mprf mode and aircraft radar using the same
KR20190014917A (en) * 2017-08-04 2019-02-13 국방과학연구소 Tracking prf selection method of mprf mode and aircraft radar using the same
JP2020051752A (en) * 2018-09-21 2020-04-02 株式会社東芝 Radar device and radar device pri control method
JP7242227B2 (en) 2018-09-21 2023-03-20 株式会社東芝 RADAR SYSTEM AND PRI CONTROL METHOD OF RADAR SYSTEM
CN113330324A (en) * 2019-01-31 2021-08-31 三菱电机株式会社 Radar apparatus and signal processing method
CN113330324B (en) * 2019-01-31 2023-06-06 三菱电机株式会社 Radar apparatus and signal processing method
CN113341404A (en) * 2020-03-02 2021-09-03 加特兰微电子科技(上海)有限公司 Constant false alarm detection method and device, integrated circuit and radio device
JP7170948B2 (en) 2020-08-04 2022-11-14 三菱電機株式会社 SIGNAL PROCESSING DEVICE, RADAR AND RADAR SIGNAL PROCESSING METHOD
JPWO2022029865A1 (en) * 2020-08-04 2022-02-10
CN112462341A (en) * 2020-10-20 2021-03-09 西南石油大学 Small rotor unmanned aerial vehicle target detection method based on multi-pulse accumulation
CN112462341B (en) * 2020-10-20 2022-06-17 西南石油大学 Small rotor unmanned aerial vehicle target detection method based on multi-pulse accumulation
JP2022147383A (en) * 2021-03-23 2022-10-06 株式会社東芝 Target extraction device and target extraction method
JP7055915B1 (en) 2021-03-23 2022-04-18 株式会社東芝 Target extraction device and target extraction method
CN117572426A (en) * 2024-01-17 2024-02-20 长沙莫之比智能科技有限公司 Vehicle-mounted millimeter wave radar static target detection method

Also Published As

Publication number Publication date
JP6352688B2 (en) 2018-07-04

Similar Documents

Publication Publication Date Title
JP6352688B2 (en) Radar apparatus and radar signal processing method thereof
CN108051809B (en) Moving target imaging method and device based on Radon transformation and electronic equipment
US9335409B2 (en) Bistatic inverse synthetic aperture radar imaging
Munoz-Ferreras et al. Uniform rotational motion compensation for inverse synthetic aperture radar with non-cooperative targets
Huang et al. GMTI and parameter estimation via time-Doppler chirp-varying approach for single-channel airborne SAR system
JP6301749B2 (en) Doppler radar apparatus and radar signal processing method thereof
Zhu et al. A new slant-range velocity ambiguity resolving approach of fast moving targets for SAR system
JP6470152B2 (en) Radar apparatus and radar signal processing method
JP2014182010A (en) Radar apparatus
JP6367143B2 (en) Synthetic aperture radar apparatus and radar signal processing method thereof
JP6466263B2 (en) Radar apparatus and radar signal processing method
JP6546109B2 (en) Radar equipment
JP6178244B2 (en) Synthetic aperture radar apparatus and image processing method thereof
JP6367140B2 (en) Radar apparatus and radar signal processing method
Uysal et al. The effect of moving target on range-Doppler map and backprojection algorithm for focusing
Vu et al. Space time adaptive processing for moving target detection and imaging in bistatic SAR
JP6612728B2 (en) Radar apparatus and radar signal processing method thereof
Sanz-Marcos et al. Bistatic fixed-receiver parasitic SAR processor based on the back-propagation algorithm
Vu et al. Two-dimensional spectrum for BiSAR derivation based on Lagrange inversion theorem
Chen et al. Parametric Iterative Soft Thresholding Algorithm for Refocusing of Moving Targets in SAR Images
JP6400340B2 (en) Pulse compression radar apparatus and radar signal processing method thereof
JP6363524B2 (en) Radar apparatus and radar signal processing method
JP2005083814A (en) Radar system
Trischman et al. Inverse synthetic aperture radar imaging at 580 GHz
Gromek et al. A concept of the parametric autofocus method for passive ISAR imaging

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20170220

A711 Notification of change in applicant

Free format text: JAPANESE INTERMEDIATE CODE: A712

Effective date: 20170911

Free format text: JAPANESE INTERMEDIATE CODE: A711

Effective date: 20170911

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20180112

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20180123

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20180322

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: 20180508

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20180607

R150 Certificate of patent or registration of utility model

Ref document number: 6352688

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150