JP3223697B2 - Radar signal processing equipment - Google Patents

Radar signal processing equipment

Info

Publication number
JP3223697B2
JP3223697B2 JP06979694A JP6979694A JP3223697B2 JP 3223697 B2 JP3223697 B2 JP 3223697B2 JP 06979694 A JP06979694 A JP 06979694A JP 6979694 A JP6979694 A JP 6979694A JP 3223697 B2 JP3223697 B2 JP 3223697B2
Authority
JP
Japan
Prior art keywords
frequency
circuit
amplitude
waveform
area
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.)
Expired - Fee Related
Application number
JP06979694A
Other languages
Japanese (ja)
Other versions
JPH07280923A (en
Inventor
隆三郎 臼井
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.)
Mitsubishi Electric Corp
Original Assignee
Mitsubishi Electric 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 Mitsubishi Electric Corp filed Critical Mitsubishi Electric Corp
Priority to JP06979694A priority Critical patent/JP3223697B2/en
Publication of JPH07280923A publication Critical patent/JPH07280923A/en
Application granted granted Critical
Publication of JP3223697B2 publication Critical patent/JP3223697B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Description

【発明の詳細な説明】DETAILED DESCRIPTION OF THE INVENTION

【0001】[0001]

【産業上の利用分野】この発明は、例えば、高分解能レ
ーダ装置のレーダ信号処理装置における目標受信信号の
位相ずれの補正に関するものである。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to, for example, correction of a phase shift of a target reception signal in a radar signal processing device of a high resolution radar device.

【0002】[0002]

【従来の技術】図6は高分解能レーダ装置のレーダ信号
処理装置の構成図であり、図中、1はレーダ装置から入
力された目標受信信号を内部で処理できるデータ形式に
変換するデータインタフェース部、2はデータインタフ
ェース部1で変換された目標受信信号をパルス圧縮する
パルス圧縮部、3はパルス圧縮部2でパルス圧縮された
目標受信信号の時間による距離ずれを補正する距離補正
部、4は距離補正部3で補正された目標受信信号の時間
による位相ずれを補正する位相補正部、5は位相補正部
4で補正された目標受信信号のドップラ周波数を分離す
る周波数分析部、6は周波数分析部5で周波数分析され
た目標受信信号の周波数スペクトルを画像データに変換
する検波部、7は検波部6で得られた画像データと表示
器のインタフェースを調整し表示画像データを生成する
表示器インタフェース部、SMはレーダ装置から入力さ
れた目標受信信号、RSは距離補正部3で時間による距
離ずれが補正された目標受信信号、RGはレーダと目標
重心との初期距離、RDは位相補正部4で時間による位
相ずれが補正された目標受信信号、Dは表示画像データ
である。
2. Description of the Related Art FIG. 6 is a block diagram of a radar signal processing device of a high-resolution radar device. In the drawing, reference numeral 1 denotes a data interface unit for converting a target reception signal input from the radar device into a data format that can be processed internally. Reference numeral 2 denotes a pulse compression unit for pulse-compressing the target reception signal converted by the data interface unit 1. Reference numeral 3 denotes a distance correction unit for correcting a time-dependent distance shift of the target reception signal pulse-compressed by the pulse compression unit 2. A phase correction unit for correcting a phase shift due to time of the target reception signal corrected by the distance correction unit 3; a frequency analysis unit 5 for separating the Doppler frequency of the target reception signal corrected by the phase correction unit 4; A detector 7 converts the frequency spectrum of the target reception signal, which has been frequency-analyzed by the unit 5, into image data. Reference numeral 7 denotes an interface between the image data obtained by the detector 6 and the display. A display interface unit for adjusting display time and generating display image data, SM is a target reception signal input from the radar device, RS is a target reception signal whose distance deviation due to time has been corrected by the distance correction unit 3, and RG is a radar and target. An initial distance from the center of gravity, RD is a target reception signal whose phase shift due to time has been corrected by the phase correction unit 4, and D is display image data.

【0003】図7は図6のレーダ信号処理装置における
従来の位相補正部4の構成図であり、図中、RS,R
G,RD及び4は図6と同じであり、8は距離補正部3
より出力された時間による距離ずれが補正された目標受
信信号RS及びレーダと目標重心との初期距離RGを格
納するバッファ回路、9はバッファ回路8より出力され
たレーダと目標重心との初期距離RGにおける目標受信
信号を時間方向に小区間で周波数分析する区分周波数分
析回路、10は区分周波数分析回路9で得られた周波数
と振幅の波形に対して振幅値が最大となる周波数を基準
点周波数として検出する振幅値最大検出回路、11は振
幅値最大検出回路10で検出された基準点周波数の時間
方向に対する軌跡を平滑化する平滑化回路、12は平滑
化回路11で平滑化された基準点周波数の時間方向に対
する軌跡から位相補正量を算出する位相補正量算出回
路、13は位相補正量算出回路12で算出された位相補
正量を用いてバッファ回路8より出力された目標受信信
号の位相を補正する位相補正回路、GSはレーダと目標
重心との初期距離RGにおける目標受信信号である。
FIG. 7 is a block diagram of a conventional phase correction unit 4 in the radar signal processing apparatus of FIG.
G, RD and 4 are the same as in FIG.
A buffer circuit for storing the target reception signal RS and the initial distance RG between the radar and the target center of gravity corrected by the output time-dependent distance shift, and an initial distance RG between the radar output from the buffer circuit 8 and the target center of gravity. A frequency analysis circuit for frequency-analyzing the target reception signal in a small section in the time direction at 10 in the frequency and amplitude waveforms obtained by the frequency division analysis circuit 9 as a reference point frequency. A maximum amplitude value detection circuit to be detected, 11 is a smoothing circuit for smoothing the locus of the reference point frequency detected by the maximum amplitude value detection circuit 10 in the time direction, and 12 is a reference point frequency smoothed by the smoothing circuit 11. A phase correction amount calculating circuit for calculating a phase correction amount from a trajectory in the time direction of the phase correction amount calculating circuit; Phase correcting circuit for correcting the phase of the target received signal outputted from the circuit 8, GS is the target received signal in the initial distance RG between the radar and the target centroid.

【0004】次に、動作について説明する。レーダ装置
から入力された目標受信信号SMはデータインタフェー
ス部1で内部で処理できるデータ形式に変換され、パル
ス圧縮部2でパルス圧縮された後、距離補正部3で時間
による距離ずれを補正し、目標受信信号RSとして位相
補正部4に出力する。また、距離補正部3ではレーダと
目標重心との初期距離RGを算出し、位相補正部4に出
力する。位相補正部4では時間による距離ずれが補正さ
れた目標受信信号RS及びレーダと目標重心との初期距
離RGを用いて目標受信信号RSの時間による位相ずれ
を補正し、目標受信信号RDとして周波数分析部5に出
力する。この目標受信信号RDは周波数分析部5で周波
数分析されることによって周波数スペクトルに変換さ
れ、検波部6で画像データに変換された後、表示器イン
タフェース部7で表示器とのインタフェースを調整し、
表示画像データDとして出力される。
Next, the operation will be described. The target reception signal SM input from the radar device is converted into a data format that can be processed internally by the data interface unit 1, pulse-compressed by the pulse compression unit 2, and corrected by the distance correction unit 3 for a distance shift due to time. The signal is output to the phase correction unit 4 as the target reception signal RS. Further, the distance correction unit 3 calculates an initial distance RG between the radar and the target center of gravity, and outputs the calculated initial distance RG to the phase correction unit 4. The phase correction unit 4 corrects the time-dependent phase shift of the target reception signal RS using the target reception signal RS whose distance shift due to time has been corrected and the initial distance RG between the radar and the target center of gravity, and performs frequency analysis as the target reception signal RD. Output to section 5. The target reception signal RD is converted into a frequency spectrum by being frequency-analyzed by the frequency analysis unit 5 and converted into image data by the detection unit 6, and then the display interface unit 7 adjusts the interface with the display device.
It is output as display image data D.

【0005】次に、位相補正部4の動作について説明す
る。距離補正部3から入力された時間による距離ずれが
補正された目標受信信号RS及びレーダと目標重心との
初期距離RGはバッファ回路8に格納され、目標受信信
号RS及びレーダと目標重心との初期距離RGにおける
目標受信信号GSとして出力される。このレーダと目標
重心との初期距離RGにおける目標受信信号GSは区分
周波数分析回路9で時間方向に小区間で周波数分析さ
れ、得られた周波数と振幅の波形に対して振幅値最大検
出回路10で振幅値が最大となる周波数を基準点周波数
として検出した後、平滑化回路11に出力される。平滑
化回路11では基準点周波数の時間方向に対する軌跡を
平滑化し、平滑化された軌跡から位相補正量算出回路1
2で位相補正量を算出する。位相補正回路13は位相補
正量算出回路12で算出した位相補正量を用いてバッフ
ァ回路8から出力された時間による距離ずれが補正され
た目標受信信号RSの位相を補正し、時間による位相ず
れが補正された目標受信信号RDとして周波数分析部5
へ出力する。
Next, the operation of the phase correction section 4 will be described. The target reception signal RS and the initial distance RG between the radar and the target center of gravity, the distance of which has been corrected by the time input from the distance correction unit 3, are stored in the buffer circuit 8, and the initial reception signal RS and the initial distance between the radar and the target center of gravity are stored. It is output as the target reception signal GS at the distance RG. The target reception signal GS at the initial distance RG between the radar and the target center of gravity is frequency-analyzed in a small section in the time direction by the divided frequency analysis circuit 9, and the obtained maximum frequency value and amplitude waveform are detected by the maximum amplitude value detection circuit 10. After detecting the frequency at which the amplitude value becomes maximum as the reference point frequency, it is output to the smoothing circuit 11. The smoothing circuit 11 smoothes the locus of the reference point frequency in the time direction, and calculates the phase correction amount calculating circuit 1 from the smoothed locus.
In step 2, the amount of phase correction is calculated. The phase correction circuit 13 uses the phase correction amount calculated by the phase correction amount calculation circuit 12 to correct the phase of the target reception signal RS output from the buffer circuit 8 and corrected for the distance shift due to time. The frequency analysis unit 5 determines the corrected target reception signal RD as
Output to

【0006】更に、位相補正部4を図8を用いて説明す
る。図8は位相補正部4の処理方法を示した図である。
時間による距離ずれが補正された目標受信信号RSをS
i,j(ここでiはレンジビン番号、jはパルスヒット番
号、i,jは自然数である。)、レーダと目標重心との
初期距離RGの存在するレンジビン番号をrと定義する
とレーダと目標重心との初期距離RGにおける目標受信
信号GSはSr,j と表され、図8(a)のような波形が
得られる。Sr,j に対し、区分周波数分析回路9で時間
方向(パルスヒット方向)に小区間で周波数分析を行う
と図8(b)のような波形が得られ、周波数fm と振幅
k m の関係(ここでkは区分周波数分析番号、mは周
波数ビン番号、k,mは自然数である。)は“数1”で
表される。
Further, the phase correction section 4 will be described with reference to FIG. FIG. 8 is a diagram illustrating a processing method of the phase correction unit 4.
The target reception signal RS whose distance deviation due to time has been corrected is represented by S
i, j (where i is the range bin number, j is the pulse hit number, and i, j are natural numbers), and the range bin number where the initial distance RG between the radar and the target center of gravity is defined as r, the radar and the target center of gravity are defined. The target reception signal GS at an initial distance RG from the target reception signal RG is represented as S r, j, and a waveform as shown in FIG. 8A is obtained. S r, j hand, when the frequency analysis in small sections to piecewise frequency analysis circuit 9 time direction (pulse hit direction) waveform as shown in FIG. 8 (b) is obtained, the frequency f m and an amplitude A k m (Where k is a division frequency analysis number, m is a frequency bin number, and k and m are natural numbers) are represented by “Equation 1”.

【0007】[0007]

【数1】 (Equation 1)

