JP2648110B2 - Signal detection method and device - Google Patents
Signal detection method and deviceInfo
- Publication number
- JP2648110B2 JP2648110B2 JP29649494A JP29649494A JP2648110B2 JP 2648110 B2 JP2648110 B2 JP 2648110B2 JP 29649494 A JP29649494 A JP 29649494A JP 29649494 A JP29649494 A JP 29649494A JP 2648110 B2 JP2648110 B2 JP 2648110B2
- Authority
- JP
- Japan
- Prior art keywords
- signal
- frequency
- vector
- unit
- power spectrum
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Lifetime
Links
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- 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 signal detection method and apparatus for processing a directional sonobuy received signal and detecting a wideband underwater acoustic signal.
【0002】[0002]
【従来の技術】従来より、この種の信号検出装置とし
て、図3に示すように、周波数分析部11と、これから
の出力を受ける検出部12とを有している。周波数分析
部11は、指向性ソノブイの受波信号x(t)を受けて
周波数分析を行い、パワースペクトル値P(f)を一定
時間毎に出力する。検出部12は、パワースペクトル値
P(f)を受けて、この値がスレッショルドレベルを越
える信号を検出する。2. Description of the Related Art Conventionally, as this type of signal detection device, as shown in FIG. 3, a frequency analysis unit 11 and a detection unit 12 for receiving an output from the frequency analysis unit 11 are provided. The frequency analysis unit 11 receives the received signal x (t) of the directional sonobuoy, performs frequency analysis, and outputs a power spectrum value P (f) at regular intervals. Upon receiving the power spectrum value P (f), the detecting unit 12 detects a signal whose value exceeds the threshold level.
【0003】[0003]
【発明が解決しようとする課題】しかし、この従来の信
号検出装置は、周波数分析を行うことで、狭帯域分析に
より広帯域にわたって分布する雑音のレベルを低下させ
ているから、広帯域性の信号のレベルも低下するために
S/Nの効果が小さく信号検出が困難であるという問題
がある。However, in this conventional signal detection device, the frequency analysis is performed to reduce the level of noise distributed over a wide band by narrow-band analysis. Therefore, there is a problem that the effect of S / N is small and signal detection is difficult.
【0004】本発明の目的は、広帯域性の水中音響信号
の検出が容易である信号検出方法及び装置を提供するこ
とにある。[0004] It is an object of the present invention to provide a signal detection method and apparatus capable of easily detecting a wideband underwater acoustic signal.
【0005】[0005]
【課題を解決するための手段】本発明は、前記課題を解
決するために、指向性ソノブイの受波信号から周波数毎
の信号方位値とパワースペクトル値を生成し、これらの
信号方位とパワースペクトル値により周波数毎の信号ベ
クトルを合成し、これらの信号ベクトルを全周波数で加
算して全周波数の和信号ベクトルを求め、次にこの和信
号ベクトルのスレッショルドレベルを越える信号を検出
することを特徴とする。According to the present invention, a signal azimuth value and a power spectrum value for each frequency are generated from a received signal of a directional sonobuoy. A signal vector for each frequency is synthesized according to the value, a sum signal vector of all frequencies is obtained by adding these signal vectors at all frequencies, and then a signal exceeding a threshold level of the sum signal vector is detected. I do.
【0006】また、本発明は、指向性ソノブイの受波信
号を周波数分析して周波数毎の位相スペクトル値とパワ
ースペクトル値を一定時間毎に出力する周波数分析部
と、この周波数分析部から位相スペクトル値とパワース
ペクトル値を受けてその時間の信号に対応した周波数毎
の信号到来方位を算出する方位算出部と、前記周波数分
析部からのパワースペクトル値と前記方位算出部からの
信号到来方位により形成される信号ベクトルを加算して
全周波数の和信号ベクトルを求めるベクトル加算部と、
このベクトル加算部からの和信号ベクトルのスレッショ
ルドレベルを越える信号を検出する検出部とを有するこ
とを特徴とする。Further, the present invention provides a frequency analysis unit for frequency-analyzing a received signal of a directional sonobuoy and outputting a phase spectrum value and a power spectrum value for each frequency at fixed time intervals, and a phase spectrum from the frequency analysis unit. A direction calculation unit that receives a value and a power spectrum value to calculate a signal arrival direction for each frequency corresponding to a signal at that time, and a power spectrum value from the frequency analysis unit and a signal arrival direction from the direction calculation unit. A vector addition unit for adding the signal vectors to be added to obtain a sum signal vector of all frequencies,
A detection unit for detecting a signal exceeding a threshold level of the sum signal vector from the vector addition unit.
【0007】[0007]
【実施例】次に、本発明の実施例を図面に基いて詳細に
説明する。Next, an embodiment of the present invention will be described in detail with reference to the drawings.
【0008】図1は、本発明の1実施例を示すブロック
図である。本発明の信号検出装置は、周波数分析部1
と、この周波数分析部1からの出力を受ける方位算出部
2と、前記周波数分析部1および方位算出部2からの出
力を受けるベクトル加算部3と、このベクトル加算部3
の出力を受ける検出部4とを有する。FIG. 1 is a block diagram showing one embodiment of the present invention. The signal detection device according to the present invention includes a frequency analysis unit 1
An azimuth calculation unit 2 receiving an output from the frequency analysis unit 1, a vector addition unit 3 receiving an output from the frequency analysis unit 1 and the azimuth calculation unit 2, and a vector addition unit 3
And a detection unit 4 that receives the output of
【0009】前記周波数分析部1は、指向性ソノブイの
受波信号x(t)を受けて周波数分析を行い、パワース
ペクトル値P(f)および位相スペクトル値θ(f)を
一定時間毎に出力する。前記方位算出部2は、周波数分
析部1からのパワースペクトル値P(f)および位相ス
ペクトル値θ(f)に基いて、その時間の信号に対応し
た周波数毎の信号到来方位B(f)を算出する。The frequency analysis unit 1 receives a received signal x (t) of a directional sonobuoy and performs frequency analysis, and outputs a power spectrum value P (f) and a phase spectrum value θ (f) at regular time intervals. I do. Based on the power spectrum value P (f) and the phase spectrum value θ (f) from the frequency analysis unit 1, the direction calculation unit 2 calculates the signal arrival direction B (f) for each frequency corresponding to the signal at that time. calculate.
【0010】前記ベクトル加算部3は、周波数分析部1
からのパワースペクトル値P(f)と方位算出部2から
の信号到来方位B(f)とを受けて、これらにより形成
される周波数毎の信号ベクトルV(f)を加算し、次の
数1で表される全周波数の和ベクトルVを求めて出力す
る。The vector adder 3 includes a frequency analyzer 1
, And a signal arrival direction B (f) from the direction calculation unit 2, add a signal vector V (f) for each frequency formed by these, and obtain the following equation (1). Is obtained and output.
【0011】[0011]
【数1】 前記信号ベクトルV(f)は、ベクトルの角度がB
(f)であり、ベクトルの大きさがパワースペクトル値
P(f)であるベクトルである。すなわち、この信号ベ
クトルV(f)を直交する2つの方向成分VN-S (f)
とVE-W (f)に分解すると、それぞれ次の数2および
数3で表される。(Equation 1) The signal vector V (f) has a vector angle B
(F), where the magnitude of the vector is the power spectrum value P (f). That is, the signal vector V (f) is transformed into two orthogonal components V NS (f)
And V EW (f) are expressed by the following equations 2 and 3, respectively.
【0012】[0012]
【数2】 (Equation 2)
【0013】[0013]
【数3】 前記検出部4は、全周波数の和ベクトルVを受けて、ス
レッショルドレベルを越える信号を検出する。(Equation 3) The detection unit 4 receives the sum vector V of all frequencies and detects a signal exceeding a threshold level.
【0014】次に、前記ベクトル加算部3の処理を図2
を参照して説明する。Next, the processing of the vector adder 3 will be described with reference to FIG.
This will be described with reference to FIG.
【0015】前記ベクトル加算部3で受けたパワースペ
クトル値P(f)と方位算出部2からの信号到来方位B
(f)とにより形成される周波数毎の信号ベクトルV
(f)は、信号成分Vs (f)と雑音成分Vn (f)の
和になっている。信号成分Vs(f)は一定の方位から
到来する成分であり、雑音成分Vn (f)は全方位から
到来する等方性の成分であるとする。The power spectrum value P (f) received by the vector addition unit 3 and the signal arrival direction B from the direction calculation unit 2
(F) the signal vector V for each frequency formed by
(F) is the sum of the signal component V s (f) and the noise component V n (f). It is assumed that the signal component V s (f) is a component coming from a fixed direction, and the noise component V n (f) is an isotropic component coming from all directions.
【0016】信号成分Vs (f)は全周波数で信号成分
方位Bs (f)が等しい、すなわちBs (f1 )=Bs
(f2 )=…=Bs (fK )=Bs となるから、全周波
数のベクトル加算によりレベルの総和ΣPs (f)が得
られるとともに信号成分方位Bs (f)が保存される。The signal component V s (f) has the same signal component direction B s (f) at all frequencies, that is, B s (f 1 ) = B s
Since (f 2 ) =... = B s (f K ) = B s , the sum of levels ΣP s (f) is obtained by vector addition of all frequencies, and the signal component direction B s (f) is stored. .
【0017】雑音成分Vn (f)は、各周波数で雑音成
分方位Bn (f)がランダムであるため、全周波数のベ
クトル加算によりレベルが小さくなり、ΣPn (f)=
0となる。Since the noise component direction B n (f) is random at each frequency, the level of the noise component V n (f) is reduced by vector addition of all frequencies, and ΔP n (f) =
It becomes 0.
【0018】前記全周波数の和ベクトルVは、次の数4
で表される。The sum vector V of all frequencies is given by
It is represented by
【0019】[0019]
【数4】 この数4におけるVN-S は、次の数5で表される。(Equation 4) V NS in this equation (4) is expressed by the following equation ( 5).
【0020】[0020]
【数5】 また、数4におけるVE-W は、次の数6で表される。(Equation 5) V EW in Equation 4 is expressed by Equation 6 below.
【0021】[0021]
【数6】 したがって、全周波数の和ベクトルVは、信号成分のレ
ベル及び方位だけが保存された状態に近くなり、S/N
が改善される。(Equation 6) Therefore, the sum vector V of all frequencies is close to a state where only the level and the direction of the signal component are stored, and the S / N
Is improved.
【0022】[0022]
【発明の効果】本発明は、S/Nを改善することができ
るから、広帯域性の水中音響信号の検出が容易である。According to the present invention, since the S / N can be improved, it is easy to detect a wideband underwater acoustic signal.
【図1】本発明の1実施例を示すブロック図である。FIG. 1 is a block diagram showing one embodiment of the present invention.
【図2】本発明の1実施例の作用を説明するための図で
ある。FIG. 2 is a diagram for explaining the operation of one embodiment of the present invention.
【図3】従来の信号検出装置を示すブロック図である。FIG. 3 is a block diagram showing a conventional signal detection device.
1 周波数分析部 2 方位算出部 3 ベクトル加算部 4 検出部 Reference Signs List 1 frequency analysis unit 2 direction calculation unit 3 vector addition unit 4 detection unit
───────────────────────────────────────────────────── フロントページの続き (72)発明者 佐々木 康 千葉県印旛郡白井町大山口2−4−16− 203 (72)発明者 竹井 学 東京都港区芝五丁目7番1号 日本電気 株式会社内 (56)参考文献 特開 昭63−144276(JP,A) 特開 平6−201811(JP,A) ──────────────────────────────────────────────────続 き Continuing on the front page (72) Inventor Yasushi Sasaki 2-3-4 Oyamaguchi, Shirai-machi, Inba-gun, Chiba 203 (72) Inventor Manabu Takei 5-7-1 Shiba, Minato-ku, Tokyo NEC Corporation In-company (56) References JP-A-63-144276 (JP, A) JP-A-6-201811 (JP, A)
Claims (2)
の信号方位値とパワースペクトル値を生成し、これらの
信号方位とパワースペクトル値により周波数毎の信号ベ
クトルを合成し、これらの信号ベクトルを全周波数で加
算して全周波数の和信号ベクトルを求め、次にこの和信
号ベクトルのスレッショルドレベルを越える信号を検出
することを特徴とする信号検出方法。1. A signal direction value and a power spectrum value for each frequency are generated from a received signal of a directional sonobuoy, and a signal vector for each frequency is synthesized based on the signal direction and the power spectrum value. A signal detection method comprising adding a sum signal at all frequencies to obtain a sum signal vector of all frequencies, and then detecting a signal exceeding a threshold level of the sum signal vector.
して周波数毎の位相スペクトル値とパワースペクトル値
を一定時間毎に出力する周波数分析部と、この周波数分
析部から位相スペクトル値とパワースペクトル値を受け
てその時間の信号に対応した周波数毎の信号到来方位を
算出する方位算出部と、前記周波数分析部からのパワー
スペクトル値と前記方位算出部からの信号到来方位によ
り形成される信号ベクトルを加算して全周波数の和信号
ベクトルを求めるベクトル加算部と、このベクトル加算
部からの和信号ベクトルのスレッショルドレベルを越え
る信号を検出する検出部とを有することを特徴とする信
号検出装置。2. A frequency analysis unit for frequency-analyzing a received signal of a directional sonobuoy and outputting a phase spectrum value and a power spectrum value for each frequency at regular intervals, and a phase spectrum value and a power spectrum from the frequency analysis unit. A direction calculating unit for receiving a value and calculating a signal arrival direction for each frequency corresponding to the signal at that time, a signal vector formed by the power spectrum value from the frequency analysis unit and the signal arrival direction from the direction calculation unit And a detection unit for detecting a signal exceeding a threshold level of the sum signal vector from the vector addition unit.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP29649494A JP2648110B2 (en) | 1994-11-30 | 1994-11-30 | Signal detection method and device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP29649494A JP2648110B2 (en) | 1994-11-30 | 1994-11-30 | Signal detection method and device |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH08152465A JPH08152465A (en) | 1996-06-11 |
JP2648110B2 true JP2648110B2 (en) | 1997-08-27 |
Family
ID=17834286
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP29649494A Expired - Lifetime JP2648110B2 (en) | 1994-11-30 | 1994-11-30 | Signal detection method and device |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2648110B2 (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102834707A (en) * | 2010-04-12 | 2012-12-19 | 佳能株式会社 | Acoustic wave imaging apparatus and acoustic wave imaging method |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4836497B2 (en) * | 2005-06-15 | 2011-12-14 | 株式会社東芝 | Radar apparatus and signal processing method thereof |
JP5111778B2 (en) * | 2006-04-28 | 2013-01-09 | 株式会社東芝 | Radar apparatus and signal processing method thereof |
JP4900229B2 (en) * | 2007-12-18 | 2012-03-21 | 日本電気株式会社 | Target signal detection method, apparatus and program |
JP5093616B2 (en) * | 2009-01-26 | 2012-12-12 | 日本電気株式会社 | Target signal detection apparatus, method and program |
JP5423713B2 (en) * | 2011-03-30 | 2014-02-19 | 日本電気株式会社 | Active sonar device and signal processing method thereof |
JP6361546B2 (en) * | 2015-03-24 | 2018-07-25 | 三菱電機株式会社 | Signal processing device |
CN109632071B (en) * | 2018-11-08 | 2020-11-20 | 中国船舶重工集团公司七五0试验场 | Method and device for generating underwater acoustic environment noise data based on time-frequency characteristics |
-
1994
- 1994-11-30 JP JP29649494A patent/JP2648110B2/en not_active Expired - Lifetime
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102834707A (en) * | 2010-04-12 | 2012-12-19 | 佳能株式会社 | Acoustic wave imaging apparatus and acoustic wave imaging method |
US8929174B2 (en) | 2010-04-12 | 2015-01-06 | Canon Kabushiki Kaisha | Acoustic wave imaging apparatus and acoustic wave imaging method |
CN102834707B (en) * | 2010-04-12 | 2015-07-29 | 佳能株式会社 | Acoustic imaging equipment |
Also Published As
Publication number | Publication date |
---|---|
JPH08152465A (en) | 1996-06-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8577055B2 (en) | Sound source signal filtering apparatus based on calculated distance between microphone and sound source | |
US20080247566A1 (en) | Sound source localization system and sound source localization method | |
US20150016628A1 (en) | Method and circuitry for direction of arrival estimation using microphone array with a sharp null | |
JP2648110B2 (en) | Signal detection method and device | |
US9820043B2 (en) | Sound source detection apparatus, method for detecting sound source, and program | |
EP3232219B1 (en) | Sound source detection apparatus, method for detecting sound source, and program | |
JP6862799B2 (en) | Signal processing device, directional calculation method and directional calculation program | |
WO2019227353A1 (en) | Method and device for estimating a direction of arrival | |
JP4977849B2 (en) | Radio wave arrival direction detector | |
JP5093616B2 (en) | Target signal detection apparatus, method and program | |
JP4015989B2 (en) | Wave number estimation apparatus and wave number estimation method | |
JP6933303B2 (en) | Wave source direction estimator, wave source direction estimation method, and program | |
JP3139420B2 (en) | Signal detection method | |
JPS6010190A (en) | Mutual signal processing system for passive sonar | |
JP2005017143A (en) | Meteorological radar signal processor | |
JP2717331B2 (en) | Direction estimation method using narrowband signal | |
JP2845220B2 (en) | Sound source direction detection device | |
JP2001051039A (en) | Signal processing system | |
JP7056739B2 (en) | Wave source direction estimator, wave source direction estimation method, and program | |
JPH05196717A (en) | Sonar-signal detecting and processing method | |
JP2779435B2 (en) | Sonar signal detection processor | |
JPS59104577A (en) | Correcting system of estimated error of angle | |
JPH0766050B2 (en) | Direction measurement method | |
JPH0882662A (en) | Direction detection apparatus | |
JPH10288659A (en) | Interference-wave suppressing device |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 19970408 |
|
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 |
|
S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090509 Year of fee payment: 12 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100509 Year of fee payment: 13 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110509 Year of fee payment: 14 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110509 Year of fee payment: 14 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120509 Year of fee payment: 15 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120509 Year of fee payment: 15 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130509 Year of fee payment: 16 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130509 Year of fee payment: 16 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
EXPY | Cancellation because of completion of term |