JP2717331B2 - Direction estimation method using narrowband signal - Google Patents
Direction estimation method using narrowband signalInfo
- Publication number
- JP2717331B2 JP2717331B2 JP3250429A JP25042991A JP2717331B2 JP 2717331 B2 JP2717331 B2 JP 2717331B2 JP 3250429 A JP3250429 A JP 3250429A JP 25042991 A JP25042991 A JP 25042991A JP 2717331 B2 JP2717331 B2 JP 2717331B2
- Authority
- JP
- Japan
- Prior art keywords
- frequency
- analysis
- output
- estimated
- azimuth
- 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
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
- Y02E60/30—Hydrogen technology
- Y02E60/50—Fuel cells
Landscapes
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Description
【0001】[0001]
【産業上の利用分野】本発明は、水中などに存在する信
号源からの狭帯域信号を、複数のセンサからなるアレイ
により受信し、その狭帯域信号の周波数と到来方向とを
推定する方法に関し、特に、狭帯域信号の到来方向の推
定を行う方位推定方法に関するものである。BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a method for receiving a narrow-band signal from a signal source existing in water or the like by an array of a plurality of sensors and determining the frequency and the arrival direction of the narrow-band signal.
The present invention relates to an estimation method, and more particularly, to an azimuth estimation method for estimating an arrival direction of a narrowband signal.
【0002】[0002]
【従来技術】従来技術の狭帯域による方位推定方法を説
明するための方位推定装置の一例を図3にブロック図で
示す。信号源からの狭帯域信号の到来方向を推定する方
法の一つとして、複数のセンサを直線上に配列し、各セ
ンサの受信信号を周波数分析により複数の周波数帯域に
分割に、その各帯域に対し周波数領域整相により複数の
ビ−ムを形成し、各ビ−ムの出力の狭帯域分析結果の絶
対レベルを比較することにより、到来方向を求める方法
が用いられていた。2. Description of the Related Art FIG. 3 is a block diagram showing an example of an azimuth estimating apparatus for explaining a conventional azimuth estimating method using a narrow band. As one method of estimating the arrival direction of a narrowband signal from a signal source, a plurality of sensors are arranged on a straight line, and a reception signal of each sensor is divided into a plurality of frequency bands by frequency analysis. On the other hand, a method has been used in which a plurality of beams are formed by frequency domain phasing, and the arrival direction is obtained by comparing the absolute levels of the narrow band analysis results of the output of each beam.
【0003】この方位推定装置は、狭帯域信号を対象と
し、N個のセンサからなるラインアレイ31を用いて、
信号の到来方向を推定する装置である。簡単の為、ライ
ンアレイ31のN個のセンサ311〜31Nは、図4に示
すように、y軸上に配列されているものとする。又、信
号の到来方向をラインアレイ31の法線方向からの角度
θで表す。This azimuth estimating apparatus targets a narrow band signal, and uses a line array 31 composed of N sensors.
This is a device for estimating the arrival direction of a signal. For simplicity, it is assumed that the N sensors 31 1 to 31 N of the line array 31 are arranged on the y-axis as shown in FIG. Further, the arrival direction of the signal is represented by an angle θ from the normal direction of the line array 31.
【0004】ラインアレイ31の各センサ311〜31N
には、周波数分割器32が接続され、その出力側には、
周波数領域の信号の整相を行う周波数領域整相器33が
接続されている。更に、周波数領域整相器33には、周
波数分析器34が接続され、その周波数分析器34に
は、パワ−算出器35を介して、ピ−ク検出器36が接
続され、更にピ−ク検出器36の出力側には出力端子3
7が接続されている。Each sensor 31 1 to 31 N of the line array 31
Is connected to a frequency divider 32, and on the output side,
A frequency domain phasing device 33 for phasing a signal in the frequency domain is connected. Further, a frequency analyzer 34 is connected to the frequency domain phase adjuster 33, and a peak detector 36 is connected to the frequency analyzer 34 via a power calculator 35. An output terminal 3 is provided on the output side of the detector 36.
7 is connected.
【0005】次に、この方位推定装置を用いた信号源の
方位推定方法を説明する。例えば、海水中の信号源から
の信号は、ラインアレイ31に受波される。時刻tにラ
インアレイ31に入力された信号は、センサ311〜3
1N単位に電気信号に変換された後、周波数分割器32
に入力される。周波数分割器32では、ラインアレイ3
1からの(1)式で示される出力XNext, a method for estimating the azimuth of a signal source using the azimuth estimating apparatus will be described. For example, a signal from a signal source in seawater is received by the line array 31. The signals input to the line array 31 at the time t are sensors 31 1 to 3 1
After being converted into electric signals in units of 1 N , the frequency divider 32
Is input to In the frequency divider 32, the line array 3
Output X from Equation 1 as shown in equation (1)
【数1】 に対して、各センサ出力端子ごとに周波数分析を行い、
それぞれをK個の帯域に分割する。そして、(2)式で
示されるその出力Y(Equation 1) For each sensor output terminal.
Each is divided into K bands. Then, the output Y expressed by the equation (2)
【数2】 を周波数領域整相器33に出力する。(Equation 2) To the frequency domain phase adjuster 33.
【0006】周波数領域整相器33は、位相補償回路及
び加算器等で構成されており、一周波数帯域当りθ1〜
θM方向のM個のビ−ムを形成する装置である。θ1方向
のビ−ムを形成するためには、まず位相補償回路によっ
て、周波数分割器32から出力Yの各要素yn(fj)に、
(3)式の位相補償量 Φn(θ1,fj)=exp[(-j2πfjynsinθ1)/c] (j=1〜
K,n=1〜N)・・・・(3) 但し、c:音速、j:虚数単
位、π:円周率 の乗算を行い、その後、Yに空間シェ−ディング係数ベ
クトル Ω=[ω1,・・・,ωN] を乗ずる。すなわち(4)式[0006] frequency domain phasing unit 33 is constituted by a phase compensation circuit and the adder, etc., one frequency band per theta 1 ~
theta M direction of the M-bi - a device for forming a beam. theta 1 direction bi - to form a beam, the first by a phase compensation circuit, each element y n output Y from the frequency divider 32 (f j),
(3) the phase compensation amount Φ n (θ 1, f j ) of the formula = exp [(- j2πf j y n sinθ 1) / c] (j = 1~
K, n = 1 to N) (3) where c: sound velocity, j: imaginary unit, π: pi, and then Y is a space shading coefficient vector Ω = [ω 1 ,..., Ω N ]. That is, equation (4)
【数3】 の算出を行い、これをθ1方向ビ−ム出力とする。θ2〜
θM方向に対しても同様の処理を行い、z(θ2)〜z(θM)
を計算する。そして、(5)式で示されるその結果(Equation 3) Perform the calculation, which theta 1 direction bi - a beam output. θ 2 ~
The same processing is performed for the θ M direction, and z (θ 2 ) to z (θ M )
Is calculated. Then, the result shown by equation (5)
【数4】 を周波数分析器34に出力する。(Equation 4) Is output to the frequency analyzer 34.
【0007】周波数分析器34は、周波数領域整相器3
3からの出力 Z=[z(θi,fj)]i=1〜M,j=1〜k の各要素に対して、更にFFT(高速フーリエ変換)に
よる周波数分析を行い、一周波数帯域あたりP個の周波
数ビンを抽出する。そして、その出力 S=[s(θi,fjp)]i=1〜M,j=1〜K,p=1〜P をパワ−算出器35に出力する。The frequency analyzer 34 is a frequency domain phase adjuster 3
3. Output from Z = [z (θ i , f j )] For each element of i = 1 to M, j = 1 to k , a frequency analysis by FFT (Fast Fourier Transform) is further performed to obtain one frequency band. Per P frequency bins. Then, the output S = [s (θ i , f jp )] i = 1 to M, j = 1 to K, and p = 1 to P are output to the power calculator 35.
【0008】パワ−算出器35は、周波数分析器34の
出力に対し、複素領域でのパワ−算出処理、即ち絶対値
の2乗算出処理を行い、場合によっては、その後、信号
を時間的に滑らかに検出しやすくするための積分処理を
行い、その結果 |S|2=[|s(θi,fjp)|2]i=1〜M,j=1〜k,p=1〜P をピ−ク検出器36に出力する。[0008] The power calculator 35 is provided by the frequency analyzer 34.
For the output, the power calculation process in the complex domain, that is, the absolute value
, And in some cases, the signal
Integration to make it easier to detect
And the result | S |Two= [| S (θi, fjp) |Two]i = 1-M, j = 1-k, p = 1-P Is output to the peak detector 36.
【0009】ピーク検出器36は、パワー算出器35の
出力 |S|2=[|s(θi,fjp)|2]i=1〜M,j=1〜k,p=1〜P に対して方位θ及び周波数fの2次元でのピーク検出処
理を行い、そのピーク(局所的な最大値)となる要素|s
(θi,fjp)|2を検出し、その方位θi及び周波数fjpを出
力端子37に出力する。この2段の周波数分析方法、即
ち、1段目ではK個の帯域に応じた粗い分析幅で周波数
分析を行ない、2段目の周波数分析ではK×P個の周波
数ビンに対応した細かい分析幅の周波数を行う方法によ
れば、センサ出力を2つに分岐して周波数領域整相のた
めの周波数分析と周波数推定のための周波数分析を別々
に行う場合に比べて、2段目の周波数分析を行う際に周
波数分析を行う周波数領域とビームとを制限することが
できる等のため、演算量は大幅に少なくなる。そして、
信号源の狭帯域信号の周波数分析は、2段目の周波数分
析器でのK×P個の帯域(分析幅)に対応した精度で行
われ、周波数f j に比べてP倍細かい精度の推定周波数f
jp を出力することができる。 The peak detector 36 outputs the output | S | 2 = [| s (θ i , f jp ) | 2 ] i = 1 to M, j = 1 to k, p = 1 to P of the power calculator 35. Is subjected to two-dimensional peak detection processing of the azimuth θ and the frequency f, and the element | s that becomes the peak (local maximum value)
(θ i , f jp ) | 2 and outputs the azimuth θ i and the frequency f jp to the output terminal 37. This two-stage frequency analysis method,
In the first stage, the frequency is determined by the coarse analysis width corresponding to the K bands.
The analysis is performed, and K × P frequency
A method that performs a frequency with a fine analysis width corresponding to several bins
Then, the sensor output is branched into two,
Frequency analysis for frequency estimation and frequency analysis for frequency estimation separately
When performing the second-stage frequency analysis,
It is possible to limit the frequency domain and beam for wave number analysis.
Since it is possible, the amount of calculation is significantly reduced. And
The frequency analysis of the narrow band signal of the signal source
With accuracy corresponding to K × P bands (analysis width) in analyzer
Estimated frequency f, which is P times finer than frequency f j
jp can be output.
【0010】[0010]
【発明が解決しようとする課題】しかしながら、上記構
成の装置では、次のような課題があった。例えば周波数
fの狭帯域信号の方位を推定する場合を考える。θ1方
向への周波数領域整相において、周波数fの信号に対す
るセンサnにおける位相補償量Φn(θi,f)は,理論的に
は、次の(6)式で求まる。 Φn(θi,f)=exp[(-j2πfynsinθi)/c] ・・・・・・・・・・・・・・・・・・・・・(6)However, the apparatus having the above configuration has the following problems. For example, consider the case of estimating the azimuth of a narrowband signal of frequency f. In the frequency domain the phasing of the theta 1 direction, the phase compensation amount in the sensor n to the signal of the frequency f Φ n (θ i, f ) is theoretically determined by the following equation (6). Φ n (θ i , f) = exp [(-j2πfy n sin θ i ) / c] (6)
【0011】ところが、図3に示す周波数分割器32
は、比較的粗い幅の帯域への分割を行う。このために、
周波数領域整相器33で算出される位相補償量は、
(6)式における狭帯域信号周波数fを含む帯域の中心
周波数fjに置き換えた Φn(θi,fj)=exp[(-j2πfjynsinθi)/c] となっている。従って、位相補償量Φ n (θi,fj)に
は誤差が生じ、更にそれを用いて推定された方位も誤差
を含んだものとなる。特に信号の周波数が帯域分割の境
界値付近に存在する場合においては、その推定された方
位の誤差は最も大きなものとなる。However, the frequency divider 32 shown in FIG.
Performs a division into bands of relatively coarse width. For this,
The phase compensation amount calculated by the frequency domain phase adjuster 33 is
(6) [Phi replaced with the center frequency f j of the band including a narrowband signal frequency f in formula n (θ i, f j) = exp [(- j2πf j y n sinθ i) / c] and has. Therefore, an error occurs in the phase compensation amount Φ n (θ i , f j ), and the azimuth estimated using the error also includes the error. In particular, when the frequency of the signal exists near the boundary value of the band division, the error of the estimated azimuth becomes the largest.
【0012】本発明は上述の点に鑑みてなされたもの
で、上記問題点即ち、周波数分割が比較的粗い幅で実現
されているために周波数領域整相で用いられている位相
補償量の精度が悪くなることで、精度の良い方位推定が
行えないという問題点を解決し、高精度で方位の推定を
することができる狭帯域信号による方位推定方法を提供
するものである。[0012] The present invention has been made in view of the above, the problem that is, the phase compensation amount of accuracy that is used in the frequency domain phasing to frequency division is implemented with a relatively coarse width by deteriorates, to solve the problem that can not be carried out accurate direction estimation, the direction estimation with high accuracy
And an azimuth estimating method using a narrow band signal .
【0013】[0013]
【課題を解決するための手段】上記課題を解決するため
本発明は、図1に示すように狭帯域信号の到来方向推定
方法において、各ビ−ムの出力の狭帯域分析結果から該
狭帯域信号周波数と前記の分割された周波数帯域のうち
該狭帯域信号周波数を含む帯域の周波数領域整相に用い
た代表周波数を用いて推定された到来方向の補正を行う
手段を設けたものである。以下、その原理を説明する。According to the present invention, there is provided a method for estimating a direction of arrival of a narrow band signal, as shown in FIG. 1, based on a narrow band analysis result of the output of each beam. Means is provided for correcting the estimated direction of arrival using the signal frequency and the representative frequency used in the frequency domain phasing of the band including the narrow band signal frequency among the divided frequency bands. Hereinafter, the principle will be described.
【0014】従来技術の問題点で説明したように、θ方
向への周波数領域整相において、周波数fの信号に対す
るセンサnにおける位相補償φn(θ,f)は、理論的に
は、次の(7)式によて求まる。 Φn(θ,f)=exp[(-j2πfynsinθ)/c] ・・・・・・・・・・・・・・・・・・・・・・(7)As described in the prior art, in the frequency domain phasing in the θ direction, the phase compensation φ n (θ, f) in the sensor n for the signal of the frequency f is theoretically expressed as follows. It is obtained by equation (7). Φ n (θ, f) = exp [(-j2πfy n sinθ) / c] (7)
【0015】ところが、現実の装置においては、周波数
fが、適当な代表周波数fjに置き換えられた Φn(θ,fj)=exp[(-j2πfjynsinθ)/c] ・・・・・・・・・・・・・・・・・・・・(8) に近似される。(7)式と(8)式との比較から、ピ−
ク検出器で求められた信号到来方向θと実際の方位θ´
との間には、次の(9)式の関係が成り立つことが分か
る。 fsinθ´=fjsinθ ・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・(9)However, in an actual device, the frequency f is replaced by an appropriate representative frequency f j Φ n (θ, f j ) = exp [(− j2πf j n sin θ) / c] (8) From the comparison between the equations (7) and (8), the peak
Signal arrival direction θ and actual azimuth θ '
It can be seen that the following relationship (9) is established between. fsinθ '= f j sinθ (9)
【0016】従って、ピ−ク検出方位θ及び探知周波数
fを用いて、次の(10)式により方位補正を行えばバ
イアス誤差のない方位が得られることになる。 θ´=arc sin((fj/f)sinθ) ・・・・・・・・・・・・・・・・・・・・・・・・・・(10) 以上詳細に述べた方位補正手段をピ−ク検出器6の後に
設けた。Therefore, if the azimuth is corrected by the following equation (10) using the peak detection azimuth θ and the detection frequency f, an azimuth without a bias error can be obtained. θ ′ = arc sin ((f j / f) sin θ) (10) Azimuth correction described in detail above Means were provided after the peak detector 6.
【0017】[0017]
【作用】本発明によれば、上記動作原理で詳細に説明し
たように上記補正手段7を設けることによって周波数領
域整相器3に周波数分割の代表周波数を用いるために生
じる推定方位の誤差の影響を取り除くことが出来る。According to the present invention, as described in detail in the above operation principle, the influence of the error of the estimated azimuth caused by using the representative frequency of the frequency division for the frequency domain phase adjuster 3 by providing the correction means 7 is provided. Can be removed.
【0018】[0018]
【実施例】以下本発明の一実施例を図面に基づいて詳細
に説明する。図1は、本発明の実施例を示す一例で、狭
帯域信号の方位推定装置の機能ブロック図である。この
方位推定装置は、例えば従来と同様なセンサからなるラ
インアレイ1を用い、狭帯域信号の到来方向を推定する
装置である。簡単の為、ラインアレイ1のN個のセンサ
11〜1Nはy軸上に配列されているものとする。又、信
号の到来方向をラインアレイ1の法線方向からの角度θ
で表す。DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS One embodiment of the present invention will be described below in detail with reference to the drawings. FIG. 1 is an example showing an embodiment of the present invention, and is a functional block diagram of an azimuth estimating apparatus for a narrowband signal. This azimuth estimating apparatus is an apparatus for estimating an arrival direction of a narrowband signal using, for example, a line array 1 composed of sensors similar to a conventional sensor. For simplicity, it is assumed that the N sensors 11 to 1N of the line array 1 are arranged on the y-axis. Further, the arrival direction of the signal is defined as an angle θ from the normal direction of the line array 1.
Expressed by
【0019】ラインアレイ1の各センサ11〜1Nは、周
波数分割器2が接続され、その出力側には、周波数領域
の信号の整相を行う周波数領域整相器3が接続されてい
る。更に、周波数領域整相器3には、周波数分析器4が
接続され、その周波数分析器4には、パワ−算出器5を
介して、ピ−ク検出器6が接続されている。更にピ−ク
検出器6の出力側には方位補正手段7を設け、方位補正
手段7の出力側には出力端子8が接続されている。Each of the sensors 1 1 to 1 N of the line array 1 is connected to a frequency divider 2, and the output side thereof is connected to a frequency-domain phase adjuster 3 for phasing a frequency-domain signal. . Further, a frequency analyzer 4 is connected to the frequency domain phase adjuster 3, and a peak detector 6 is connected to the frequency analyzer 4 via a power calculator 5. Further, an azimuth correction means 7 is provided on the output side of the peak detector 6, and an output terminal 8 is connected to the output side of the azimuth correction means 7.
【0020】図2は、図1中の方位補正手段7の動作手
順を示す図である。sin関数計算手段21、除算手段2
2、乗算手段23、arc sin関数計算手段24で構成さ
れている。次に、以上のように構成される方位推定装置
を用いた方位推定方法を説明する。FIG. 2 is a diagram showing an operation procedure of the azimuth correction means 7 in FIG. sin function calculating means 21, dividing means 2
2. Multiplication means 23 and arc sin function calculation means 24. Next, an azimuth estimation method using the azimuth estimation device configured as described above will be described.
【0021】時刻tにおいてラインアレイ1に入力され
た信号は、センサ11〜1N単位に電気信号に変換された
後、周波数分割器2に入力される。周波数分割器2で
は、ラインアレイ1からの(11)次式で示される出力
XThe time input signal to the line array 1 at t is converted into an electric signal to the sensor 1 1 to 1 N units, is input to the frequency divider 2. In the frequency divider 2, the output X from the line array 1 represented by (11)
【数5】 に対して、各センサ出力ごとに周波数分析を行い、それ
ぞれK個の帯域に分割する。そして、(12)で示され
るその出力(Equation 5) Is subjected to frequency analysis for each sensor output, and is divided into K bands. And its output shown in (12)
【数6】 を周波数領域整相器4に出力する。(Equation 6) To the frequency domain phase adjuster 4.
【0022】周波数領域位相器4は、位相補償回路及び
加算器等で構成されており、一周波数帯域当りθ1〜θM
方向のM個のビ−ムを形成する装置である。θ1方向の
ビ−ムを形成するためには、まず周波数分割器2からの
出力Yに対して、まず位相補償回路によって、各周波数
帯域fiごとに位相補償量 Φn(θ1,fj)=exp[(-j2πfjynsinθ1)/c] (j=1〜K) ・・・・・・・・(13) 但し、c:音速、j:虚数単位、π:円周率 の乗算を行い、その後、Yに式で示されるシェ−ンディ
ング係数ベクトル Ω=[ω1,・・・,ωN] を乗ずる。すなわちThe frequency domain phase shifter 4 is constituted by a phase compensation circuit and the adder, etc., one frequency band per theta 1 through? M
This is an apparatus for forming M beams in one direction. theta 1 direction bi - to form a beam, first the output Y from the frequency divider 2, firstly by the phase compensation circuit, the frequency band f phase compensation amount for each i Φ n (θ 1, f j) = exp [(- j2πf j y n sinθ 1) / c] (j = 1~K) ········ (13) where, c: sound velocity, j: imaginary unit, [pi: circle Then, Y is multiplied by a shading coefficient vector Ω = [ω 1 ,..., Ω N ] shown in the equation. Ie
【数7】 を算出し、これをθ1方向のビ−ム出力とする。θ2〜θ
M方向に対しても同様の処理を行い、z(θ2)〜z
(θM)を算出する。そして、(15)次式で示されるそ
の結果Z(Equation 7) It is calculated and this theta 1 direction bi - a beam output. θ 2 to θ
The same processing is performed for the M direction, and z (θ 2 ) to z
(θ M ) is calculated. Then, (15) the result Z expressed by the following equation
【数8】 を周波数分析器4に出力する。(Equation 8) Is output to the frequency analyzer 4.
【0023】周波数分析器4は、周波数領域整相器3か
らの(16)式で示される出力 Z=[z(θi,fj)]i=1〜M,j=1〜k ・・・・・・・・・・・・・・・・・(16) の各要素に対して、更にFFTによる周波数解析を行
い、一周波数帯域あたりP個の周波数ビンを抽出する。
そして、その出力 S=[s(θi,fjp)]i=1〜M,j=1〜K,p=1〜P ・・・・・・・・・・・・・・・・・(17) をパワー算出器5に出力する。ここで、f j は、1段目の
周波数分析の代表周波数に相当し、f jp は2段目の周波
数の代表周波数に相当し、f j よりも精度の高いものであ
る。なお、2段目の周波数分析は、1段目の周波数分析
で帯域分割した信号が対象となるため、必然的に長い分
析時間をとり、例えば、1段目の周波数分析を1/10
0秒毎に時間フレームで行い、2段目の周波数分析を5
秒のフレームで行う。 The frequency analyzer 4 outputs an output Z = [z (θ i , f j )] i = 1 to M, j = 1 to k ... ... For each element of (16), a frequency analysis is further performed by FFT to extract P frequency bins per frequency band.
And its output S = [s (θ i , f jp )] i = 1 to M, j = 1 to K, p = 1 to P ... (17) is output to the power calculator 5. Here, f j is the first stage
This corresponds to the representative frequency of frequency analysis, and fjp is the frequency of the second stage.
Frequency, which is more accurate than f j
You. The frequency analysis in the second stage is the frequency analysis in the first stage.
Because the signal divided in
Take the analysis time, for example, the first-stage frequency analysis
The frequency analysis of the second stage is performed every 5 seconds in a time frame.
Perform in seconds frames.
【0024】パワ−算出器5は、周波数分析器4の出力
に対し、複素領域でのパワ−算出処理、即ち絶対値の2
乗算出処理を行い、場合によっては、その後、信号を時
間的に滑らかに検出しやすくするために積分等の処理を
行い、その結果 |S|2=[|s(θi,fjp)|2]i=1〜M,j=1〜k,p=1〜P ・・・・・・・・・・・(18) をピ−ク検出器6に出力する。The power calculator 5 performs a power calculation process in a complex domain on the output of the frequency analyzer 4, that is, 2 of the absolute value.
A multiplication calculation process is performed, and in some cases, a process such as integration is then performed to make it easier to detect the signal temporally smoothly. As a result, | S | 2 = [| s (θ i , f jp ) | 2 ] i = 1 to M, j = 1 to k, p = 1 to P (18) are output to the peak detector 6.
【0025】ピ−ク検出器6は、パワ−算出器5の(1
9)式で示される出力 |S|2=[|s(θi,fjp)|2]i=1〜M,j=1〜k,p=1〜P ・・・・・・・・・・・(19) に対して方位θ及び周波数fの2次元でのピ−ク検出処
理を行い、そのピ−クとなる要素|s(θi,fjp)|2の方向
θi及び周波数fjp 、更に、周波数fjpを含む周波数領域
整相に用いた分割帯域の代表周波数fjを方位補正手段7
に出力する。The peak detector 6 is provided by (1) of the power calculator 5.
9) Output | S | 2 = [| s (θ i , f jp ) | 2 ] i = 1 to M, j = 1 to k, p = 1 to P ... peak of a two-dimensional orientation theta and the frequency f with respect ... (19) - perform click detection process, the peak - click become elements | s (θ i, f jp ) | 2 direction theta i and The frequency f jp and the representative frequency f j of the divided band used for the frequency domain phasing including the frequency f jp are corrected by the azimuth correction means 7.
Output to
【0026】方位補正手段7を図2で説明する。ピ−ク
検出器6の出力のうちθiをsin関数計算手段21に、f
jp及びfjを除算手段22に入力する。sin関数計算手段
21はθiのsin関数値を算出し、その結果を乗算手段2
3に送る。除算手段22は、fj/fjpを算出し、その結果
を乗算手段23に送る。乗算手段23は、sin関数計算
手段21の出力sinθiと除算手段22の出力fj/fjpの積
を算出し、その結果(fi/fjp)sinθiをarc sin関数計算
手段24に送る。arc sin関数計算手段24は、乗算手
段23の出力(fj/fjp)sinθiのarc sin関数値を算出
し、その結果arc sin(fj/fjp)sinθiを出力する。方位
補正手段7は、上記手段によって算出されたθを信号源
の方位推定値として出力端子8に出力する。The azimuth correcting means 7 will be described with reference to FIG. Among the outputs of the peak detector 6, θ i is given to the sin function calculating means 21 and f
jp and f j are input to the dividing means 22. The sin function calculating means 21 calculates the sin function value of θ i and multiplies the result by the multiplying means 2
Send to 3. The divider 22 calculates f j / f jp and sends the result to the multiplier 23. The multiplication means 23 calculates the product of the output sin θ i of the sin function calculation means 21 and the output f j / f jp of the division means 22, and outputs the result (f i / f jp ) sin θ i to the arc sin function calculation means 24. send. The arc sin function calculating means 24 calculates an arc sin function value of the output (f j / f jp ) sin θ i of the multiplying means 23, and outputs arc sin (f j / f jp ) sin θ i as a result. The azimuth correcting means 7 outputs the θ calculated by the above means to the output terminal 8 as the azimuth estimated value of the signal source.
【0027】[0027]
【発明の効果】以上、詳細に説明したように本発明によ
れば、下記のような優れた効果が期待できる。ピーク検
出処理によって推定した推定方位θ i を、ピーク検出処
理によって推定した精度の高い推定周波数f jp と該推定
周波数f jp に対応した第1の周波数分析における代表周
波数f j との比を用いて修正するので、周波数領域整相
に周波数分割の代表周波数を用いるために生じる推定方
位の誤差の影響を取り除くことが出来、目標とする信号
源の方向を高精度に推定出来る。Effect of the Invention] According to the present invention as described in detail, it can be expected excellent effects as follows. Peak detection
The estimated orientation theta i estimated by output processing, the peak detection processing
Frequency f jp with high accuracy estimated by
Representative frequency in the first frequency analysis corresponding to frequency fjp
Since corrected using the ratio of the wave number f j, you can eliminate the influence of the estimated heading error that occur in order to use the typical frequency of the frequency division to the frequency domain phasing, the direction of the signal source for the target precision Can be estimated.
【図1】本発明の実施例の狭帯域信号による方位推定装
置の機能ブロック図である。FIG. 1 is a functional block diagram of an azimuth estimating apparatus based on a narrowband signal according to an embodiment of the present invention.
【図2】方位補正手段の手順を示すブロック図である。FIG. 2 is a block diagram illustrating a procedure of an azimuth correction unit.
【図3】従来の狭帯域信号による方位推定装置の機能ブ
ロック図である。FIG. 3 is a functional block diagram of a conventional azimuth estimating apparatus based on a narrowband signal.
【図4】ラインアレイを説明する図である。FIG. 4 is a diagram illustrating a line array.
1 ラインアレイ 2 周波数分析器 3 周波数領域整相器 4 周波数分析器 5 パワ−算出器 6 ピ−ク検出器 7 方位補正手段 DESCRIPTION OF SYMBOLS 1 Line array 2 Frequency analyzer 3 Frequency domain phase adjuster 4 Frequency analyzer 5 Power calculator 6 Peak detector 7 Direction correction means
───────────────────────────────────────────────────── フロントページの続き (56)参考文献 特開 昭59−104578(JP,A) 特公 平6−93014(JP,B2) 特公 平6−105289(JP,B2) ──────────────────────────────────────────────────続 き Continuation of the front page (56) References JP-A-59-104578 (JP, A) JP 6-93014 (JP, B2) JP 6-105289 (JP, B2)
Claims (1)
析により複数の周波数帯域に分割し、該各周波数帯域に
対し周波数領域整相により複数のビームを形成し、該各
ビームの各周波数帯域の出力を更に第2の周波数分析に
より周波数分析し、該分析結果に対して方位と周波数と
の2次元でのピーク検出処理を行って、狭帯域信号の到
来方向の方位θ i と狭帯域信号の周波数f jp を推定し、 前記ピーク検出処理によって推定した前記推定方位θ i
を、前記ピーク検出処理によって推定した前記推定周波
数f jp と該推定周波数f jp に対応した前記第1の周波数
分析における 代表周波数f j との比を用いて、修正する
ことを特徴とする狭帯域信号による方位推定方法。1. A method for receiving signals from a plurality of sensors.FirstFrequency
Analysis into multiple frequency bands,TheeachfrequencyTo the band
On the other hand, multiple beams are formed by frequency domain phasing,Theeach
BeamOf each frequency bandoutputTo the second frequency analysis
More frequency analysis, the direction and frequency with respect to the analysis result
Perform the two-dimensional peak detection process ofNarrowband signal arrival
Coming directionOrientation θ i And the frequency f of the narrowband signal jp EstimateAnd The estimated orientation θ estimated by the peak detection processing i
Is the estimated frequency estimated by the peak detection process.
Number f jp And the estimated frequency f jp The first frequency corresponding to
In the analysis Representative frequencyf j Using the ratio of
Characterized byBy narrowband signalBearing estimation method.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP3250429A JP2717331B2 (en) | 1991-09-03 | 1991-09-03 | Direction estimation method using narrowband signal |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP3250429A JP2717331B2 (en) | 1991-09-03 | 1991-09-03 | Direction estimation method using narrowband signal |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH0560851A JPH0560851A (en) | 1993-03-12 |
JP2717331B2 true JP2717331B2 (en) | 1998-02-18 |
Family
ID=17207754
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP3250429A Expired - Fee Related JP2717331B2 (en) | 1991-09-03 | 1991-09-03 | Direction estimation method using narrowband signal |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2717331B2 (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5691753B2 (en) * | 2011-04-04 | 2015-04-01 | 沖電気工業株式会社 | Wideband signal processing system |
-
1991
- 1991-09-03 JP JP3250429A patent/JP2717331B2/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JPH0560851A (en) | 1993-03-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wong et al. | A combined SAR Doppler centroid estimation scheme based upon signal phase | |
US5874916A (en) | Frequency selective TDOA/FDOA cross-correlation | |
CA2054995C (en) | System and method for compensation of in-phase and quadrature phase and gain imbalance | |
KR960011786B1 (en) | Fast phase difference autofocusing method | |
US5495256A (en) | Method of efficiently computing maximum likelihood bearing estimator | |
US20060248138A1 (en) | Signal processing device for computing phase difference between alternating current signals | |
JP4681599B2 (en) | Direction of arrival estimation device | |
JP5110531B2 (en) | FM-CW radar equipment | |
JP2717331B2 (en) | Direction estimation method using narrowband signal | |
US5610612A (en) | Method for maximum likelihood estimations of bearings | |
Maskell et al. | Adaptive subsample delay estimation using a modified quadrature phase detector | |
JP2648110B2 (en) | Signal detection method and device | |
JPS6349193B2 (en) | ||
US20200142020A1 (en) | Method and device for estimating an angle of arrival of an incident radio signal | |
Zamlutti et al. | Incoherent scatter multiple‐pulse measurements at Arecibo | |
JP3269965B2 (en) | Orientation detection device | |
JPH05188145A (en) | Method and apparatus for extracting phase error caused by atmosphere in backward scattering signal of interference-type image radar system mounted on carrier out of raw radar data | |
JP3916952B2 (en) | Synthetic aperture radar apparatus and synthetic aperture radar image reproduction method | |
JP3247549B2 (en) | Monopulse radar device | |
JPH0882664A (en) | Estimation method of signal arrival direction | |
JP2003232849A (en) | Adaptive phasing method | |
JPH0138270B2 (en) | ||
JPH0257993A (en) | Radar target detecting apparatus | |
JP2758871B2 (en) | Frequency analysis method by maximum entropy method | |
JP2000101327A (en) | Radio wave receiver with improved antenna resolution |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20081114 Year of fee payment: 11 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20081114 Year of fee payment: 11 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091114 Year of fee payment: 12 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101114 Year of fee payment: 13 |
|
LAPS | Cancellation because of no payment of annual fees |