【0008】振幅値最大検出回路10で各区分周波数分
析番号kに対し、振幅Ak m が最大値をとる時の周波数
を検出し、それを基準点周波数fk とすると、時間tk
と基準点周波数fk の関係は図8(c)のプロットのよ
うになる。平滑化回路11で図8(c)のプロットに対
し、平滑化を行うと図8(c)の実線のような波形が得
られ、時間tk と周波数f’k の関係は“数2”で表さ
れる。
[0008] In the amplitude maximum detection circuit 10 for each division frequency analysis number k, the amplitude A k m detects the frequency at which the maximum value to it as a reference point frequency f k, the time t k
The relationship between and the reference point frequency f k is as shown in the plot of FIG. To plot in Figure smoothing circuit 11 8 (c), the smoothing waveform shown by the solid line shown in FIG. 8 (c) is obtained, the relationship between the time t k and the frequency f 'k' number 2 ' It is represented by

【0009】[0009]

【数2】 (Equation 2)

【0010】位相補正量算出回路12では、位相補正量
j を“数3”で算出する。
The phase correction amount calculation circuit 12 calculates the phase correction amount Wj by "Equation 3".

【0011】[0011]

【数3】 (Equation 3)

【0012】位相補正回路13では、Si,j の位相を位
相補正量Wj を用いて“数4”で補正する。但し、時間
による位相ずれが補正された目標受信信号RDをS’
i,j と定義する。
The phase correction circuit 13 corrects the phase of S i, j by “Equation 4” using the phase correction amount W j . However, the target reception signal RD whose phase shift due to time has been corrected is represented by S ′.
Define as i, j .

【0013】[0013]

【数4】 (Equation 4)

【0014】[0014]

【発明が解決しょうとする課題】上記のような従来のレ
ーダ信号処理装置では、区分周波数分析後の周波数と振
幅の波形が多峰で、かつ振幅値が最大となる峰の位置が
時間毎に大きく変動するような場合、検出した基準点周
波数がふらついて正確な位相補量を算出することがで
きず、画像がぼけたり、にじんだりするという課題があ
った。
In the conventional radar signal processing apparatus as described above, the frequency and amplitude waveforms after the divisional frequency analysis have multiple peaks, and the position of the peak at which the amplitude value is maximum is located at every time. If such fluctuates greatly, it is impossible to calculate an accurate phase compensation amount unsteady detected reference point frequency, the image is blurred, there is a problem that blurred.

【0015】この発明はかかる課題を解決するためにな
されたもので、検出した基準点周波数がふらついて正確
な位相補量を算出することができず、画像がぼけた
り、にじんだりするということを防止するレーダ信号処
理装置を得ることを目的とする。
[0015] The present invention has been made to solve the above problem, it is impossible to calculate an accurate phase compensation amount unsteady detected reference point frequency, the image is blurred, that is disappear, It is an object of the present invention to obtain a radar signal processing device that prevents the occurrence of a radar signal.

【0016】[0016]

【課題を解決するための手段】この発明に係るレーダ信
号処理装置では、上記位相補正部を距離補正部で補正さ
れた目標受信信号及びレーダと目標重心との初期距離を
格納するバッファ回路と、上記バッファ回路より出力さ
れたレーダと目標重心との初期距離における目標受信信
号を時間方向に小区間で周波数分析する区分周波数分析
回路と、上記区分周波数分析回路で得られた周波数と振
幅の波形を周波数方向に補間する補間回路と、周波数と
振幅の波形について許容される周波数領域内に基準とな
る周波数を設定する周波数設定回路と、周波数と振幅の
波形を許容される周波数領域の開始点から上記周波数設
定回路で設定した周波数まで積分して面積を算出する第
1の積分回路と、上記第1の積分回路で算出された面積
が半分になる時の周波数を算出する第1の面積中心算出
回路と、周波数と振幅の波形を上記周波数設定回路で設
定した周波数から許容される周波数領域の終了点まで積
分して面積を算出する第2の積分回路と、上記第2の積
分回路で算出された面積が半分になる時の周波数を算出
する第2の面積中心算出回路と、上記第1の面積中心算
出回路で算出された周波数と上記第2の面積中心算出回
路で算出された周波数の中心周波数を基準点周波数とし
て検出する中心周波数検出回路と、上記中心周波数検出
回路で検出された基準点周波数の時間方向に対する軌跡
を平滑化する平滑化回路と、上記平滑化回路で平滑化さ
れた基準点周波数の時間方向に対する軌跡から位相補正
量を算出する位相補正量算出回路と、上記位相補正量算
出回路で算出された位相補正量を用いて上記バッファ回
路から出力された目標受信信号の位相を補正する位相補
正回路とによって構成したものである。
In the radar signal processing apparatus according to the present invention, the phase correction unit stores a target reception signal corrected by the distance correction unit and an initial distance between the radar and the target center of gravity. A divided frequency analysis circuit that analyzes the frequency of a target reception signal at an initial distance between the radar output from the buffer circuit and the target center of gravity in a small section in the time direction, and a waveform of the frequency and amplitude obtained by the divided frequency analysis circuit. An interpolation circuit for interpolating in the frequency direction, a frequency setting circuit for setting a reference frequency within a frequency region allowed for the frequency and amplitude waveforms, and a frequency and amplitude waveform starting from the start point of the allowed frequency region. A first integration circuit for calculating an area by integrating up to a frequency set by the frequency setting circuit, and a first integration circuit for calculating the area when the area calculated by the first integration circuit is halved. A first area center calculating circuit for calculating a wave number, a second integrating circuit for calculating an area by integrating the frequency and amplitude waveforms from the frequency set by the frequency setting circuit to an end point of an allowable frequency region; A second area center calculation circuit for calculating a frequency when the area calculated by the second integration circuit is halved, a frequency calculated by the first area center calculation circuit, and the second area A center frequency detection circuit that detects the center frequency of the frequency calculated by the center calculation circuit as a reference point frequency, and a smoothing circuit that smoothes a locus in the time direction of the reference point frequency detected by the center frequency detection circuit, A phase correction amount calculating circuit for calculating a phase correction amount from a trajectory of the reference point frequency smoothed by the smoothing circuit in the time direction, and a phase correction amount calculated by the phase correction amount calculating circuit. Te is obtained constituted by a phase correction circuit for correcting the phase of the target received signal output from the buffer circuit.

【0017】また、上記位相補正部に上記補間回路で補
間された周波数と振幅の波形についてデータを切り出す
データ切り出し回路と、上記データ切り出し回路でデー
タを切り出す際に必要な閾値を設定する閾値設定回路を
設ける。
Further, a data extracting circuit for extracting data on the frequency and amplitude waveforms interpolated by the interpolating circuit in the phase correction section, and a threshold setting circuit for setting a threshold necessary for extracting data by the data extracting circuit Is provided.

【0018】また、上記位相補正部に上記補間回路で補
間された周波数と振幅の波形について各振幅値のべき乗
を計算するべき乗計算回路を設ける。
Further, the phase correction section is provided with a power calculation circuit for calculating a power of each amplitude value for the frequency and amplitude waveforms interpolated by the interpolation circuit.

【0019】また、上記位相補正部に上記補間回路で補
間された周波数と振幅の波形について各振幅値のべき乗
を計算するべき乗計算回路と、上記べき乗計算回路で計
算された周波数と振幅のべき乗の波形についてデータを
切り出すデータ切り出し回路と、上記データ切り出し回
路でデータを切り出す際に必要な閾値を設定する閾値設
定回路を設ける。
The phase correction unit calculates a power of each amplitude value for the frequency and amplitude waveforms interpolated by the interpolation circuit, and a power of the frequency and amplitude calculated by the power calculation circuit. A data extracting circuit for extracting data on a waveform and a threshold setting circuit for setting a threshold necessary for extracting data by the data extracting circuit are provided.

【0020】また、上記位相補正部に上記補間回路で補
間された周波数と振幅の波形についてデータを切り出す
データ切り出し回路と、上記データ切り出し回路でデー
タを切り出す際に必要な閾値を設定する閾値設定回路
と、上記データ切り出し回路で切り出された周波数と振
幅の波形について各振幅値のべき乗を計算するべき乗計
算回路を設ける。
Further, a data extraction circuit for extracting data on the frequency and amplitude waveforms interpolated by the interpolation circuit in the phase correction section, and a threshold setting circuit for setting a threshold necessary for extracting data by the data extraction circuit And a power calculation circuit for calculating a power of each amplitude value for the frequency and amplitude waveforms cut out by the data cutout circuit.

【0021】[0021]

【作用】この発明においては、区分周波数分析後の周波
数と振幅の波形を周波数方向に補間し、補間された波形
を許容される周波数領域の開始点から周波数設定回路で
設定した周波数まで積分して得られた面積が半分になる
時の周波数と、周波数設定回路で設定した周波数から許
容される周波数領域の終了点まで積分して得られた面積
が半分になる時の周波数との中心周波数を基準点周波数
として検出するので、周波数と振幅の波形が多峰で、か
つ振幅値が最大となる峰の位置が時間毎に大きく変動す
るような場合でも変動の影響を受けにくく、安定して基
準点周波数を検出することができ、画像のぼけやにじみ
を除去することができる。
According to the present invention, the frequency and amplitude waveforms after the divided frequency analysis are interpolated in the frequency direction, and the interpolated waveform is integrated from the start point of the allowable frequency region to the frequency set by the frequency setting circuit. The reference frequency is the center frequency between the frequency at which the obtained area is halved and the frequency at which the area obtained by integrating from the frequency set by the frequency setting circuit to the end point of the allowable frequency region is halved Since it is detected as a point frequency, even if the frequency and amplitude waveforms are multi-peaked and the position of the peak with the largest amplitude value fluctuates greatly with time, it is hardly affected by the fluctuation, and the reference point is stable. The frequency can be detected, and blur and blur of the image can be removed.

【0022】また、区分周波数分析後の周波数と振幅の
波形を周波数方向に補間し、補間された波形を設定した
閾値に基づいて切り出すので、閾値以下の不要なノイズ
成分を除去することができる。このようにして得られた
波形を切り出した周波数領域の開始点から周波数設定回
路で設定した周波数まで積分して得られた面積が半分に
なる時の周波数と、周波数設定回路で設定した周波数か
ら切り出した周波数領域の終了点まで積分して得られた
面積が半分になる時の周波数との中心周波数を基準点周
波数として検出するので、周波数と振幅の波形が多峰
で、かつ振幅値が最大となる峰の位置が時間毎に大きく
変動するような場合でも変動の影響を受けにくく、安定
して基準点周波数を検出することができ、画像のぼけや
にじみを除去することができる。
Further, since the frequency and amplitude waveforms after the division frequency analysis are interpolated in the frequency direction and the interpolated waveform is cut out based on the set threshold value, unnecessary noise components below the threshold value can be removed. The waveform obtained in this way is cut out from the frequency at which the area obtained by integrating from the starting point of the cut-out frequency domain to the frequency set by the frequency setting circuit becomes half, and the frequency set by the frequency setting circuit The center frequency of the frequency at which the area obtained by integrating up to the end point of the frequency domain that is obtained by half is detected as the reference point frequency, so that the frequency and amplitude waveforms are multi-peaked and the amplitude value is the maximum. Even when the position of a peak fluctuates greatly with time, it is hardly affected by the fluctuation, the reference point frequency can be detected stably, and blur and blur of an image can be removed.

【0023】また、区分周波数分析後の周波数と振幅の
波形を周波数方向に補間し、補間された波形について各
振幅値のべき乗を計算するので、S/N比を向上させる
ことができる。このようにして得られた波形を許容され
る周波数領域の開始点から周波数設定回路で設定した周
波数まで積分して得られた面積が半分になる時の周波数
と、周波数設定回路で設定した周波数から許容される周
波数領域の終了点まで積分して得られた面積が半分にな
る時の周波数との中心周波数を基準点周波数として検出
するので、周波数と振幅の波形が多峰で、かつ振幅値が
最大となる峰の位置が時間毎に大きく変動するような場
合でも変動の影響を受けにくく、安定して基準点周波数
を検出することができ、画像のぼけやにじみを除去する
ことができる。
Further, since the waveform of the frequency and the amplitude after the division frequency analysis is interpolated in the frequency direction and the power of each amplitude value is calculated for the interpolated waveform, the S / N ratio can be improved. The frequency obtained when the area obtained by integrating the waveform obtained in this way from the start point of the allowable frequency region to the frequency set by the frequency setting circuit is halved, and the frequency set by the frequency setting circuit Since the center frequency with the frequency at which the area obtained by integrating up to the end point of the allowable frequency region is halved is detected as the reference point frequency, the waveform of the frequency and amplitude is multi-peak, and the amplitude value is Even when the position of the maximum peak fluctuates greatly with time, it is hardly affected by the fluctuation, the reference point frequency can be detected stably, and blur and blur of the image can be removed.

