JP5863443B2 - Radar equipment - Google Patents
Radar equipment Download PDFInfo
- Publication number
- JP5863443B2 JP5863443B2 JP2011281614A JP2011281614A JP5863443B2 JP 5863443 B2 JP5863443 B2 JP 5863443B2 JP 2011281614 A JP2011281614 A JP 2011281614A JP 2011281614 A JP2011281614 A JP 2011281614A JP 5863443 B2 JP5863443 B2 JP 5863443B2
- Authority
- JP
- Japan
- Prior art keywords
- pulse
- frequency spectrum
- pulse compression
- calculation unit
- calculated
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000007906 compression Methods 0.000 claims description 178
- 230000006835 compression Effects 0.000 claims description 176
- 239000011159 matrix material Substances 0.000 claims description 175
- 238000001228 spectrum Methods 0.000 claims description 107
- 239000013598 vector Substances 0.000 claims description 65
- 230000003595 spectral effect Effects 0.000 claims description 23
- 230000010354 integration Effects 0.000 claims description 21
- 230000001427 coherent effect Effects 0.000 claims description 17
- 230000008054 signal transmission Effects 0.000 claims description 14
- 238000001514 detection method Methods 0.000 claims description 13
- 230000035485 pulse pressure Effects 0.000 claims 2
- 230000002194 synthesizing effect Effects 0.000 claims 1
- 230000005540 biological transmission Effects 0.000 description 60
- 238000000034 method Methods 0.000 description 17
- 238000010586 diagram Methods 0.000 description 16
- 230000006870 function Effects 0.000 description 12
- 230000000694 effects Effects 0.000 description 6
- 238000006243 chemical reaction Methods 0.000 description 3
- 238000010606 normalization Methods 0.000 description 3
- 239000004065 semiconductor Substances 0.000 description 3
- 230000007850 degeneration Effects 0.000 description 2
- 230000006866 deterioration Effects 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000001174 ascending effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 230000000593 degrading effect Effects 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Description
この発明は、受信信号をパルス圧縮して、レンジサイドローブのレベルを低減するレーダ装置に関するものである。 The present invention relates to a radar apparatus that reduces the level of a range side lobe by pulse-compressing a received signal.
レーダ装置では、周波数変調や位相変調などを施しているパルス信号を空中に放射し、目標又はクラッタに反射されたパルス信号の散乱波を受信する際、送信したパルス信号と相関が高いレプリカを用いた相関処理を行う受信信号のパルス圧縮が広く実用化されている。
受信信号のパルス圧縮を行うことで、送信パルス幅の増長に伴う送信エネルギー増大と、相関処理による圧縮後のパルス幅の狭小化が可能になり、レーダ装置の目標探知距離と距離分解能を同時に高めることができる。
In radar equipment, a pulse signal that has been subjected to frequency modulation or phase modulation is emitted into the air, and when a scattered wave of the pulse signal reflected by the target or clutter is received, a replica that has a high correlation with the transmitted pulse signal is used. The pulse compression of received signals for performing correlation processing has been widely put into practical use.
By performing pulse compression of the received signal, it is possible to increase the transmission energy accompanying the increase of the transmission pulse width and to narrow the pulse width after compression by correlation processing, and simultaneously increase the target detection distance and distance resolution of the radar device. be able to.
パルス圧縮後の出力波形は、目標等の距離に対応する遅延時間を中心とするメインローブと、そのメインローブの前後の送信パルス幅相当の範囲に分布するレンジサイドローブとから構成される。
このとき、レンジサイドローブ領域に存在する別目標のSIR(Signal to Interference Ratio)を劣化させる問題がある。
The output waveform after the pulse compression is composed of a main lobe centered on a delay time corresponding to a distance such as a target, and range side lobes distributed in a range corresponding to the transmission pulse width before and after the main lobe.
At this time, there is a problem of degrading another target SIR (Signal to Interference Ratio) existing in the range side lobe region.
この問題を解決する方法として、窓関数を利用するパルス圧縮によってレンジサイドローブを低減する方法がある。
例えば、Taylor窓関数を用いるパルス圧縮を行うことで、レンジサイドローブのピークレベルを−60dBまで低減することが可能である。窓関数を利用しないでパルス圧縮を行う場合のピークレベルが−13.2dBであるため、ピークレベルを46.8dB改善できることになる。
As a method of solving this problem, there is a method of reducing range side lobes by pulse compression using a window function.
For example, it is possible to reduce the peak level of the range side lobe to −60 dB by performing pulse compression using a Taylor window function. Since the peak level when performing pulse compression without using a window function is −13.2 dB, the peak level can be improved by 46.8 dB.
しかし、窓関数を乗じることで、距離分解能を劣化させる圧縮後パルス幅の拡大や、探知距離を縮退させるミスマッチ損失が原理的に発生することが以下の非特許文献1に開示されている。
レンジサイドローブを低減する特性と、距離分解能の劣化や探知距離の縮退させる特性とは、トレードオフの関係にあるため、レーダ装置を設計する際、レーダの要求性能を満たす条件で最も適する窓関数を選択する。
However, the following non-patent document 1 discloses that by multiplying by a window function, in principle, a post-compression pulse width that degrades the distance resolution and a mismatch loss that degenerates the detection distance occur.
There is a trade-off between the characteristics that reduce the range side lobe and the characteristics that degrade the range resolution and reduce the detection distance. Therefore, when designing a radar system, the window function that is most suitable under the conditions that satisfy the required performance of the radar. Select.
従来のレーダ装置は以上のように構成されているので、パルス圧縮によってレンジサイドローブの低減を図ると、圧縮後パルス幅の拡大やミスマッチ損失が発生して、距離分解能の劣化や探知距離の縮退を招いてしまうなどの課題があった。 Since the conventional radar device is configured as described above, if the range side lobe is reduced by pulse compression, the pulse width after compression and mismatch loss occur, resulting in degradation of distance resolution and reduction of detection distance. There was a problem such as inviting.
この発明は上記のような課題を解決するためになされたもので、距離分解能の劣化や探知距離の縮退を招くことなく、低レンジサイドローブ特性を実現することができるレーダ装置を得ることを目的とする。 The present invention has been made to solve the above-described problems, and an object of the present invention is to obtain a radar apparatus capable of realizing a low-range sidelobe characteristic without causing deterioration in distance resolution and degeneration of detection distance. And
この発明に係るレーダ装置は、パルス信号を空中に放射するパルス信号送信手段と、パルス信号送信手段から空中に放射されたのち、目標又はクラッタに反射された上記パルス信号の散乱波を受信して、その散乱波の受信信号を出力する散乱波受信手段と、パルス信号送信手段により放射されたパルス信号の周波数スペクトルから散乱波ステアリングベクトルを算出し、パルス圧縮に伴うピークレベルの低下を防ぐガード対象のレンジビンの前後のレンジサイドローブ領域における散乱波ステアリングベクトルから構成されるアレーマニフォルドを定義して、そのアレーマニフォルドからサイドローブ干渉波の相関行列を算出し、その相関行列を用いて、散乱波受信手段から出力された受信信号をパルス圧縮するパルス圧縮手段とを設け、目標検出手段が、パルス圧縮手段によりパルス圧縮された受信信号から目標を検出するようにしたものである。 A radar apparatus according to the present invention receives a pulse signal transmission unit that radiates a pulse signal in the air, and receives a scattered wave of the pulse signal reflected from a target or a clutter after being radiated from the pulse signal transmission unit to the air. A scattered wave receiving means for outputting the received signal of the scattered wave, and a guard object for calculating a scattered wave steering vector from the frequency spectrum of the pulse signal radiated by the pulse signal transmitting means to prevent a decrease in peak level due to pulse compression. Define an array manifold composed of scattered wave steering vectors in the range sidelobe region before and after the range bin, calculate the correlation matrix of the sidelobe interference wave from the array manifold, and use the correlation matrix to receive the scattered wave And a pulse compression means for pulse-compressing the received signal output from the means. Detecting means, in which to detect the target from a received signal pulse compression by the pulse compression means.
この発明によれば、パルス信号送信手段により放射されたパルス信号の周波数スペクトルから散乱波ステアリングベクトルを算出し、パルス圧縮に伴うピークレベルの低下を防ぐガード対象のレンジビンの前後のレンジサイドローブ領域における散乱波ステアリングベクトルから構成されるアレーマニフォルドを定義して、そのアレーマニフォルドからサイドローブ干渉波の相関行列を算出し、その相関行列を用いて、散乱波受信手段から出力された受信信号をパルス圧縮するように構成したので、距離分解能の劣化や探知距離の縮退を招くことなく、低レンジサイドローブ特性を実現することができる効果がある。 According to the present invention, the scattered wave steering vector is calculated from the frequency spectrum of the pulse signal radiated by the pulse signal transmitting means, and in the range sidelobe region before and after the guard target range bin to prevent the peak level from being lowered due to the pulse compression. Define an array manifold composed of scattered wave steering vectors, calculate the correlation matrix of the sidelobe interference wave from the array manifold, and use this correlation matrix to pulse-compress the received signal output from the scattered wave receiver Thus, the low-range side lobe characteristic can be realized without causing deterioration of the distance resolution and degeneration of the detection distance.
実施の形態1.
図1はこの発明の実施の形態1によるレーダ装置を示す構成図である。
図1において、送信機1は所定の送信周波数信号を生成し、所定のPRI(Pulse Repetition Interval)間隔で、その送信周波数信号を用いて、中心周波数fc,変調帯域幅B,パルス幅τucの送信パルスr(t)を生成し、その送信パルスr(t)をデュプレクサ2に出力する処理を実施する。
デュプレクサ2は送信機1から出力された送信パルスr(t)をアンテナ3に出力する一方、アンテナ3により受信された散乱波(アンテナ3から空中に放射されたのち、目標又はクラッタに反射された送信パルスr(t)の散乱波)を受信機4に出力する切替器である。
Embodiment 1 FIG.
1 is a block diagram showing a radar apparatus according to Embodiment 1 of the present invention.
In FIG. 1, a transmitter 1 generates a predetermined transmission frequency signal, and uses the transmission frequency signal at a predetermined PRI (Pulse Repetition Interval) interval to use a center frequency f c , a modulation bandwidth B, and a pulse width τ uc. The transmission pulse r (t) is generated and the transmission pulse r (t) is output to the
The
アンテナ3はデュプレクサ2から出力された送信パルスr(t)を所定のビーム指向方向に向けて空中に放射する一方、目標又はクラッタに反射された送信パルスr(t)の散乱波を受信する。
なお、送信機1、デュプレクサ2及びアンテナ3からパルス信号送信手段が構成されている。
The antenna 3 radiates the transmission pulse r (t) output from the
The transmitter 1, the
受信機4はアンテナ3により受信された散乱波をデュプレクサ2経由で受けると、その散乱波を示す信号の周波数をベースバンド帯に変換して、ベースバンド帯のアナログ受信信号を生成する処理を実施する。
AD変換器5は受信機4により生成されたアナログ受信信号をAD変換して、ディジタルの受信信号z(t)をパルス圧縮器6に出力する処理を実施する。
なお、デュプレクサ2、アンテナ3、受信機4及びAD変換器5から散乱波受信手段が構成されている。
When the receiver 4 receives the scattered wave received by the antenna 3 via the
The
The
パルス圧縮器6はPRI周期毎に、AD変換器5から出力された受信信号z(t)をパルス圧縮する処理器である。
即ち、パルス圧縮器6は送信パルスr(t)の周波数スペクトルR(f)から、パルス圧縮に伴うピークレベルの低下を防ぐガード対象のレンジビンの前後のレンジサイドローブ領域における散乱波ステアリングベクトルから構成されるアレーマニフォルドAi (AM)を算出して、そのアレーマニフォルドAi (AM)からサイドローブ干渉波の相関行列Riを算出し、その相関行列Riの逆行列Ri −1を用いて、AD変換器5から出力された受信信号z(t)をパルス圧縮する処理を実施する。なお、パルス圧縮器6はパルス圧縮手段を構成している。
目標検出器7はパルス圧縮器6によりパルス圧縮された受信信号z’(t)から目標を検出する処理を実施する。なお、目標検出器7は目標検出手段を構成している。
The
That is, the
The
図2はこの発明の実施の形態1によるレーダ装置のパルス圧縮器6を示す構成図である。
図2において、第1の周波数スペクトル算出部である離散フーリエ変換部11はAD変換器5から出力された受信信号z(t)を離散フーリエ変換して、その受信信号z(t)の周波数スペクトルZ(f)を算出する処理を実施する。
第2の周波数スペクトル算出部である離散フーリエ変換部12は送信機1から出力された送信パルスr(t)をレプリカとして、そのレプリカを離散フーリエ変換して、そのレプリカの周波数スペクトルR(f)を算出する処理を実施する。
FIG. 2 is a block diagram showing a
In FIG. 2, a discrete Fourier
The discrete Fourier
スペクトル積算出部13は離散フーリエ変換部12により算出されたレプリカの周波数スペクトルR(f)の複素共役R*(f)を離散フーリエ変換部11により算出された受信信号z(t)の周波数スペクトルZ(f)に乗算してスペクトル積X(f)を算出する処理を実施する。
アレーマニフォルド定義部14は離散フーリエ変換部12により算出されたレプリカの周波数スペクトルR(f)から散乱波ステアリングベクトルを算出して、パルス圧縮に伴うピークレベルの低下を防ぐガード対象のレンジビン(ガード領域)の前後のレンジサイドローブ領域における散乱波ステアリングベクトルから構成されるアレーマニフォルドAi (AM)を定義する処理を実施する。
The spectrum
The array
相関行列算出部15はアレーマニフォルド定義部14により定義されたアレーマニフォルドAi (AM)と対角荷重であるDLL(Diagnal Loading Level)からサイドローブ干渉波の相関行列Riを算出する処理を実施する。
パルス圧縮ウェイト算出部16は相関行列算出部15により算出された相関行列Riの逆行列Ri −1を用いて、各レンジビンに対応するパルス圧縮ウェイトw(pro)を算出する処理を実施する。
コヒーレント積分処理部17はパルス圧縮ウェイト算出部16により算出されたパルス圧縮ウェイトw(pro)を用いて、スペクトル積算出部13により算出されたスペクトル積X(f)をコヒーレント積分し、その積分結果をパルス圧縮後の受信信号z’(t)として目標検出器7に出力する処理を実施する。
The correlation
The pulse compression
The coherent
図2の例では、パルス圧縮器6の構成要素である離散フーリエ変換部11,12、スペクトル積算出部13、アレーマニフォルド定義部14、相関行列算出部15、パルス圧縮ウェイト算出部16及びコヒーレント積分処理部17のそれぞれが専用のハードウェア(例えば、CPUを実装している半導体集積回路、あるいは、ワンチップマイコン)で構成されているものを想定しているが、パルス圧縮器6がコンピュータで構成されていてもよい。
パルス圧縮器6がコンピュータで構成されている場合、離散フーリエ変換部11,12、スペクトル積算出部13、アレーマニフォルド定義部14、相関行列算出部15、パルス圧縮ウェイト算出部16及びコヒーレント積分処理部17の処理内容を記述しているプログラムをコンピュータのメモリに格納し、当該コンピュータのCPUが当該メモリに格納されているプログラムを実行するようにすればよい。
In the example of FIG. 2, discrete Fourier
When the
次に動作について説明する。
送信機1は、所定の送信周波数信号を生成する。
そして、送信機1は、所定のPRI間隔で、その送信周波数信号を用いて、中心周波数fc,変調帯域幅B,パルス幅τucの送信パルスr(t)を生成し、その送信パルスr(t)をデュプレクサ2に出力する。
デュプレクサ2は、送信機1から送信パルスr(t)を受けると、その送信パルスr(t)をアンテナ3に出力する。
これにより、アンテナ3から送信パルスr(t)が所定のビーム指向方向に向けて空中に放射される。
その後、目標又はクラッタに反射された送信パルスr(t)の散乱波がアンテナ3に到来し、アンテナ3により送信パルスr(t)の散乱波が受信される。
Next, the operation will be described.
The transmitter 1 generates a predetermined transmission frequency signal.
Then, the transmitter 1 generates a transmission pulse r (t) having a center frequency f c , a modulation bandwidth B, and a pulse width τ uc using the transmission frequency signal at a predetermined PRI interval, and the transmission pulse r (T) is output to the
When receiving the transmission pulse r (t) from the transmitter 1, the
As a result, the transmission pulse r (t) is radiated from the antenna 3 toward the predetermined beam directing direction.
Thereafter, the scattered wave of the transmission pulse r (t) reflected by the target or the clutter arrives at the antenna 3, and the scattered wave of the transmission pulse r (t) is received by the antenna 3.
デュプレクサ2は、アンテナ3により受信された散乱波を受信機4に出力する。
受信機4は、アンテナ3により受信された散乱波をデュプレクサ2経由で受けると、その散乱波を示す信号の周波数をベースバンド帯に変換して、ベースバンド帯のアナログ受信信号を生成する。
AD変換器5は、受信機4がアナログ受信信号を生成すると、そのアナログ受信信号をAD変換して、ディジタルの受信信号z(t)をパルス圧縮器6に出力する。
The
When the receiver 4 receives the scattered wave received by the antenna 3 via the
When the receiver 4 generates an analog reception signal, the
パルス圧縮器6は、PRI周期毎に、AD変換器5から出力された受信信号z(t)をパルス圧縮し、パルス圧縮後の受信信号z’(t)を目標検出器7に出力する。
即ち、パルス圧縮器6は、送信パルスr(t)の周波数スペクトルR(f)から、パルス圧縮に伴うピークレベルの低下を防ぐガード対象のレンジビンの前後のレンジサイドローブ領域における散乱波ステアリングベクトルから構成されるアレーマニフォルドAi (AM)を算出して、そのアレーマニフォルドAi (AM)からサイドローブ干渉波の相関行列Riを算出し、その相関行列Riの逆行列Ri −1を用いて、AD変換器5から出力された受信信号z(t)をパルス圧縮する。
The
That is, the
目標検出器7は、パルス圧縮器6からパルス圧縮後の受信信号z’(t)を受けると、その受信信号z’(t)から目標を検出する。
目標検出器7における目標の検出処理は公知の技術であるため、ここでは詳細な説明を省略する。
When the
Since the target detection process in the
以下、パルス圧縮器6の処理内容を具体的に説明する。
ただし、以下の説明では、説明の簡単化のため、送信と受信の切替タイミング付近で発生する受信パルスのエクリプスについては考えない。
まず、AD変換器5からPRI間隔で出力される受信信号z(t)は、K個の散乱波を含んでおり、下記の式(1)のように表される。
式(1)において、tkは第k番目の散乱波の遅延時間、fk (D)は第k番目の散乱波のドップラ周波数である。
また、αkは第k番目の散乱波の距離等による減衰係数、n(t)は受信機雑音である。
Hereinafter, the processing content of the
However, in the following description, for simplification of explanation, Eclipse of received pulses generated near the transmission and reception switching timing is not considered.
First, the received signal z (t) output from the
In the formula (1), t k is the delay time of the k-th scattered wave, f k (D) is the Doppler frequency of the k-th scattered wave.
Α k is an attenuation coefficient due to the distance of the k-th scattered wave, and n (t) is receiver noise.
パルス圧縮器6の離散フーリエ変換部11は、AD変換器5から受信信号z(t)を受けると、その受信信号z(t)を離散フーリエ変換して、その受信信号z(t)の周波数スペクトルZ(f)を算出する。受信信号z(t)の周波数スペクトルZ(f)は、下記の式(2)のように表される。
式(2)において、R(f)は送信パルスr(t)の周波数スペクトル(離散フーリエ変換部12により算出される周波数スペクトル)、N(f)は受信機雑音の周波数スペクトルである。
When receiving the received signal z (t) from the
In Equation (2), R (f) is the frequency spectrum of the transmission pulse r (t) (frequency spectrum calculated by the discrete Fourier transform unit 12), and N (f) is the frequency spectrum of the receiver noise.
離散フーリエ変換部12は、送信機1から出力された送信パルスr(t)をレプリカとし、そのレプリカを離散フーリエ変換して、そのレプリカの周波数スペクトルR(f)を算出する。
The discrete
スペクトル積算出部13は、離散フーリエ変換部11が受信信号z(t)の周波数スペクトルZ(f)を算出し、離散フーリエ変換部12がレプリカの周波数スペクトルR(f)を算出すると、下記の式(3)に示すように、そのレプリカの周波数スペクトルR(f)の複素共役R*(f)を受信信号z(t)の周波数スペクトルZ(f)に乗算して、スペクトル積X(f)を算出する。
When the discrete
ここでは、周波数スペクトルZ(f),R(f)がM点のサンプルより構成されているものとして、下記の式(4)のような受信信号ベクトルxを定義する。
なお、Mはレンジビン数でもあり、z(t),r(t)がM点のサンプル点で表されることに相当する。
ただし、第m番目の周波数fmは、下記の式(5)のように与えるものとする。
また、Δfはリファレンス信号の周波数領域でのサンプリング間隔であり、B≦MΔfを満たすものとする。
Here, it is assumed that the frequency spectrums Z (f) and R (f) are composed of M points of samples, and a received signal vector x as shown in the following equation (4) is defined.
Note that M is also the number of range bins, which corresponds to z (t) and r (t) being represented by M sample points.
However, it is assumed that the mth frequency fm is given by the following equation (5).
Δf is a sampling interval in the frequency domain of the reference signal, and B ≦ MΔf is satisfied.
ここで、式(4)に式(3)のスペクトル積X(f)を代入すると、以下のように変形することができる。
式(6)において、akは第k番目の目標に対する遅延時間tk及びドップラ周波数fk (D)から決まるステアリングベクトルである。
また、AはK個のステアリングベクトルを列ベクトルとする行列、sは振幅ベクトル、nは受信機雑音ベクトルである。
Here, when the spectral product X (f) of the equation (3) is substituted into the equation (4), the following modification can be made.
In Equation (6), a k is a steering vector determined from the delay time t k and the Doppler frequency f k (D) for the k-th target.
A is a matrix having K steering vectors as column vectors, s is an amplitude vector, and n is a receiver noise vector.
アレーマニフォルド定義部14は、離散フーリエ変換部12がレプリカの周波数スペクトルR(f)を算出すると、そのレプリカの周波数スペクトルR(f)から散乱波ステアリングベクトルを算出して、パルス圧縮に伴うピークレベルの低下を防ぐガード対象のレンジビン(ガード領域)の前後のレンジサイドローブ領域における散乱波ステアリングベクトルから構成されるアレーマニフォルドAi (AM)を定義する。
具体的には、以下の通りである。
When the discrete
Specifically, it is as follows.
ここでは、PRI内における遅延時間ti (AM)に対応する散乱波の受信信号ベクトルci (AM)を考える。
ただし、ci (AM)はパルス圧縮ウェイトの計算に用いるものであり、遅延時間ti (AM)は任意ではなく、PRI内を時間間隔dtで区切った場合の第i番目のサンプル点に対応するものとする。
このとき、遅延時間ti (AM)とレンジビン数Mから決まるtmとの間には、以下の関係が成り立つものとする。なお、ΔtはAD変換器5のサンプリング時間間隔である。
Here, the received signal vector c i (AM) of the scattered wave corresponding to the delay time t i (AM) in the PRI is considered.
However, c i (AM) is used for calculating the pulse compression weight, and the delay time t i (AM) is not arbitrary, and corresponds to the i-th sample point when the PRI is divided by the time interval dt. It shall be.
At this time, it is assumed that the following relationship holds between the delay time t i (AM) and t m determined from the number of range bins M. Δt is a sampling time interval of the
例えば、散乱波の強度を1とすると、受信信号ベクトルci (AM)は散乱波のステアリングベクトルai (AM)と等しくなり、下記の式(13)のように与えられる。
ただし、式(13)のDiは遅延時間ti (AM)に関する位相回転行列であり、下記の式(14)のように与えられる。
また、a1 (AM)は下記の式(15)のように与えられる。
For example, if the intensity of the scattered wave is 1, the received signal vector c i (AM) is equal to the steering vector a i (AM) of the scattered wave, and is given by the following equation (13).
However, D i in the equation (13) is a phase rotation matrix with respect to the delay time t i (AM) , and is given by the following equation (14).
Further, a 1 (AM) is given by the following equation (15).
遅延時間ti (AM)の散乱点に対応するパルス圧縮出力は、遅延時間ti (AM)に現れるメインローブを中心とすると、ti (AM)−τucからti (AM)+τucの範囲にサイドローブを伴う波形となる。
このとき、同範囲に2L+1個の散乱点(ただし、強度は1とする)が時間間隔dtで分布するものとする。
ただし、Lは遅延時間ti (AM)の前側あるいは後側の送信パルス幅τucに相当する時間領域に分布する散乱波点数である。
遅延時間ti (AM)の前側あるいは後側のti−l (AM),ti+l (AM)(ただし、l=1,…,L)は、下記の式(16)(17)で表される。
The pulse compression output corresponding to the scattering point of the delay time t i (AM), when the center of the main lobe appearing in the delay time t i (AM), from t i (AM) -τ uc t i (AM) + τ uc A waveform with side lobes in the range.
At this time, it is assumed that 2L + 1 scattering points (however, the intensity is 1) are distributed at the time interval dt in the same range.
Here, L is the number of scattered wave points distributed in the time region corresponding to the transmission pulse width τ uc on the front side or the rear side of the delay time t i (AM) .
T i−l (AM) , t i + l (AM) on the front side or rear side of the delay time t i (AM) (where l = 1,..., L) is expressed by the following equations (16) and (17). Is done.
遅延時間ti−l (AM),ti+l (AM)に対応する散乱波の受信信号ベクトルci−l (AM),ci+l (AM)は、式(13)を利用すると、下記の式(18)(19)のように与えられる。
The received signal vectors c i−l (AM) and c i + l (AM) of the scattered wave corresponding to the delay times t i−l (AM) and t i + l (AM) are expressed by the following equation using the equation (13): (18) It is given as (19).
次に、2L+1個の散乱点のステアリングベクトルを並べたアレーマニフォルドを考える。
さらに、散乱波のステアリングベクトルai (AM)と、その前後のLG個のステアリングベクトルをゼロベクトル0にするため、下記の式(21)に示すようなガード行列Gを定義する。ただし、l’=1,…,2L+1である。
ゼロベクトル0となる範囲を「ガード領域」と称する。
Next, consider an array manifold in which steering vectors of 2L + 1 scattering points are arranged.
Further, a steering vector a i of the scattered wave (AM), for the L G number of steering vectors before and after the zero
A range where the zero vector is 0 is referred to as a “guard region”.
アレーマニフォルド定義部14は、式(21)のガード行列Gによって、パルス圧縮のメインローブ付近の遅延時間に対応するステアリングベクトルをゼロベクトルにするものとし、下記の式(22)に示すように、主にサイドローブ領域の散乱点ステアリングベクトルから構成されるアレーマニフォルドAi (AM)を定義する。
The array
相関行列算出部15は、アレーマニフォルド定義部14がアレーマニフォルドAi (AM)を定義すると、そのアレーマニフォルドAi (AM)とDLLからサイドローブ干渉波の相関行列Riを算出する。
即ち、相関行列算出部15は、下記の式(23)に示すように、遅延時間ti (AM)に関するパルス圧縮におけるサイドローブ干渉波の相関行列Riを定義する。
式(23)において、ρはDLL、ρIは「Diagonal Loading」のための行列、Diは遅延時間ti (AM)に関する位相回転行列である。
When the array
That is, the correlation
In Expression (23), ρ is a DLL, ρI is a matrix for “Diagonal Loading”, and D i is a phase rotation matrix related to the delay time t i (AM) .
式(23)より、遅延時間ti (AM)に関する相関行列Riは、遅延時間t1 (AM)に関する相関行列R1を求めておけば、遅延時間ti (AM)に関する位相回転行列Diを利用して簡単に求められることがわかる。
ここで、サイドローブの干渉波を抑圧するための行列を相関行列Riの逆行列Ri −1とすると、相関行列Riの逆行列Ri −1は、下記の式(24)のように与えられる。
From equation (23), the correlation matrix R i relating to a delay time t i (AM) is, if the correlation matrix R 1 concerning the delay time t 1 (AM), phase rotation matrix about a delay time t i (AM) D It can be seen that i can be easily obtained.
Here, when the matrix for suppressing interference waves sidelobe and inverse matrix R i -1 of the correlation matrix R i, the inverse matrix R i -1 of the correlation matrix R i is as the following equation (24) Given to.
パルス圧縮ウェイト算出部16は、相関行列算出部15がサイドローブ干渉波の相関行列Riを算出すると、その相関行列Riの逆行列Ri −1を用いて、各レンジビンに対応するパルス圧縮ウェイトw(pro)を算出する。
ここで、遅延時間ti (AM)に関する相関行列Riの逆行列Ri −1は、式(24)より、遅延時間ti (AM)から直接求めずに、遅延時間t1 (AM)に関する逆行列R1 −1を求めておき、遅延時間ti (AM)に関する位相回転行列Diを逆行列R1 −1の両側から乗じることで求められることがわかる。
When the correlation
Here, the inverse matrix R i -1 of the correlation matrix R i relating to a delay time t i (AM), from equation (24), without requiring direct from the delay time t i (AM), the delay time t 1 (AM) It can be seen that the inverse matrix R 1 −1 is obtained and the phase rotation matrix D i related to the delay time t i (AM) is multiplied from both sides of the inverse matrix R 1 −1 .
そのため、パルス圧縮ウェイト算出部16は、第m番目のレンジビンに対応するパルス圧縮ウェイトwm (pro)を下記の式(25)のように与える。
ただし、αmは規格化係数、wmは第m番目のレンジビンに対応する時間遅延ベクトルである。
Therefore, the pulse compression
Here, α m is a normalization coefficient, and w m is a time delay vector corresponding to the m-th range bin.
コヒーレント積分処理部17は、パルス圧縮ウェイト算出部16が各レンジビンに対応するパルス圧縮ウェイトw(pro)を算出すると、そのパルス圧縮ウェイトw(pro)を用いて、スペクトル積算出部13により算出されたスペクトル積X(f)をコヒーレント積分し、その積分結果をパルス圧縮後の受信信号z’(t)として、目標検出器7に出力する。
When the pulse compression
以上で明らかなように、この実施の形態1によれば、アンテナ3より放射された送信パルスr(t)の周波数スペクトルR(f)から、パルス圧縮に伴うピークレベルの低下を防ぐガード対象のレンジビンの前後のレンジサイドローブ領域における散乱波ステアリングベクトルから構成されるアレーマニフォルドAi (AM)を算出して、そのアレーマニフォルドAi (AM)からサイドローブ干渉波の相関行列Riを算出し、その相関行列Riの逆行列Ri −1を用いて、AD変換器5から出力された受信信号z(t)をパルス圧縮するように構成したので、距離分解能の劣化や探知距離の縮退を招くことなく、低レンジサイドローブ特性を実現することができる効果を奏する。
As is apparent from the above, according to the first embodiment, the frequency spectrum R (f) of the transmission pulse r (t) radiated from the antenna 3 is used as a guard target for preventing the peak level from being lowered due to pulse compression. An array manifold A i (AM) composed of scattered wave steering vectors in the range side lobe region before and after the range bin is calculated, and a correlation matrix R i of the side lobe interference wave is calculated from the array manifold A i (AM). Since the reception signal z (t) output from the
即ち、この実施の形態1によれば、着目レンジビンの前後のサイドローブ干渉波を考え、そのサイドローブ干渉波の受信電力を最小化するようなパルス圧縮ウェイトw(pro)を利用して、受信信号z(t)のパルス圧縮を行うものである。特にパルス圧縮ウェイトw(pro)の計算では、Diagonal Loadingを行ったサイドローブ干渉波の相関行列Riの逆行列Ri −1を利用するものである。このため、所望の低レンジサイドローブ特性を達成しつつ、窓関数を用いた場合と比べて、圧縮後パルス幅の拡大やミスマッチ損失を小さくすることができるという効果を奏する。 That is, according to the first embodiment, the side lobe interference wave before and after the target range bin is considered, and the pulse compression weight w (pro) that minimizes the reception power of the side lobe interference wave is used to receive the signal. The signal z (t) is subjected to pulse compression. In particular, in the calculation of the pulse compression weight w (pro) , the inverse matrix R i −1 of the correlation matrix R i of the sidelobe interference wave subjected to the diagonal loading is used. For this reason, there is an effect that the pulse width after compression and mismatch loss can be reduced as compared with the case where the window function is used while achieving the desired low range sidelobe characteristics.
実施の形態2.
図3はこの発明の実施の形態2によるレーダ装置のパルス圧縮器6を示す構成図であり、図において、図2と同一符号は同一または相当部分を示すので説明を省略する。
相関行列算出部21はアレーマニフォルド定義部14により定義されたアレーマニフォルドAi (AM)からサイドローブ干渉波の相関行列Ri(ただし、ρ=0)を算出する処理を実施する。
パルス圧縮ウェイト算出部22は相関行列算出部21により算出された相関行列Riの一部を構成する固有ベクトルに関する行列Eiからサイドローブ干渉波ランクrankの数(指定数)の固有ベクトルが取り出された行列Ei (rank)を用いて、その相関行列Riに関する射影行列Pnulliを算出し、その射影行列Pnulliを用いて、各レンジビンに対応するパルス圧縮ウェイトw(pro)を算出する処理を実施する。
3 is a block diagram showing a
The correlation
The pulse compression
図3の例では、パルス圧縮器6の構成要素である離散フーリエ変換部11,12、スペクトル積算出部13、アレーマニフォルド定義部14、相関行列算出部21、パルス圧縮ウェイト算出部22及びコヒーレント積分処理部17のそれぞれが専用のハードウェア(例えば、CPUを実装している半導体集積回路、あるいは、ワンチップマイコン)で構成されているものを想定しているが、パルス圧縮器6がコンピュータで構成されていてもよい。
パルス圧縮器6がコンピュータで構成されている場合、離散フーリエ変換部11,12、スペクトル積算出部13、アレーマニフォルド定義部14、相関行列算出部21、パルス圧縮ウェイト算出部22及びコヒーレント積分処理部17の処理内容を記述しているプログラムをコンピュータのメモリに格納し、当該コンピュータのCPUが当該メモリに格納されているプログラムを実行するようにすればよい。
In the example of FIG. 3, discrete
When the
次に動作について説明する。
上記実施の形態1と比べて、パルス圧縮器6の相関行列算出部15及びパルス圧縮ウェイト算出部16が、相関行列算出部21及びパルス圧縮ウェイト算出部22に代わっている以外は同様であるため、相関行列算出部21及びパルス圧縮ウェイト算出部22の処理内容だけを説明する。
相関行列算出部21は、アレーマニフォルド定義部14がアレーマニフォルドAi (AM)を定義すると、そのアレーマニフォルドAi (AM)からサイドローブ干渉波の相関行列Riを算出する。
Next, the operation will be described.
Compared to the first embodiment, the correlation
When the array
この実施の形態2では、Diagonal Loadingを行わないので(ρ=0)、遅延時間t1 (AM)に関する相関行列R1に着目すると、相関行列R1は固有値・固有ベクトルによって、下記の式(28)のように表される。ただし、GEは昇順に並べられた固有値を対角項に有する対角行列であり、E1は固有値に対応する固有ベクトルが列ベクトルとして配されている行列である。
このとき、固有ベクトルは互いに直交するので、下記の式(29)が成立する。
In the second embodiment, since the diagonal loading is not performed (ρ = 0), when attention is paid to the correlation matrix R 1 regarding the delay time t 1 (AM) , the correlation matrix R 1 is expressed by the following equation (28 ). However, G E is a diagonal matrix with eigenvalues ordered in ascending order in diagonal, E 1 is a matrix eigenvector corresponding to the eigenvalue is arranged as a column vector.
At this time, since the eigenvectors are orthogonal to each other, the following equation (29) is established.
パルス圧縮ウェイト算出部22は、相関行列算出部15が相関行列Riを算出すると、その相関行列Riを構成する行列Eiからサイドローブ干渉波ランクrankの数の固有ベクトルが取り出された行列Ei (rank)を用いて、相関行列Ri(ただし、ρ=0)に関する射影行列Pnulliを算出し、その射影行列Pnulliを用いて、各レンジビンに対応するパルス圧縮ウェイトw(pro)を算出する。
具体的には、以下のようにして、各レンジビンに対応するパルス圧縮ウェイトw(pro)を算出する。
When the correlation
Specifically, the pulse compression weight w (pro) corresponding to each range bin is calculated as follows.
まず、相関行列R1(ただし、ρ=0)に関する射影行列Pnull1は、下記の式(30)のように求まる。
ただし、E1 (rank)はE1よりサイドローブ干渉波ランクrankの数の固有ベクトルが取り出された行列である。
また、式(28)を用いると、遅延時間ti (AM)に関する相関行列Ri(ρ=0)は、下記の式(31)のように与えられる。
First, the projection matrix P null1 regarding the correlation matrix R 1 (where ρ = 0) is obtained as in the following Expression (30).
However, E 1 (rank) is a matrix in which eigenvectors of the number of side lobe interference wave ranks rank are extracted from E 1 .
Further, when Expression (28) is used, a correlation matrix R i (ρ = 0) related to the delay time t i (AM ) is given as Expression (31) below.
ここで、式(31)のDiE1に関して、Di HDi=Iと式(29)より、下記の式(32)が成立する。
式(32)より、DiE1の列ベクトルが互いに直交していることがわかる。即ち、式(31)において、DiE1の列ベクトルは相関行列Ri(ρ=0)の固有ベクトルとなっており、下記の式(33)が成立する。ただし、Eiは相関行列Ri(ρ=0)の固有ベクトルが列ベクトルとして配されている行列である。
Here, with respect to D i E 1 in equation (31), the following equation (32) is established from D i H D i = I and equation (29).
From Equation (32), it can be seen that the column vectors of D i E 1 are orthogonal to each other. That is, in Expression (31), the column vector of D i E 1 is an eigenvector of the correlation matrix R i (ρ = 0), and the following Expression (33) is established. However, E i is a matrix in which eigenvectors of the correlation matrix R i (ρ = 0) are arranged as column vectors.
よって、相関行列Ri(ただし、ρ=0)に関する射影行列Pnulliは、下記の式(34)のように求まる。
Therefore, the projection matrix P null regarding the correlation matrix R i (where ρ = 0) is obtained as in the following Expression (34).
したがって、相関行列Ri(ただし、ρ=0)に関する射影行列Pnulliは、遅延時間ti (AM)から直接求めずに、遅延時間t1 (AM)に関する射影行列Pnull1を求めておき、遅延時間ti (AM)に関する位相回転行列Diを射影行列Pnull1の両側から乗じることで求められることがわかる。
そのため、パルス圧縮ウェイト算出部22は、第m番目のレンジビンに対応するパルス圧縮ウェイトwm (pro)を下記の式(35)のように与える。
ただし、βmは規格化係数である。
Therefore, the projection matrix P null1 relating to the correlation matrix R i (where ρ = 0) is not obtained directly from the delay time t i (AM) , but the projection matrix P null1 relating to the delay time t 1 (AM) is obtained. It can be seen that the phase rotation matrix D i related to the delay time t i (AM) is obtained by multiplying from both sides of the projection matrix P null1 .
Therefore, the pulse compression
Where β m is a normalization factor.
この実施の形態2によれば、着目レンジビンの前後のサイドローブ干渉波を考え、そのサイドローブ干渉波の受信電力を最小化するようなパルス圧縮ウェイトw(pro)を利用して、受信信号z(t)のパルス圧縮を行うものである。特にパルス圧縮ウェイトw(pro)の計算では、サイドローブ干渉波の相関行列Riとサイドローブ干渉波ランクrankによる射影行列Pnulliを利用するものである。このため、所望の低レンジサイドローブ特性を達成しつつ、窓関数を用いた場合と比べて、圧縮後パルス幅の拡大やミスマッチ損失を小さくすることができるという効果を奏する。 According to the second embodiment, the side lobe interference wave before and after the target range bin is considered, and the received signal z is used by using the pulse compression weight w (pro) that minimizes the reception power of the side lobe interference wave. The pulse compression of (t) is performed. In particular, in the calculation of the pulse compression weight w (pro) , a correlation matrix R i of the side lobe interference wave and a projection matrix P null by the side lobe interference wave rank are used. For this reason, there is an effect that the pulse width after compression and mismatch loss can be reduced as compared with the case where the window function is used while achieving the desired low range sidelobe characteristics.
実施の形態3.
図4はこの発明の実施の形態3によるレーダ装置を示す構成図であり、図において、図1と同一符号は同一または相当部分を示すので説明を省略する。
パルス圧縮器31はパルス圧縮器6と異なる手法で、AD変換器5から出力された受信信号z(t)をパルス圧縮する処理器である。
パルス圧縮器31としては、例えば、窓関数を利用してパルス圧縮を行う従来のパルス圧縮器などが用いられる。
なお、パルス圧縮器31は第2のパルス圧縮手段を構成している。
Embodiment 3 FIG.
FIG. 4 is a block diagram showing a radar apparatus according to Embodiment 3 of the present invention. In the figure, the same reference numerals as those in FIG.
The
As the
The
パルス圧縮出力選択器32はレンジビン毎に、パルス圧縮器6によりパルス圧縮された受信信号z’(t)又はパルス圧縮器31によりパルス圧縮された受信信号z’(t)を選択し、その選択した受信信号z’(t)を目標検出器7に出力する処理を実施する。
例えば、レンジビンが0〜100まであるとき、レンジビン0〜10についてはパルス圧縮器6によりパルス圧縮された受信信号z’(t)を選択し、レンジビン11〜100についてはパルス圧縮器31によりパルス圧縮された受信信号z’(t)を選択するなどの態様が考えられる。レンジビンとパルス圧縮器6,31の対応関係は予め設定されているものとする。
なお、パルス圧縮出力選択器32は選択手段を構成している。
The pulse
For example, when there are range bins from 0 to 100, the received signal z ′ (t) pulse-compressed by the
The pulse
次に動作について説明する。
この実施の形態3が想定するシナリオとして、例えば、艦船搭載型レーダのパルス送信方向に島嶼があるような伝搬環境が考えられる。
島嶼によるクラッタ電力は一般に高い電力で受信されるため、そのレンジサイドローブ領域に存在する目標のSIRが劣化してしまう場合がある。
一方、島嶼によるクラッタのレンジサイドローブ領域外では、同クラッタのような高いクラッタ電力が受信される可能性は小さい。
以上のようなシナリオで適用されるのが実施の形態3であり、島嶼周辺のレンジビンでは、従来のパルス圧縮と比べて、演算負荷が高くても、高性能なパルス圧縮特性が達成されるパルス圧縮器6を使用し、島嶼から離れたレンジビンでは、演算負荷が低いパルス圧縮を行うパルス圧縮器31を使用するものとする。
Next, the operation will be described.
As a scenario assumed by the third embodiment, for example, a propagation environment in which there is an island in the pulse transmission direction of a ship-borne radar is conceivable.
Since clutter power due to islands is generally received at high power, the target SIR existing in the range sidelobe region may deteriorate.
On the other hand, outside the range side lobe region of the clutter by the island, there is little possibility of receiving high clutter power like the clutter.
The third embodiment is applied in the above-described scenario, and in the range bin around the island, a pulse that achieves high-performance pulse compression characteristics even when the calculation load is high compared to the conventional pulse compression. In the range bin that uses the
パルス圧縮器31は、上述したように、パルス圧縮器6と異なる手法で、AD変換器5から出力された受信信号z(t)をパルス圧縮する。例えば、窓関数を利用するパルス圧縮を行う。
このパルス圧縮では、パルス圧縮器6におけるコヒーレント積分処理部17に相当する手段として、高速な離散逆フーリエ変換を用いることができる。
As described above, the
In this pulse compression, high-speed discrete inverse Fourier transform can be used as means corresponding to the coherent
パルス圧縮出力選択器32は、パルス圧縮器6によりパルス圧縮された受信信号z’(t)と、パルス圧縮器31によりパルス圧縮された受信信号z’(t)とを入力し、レンジビン毎に、いずれか一方の受信信号z’(t)を選択して目標検出器7に出力する。
上記のシナリオであれば、島嶼周辺のレンジビンでは、パルス圧縮器6によりパルス圧縮された受信信号z’(t)を選択し、それ以外のレンジビンでは、パルス圧縮器31によりパルス圧縮された受信信号z’(t)を選択する。
The pulse
In the above scenario, the received signal z ′ (t) pulse-compressed by the
以上で明らかなように、この実施の形態3によれば、パルス圧縮出力選択器32が、レンジビン毎に、パルス圧縮器6によりパルス圧縮された受信信号z’(t)又はパルス圧縮器31によりパルス圧縮された受信信号z’(t)を選択して、その選択した受信信号z’(t)を目標検出器7に出力するように構成したので、演算負荷の高いパルス圧縮器6を全レンジビンではなく、所定範囲のレンジビンのみでパルス圧縮を行い、残りの範囲では演算負荷の低いパルス圧縮器31でパルス圧縮を行うことが可能になり、全レンジビンでパルス圧縮器6を使用する場合よりも、演算負荷を低減することができる効果を奏する。
As is apparent from the above, according to the third embodiment, the pulse
実施の形態4.
この実施の形態4では、Mt(Mtは2以上の整数)個のMIMO送信パルスを空中に同時に放射するMIMOレーダ送信機を搭載し、M個のパルス圧縮器が受信信号のパルス圧縮を行うレーダ装置について説明する。
Embodiment 4 FIG.
In the fourth embodiment, a MIMO radar transmitter that simultaneously emits Mt (Mt is an integer of 2 or more) MIMO transmission pulses in the air is mounted, and a radar in which M pulse compressors perform pulse compression of a received signal. The apparatus will be described.
図5はこの発明の実施の形態4によるレーダ装置を示す構成図であり、図において、図1と同一符号は同一または相当部分を示すので説明を省略する。
MIMOレーダ送信機41は所定の送信周波数信号を生成し、所定のPRI間隔で、その送信周波数信号を用いて、Mt個(アンテナ数と一致する個数)のMIMO送信パルスrj(t)(j=1,…,Mt)を生成する処理を実施する。
なお、MIMOレーダ送信機41は変調パラメータを変えることで、Mt個のMIMO送信パルスを生成する。
5 is a block diagram showing a radar apparatus according to Embodiment 4 of the present invention. In the figure, the same reference numerals as those in FIG.
The
The
MIMOレーダアンテナ42はMt個のアンテナから構成されており、MIMOレーダ送信機41により生成されたMt個のMIMO送信パルスrj(t)を所定のビーム指向方向に向けて空中に同時に放射する一方、目標又はクラッタに反射されたMIMO送信パルスの散乱波を受信する。
なお、MIMOレーダ送信機41及びMIMOレーダアンテナ42はパルス信号送信手段を構成しており、MIMOレーダアンテナ42は散乱波受信手段を構成している。
The
The
パルス圧縮器43−1〜43−MtはPRI周期毎に、AD変換器5から出力された受信信号z(t)をパルス圧縮する処理器である。
即ち、パルス圧縮器43−1〜43−MtはMt個のMIMO送信パルスrj(t)の周波数スペクトルRj(f)から、パルス圧縮に伴うピークレベルの低下を防ぐガード対象のレンジビンの前後のレンジサイドローブ領域における散乱波ステアリングベクトルから構成されるアレーマニフォルドBi (AM)を算出して、そのアレーマニフォルドBi (AM)からサイドローブ干渉波の相関行列Riを算出し、その相関行列Riの逆行列Ri −1を用いて、AD変換器5から出力された受信信号z(t)をパルス圧縮する処理を実施する。
なお、パルス圧縮器43−1〜43−Mtはパルス圧縮手段を構成している。
The pulse compressors 43-1 to 43-Mt are processors that pulse-compress the received signal z (t) output from the
That is, the pulse compressors 43-1 to 43-Mt use the frequency spectrum R j (f) of the Mt MIMO transmission pulses r j (t) before and after the range bins to be guarded to prevent the peak level from being lowered due to pulse compression. An array manifold B i (AM) composed of scattered wave steering vectors in the range side lobe region is calculated, a correlation matrix R i of the side lobe interference wave is calculated from the array manifold B i (AM) , and the correlation is calculated. matrix using the inverse matrix R i -1 of R i, carries out a process of pulse compression of the received signal z outputted from the AD converter 5 (t).
The pulse compressors 43-1 to 43-Mt constitute pulse compression means.
MIMOレーダ送信ビーム合成器44はパルス圧縮器43−1〜43−Mtによりパルス圧縮された受信信号z’(t)を合成し、合成後の受信信号を目標検出器7に出力する処理を実施する。なお、MIMOレーダ送信ビーム合成器44は合成手段を構成している。
The MIMO radar
図6はこの発明の実施の形態4によるレーダ装置のパルス圧縮器43−mtを示す構成図であり、図において、図2と同一符号は同一または相当部分を示すので説明を省略する。
第2の周波数スペクトル算出部である離散フーリエ変換部51−1〜51−MtはMIMOレーダ送信機41から出力されたMIMO送信パルスr1(t)〜rMt(t)をそれぞれレプリカ#1〜#Mtとして、そのレプリカ#1〜#Mtを離散フーリエ変換して、そのレプリカ#1〜#Mtの周波数スペクトルR1(f)〜RMt(f)を算出する処理を実施する。
FIG. 6 is a block diagram showing a pulse compressor 43-mt of a radar apparatus according to Embodiment 4 of the present invention. In the figure, the same reference numerals as those in FIG.
Discrete Fourier transform units 51-1 to 51-Mt, which are second frequency spectrum calculation units, convert the MIMO transmission pulses r 1 (t) to r Mt (t) output from the
スペクトル積算出部52は離散フーリエ変換部51−1〜51−Mtにより算出されたレプリカ#1〜#Mtの周波数スペクトルR1(f)〜RMt(f)のうち、第mt番目のレプリカ#mtの周波数スペクトルRmt(f)の複素共役R* mt(f)を離散フーリエ変換部11により算出された受信信号z(t)の周波数スペクトルZ(f)に乗算してスペクトル積X(f)を算出する処理を実施する。
アレーマニフォルド定義部53は離散フーリエ変換部51−1〜51−Mtにより算出されたレプリカ#1〜#Mtの周波数スペクトルR1(f)〜RMt(f)から散乱波ステアリングベクトルを算出して、パルス圧縮に伴うピークレベルの低下を防ぐガード対象のレンジビン(ガード領域)の前後のレンジサイドローブ領域における散乱波ステアリングベクトルから構成されるアレーマニフォルドBi (AM)を定義する処理を実施する。
The spectrum
The array
図6の例では、パルス圧縮器43−mtの構成要素である離散フーリエ変換部11,51−1〜51−Mt、スペクトル積算出部52、アレーマニフォルド定義部53、相関行列算出部15、パルス圧縮ウェイト算出部16及びコヒーレント積分処理部17のそれぞれが専用のハードウェア(例えば、CPUを実装している半導体集積回路、あるいは、ワンチップマイコン)で構成されているものを想定しているが、パルス圧縮器43−mtがコンピュータで構成されていてもよい。
パルス圧縮器43−mtがコンピュータで構成されている場合、離散フーリエ変換部11,51−1〜51−Mt、スペクトル積算出部52、アレーマニフォルド定義部53、相関行列算出部15、パルス圧縮ウェイト算出部16及びコヒーレント積分処理部17の処理内容を記述しているプログラムをコンピュータのメモリに格納し、当該コンピュータのCPUが当該メモリに格納されているプログラムを実行するようにすればよい。
In the example of FIG. 6, discrete
When the pulse compressor 43-mt is configured by a computer, the discrete
次に動作について説明する。
MIMOレーダ送信機41は、所定の送信周波数信号を生成し、所定のPRI間隔で、その送信周波数信号を用いて、互いに波形が異なるMt個のMIMO送信パルスrj(t)(j=1,…,Mt)を生成する。
Mt個のMIMO送信パルスrj(t)は、いずれも中心周波数がfcであり、一定のパルス幅τucである。
この実施の形態4では、第mt番目のMIMO送信パルスrmt(t)の変調帯域幅がBmtであるとする。
Next, the operation will be described.
The
Mt number of MIMO transmit pulse r j (t) are both center frequency f c, a constant pulse width tau uc.
In the fourth embodiment, it is assumed that the modulation bandwidth of the mtth MIMO transmission pulse r mt (t) is B mt .
デュプレクサ2は、MIMOレーダ送信機41からMt個のMIMO送信パルスrj(t)を受けると、Mt個のMIMO送信パルスrj(t)をMIMOレーダアンテナ42に出力する。
これにより、MIMOレーダアンテナ42を構成しているMt個のアンテナから、Mt個のMIMO送信パルスrj(t)が所定のビーム指向方向に向けて空中に同時に放射される。
その後、目標又はクラッタに反射されたMIMO送信パルスrj(t)の散乱波がMIMOレーダアンテナ42に到来し、MIMOレーダアンテナ42によりMIMO送信パルスrj(t)の散乱波が受信される。
Upon receiving Mt MIMO transmission pulses r j (t) from the
As a result, Mt MIMO transmission pulses r j (t) are simultaneously radiated from the Mt antennas constituting the
Thereafter, the scattered wave of the MIMO transmission pulse r j (t) reflected by the target or the clutter arrives at the
デュプレクサ2は、MIMOレーダアンテナ42により受信された散乱波を受信機4に出力する。
受信機4は、MIMOレーダアンテナ42により受信された散乱波をデュプレクサ2経由で受けると、その散乱波を示す信号の周波数をベースバンド帯に変換して、ベースバンド帯のアナログ受信信号を生成する。
AD変換器5は、受信機4がアナログ受信信号を生成すると、そのアナログ受信信号をAD変換して、ディジタルの受信信号z(t)をパルス圧縮器43−1〜43−Mtに出力する。
The
When the receiver 4 receives the scattered wave received by the
When the receiver 4 generates an analog reception signal, the
パルス圧縮器43−1〜43−Mtは、PRI周期毎に、AD変換器5から出力された受信信号z(t)をパルス圧縮し、パルス圧縮後の受信信号z’(t)をMIMOレーダ送信ビーム合成器44に出力する。
即ち、パルス圧縮器43−1〜43−Mtは、Mt個のMIMO送信パルスrj(t)の周波数スペクトルRj(f)から、パルス圧縮に伴うピークレベルの低下を防ぐガード対象のレンジビンの前後のレンジサイドローブ領域における散乱波ステアリングベクトルから構成されるアレーマニフォルドBi (AM)を算出して、そのアレーマニフォルドBi (AM)からサイドローブ干渉波の相関行列Riを算出し、その相関行列Riの逆行列Ri −1を用いて、AD変換器5から出力された受信信号z(t)をパルス圧縮する。
The pulse compressors 43-1 to 43-Mt pulse-compress the received signal z (t) output from the
That is, the pulse compressors 43-1 to 43-Mt use the frequency spectrum R j (f) of the Mt MIMO transmission pulses r j (t) to detect the range bins to be guarded that prevent the peak level from being lowered due to the pulse compression. An array manifold B i (AM) composed of scattered wave steering vectors in the front and rear range side lobe regions is calculated, and a correlation matrix R i of the side lobe interference wave is calculated from the array manifold B i (AM). The received signal z (t) output from the
MIMOレーダ送信ビーム合成器44は、パルス圧縮器43−1〜43−Mtからパルス圧縮後の受信信号z’(t)を受けると、それらの受信信号z’(t)を合成し、合成後の受信信号を目標検出器7に出力する。
目標検出器7は、MIMOレーダ送信ビーム合成器44から合成後の受信信号z’(t)を受けると、その受信信号z’(t)から目標を検出する。
When receiving a received signal z ′ (t) after pulse compression from the pulse compressors 43-1 to 43-Mt, the MIMO radar
Upon receiving the combined received signal z ′ (t) from the MIMO radar
以下、第mt番目のパルス圧縮器43−mtの処理内容を具体的に説明する。
ただし、以下の説明では、説明の簡単化のため、送信と受信の切替タイミング付近で発生する受信パルスのエクリプスについては考えない。
まず、AD変換器5からPRI間隔で出力される受信信号z(t)は、K個の散乱波を含んでおり、下記の式(36)のように表される。
式(36)において、δtjはMIMOレーダアンテナ42を構成している第j番目のアンテナの送信行路差による時間遅延である。
The processing contents of the mtth pulse compressor 43-mt will be specifically described below.
However, in the following description, for simplification of explanation, Eclipse of received pulses generated near the transmission and reception switching timing is not considered.
First, the received signal z (t) output from the
In Expression (36), δt j is a time delay due to the transmission path difference of the j-th antenna constituting the
パルス圧縮器43−mtの離散フーリエ変換部11は、AD変換器5から受信信号z(t)を受けると、その受信信号z(t)を離散フーリエ変換して、その受信信号z(t)の周波数スペクトルZ(f)を算出する。受信信号z(t)の周波数スペクトルZ(f)は、下記の式(37)のように表される。
式(37)において、Rj(f)はMIMO送信パルスrj(t)の周波数スペクトルである。
When receiving the received signal z (t) from the
In Expression (37), R j (f) is a frequency spectrum of the MIMO transmission pulse r j (t).
離散フーリエ変換部51−1〜51−Mtは、MIMOレーダ送信機41から出力されたMIMO送信パルスr1(t)〜rMt(t)をそれぞれレプリカ#1〜#Mtとし、レプリカ#1〜#Mtを離散フーリエ変換して、そのレプリカ#1〜#Mtの周波数スペクトルR1(f)〜RMt(f)を算出する。
The discrete Fourier transform units 51-1 to 51-Mt use the MIMO transmission pulses r 1 (t) to r Mt (t) output from the
スペクトル積算出部52は、離散フーリエ変換部11が受信信号z(t)の周波数スペクトルZ(f)を算出し、離散フーリエ変換部51−mtがレプリカ#mtの周波数スペクトルRmt(f)を算出すると、下記の式(38)に示すように、レプリカ#mtの周波数スペクトルRmt(f)の複素共役R* mt(f)を離散フーリエ変換部11により算出された受信信号z(t)の周波数スペクトルZ(f)に乗算してスペクトル積X(f)を算出する。
In the spectral
ここでは、周波数スペクトルZ(f),Rmt(f)がM点のサンプルより構成されているものとして、下記の式(39)のような受信信号ベクトルx(mt)を定義する。
式(38)を式(39)に代入すると、下記の式(40)のように変形することができる。
式(40)において、bk (mt)はレプリカ#mtにおける第k番目の目標に対する遅延時間tk及びドップラ周波数fk (D)から決まるステアリングベクトル、B(mt)はK個のステアリングベクトルを列ベクトルとする行列、n(mt)は受信機雑音ベクトルである。
Here, assuming that the frequency spectrums Z (f) and R mt (f) are composed of M-point samples, a received signal vector x (mt) is defined as in the following equation (39).
By substituting equation (38) into equation (39), it can be transformed into the following equation (40).
In Equation (40), b k (mt) is a steering vector determined from the delay time t k and the Doppler frequency f k (D) for the kth target in the replica #mt, and B (mt) is K steering vectors. A matrix n (mt) as a column vector is a receiver noise vector.
ここで、上記実施の形態1で示している式(6)と、上記の式(40)を比較すると、MIMOレーダ送信機41を用いる場合には、akをbk (mt)に置き換え、AをB(mt)に置き換え、nをn(mt)に置き換えた形となっていることが明らかである。
このため、この実施の形態4におけるパルス圧縮ウェイト計算の処理のながれも、上記の置き換えを念頭におくと、上記実施の形態1と同様に説明することができる。
Here, when the equation (6) shown in the first embodiment is compared with the above equation (40), when the
For this reason, the flow of the pulse compression weight calculation process in the fourth embodiment can be described in the same manner as in the first embodiment, taking the above replacement into account.
アレーマニフォルド定義部53は、離散フーリエ変換部51−1〜51−Mtがレプリカ#1〜#Mtの周波数スペクトルR1(f)〜RMt(f)を算出すると、レプリカ#1〜#Mtの周波数スペクトルR1(f)〜RMt(f)から散乱波ステアリングベクトルを算出して、パルス圧縮に伴うピークレベルの低下を防ぐガード対象のレンジビン(ガード領域)の前後のレンジサイドローブ領域における散乱波ステアリングベクトルから構成されるアレーマニフォルドBi (AM)を定義する。
具体的には、以下の通りである。
The array
Specifically, it is as follows.
ここでは、PRI内における遅延時間ti (AM)に対応する散乱波の受信信号ベクトルci (AM|mt)を考える。
ci (AM|mt)は散乱波のステアリングベクトルbi (AM|mt)と等しく、下記の式(44)のように与えられる。
ただし、b1 (AM|mt)は下記の式(45)で示されるものであり、式(45)中のδtj (TBR)はMIMOレーダ送信ビームを形成するための時間遅延の設定値である。
Here, the received signal vector c i (AM | mt) of the scattered wave corresponding to the delay time t i (AM) in the PRI is considered.
c i (AM | mt) is the scattered wave steering vectors b i | equal (AM mt), given as the following equation (44).
However, b 1 (AM | mt) is expressed by the following equation (45), and δ tj (TBR) in equation (45 ) is a set value of a time delay for forming a MIMO radar transmission beam. is there.
式(44)より、遅延時間ti−l (AM),ti+l (AM)に対応する散乱波の受信信号ベクトルci−l (AM|mt),ci+l (AM|mt)は、下記の式(46)(47)のように与えられる。
From Expression (44), the received signal vectors c i−l (AM | mt) and c i + l (AM | mt) of the scattered wave corresponding to the delay times t i−l (AM) and t i + l (AM) are (46) and (47).
アレーマニフォルド定義部53は、上記の式(21)のガード行列Gによって、パルス圧縮のメインローブ付近の遅延時間に対応するステアリングベクトルをゼロベクトルにするものとし、下記の式(48)に示すように、サイドローブ領域の散乱点ステアリングベクトルから構成されるアレーマニフォルドBi (AM|mt)を定義する。
The array
相関行列算出部15は、アレーマニフォルド定義部53がアレーマニフォルドBi (AM|mt)を定義すると、そのアレーマニフォルドBi (AM|mt)とDLLからサイドローブ干渉波の相関行列Riを算出する。
即ち、相関行列算出部15は、下記の式(49)に示すように、遅延時間ti (AM)に関するパルス圧縮におけるサイドローブ干渉波の相関行列Riを定義する。
When the array
That is, the correlation
式(49)より、遅延時間ti (AM)に関する相関行列Ri (mt)は、遅延時間t1 (AM)に関する相関行列R1 (mt)を求めておけば、遅延時間ti (AM)に関する位相回転行列Diを利用して簡単に求められることがわかる。
ここで、サイドローブの干渉波を抑圧するための行列を相関行列Ri (mt)の逆行列(Ri (mt))−1とすると、相関行列Ri (mt)の逆行列(Ri (mt))−1は、下記の式(50)のように与えられる。
From equation (49), the correlation matrix for the delay time t i (AM) R i ( mt) is, if the correlation matrix for the delay time t 1 (AM) R 1 ( mt), the delay time t i (AM It can be seen that it can be easily obtained using the phase rotation matrix D i for ) .
Here, when the matrix for suppressing interference waves sidelobe and inverse matrix (R i (mt)) -1 of the correlation matrix R i (mt), the inverse matrix of the correlation matrix R i (mt) (R i (Mt) ) -1 is given by the following equation (50).
パルス圧縮ウェイト算出部16は、相関行列算出部15がサイドローブ干渉波の相関行列Ri (mt)を算出すると、その相関行列Ri (mt)の逆行列(Ri (mt))−1を用いて、各レンジビンに対応するパルス圧縮ウェイトw(pro)を算出する。
ここで、遅延時間ti (AM)に関する相関行列Ri (mt)の逆行列(Ri (mt))−1は、式(50)より、遅延時間ti (AM)から直接求めずに、遅延時間t1 (AM)に関する逆行列(R1 (mt))−1を求めておき、遅延時間ti (AM)に関する位相回転行列Diを逆行列(R1 (mt))−1の両側から乗じることで求められることがわかる。
When the correlation
Here, the inverse matrix of the correlation matrix for the delay time t i (AM) R i ( mt) (R i (mt)) -1 , from expression (50), without requiring direct from the delay time t i (AM) Inverse matrix (R 1 (mt) ) −1 relating to delay time t 1 (AM) is obtained, and phase rotation matrix D i relating to delay time t i (AM) is obtained as inverse matrix (R 1 (mt) ) −1. It turns out that it is calculated by multiplying from both sides.
そのため、パルス圧縮ウェイト算出部16は、第m番目のレンジビンに対応するパルス圧縮ウェイトwm (pro|mt)を下記の式(51)のように与える。0044
Therefore, the pulse compression
コヒーレント積分処理部17は、パルス圧縮ウェイト算出部16が各レンジビンに対応するパルス圧縮ウェイトw(pro|mt)を算出すると、そのパルス圧縮ウェイトw(pro|mt)を用いて、スペクトル積算出部13により算出されたスペクトル積X(f)をコヒーレント積分し、その積分結果をパルス圧縮後の受信信号z’(t)として、目標検出器7に出力する。
When the pulse compression
この実施の形態4によれば、MIMOレーダ送信機41を用いる際、着目レンジビンの前後のサイドローブ干渉波を考え、そのサイドローブ干渉波の受信電力を最小化するようなパルス圧縮ウェイトw(pro|mt)を利用して、受信信号z(t)のパルス圧縮を行うものである。特にパルス圧縮ウェイトw(pro|mt)の計算では、Diagonal Loadingを行ったサイドローブ干渉波の相関行列Ri (mt)の逆行列(Ri (mt))−1を利用するものである。このため、自己相関のみならず、相互相関によるレンジサイドローブも抑圧して、所望の低レンジサイドローブ特性を達成しつつ、窓関数を用いた場合と比べて、圧縮後パルス幅の拡大やミスマッチ損失を小さくすることができるという効果を奏する。
According to the fourth embodiment, when the
実施の形態5.
図7はこの発明の実施の形態5によるレーダ装置のパルス圧縮器43−mtを示す構成図であり、図において、図3及び図6と同一符号は同一または相当部分を示すので説明を省略する。
上記実施の形態4と比べて、パルス圧縮器43−mtの相関行列算出部15及びパルス圧縮ウェイト算出部16が、相関行列算出部21及びパルス圧縮ウェイト算出部22に代わっている以外は同様であるため、相関行列算出部21及びパルス圧縮ウェイト算出部22の処理内容だけを説明する。
相関行列算出部21は、アレーマニフォルド定義部53がアレーマニフォルドBi (AM|mt)を定義すると、そのアレーマニフォルドBi (AM|mt)からサイドローブ干渉波の相関行列Riを算出する。
7 is a block diagram showing a pulse compressor 43-mt of a radar apparatus according to
Compared to the fourth embodiment, except that the correlation
When the array
この実施の形態5では、Diagonal Loadingを行わないので(ρ=0)、上記実施の形態2と同様に、遅延時間ti (AM)に関する相関行列Ri (mt)から、サイドローブ干渉波の抑圧に用いる射影行列Pnulli (mt)を下記の式(52)のように求めることができる。
In the fifth embodiment, since the diagonal loading is not performed (ρ = 0), as in the second embodiment, from the correlation matrix R i (mt) regarding the delay time t i (AM) , A projection matrix P null (mt) used for suppression can be obtained as in the following equation (52).
したがって、相関行列Ri (mt)(ただし、ρ=0)に関する射影行列Pnulli (mt)は、遅延時間ti (AM)から直接求めずに、遅延時間t1 (AM)に関する射影行列Pnull1 (mt)を求めておき、遅延時間ti (AM)に関する位相回転行列Diを射影行列Pnull1 (mt)の両側から乗じることで求められることがわかる。
そのため、パルス圧縮ウェイト算出部22は、第m番目のレンジビンに対応するパルス圧縮ウェイトwm (pro|mt)を下記の式(53)のように与える。
ただし、βm (mt)は規格化係数である。
Therefore, the projection matrix P null (mt ) relating to the correlation matrix R i (mt) (where ρ = 0) is not directly obtained from the delay time t i (AM) , but the projection matrix P relating to the delay time t 1 (AM). to previously obtain the null1 (mt), it can be seen that the obtained by multiplying the phase rotation matrix D i about a delay time t i (AM) from both sides of the projection matrix P null1 (mt).
Therefore, the pulse compression
Where β m (mt) is a normalization coefficient.
この実施の形態5によれば、MIMOレーダ送信機41を用いる際、着目レンジビンの前後のサイドローブ干渉波を考え、そのサイドローブ干渉波の受信電力を最小化するようなパルス圧縮ウェイトw(pro|mt)を利用して、受信信号z(t)のパルス圧縮を行うものである。特にパルス圧縮ウェイトw(pro|mt)の計算では、サイドローブ干渉波の相関行列Riとサイドローブ干渉波ランクrankによる射影行列Pnulli (mt)を利用するものである。このため、自己相関のみならず、相互相関によるレンジサイドローブも抑圧して、所望の低レンジサイドローブ特性を達成しつつ、窓関数を用いた場合と比べて、圧縮後パルス幅の拡大やミスマッチ損失を小さくすることができるという効果を奏する。
According to the fifth embodiment, when the
実施の形態6.
図8はこの発明の実施の形態6によるレーダ装置を示す構成図であり、図において、図5と同一符号は同一または相当部分を示すので説明を省略する。
パルス圧縮器61−1〜61−Mtは図4のパルス圧縮器31に相当する処理器であり、パルス圧縮器43−1〜43−Mt(第1のパルス圧縮手段)と異なる手法で、AD変換器5から出力された受信信号z(t)をパルス圧縮する。
なお、パルス圧縮器61−1〜61−Mtは第2のパルス圧縮手段を構成している。
FIG. 8 is a block diagram showing a radar apparatus according to
The pulse compressors 61-1 to 61-Mt are processors corresponding to the
The pulse compressors 61-1 to 61-Mt constitute second pulse compression means.
パルス圧縮出力選択器62−1〜62−Mt/2は図4のパルス圧縮出力選択器32に相当する選択器であり、レンジビン毎に、パルス圧縮器43−mtによりパルス圧縮された受信信号z’(t)又はパルス圧縮器61−mtによりパルス圧縮された受信信号z’(t)を選択し、その選択した受信信号z’(t)をMIMOレーダ送信ビーム合成器44に出力する処理を実施する。
なお、パルス圧縮出力選択器62−1〜62−Mt/2は選択手段を構成している。
The pulse compression output selectors 62-1 to 62-Mt / 2 are selectors corresponding to the pulse
The pulse compression output selectors 62-1 to 62-Mt / 2 constitute selection means.
次に動作について説明する。
この実施の形態6が想定するシナリオは、上記実施の形態3と同様である。
このため、島嶼周辺のレンジビンでは、従来のパルス圧縮と比べて、演算負荷が高くても、高性能なパルス圧縮特性が達成されるパルス圧縮器43−1〜43−Mtを使用し、島嶼から離れたレンジビンでは、演算負荷が低いパルス圧縮を行うパルス圧縮器61−1〜61−Mtを使用するものとする。
Next, the operation will be described.
The scenario assumed in the sixth embodiment is the same as that in the third embodiment.
For this reason, the range bins around the island use pulse compressors 43-1 to 43-Mt that achieve high-performance pulse compression characteristics even when the calculation load is high compared to conventional pulse compression. In the remote range bin, pulse compressors 61-1 to 61-Mt that perform pulse compression with a low calculation load are used.
パルス圧縮器61−1〜61−Mtは、上述したように、パルス圧縮器43−1〜43−Mtと異なる手法で、AD変換器5から出力された受信信号z(t)をパルス圧縮する。例えば、窓関数を利用するパルス圧縮を行う。
このパルス圧縮では、パルス圧縮器43−1〜43−Mtにおけるコヒーレント積分処理部17に相当する手段として、高速な離散逆フーリエ変換を用いることができる。
As described above, the pulse compressors 61-1 to 61-Mt pulse-compress the received signal z (t) output from the
In this pulse compression, high-speed discrete inverse Fourier transform can be used as means corresponding to the coherent
パルス圧縮出力選択器62−1〜62−Mt/2は、パルス圧縮器43−1〜43−Mtによりパルス圧縮された受信信号z’(t)と、パルス圧縮器61−1〜61−Mtによりパルス圧縮された受信信号z’(t)とを入力し、レンジビン毎に、いずれか一方の受信信号z’(t)を選択してMIMOレーダ送信ビーム合成器44に出力する。
上記のシナリオであれば、島嶼周辺のレンジビンでは、パルス圧縮器43−1〜43−Mtによりパルス圧縮された受信信号z’(t)を選択し、それ以外のレンジビンでは、パルス圧縮器61−1〜61−Mtによりパルス圧縮された受信信号z’(t)を選択する。
The pulse compression output selectors 62-1 to 62-Mt / 2 include the received signal z ′ (t) pulse-compressed by the pulse compressors 43-1 to 43-Mt and the pulse compressors 61-1 to 61-Mt. The received signal z ′ (t) pulse-compressed by is input, and one of the received signals z ′ (t) is selected and output to the MIMO radar
In the above scenario, the received signal z ′ (t) pulse-compressed by the pulse compressors 43-1 to 43-Mt is selected in the range bin around the island, and the pulse compressor 61- in the other range bins. The received signal z ′ (t) pulse-compressed by 1 to 61-Mt is selected.
以上で明らかなように、この実施の形態6によれば、パルス圧縮出力選択器62−1〜62−Mt/2が、レンジビン毎に、パルス圧縮器43−1〜43−Mtによりパルス圧縮された受信信号z’(t)又はパルス圧縮器61−1〜61−Mtによりパルス圧縮された受信信号z’(t)を選択して、その選択した受信信号z’(t)をMIMOレーダ送信ビーム合成器44に出力するように構成したので、演算負荷の高いパルス圧縮器43−1〜43−Mtを全レンジビンではなく、所定範囲のレンジビンのみでパルス圧縮を行い、残りの範囲では演算負荷の低いパルス圧縮器61−1〜61−Mtでパルス圧縮を行うことが可能になり、全レンジビンでパルス圧縮器43−1〜43−Mtを使用する場合よりも、演算負荷を低減することができる効果を奏する。
As apparent from the above, according to the sixth embodiment, the pulse compression output selectors 62-1 to 62-Mt / 2 are pulse-compressed by the pulse compressors 43-1 to 43-Mt for each range bin. The received signal z ′ (t) or the received signal z ′ (t) pulse-compressed by the pulse compressors 61-1 to 61-Mt is selected, and the selected received signal z ′ (t) is transmitted by MIMO radar. Since it is configured to output to the
なお、本願発明はその発明の範囲内において、各実施の形態の自由な組み合わせ、あるいは各実施の形態の任意の構成要素の変形、もしくは各実施の形態において任意の構成要素の省略が可能である。 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 AD変換器(散乱波受信手段)、6 パルス圧縮器(パルス圧縮手段)、7 目標検出器(目標検出手段)、11 離散フーリエ変換部(第1の周波数スペクトル算出部)、12 離散フーリエ変換部(第2の周波数スペクトル算出部)、13 スペクトル積算出部、14 アレーマニフォルド定義部、15,21 相関行列算出部、16,22 パルス圧縮ウェイト算出部、17 コヒーレント積分処理部、31 パルス圧縮器(第2のパルス圧縮手段)、32 パルス圧縮出力選択器(選択手段)、41 MIMOレーダ送信機(パルス信号送信手段)、42 MIMOレーダアンテナ(パルス信号送信手段、散乱波受信手段)、43−1〜43−Mt パルス圧縮器(パルス圧縮手段、第1のパルス圧縮手段)、44 MIMOレーダ送信ビーム合成器(合成手段)、51−1〜51−Mt 離散フーリエ変換部(第2の周波数スペクトル算出部)、52 スペクトル積算出部、53 アレーマニフォルド定義部、61−1〜61−Mt パルス圧縮器(第2のパルス圧縮手段)、62−1〜62−Mt/2 パルス圧縮出力選択器(選択手段)。 DESCRIPTION OF SYMBOLS 1 Transmitter (pulse signal transmission means), 2 Duplexer (pulse signal transmission means, scattered wave reception means), 3 Antenna (pulse signal transmission means, scattered wave reception means), 4 Receiver (scattered wave reception means), 5 AD Converter (scattered wave reception means), 6 pulse compressor (pulse compression means), 7 target detector (target detection means), 11 discrete Fourier transform section (first frequency spectrum calculation section), 12 discrete Fourier transform section ( (Second frequency spectrum calculation unit), 13 spectrum product calculation unit, 14 array manifold definition unit, 15, 21 correlation matrix calculation unit, 16, 22 pulse compression weight calculation unit, 17 coherent integration processing unit, 31 pulse compressor (first) 2 pulse compression means), 32 pulse compression output selector (selection means), 41 MIMO radar transmitter (pulse signal transmission means), 42 MIMO radar antenna (pulse signal transmission means, scattered wave reception means), 43-1 to 43-Mt pulse compressor (pulse compression means, first pulse compression means), 44 MIMO radar transmission beam synthesizer (synthesis means) , 51-1 to 51-Mt discrete Fourier transform unit (second frequency spectrum calculation unit), 52 spectral product calculation unit, 53 array manifold definition unit, 61-1 to 61-Mt pulse compressor (second pulse compression) Means), 62-1 to 62-Mt / 2 pulse compression output selector (selection means).
Claims (8)
上記複数のパルス圧縮手段は、上記パルス信号送信手段により放射された複数のMIMOパルス信号の周波数スペクトルから散乱波ステアリングベクトルを算出し、パルス圧縮に伴うピークレベルの低下を防ぐガード対象のレンジビンの前後のレンジサイドローブ領域における散乱波ステアリングベクトルから構成されるアレーマニフォルドを定義して、上記アレーマニフォルドからサイドローブ干渉波の相関行列を算出し、上記相関行列を用いて、散乱波受信手段から出力された受信信号をパルス圧縮する処理を実施し、
上記複数のパルス信号送信手段によりパルス圧縮された受信信号を合成し、合成後の受信信号を目標検出手段に出力する合成手段
を設けたことを特徴とする請求項1記載のレーダ装置。 The pulse signal transmission means is composed of a MIMO radar transmitter, when a plurality of MIMO pulse signal generated by the pulse signal sending unit is simultaneously emitted in the air, the MIMO pulse signal the same number of pulse compression means Installed,
The plurality of pulse compression means calculate scattered wave steering vectors from the frequency spectra of the plurality of MIMO pulse signals radiated by the pulse signal transmission means, and before and after the guard target range bin that prevents the peak level from being lowered due to pulse compression. Define an array manifold composed of scattered wave steering vectors in the range sidelobe region, calculate the correlation matrix of the sidelobe interference wave from the array manifold, and output from the scattered wave receiving means using the correlation matrix The received signal is pulse-compressed,
2. The radar apparatus according to claim 1, further comprising a combining unit that combines the reception signals compressed by the plurality of pulse signal transmission units and outputs the combined reception signal to the target detection unit.
第m(m=1,…,M)番目のパルス圧縮手段は、
上記散乱波受信手段から出力された受信信号の周波数スペクトルを算出する第1の周波数スペクトル算出部と、上記パルス信号送信手段により放射された複数のMIMOパルス信号をそれぞれレプリカとして、そのレプリカの周波数スペクトルを算出する複数の第2の周波数スペクトル算出部と、上記複数の第2の周波数スペクトル算出部により算出された複数のレプリカの周波数スペクトルのうち、第m番目のレプリカの周波数スペクトルの複素共役を上記第1の周波数スペクトル算出部により算出された受信信号の周波数スペクトルに乗算してスペクトル積を算出するスペクトル積算出部と、上記複数の第2の周波数スペクトル算出部により算出されたレプリカの周波数スペクトルから散乱波ステアリングベクトルを算出して、パルス圧縮に伴うピークレベルの低下を防ぐガード対象のレンジビンの前後のレンジサイドローブ領域における散乱波ステアリングベクトルから構成されるアレーマニフォルドを定義するアレーマニフォルド定義部と、上記アレーマニフォルド定義部により定義されたアレーマニフォルドと対角荷重であるDLLからサイドローブ干渉波の相関行列を算出する相関行列算出部と、上記相関行列算出部により算出された相関行列の逆行列を用いて、各レンジビンに対応するパルス圧縮ウェイトを算出するパルス圧縮ウェイト算出部と、上記パルス圧縮ウェイト算出部により算出されたパルス圧縮ウェイトを用いて、上記スペクトル積算出部により算出されたスペクトル積をコヒーレント積分するコヒーレント積分処理部とから構成されていることを特徴とする請求項5記載のレーダ装置。 When M (M is an integer of 1 or more) pulse compression means is mounted,
The mth (m = 1,..., M) th pulse compression means is
A first frequency spectrum calculation section that calculates a frequency spectrum of the received signal output from the scattered wave receiving means, a plurality of MIMO pulse signal emitted as a replica, respectively by the pulse signal sending unit, the frequency spectrum of the replica And a complex conjugate of the frequency spectrum of the m-th replica among the plurality of replica frequency spectra calculated by the plurality of second frequency spectrum calculation units. From the spectrum product calculation unit that calculates the spectrum product by multiplying the frequency spectrum of the reception signal calculated by the first frequency spectrum calculation unit, and the frequency spectrum of the replica calculated by the plurality of second frequency spectrum calculation units Calculate scattered wave steering vector and calculate pulse pressure An array manifold definition section for defining an array manifold composed of scattered wave steering vectors in the range sidelobe region before and after the guarded range bin to prevent a drop in peak level due to the above, and an array manifold defined by the above array manifold definition section And a correlation matrix calculation unit that calculates a correlation matrix of a sidelobe interference wave from a DLL that is a diagonal load, and a pulse compression weight corresponding to each range bin using an inverse matrix of the correlation matrix calculated by the correlation matrix calculation unit And a coherent integration processing unit that coherently integrates the spectral product calculated by the spectral product calculation unit using the pulse compression weight calculated by the pulse compression weight calculation unit. It is characterized by The radar apparatus according to claim 5.
第m(m=1,…,M)番目のパルス圧縮手段は、
上記散乱波受信手段から出力された受信信号の周波数スペクトルを算出する第1の周波数スペクトル算出部と、上記パルス信号送信手段により放射された複数のMIMOパルス信号をそれぞれレプリカとして、そのレプリカの周波数スペクトルを算出する複数の第2の周波数スペクトル算出部と、上記複数の第2の周波数スペクトル算出部により算出されたレプリカの周波数スペクトルのうち、第m番目のレプリカの周波数スペクトルの複素共役を上記第1の周波数スペクトル算出部により算出された受信信号の周波数スペクトルに乗算してスペクトル積を算出するスペクトル積算出部と、上記複数の第2の周波数スペクトル算出部により算出されたレプリカの周波数スペクトルから散乱波ステアリングベクトルを算出して、パルス圧縮に伴うピークレベルの低下を防ぐガード対象のレンジビンの前後のレンジサイドローブ領域における散乱波ステアリングベクトルから構成されるアレーマニフォルドを定義するアレーマニフォルド定義部と、上記アレーマニフォルド定義部により定義されたアレーマニフォルドからサイドローブ干渉波の相関行列を算出する相関行列算出部と、上記相関行列算出部により算出された相関行列の一部を構成する固有ベクトルに関する行列から指定数の固有ベクトルが取り出された行列を用いて、上記相関行列に関する射影行列を算出し、上記射影行列を用いて、各レンジビンに対応するパルス圧縮ウェイトを算出するパルス圧縮ウェイト算出部と、上記パルス圧縮ウェイト算出部により算出されたパルス圧縮ウェイトを用いて、上記スペクトル積算出部により算出されたスペクトル積をコヒーレント積分するコヒーレント積分処理部とから構成されていることを特徴とする請求項5記載のレーダ装置。 When M (M is an integer of 1 or more) pulse compression means is mounted,
The mth (m = 1,..., M) th pulse compression means is
A first frequency spectrum calculation section that calculates a frequency spectrum of the received signal output from the scattered wave receiving means, a plurality of MIMO pulse signal emitted as a replica, respectively by the pulse signal sending unit, the frequency spectrum of the replica And the complex conjugate of the frequency spectrum of the m-th replica among the frequency spectra of the replicas calculated by the plurality of second frequency spectrum calculation units and the plurality of second frequency spectrum calculation units. A spectrum product calculating unit that calculates a spectrum product by multiplying the frequency spectrum of the received signal calculated by the frequency spectrum calculating unit of the received signal, and a scattered wave from the frequency spectrum of the replica calculated by the plurality of second frequency spectrum calculating units. Steering vector is calculated and An array manifold definition section that defines an array manifold composed of scattered wave steering vectors in the range sidelobe region before and after the guarded range bin to prevent a decrease in peak level, and the side from the array manifold defined by the array manifold definition section. A correlation matrix calculation unit for calculating a correlation matrix of a lobe interference wave, and a matrix obtained by extracting a specified number of eigenvectors from a matrix related to eigenvectors constituting a part of the correlation matrix calculated by the correlation matrix calculation unit, A projection matrix related to the correlation matrix is calculated, and a pulse compression weight calculation unit that calculates a pulse compression weight corresponding to each range bin using the projection matrix, and a pulse compression weight calculated by the pulse compression weight calculation unit. , Spectral integration above The radar device according to claim 5, characterized in that it is constituted spectral product calculated from the coherent integration processor for coherent integration by parts.
上記第1のパルス圧縮手段と異なる手法で、上記散乱波受信手段から出力された受信信号をパルス圧縮するM個の第2のパルス圧縮手段と、レンジビン毎に、上記第1のパルス圧縮手段における第m(m=1,…,M)番目のパルス信号送信手段によりパルス圧縮された受信信号又は上記第2のパルス圧縮手段における第m番目のパルス信号送信手段によりパルス圧縮された受信信号を選択して、その選択した受信信号を合成手段に出力する複数の選択手段とを設けたことを特徴とする請求項5記載のレーダ装置。 M (M is an integer of 1 or more) pulse compression means are used as the first pulse compression means,
In a manner different from the first pulse compression means, and M second pulse compression means for pulse compression of the received signal output from the scattered wave receiving means, for each range bin, in the first pulse compression means Select a reception signal pulse-compressed by the m-th (m = 1,..., M) -th pulse signal transmission means or a reception signal pulse-compressed by the m-th pulse signal transmission means in the second pulse compression means. 6. A radar apparatus according to claim 5, further comprising a plurality of selecting means for outputting the selected received signal to the synthesizing means.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2011281614A JP5863443B2 (en) | 2011-12-22 | 2011-12-22 | Radar equipment |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2011281614A JP5863443B2 (en) | 2011-12-22 | 2011-12-22 | Radar equipment |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2013130527A JP2013130527A (en) | 2013-07-04 |
JP5863443B2 true JP5863443B2 (en) | 2016-02-16 |
Family
ID=48908179
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2011281614A Active JP5863443B2 (en) | 2011-12-22 | 2011-12-22 | Radar equipment |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP5863443B2 (en) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6289252B2 (en) * | 2014-05-02 | 2018-03-07 | 三菱電機株式会社 | Radar equipment |
JP6296907B2 (en) * | 2014-06-06 | 2018-03-20 | 株式会社東芝 | Radar apparatus and radar signal processing method thereof |
JP6079709B2 (en) * | 2014-07-08 | 2017-02-15 | 三菱電機株式会社 | Radar equipment |
JP6462365B2 (en) * | 2015-01-13 | 2019-01-30 | 株式会社東芝 | Radar apparatus and radar signal processing method thereof |
JP6400494B2 (en) * | 2015-01-23 | 2018-10-03 | 株式会社東芝 | Radar apparatus and radar signal processing method thereof |
JP6479602B2 (en) * | 2015-08-06 | 2019-03-06 | 株式会社東芝 | Radar apparatus and radar signal processing method thereof |
CN111175705B (en) * | 2019-12-24 | 2023-05-16 | 西安电子科技大学 | Radar communication integrated design method based on spatial waveform two-dimensional weighting |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS6474479A (en) * | 1987-09-16 | 1989-03-20 | Nec Corp | Transmitter-receiver for pulse-compression radar |
US5414428A (en) * | 1993-08-06 | 1995-05-09 | Martin Marietta Corp. | Radar system with pulse compression and range sidelobe suppression preceding doppler filtering |
JP4444057B2 (en) * | 2003-09-30 | 2010-03-31 | 株式会社東芝 | Pulse compression processor |
JP5075475B2 (en) * | 2006-05-25 | 2012-11-21 | 株式会社東芝 | Transmission signal generator and radar transmitter using the same |
JP2010203965A (en) * | 2009-03-04 | 2010-09-16 | Toshiba Corp | Radar apparatus, receiver, and correlation component detector |
JP5369350B2 (en) * | 2009-08-14 | 2013-12-18 | 東京計器株式会社 | Transmission waveform generation method in pulse compression, transmission waveform generation program, and pulse compression apparatus manufactured by transmission waveform generation method |
JP5417187B2 (en) * | 2010-01-08 | 2014-02-12 | 株式会社東芝 | Correlation reception processing device |
-
2011
- 2011-12-22 JP JP2011281614A patent/JP5863443B2/en active Active
Also Published As
Publication number | Publication date |
---|---|
JP2013130527A (en) | 2013-07-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5863443B2 (en) | Radar equipment | |
US9070972B2 (en) | Wideband beam forming device; wideband beam steering device and corresponding methods | |
US7038618B2 (en) | Method and apparatus for performing bistatic radar functions | |
JP4987456B2 (en) | Radar equipment | |
JP7108929B2 (en) | Radar device and target determination method | |
JP6391888B2 (en) | Radar equipment | |
JP5554018B2 (en) | Angle measuring device, monopulse angle measuring device, monopulse radar, multi-static radar | |
JPWO2014128995A1 (en) | Passive radar device | |
JP2007064671A (en) | Radar system | |
US20220003835A1 (en) | Radar apparatus and method for determining range side lobe | |
JP2010175457A (en) | Radar apparatus | |
JP6573748B2 (en) | Radar equipment | |
Lesturgie | Some relevant applications of MIMO to radar | |
JP6164936B2 (en) | Radar equipment | |
JP7143146B2 (en) | Radar system and its radar signal processing method | |
JP5705066B2 (en) | Passive radar device | |
JP2009250952A (en) | Radar system | |
JP6271834B2 (en) | Radar equipment | |
JP5618494B2 (en) | Radar equipment | |
JP6373528B2 (en) | Radar equipment | |
JP7117557B2 (en) | radar equipment | |
JP5933330B2 (en) | Pulse compressor | |
JP7266258B2 (en) | radar equipment | |
JP2009180538A (en) | Radar system | |
JP2005195339A (en) | Radar signal processing apparatus |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20141030 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20150817 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20150915 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20151029 |
|
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: 20151124 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20151222 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 5863443 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |