JP5865794B2 - Radar equipment - Google Patents
Radar equipment Download PDFInfo
- Publication number
- JP5865794B2 JP5865794B2 JP2012157392A JP2012157392A JP5865794B2 JP 5865794 B2 JP5865794 B2 JP 5865794B2 JP 2012157392 A JP2012157392 A JP 2012157392A JP 2012157392 A JP2012157392 A JP 2012157392A JP 5865794 B2 JP5865794 B2 JP 5865794B2
- Authority
- JP
- Japan
- Prior art keywords
- block
- pulse
- range
- fourier transform
- frequency spectrum
- 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
Links
Images
Description
この発明は、目標のレーダ画像を観測するに際して、観測時間を長くすることでドップラー周波数の分解能と信号対雑音電力比を高めるレーダ装置に関するものである。 The present invention relates to a radar apparatus that increases the resolution of the Doppler frequency and the signal-to-noise power ratio by extending the observation time when observing a target radar image.
レーダ装置として、既存の電波発信源から発信される電波を利用して、目標のレーダ画像を観測するパッシブレーダがある。
パッシブレーダは、自ら電波を放射しないため、省電力・省電波資源に資する方式として注目されている。電波発信源としては、テレビやラジオなどの電波発信源に加えて、GNSS(Global Navigation Satellite System) などが検討されている。
パッシブレーダでは、図1に示すように、電波発信源から送信された電波のうち、その電波発信源から直接受信局に到来する直接波と、目標に散乱された後に受信局に到来する散乱波とを受信局で受信することで、その直接波の経路と散乱波の経路との差(経路長差)や、散乱波のドップラー周波数シフトを計測する方式を採用している。
As a radar apparatus, there is a passive radar that observes a target radar image using radio waves transmitted from an existing radio wave transmission source.
Since passive radar does not emit radio waves by itself, it is attracting attention as a method that contributes to power saving and radio resource saving. As radio wave sources, in addition to radio wave sources such as television and radio, GNSS (Global Navigation Satellite System) and the like are being studied.
In the passive radar, as shown in FIG. 1, among the radio waves transmitted from the radio wave transmission source, the direct wave that directly arrives at the receiving station from the radio wave source and the scattered wave that arrives at the receiving station after being scattered by the target. Is received by the receiving station, and a method of measuring the difference between the direct wave path and the scattered wave path (path length difference) and the Doppler frequency shift of the scattered wave is adopted.
例えば、以下の非特許文献1には、パッシブレーダに関する従来の開発成果や、パッシブレーダの利点・欠点などが体系的に記載されている。
パッシブレーダの最大の課題は、探知距離の延伸である。
探知距離を延伸させるために、信号の積分時間を延長して、SNR(Signal to Noise Ratio)を改善する方式が、例えば、以下の特許文献1〜3に開示されている。
信号の積分時間を延長すると、目標が積分時間内にレンジセルを移動してしまう問題が発生するが、以下の特許文献1〜3には、この問題に対処する方法が開示されている。
For example, Non-Patent
The biggest problem with passive radar is extending the detection distance.
In order to extend the detection distance, a method of extending the signal integration time to improve SNR (Signal to Noise Ratio) is disclosed in, for example,
If the integration time of the signal is extended, there is a problem that the target moves the range cell within the integration time. However,
以下の特許文献1〜3に開示されている対処方法は、初めに、比較的短い積分時間でドップラー処理を実施し、その処理結果を用いて目標候補を検出する。
次に、目標候補の検出信号のドップラー周波数をもって、目標信号の補償処理を実施する方式である。
したがって、この対処方法では、比較的短い積分時間で、ある程度、目標が検出されていることが前提となる。
In the coping methods disclosed in
Next, a target signal compensation process is performed with the Doppler frequency of the target candidate detection signal.
Therefore, in this coping method, it is assumed that the target is detected to some extent in a relatively short integration time.
また、以下の非特許文献2,3には、電波発信源から直接受信局に到来する直接波の信号に対して、“Stretch Processing”という処理を施すことにより、直接波の信号から、目標の速度に応じたレンジセルの移動を予め織り込んだ参照信号を生成し、この参照信号と散乱波の信号との相互相関を求めることによって、レンジセルの移動を補償し、積分可能な時間を延長する方式が開示されている。
Further, in the following
以下の非特許文献4には、ドップラー周波数の変化による加速度的な影響を考慮して、目標の加速度に応じたレンジセルの移動を予め織り込んだ参照信号を生成し、この参照信号と散乱波の信号との相互相関を求めることによって、レンジセルの移動を補償し、非特許文献2,3に開示されている方式よりも、さらに積分可能な時間を延長する方式が開示されている。
しかし、非特許文献2〜4に開示されている方式では、演算量が多くなってしまう問題がある。
In
However, the methods disclosed in
パッシブレーダの最大の課題は、上述したように、探知距離の延伸であるが、パッシブレーダにおいて、探知距離の延伸が困難である原因は以下の2つである。
(1)警戒管制レーダなどと比較して、電波発信源から発信される電波(放送波)の送信電力が微弱であるため、SNR(Signal to Noise Ratio)が極めて低い。
(2)散乱波の観測チャネルにおいては、直接波が干渉波として振舞うが、直接波の信号レベルが散乱波に対してはるかに大きいため、SIR(Signal to Interference Ratio)が極めて低い。
したがって、探知距離を延伸させるには、SNRの向上とSIRの向上が課題となる。
As described above, the biggest problem of the passive radar is the extension of the detection distance. However, in the passive radar, the extension of the detection distance is difficult for the following two reasons.
(1) Since the transmission power of the radio wave (broadcast wave) transmitted from the radio wave transmission source is weak compared to a warning control radar or the like, the SNR (Signal to Noise Ratio) is extremely low.
(2) In the observation channel of the scattered wave, the direct wave behaves as an interference wave, but since the signal level of the direct wave is much higher than that of the scattered wave, the SIR (Signal to Interference Ratio) is extremely low.
Therefore, in order to extend the detection distance, improvement of SNR and improvement of SIR are problems.
従来のレーダ装置は以上のように構成されているので、信号の積分時間を延長すれば、探知距離を延伸させることができるが、信号の積分時間を延長すると、目標が積分時間内にレンジセルを移動してしまう問題が発生する。
特許文献1〜3では、この問題に対処する方法を開示しているが、初めに、比較的短い積分時間で、ある程度、目標を検出できていることが前提であるため、目標を検出できていなければ、この問題に対処することができない課題があった。
また、非特許文献2〜4にも、積分可能な時間を延長する方式が開示されているが、演算量が多くなってしまう課題があった。
Since the conventional radar apparatus is configured as described above, the detection distance can be extended by extending the signal integration time. However, if the signal integration time is extended, the target can set the range cell within the integration time. The problem of moving will occur.
この発明は上記のような課題を解決するためになされたもので、積分時間内の目標のレンジセルの移動や演算量の増大を抑えつつ、目標の探知距離を延伸させることができるレーダ装置を得ることを目的とする。 The present invention has been made to solve the above-described problems, and provides a radar apparatus capable of extending a target detection distance while suppressing movement of a target range cell within an integration time and an increase in calculation amount. For the purpose.
この発明に係るレーダ装置は、電波発信源から送信された電波の直接波を受信するとともに、電波発信源から送信された電波が目標によって散乱された散乱波を受信する電波受信手段と、電波受信手段により受信された直接波の信号及び散乱波の信号をパルスに分割し、各々のパルスについて、直接波の信号と散乱波の信号との相互相関を求めることで、パルス毎のレンジプロフィールを生成するパルス毎レンジ圧縮手段と、パルス毎レンジ圧縮手段により生成されたパルス毎のレンジプロフィールをブロック単位にまとめ、複数のパルスからなるブロック単位のレンジプロフィールをヒット方向にフーリエ変換することで、ブロック毎のドップラー周波数スペクトルを算出するブロック毎ドップラー処理手段と、ブロック毎ドップラー処理手段により算出されたドップラー周波数スペクトルの中で、ブロック方向に発生している1次のレンジマイグレーションを補償する1次レンジマイグレーション補償手段と、1次レンジマイグレーション補償手段により1次のレンジマイグレーションが補償されたドップラー周波数スペクトルの中で、ブロック方向に発生している2次以上のレンジマイグレーションを補償する高次レンジマイグレーション補償手段と、高次レンジマイグレーション補償手段により2次以上のレンジマイグレーションが補償されたドップラー周波数スペクトルをブロック方向にフーリエ変換することで、ドップラー周波数スペクトルを算出するブロック方向ドップラー処理手段とを設け、目標検出手段が、ブロック方向ドップラー処理手段により算出されたドップラー周波数スペクトルから目標を検出するようにしたものである。
また、1次レンジマイグレーション補償手段が、ブロック毎ドップラー処理手段により算出されたドップラー周波数スペクトルをレンジ方向にフーリエ変換するパルス毎フーリエ変換手段と、ブロック毎ドップラー処理手段により算出されたドップラー周波数スペクトルにおける各々のドップラー周波数セルに対応するブロック方向の1次のレンジマイグレーションの量に相当する位相変化を補償するための1次レンジマイグレーション位相補償関数をパルス毎フーリエ変換手段のフーリエ変換結果に乗算する位相補償関数乗算手段と、位相補償関数乗算手段の乗算結果をレンジ方向に逆フーリエ変換するパルス毎逆フーリエ変換手段とから構成されているようにしたものである。
A radar apparatus according to the present invention includes a radio wave receiving means for receiving a direct wave of a radio wave transmitted from a radio wave source and receiving a scattered wave obtained by scattering the radio wave transmitted from the radio wave source from a target; Generates a range profile for each pulse by dividing the direct wave signal and scattered wave signal received by the means into pulses and obtaining the cross-correlation between the direct wave signal and the scattered wave signal for each pulse. The pulse-by-pulse range compressing means and the range profile for each pulse generated by the pulse-by-pulse range compressing means are grouped into blocks, and the block-by-block range profile consisting of a plurality of pulses is Fourier-transformed in the hit direction. Block Doppler processing means for calculating the Doppler frequency spectrum of the block, and block Doppler processing In the Doppler frequency spectrum calculated by the stage, the primary range migration compensation means for compensating the primary range migration occurring in the block direction and the primary range migration compensation means compensate the primary range migration. In the Doppler frequency spectrum, higher-order range migration compensation means for compensating second-order or higher range migration generated in the block direction, and second-order or higher range migration compensated by the higher-order range migration compensation means. Block direction Doppler processing means for calculating the Doppler frequency spectrum by Fourier transforming the frequency spectrum in the block direction, and the target detecting means is configured to calculate the Doppler processing means by the block direction Doppler processing means. It is obtained by the puller frequency spectrum to detect a target.
Further, the primary range migration compensation unit includes a Fourier transform unit for each pulse Fourier transforming the Doppler frequency spectrum calculated by the Doppler processing unit for each block in the range direction, and a Doppler frequency spectrum calculated by the Doppler processing unit for each block. Compensation function for multiplying the Fourier transform result of the pulse-by-pulse Fourier transform means by the primary range migration phase compensation function for compensating for the phase change corresponding to the amount of the first-order range migration in the block direction corresponding to the Doppler frequency cell. The multiplication means and a pulse-by-pulse inverse Fourier transform means for performing inverse Fourier transform on the multiplication result of the phase compensation function multiplying means in the range direction.
この発明によれば、高次レンジマイグレーション補償手段が、1次レンジマイグレーション補償手段により1次のレンジマイグレーションが補償されたドップラー周波数スペクトルの中で、ブロック方向に発生している2次以上のレンジマイグレーションを補償するように構成したので、積分時間内の目標のレンジセルの移動や演算量の増大を抑えつつ、目標の探知距離を延伸させることができる効果がある。 According to the present invention, the higher-order range migration compensating means has a second or higher order range migration occurring in the block direction in the Doppler frequency spectrum in which the first-order range migration is compensated by the first-order range migration compensating means. Therefore, there is an effect that the target detection distance can be extended while suppressing the movement of the target range cell within the integration time and the increase in the calculation amount.
実施の形態1.
この実施の形態1では、パッシブレーダ方式を採用しているレーダ装置について説明するが、信号を得るまでの過程において、電波発信源がレーダ装置の制御下にあれば、直接波受信用のアンテナを省略することが可能であり、アクティブレーダとして扱えることが可能である。
図1はパッシブレーダ方式のレーダ装置を示す概念図である。
図1において、電波発信源1は搬送波周波数がfc、信号帯域がBの信号(電波)を継続的に発信(放送)している発信源である。
In the first embodiment, a radar apparatus adopting a passive radar system will be described. However, if a radio wave source is under the control of the radar apparatus in the process until a signal is obtained, an antenna for direct wave reception is used. It can be omitted and can be treated as an active radar.
FIG. 1 is a conceptual diagram showing a passive radar type radar apparatus.
In FIG. 1, a radio
直接波受信用アンテナ2は電波発信源1を指向するように配置されており、電波発信源1から送信された電波の直接波を受信する。
散乱波受信用アンテナ3は、目標の存在する観測領域を指向するように配置されており、電波発信源1から送信された電波が目標によって散乱された散乱波を受信する。
受信局4は直接波受信用アンテナ2及び散乱波受信用アンテナ3の受信信号を増幅し、増幅後の受信信号を帯域フィルタに通して、所望の帯域の信号を取り出した後、所望の帯域の信号をダウンコンバートしてサンプリングする処理を実施する。
The direct
The scattered
The receiving
ただし、信号帯域が既知又は可変であり、かつ、サンプリング周波数が可変又は既知であれば、信号帯域やサンプリング周波数を選択することで、所望の帯域の信号をダウンコンバートせずにサンプリングすることも可能である。
また、電波発信源1の信号が既知であれば、直接波受信用アンテナ2を省略することが可能であることは言うまでもない。
この実施の形態1では、直接波受信用アンテナ2と散乱波受信用アンテナ3を別の実体として説明するが、2つ以上のアンテナで受信した信号を用いて、デジタルビームフォーミングによって直接波と散乱波を分離するように構成しても構わない。
However, if the signal band is known or variable and the sampling frequency is variable or known, the signal in the desired band can be sampled without down-conversion by selecting the signal band or sampling frequency. It is.
Needless to say, if the signal from the
In the first embodiment, the direct
また、散乱波受信用アンテナ3として、互いに直交する偏波特性を有する2つの受信アンテナを用いることによって、目標によって散乱された散乱波の偏波特性を計測することも可能である。
この場合、以下に説明する各処理を、2つの散乱波受信用アンテナで得られた信号に対して各々適用することによって、偏波特性の異なる2つのレーダ画像を生成することが可能である。
Further, by using two receiving antennas having polarization characteristics orthogonal to each other as the scattered
In this case, it is possible to generate two radar images having different polarization characteristics by applying each processing described below to signals obtained by the two scattered wave receiving antennas. .
図1において、psバー(明細書の文書中では、電子出願の関係上、文字の上部に“−”の記号を付することができないため、「psバー」のように表記している)は固定の電波発信源1の位置を表す位置ベクトルである。
prバーは受信局4の位置を表す位置ベクトルである。
ptバーは目標の重心位置を表す位置ベクトルであり、vバーは目標の速度を表す速度ベクトルである。
図1では、電波発信源1と受信局4が異なる位置に設置されている例を示しているが、電波発信源1と受信局4が同じ位置に設置されていてもよい。
In FIG. 1, p s bar (in the specification document, the symbol “−” cannot be attached to the upper part of the letter because of electronic application, so it is represented as “ ps bar”. ) Is a position vector representing the position of the fixed
The pr bar is a position vector representing the position of the receiving
The pt bar is a position vector that represents the center of gravity of the target, and the v bar is a speed vector that represents the target speed.
Although FIG. 1 shows an example in which the
isドットハット(明細書の文書中では、電子出願の関係上、文字の上部に“・”や“^”の記号を付することができないため、「isドットハット」のように表記している)は目標から電波発信源1への向きを表す単位ベクトルである。
irハットは目標から受信局4への向きを表す単位ベクトルである。
i s dot hat (In the document of the description, because of electronic application, the symbol “・” or “^” cannot be added to the upper part of the letter, so it is expressed as “i s dot hat”. Is a unit vector representing the direction from the target to the
The ir hat is a unit vector representing the direction from the target to the receiving
このとき、目標と電波発信源1の距離rs、目標と受信局4の距離rr、電波発信源1と受信局4の距離rdは、下記の式(2)で表すことができる。
ただし、目標は移動しているので、ptバー、vバー、isドットハット、irドットハット及びrrは時刻tの関数であり、以下の説明では必要に応じて明示的にptバー(t)などのように表記する。
At this time, the target
However, since the target is moving, p t bar, v bar, i s dot hat, i r dot hat and r r is a function of time t, explicitly p t, if necessary in the following description Notation such as bar (t) is used.
図2はこの発明の実施の形態1によるレーダ装置を示す構成図であり、図2において、図1と同一符号は同一または相当部分を示すので説明を省略する。
直接波受信機5は受信局4に搭載されており、直接波受信用アンテナ2の受信信号を増幅し、増幅後の受信信号を帯域フィルタに通して、所望の帯域の信号を取り出した後、所望の帯域の信号をダウンコンバートしてサンプリングする。
散乱波受信機6は受信局4に搭載されており、散乱波受信用アンテナ3の受信信号を増幅し、増幅後の受信信号を帯域フィルタに通して、所望の帯域の信号を取り出した後、所望の帯域の信号をダウンコンバートしてサンプリングする。
なお、直接波受信用アンテナ2、散乱波受信用アンテナ3及び受信局4から電波受信手段が構成されている。
2 is a block diagram showing a radar apparatus according to
The
The scattered
The direct
パルス毎レンジ圧縮部7は例えばCPUを実装している半導体集積回路、あるいは、ワンチップマイコンなどから構成されており、受信局4の直接波受信機5から出力された直接波信号(ダウンコンバートしてサンプリングされた信号)及び散乱波受信機6から出力された散乱波信号(ダウンコンバートしてサンプリングされた信号)を短時間のパルスに分割し、各々のパルスについて、その直接波信号と散乱波信号との相互相関を求めることで、パルス毎のレンジプロフィールを生成する処理を実施する。なお、パルス毎レンジ圧縮部7はパルス毎レンジ圧縮手段を構成している。
The pulse-by-pulse
クラッタ抑圧部8は例えばCPUを実装している半導体集積回路、あるいは、ワンチップマイコンなどから構成されており、パルス毎レンジ圧縮部7により生成されたパルス毎のレンジプロフィールから、各パルスのレンジプロフィールのパルス方向平均を差し引く処理を実施することで、背景の静止物からの反射信号(クラッタ)を抑圧する。
ブロック毎ドップラー処理部9は例えばCPUを実装している半導体集積回路、あるいは、ワンチップマイコンなどから構成されており、クラッタ抑圧部8によりクラッタが抑圧されたパルス毎のレンジプロフィールをブロック単位にまとめ、ブロック単位のレンジプロフィールをヒット方向にフーリエ変換することで、ブロック毎のドップラー周波数スペクトルを算出する処理を実施する。なお、ブロック毎ドップラー処理部9はブロック毎ドップラー処理手段を構成している。
The
The block-by-block
ドップラー周波数セル対応1次レンジマイグレーション・位相補償部10は例えばCPUを実装している半導体集積回路、あるいは、ワンチップマイコンなどから構成されており、ブロック毎ドップラー処理部9により算出されたドップラー周波数スペクトルの中で、ブロック方向に発生している1次のレンジマイグレーションをドップラー周波数セル毎に補償する処理を実施する。なお、ドップラー周波数セル対応1次レンジマイグレーション・位相補償部10は1次レンジマイグレーション補償手段を構成している。
The Doppler frequency cell compatible primary range migration /
高次レンジマイグレーション・位相補償部11は例えばCPUを実装している半導体集積回路、あるいは、ワンチップマイコンなどから構成されており、ドップラー周波数セル対応1次レンジマイグレーション・位相補償部10により1次のレンジマイグレーションが補償されたドップラー周波数スペクトルの中で、ブロック方向に発生している2次のレンジマイグレーションを補償する処理を実施する。なお、高次レンジマイグレーション・位相補償部11は高次レンジマイグレーション補償手段を構成している。
この実施の形態1では、説明の簡単化のため、高次レンジマイグレーション・位相補償部11がブロック方向に発生している2次のレンジマイグレーションを補償するものを示しているが、3次以上のレンジマイグレーションを補償するようにしてもよい。
The high-order range migration /
In the first embodiment, for simplification of explanation, the high-order range migration /
ブロック方向ドップラー処理部12は例えばCPUを実装している半導体集積回路、あるいは、ワンチップマイコンなどから構成されており、高次レンジマイグレーション・位相補償部11により2次のレンジマイグレーションが補償されたドップラー周波数スペクトルをブロック方向にフーリエ変換することで、ドップラー周波数スペクトルを算出する処理を実施する。なお、ブロック方向ドップラー処理部12はブロック方向ドップラー処理手段を構成している。
The block direction
目標検出部13は例えばCPUを実装している半導体集積回路、あるいは、ワンチップマイコンなどから構成されており、ブロック方向ドップラー処理部12により算出されたドップラー周波数スペクトルから目標を検出する処理を実施する。なお、目標検出部13は目標検出手段を構成している。
The
図2の例では、レーダ装置の構成要素である直接波受信用アンテナ2、散乱波受信用アンテナ3、受信局4、パルス毎レンジ圧縮部7、クラッタ抑圧部8、ブロック毎ドップラー処理部9、ドップラー周波数セル対応1次レンジマイグレーション・位相補償部10、高次レンジマイグレーション・位相補償部11、ブロック方向ドップラー処理部12及び目標検出部13のそれぞれが専用のハードウェアで構成されているものを想定しているが、レーダ装置の全部又は一部がコンピュータで構成されていてもよい。
例えば、レーダ装置の一部(直接波受信用アンテナ2、散乱波受信用アンテナ3及び受信局4を除く部分)がコンピュータで構成されている場合、パルス毎レンジ圧縮部7、クラッタ抑圧部8、ブロック毎ドップラー処理部9、ドップラー周波数セル対応1次レンジマイグレーション・位相補償部10、高次レンジマイグレーション・位相補償部11、ブロック方向ドップラー処理部12及び目標検出部13の処理内容を記述しているプログラムをコンピュータのメモリに格納し、当該コンピュータのCPUが当該メモリに格納されているプログラムを実行するようにすればよい。
In the example of FIG. 2, the direct
For example, when a part of the radar device (a part excluding the direct
図3はこの発明の実施の形態1によるレーダ装置のパルス毎レンジ圧縮部7を示す構成図である。
図3において、パルス毎FFT部21は受信局4の直接波受信機5から出力された直接波信号をレンジ方向にフーリエ変換することで、その直接波信号を短時間のパルスに分割する処理を実施する。
パルス毎FFT部22は受信局4の散乱波受信機6から出力された散乱波信号をレンジ方向にフーリエ変換することで、その散乱波信号を短時間のパルスに分割する処理を実施する。
FIG. 3 is a block diagram showing the pulse-by-pulse
In FIG. 3, the
The per-
複素共役乗算部23はパルス毎FFT部21のフーリエ変換結果とパルス毎FFT部22のフーリエ変換結果との複素共役乗算を行うことで、その直接波信号と散乱波信号との相互相関を求める処理を実施する。
パルス毎IFFT部24は複素共役乗算部23の乗算結果をレンジ方向に逆フーリエ変換する処理を実施する。
The
The pulse-by-
図4はこの発明の実施の形態1によるレーダ装置のドップラー周波数セル対応1次レンジマイグレーション・位相補償部10を示す構成図である。
図4において、パルス毎FFT部31はブロック毎ドップラー処理部9により算出されたドップラー周波数スペクトルをレンジ方向にフーリエ変換する処理を実施する。なお、パルス毎FFT部31はパルス毎フーリエ変換手段を構成している。
1次レンジマイグレーション・位相補償関数格納部32はブロック毎ドップラー処理部9により算出されたドップラー周波数スペクトルにおける各々のドップラー周波数セルに対応するブロック方向の1次のレンジマイグレーションの量に相当する位相変化を補償するための1次レンジマイグレーション位相補償関数を格納しているメモリである。
FIG. 4 is a block diagram showing the primary range migration /
In FIG. 4, the pulse-by-
The primary range migration / phase compensation
位相補償関数乗算部33は1次レンジマイグレーション・位相補償関数格納部32により格納されている1次レンジマイグレーション位相補償関数をパルス毎FFT部31のフーリエ変換結果に乗算する処理を実施する。なお、位相補償関数乗算部33は位相補償関数乗算手段を構成している。
パルス毎IFFT部34は位相補償関数乗算部33の乗算結果をレンジ方向に逆フーリエ変換する処理を実施する。なお、パルス毎IFFT部34はパルス毎逆フーリエ変換手段を構成している。
The phase compensation
The pulse-by-
図5はこの発明の実施の形態1によるレーダ装置の高次レンジマイグレーション・位相補償部11を示す構成図である。
図5において、ブロック毎IFFT部41はドップラー周波数セル対応1次レンジマイグレーション・位相補償部10により1次のレンジマイグレーションが補償されたドップラー周波数スペクトルをヒット方向に逆フーリエ変換する処理を実施する。なお、ブロック毎IFFT部41はブロック毎逆フーリエ変換手段を構成している。
FIG. 5 is a block diagram showing the high-order range migration /
In FIG. 5, the block-by-
高次項設定部42はレンジマイグレーションの2次の項が取り得る範囲を離散量で設定する処理を実施する。
ドップラーセルマイグレーション補償関数格納部43はレンジマイグレーションにおける2次の項に影響されるドップラーセルマイグレーションを補償するためのドップラーセルマイグレーション補償関数(高次項設定部42の設定内容である後述のαが反映されたドップラーセルマイグレーション補償関数)を格納しているメモリである。
高次レンジマイグレーション・位相補償関数格納部44は2次のレンジマイグレーションを補償するための高次レンジマイグレーション位相補償関数(高次項設定部42の設定内容である後述のαが反映された高次レンジマイグレーション位相補償関数)を格納しているメモリである。
The high-order
The Doppler cell migration compensation
The high-order range migration / phase compensation
ドップラーセルマイグレーション補償関数乗算部45はドップラーセルマイグレーション補償関数格納部43により格納されているドップラーセルマイグレーション補償関数をブロック毎IFFT部41の逆フーリエ変換結果に乗算する処理を実施する。なお、ドップラーセルマイグレーション補償関数乗算部45は第1の補償関数乗算手段を構成している。
位相補償関数乗算部46は高次レンジマイグレーション・位相補償関数格納部44により格納されている高次レンジマイグレーション位相補償関数をドップラーセルマイグレーション補償関数乗算部45の乗算結果に乗算する処理を実施する。なお、位相補償関数乗算部46は第2の補償関数乗算手段を構成している。
The Doppler cell migration compensation
The phase compensation
ブロック毎FFT部47は位相補償関数乗算部46の乗算結果をヒット方向にフーリエ変換する処理を実施する。なお、ブロック毎FFT部47はブロック毎フーリエ変換手段を構成している。
パルス毎IFFT部48はブロック毎FFT部47のフーリエ変換結果をレンジ方向に逆フーリエ変換する処理を実施する。なお、パルス毎IFFT部48はパルス毎逆フーリエ変換手段を構成している。
The block-by-
The pulse-by-
次に動作について説明する。
図6はこの発明の実施の形態1によるレーダ装置の処理における受信信号の取り扱いを示す説明図である。
図6の横軸は時刻であり、図6に示した観測時間T[sec]の間、直接波受信機5及び散乱波受信機6が信号を受信する。
この実施の形態1では、受信信号が時間幅Tb[sec]の ブロックに分割される(N個のブロックに分割される)。
また、各々のブロックが時間幅T0[sec]のパルスに分割される(M個のパルスに分割される)。
Next, the operation will be described.
FIG. 6 is an explanatory diagram showing how received signals are handled in the processing of the radar apparatus according to
The horizontal axis of FIG. 6 is time, and the
In the first embodiment, the received signal is divided into blocks of time width T b [sec] (divided into N blocks).
Each block is divided into pulses having a time width T 0 [sec] (divided into M pulses).
以降、時間幅Tbをブロック幅と称し、時間幅T0をパルス幅と称する。
なお、パルスとパルスは隣接しているため、パルスの繰り返し周期はパルス幅T0と一致する。そのため、必要に応じてT0[sec]をパルス繰返し周期(PRI:Pulse Repetition Interval)と呼ぶ場合がある。
定義より、T、Tb、T0は、下記の次の関係を満足する。
Hereinafter, the time width T b is referred to as a block width, and the time width T 0 is referred to as a pulse width.
Since the pulses are adjacent to each other, the pulse repetition period matches the pulse width T 0 . For this reason, T 0 [sec] may be referred to as a pulse repetition period (PRI) as necessary.
By definition, T, T b , and T 0 satisfy the following relationship.
この実施の形態1において、信号をパルスに分割したり、ブロックにまとめたりするのは、処理の高速化を図るためである。
詳細は後述するが、長時間観測したデータを短いパルスに区切ることにより、直接波と散乱波の長時間の相関処理を全てFFT(Fast Fourier Transform)によって構成することができる。
また、いくつかのパルスをまとめたブロック単位の処理を導入することで、観測時間中の1次レンジマイグレーションへの対処を高速化することができる。
In the first embodiment, the signal is divided into pulses or grouped into blocks in order to increase the processing speed.
Although details will be described later, all long-time correlation processing of the direct wave and the scattered wave can be configured by FFT (Fast Fourier Transform) by dividing the data observed for a long time into short pulses.
In addition, by introducing processing in units of blocks in which several pulses are combined, it is possible to speed up the response to primary range migration during the observation time.
続いて、パルス毎レンジ圧縮部7及びその後の処理ブロックの動作を説明するために、電波発信源1から発信された直接波の信号と、目標によって散乱された散乱波の信号とを定式化する。
電波発信源1から発信される信号は、信号帯域幅がB、中心周波数がfc、ベースバンド信号がw(t)の狭帯域信号であるから、直接波の信号をsd(t)、散乱波の信号をss(t)とすると、下記のように表すことができる。
以下では、広義定常性(WSS:wide−sense stationary)が成立するものとする。
ただし、adは直接波受信時の信号の振幅、asは散乱波受信時の信号の振幅、τdは直接波が受信局4に到達するまでの遅延時間、R0(t)は散乱波の遅延距離、cは伝播線路中の電磁波の速さ、tは時間変数、fdはドップラー周波数である。
Subsequently, in order to explain the operation of the pulse-by-pulse
Since the signal transmitted from the
In the following, it is assumed that wide-sense stationarity (WSS: wide-sense stationery) is established.
However, a d is the amplitude of the signal at the direct wave reception, a s the delay time until the amplitude of the signal at the time of the scattered wave reception, tau d is direct wave reaches the receiving
ここで、直接波の信号sd(t)及び散乱波の信号ss(t)の位相検波を実施することにより、次のビデオ信号sd(t)チルダ,ss(t)チルダが得られる。
明細書の文書中では、電子出願の関係上、文字の上部に“〜”の記号を付することができないため、「sd(t)チルダ」のように表記している。
Here, the next video signal s d (t) tilde and s s (t) tilde are obtained by performing phase detection of the direct wave signal s d (t) and the scattered wave signal s s (t). It is done.
In the document of the specification, the symbol “˜” cannot be added to the upper part of the character because of the electronic application, and therefore, “s d (t) tilde” is used.
直接波信号と散乱波信号との遅延距離差R(t)は、下記のように定義する。
この実施の形態1では、目標のドップラー周波数fd(t)が1次の時間変化をする場合を想定する。
そこで、ドップラー周波数の変化率α[Hz/s]を導入し(αは、高次レンジマイグレーション・位相補償部11の高次項設定部42により設定されるレンジマイグレーションの2次の項が取り得る範囲に相当する)、目標のドップラー周波数fd(t)を次式でモデル化する。
The delay distance difference R (t) between the direct wave signal and the scattered wave signal is defined as follows.
In the first embodiment, it is assumed that the target Doppler frequency f d (t) undergoes a first-order time change.
Therefore, a change rate α [Hz / s] of the Doppler frequency is introduced (α is a range that can be taken by the secondary term of the range migration set by the high-order
よって、
であることを考慮すると、式(9)及び式(10)より、遅延距離差R(t)とラジアル速度v(t)は次式で表される。
Therefore,
Considering that, the delay distance difference R (t) and the radial velocity v (t) are expressed by the following equations from the equations (9) and (10).
なお、アルゴリズムの導出の過程で、ブロック幅Tb[sec]の間は、ドップラー周波数が一定であると見なせると仮定する。
あるブロックにおけるm番目のパルスの中心時刻をtmとし、n番目のブロックの中心時刻をtnとすると、tm,tnは次式によって定義される。
R(t)はtmの定義域において、次の1次関数Rn(tm)とする。
In the process of deriving the algorithm, it is assumed that the Doppler frequency can be assumed to be constant during the block width T b [sec].
The central time of the m-th pulse in a given block and t m, when the center time of the n th block to t n, t m, t n is defined by the following equation.
R (t) in the domain of t m, the following linear function R n (t m).
以下、レーダ装置の各ブロックの処理内容を具体的に説明する。
パルス毎レンジ圧縮部7は、受信局4の直接波受信機5から直接波信号を受け、散乱波受信機6から散乱波信号を受けると、その直接波信号及び散乱波信号を短時間のパルスに分割し、パルス毎に、その直接波信号と散乱波信号との相互相関を求めることで、パルス毎のレンジプロフィールを生成する。
即ち、パルス毎レンジ圧縮部7のパルス毎FFT部21は、受信局4の直接波受信機5から出力された直接波信号をレンジ方向にフーリエ変換することで、その直接波信号を短時間のパルスに分割する。
また、パルス毎レンジ圧縮部7のパルス毎FFT部22は、受信局4の散乱波受信機6から出力された散乱波信号をレンジ方向にフーリエ変換することで、その散乱波信号を短時間のパルスに分割する。
Hereinafter, the processing contents of each block of the radar apparatus will be specifically described.
When receiving the direct wave signal from the
That is, the pulse-by-
The pulse-by-
パルス毎レンジ圧縮部7の複素共役乗算部23は、パルス毎FFT部21のフーリエ変換結果とパルス毎FFT部22のフーリエ変換結果との複素共役乗算を行うことで、その直接波信号と散乱波信号との相互相関を求める。
ここで、複素共役乗算部23により求められる相互相関を示す相互相関関数x(τ,tm,tn)は、下記のように表され、レンジプロフィールと呼ばれる。
ただし、x(τ,tm,tn)はn番目のブロックのm番目のパルスにおけるレンジプロフィールである。
The
Here, the cross-correlation function x (τ, t m , t n ) indicating the cross-correlation obtained by the
Here, x (τ, t m , t n ) is a range profile in the mth pulse of the nth block.
ここで、ベースバンド信号w(t)については、広義定常性(WSS)が成立すると仮定しているので、自己相関関数の期待値C(τ)は次式で表される。
期待値について議論すると、次のようになる。
ただし、τはレンジ圧縮処理後の時間変数、*は複素共役、E[・]は期待値、C(・)は自己相関関数である。
Here, for the baseband signal w (t), since it is assumed that wide-sense stationarity (WSS) is established, the expected value C (τ) of the autocorrelation function is expressed by the following equation.
The expected value is discussed as follows.
Where τ is a time variable after range compression processing, * is a complex conjugate, E [•] is an expected value, and C (•) is an autocorrelation function.
周波数空間上でのレンジ圧縮は、次式で表現される。
ただし、Cf(f)は自己相関関数C(τ)のフーリエ変換である。F[g(τ)]τは関数g(τ)の変数τについてのフーリエ変換を表している。
パルス毎IFFT部24は、複素共役乗算部23の乗算結果をレンジ方向に逆フーリエ変換し、その逆フーリエ変換結果をパルス毎のレンジプロフィールとして、クラッタ抑圧部8に出力する。
Range compression in the frequency space is expressed by the following equation.
However, C f (f) is a Fourier transform of the autocorrelation function C (τ). F [g (τ)] τ represents the Fourier transform for the variable τ of the function g (τ).
The pulse IFFT
クラッタ抑圧部8は、パルス毎レンジ圧縮部7からパルス毎のレンジプロフィールを受けると、パルス毎のレンジプロフィールから、各パルスのレンジプロフィールのパルス方向平均を差し引く処理を実施することで、背景の静止物からの反射信号(クラッタ)を抑圧する。
即ち、クラッタ抑圧部8は、下記の式(20)に示す処理を実施することで、クラッタを抑圧する。
When the
That is, the
これにより、式(20)に示すクラッタ抑圧後の信号X逆弧(f,tm,tn)(明細書の文書中では、電子出願の関係上、文字の上部に“∪”の記号を付することができないため、「X逆弧」のように表記している)を出力するが、静止物からの反射波の信号については、遅延時間差がパルス毎に変化しない。即ち、ドップラー周波数シフトを受けていない。
このことは、式(19)において、目標が固定であれば、Rn(tm)=constとなり、パルス番号mによらず信号の位相が一定になることからも確認できる。
式(20)において、Nパルスのレンジプロフィールを平均することによって、ドップラー周波数がゼロの信号を抽出し、これを各レンジプロフィールから差し引くことで、ドップラー周波数がゼロの信号を抑圧している。これにより、背景の静止物からの反射信号を抑圧することができる。
なお、このクラッタ抑圧部8は、この実施の形態1を構成する必須の機能ではない。仮にクラッタ抑圧部8の処理を省いても、その他については、同様の処理を実施することが可能である。
As a result, the signal X reverse arc (f, t m , t n ) after clutter suppression shown in equation (20) (in the document of the specification, the symbol “∪” is placed above the characters for the purpose of electronic application. However, the delay time difference does not change for each pulse of a reflected wave signal from a stationary object. That is, it has not undergone a Doppler frequency shift.
This can be confirmed from the fact that, in the equation (19), if the target is fixed, R n (t m ) = const, and the phase of the signal is constant regardless of the pulse number m.
In Equation (20), a signal with a zero Doppler frequency is extracted by averaging the range profile of N pulses, and the signal with a zero Doppler frequency is suppressed by subtracting it from each range profile. Thereby, the reflected signal from the stationary object in the background can be suppressed.
The
ブロック毎ドップラー処理部9は、クラッタ抑圧部8がクラッタの抑圧処理を行うと、クラッタが抑圧されたパルス毎のレンジプロフィールをブロック単位にまとめ、ブロック単位のレンジプロフィールをヒット方向にフーリエ変換することで、ブロック毎のドップラー周波数スペクトルを算出する。
以下、ブロック毎ドップラー処理部9の処理内容を具体的に説明する。
まず、n番目のブロックにおけるレンジ圧縮後の信号X(f,tm,tn)は、式(15)を式(19)に代入することで、下式のように表すことができる。
ここでは、レンジ圧縮後の信号X(f,tm,tn)を式(21)で表しているが、式(15)を式(20)に代入することで、クラッタが抑圧されたレンジ圧縮後の信号X逆弧(f,tm,tn)を表すことができる。
When the
Hereinafter, the processing contents of the block-by-block
First, the range-compressed signal X (f, t m , t n ) in the nth block can be expressed by the following equation by substituting equation (15) into equation (19).
Here, the signal X (f, t m , t n ) after the range compression is represented by Expression (21), but the range in which clutter is suppressed by substituting Expression (15) into Expression (20). The compressed signal X reverse arc (f, t m , t n ) can be represented.
レンジ圧縮後の信号X(f,tm,tn)を、mについて離散フーリエ変換すると、次式が得られる。
ただし、fdkはk番目のドップラー周波数セルの中心周波数であり、次式で定義される。
When the range-compressed signal X (f, t m , t n ) is discrete Fourier transformed with respect to m, the following equation is obtained.
Here, f dk is the center frequency of the kth Doppler frequency cell and is defined by the following equation.
式(22)における近似は、次の関係式による。
fc >> f∈[−B/2,B/2] (24)
式(24)は、信号の比帯域が小さければ、式(22)の近似が成り立つことを示している。
ここでは、ブロック幅Tbにおける目標の移動量がレンジ分解能以下であることを示す下記の式(25)の関係を満たしている必要がある。
ただし、fdmaxはPRIで決まる値であり、fdmax=1/(2T0)である。したがって、式(25)は次のように整理することができる。
The approximation in equation (22) is based on the following relational expression.
f c >> f∈ [−B / 2, B / 2] (24)
Equation (24) indicates that the approximation of Equation (22) holds if the signal bandwidth is small.
Here, it is necessary to satisfy the relationship of formula (25) below which show that the amount of movement of the target in the block width T b is less range resolution.
However, f dmax is a value determined by PRI, and f dmax = 1 / (2T 0 ). Therefore, equation (25) can be organized as follows.
ドップラー周波数セル対応1次レンジマイグレーション・位相補償部10は、ブロック毎ドップラー処理部9がブロック毎のドップラー周波数スペクトルを算出すると、そのドップラー周波数スペクトルの中で、ブロック方向に発生している1次のレンジマイグレーションをドップラー周波数セル毎に補償する。
以下、ドップラー周波数セル対応1次レンジマイグレーション・位相補償部10の処理内容を具体的に説明する。
When the Doppler frequency cell compatible primary range migration /
The processing contents of the primary range migration /
ブロック毎ドップラー処理部9では、ブロック幅Tbの間は、ドップラー周波数が一定であることを仮定して処理を実施しているが、観測時間Tの全体を考える場合は、ドップラー周波数fd(t)が1次の変化をする式(9)のモデルを考える。
このことを踏まえ、式(22)のR(tn),v(tn)に式(11),式(12)を代入することで次式を得る。
The block-by-block
Based on this, the following equation is obtained by substituting the equations (11) and (12) into R (t n ) and v (t n ) of the equation (22).
この実施の形態1では、ブロック方向の処理を実施するために、目標信号のレンジ方向の移動(マイグレーション)を補償する。
ここでは、ドップラー周波数が1次の変化をするモデルを考えているので、目標信号のレンジは1次と2次の変化をする。ドップラー周波数セル対応1次レンジマイグレーション・位相補償部10では、1次の変化量を補償する。
k番目のドップラー周波数セルに含まれる信号のドップラー周波数はfdkであるから、n番目のブロックにおける遅延量の1次成分ΔRn(tn)は次式で表される。
In the first embodiment, in order to perform the processing in the block direction, the movement (migration) of the target signal in the range direction is compensated.
Here, since a model in which the Doppler frequency changes linearly is considered, the range of the target signal changes between primary and secondary. The primary range migration /
Since the Doppler frequency of the signal included in the k-th Doppler frequency cell is f dk , the primary component ΔR n (t n ) of the delay amount in the n-th block is expressed by the following equation.
式(27)に対する1次レンジマイグレーション位相補償関数Ψcmp1(f,fdm,tn)は次式で表され、予め、1次レンジマイグレーション位相補償関数Ψcmp1(f,fdm,tn)が1次レンジマイグレーション・位相補償関数格納部32に格納される。
ただし、*は複素共役である。
Equation (27) primary range for the migration phase compensation function Ψ cmp1 (f, f dm, t n) is expressed by the following equation, in advance, the primary-range migration phase compensation function Ψ cmp1 (f, f dm, t n) Is stored in the primary range migration / phase
However, * is a complex conjugate.
ドップラー周波数セル対応1次レンジマイグレーション・位相補償部10のパルス毎FFT部31は、ブロック毎ドップラー処理部9がブロック毎のドップラー周波数スペクトルを算出すると、そのドップラー周波数スペクトルをレンジ方向にフーリエ変換する。
ドップラー周波数セル対応1次レンジマイグレーション・位相補償部10の位相補償関数乗算部33は、パルス毎FFT部31がブロック毎のドップラー周波数スペクトルを算出すると、1次レンジマイグレーション・位相補償関数格納部32により格納されている1次レンジマイグレーション位相補償関数Ψcmp1(f,fdm,tn)をブロック毎のドップラー周波数スペクトルに乗算することで、1次のレンジマイグレーションを補償する。
When the block-by-block
The phase compensation
1次レンジマイグレーション位相補償関数Ψcmp1(f,fdm,tn)を用いた関数y(τ,fdk,tn)の1次のレンジマイグレーション補償は次のようになる。これは周波数空間における畳み込み演算に相当する。
ただし、δfdk=β−fdkである。
ドップラー周波数セル対応1次レンジマイグレーション・位相補償部10のパルス毎IFFT部34は、位相補償関数乗算部33の乗算結果をレンジ方向に逆フーリエ変換する。
The primary range migration compensation of the function y (τ, f dk , t n ) using the primary range migration phase compensation function ψ cmp1 (f, f dm , t n ) is as follows. This corresponds to a convolution operation in the frequency space.
However, δf dk = β−f dk .
The pulse IFFT
高次レンジマイグレーション・位相補償部11は、ドップラー周波数セル対応1次レンジマイグレーション・位相補償部10が1次のレンジマイグレーションを補償すると、1次のレンジマイグレーションが補償されたドップラー周波数スペクトルの中で、ブロック方向に発生している2次のレンジマイグレーションを補償する。
以下、高次レンジマイグレーション・位相補償部11の処理内容を具体的に説明する。
When the primary range migration /
The processing contents of the high-order range migration /
まず、高次レンジマイグレーション・位相補償部11のブロック毎IFFT部41は、ドップラー周波数セル対応1次レンジマイグレーション・位相補償部10により1次のレンジマイグレーションが補償されたドップラー周波数スペクトルをヒット方向に逆フーリエ変換する。
即ち、式(30)のドップラー周波数fdkについて、下式に示すように、逆離散フーリエ変換を実施する。
First, the block-by-
That is, the inverse discrete Fourier transform is performed on the Doppler frequency f dk of Expression (30) as shown in the following expression.
ここで、ドップラーセルマイグレーション補償関数格納部43には、レンジマイグレーションにおける2次の項に影響されるドップラーセルマイグレーションを補償するためのドップラーセルマイグレーション補償関数Ψcmp2(tm,tn,αハット)が格納されている。αハットは、探索的に代入する必要がある変数であるが、高次項設定部42により設定されたドップラー周波数の変化率α[Hz/s]と等しいものとする。
また、高次レンジマイグレーション・位相補償関数格納部44には、2次のレンジマイグレーションを補償するための高次レンジマイグレーション位相補償関数Ψcmp2(f,tn,αハット)が格納されている。高次レンジマイグレーション位相補償関数Ψcmp2(f,tn,αハット)内のαハットについても、探索的に代入する必要がある変数であるが、高次項設定部42により設定されたドップラー周波数の変化率α[Hz/s]と等しいものとする。
Here, the Doppler cell migration compensation
The high-order range migration / phase
ドップラーセルマイグレーション補償関数乗算部45は、ブロック毎IFFT部41が、1次のレンジマイグレーションが補償されたドップラー周波数スペクトルをヒット方向に逆フーリエ変換すると、ドップラーセルマイグレーション補償関数格納部43により格納されているドップラーセルマイグレーション補償関数Ψcmp2(tm,tn,αハット)をブロック毎IFFT部41の逆フーリエ変換結果に乗算することで、ドップラーセルマイグレーションを補償する。図7はドップラーセルマイグレーションの補償処理を示す説明図である。
また、位相補償関数乗算部46は、高次レンジマイグレーション・位相補償関数格納部44により格納されている高次レンジマイグレーション位相補償関数Ψcmp2(f,tn,αハット)をドップラーセルマイグレーション補償関数乗算部45の乗算結果に乗算する。
The Doppler cell migration compensation
Further, the phase compensation
即ち、式(31)に示す逆フーリエ変換結果に対して、式(32)に示すドップラーセルマイグレーション補償関数Ψcmp2(tm,tn,αハット)と、式(33)に示す高次レンジマイグレーション位相補償関数Ψcmp2(f,tn,αハット)とを乗算する処理を実施する。
αハットは、探索的に代入する必要があるが、以下、αハット=αが成立する場合について議論を進める。
That is, the Doppler cell migration compensation function Ψ cmp2 (t m , t n , α hat) shown in Expression (32) and the higher-order range shown in Expression (33) with respect to the inverse Fourier transform result shown in Expression (31). A process of multiplying the migration phase compensation function Ψ cmp2 (f, t n , α hat) is performed.
The α hat needs to be substituted in a search manner, but the discussion will proceed below for the case where α hat = α holds.
ブロック毎FFT部47は、位相補償関数乗算部46の乗算結果をヒット方向にフーリエ変換する。
パルス毎IFFT部48は、ブロック毎FFT部47のフーリエ変換結果をレンジ方向に逆フーリエ変換する。
式(34)に対して、fについて逆フーリエ変換、tm,tnについて離散フーリエ変換を実施すると、次のようになる。
The block-by-
The pulse-by-
When the inverse Fourier transform is performed on f and the discrete Fourier transform is performed on t m and t n with respect to the equation (34), the result is as follows.
ブロック方向ドップラー処理部12は、高次レンジマイグレーション・位相補償部11が2次のレンジマイグレーションを補償すると、2次のレンジマイグレーションが補償されたドップラー周波数スペクトルをブロック方向にフーリエ変換することで、ドップラー周波数スペクトルfdiを算出する。
以下、ブロック方向ドップラー処理部12の処理内容を具体的に説明する。
When the high-order range migration /
The processing contents of the block direction
式(35)の近似については、レンジマイグレーション補償処理後は、信号がブロック番号によらず、同一のレンジセルに存在していると見なせる。以降では、fdnを次式で表すようにする。ただし、fdkはk番目のドップラー周波数セルの中心周波数であり、次式で定義される。
また、新たなドップラー周波数fdiを次式のように定義する。
As for the approximation of Expression (35), after the range migration compensation process, it can be considered that the signal exists in the same range cell regardless of the block number. Hereinafter, f dn is expressed by the following equation. Here, f dk is the center frequency of the kth Doppler frequency cell and is defined by the following equation.
Also, a new Doppler frequency f di is defined as follows:
ただし、m(i),n(i)は次のように定義される。
floor(i)は、iの小数点以下を切り捨てる演算子であり、mod(i,N)はiをNで割った余りを得る演算子である。
式(37),(38)より、式(35)は次のようになる。
式(39)において、実質の周波数分解能Δfdは、1つ目のsinc関数よって、次式のようになる。
However, m (i) and n (i) are defined as follows.
floor (i) is an operator that truncates the decimal part of i, and mod (i, N) is an operator that obtains a remainder obtained by dividing i by N.
From the equations (37) and (38), the equation (35) is as follows.
In the equation (39), the actual frequency resolution Δf d is expressed by the following equation using the first sinc function.
目標検出部13は、ブロック方向ドップラー処理部12がドップラー周波数スペクトルfdiを算出すると、そのドップラー周波数スペクトルfdiを用いて、目標を検出する。
即ち、目標検出部13は、式(35)によって算出される信号xcmp2(τ,fdk,fdl)の強度P1(τ,fdk,fdl)、あるいは、式(39)によって算出される信号xcmp2(τ,fdi)の強度P2(τ,fdi)を算出する。
目標検出部13は、信号強度P1(τ,fdk,fdl)、あるいは、信号強度P2(τ,fdi)に対して、CFAR(Constant False Alarm Rate)処理などの検出処理を適用することによって、目標信号を検出する。CFAR処理については公知の技術であるため詳細な説明を省略する。
When the block direction
That is, the
The
なお、CFAR処理を実施する前に、信号強度P1(τ,fdk,fdl)、あるいは、信号強度P2(τ,fdi)について、レンジ方向又はドップラー周波数方向、あるいは、その両方向に互いに隣接する複数セルにまたがってインコヒーレントに積分を行うようにしてもよい。
この処理により、雑音成分の標準偏差を低減することができるほか、目標信号がドップラーレンジ方向又はドップラー周波数方向、あるいは、その両方向に広がっている場合、信号成分をインコヒーレントに積み上げることができるため、SNR(Signal to Noise Ratio)を改善することが可能になる。
Before performing the CFAR process, the signal intensity P 1 (τ, f dk , f dl ) or the signal intensity P 2 (τ, f di ) is set in the range direction, the Doppler frequency direction, or both directions. The integration may be performed incoherently across a plurality of adjacent cells.
This process can reduce the standard deviation of the noise component, and when the target signal spreads in the Doppler range direction or Doppler frequency direction, or both directions, the signal component can be accumulated incoherently, The SNR (Signal to Noise Ratio) can be improved.
以上の処理が、図1のレーダ装置の処理内容であるが、ドップラー周波数の分解能と、ドップラー周波数変化率のふり幅の全観測時間については、次の関係が成り立てばよい。
探索されるαは,次のようになる。
The above processing is the processing content of the radar apparatus of FIG. 1, but the following relationship may be established between the resolution of the Doppler frequency and the total observation time of the fluctuation width of the Doppler frequency change rate.
Α to be searched is as follows.
以上で明らかなように、この実施の形態1によれば、高次レンジマイグレーション・位相補償部11が、ドップラー周波数セル対応1次レンジマイグレーション・位相補償部10により1次のレンジマイグレーションが補償されたドップラー周波数スペクトルの中で、ブロック方向に発生している2次のレンジマイグレーションを補償するように構成したので、積分時間内の目標のレンジセルの移動や演算量の増大を抑えつつ、目標の探知距離を延伸させることができる効果を奏する。
即ち、1次のレンジマイグレーションを補償した後に、2次のレンジマイグレーションを補償することで、精度・効率よく、長時間コヒーレント積分が可能になるため、積分時間内の目標のレンジセルの移動や演算量の増大を抑えつつ、目標の探知距離を延伸させることができる。
As is apparent from the above, according to the first embodiment, the high-order range migration /
In other words, by compensating the secondary range migration after compensating the primary range migration, it is possible to perform coherent integration over a long period of time with high accuracy. It is possible to extend the target detection distance while suppressing the increase.
実施の形態2.
図8はこの発明の実施の形態2によるレーダ装置の高次レンジマイグレーション・位相補償部11を示す構成図であり、図において、図5と同一符号は同一または相当部分を示すので説明を省略する。
上記実施の形態1では、高次レンジマイグレーション・位相補償部11が図5に示すように、ブロック毎IFFT部41、ドップラーセルマイグレーション補償関数格納部43及びドップラーセルマイグレーション補償関数乗算部45を実装しているものを示したが、図8に示すように、ブロック毎IFFT部41、ドップラーセルマイグレーション補償関数格納部43及びドップラーセルマイグレーション補償関数乗算部45を省略するようにしてもよい。
この場合、ドップラーセルマイグレーションの補償が行われないため(ドップラーセルマイグレーションの補償処理を省略するための条件については、下記の実施の形態5で説明する)、上記実施の形態1と比べて、レンジマイグレーションの補償精度が若干低下するが、処理速度が速くなる効果が得られる。
FIG. 8 is a block diagram showing the high-order range migration /
In the first embodiment, the high-order range migration /
In this case, since compensation for Doppler cell migration is not performed (conditions for omitting compensation processing for Doppler cell migration will be described in the fifth embodiment below), the range is different from that in the first embodiment. Although the migration compensation accuracy is slightly lowered, the effect of increasing the processing speed can be obtained.
実施の形態3.
図9はこの発明の実施の形態3によるレーダ装置の高次レンジマイグレーション・位相補償部11を示す構成図であり、図において、図5と同一符号は同一または相当部分を示すので説明を省略する。
パルス毎IFFT部51はドップラー周波数セル対応1次レンジマイグレーション・位相補償部10により1次のレンジマイグレーションが補償されたドップラー周波数スペクトルをレンジ方向に逆フーリエ変換する処理を実施する。なお、パルス毎IFFT部51はパルス毎逆フーリエ変換手段を構成している。
FIG. 9 is a block diagram showing a high-order range migration /
The pulse-by-
高次位相補償関数格納部52はレンジマイグレーションにおける2次の項に影響される2次の位相変調を補償するための高次位相補償関数(高次項設定部42の設定内容であるαが反映された高次位相補償関数)を格納しているメモリである。
高次位相補償関数格納部52により格納されている高次位相補償関数の詳細については、下記の実施の形態6で述べる。
高次位相補償関数乗算部53は高次位相補償関数格納部52により格納されている高次位相補償関数をパルス毎IFFT部51の逆フーリエ変換結果に乗算する処理を実施する。なお、高次位相補償関数乗算部53は補償関数乗算手段を構成している。
The high-order phase compensation
Details of the high-order phase compensation function stored by the high-order phase compensation
The high-order phase compensation
上記実施の形態1では、ブロック毎IFFT部41が、ドップラー周波数セル対応1次レンジマイグレーション・位相補償部10により1次のレンジマイグレーションが補償されたドップラー周波数スペクトルをヒット方向に逆フーリエ変換してから、ドップラーセルマイグレーション補償関数Ψcmp2(tm,tn,αハット)及び高次レンジマイグレーション位相補償関数Ψcmp2(f,tn,αハット)をブロック毎IFFT部41の逆フーリエ変換結果(ヒット方向に逆フーリエ変換されたドップラー周波数スペクトル)に乗算するものを示したが、図9に示すように、パルス毎IFFT部51が、ドップラー周波数セル対応1次レンジマイグレーション・位相補償部10により1次のレンジマイグレーションが補償されたドップラー周波数スペクトルをレンジ方向に逆フーリエ変換してから、高次位相補償関数格納部52により格納されている高次位相補償関数をパルス毎IFFT部51の逆フーリエ変換結果(レンジ方向に逆フーリエ変換されたドップラー周波数スペクトル)に乗算するようにしても、上記実施の形態1と同様に、2次のレンジマイグレーションを補償することができる。
この場合、補償関数の乗算処理が1回で済むため、上記実施の形態1と比べて、処理速度が速くなる効果が得られる。
In the first embodiment, the block-by-
In this case, since the compensation function is multiplied once, an effect of increasing the processing speed can be obtained as compared with the first embodiment.
この実施の形態3では、ドップラーセルマイグレーションの補償処理及び2次のレンジマイグレーションの補償処理を省略して、高次位相補償関数乗算部53が、2次の位相補償処理だけを実施しているが、ドップラーセルマイグレーションの補償処理及び2次のレンジマイグレーションの補償処理を省略して、2次の位相補償処理を実施する条件については、下記の実施の形態6で説明する。
In the third embodiment, the Doppler cell migration compensation process and the secondary range migration compensation process are omitted, and the high-order phase compensation
実施の形態4.
図10はこの発明の実施の形態4によるレーダ装置の高次レンジマイグレーション・位相補償部11を示す構成図であり、図において、図5及び図9と同一符号は同一または相当部分を示すので説明を省略する。
上記実施の形態1では、ドップラーセルマイグレーション補償関数Ψcmp2(tm,tn,αハット)及び高次レンジマイグレーション位相補償関数Ψcmp2(f,tn,αハット)をブロック毎IFFT部41の逆フーリエ変換結果(ヒット方向に逆フーリエ変換されたドップラー周波数スペクトル)に乗算し、上記実施の形態3では、高次位相補償関数をパルス毎IFFT部51の逆フーリエ変換結果(レンジ方向に逆フーリエ変換されたドップラー周波数スペクトル)に乗算するものを示したが、図10に示すように、1次のレンジマイグレーションが補償されたドップラー周波数スペクトルをレンジ方向及びヒット方向に逆フーリエ変換してから、ドップラーセルマイグレーション補償関数Ψcmp2(tm,tn,αハット)及び高次位相補償関数をレンジ方向及びヒット方向の逆フーリエ変換に乗算するようにしても、上記実施の形態1と同様に、2次のレンジマイグレーションを補償することができる。
10 is a block diagram showing a high-order range migration /
In the first embodiment, the Doppler cell migration compensation function Ψ cmp2 (t m , t n , α hat) and the higher-order range migration phase compensation function Ψ cmp2 (f, t n , α hat) are stored in the block-by-
この実施の形態4では、2次のレンジマイグレーションの補償処理を省略して、ドップラーセルマイグレーションの補償処理と2次の位相補償処理を実施するようにしているが、2次のレンジマイグレーションの補償処理を省略して、ドップラーセルマイグレーションの補償処理と2次の位相補償処理を実施する条件については、下記の実施の形態7で説明する。
In the fourth embodiment, the compensation process for the secondary range migration is omitted and the compensation process for the Doppler cell migration and the secondary phase compensation process are performed. However, the compensation process for the secondary range migration is performed. The conditions for performing the Doppler cell migration compensation process and the secondary phase compensation process will be described in
実施の形態5.
図11はこの発明の実施の形態5によるレーダ装置を示す構成図であり、図において、図2と同一符号は同一または相当部分を示すので説明を省略する。
この実施の形態5のレーダ装置では、パルス毎レンジ圧縮部7の構成要素が、パルス毎FFT部21,22及び複素共役乗算部23であり、図3に示すパルス毎IFFT部24が省略されている。
また、ドップラー周波数セル対応1次レンジマイグレーション・位相補償部10の構成要素が、1次レンジマイグレーション・位相補償関数格納部32及び位相補償関数乗算部33であり、図4に示すパルス毎FFT部31及びパルス毎IFFT部34が省略されている。
高次レンジマイグレーション・位相補償部11の構成は、図5の構成と同一である。
11 is a block diagram showing a radar apparatus according to
In the radar apparatus according to the fifth embodiment, the components of the pulse-by-
Further, the constituent elements of the primary range migration /
The configuration of the high-order range migration /
図11のレーダ装置では、ドップラーセルマイグレーションの補償処理を導入しているが、ドップラーセルマイグレーションの補償処理を省略するようにパラメータ設定することで、演算量の削減を図ることができる。
ドップラーセルマイグレーションを無視する状況下では、2次のレンジマイグレーションも、無視できる程度に小さいと考えられる。
図12はドップラーセルマイグレーションの補償処理を省略しているレーダ装置を示す構成図である。
図12のレーダ装置では、ドップラーセルマイグレーションの補償を省略しているので、図11のレーダ装置におけるブロック毎IFFT部41、ドップラーセルマイグレーション補償関数格納部43、ドップラーセルマイグレーション補償関数乗算部45及びブロック毎FFT部47を削減することができる。
In the radar apparatus of FIG. 11, compensation processing for Doppler cell migration is introduced, but by setting parameters so as to omit compensation processing for Doppler cell migration, the amount of calculation can be reduced.
Under the situation where Doppler cell migration is ignored, the secondary range migration is also considered to be negligibly small.
FIG. 12 is a block diagram showing a radar apparatus in which the compensation process for Doppler cell migration is omitted.
12, the Doppler cell migration compensation is omitted, so the block-by-
ただし、ドップラーセルマイグレーションは、常に無視できる訳ではなく、所定の条件を満足する場合に限り、無視することができる。
以下、ドップラーセルマイグレーションの補償処理を省略するための条件を明示する。
1次のレンジマイグレーションの補償処理が行われた後、ブロック幅Tbについては、再び自由に決定することが可能である。
ただし、パルス幅T0は、最大のドップラー周波数や観測領域を支配するため、ここでは不変の定数であるとする。
However, the Doppler cell migration is not always negligible and can be ignored only when a predetermined condition is satisfied.
In the following, conditions for omitting the Doppler cell migration compensation process will be clarified.
After compensation of the primary range migration occurs, the block width T b, it is possible to freely determine again.
However, since the pulse width T 0 dominates the maximum Doppler frequency and the observation region, it is assumed here to be an invariant constant.
この実施の形態5では、位相補償関数乗算部33によって、1次レンジマイグレーション位相補償関数Ψcmp1(f,fdm,tn)が乗算されたのち、パルスの分割数であるMや、ブロックの分割数であるNを変更できるものとし、変更後のM,NをM’,N’と表記する。
M’を可変とすることによるドップラーセルマイグレーションの補償処理を省略するための条件としては、ブロックサイズ変更後の周波数分解能1/M’T0について、式(43)と同様に、次の関係が成り立てばよい。
In the fifth embodiment, after being multiplied by the primary range migration phase compensation function Ψ cmp1 (f, f dm , t n ) by the phase compensation
As a condition for omitting the Doppler cell migration compensation process by making M ′ variable, the following relationship is obtained for the
以上で明らかなように、ブロックサイズを可変として、パラメータを最適に設定することにより、ドップラーセルマイグレーションを省略することができ、その結果、演算量を削減した高次のレンジマイグレーション補償が可能になる効果がある。 As is apparent from the above, Doppler cell migration can be omitted by making the block size variable and optimally setting the parameters, and as a result, higher-order range migration compensation with a reduced amount of computation becomes possible. effective.
実施の形態6.
上記実施の形態5では、ドップラーセルマイグレーションの補償処理を省略して、位相補償関数乗算部46が、高次レンジマイグレーション位相補償関数Ψcmp2(f,tn,αハット)を乗算することで、2次のレンジマイグレーションを補償するものを示したが、ドップラーセルマイグレーションの補償処理及び2次のレンジマイグレーションの補償処理を省略して、高次位相補償関数乗算部53が、高次位相補償関数Ψcmp2(tn,αハット)を乗算することで、2次の位相補償処理だけを実施するようにしてもよい。
図13はこの発明の実施の形態6によるレーダ装置を示す構成図である。
In the fifth embodiment, the Doppler cell migration compensation process is omitted, and the phase compensation
13 is a block diagram showing a radar apparatus according to
まず、レンジ方向のデータ量は、ヒット方向やブロック方向と比べて、高いオーダである。
レンジ方向に逆フーリエ変換するパルス毎IFFT部48が、高次項設定部42の後段に設置されているのは(例えば、図11を参照)、2次のレンジマイグレーションを補償するためであるが、信号の比帯域が十分に小さければ、2次のレンジマイグレーションは無視できる程度に小さいと考えられる。
2次のレンジマイグレーションの補償処理を省略することができる条件は、次に示す通りである。
First, the amount of data in the range direction is higher than that in the hit direction or the block direction.
The reason why the pulse-by-
The conditions under which the secondary range migration compensation process can be omitted are as follows.
したがって、式(45)及び式(46)に示す条件の双方を満足する場合、図13に示すように、ドップラーセルマイグレーションの補償処理及び2次のレンジマイグレーションの補償処理を省略して、2次の位相補償処理だけを実施する(高次位相補償関数乗算部53が、高次位相補償関数Ψcmp2(tn,αハット)を乗算する処理を実施する)ようにすることができる。
これにより、2次のレンジマイグレーションの補償処理に要する畳み込み積分の演算量を削減することができる。
Therefore, when both of the conditions shown in Expression (45) and Expression (46) are satisfied, the Doppler cell migration compensation process and the secondary range migration compensation process are omitted as shown in FIG. (The high-order phase compensation
Thereby, the amount of calculation of convolution integral required for the compensation process of the secondary range migration can be reduced.
図13のレーダ装置では、ドップラーセルマイグレーションの補償処理及び2次のレンジマイグレーションの補償処理を省略することで、反復処理における更なる演算量の削減を考慮している。
ただし、経験的には、1次のレンジマイグレーションの補償処理のみで十分に補償されていることを確認しているが、ドップラーセルマイグレーションの補償処理を省略することによる問題もある。
δfdの同定まではできないため、ドップラーセルマイグレーションが2セルにまたがって発生する可能性がある。言い換えると、ドップラーセルマイグレーションを最大でも2セルまでに抑えることは可能である。
In the radar apparatus of FIG. 13, further reduction of the calculation amount in the iterative process is considered by omitting the Doppler cell migration compensation process and the secondary range migration compensation process.
However, empirically, it has been confirmed that the compensation is sufficiently performed only by the compensation process for the first-order range migration, but there is a problem that the compensation process for the Doppler cell migration is omitted.
It can not be until the identification of delta] f d, which may doppler cell migration occurs across the two cells. In other words, it is possible to suppress the Doppler cell migration to a maximum of 2 cells.
これにより、目標電力Prevが不変(レーダ方程式について変化がない)という仮定の下では、最悪の場合で3dBの損失が生じる。
しかし、実際のところ、目標移動によるRCS(Radar Cross Section)揺らぎ、距離変化によって、受信される電力にも増減(シンチレーション)を伴うため、データ長に対して実効的な積分時間が短くなる傾向がある。
そのため、ドップラーセルマイグレーションの補償処理を省略した場合でも、3dBまでの積分ロスは発生しないと考えられる。問題があるとすれば、隣のドップラーセルに偽像が発生する可能性を含むことである。
演算量の削減に伴ってデメリットもあるが、反復処理内の処理量の削減が、全体の演算量を大きく左右するため、演算量の観点から言えば、有効性が高いことが期待される。
This causes a 3 dB loss in the worst case under the assumption that the target power Prev is unchanged (no change with respect to the radar equation).
However, in reality, the RCS (Radar Cross Section) fluctuation due to the target movement and the received power are accompanied by an increase / decrease (scintillation) in the received power, so that the effective integration time tends to be shortened with respect to the data length. is there.
Therefore, even when the compensation process for Doppler cell migration is omitted, it is considered that an integration loss up to 3 dB does not occur. If there is a problem, it includes the possibility of a false image occurring in the adjacent Doppler cell.
Although there is a demerit associated with the reduction in the amount of computation, since the reduction in the amount of processing in the iterative process greatly affects the overall amount of computation, it is expected to be highly effective from the viewpoint of the amount of computation.
実施の形態7.
上記実施の形態5では、ドップラーセルマイグレーションの補償処理を省略して、位相補償関数乗算部46が、高次レンジマイグレーション位相補償関数Ψcmp2(f,tn,αハット)を乗算することで、2次のレンジマイグレーションを補償し、上記実施の形態6では、ドップラーセルマイグレーションの補償処理及び2次のレンジマイグレーションの補償処理を省略して、高次位相補償関数乗算部53が、高次位相補償関数Ψcmp2(tn,αハット)を乗算することで、2次の位相補償処理だけを実施するものを示したが、2次のレンジマイグレーションの補償処理を省略して、ドップラーセルマイグレーションの補償処理と2次の位相補償処理を実施するようにしてもよい。
図14はこの発明の実施の形態6によるレーダ装置を示す構成図である。
In the fifth embodiment, the Doppler cell migration compensation process is omitted, and the phase compensation
FIG. 14 is a block diagram showing a radar apparatus according to
2次のレンジマイグレーションの補償処理を省略して、ドップラーセルマイグレーションの補償処理と2次の位相補償処理を実施する条件は、上記の式(46)を満足することである。
図14のレーダ装置は、探知レーダとして扱う場合に不可避であるレンジ方向のデータ点数の増大を考慮した演算量削減方式を適用するものである。
A condition for performing the Doppler cell migration compensation process and the secondary phase compensation process while omitting the compensation process for the secondary range migration is to satisfy the above equation (46).
The radar apparatus of FIG. 14 applies a calculation amount reduction method that takes into account an increase in the number of data points in the range direction, which is unavoidable when handled as a detection radar.
なお、本願発明はその発明の範囲内において、各実施の形態の自由な組み合わせ、あるいは各実施の形態の任意の構成要素の変形、もしくは各実施の形態において任意の構成要素の省略が可能である。 In the present invention, within the scope of the invention, any combination of the embodiments, or any modification of any component in each embodiment, or omission of any component in each embodiment is possible. .
1 電波発信源、2 直接波受信用アンテナ(電波受信手段)、3 散乱波受信用アンテナ(電波受信手段)、4 受信局(電波受信手段)、5 直接波受信機、6 散乱波受信機、7 パルス毎レンジ圧縮部(パルス毎レンジ圧縮手段)、8 クラッタ抑圧部、9 ブロック毎ドップラー処理部(ブロック毎ドップラー処理手段)、10 ドップラー周波数セル対応1次レンジマイグレーション・位相補償部(1次レンジマイグレーション補償手段)、11 高次レンジマイグレーション・位相補償部(高次レンジマイグレーション補償手段)、12 ブロック方向ドップラー処理部(ブロック方向ドップラー処理手段)、13 目標検出部(目標検出手段)、21 パルス毎FFT部、22 パルス毎FFT部、23 複素共役乗算部、24 パルス毎IFFT部、31 パルス毎FFT部(パルス毎フーリエ変換手段)、32 1次レンジマイグレーション・位相補償関数格納部、33 位相補償関数乗算部(位相補償関数乗算手段)、34 パルス毎IFFT部(パルス毎逆フーリエ変換手段)、41 ブロック毎IFFT部(ブロック毎逆フーリエ変換手段)、42 高次項設定部、43 ドップラーセルマイグレーション補償関数格納部、44 高次レンジマイグレーション・位相補償関数格納部、45 ドップラーセルマイグレーション補償関数乗算部(第1の補償関数乗算手段)、46 位相補償関数乗算部(第2の補償関数乗算手段)、47 ブロック毎FFT部(ブロック毎フーリエ変換手段)、48 パルス毎IFFT部(パルス毎逆フーリエ変換手段)、51 パルス毎IFFT部(パルス毎逆フーリエ変換手段)、52 高次位相補償関数格納部、53 高次位相補償関数乗算部(補償関数乗算手段)。 1 radio wave source, 2 direct wave receiving antenna (radio wave receiving means), 3 scattered wave receiving antenna (radio wave receiving means), 4 receiving station (radio wave receiving means), 5 direct wave receiver, 6 scattered wave receiver, 7 Pulse-by-pulse range compression section (pulse-by-pulse range compression means) 8 Clutter suppression section 9 Block-by-block Doppler processing section (block-by-block Doppler processing means) Migration compensation means), 11 high-order range migration / phase compensation section (high-order range migration compensation means), 12 block direction Doppler processing section (block direction Doppler processing means), 13 target detection section (target detection means), 21 per pulse FFT unit, 22 pulse-by-pulse FFT unit, 23 complex conjugate multiplier, 24 pulses IFFT unit, 31 FFT unit per pulse (Fourier transform unit per pulse), 32 Primary range migration / phase compensation function storage unit, 33 Phase compensation function multiplication unit (phase compensation function multiplication unit), 34 IFFT unit per pulse (pulse Every inverse Fourier transform means), 41 each block IFFT section (each block inverse Fourier transform means), 42 higher order term setting section, 43 Doppler cell migration compensation function storage section, 44 higher order range migration / phase compensation function storage section, 45 Doppler Cell migration compensation function multiplication unit (first compensation function multiplication unit), 46 Phase compensation function multiplication unit (second compensation function multiplication unit), 47 FFT unit for each block (Fourier transform unit for each block), 48 IFFT unit for each pulse (Inverse Fourier transform means per pulse), 51 IFF per pulse Part (pulse per inverse Fourier transform unit) 52 high-order phase compensation function storage section, 53 high-order phase compensation function multiplier (compensation function multiplication means).
Claims (10)
上記電波受信手段により受信された直接波の信号及び散乱波の信号をパルスに分割し、各々のパルスについて、上記直接波の信号と上記散乱波の信号との相互相関を求めることで、パルス毎のレンジプロフィールを生成するパルス毎レンジ圧縮手段と、
上記パルス毎レンジ圧縮手段により生成されたパルス毎のレンジプロフィールをブロック単位にまとめ、複数のパルスからなるブロック単位のレンジプロフィールをヒット方向にフーリエ変換することで、ブロック毎のドップラー周波数スペクトルを算出するブロック毎ドップラー処理手段と、
上記ブロック毎ドップラー処理手段により算出されたドップラー周波数スペクトルの中で、ブロック方向に発生している1次のレンジマイグレーションを補償する1次レンジマイグレーション補償手段と、
上記1次レンジマイグレーション補償手段により1次のレンジマイグレーションが補償されたドップラー周波数スペクトルの中で、ブロック方向に発生している2次以上のレンジマイグレーションを補償する高次レンジマイグレーション補償手段と、
上記高次レンジマイグレーション補償手段により2次以上のレンジマイグレーションが補償されたドップラー周波数スペクトルをブロック方向にフーリエ変換することで、ドップラー周波数スペクトルを算出するブロック方向ドップラー処理手段と、
上記ブロック方向ドップラー処理手段により算出されたドップラー周波数スペクトルから目標を検出する目標検出手段とを備え、
上記1次レンジマイグレーション補償手段は、
上記ブロック毎ドップラー処理手段により算出されたドップラー周波数スペクトルをレンジ方向にフーリエ変換するパルス毎フーリエ変換手段と、
上記ブロック毎ドップラー処理手段により算出されたドップラー周波数スペクトルにおける各々のドップラー周波数セルに対応するブロック方向の1次のレンジマイグレーションの量に相当する位相変化を補償するための1次レンジマイグレーション位相補償関数を上記パルス毎フーリエ変換手段のフーリエ変換結果に乗算する位相補償関数乗算手段と、
上記位相補償関数乗算手段の乗算結果をレンジ方向に逆フーリエ変換するパルス毎逆フーリエ変換手段と
から構成されていることを特徴とするレーダ装置。 A radio wave receiving means for receiving a direct wave of a radio wave transmitted from a radio wave source and receiving a scattered wave in which the radio wave transmitted from the radio wave source is scattered by a target;
The direct wave signal and the scattered wave signal received by the radio wave receiving means are divided into pulses, and for each pulse, the cross correlation between the direct wave signal and the scattered wave signal is obtained for each pulse. A pulse-by-pulse range compression means for generating a range profile of
The Doppler frequency spectrum for each block is calculated by summarizing the range profile for each pulse generated by the above-mentioned pulse-by-pulse range compressing unit into block units, and Fourier-transforming the range profile of the block unit consisting of a plurality of pulses in the hit direction. Block-by-block Doppler processing means;
Primary range migration compensation means for compensating primary range migration occurring in the block direction in the Doppler frequency spectrum calculated by the block-by-block Doppler processing means;
High-order range migration compensation means for compensating for secondary or higher range migration occurring in the block direction in the Doppler frequency spectrum in which the primary range migration is compensated by the primary range migration compensation means;
A block direction Doppler processing means for calculating a Doppler frequency spectrum by Fourier-transforming the Doppler frequency spectrum in which second-order or higher range migration has been compensated by the high-order range migration compensation means in the block direction;
E Bei a target detection means for detecting a target from the Doppler frequency spectrum calculated by the block direction Doppler processing means,
The primary range migration compensation means is:
A pulse-by-pulse Fourier transform means for Fourier-transforming the Doppler frequency spectrum calculated by the block-by-block Doppler processing means in the range direction;
A primary range migration phase compensation function for compensating for a phase change corresponding to the amount of primary range migration in the block direction corresponding to each Doppler frequency cell in the Doppler frequency spectrum calculated by the block-by-block Doppler processing means. Phase compensation function multiplying means for multiplying the Fourier transform result of the pulse-by-pulse Fourier transform means;
A pulse-by-pulse inverse Fourier transform means for performing an inverse Fourier transform on the multiplication result of the phase compensation function multiplying means in the range direction;
A radar apparatus comprising:
上記電波受信手段により受信された直接波の信号及び散乱波の信号をパルスに分割し、各々のパルスについて、上記直接波の信号と上記散乱波の信号との相互相関を求めることで、パルス毎のレンジプロフィールを生成するパルス毎レンジ圧縮手段と、
上記パルス毎レンジ圧縮手段により生成されたパルス毎のレンジプロフィールをブロック単位にまとめ、複数のパルスからなるブロック単位のレンジプロフィールをヒット方向にフーリエ変換することで、ブロック毎のドップラー周波数スペクトルを算出するブロック毎ドップラー処理手段と、
上記ブロック毎ドップラー処理手段により算出されたドップラー周波数スペクトルの中で、ブロック方向に発生している1次のレンジマイグレーションを補償する1次レンジマイグレーション補償手段と、
上記1次レンジマイグレーション補償手段により1次のレンジマイグレーションが補償されたドップラー周波数スペクトルの中で、ブロック方向に発生している2次以上のレンジマイグレーションを補償する高次レンジマイグレーション補償手段と、
上記高次レンジマイグレーション補償手段により2次以上のレンジマイグレーションが補償されたドップラー周波数スペクトルをブロック方向にフーリエ変換することで、ドップラー周波数スペクトルを算出するブロック方向ドップラー処理手段と、
上記ブロック方向ドップラー処理手段により算出されたドップラー周波数スペクトルから目標を検出する目標検出手段とを備え、
上記高次レンジマイグレーション補償手段は、
上記1次レンジマイグレーション補償手段により補償されたドップラー周波数スペクトルをヒット方向に逆フーリエ変換するブロック毎逆フーリエ変換手段と、
レンジマイグレーションにおける2次以上の項に影響されるドップラーセルマイグレーションを補償するためのドップラーセルマイグレーション補償関数を上記ブロック毎逆フーリエ変換手段の逆フーリエ変換結果に乗算する第1の補償関数乗算手段と、
2次以上のレンジマイグレーションを補償するための高次レンジマイグレーション位相補償関数を上記第1の補償関数乗算手段の乗算結果に乗算する第2の補償関数乗算手段と、
上記第2の補償関数乗算手段の乗算結果をヒット方向にフーリエ変換するブロック毎フーリエ変換手段と、
上記ブロック毎フーリエ変換手段のフーリエ変換結果をレンジ方向に逆フーリエ変換するパルス毎逆フーリエ変換手段と
から構成されていることを特徴とするレーダ装置。 A radio wave receiving means for receiving a direct wave of a radio wave transmitted from a radio wave source and receiving a scattered wave in which the radio wave transmitted from the radio wave source is scattered by a target;
The direct wave signal and the scattered wave signal received by the radio wave receiving means are divided into pulses, and for each pulse, the cross correlation between the direct wave signal and the scattered wave signal is obtained for each pulse. A pulse-by-pulse range compression means for generating a range profile of
The Doppler frequency spectrum for each block is calculated by summarizing the range profile for each pulse generated by the above-mentioned pulse-by-pulse range compressing unit into block units, and Fourier-transforming the range profile of the block unit consisting of a plurality of pulses in the hit direction. Block-by-block Doppler processing means;
Primary range migration compensation means for compensating primary range migration occurring in the block direction in the Doppler frequency spectrum calculated by the block-by-block Doppler processing means;
High-order range migration compensation means for compensating for secondary or higher range migration occurring in the block direction in the Doppler frequency spectrum in which the primary range migration is compensated by the primary range migration compensation means;
A block direction Doppler processing means for calculating a Doppler frequency spectrum by Fourier-transforming the Doppler frequency spectrum in which second-order or higher range migration has been compensated by the high-order range migration compensation means in the block direction;
Target detection means for detecting a target from the Doppler frequency spectrum calculated by the block direction Doppler processing means,
The higher range migration compensation means is
And inverse Fourier transform means for each block to inverse Fourier transform the Doppler frequency spectrum to the hit direction, which is compensated by the primary-range migration compensation means,
First compensation function multiplying means for multiplying the inverse Fourier transform result of the block-by-block inverse Fourier transform means by a Doppler cell migration compensation function for compensating for Doppler cell migration affected by a second-order or higher term in range migration;
Second compensation function multiplying means for multiplying a multiplication result of the first compensation function multiplying means by a higher order range migration phase compensation function for compensating for second or higher order range migration;
Block-by-block Fourier transform means for Fourier transforming the multiplication result of the second compensation function multiplying means in the hit direction;
A radar apparatus comprising: a pulse-by-pulse inverse Fourier transform means for performing an inverse Fourier transform on a Fourier transform result of the block-by-block Fourier transform means in a range direction.
上記電波受信手段により受信された直接波の信号及び散乱波の信号をパルスに分割し、各々のパルスについて、上記直接波の信号と上記散乱波の信号との相互相関を求めることで、パルス毎のレンジプロフィールを生成するパルス毎レンジ圧縮手段と、
上記パルス毎レンジ圧縮手段により生成されたパルス毎のレンジプロフィールをブロック単位にまとめ、複数のパルスからなるブロック単位のレンジプロフィールをヒット方向にフーリエ変換することで、ブロック毎のドップラー周波数スペクトルを算出するブロック毎ドップラー処理手段と、
上記ブロック毎ドップラー処理手段により算出されたドップラー周波数スペクトルの中で、ブロック方向に発生している1次のレンジマイグレーションを補償する1次レンジマイグレーション補償手段と、
上記1次レンジマイグレーション補償手段により1次のレンジマイグレーションが補償されたドップラー周波数スペクトルの中で、ブロック方向に発生している2次以上のレンジマイグレーションを補償する高次レンジマイグレーション補償手段と、
上記高次レンジマイグレーション補償手段により2次以上のレンジマイグレーションが補償されたドップラー周波数スペクトルをブロック方向にフーリエ変換することで、ドップラー周波数スペクトルを算出するブロック方向ドップラー処理手段と、
上記ブロック方向ドップラー処理手段により算出されたドップラー周波数スペクトルから目標を検出する目標検出手段とを備え、
上記高次レンジマイグレーション補償手段は、
2次以上のレンジマイグレーションを補償するための高次レンジマイグレーション位相補償関数を1次レンジマイグレーション補償手段により補償されたドップラー周波数スペクトルに乗算する補償関数乗算手段と、
上記補償関数乗算手段の乗算結果をヒット方向にフーリエ変換するブロック毎フーリエ変換手段と、
上記ブロック毎フーリエ変換手段のフーリエ変換結果をレンジ方向に逆フーリエ変換するパルス毎逆フーリエ変換手段と
から構成されていることを特徴とするレーダ装置。 A radio wave receiving means for receiving a direct wave of a radio wave transmitted from a radio wave source and receiving a scattered wave in which the radio wave transmitted from the radio wave source is scattered by a target;
The direct wave signal and the scattered wave signal received by the radio wave receiving means are divided into pulses, and for each pulse, the cross correlation between the direct wave signal and the scattered wave signal is obtained for each pulse. A pulse-by-pulse range compression means for generating a range profile of
The Doppler frequency spectrum for each block is calculated by summarizing the range profile for each pulse generated by the above-mentioned pulse-by-pulse range compressing unit into block units, and Fourier-transforming the range profile of the block unit consisting of a plurality of pulses in the hit direction. Block-by-block Doppler processing means;
Primary range migration compensation means for compensating primary range migration occurring in the block direction in the Doppler frequency spectrum calculated by the block-by-block Doppler processing means;
High-order range migration compensation means for compensating for secondary or higher range migration occurring in the block direction in the Doppler frequency spectrum in which the primary range migration is compensated by the primary range migration compensation means;
A block direction Doppler processing means for calculating a Doppler frequency spectrum by Fourier-transforming the Doppler frequency spectrum in which second-order or higher range migration has been compensated by the high-order range migration compensation means in the block direction;
Target detection means for detecting a target from the Doppler frequency spectrum calculated by the block direction Doppler processing means,
The higher range migration compensation means is
Compensation function multiplication means for multiplying the Doppler frequency spectrum compensated by the primary range migration compensation means by a higher order range migration phase compensation function for compensating for the second or higher order range migration;
Block-by-block Fourier transform means for Fourier transforming the multiplication result of the compensation function multiplying means in the hit direction;
A radar apparatus comprising: a pulse-by-pulse inverse Fourier transform means for performing an inverse Fourier transform on a Fourier transform result of the block-by-block Fourier transform means in a range direction.
上記電波受信手段により受信された直接波の信号及び散乱波の信号をパルスに分割し、各々のパルスについて、上記直接波の信号と上記散乱波の信号との相互相関を求めることで、パルス毎のレンジプロフィールを生成するパルス毎レンジ圧縮手段と、
上記パルス毎レンジ圧縮手段により生成されたパルス毎のレンジプロフィールをブロック単位にまとめ、複数のパルスからなるブロック単位のレンジプロフィールをヒット方向にフーリエ変換することで、ブロック毎のドップラー周波数スペクトルを算出するブロック毎ドップラー処理手段と、
上記ブロック毎ドップラー処理手段により算出されたドップラー周波数スペクトルの中で、ブロック方向に発生している1次のレンジマイグレーションを補償する1次レンジマイグレーション補償手段と、
上記1次レンジマイグレーション補償手段により1次のレンジマイグレーションが補償されたドップラー周波数スペクトルの中で、ブロック方向に発生している2次以上のレンジマイグレーションを補償する高次レンジマイグレーション補償手段と、
上記高次レンジマイグレーション補償手段により2次以上のレンジマイグレーションが補償されたドップラー周波数スペクトルをブロック方向にフーリエ変換することで、ドップラー周波数スペクトルを算出するブロック方向ドップラー処理手段と、
上記ブロック方向ドップラー処理手段により算出されたドップラー周波数スペクトルから目標を検出する目標検出手段とを備え、
上記高次レンジマイグレーション補償手段は、
上記1次レンジマイグレーション補償手段により補償されたドップラー周波数スペクトルをレンジ方向に逆フーリエ変換するパルス毎逆フーリエ変換手段と、
レンジマイグレーションにおける2次以上の項に影響される2次以上の位相変調を補償するための高次位相補償関数を上記パルス毎逆フーリエ変換手段の逆フーリエ変換結果に乗算する補償関数乗算手段と、
上記補償関数乗算手段の乗算結果をヒット方向にフーリエ変換するブロック毎フーリエ変換手段と
から構成されていることを特徴とするレーダ装置。 A radio wave receiving means for receiving a direct wave of a radio wave transmitted from a radio wave source and receiving a scattered wave in which the radio wave transmitted from the radio wave source is scattered by a target;
The direct wave signal and the scattered wave signal received by the radio wave receiving means are divided into pulses, and for each pulse, the cross correlation between the direct wave signal and the scattered wave signal is obtained for each pulse. A pulse-by-pulse range compression means for generating a range profile of
The Doppler frequency spectrum for each block is calculated by summarizing the range profile for each pulse generated by the above-mentioned pulse-by-pulse range compressing unit into block units, and Fourier-transforming the range profile of the block unit consisting of a plurality of pulses in the hit direction. Block-by-block Doppler processing means;
Primary range migration compensation means for compensating primary range migration occurring in the block direction in the Doppler frequency spectrum calculated by the block-by-block Doppler processing means;
High-order range migration compensation means for compensating for secondary or higher range migration occurring in the block direction in the Doppler frequency spectrum in which the primary range migration is compensated by the primary range migration compensation means;
A block direction Doppler processing means for calculating a Doppler frequency spectrum by Fourier-transforming the Doppler frequency spectrum in which second-order or higher range migration has been compensated by the high-order range migration compensation means in the block direction;
Target detection means for detecting a target from the Doppler frequency spectrum calculated by the block direction Doppler processing means,
The higher range migration compensation means is
A pulse each inverse Fourier transform means for inverse Fourier transform of the Doppler frequency spectrum to the range direction which is compensated by the primary-range migration compensation means,
Compensation function multiplication means for multiplying the inverse Fourier transform result of the inverse Fourier transform means for each pulse by a higher-order phase compensation function for compensating for the second or higher order phase modulation influenced by the second or higher order term in the range migration;
A radar apparatus comprising: a block-by-block Fourier transform unit that Fourier-transforms a multiplication result of the compensation function multiplying unit in a hit direction.
上記電波受信手段により受信された直接波の信号及び散乱波の信号をパルスに分割し、各々のパルスについて、上記直接波の信号と上記散乱波の信号との相互相関を求めることで、パルス毎のレンジプロフィールを生成するパルス毎レンジ圧縮手段と、
上記パルス毎レンジ圧縮手段により生成されたパルス毎のレンジプロフィールをブロック単位にまとめ、複数のパルスからなるブロック単位のレンジプロフィールをヒット方向にフーリエ変換することで、ブロック毎のドップラー周波数スペクトルを算出するブロック毎ドップラー処理手段と、
上記ブロック毎ドップラー処理手段により算出されたドップラー周波数スペクトルの中で、ブロック方向に発生している1次のレンジマイグレーションを補償する1次レンジマイグレーション補償手段と、
上記1次レンジマイグレーション補償手段により1次のレンジマイグレーションが補償されたドップラー周波数スペクトルの中で、ブロック方向に発生している2次以上のレンジマイグレーションを補償する高次レンジマイグレーション補償手段と、
上記高次レンジマイグレーション補償手段により2次以上のレンジマイグレーションが補償されたドップラー周波数スペクトルをブロック方向にフーリエ変換することで、ドップラー周波数スペクトルを算出するブロック方向ドップラー処理手段と、
上記ブロック方向ドップラー処理手段により算出されたドップラー周波数スペクトルから目標を検出する目標検出手段とを備え、
上記高次レンジマイグレーション補償手段は、
上記1次レンジマイグレーション補償手段により補償されたドップラー周波数スペクトルをレンジ方向に逆フーリエ変換するパルス毎逆フーリエ変換手段と、
上記パルス毎逆フーリエ変換手段の逆フーリエ変換結果をヒット方向に逆フーリエ変換するブロック毎逆フーリエ変換手段と、
レンジマイグレーションにおける2次以上の項に影響されるドップラーセルマイグレーションを補償するためのドップラーセルマイグレーション補償関数を上記ブロック毎逆フーリエ変換手段の逆フーリエ変換結果に乗算する第1の補償関数乗算手段と、
レンジマイグレーションにおける2次以上の項に影響される2次以上の位相変調を補償するための高次位相補償関数を上記第1の補償関数乗算手段の乗算結果に乗算する第2の補償関数乗算手段と、
上記第2の補償関数乗算手段をヒット方向にフーリエ変換するブロック毎フーリエ変換手段と
から構成されていることを特徴とするレーダ装置。 A radio wave receiving means for receiving a direct wave of a radio wave transmitted from a radio wave source and receiving a scattered wave in which the radio wave transmitted from the radio wave source is scattered by a target;
The direct wave signal and the scattered wave signal received by the radio wave receiving means are divided into pulses, and for each pulse, the cross correlation between the direct wave signal and the scattered wave signal is obtained for each pulse. A pulse-by-pulse range compression means for generating a range profile of
The Doppler frequency spectrum for each block is calculated by summarizing the range profile for each pulse generated by the above-mentioned pulse-by-pulse range compressing unit into block units, and Fourier-transforming the range profile of the block unit consisting of a plurality of pulses in the hit direction. Block-by-block Doppler processing means;
Primary range migration compensation means for compensating primary range migration occurring in the block direction in the Doppler frequency spectrum calculated by the block-by-block Doppler processing means;
High-order range migration compensation means for compensating for secondary or higher range migration occurring in the block direction in the Doppler frequency spectrum in which the primary range migration is compensated by the primary range migration compensation means;
A block direction Doppler processing means for calculating a Doppler frequency spectrum by Fourier-transforming the Doppler frequency spectrum in which second-order or higher range migration has been compensated by the high-order range migration compensation means in the block direction;
Target detection means for detecting a target from the Doppler frequency spectrum calculated by the block direction Doppler processing means,
The higher range migration compensation means is
A pulse each inverse Fourier transform means for inverse Fourier transform of the Doppler frequency spectrum to the range direction which is compensated by the primary-range migration compensation means,
Block inverse Fourier transform means for inverse Fourier transforming the inverse Fourier transform result of the pulse-by-pulse inverse Fourier transform means in the hit direction;
First compensation function multiplying means for multiplying the inverse Fourier transform result of the block-by-block inverse Fourier transform means by a Doppler cell migration compensation function for compensating for Doppler cell migration affected by a second-order or higher term in range migration;
Second compensation function multiplying means for multiplying the multiplication result of the first compensation function multiplying means by a higher order phase compensation function for compensating for the second or higher order phase modulation affected by the second or higher order term in the range migration. When,
A radar apparatus comprising: a block-by-block Fourier transform unit that Fourier-transforms the second compensation function multiplying unit in a hit direction.
上記電波受信手段により受信された直接波の信号及び散乱波の信号をパルスに分割し、各々のパルスについて、上記直接波の信号と上記散乱波の信号との相互相関を求めることで、パルス毎のレンジプロフィールを生成するパルス毎レンジ圧縮手段と、
上記パルス毎レンジ圧縮手段により生成されたパルス毎のレンジプロフィールをブロック単位にまとめ、複数のパルスからなるブロック単位のレンジプロフィールをヒット方向にフーリエ変換することで、ブロック毎のドップラー周波数スペクトルを算出するブロック毎ドップラー処理手段と、
上記ブロック毎ドップラー処理手段により算出されたドップラー周波数スペクトルにおける各々のドップラー周波数セルに対応するブロック方向の1次のレンジマイグレーションの量に相当する位相変化を補償するための1次レンジマイグレーション位相補償関数を上記ドップラー周波数スペクトルに乗算する位相補償関数乗算手段と、
上記位相補償関数乗算手段の乗算結果をヒット方向に逆フーリエ変換するブロック毎逆フーリエ変換手段と、
レンジマイグレーションにおける2次以上の項に影響されるドップラーセルマイグレーションを補償するためのドップラーセルマイグレーション補償関数を上記ブロック毎逆フーリエ変換手段の逆フーリエ変換結果に乗算する第1の補償関数乗算手段と、
2次以上のレンジマイグレーションを補償するための高次レンジマイグレーション位相補償関数を上記第1の補償関数乗算手段の乗算結果に乗算する第2の補償関数乗算手段と、
上記第2の補償関数乗算手段の乗算結果をヒット方向にフーリエ変換するブロック毎フーリエ変換手段と、
上記ブロック毎フーリエ変換手段のフーリエ変換結果をレンジ方向に逆フーリエ変換するパルス毎逆フーリエ変換手段と
上記パルス毎逆フーリエ変換手段の逆フーリエ変換結果をブロック方向にフーリエ変換することで、ドップラー周波数スペクトルを算出するブロック方向ドップラー処理手段と、
上記ブロック方向ドップラー処理手段により算出されたドップラー周波数スペクトルから目標を検出する目標検出手段と
を備えたレーダ装置。 A radio wave receiving means for receiving a direct wave of a radio wave transmitted from a radio wave source and receiving a scattered wave in which the radio wave transmitted from the radio wave source is scattered by a target;
The direct wave signal and the scattered wave signal received by the radio wave receiving means are divided into pulses, and for each pulse, the cross correlation between the direct wave signal and the scattered wave signal is obtained for each pulse. A pulse-by-pulse range compression means for generating a range profile of
The Doppler frequency spectrum for each block is calculated by summarizing the range profile for each pulse generated by the above-mentioned pulse-by-pulse range compressing unit into block units, and Fourier-transforming the range profile of the block unit consisting of a plurality of pulses in the hit direction. Block-by-block Doppler processing means;
A primary range migration phase compensation function for compensating for a phase change corresponding to the amount of primary range migration in the block direction corresponding to each Doppler frequency cell in the Doppler frequency spectrum calculated by the block-by-block Doppler processing means. Phase compensation function multiplication means for multiplying the Doppler frequency spectrum;
Block-by-block inverse Fourier transform means for inverse Fourier transforming the multiplication result of the phase compensation function multiplying means in the hit direction;
First compensation function multiplying means for multiplying the inverse Fourier transform result of the block-by-block inverse Fourier transform means by a Doppler cell migration compensation function for compensating for Doppler cell migration affected by a second-order or higher term in range migration;
Second compensation function multiplying means for multiplying a multiplication result of the first compensation function multiplying means by a higher order range migration phase compensation function for compensating for second or higher order range migration;
Block-by-block Fourier transform means for Fourier transforming the multiplication result of the second compensation function multiplying means in the hit direction;
A pulse-by-pulse inverse Fourier transform means for performing an inverse Fourier transform on the Fourier transform result of the block-by-block Fourier transform means in a range direction, and a Doppler frequency spectrum by Fourier-transforming the inverse Fourier transform result of the pulse-by-pulse inverse Fourier transform means in a block direction. Block direction Doppler processing means for calculating
A radar apparatus comprising: target detection means for detecting a target from the Doppler frequency spectrum calculated by the block direction Doppler processing means.
上記電波受信手段により受信された直接波の信号及び散乱波の信号をパルスに分割し、各々のパルスについて、上記直接波の信号と上記散乱波の信号との相互相関を求めることで、パルス毎のレンジプロフィールを生成するパルス毎レンジ圧縮手段と、
上記パルス毎レンジ圧縮手段により生成されたパルス毎のレンジプロフィールをブロック単位にまとめ、複数のパルスからなるブロック単位のレンジプロフィールをヒット方向にフーリエ変換することで、ブロック毎のドップラー周波数スペクトルを算出するブロック毎ドップラー処理手段と、
上記ブロック毎ドップラー処理手段により算出されたドップラー周波数スペクトルにおける各々のドップラー周波数セルに対応するブロック方向の1次のレンジマイグレーションの量に相当する位相変化を補償するための1次レンジマイグレーション位相補償関数を上記ドップラー周波数スペクトルに乗算する位相補償関数乗算手段と、
2次以上のレンジマイグレーションを補償するための高次レンジマイグレーション位相補償関数を上記位相補償関数乗算手段の乗算結果に乗算する補償関数乗算手段と、
上記補償関数乗算手段の乗算結果をレンジ方向に逆フーリエ変換するパルス毎逆フーリエ変換手段と、
上記パルス毎逆フーリエ変換手段の逆フーリエ変換結果をブロック方向にフーリエ変換することで、ドップラー周波数スペクトルを算出するブロック方向ドップラー処理手段と、
上記ブロック方向ドップラー処理手段により算出されたドップラー周波数スペクトルから目標を検出する目標検出手段と
を備えたレーダ装置。 A radio wave receiving means for receiving a direct wave of a radio wave transmitted from a radio wave source and receiving a scattered wave in which the radio wave transmitted from the radio wave source is scattered by a target;
The direct wave signal and the scattered wave signal received by the radio wave receiving means are divided into pulses, and for each pulse, the cross correlation between the direct wave signal and the scattered wave signal is obtained for each pulse. A pulse-by-pulse range compression means for generating a range profile of
The Doppler frequency spectrum for each block is calculated by summarizing the range profile for each pulse generated by the above-mentioned pulse-by-pulse range compressing unit into block units, and Fourier-transforming the range profile of the block unit consisting of a plurality of pulses in the hit direction. Block-by-block Doppler processing means;
A primary range migration phase compensation function for compensating for a phase change corresponding to the amount of primary range migration in the block direction corresponding to each Doppler frequency cell in the Doppler frequency spectrum calculated by the block-by-block Doppler processing means. Phase compensation function multiplication means for multiplying the Doppler frequency spectrum;
Compensation function multiplying means for multiplying the multiplication result of the phase compensation function multiplying means by a higher order range migration phase compensation function for compensating for the second or higher order range migration;
Per-pulse inverse Fourier transform means for inverse Fourier transforming the multiplication result of the compensation function multiplying means in the range direction;
A block direction Doppler processing means for calculating a Doppler frequency spectrum by Fourier transforming the inverse Fourier transform result of the inverse per-pulse Fourier transform means in the block direction;
A radar apparatus comprising: target detection means for detecting a target from the Doppler frequency spectrum calculated by the block direction Doppler processing means.
上記電波受信手段により受信された直接波の信号及び散乱波の信号をパルスに分割し、各々のパルスについて、上記直接波の信号と上記散乱波の信号との相互相関を求めることで、パルス毎のレンジプロフィールを生成するパルス毎レンジ圧縮手段と、
上記パルス毎レンジ圧縮手段により生成されたパルス毎のレンジプロフィールをブロック単位にまとめ、複数のパルスからなるブロック単位のレンジプロフィールをヒット方向にフーリエ変換することで、ブロック毎のドップラー周波数スペクトルを算出するブロック毎ドップラー処理手段と、
上記ブロック毎ドップラー処理手段により算出されたドップラー周波数スペクトルにおける各々のドップラー周波数セルに対応するブロック方向の1次のレンジマイグレーションの量に相当する位相変化を補償するための1次レンジマイグレーション位相補償関数を上記ドップラー周波数スペクトルに乗算する位相補償関数乗算手段と、
上記位相補償関数乗算手段の乗算結果をレンジ方向に逆フーリエ変換するパルス毎逆フーリエ変換手段と、
レンジマイグレーションにおける2次以上の項に影響される2次以上の位相変調を補償するための高次位相補償関数を上記パルス毎逆フーリエ変換手段の逆フーリエ変換結果に乗算する補償関数乗算手段と、
上記補償関数乗算手段の乗算結果をブロック方向にフーリエ変換することで、ドップラー周波数スペクトルを算出するブロック方向ドップラー処理手段と、
上記ブロック方向ドップラー処理手段により算出されたドップラー周波数スペクトルから目標を検出する目標検出手段と
を備えたレーダ装置。 A radio wave receiving means for receiving a direct wave of a radio wave transmitted from a radio wave source and receiving a scattered wave in which the radio wave transmitted from the radio wave source is scattered by a target;
The direct wave signal and the scattered wave signal received by the radio wave receiving means are divided into pulses, and for each pulse, the cross correlation between the direct wave signal and the scattered wave signal is obtained for each pulse. A pulse-by-pulse range compression means for generating a range profile of
The Doppler frequency spectrum for each block is calculated by summarizing the range profile for each pulse generated by the above-mentioned pulse-by-pulse range compressing unit into block units, and Fourier-transforming the range profile of the block unit consisting of a plurality of pulses in the hit direction. Block-by-block Doppler processing means;
A primary range migration phase compensation function for compensating for a phase change corresponding to the amount of primary range migration in the block direction corresponding to each Doppler frequency cell in the Doppler frequency spectrum calculated by the block-by-block Doppler processing means. Phase compensation function multiplication means for multiplying the Doppler frequency spectrum;
A pulse-by-pulse inverse Fourier transform means for inverse Fourier transforming the multiplication result of the phase compensation function multiplying means in the range direction;
Compensation function multiplication means for multiplying the inverse Fourier transform result of the inverse Fourier transform means for each pulse by a higher-order phase compensation function for compensating for the second or higher order phase modulation influenced by the second or higher order term in the range migration;
Block direction Doppler processing means for calculating a Doppler frequency spectrum by Fourier-transforming the multiplication result of the compensation function multiplication means in the block direction;
A radar apparatus comprising: target detection means for detecting a target from the Doppler frequency spectrum calculated by the block direction Doppler processing means.
上記電波受信手段により受信された直接波の信号及び散乱波の信号をパルスに分割し、各々のパルスについて、上記直接波の信号と上記散乱波の信号との相互相関を求めることで、パルス毎のレンジプロフィールを生成するパルス毎レンジ圧縮手段と、
上記パルス毎レンジ圧縮手段により生成されたパルス毎のレンジプロフィールをブロック単位にまとめ、複数のパルスからなるブロック単位のレンジプロフィールをヒット方向にフーリエ変換することで、ブロック毎のドップラー周波数スペクトルを算出するブロック毎ドップラー処理手段と、
上記ブロック毎ドップラー処理手段により算出されたドップラー周波数スペクトルにおける各々のドップラー周波数セルに対応するブロック方向の1次のレンジマイグレーションの量に相当する位相変化を補償するための1次レンジマイグレーション位相補償関数を上記ドップラー周波数スペクトルに乗算する位相補償関数乗算手段と、
上記位相補償関数乗算手段の乗算結果をレンジ方向に逆フーリエ変換するパルス毎逆フーリエ変換手段と、
上記パルス毎逆フーリエ変換手段の逆フーリエ変換結果をヒット方向に逆フーリエ変換するブロック毎逆フーリエ変換手段と、
レンジマイグレーションにおける2次以上の項に影響されるドップラーセルマイグレーションを補償するためのドップラーセルマイグレーション補償関数を上記ブロック毎逆フーリエ変換手段の逆フーリエ変換結果に乗算する第1の補償関数乗算手段と、
レンジマイグレーションにおける2次以上の項に影響される2次以上の位相変調を補償するための高次位相補償関数を上記第1の補償関数乗算手段の乗算結果に乗算する第2の補償関数乗算手段と、
上記第2の補償関数乗算手段をヒット方向にフーリエ変換するブロック毎フーリエ変換手段と、
上記ブロック毎フーリエ変換手段のフーリエ変換結果をブロック方向にフーリエ変換することで、ドップラー周波数スペクトルを算出するブロック方向ドップラー処理手段と、
上記ブロック方向ドップラー処理手段により算出されたドップラー周波数スペクトルから目標を検出する目標検出手段と
を備えたレーダ装置。 A radio wave receiving means for receiving a direct wave of a radio wave transmitted from a radio wave source and receiving a scattered wave in which the radio wave transmitted from the radio wave source is scattered by a target;
The direct wave signal and the scattered wave signal received by the radio wave receiving means are divided into pulses, and for each pulse, the cross correlation between the direct wave signal and the scattered wave signal is obtained for each pulse. A pulse-by-pulse range compression means for generating a range profile of
The Doppler frequency spectrum for each block is calculated by summarizing the range profile for each pulse generated by the above-mentioned pulse-by-pulse range compressing unit into block units, and Fourier-transforming the range profile of the block unit consisting of a plurality of pulses in the hit direction. Block-by-block Doppler processing means;
A primary range migration phase compensation function for compensating for a phase change corresponding to the amount of primary range migration in the block direction corresponding to each Doppler frequency cell in the Doppler frequency spectrum calculated by the block-by-block Doppler processing means. Phase compensation function multiplication means for multiplying the Doppler frequency spectrum;
A pulse-by-pulse inverse Fourier transform means for inverse Fourier transforming the multiplication result of the phase compensation function multiplying means in the range direction;
Block inverse Fourier transform means for inverse Fourier transforming the inverse Fourier transform result of the pulse-by-pulse inverse Fourier transform means in the hit direction;
First compensation function multiplying means for multiplying the inverse Fourier transform result of the block-by-block inverse Fourier transform means by a Doppler cell migration compensation function for compensating for Doppler cell migration affected by a second-order or higher term in range migration;
Second compensation function multiplying means for multiplying the multiplication result of the first compensation function multiplying means by a higher order phase compensation function for compensating for the second or higher order phase modulation affected by the second or higher order term in the range migration. When,
Block-by-block Fourier transform means for Fourier transforming the second compensation function multiplication means in the hit direction;
Block direction Doppler processing means for calculating a Doppler frequency spectrum by Fourier transforming the Fourier transform result of the block-by-block Fourier transform means in the block direction; and
A radar apparatus comprising: target detection means for detecting a target from the Doppler frequency spectrum calculated by the block direction Doppler processing means.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2012157392A JP5865794B2 (en) | 2012-07-13 | 2012-07-13 | Radar equipment |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2012157392A JP5865794B2 (en) | 2012-07-13 | 2012-07-13 | Radar equipment |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2014020820A JP2014020820A (en) | 2014-02-03 |
JP5865794B2 true JP5865794B2 (en) | 2016-02-17 |
Family
ID=50195889
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2012157392A Expired - Fee Related JP5865794B2 (en) | 2012-07-13 | 2012-07-13 | Radar equipment |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP5865794B2 (en) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5823062B2 (en) * | 2012-12-14 | 2015-11-25 | 三菱電機株式会社 | Radar equipment |
JP6081394B2 (en) * | 2014-02-28 | 2017-02-15 | 株式会社東芝 | Positioning system and positioning method |
CN103901419B (en) * | 2014-04-22 | 2016-03-30 | 北京理工大学 | A kind of external illuminators-based radar range migration compensation method based on frequency domain phase correction |
CN104062640B (en) * | 2014-06-30 | 2016-03-30 | 北京理工大学 | The Fast implementation that a kind of external illuminators-based radar range migration compensates |
RU2583537C1 (en) * | 2015-03-26 | 2016-05-10 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Рязанский государственный радиотехнический университет" | Auto-compensator for doppler phase of passive interference |
RU2624795C1 (en) * | 2016-05-10 | 2017-07-06 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Рязанский государственный радиотехнический университет" | Autocompensor of doppler shifts of phase of interference |
RU2642808C1 (en) * | 2016-12-20 | 2018-01-26 | Федеральное государственное бюджетное образовательное учреждение высшего образования "Рязанский государственный радиотехнический университет" | Interference suppressor |
WO2023089713A1 (en) | 2021-11-18 | 2023-05-25 | 三菱電機株式会社 | Radar signal processing device, radar signal processing method, and target observation system |
Family Cites Families (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0664142B2 (en) * | 1985-08-07 | 1994-08-22 | 株式会社東芝 | Radar device |
US5235338A (en) * | 1990-10-31 | 1993-08-10 | Hsiao Stephen S | Moving target detection through range cell migration radar |
DE4122592C1 (en) * | 1991-07-08 | 1993-01-28 | Deutsche Forschungsanstalt Fuer Luft- Und Raumfahrt Ev, 5300 Bonn, De | |
US5200754A (en) * | 1991-11-27 | 1993-04-06 | Hughes Aircraft Company | Fourth-order-product phase difference autofocus |
DE4427657C2 (en) * | 1994-08-05 | 1996-10-24 | Deutsche Forsch Luft Raumfahrt | Process for image generation by means of two-dimensional data processing on a radar with a synthetic aperture |
JPH08179037A (en) * | 1994-12-22 | 1996-07-12 | Mitsubishi Electric Corp | Radar device |
JP4046422B2 (en) * | 1998-09-10 | 2008-02-13 | 三菱電機株式会社 | Motion compensation circuit and radar device |
JP4597444B2 (en) * | 2001-09-19 | 2010-12-15 | 三菱電機株式会社 | Synthetic aperture radar apparatus and image reproduction method in synthetic aperture radar apparatus |
JP4110896B2 (en) * | 2002-09-10 | 2008-07-02 | 三菱電機株式会社 | Radar equipment |
JP4377772B2 (en) * | 2004-08-11 | 2009-12-02 | 株式会社東芝 | Radar equipment |
JP2006071363A (en) * | 2004-08-31 | 2006-03-16 | Toshiba Corp | Radar installation |
JP4602179B2 (en) * | 2005-02-15 | 2010-12-22 | 三菱電機株式会社 | Radar equipment |
JP4943065B2 (en) * | 2006-06-15 | 2012-05-30 | 三菱電機株式会社 | Image radar device |
JP5197138B2 (en) * | 2008-04-30 | 2013-05-15 | 三菱電機株式会社 | Multi static radar device |
JP2010019623A (en) * | 2008-07-09 | 2010-01-28 | Toshiba Corp | Signal processing circuit |
JP5606097B2 (en) * | 2010-02-25 | 2014-10-15 | 三菱電機株式会社 | Passive radar device |
JP2011247593A (en) * | 2010-05-21 | 2011-12-08 | Mitsubishi Electric Corp | Image radar device |
EP2677343B1 (en) * | 2011-02-18 | 2021-03-24 | Mitsubishi Electric Corporation | Passive radar device |
JP5679904B2 (en) * | 2011-05-18 | 2015-03-04 | 三菱電機株式会社 | Image radar device |
JP5979868B2 (en) * | 2011-12-21 | 2016-08-31 | 三菱電機株式会社 | Image radar device |
-
2012
- 2012-07-13 JP JP2012157392A patent/JP5865794B2/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JP2014020820A (en) | 2014-02-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5865794B2 (en) | Radar equipment | |
JP5665893B2 (en) | Passive radar device | |
JP5823062B2 (en) | Radar equipment | |
US7375676B1 (en) | Mutual interference processor for pulse radar | |
JP5606097B2 (en) | Passive radar device | |
Tian et al. | Multiband radar signal coherent fusion processing with IAA and apFFT | |
KR101703773B1 (en) | Method for acquiring radar image using serveillance radar and System thereof | |
JP5251591B2 (en) | Pulse compressor | |
US8760340B2 (en) | Processing radar return signals to detect targets | |
CN102520396A (en) | Time-varying narrow-band interference suppression method based on complex empirical mode decomposition | |
JP2009257907A (en) | Target detector | |
CN110568434A (en) | multi-channel uniform acceleration trajectory SAR moving target two-dimensional speed estimation method | |
US10859691B2 (en) | Radar range accuracy improvement method | |
EP3399334B1 (en) | Object detecting device and sensor device | |
Wan et al. | A Doppler ambiguity tolerated method for radar sensor maneuvering target focusing and detection | |
US20170269194A1 (en) | Signal processing device, radar apparatus and signal processing method | |
JP4793240B2 (en) | High resolution radar equipment | |
CN107976663A (en) | It is a kind of based on the external illuminators-based radar of subspace projection to targeted compression detection method | |
Liu et al. | Novel range profile synthesis algorithm for linearly stepped-frequency modulated inversed synthetic aperture radar imaging of remote manoeuvring target | |
JP4982229B2 (en) | Moving target detection device | |
JP6164936B2 (en) | Radar equipment | |
KR101007389B1 (en) | Apparatus and method for generating image | |
JP5625326B2 (en) | Radar apparatus and distance measurement method thereof | |
JP2016038267A (en) | Delay time estimation device and range-height device | |
CN115932921B (en) | Multi-satellite multi-carrier non-coherent combined accumulation processing method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20141204 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20150907 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20150915 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20151113 |
|
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: 20151201 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20151228 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 5865794 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
LAPS | Cancellation because of no payment of annual fees |