【0024】また、区分周波数分析後の周波数と振幅の
波形を周波数方向に補間し、補間された波形について各
振幅値のべき乗を計算するので、S/N比を向上させる
ことができる。更に、べき乗計算された波形を設定した
閾値に基づいて切り出すので、閾値以下の不要なノイズ
成分を除去することができる。このようにして得られた
波形を切り出した周波数領域の開始点から周波数設定回
路で設定した周波数まで積分して得られた面積が半分に
なる時の周波数と、周波数設定回路で設定した周波数か
ら切り出した周波数領域の終了点まで積分して得られた
面積が半分になる時の周波数との中心周波数を基準点周
波数として検出するので、周波数と振幅の波形が多峰
で、かつ振幅値が最大となる峰の位置が時間毎に大きく
変動するような場合でも変動の影響を受けにくく、安定
して基準点周波数を検出することができ、画像のぼけや
にじみを除去することができる。
Further, since the waveform of the frequency and the amplitude after the divided frequency analysis is interpolated in the frequency direction and the power of each amplitude value is calculated for the interpolated waveform, the S / N ratio can be improved. Furthermore, since the power-calculated waveform is cut out based on the set threshold value, unnecessary noise components below the threshold value can be removed. The waveform obtained in this way is cut out from the frequency at which the area obtained by integrating from the starting point of the cut-out frequency domain to the frequency set by the frequency setting circuit becomes half, and the frequency set by the frequency setting circuit The center frequency of the frequency at which the area obtained by integrating up to the end point of the frequency domain that is obtained by half is detected as the reference point frequency, so that the frequency and amplitude waveforms are multi-peaked and the amplitude value is the maximum. Even when the position of a peak fluctuates greatly with time, it is hardly affected by the fluctuation, the reference point frequency can be detected stably, and blur and blur of an image can be removed.

【0025】また、区分周波数分析後の周波数と振幅の
波形を周波数方向に補間し、補間された波形を設定した
閾値に基づいて切り出すので、閾値以下の不要なノイズ
成分を除去することができる。更に、切り出された波形
について各振幅値のべき乗を計算するので、計算量の削
減とS/N比の向上が可能である。このようにして得ら
れた波形を切り出した周波数領域の開始点から周波数設
定回路で設定した周波数まで積分して得られた面積が半
分になる時の周波数と、周波数設定回路で設定した周波
数から切り出した周波数領域の終了点まで積分して得ら
れた面積が半分になる時の周波数との中心周波数を基準
点周波数として検出するので、周波数と振幅の波形が多
峰で、かつ振幅値が最大となる峰の位置が時間毎に大き
く変動するような場合でも変動の影響を受けにくく、安
定して基準点周波数を検出することができ、画像のぼけ
やにじみを除去することができる。
Further, since the frequency and amplitude waveforms after the division frequency analysis are interpolated in the frequency direction and the interpolated waveform is cut out based on the set threshold value, unnecessary noise components below the threshold value can be removed. Furthermore, since the power of each amplitude value is calculated for the cut-out waveform, the calculation amount can be reduced and the S / N ratio can be improved. The waveform obtained in this way is cut out from the frequency at which the area obtained by integrating from the starting point of the cut-out frequency domain to the frequency set by the frequency setting circuit becomes half, and the frequency set by the frequency setting circuit The center frequency of the frequency at which the area obtained by integrating up to the end point of the frequency domain that is obtained by half is detected as the reference point frequency, so that the frequency and amplitude waveforms are multi-peaked and the amplitude value is the maximum. Even when the position of a peak fluctuates greatly with time, it is hardly affected by the fluctuation, the reference point frequency can be detected stably, and blur and blur of an image can be removed.

【0026】[0026]

【実施例】【Example】

実施例1 図1は、図6に示すレーダ信号処理装置におけるこの発
明の位相補正部の一実施例を示すものである。図におい
て、4,RS,RG及びRDは図6及び図7と同じであ
る。また、図において、8,9,11,12,13及び
GSは図7と同じである。14は区分周波数分析回路9
で得られた周波数と振幅の波形を周波数方向に補間する
補間回路、15は周波数と振幅の波形について許容され
る周波数領域内に基準となる周波数を設定する周波数設
定回路、16は周波数と振幅の波形を許容される周波数
領域の開始点から周波数設定回路15で設定した周波数
まで積分して面積を算出する第1の積分回路、17は第
1の積分回路16で算出された面積が半分になる時の周
波数を算出する第1の面積中心算出回路、18は周波数
と振幅の波形を周波数設定回路15で設定した周波数か
ら許容される周波数領域の終了点まで積分して面積を算
出する第2の積分回路、19は第2の積分回路18で算
出された面積が半分になる時の周波数を算出する第2の
面積中心算出回路、20は第1の面積中心算出回路17
で算出された周波数と第2の面積中心算出回路19で算
出された周波数の中心周波数を検出する中心周波数検出
回路である。
Embodiment 1 FIG. 1 shows an embodiment of the phase correction section of the present invention in the radar signal processing device shown in FIG. In the figure, 4, RS, RG and RD are the same as in FIGS. In the drawing, 8, 9, 11, 12, 13 and GS are the same as those in FIG. 14 is a division frequency analysis circuit 9
Interpolation circuit for interpolating the frequency and amplitude waveforms obtained in the above in the frequency direction, 15 is a frequency setting circuit for setting a reference frequency within a frequency region allowed for the frequency and amplitude waveforms, and 16 is a frequency and amplitude The first integration circuit 17 calculates the area by integrating the waveform from the start point of the allowable frequency region to the frequency set by the frequency setting circuit 15, and the area calculated by the first integration circuit 16 is halved. A first area center calculating circuit for calculating the frequency at the time, a second for calculating the area by integrating the frequency and amplitude waveforms from the frequency set by the frequency setting circuit 15 to the end point of the allowable frequency region. An integrating circuit 19 is a second area center calculating circuit for calculating a frequency when the area calculated by the second integrating circuit 18 is halved, and 20 is a first area center calculating circuit 17
Is a center frequency detection circuit that detects the center frequency between the frequency calculated in step (a) and the frequency calculated in the second area center calculation circuit 19.

【0027】次に、上記図1のように構成された位相補
正部4の動作について説明する。距離補正部3から入力
された時間による距離ずれが補正された目標受信信号R
S及びレーダと目標重心との初期距離RGはバッファ回
路8に格納され、目標受信信号RS及びレーダと目標重
心との初期距離RGにおける目標受信信号GSとして出
力される。このレーダと目標重心との初期距離RGにお
ける目標受信信号GSは区分周波数分析回路9で時間方
向に小区間で周波数分析され、補間回路14で周波数方
向に補間される。第1の積分回路16は補間回路14で
補間された波形を許容される周波数領域の開始点から周
波数設定回路15で設定した周波数まで積分して面積を
算出する。第2の積分回路18は補間回路14で補間さ
れた波形を周波数設定回路15で設定した周波数から許
容される周波数領域の終了点まで積分して面積を算出す
る。周波数設定回路15では、例えば、周波数と振幅の
波形について振幅値が最大になる時の周波数や周波数領
域の中心周波数など、許容される周波数領域内に基準と
なる周波数を設定する。第1の面積中心算出回路17は
第1の積分回路16で算出された面積が半分になる時の
周波数を算出し、第2の面積中心算出回路19は第2の
積分回路18で算出された面積が半分になる時の周波数
を算出する。
Next, the operation of the phase corrector 4 configured as shown in FIG. 1 will be described. The target reception signal R in which the distance shift due to the time input from the distance correction unit 3 has been corrected.
S and the initial distance RG between the radar and the target center of gravity are stored in the buffer circuit 8, and output as the target reception signal RS and the target reception signal GS at the initial distance RG between the radar and the target center of gravity. The target reception signal GS at the initial distance RG between the radar and the target center of gravity is frequency-analyzed in a small section in the time direction by the divided frequency analysis circuit 9, and is interpolated in the frequency direction by the interpolation circuit 14. The first integration circuit 16 calculates the area by integrating the waveform interpolated by the interpolation circuit 14 from the start point of the allowable frequency region to the frequency set by the frequency setting circuit 15. The second integration circuit 18 calculates the area by integrating the waveform interpolated by the interpolation circuit 14 from the frequency set by the frequency setting circuit 15 to the end point of the allowable frequency region. The frequency setting circuit 15 sets a reference frequency in an allowable frequency range, such as a frequency at which the amplitude value of the waveform of the frequency and the amplitude is maximized and a center frequency of the frequency range. The first area center calculation circuit 17 calculates the frequency when the area calculated by the first integration circuit 16 becomes half, and the second area center calculation circuit 19 calculates the frequency by the second integration circuit 18. Calculate the frequency when the area becomes half.

【0028】中心周波数検出回路20は第1の面積中心
算出回路17で算出された周波数と第2の面積中心算出
回路19で算出された周波数の中心周波数を基準点周波
数として検出した後、平滑化回路11に出力する。平滑
化回路11は基準点周波数の時間方向に対する軌跡を平
滑化し、平滑化された軌跡から位相補正量算出回路12
で位相補正量を算出する。位相補正回路13は位相補正
量算出回路12で算出した位相補正量を用いてバッファ
回路8から出力された時間による距離ずれが補正された
目標受信信号RSの位相を補正し、時間による位相ずれ
が補正された目標受信信号RDとして周波数分析部5へ
出力する。
The center frequency detection circuit 20 detects the center frequency of the frequency calculated by the first area center calculation circuit 17 and the center frequency calculated by the second area center calculation circuit 19 as a reference point frequency, and then performs smoothing. Output to the circuit 11. The smoothing circuit 11 smoothes the locus of the reference point frequency in the time direction, and calculates a phase correction amount calculating circuit 12 from the smoothed locus.
Is used to calculate the phase correction amount. The phase correction circuit 13 uses the phase correction amount calculated by the phase correction amount calculation circuit 12 to correct the phase of the target reception signal RS output from the buffer circuit 8 and corrected for the distance shift due to time. The signal is output to the frequency analysis unit 5 as the corrected target reception signal RD.

【0029】次に、上記図1のように構成された位相補
正部4を図8及び図9を用いて説明する。図8は位相補
正部4の処理方法、図9は区分周波数分析後の波形及び
データ補間後の波形を示した図である。時間による距離
ずれが補正された目標受信信号RSをSi,j (ここでi
はレンジビン番号、jはパルスヒット番号、i,jは自
然数である。)、レーダと目標重心との初期距離RGの
存在するレンジビン番号をrと定義するとレーダと目標
重心との初期距離RGにおける目標受信信号GSはS
r,j と表され、図8(a)のような波形が得られる。S
r,j に対し、区分周波数分析回路9で時間方向(パルス
ヒット方向)に小区間で周波数分析を行うと図8(b)
のような波形が得られ、周波数fm と振幅Ak m の関係
(ここでkは区分周波数分析番号、mは周波数ビン番
号、k,mは自然数である。)は“数1”で表される。
Next, the phase corrector 4 configured as shown in FIG. 1 will be described with reference to FIGS. FIG. 8 is a diagram illustrating a processing method of the phase correction unit 4, and FIG. 9 is a diagram illustrating a waveform after the segmented frequency analysis and a waveform after the data interpolation. The target reception signal RS whose distance deviation due to time has been corrected is represented by S i, j (where i
Is a range bin number, j is a pulse hit number, and i and j are natural numbers. If the range bin number where the initial distance RG between the radar and the target center of gravity exists is defined as r, the target reception signal GS at the initial distance RG between the radar and the target center of gravity becomes S
r, j, and a waveform as shown in FIG. 8A is obtained. S
FIG. 8 (b) shows that the divided frequency analysis circuit 9 performs frequency analysis on small sections in the time direction (pulse hit direction) for r and j .
Waveform is obtained as the relation (k is here divided frequency analysis number, m is a frequency bin number, k, m is a natural number.) Of the frequency f m and an amplitude A k m tables in "Number 1" Is done.

【0030】各区分周波数分析番号kにおいて周波数f
m と振幅Ak m の波形は図9(a)のように表され、補
間回路14で周波数方向に補間すると図9(b)のよう
な波形が得られる。ここで、得られた波形を第1の積分
回路16で許容される周波数領域の開始点f1 から周波
数設定回路15で設定した周波数ck まで積分すると面
積Uk は“数5”で算出される。
The frequency f at each divisional frequency analysis number k
m and the waveform of the amplitude A k m is expressed as shown in FIG. 9 (a), interpolated in the frequency direction is waveform as shown in FIG. 9 (b) is obtained by the interpolation circuit 14. Here, when the obtained waveform is integrated from the start point f 1 of the frequency region permitted by the first integration circuit 16 to the frequency c k set by the frequency setting circuit 15, the area U k is calculated by “Equation 5”. You.

【0031】[0031]

【数5】 (Equation 5)

【0032】また、得られた波形を第2の積分回路18
で周波数設定回路15で設定した周波数ck から許容さ
れる周波数領域の終了点f1 まで積分すると面積Vk
“数6”で算出される。
The obtained waveform is converted to a second integrating circuit 18.
By integrating from the frequency ck set by the frequency setting circuit 15 to the end point f 1 of the allowable frequency region, the area V k is calculated by “Equation 6”.

【0033】[0033]

【数6】 (Equation 6)

【0034】第1の面積中心算出回路17は、面積Uk
が半分になる時の周波数を“数7”を満足するように算
出し、その時の周波数をXk とする。
The first area center calculation circuit 17 calculates the area U k
Is calculated so as to satisfy “Equation 7”, and the frequency at that time is defined as X k .

【0035】[0035]

【数7】 (Equation 7)

【0036】第2の面積中心算出回路19は、面積Vk
が半分になる時の周波数を“数8”を満足するように算
出し、その時の周波数をYk とする。
The second area center calculating circuit 19 calculates the area V k
Is calculated so as to satisfy “Equation 8”, and the frequency at that time is defined as Y k .

【0037】[0037]

【数8】 (Equation 8)

【0038】中心周波数検出回路20は、周波数Xk
び周波数Yk を用いて、“数9”で基準点周波数fk
算出する。
The center frequency detection circuit 20 calculates the reference point frequency f k by using “Expression 9” using the frequency X k and the frequency Y k .

【0039】[0039]

【数9】 (Equation 9)

【0040】時間tk と基準点周波数fk の関係は図8
(c)のプロットのようになり、平滑化回路11で平滑
化を行うと図8(c)の実線のような波形が得られ、時
間tk と周波数f’k の関係は“数2”で表される。
The relationship between the time t k and the reference point frequency f k is shown in FIG.
Look like the plot of (c), the smoothing circuit 11 smoothes the waveform shown by the solid line shown in FIG. 8 (c) is obtained, the relationship between the time t k and the frequency f 'k' number 2 ' It is represented by

【0041】位相補正量算出回路12は、位相補正量W
j を“数3”で算出する。位相補正回路13は、Si,j
の位相を位相補正量Wj を用いて“数4”で補正する。
但し、時間による位相ずれが補正された目標受信信号R
DをS’i,j と定義する。
The phase correction amount calculation circuit 12 calculates the phase correction amount W
j is calculated by “Equation 3”. The phase correction circuit 13 calculates S i, j
It is corrected by "Number 4" by using a phase correction amount W j of the phase.
However, the target reception signal R in which the phase shift due to time is corrected
D is defined as S ′ i, j .

【0042】実施例2 図2に示される実施例では、上記実施例1における補間
回路14の直後に、補間回路14で補間された周波数と
振幅の波形についてデータを切り出すデータ切り出し回
路21と、データを切り出す際に必要な閾値を設定する
閾値設定回路22を設けている。このような実施態様に
よれば、補間された周波数と振幅の波形を設定した閾値
に基づいて切り出すので、閾値以下の不要なノイズ成分
を除去することができる。
Embodiment 2 In the embodiment shown in FIG. 2, immediately after the interpolation circuit 14 in the first embodiment, a data extraction circuit 21 for extracting data on the frequency and amplitude waveforms interpolated by the interpolation circuit 14, Is provided with a threshold setting circuit 22 for setting a threshold necessary for cutting out. According to such an embodiment, since the interpolated frequency and amplitude waveforms are cut out based on the set threshold value, unnecessary noise components below the threshold value can be removed.

【0043】次に、上記図2のように構成された位相補
正部4の動作について説明する。距離補正部3から入力
された時間による距離ずれが補正された目標受信信号R
S及びレーダと目標重心との初期距離RGはバッファ回
路8に格納され、目標受信信号RS及びレーダと目標重
心との初期距離RGにおける目標受信信号GSとして出
力される。このレーダと目標重心との初期距離RGにお
ける目標受信信号GSは区分周波数分析回路9で時間方
向に小区間で周波数分析され、補間回路14で周波数方
向に補間された後、データ切り出し回路21で各振幅値
が閾値設定回路22で設定した閾値を越える領域につい
てデータが切り出される。閾値設定回路22では、例え
ば、メインローブレベルとサイドローブレベルの間に一
定の閾値を設定する固定スレッショルドやアダプティブ
に閾値を設定するCFAR(Constant Fal
se Alarm Rate)等を用いて、各区分周波
数分析後の波形毎に閾値を設定する。第1の積分回路1
6ではデータ切り出し回路21で切り出された波形を切
り出した周波数領域の開始点から周波数設定回路15で
設定した周波数まで積分して面積を算出する。第2の積
分回路18ではデータ切り出し回路21で切り出された
波形を周波数設定回路15で設定した周波数から切り出
した周波数領域の終了点まで積分して面積を算出する。
周波数設定回路15では、例えば、周波数と振幅の波形
について振幅値が最大になる時の周波数や周波数領域の
中心周波数など、許容される周波数領域内に基準となる
周波数を設定する。第1の面積中心算出回路17では第
1の積分回路16で算出された面積が半分になる時の周
波数を算出、第2の面積中心算出回路19では第2の積
分回路18で算出された面積が半分になる時の周波数を
算出する。中心周波数検出回路20では第1の面積中心
算出回路17で算出された周波数と第2の面積中心算出
回路19で算出された周波数の中心周波数を基準点周波
数として検出した後、平滑化回路11に出力する。平滑
化回路11では基準点周波数の時間方向に対する軌跡を
平滑化し、平滑化された軌跡から位相補正量算出回路1
2で位相補正量を算出する。位相補正回路13は位相補
正量算出回路12で算出した位相補正量を用いてバッフ
ァ回路8から出力された時間による距離ずれが補正され
た目標受信信号RSの位相を補正し、時間による位相ず
れが補正された目標受信信号RDとして周波数分析部5
へ出力する。
Next, the operation of the phase corrector 4 configured as shown in FIG. 2 will be described. The target reception signal R in which the distance shift due to the time input from the distance correction unit 3 has been corrected.
S and the initial distance RG between the radar and the target center of gravity are stored in the buffer circuit 8, and output as the target reception signal RS and the target reception signal GS at the initial distance RG between the radar and the target center of gravity. The target reception signal GS at the initial distance RG between the radar and the target center of gravity is frequency-analyzed in a small section in the time direction by the divided frequency analysis circuit 9 and interpolated in the frequency direction by the interpolation circuit 14. Data is cut out for a region where the amplitude value exceeds the threshold set by the threshold setting circuit 22. In the threshold setting circuit 22, for example, a fixed threshold for setting a constant threshold between the main lobe level and the side lobe level, or a CFAR (Constant Fal) for adaptively setting a threshold.
For example, a threshold is set for each waveform after the division frequency analysis using “se Alarm Rate” or the like. First integration circuit 1
In step S6, an area is calculated by integrating the waveform extracted by the data extraction circuit 21 from the start point of the extracted frequency region to the frequency set by the frequency setting circuit 15. The second integration circuit 18 calculates the area by integrating the waveform cut out by the data cut-out circuit 21 from the frequency set by the frequency setting circuit 15 to the end point of the cut-out frequency region.
The frequency setting circuit 15 sets a reference frequency in an allowable frequency range, such as a frequency at which the amplitude value of the waveform of the frequency and the amplitude is maximized and a center frequency of the frequency range. The first area center calculation circuit 17 calculates the frequency when the area calculated by the first integration circuit 16 is reduced to half, and the second area center calculation circuit 19 calculates the area calculated by the second integration circuit 18. Calculate the frequency at which is reduced by half. The center frequency detection circuit 20 detects the center frequency of the frequency calculated by the first area center calculation circuit 17 and the center frequency of the frequency calculated by the second area center calculation circuit 19 as a reference point frequency. Output. The smoothing circuit 11 smoothes the locus of the reference point frequency in the time direction, and calculates the phase correction amount calculating circuit 1 from the smoothed locus.
In step 2, the amount of phase correction is calculated. The phase correction circuit 13 uses the phase correction amount calculated by the phase correction amount calculation circuit 12 to correct the phase of the target reception signal RS output from the buffer circuit 8 and corrected for the distance shift due to time. The frequency analysis unit 5 determines the corrected target reception signal RD as
Output to

【0044】次に、上記図2のように構成された位相補
正部4を図9及び図10を用いて説明する。図10はデ
ータ補間した波形の切り出し後の波形を示した図であ
る。各区分周波数分析番号kにおいて周波数fm と振幅
k m の波形は図9(a)のように表され、補間回路1
4で周波数方向に補間すると図9(b)のような波形が
得られ、データ切り出し回路21で各振幅値が閾値設定
回路22で設定した閾値(例えば、図9(b)中の一点
破線u)を越える周波数領域[aK ,bk ](ここでa
k は周波数領域の開始点、bk は周波数領域の終了点を
表す。)についてデータの切り出しを行うと図10のよ
うな波形が得られる。ここで、得られた波形を第1の積
分回路16で周波数領域の開始点ak から周波数設定回
路15で設定した周波数ck まで積分すると面積Uk
“数10”で算出される。
Next, the phase corrector 4 configured as shown in FIG. 2 will be described with reference to FIGS. FIG. 10 is a diagram illustrating a waveform after the data interpolation is performed on the waveform. Waveform of the frequency f m and an amplitude A k m at each division frequency analysis number k is expressed as in FIG. 9 (a), the interpolation circuit 1
4, a waveform as shown in FIG. 9B is obtained, and each amplitude value is set by the data cutout circuit 21 to the threshold value set by the threshold value setting circuit 22 (for example, the dashed line u in FIG. 9B). ) Beyond the frequency range [a K , b k ] (where a
k represents the start point of the frequency domain, and b k represents the end point of the frequency domain. When the data is cut out for ()), a waveform as shown in FIG. 10 is obtained. Here, when the obtained waveform is integrated by the first integration circuit 16 from the start point a k of the frequency domain to the frequency c k set by the frequency setting circuit 15, the area U k is calculated by “Equation 10”.

【0045】[0045]

【数10】 (Equation 10)

【0046】また、得られた波形を第2の積分回路18
で周波数設定回路15で設定した周波数ck から周波数
領域の終了点bk まで積分すると面積Vk は“数11”
で算出される。
The obtained waveform is converted to a second integrating circuit 18.
When integrating from the frequency c k set by the frequency setting circuit 15 to the end point b k of the frequency domain, the area V k becomes “Equation 11”.
Is calculated.

【0047】[0047]

【数11】 [Equation 11]

【0048】第1の面積中心算出回路17は、面積Uk
が半分になる時の周波数を“数12”を満足するように
算出し、その時の周波数をXk とする。
The first area center calculating circuit 17 calculates the area U k
Is calculated so as to satisfy “Equation 12”, and the frequency at that time is defined as X k .

【0049】[0049]

【数12】 (Equation 12)

【0050】第2の面積中心算出回路19は、面積Vk
が半分になる時の周波数を“数8”を満足するように算
出し、その時の周波数をYk とする。
The second area center calculating circuit 19 calculates the area V k
Is calculated so as to satisfy “Equation 8”, and the frequency at that time is defined as Y k .

【0051】中心周波数検出回路20は、周波数Xk
び周波数Yk を用いて、“数9”で基準点周波数fk
算出する。
The center frequency detecting circuit 20 calculates the reference point frequency f k by using “Equation 9” using the frequency X k and the frequency Y k .

【0052】実施例3 図3に示される実施例では、上記実施例1における補間
回路14の直後に、補間回路14で補間された周波数と
振幅の波形について各振幅値のべき乗を計算するべき乗
計算回路23を設けている。このような実施態様によれ
ば、補間された周波数と振幅の波形について各振幅値の
べき乗を計算するので、S/N比を向上させることがで
きる。
Embodiment 3 In the embodiment shown in FIG. 3, immediately after the interpolation circuit 14 in the first embodiment, a power calculation for calculating the power of each amplitude value for the frequency and amplitude waveforms interpolated by the interpolation circuit 14 is performed. A circuit 23 is provided. According to such an embodiment, since the power of each amplitude value is calculated for the interpolated frequency and amplitude waveforms, the S / N ratio can be improved.

【0053】次に、上記図3のように構成された位相補
正部4の動作について説明する。距離補正部3から入力
された時間による距離ずれが補正された目標受信信号R
S及びレーダと目標重心との初期距離RGはバッファ回
路8に格納され、目標受信信号RS及びレーダと目標重
心との初期距離RGにおける目標受信信号GSとして出
力される。このレーダと目標重心との初期距離RGにお
ける目標受信信号GSは区分周波数分析回路9で時間方
向に小区間で周波数分析され、補間回路14で周波数方
向に補間された後、べき乗計算回路23で各振幅値のべ
き乗が計算される。第1の積分回路16ではべき乗計算
回路23で計算された波形を許容される周波数領域の開
始点から周波数設定回路15で設定した周波数まで積分
して面積を算出する。第2の積分回路18ではべき乗計
算回路23で計算された波形を周波数設定回路15で設
定した周波数から許容される周波数領域の終了点まで積
分して面積を算出する。周波数設定回路15では、例え
ば、周波数と振幅の波形について振幅値が最大になる時
の周波数や周波数領域の中心周波数など、許容される周
波数領域内に基準となる周波数を設定する。第1の面積
中心算出回路17では第1の積分回路16で算出された
面積が半分になる時の周波数を算出、第2の面積中心算
出回路19では第2の積分回路18で算出された面積が
半分になる時の周波数を算出する。中心周波数検出回路
20では第1の面積中心算出回路17で算出された周波
数と第2の面積中心算出回路19で算出された周波数の
中心周波数を基準点周波数として検出した後、平滑化回
路11に出力する。平滑化回路11では基準点周波数の
時間方向に対する軌跡を平滑化し、平滑化された軌跡か
ら位相補正量算出回路12で位相補正量を算出する。位
相補正回路13は位相補正量算出回路12で算出した位
相補正量を用いてバッファ回路8から出力された時間に
よる距離ずれが補正された目標受信信号RSの位相を補
正し、時間による位相ずれが補正された目標受信信号R
Dとして周波数分析部5へ出力する。
Next, the operation of the phase corrector 4 configured as shown in FIG. 3 will be described. The target reception signal R in which the distance shift due to the time input from the distance correction unit 3 has been corrected.
S and the initial distance RG between the radar and the target center of gravity are stored in the buffer circuit 8, and output as the target reception signal RS and the target reception signal GS at the initial distance RG between the radar and the target center of gravity. The target reception signal GS at the initial distance RG between the radar and the target center of gravity is subjected to frequency analysis in a small section in the time direction by the divided frequency analysis circuit 9, interpolated in the frequency direction by the interpolation circuit 14, and A power of the amplitude value is calculated. The first integration circuit 16 integrates the waveform calculated by the exponentiation calculation circuit 23 from the start point of the allowable frequency region to the frequency set by the frequency setting circuit 15 to calculate the area. The second integrating circuit 18 calculates the area by integrating the waveform calculated by the exponentiation calculating circuit 23 from the frequency set by the frequency setting circuit 15 to the end point of the allowable frequency region. The frequency setting circuit 15 sets a reference frequency in an allowable frequency range, such as a frequency at which the amplitude value of the waveform of the frequency and the amplitude is maximized and a center frequency of the frequency range. The first area center calculation circuit 17 calculates the frequency when the area calculated by the first integration circuit 16 is reduced to half, and the second area center calculation circuit 19 calculates the area calculated by the second integration circuit 18. Calculate the frequency at which is reduced by half. The center frequency detection circuit 20 detects the center frequency of the frequency calculated by the first area center calculation circuit 17 and the center frequency of the frequency calculated by the second area center calculation circuit 19 as a reference point frequency. Output. The smoothing circuit 11 smoothes the locus of the reference point frequency in the time direction, and the phase correction amount calculating circuit 12 calculates the phase correction amount from the smoothed locus. The phase correction circuit 13 uses the phase correction amount calculated by the phase correction amount calculation circuit 12 to correct the phase of the target reception signal RS output from the buffer circuit 8 and corrected for the distance shift due to time. Corrected target received signal R
D is output to the frequency analysis unit 5.

【0054】次に、上記図3のように構成された位相補
正部4を図9及び図11を用いて説明する。図11はべ
き乗計算後の波形を示した図である。各区分周波数分析
番号kにおいて周波数fm と振幅Ak m の波形は図9
(a)のように表され、補間回路14で周波数方向に補
間すると図9(b)のような波形が得られ、べき乗計算
回路23で各振幅値のべき乗を計算すると図11のよう
な波形が得られる。ここで、得られた波形を第1の積分
回路16で許容される周波数領域の開始点f1 から周波
数設定回路15で設定した周波数ck まで積分すると面
積Uk は“数13”で算出される。
Next, the phase correcting section 4 configured as shown in FIG. 3 will be described with reference to FIGS. FIG. 11 is a diagram showing the waveform after the exponentiation calculation. Waveform of the frequency f m and an amplitude A k m at each division frequency analysis number k is 9
When the interpolation circuit 14 interpolates in the frequency direction, a waveform as shown in FIG. 9B is obtained. When the exponentiation circuit 23 calculates the power of each amplitude value, a waveform as shown in FIG. 11 is obtained. Is obtained. Here, when the obtained waveform is integrated from the start point f 1 of the frequency region permitted by the first integration circuit 16 to the frequency c k set by the frequency setting circuit 15, the area U k is calculated by “Equation 13”. You.

【0055】[0055]

【数13】 (Equation 13)

【0056】また、得られた波形を第2の積分回路18
で周波数設定回路15で設定した周波数ck から許容さ
れる周波数領域の終了点f1 まで積分すると面積Vk
“数14”で算出される。
The obtained waveform is converted to a second integrating circuit 18.
By integrating from the frequency c k set by the frequency setting circuit 15 to the end point f 1 of the allowable frequency region, the area V k is calculated by “Equation 14”.

【0057】[0057]

【数14】 [Equation 14]

【0058】第1の面積中心算出回路17は、面積Uk
が半分になる時の周波数を“数15”を満足するように
算出し、その時の周波数をXk とする。
The first area center calculating circuit 17 calculates the area U k
Is calculated so as to satisfy “Equation 15”, and the frequency at that time is defined as X k .

【0059】[0059]

【数15】 (Equation 15)

【0060】第2の面積中心算出回路19は、面積Vk
が半分になる時の周波数を“数16”を満足するように
算出し、その時の周波数をYk とする。
The second area center calculating circuit 19 calculates the area V k
Is calculated to satisfy “Equation 16”, and the frequency at that time is defined as Y k .

【0061】[0061]

【数16】 (Equation 16)

【0062】中心周波数検出回路20は、周波数Xk
び周波数Yk を用いて、“数9”で基準点周波数fk
算出する。
The center frequency detecting circuit 20 calculates the reference point frequency f k by using “Equation 9” using the frequency X k and the frequency Y k .

【0063】実施例4 図4に示される実施例では、上記実施例1における補間
回路14の直後に、補間回路14で補間された周波数と
振幅の波形について各振幅値のべき乗を計算するべき乗
計算回路23と、べき乗計算回路23で計算された周波
数と振幅のべき乗の波形についてデータを切り出すデー
タ切り出し回路21と、データを切り出す際に必要な閾
値を設定する閾値設定回路22を設けている。このよう
な実施態様によれば、補間された周波数と振幅の波形に
ついて各振幅値のべき乗を計算するので、S/N比を向
上させることができる。更に、べき乗計算された波形を
設定した閾値に基づいて切り出すので、閾値以下の不要
なノイズ成分を除去することができる。
Fourth Embodiment In the embodiment shown in FIG. 4, immediately after the interpolation circuit 14 in the first embodiment, a power calculation is performed to calculate the power of each amplitude value for the frequency and amplitude waveforms interpolated by the interpolation circuit 14. A circuit 23, a data extracting circuit 21 for extracting data on the power and the waveform of the frequency and the amplitude calculated by the exponentiation calculating circuit 23, and a threshold setting circuit 22 for setting a threshold necessary for extracting the data are provided. According to such an embodiment, since the power of each amplitude value is calculated for the interpolated frequency and amplitude waveforms, the S / N ratio can be improved. Furthermore, since the power-calculated waveform is cut out based on the set threshold value, unnecessary noise components below the threshold value can be removed.

【0064】次に、上記図4のように構成された位相補
正部4の動作について説明する。距離補正部3から入力
された時間による距離ずれが補正された目標受信信号R
S及びレーダと目標重心との初期距離RGはバッファ回
路8に格納され、目標受信信号RS及びレーダと目標重
心との初期距離RGにおける目標受信信号GSとして出
力される。このレーダと目標重心との初期距離RGにお
ける目標受信信号GSは区分周波数分析回路9で時間方
向に小区間で周波数分析され、補間回路14で周波数方
向に補間された後、べき乗計算回路23で各振幅値のべ
き乗が計算され、データ切り出し回路21で各振幅値が
閾値設定回路22で設定した閾値を越える領域について
データが切り出される。閾値設定回路22では、例え
ば、メインローブレベルとサイドローブレベルの間に一
定の閾値を設定する固定スレッショルドやアダプティブ
に閾値を設定するCFAR(Constant Fal
seAlarm Rate)等を用いて、各区分周波数
分析後の波形毎に閾値を設定する。第1の積分回路16
ではデータ切り出し回路21で切り出された波形を切り
出した周波数領域の開始点から周波数設定回路15で設
定した周波数まで積分して面積を算出する。第2の積分
回路18ではデータ切り出し回路21で切り出された波
形を周波数設定回路15で設定した周波数から切り出し
た周波数領域の終了点まで積分して面積を算出する。周
波数設定回路15では、例えば、周波数と振幅の波形に
ついて振幅値が最大になる時の周波数や周波数領域の中
心周波数など、許容される周波数領域内に基準となる周
波数を設定する。第1の面積中心算出回路17では第1
の積分回路16で算出された面積が半分になる時の周波
数を算出、第2の面積中心算出回路19では第2の積分
回路18で算出された面積が半分になる時の周波数を算
出する。
Next, the operation of the phase corrector 4 configured as shown in FIG. 4 will be described. The target reception signal R in which the distance shift due to the time input from the distance correction unit 3 has been corrected.
S and the initial distance RG between the radar and the target center of gravity are stored in the buffer circuit 8, and output as the target reception signal RS and the target reception signal GS at the initial distance RG between the radar and the target center of gravity. The target reception signal GS at the initial distance RG between the radar and the target center of gravity is subjected to frequency analysis in a small section in the time direction by the divided frequency analysis circuit 9, interpolated in the frequency direction by the interpolation circuit 14, and The power of the amplitude value is calculated, and data is cut out by the data cutout circuit 21 for an area where each amplitude value exceeds the threshold set by the threshold setting circuit 22. In the threshold setting circuit 22, for example, a fixed threshold for setting a constant threshold between the main lobe level and the side lobe level, or a CFAR (Constant Fal) for adaptively setting a threshold.
For example, a threshold is set for each waveform after each section frequency analysis by using (seAlarm Rate) or the like. First integration circuit 16
In the above, the area cut out from the start point of the frequency region cut out from the waveform cut out by the data cutout circuit 21 to the frequency set by the frequency setting circuit 15 is calculated. The second integration circuit 18 calculates the area by integrating the waveform cut out by the data cut-out circuit 21 from the frequency set by the frequency setting circuit 15 to the end point of the cut-out frequency region. The frequency setting circuit 15 sets a reference frequency in an allowable frequency range, such as a frequency at which the amplitude value of the waveform of the frequency and the amplitude is maximized and a center frequency of the frequency range. In the first area center calculation circuit 17, the first
The frequency when the area calculated by the integrating circuit 16 is halved is calculated, and the second area center calculating circuit 19 calculates the frequency when the area calculated by the second integrating circuit 18 is halved.

【0065】中心周波数検出回路20では第1の面積中
心算出回路17で算出された周波数と第2の面積中心算
出回路19で算出された周波数の中心周波数を基準点周
波数として検出した後、平滑化回路11に出力する。平
滑化回路11では基準点周波数の時間方向に対する軌跡
を平滑化し、平滑化された軌跡から位相補正量算出回路
12で位相補正量を算出する。位相補正回路13は位相
補正量算出回路12で算出した位相補正量を用いてバッ
ファ回路8から出力された時間による距離ずれが補正さ
れた目標受信信号RSの位相を補正し、時間による位相
ずれが補正された目標受信信号RDとして周波数分析部
5へ出力する。
The center frequency detection circuit 20 detects the center frequency of the frequency calculated by the first area center calculation circuit 17 and the center frequency calculated by the second area center calculation circuit 19 as a reference point frequency, and then performs smoothing. Output to the circuit 11. The smoothing circuit 11 smoothes the locus of the reference point frequency in the time direction, and the phase correction amount calculating circuit 12 calculates the phase correction amount from the smoothed locus. The phase correction circuit 13 uses the phase correction amount calculated by the phase correction amount calculation circuit 12 to correct the phase of the target reception signal RS output from the buffer circuit 8 and corrected for the distance shift due to time. The signal is output to the frequency analysis unit 5 as the corrected target reception signal RD.

【0066】次に、上記図4のように構成された位相補
正部4を図9,図11及び図12を用いて説明する。図
12はべき乗計算した波形の切り出し後の波形を示した
図である。各区分周波数分析番号kにおいて周波数fm
と振幅Ak m の波形は図9(a)のように表され、補間
回路14で周波数方向に補間すると図9(b)のような
波形が得られ、べき乗計算回路23で各振幅値のべき乗
を計算すると図11のような波形が得られる。更に、デ
ータ切り出し回路21で各振幅値が閾値設定回路22で
設定した閾値(例えば、図11中の一点破線v)を越え
る周波数領域[aK ,bk ](ここでak は周波数領域
の開始点、bk は周波数領域の終了点を表す。)につい
てデータの切り出しを行うと図12のような波形が得ら
れる。ここで、得られた波形を第1の積分回路16で周
波数領域の開始点ak から周波数設定回路15で設定し
た周波数ck まで積分すると面積Uk は“数17”で算
出される。
Next, the phase corrector 4 configured as shown in FIG. 4 will be described with reference to FIGS. 9, 11 and 12. FIG. 12 is a diagram showing a waveform after the extraction of the power-calculated waveform. Frequency f m at each division frequency analysis numbers k
The waveform of the amplitude A k m is expressed as FIG. 9 (a), the waveform obtained as the interpolated in the frequency direction in the interpolation circuit 14 FIG. 9 (b), the a power calculating circuit 23 for amplitude values When a power is calculated, a waveform as shown in FIG. 11 is obtained. Further, a frequency region [a K , b k ] in which each amplitude value exceeds a threshold value (for example, a dashed line v in FIG. 11) set by the threshold value setting circuit 22 in the data extraction circuit 21 (where a k is a frequency region) When the data is cut out at the start point and b k indicates the end point of the frequency domain, a waveform as shown in FIG. 12 is obtained. Here, when the obtained waveform is integrated by the first integration circuit 16 from the start point a k of the frequency domain to the frequency c k set by the frequency setting circuit 15, the area U k is calculated by “Equation 17”.

【0067】[0067]

【数17】 [Equation 17]

【0068】また、得られた波形を第2の積分回路18
で周波数設定回路15で設定した周波数ck から周波数
領域の終了点bk まで積分すると面積Vk は“数18”
で算出される。
The obtained waveform is converted to a second integrating circuit 18.
When integrating from the frequency c k set by the frequency setting circuit 15 to the end point b k of the frequency domain, the area V k becomes “Equation 18”.
Is calculated.

【0069】[0069]

【数18】 (Equation 18)

【0070】第1の面積中心算出回路17は、面積Uk
が半分になる時の周波数を“数19”を満足するように
算出し、その時の周波数をXk とする。
The first area center calculation circuit 17 calculates the area U k
Is calculated so as to satisfy “Equation 19”, and the frequency at that time is defined as X k .

【0071】[0071]

【数19】 [Equation 19]

【0072】第2の面積中心算出回路19は、面積Vk
が半分になる時の周波数を“数16”を満足するように
算出し、その時の周波数をYk とする。中心周波数検出
回路20は、周波数Xk 及び周波数Yk を用いて、“数
9”で基準点周波数fk を算出する。
The second area center calculating circuit 19 calculates the area V k
Is calculated to satisfy “Equation 16”, and the frequency at that time is defined as Y k . The center frequency detection circuit 20 calculates the reference point frequency f k using “Expression 9” using the frequency X k and the frequency Y k .

【0073】実施例5 図5に示される実施例では、上記実施例1における補間
回路14の直後に、補間回路14で補間された周波数と
振幅の波形についてデータを切り出すデータ切り出し回
路21と、データを切り出す際に必要な閾値を設定する
閾値設定回路22と、データ切り出し回路21で切り出
された周波数と振幅の波形について各振幅値のべき乗を
計算するべき乗計算回路23を設けている。このような
実施態様によれば、補間された周波数と振幅の波形を設
定した閾値に基づいて切り出すので、閾値以下の不要な
ノイズ成分を除去することができる。更に、切り出され
た波形について各振幅値のべき乗を計算するので、計算
量の削減とS/N比の向上が可能である。
Fifth Embodiment In the embodiment shown in FIG. 5, immediately after the interpolation circuit 14 in the first embodiment, a data extraction circuit 21 for extracting data on the frequency and amplitude waveforms interpolated by the interpolation circuit 14, A threshold setting circuit 22 for setting a threshold necessary for extracting the data, and a power calculation circuit 23 for calculating the power of each amplitude value for the frequency and amplitude waveforms extracted by the data extraction circuit 21. According to such an embodiment, since the interpolated frequency and amplitude waveforms are cut out based on the set threshold value, unnecessary noise components below the threshold value can be removed. Furthermore, since the power of each amplitude value is calculated for the cut-out waveform, the calculation amount can be reduced and the S / N ratio can be improved.

【0074】次に、上記図5のように構成された位相補
正部4の動作について説明する。距離補正部3から入力
された時間による距離ずれが補正された目標受信信号R
S及びレーダと目標重心との初期距離RGはバッファ回
路8に格納され、目標受信信号RS及びレーダと目標重
心との初期距離RGにおける目標受信信号GSとして出
力される。このレーダと目標重心との初期距離RGにお
ける目標受信信号GSは区分周波数分析回路9で時間方
向に小区間で周波数分析され、補間回路14で周波数方
向に補間された後、データ切り出し回路21で各振幅値
が閾値設定回路22で設定した閾値を越える領域につい
てデータが切り出され、べき乗計算回路23で各振幅値
のべき乗が計算される。閾値設定回路22では、例え
ば、メインローブレベルとサイドローブレベルの間に一
定の閾値を設定する固定スレッショルドやアダプティブ
に閾値を設定するCFAR(Constant Fal
seAlarm Rate)等を用いて、各区分周波数
分析後の波形毎に閾値を設定する。第1の積分回路16
ではべき乗計算回路23で計算された波形を切り出した
周波数領域の開始点から周波数設定回路15で設定した
周波数まで積分して面積を算出する。第2の積分回路1
8ではべき乗計算回路23で計算された波形を周波数設
定回路15で設定した周波数から切り出した周波数領域
の終了点まで積分して面積を算出する。周波数設定回路
15では、例えば、周波数と振幅の波形について振幅値
が最大になる時の周波数や周波数領域の中心周波数な
ど、許容される周波数領域内に基準となる周波数を設定
する。
Next, the operation of the phase corrector 4 configured as shown in FIG. 5 will be described. The target reception signal R in which the distance shift due to the time input from the distance correction unit 3 has been corrected.
S and the initial distance RG between the radar and the target center of gravity are stored in the buffer circuit 8, and output as the target reception signal RS and the target reception signal GS at the initial distance RG between the radar and the target center of gravity. The target reception signal GS at the initial distance RG between the radar and the target center of gravity is frequency-analyzed in a small section in the time direction by the divided frequency analysis circuit 9 and interpolated in the frequency direction by the interpolation circuit 14. Data is cut out for a region where the amplitude value exceeds the threshold value set by the threshold value setting circuit 22, and the power calculation circuit 23 calculates the power of each amplitude value. In the threshold setting circuit 22, for example, a fixed threshold for setting a constant threshold between the main lobe level and the side lobe level, or a CFAR (Constant Fal) for adaptively setting a threshold.
For example, a threshold is set for each waveform after each section frequency analysis by using (seAlarm Rate) or the like. First integration circuit 16
Then, the waveform calculated by the exponentiation calculation circuit 23 is integrated from the start point of the cut-out frequency region to the frequency set by the frequency setting circuit 15 to calculate the area. Second integration circuit 1
In step 8, the area calculated by integrating the waveform calculated by the power calculation circuit 23 from the frequency set by the frequency setting circuit 15 to the end point of the frequency region cut out. The frequency setting circuit 15 sets a reference frequency in an allowable frequency range, such as a frequency at which the amplitude value of the waveform of the frequency and the amplitude is maximized and a center frequency of the frequency range.

【0075】第1の面積中心算出回路17では第1の積
分回路16で算出された面積が半分になる時の周波数を
算出、第2の面積中心算出回路19では第2の積分回路
18で算出された面積が半分になる時の周波数を算出す
る。中心周波数検出回路20では第1の面積中心算出回
路17で算出された周波数と第2の面積中心算出回路1
9で算出された周波数の中心周波数を基準点周波数とし
て検出した後、平滑化回路11に出力する。平滑化回路
11では基準点周波数の時間方向に対する軌跡を平滑化
し、平滑化された軌跡から位相補正量算出回路12で位
相補正量を算出する。位相補正回路13は位相補正量算
出回路12で算出した位相補正量を用いてバッファ回路
8から出力された時間による距離ずれが補正された目標
受信信号RSの位相を補正し、時間による位相ずれが補
正された目標受信信号RDとして周波数分析部5へ出力
する。
The first area center calculating circuit 17 calculates the frequency when the area calculated by the first integrating circuit 16 is reduced to half, and the second area center calculating circuit 19 calculates the frequency by the second integrating circuit 18. The frequency at which the calculated area is halved is calculated. In the center frequency detection circuit 20, the frequency calculated by the first area center calculation circuit 17 and the second area center calculation circuit 1
After detecting the center frequency of the frequency calculated in step 9 as a reference point frequency, the frequency is output to the smoothing circuit 11. The smoothing circuit 11 smoothes the locus of the reference point frequency in the time direction, and the phase correction amount calculating circuit 12 calculates the phase correction amount from the smoothed locus. The phase correction circuit 13 uses the phase correction amount calculated by the phase correction amount calculation circuit 12 to correct the phase of the target reception signal RS output from the buffer circuit 8 and corrected for the distance shift due to time. The signal is output to the frequency analysis unit 5 as the corrected target reception signal RD.

【0076】次に、上記図5のように構成された位相補
正部4を図9,図10及び図12を用いて説明する。図
12は切り出した波形のべき乗計算後の波形を示した図
である。各区分周波数分析番号kにおいて周波数fm
振幅Ak m の波形は図9(a)のように表され、補間回
路14で周波数方向に補間すると図9(b)のような波
形が得られ、データ切り出し回路21で各振幅値が閾値
設定回路22で設定した閾値(例えば、図9(b)中の
一点破線u)を越える周波数領域[aK ,bk](ここ
でak は周波数領域の開始点、bk は周波数領域の終了
点を表す。)についてデータの切り出しを行うと図10
のような波形が得られる。更に、べき乗計算回路23で
各振幅値のべき乗を計算すると図12のような波形が得
られる。ここで、得られた波形を第1の積分回路16で
周波数領域の開始点ak から周波数設定回路15で設定
した周波数ck まで積分すると面積Uk は“数17”で
算出される。
Next, the phase correcting section 4 configured as shown in FIG. 5 will be described with reference to FIGS. 9, 10 and 12. FIG. 12 is a diagram showing a waveform after exponentiation calculation of the extracted waveform. Waveform of the frequency f m and an amplitude A k m at each division frequency analysis number k is expressed as in FIG. 9 (a), when interpolation in the frequency direction is waveform as shown in FIG. 9 (b) obtained by the interpolation circuit 14 The frequency region [a K , b k ] (where a k is the frequency) where each amplitude value exceeds the threshold set by the threshold setting circuit 22 (for example, the one-dot broken line u in FIG. 9B) in the data cutout circuit 21 When the data is cut out for the start point of the area and b k indicates the end point of the frequency area, FIG.
Is obtained. Further, when the power calculation circuit 23 calculates the power of each amplitude value, a waveform as shown in FIG. 12 is obtained. Here, when the obtained waveform is integrated by the first integration circuit 16 from the start point a k of the frequency domain to the frequency c k set by the frequency setting circuit 15, the area U k is calculated by “Equation 17”.

【0077】また、得られた波形を第2の積分回路18
で周波数設定回路15で設定した周波数ck から周波数
領域の終了点bk まで積分すると面積Vk は“数18”
で算出される。第1の面積中心算出回路17は、面積U
k が半分になる時の周波数を“数19”を満足するよう
に算出し、その時の周波数をXk とする。
The obtained waveform is converted to a second integrating circuit 18.
When integrating from the frequency c k set by the frequency setting circuit 15 to the end point b k of the frequency domain, the area V k becomes “Equation 18”.
Is calculated. The first area center calculation circuit 17 calculates the area U
The frequency when k becomes half is calculated so as to satisfy “Equation 19”, and the frequency at that time is set to X k .

【0078】第2の面積中心算出回路19は、面積Vk
が半分になる時の周波数を“数16”を満足するように
算出し、その時の周波数をYk とする。中心周波数検出
回路20は、周波数Xk 及び周波数Yk を用いて、“数
9”で基準点周波数fk を算出する。
The second area center calculating circuit 19 calculates the area V k
Is calculated to satisfy “Equation 16”, and the frequency at that time is defined as Y k . The center frequency detection circuit 20 calculates the reference point frequency f k using “Expression 9” using the frequency X k and the frequency Y k .

【0079】[0079]

【発明の効果】この発明は以上説明した通り、区分周波
数分析後の波形の面積に着目し、周波数領域の開始点か
ら周波数設定回路で設定した周波数まで積分して得られ
た面積が半分になる時の周波数と、周波数設定回路で設
定した周波数から周波数領域の終了点まで積分して得ら
れた面積が半分になる時の周波数との中心周波数を基準
点周波数として検出するので、周波数と振幅の波形が多
峰で、かつ振幅値が最大となる峰の位置が時間毎に大き
く変動するような場合でも変動の影響を受けにくく、安
定して基準点周波数を検出することができ、画像のぼけ
やにじみを除去することができるという効果がある。
As described above, according to the present invention, the area obtained by integrating from the starting point of the frequency domain to the frequency set by the frequency setting circuit is reduced by focusing on the area of the waveform after the divided frequency analysis. The center frequency between the time frequency and the frequency when the area obtained by integrating from the frequency set by the frequency setting circuit to the end point of the frequency domain is halved is detected as the reference point frequency. Even when the position of the peak where the waveform is multi-peak and the amplitude value is maximum fluctuates greatly with time, it is hardly affected by the fluctuation, the reference point frequency can be detected stably, and the image blur There is an effect that bleeding can be removed.

【0080】また、この発明は区分周波数分析後の周波
数と振幅の波形を周波数方向に補間し、補間された波形
を設定した閾値に基づいて切り出すので、閾値以下の不
要なノイズ成分を除去することができる。このようにし
て得られた波形を切り出した周波数領域の開始点から周
波数設定回路で設定した周波数まで積分して得られた面
積が半分になる時の周波数と、周波数設定回路で設定し
た周波数から切り出した周波数領域の終了点まで積分し
て得られた面積が半分になる時の周波数との中心周波数
を基準点周波数として検出するので、周波数と振幅の波
形が多峰で、かつ振幅値が最大となる峰の位置が時間毎
に大きく変動するような場合でも変動の影響を受けにく
く、安定して基準点周波数を検出することができ、画像
のぼけやにじみを除去することができる。
Further, according to the present invention, the frequency and amplitude waveforms after the divided frequency analysis are interpolated in the frequency direction, and the interpolated waveform is cut out based on the set threshold value, so that unnecessary noise components below the threshold value are removed. Can be. The waveform obtained in this way is cut out from the frequency at which the area obtained by integrating from the starting point of the cut-out frequency domain to the frequency set by the frequency setting circuit becomes half, and the frequency set by the frequency setting circuit The center frequency of the frequency at which the area obtained by integrating up to the end point of the frequency domain that is obtained by half is detected as the reference point frequency, so that the frequency and amplitude waveforms are multi-peaked and the amplitude value is the maximum. Even when the position of a peak fluctuates greatly with time, it is hardly affected by the fluctuation, the reference point frequency can be detected stably, and blur and blur of an image can be removed.

【0081】この発明は、区分周波数分析後の周波数と
振幅の波形を周波数方向に補間し、補間された波形につ
いて各振幅値のべき乗を計算するので、S/N比を向上
させることができる。このようにして得られた波形を許
容される周波数領域の開始点から周波数設定回路で設定
した周波数まで積分して得られた面積が半分になる時の
周波数と、周波数設定回路で設定した周波数から許容さ
れる周波数領域の終了点まで積分して得られた面積が半
分になる時の周波数との中心周波数を基準点周波数とし
て検出するので、周波数と振幅の波形が多峰で、かつ振
幅値が最大となる峰の位置が時間毎に大きく変動するよ
うな場合でも変動の影響を受けにくく、安定して基準点
周波数を検出することができ、画像のぼけやにじみを除
去することができる。
According to the present invention, the frequency and amplitude waveforms after the division frequency analysis are interpolated in the frequency direction, and the power of each amplitude value is calculated for the interpolated waveform, so that the S / N ratio can be improved. The frequency obtained when the area obtained by integrating the waveform obtained in this way from the start point of the allowable frequency region to the frequency set by the frequency setting circuit is halved, and the frequency set by the frequency setting circuit Since the center frequency with the frequency at which the area obtained by integrating up to the end point of the allowable frequency region is halved is detected as the reference point frequency, the waveform of the frequency and amplitude is multi-peak, and the amplitude value is Even when the position of the maximum peak fluctuates greatly with time, it is hardly affected by the fluctuation, the reference point frequency can be detected stably, and blur and blur of the image can be removed.

【0082】また、この発明は区分周波数分析後の周波
数と振幅の波形を周波数方向に補間し、補間された波形
について各振幅値のべき乗を計算するので、S/N比を
向上させることができる。更に、べき乗計算された波形
を設定した閾値に基づいて切り出すので、閾値以下の不
要なノイズ成分を除去することができる。このようにし
て得られた波形を切り出した周波数領域の開始点から周
波数設定回路で設定した周波数まで積分して得られた面
積が半分になる時の周波数と、周波数設定回路で設定し
た周波数から切り出した周波数領域の終了点まで積分し
て得られた面積が半分になる時の周波数との中心周波数
を基準点周波数として検出するので、周波数と振幅の波
形が多峰で、かつ振幅値が最大となる峰の位置が時間毎
に大きく変動するような場合でも変動の影響を受けにく
く、安定して基準点周波数を検出することができ、画像
のぼけやにじみを除去することができる。
According to the present invention, the frequency and amplitude waveforms after the division frequency analysis are interpolated in the frequency direction, and the power of each amplitude value is calculated for the interpolated waveform, so that the S / N ratio can be improved. . Furthermore, since the power-calculated waveform is cut out based on the set threshold value, unnecessary noise components below the threshold value can be removed. The waveform obtained in this way is cut out from the frequency at which the area obtained by integrating from the starting point of the cut-out frequency domain to the frequency set by the frequency setting circuit becomes half, and the frequency set by the frequency setting circuit The center frequency of the frequency at which the area obtained by integrating up to the end point of the frequency domain that is obtained by half is detected as the reference point frequency, so that the frequency and amplitude waveforms are multi-peaked and the amplitude value is the maximum. Even when the position of a peak fluctuates greatly with time, it is hardly affected by the fluctuation, the reference point frequency can be detected stably, and blur and blur of an image can be removed.

【0083】この発明は、区分周波数分析後の周波数と
振幅の波形を周波数方向に補間し、補間された波形を設
定した閾値に基づいて切り出すので、閾値以下の不要な
ノイズ成分を除去することができる。更に、切り出され
た波形について各振幅値のべき乗を計算するので、計算
量の削減とS/N比の向上が可能である。このようにし
て得られた波形を切り出した周波数領域の開始点から周
波数設定回路で設定した周波数まで積分して得られた面
積が半分になる時の周波数と、周波数設定回路で設定し
た周波数から切り出した周波数領域の終了点まで積分し
て得られた面積が半分になる時の周波数との中心周波数
を基準点周波数として検出するので、周波数と振幅の波
形が多峰で、かつ振幅値が最大となる峰の位置が時間毎
に大きく変動するような場合でも変動の影響を受けにく
く、安定して基準点周波数を検出することができ、画像
のぼけやにじみを除去することができる。
According to the present invention, since the frequency and amplitude waveforms after the division frequency analysis are interpolated in the frequency direction and the interpolated waveform is cut out based on the set threshold, unnecessary noise components below the threshold can be removed. it can. Furthermore, since the power of each amplitude value is calculated for the cut-out waveform, the calculation amount can be reduced and the S / N ratio can be improved. The waveform obtained in this way is cut out from the frequency at which the area obtained by integrating from the starting point of the cut-out frequency domain to the frequency set by the frequency setting circuit becomes half, and the frequency set by the frequency setting circuit The center frequency of the frequency at which the area obtained by integrating up to the end point of the frequency domain that is obtained by half is detected as the reference point frequency, so that the frequency and amplitude waveforms are multi-peaked and the amplitude value is the maximum. Even when the position of a peak fluctuates greatly with time, it is hardly affected by the fluctuation, the reference point frequency can be detected stably, and blur and blur of an image can be removed.

【図面の簡単な説明】[Brief description of the drawings]

【図1】この発明の実施例1を示す位相補正部の構成図
である。
FIG. 1 is a configuration diagram of a phase correction unit according to a first embodiment of the present invention.

【図2】この発明の実施例2を示す位相補正部の構成図
である。
FIG. 2 is a configuration diagram of a phase correction unit according to a second embodiment of the present invention.

【図3】この発明の実施例3を示す位相補正部の構成図
である。
FIG. 3 is a configuration diagram of a phase correction unit according to a third embodiment of the present invention.

【図4】この発明の実施例4を示す位相補正部の構成図
である。
FIG. 4 is a configuration diagram of a phase correction unit according to a fourth embodiment of the present invention.

【図5】この発明の実施例5を示す位相補正部の構成図
である。
FIG. 5 is a configuration diagram of a phase correction unit according to a fifth embodiment of the present invention.

【図6】高分解能レーダ装置のレーダ信号処理装置の構
成図である。
FIG. 6 is a configuration diagram of a radar signal processing device of the high resolution radar device.

【図7】従来の位相補正部の構成図である。FIG. 7 is a configuration diagram of a conventional phase correction unit.

【図8】位相補正部の処理方法を示す図である。FIG. 8 is a diagram illustrating a processing method of a phase correction unit.

【図9】区分周波数分析後の波形及びデータ補間後の波
形を示す図である。
FIG. 9 is a diagram showing a waveform after a divided frequency analysis and a waveform after data interpolation.

【図10】データ補間した波形の切り出し後の波形を示
す図である。
FIG. 10 is a diagram showing a waveform of a data-interpolated waveform after clipping.

【図11】べき乗計算後の波形を示す図である。FIG. 11 is a diagram showing a waveform after exponentiation calculation.

【図12】べき乗計算した波形の切り出し後の波形又は
切り出した波形のべき乗計算後の波形を示す図である。
FIG. 12 is a diagram showing a waveform after a power-calculated waveform is cut out or a waveform after a power-out calculation of the cut-out waveform.

【符号の説明】[Explanation of symbols]

1 データインタフェース部 2 パルス圧縮部 3 距離補正部 4 位相補正部 5 周波数分析部 6 検波部 7 表示器インタフェース部 8 バッファ回路 9 区分周波数分析回路 10 振幅値最大検出回路 11 平滑化回路 12 位相補正量算出回路 13 位相補正回路 14 補間回路 15 周波数設定回路 16 第1の積分回路 17 第1の面積中心算出回路 18 第2の積分回路 19 第2の面積中心算出回路 20 中心周波数検出回路 21 データ切り出し回路 22 閾値設定回路 23 べき乗計算回路 SM レーダ装置から入力された目標受信信号 RS 距離補正部3で時間による距離ずれが補正された
目標受信信号 RG レーダと目標重心との初期距離 GS レーダと目標重心との初期距離RGにおける目標
受信信号 RD 位相補正部4で時間による位相ずれが補正された
目標受信信号 D 表示画像データ
Reference Signs List 1 data interface unit 2 pulse compression unit 3 distance correction unit 4 phase correction unit 5 frequency analysis unit 6 detection unit 7 display interface unit 8 buffer circuit 9 division frequency analysis circuit 10 maximum amplitude value detection circuit 11 smoothing circuit 12 phase correction amount Calculation circuit 13 Phase correction circuit 14 Interpolation circuit 15 Frequency setting circuit 16 First integration circuit 17 First area center calculation circuit 18 Second integration circuit 19 Second area center calculation circuit 20 Center frequency detection circuit 21 Data extraction circuit Reference Signs List 22 threshold setting circuit 23 exponentiation calculation circuit target reception signal input from SM radar device RS target reception signal whose distance deviation due to time has been corrected by distance correction unit 3 RG initial distance between radar and target center of gravity GS radar and target center of gravity Target reception signal RD at the initial distance RG of the Target reception signal D with phase shift corrected D Display image data

───────────────────────────────────────────────────── フロントページの続き (58)調査した分野(Int.Cl.7,DB名) G01S 7/00 - 7/42 G01S 13/00 - 13/96 ──────────────────────────────────────────────────続 き Continued on the front page (58) Field surveyed (Int. Cl. 7 , DB name) G01S 7 /00-7/42 G01S 13/00-13/96

Claims (5)

(57)【特許請求の範囲】(57) [Claims] 【請求項1】 レーダ装置から入力された目標受信信号
を内部で処理できるデータ形式に変換するデータインタ
フェース部と、上記データインタフェース部で変換され
た目標受信信号をパルス圧縮するパルス圧縮部と、上記
パルス圧縮部でパルス圧縮された目標受信信号の時間に
よる距離ずれを補正する距離補正部と、上記距離補正部
で補正された目標受信信号の時間による位相ずれを補正
する位相補正部と、上記位相補正部で補正された目標受
信信号のドップラ周波数を分離する周波数分析部と、上
記周波数分析部で周波数分析された目標受信信号の周波
数スペクトルを画像データに変換する検波部と、上記検
波部で得られた画像データと表示器のインタフェースを
調整し表示画像データを生成する表示器インタフェース
部とを備えたレーダ信号処理装置において、上記位相補
正部を上記距離補正部で補正された目標受信信号及びレ
ーダと目標重心との初期距離を格納するバッファ回路
と、上記バッファ回路より出力されたレーダと目標重心
との初期距離における目標受信信号を時間方向に小区間
で周波数分析する区分周波数分析回路と、上記区分周波
数分析回路で得られた周波数と振幅の波形を周波数方向
に補間する補間回路と、周波数と振幅の波形について許
容される周波数領域内に基準となる周波数を設定する周
波数設定回路と、周波数と振幅の波形を許容される周波
数領域の開始点から上記周波数設定回路で設定した周波
数まで積分して面積を算出する第1の積分回路と、上記
第1の積分回路で算出された面積が半分になる時の周波
数を算出する第1の面積中心算出回路と、周波数と振幅
の波形を上記周波数設定回路で設定した周波数から許容
される周波数領域の終了点まで積分して面積を算出する
第2の積分回路と、上記第2の積分回路で算出された面
積が半分になる時の周波数を算出する第2の面積中心算
出回路と、上記第1の面積中心算出回路で算出された周
波数と上記第2の面積中心算出回路で算出された周波数
の中心周波数を基準点周波数として検出する中心周波数
検出回路と、上記中心周波数検出回路で検出された基準
点周波数の時間方向に対する軌跡を平滑化する平滑化回
路と、上記平滑化回路で平滑化された基準点周波数の時
間方向に対する軌跡から位相補正量を算出する位相補正
量算出回路と、上記位相補正量算出回路で算出された位
相補正量を用いて上記バッファ回路から出力された目標
受信信号の位相を補正する位相補正回路とによって構成
したことを特徴とするレーダ信号処理装置。
A data interface unit for converting a target reception signal input from a radar device into a data format that can be internally processed; a pulse compression unit for pulse-compressing the target reception signal converted by the data interface unit; A distance correction unit that corrects a distance shift due to time of the target reception signal pulse-compressed by the pulse compression unit; a phase correction unit that corrects a phase shift due to time of the target reception signal corrected by the distance correction unit; A frequency analysis unit that separates the Doppler frequency of the target reception signal corrected by the correction unit, a detection unit that converts the frequency spectrum of the target reception signal frequency-analyzed by the frequency analysis unit into image data, and a detection unit that obtains the image data. Having a display interface unit that adjusts the interface of the displayed image data and the display to generate display image data In the signal processing device, the phase correction unit stores a target reception signal corrected by the distance correction unit and a buffer circuit that stores an initial distance between the radar and the target center of gravity, and a radar and a target center of gravity output from the buffer circuit. A division frequency analysis circuit for performing frequency analysis of a target reception signal at an initial distance in a small section in a time direction, an interpolation circuit for interpolating in the frequency direction a frequency and amplitude waveform obtained by the above-described division frequency analysis circuit, and a frequency and amplitude A frequency setting circuit for setting a reference frequency within a frequency range permitted for the waveform, and an area obtained by integrating the frequency and amplitude waveforms from the start point of the permitted frequency range to the frequency set by the frequency setting circuit. A first integration circuit for calculating, and a first area center calculation circuit for calculating a frequency when the area calculated by the first integration circuit is halved; A second integrating circuit for integrating the frequency and amplitude waveforms from the frequency set by the frequency setting circuit to the end point of the allowable frequency region to calculate an area; and calculating the area calculated by the second integrating circuit. A second area center calculating circuit for calculating the frequency when the frequency is reduced to half; a frequency calculated by the first area center calculating circuit and a center frequency of the frequency calculated by the second area center calculating circuit as a reference; A center frequency detection circuit that detects the point frequency, a smoothing circuit that smoothes a trajectory in the time direction of the reference point frequency detected by the center frequency detection circuit, and a reference point frequency that is smoothed by the smoothing circuit. A phase correction amount calculation circuit for calculating a phase correction amount from a trajectory in the time direction, and a target reception amount output from the buffer circuit using the phase correction amount calculated by the phase correction amount calculation circuit. The radar signal processing apparatus characterized by being configured by a phase correction circuit for correcting the phase of the signal.
【請求項2】 上記位相補正部に上記補間回路で補間さ
れた周波数と振幅の波形についてデータを切り出すデー
タ切り出し回路と、上記データ切り出し回路でデータを
切り出す際に必要な閾値を設定する閾値設定回路を設け
たことを特徴とする請求項1記載のレーダ信号処理装
置。
2. A data extraction circuit for extracting data on a frequency and amplitude waveform interpolated by the interpolation circuit in the phase correction unit, and a threshold setting circuit for setting a threshold necessary for extracting data by the data extraction circuit The radar signal processing device according to claim 1, further comprising:
【請求項3】 上記位相補正部に上記補間回路で補間さ
れた周波数と振幅の波形について各振幅値のべき乗を計
算するべき乗計算回路を設けたことを特徴とする請求項
1記載のレーダ信号処理装置。
3. The radar signal processing according to claim 1, wherein said phase correction unit is provided with a power calculation circuit for calculating a power of each amplitude value with respect to the frequency and amplitude waveforms interpolated by said interpolation circuit. apparatus.
【請求項4】 上記位相補正部に上記補間回路で補間さ
れた周波数と振幅の波形について各振幅値のべき乗を計
算するべき乗計算回路と、上記べき乗計算回路で計算さ
れた周波数と振幅のべき乗の波形についてデータを切り
出すデータ切り出し回路と、上記データ切り出し回路で
データを切り出す際に必要な閾値を設定する閾値設定回
路を設けたことを特徴とする請求項1記載のレーダ信号
処理装置。
4. A power calculation circuit for calculating a power of each amplitude value for the frequency and amplitude waveforms interpolated by the interpolation circuit in the phase correction section, and a power of the frequency and amplitude calculated by the power calculation circuit. 2. The radar signal processing apparatus according to claim 1, further comprising a data extracting circuit for extracting data on the waveform, and a threshold setting circuit for setting a threshold necessary for extracting data by the data extracting circuit.
【請求項5】 上記位相補正部に上記補間回路で補間さ
れた周波数と振幅の波形についてデータを切り出すデー
タ切り出し回路と、上記データ切り出し回路でデータを
切り出す際に必要な閾値を設定する閾値設定回路と、上
記データ切り出し回路で切り出された周波数と振幅の波
形について各振幅値のべき乗を計算するべき乗計算回路
を設けたことを特徴とする請求項1記載のレーダ信号処
理装置。
5. A data extraction circuit for extracting data on a frequency and amplitude waveform interpolated by the interpolation circuit in the phase correction section, and a threshold setting circuit for setting a threshold necessary for extracting data by the data extraction circuit 2. The radar signal processing apparatus according to claim 1, further comprising a power calculation circuit for calculating a power of each amplitude value for the waveform of the frequency and the amplitude extracted by the data extraction circuit.
JP06979694A 1994-04-07 1994-04-07 Radar signal processing equipment Expired - Fee Related JP3223697B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP06979694A JP3223697B2 (en) 1994-04-07 1994-04-07 Radar signal processing equipment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP06979694A JP3223697B2 (en) 1994-04-07 1994-04-07 Radar signal processing equipment

Publications (2)

Publication Number Publication Date
JPH07280923A JPH07280923A (en) 1995-10-27
JP3223697B2 true JP3223697B2 (en) 2001-10-29

Family

ID=13413073

Family Applications (1)

Application Number Title Priority Date Filing Date
JP06979694A Expired - Fee Related JP3223697B2 (en) 1994-04-07 1994-04-07 Radar signal processing equipment

Country Status (1)

Country Link
JP (1) JP3223697B2 (en)

Also Published As

Publication number Publication date
JPH07280923A (en) 1995-10-27

Similar Documents

Publication Publication Date Title
JP3221291B2 (en) Image processing device, image processing method, noise elimination device, and noise elimination method
US7697067B2 (en) Digital video processing systems and methods for estimating horizontal sync in digital video signals
US7720266B2 (en) Ultrasound image enhancement and speckle mitigation method
US7420623B2 (en) Image-signal processing apparatus for enhancing specific frequency component of image signal and correcting the edge of reproduced image
EP0544546A2 (en) Multiple discrete autofocus
EP2537331B1 (en) Imaging device including target tracking function
US6278494B1 (en) Edge emphasis device, image forming apparatus using the same, image signal processing method, and image forming method using the same
JP3223697B2 (en) Radar signal processing equipment
JP2822827B2 (en) Radar signal processing equipment
JP2000342585A (en) Ultrasonic diagnostic device
JP3788117B2 (en) Radar signal processing device
US20110181777A1 (en) Automatic quantization clock phase adjustable display apparatus
JP3788079B2 (en) Radar signal processing device
JP3747856B2 (en) Radar signal processing device
JP3788157B2 (en) Radar signal processing device
JP5403450B1 (en) Image processing apparatus and image processing method
JP2002116251A (en) Radar signal processing device
JP2004061332A (en) Radar signal processor
JP2909535B1 (en) Pulse compression method and apparatus
JP2008259097A (en) Video signal processing circuit and video display device
JP4995460B2 (en) Digital video processing system and method for setting horizontal synchronization of digital video signal
JP3791442B2 (en) Radar signal processing device
JP3709701B2 (en) Radar signal processor
JP4063069B2 (en) Radar signal processing device
JP3770176B2 (en) Radar signal processing device

Legal Events

Date Code Title Description
FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20070824

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20080824

Year of fee payment: 7

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20080824

Year of fee payment: 7

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20090824

Year of fee payment: 8

LAPS Cancellation because of no payment of annual fees