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]
【産業上の利用分野】本発明は、指向性ソノブイの受波
信号を処理し広帯域性の水中音響信号を検出する信号検
出方法及び装置に関する。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]
【従来の技術】従来より、この種の信号検出装置とし
て、図３に示すように、周波数分析部１１と、これから
の出力を受ける検出部１２とを有している。周波数分析
部１１は、指向性ソノブイの受波信号ｘ（ｔ）を受けて
周波数分析を行い、パワースペクトル値Ｐ（ｆ）を一定
時間毎に出力する。検出部１２は、パワースペクトル値
Ｐ（ｆ）を受けて、この値がスレッショルドレベルを越
える信号を検出する。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]
【発明が解決しようとする課題】しかし、この従来の信
号検出装置は、周波数分析を行うことで、狭帯域分析に
より広帯域にわたって分布する雑音のレベルを低下させ
ているから、広帯域性の信号のレベルも低下するために
Ｓ／Ｎの効果が小さく信号検出が困難であるという問題
がある。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] 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]
【課題を解決するための手段】本発明は、前記課題を解
決するために、指向性ソノブイの受波信号から周波数毎
の信号方位値とパワースペクトル値を生成し、これらの
信号方位とパワースペクトル値により周波数毎の信号ベ
クトルを合成し、これらの信号ベクトルを全周波数で加
算して全周波数の和信号ベクトルを求め、次にこの和信
号ベクトルのスレッショルドレベルを越える信号を検出
することを特徴とする。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.
【０００６】また、本発明は、指向性ソノブイの受波信
号を周波数分析して周波数毎の位相スペクトル値とパワ
ースペクトル値を一定時間毎に出力する周波数分析部
と、この周波数分析部から位相スペクトル値とパワース
ペクトル値を受けてその時間の信号に対応した周波数毎
の信号到来方位を算出する方位算出部と、前記周波数分
析部からのパワースペクトル値と前記方位算出部からの
信号到来方位により形成される信号ベクトルを加算して
全周波数の和信号ベクトルを求めるベクトル加算部と、
このベクトル加算部からの和信号ベクトルのスレッショ
ルドレベルを越える信号を検出する検出部とを有するこ
とを特徴とする。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]
【実施例】次に、本発明の実施例を図面に基いて詳細に
説明する。Next, an embodiment of the present invention will be described in detail with reference to the drawings.
【０００８】図１は、本発明の１実施例を示すブロック
図である。本発明の信号検出装置は、周波数分析部１
と、この周波数分析部１からの出力を受ける方位算出部
２と、前記周波数分析部１および方位算出部２からの出
力を受けるベクトル加算部３と、このベクトル加算部３
の出力を受ける検出部４とを有する。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
【０００９】前記周波数分析部１は、指向性ソノブイの
受波信号ｘ（ｔ）を受けて周波数分析を行い、パワース
ペクトル値Ｐ（ｆ）および位相スペクトル値θ（ｆ）を
一定時間毎に出力する。前記方位算出部２は、周波数分
析部１からのパワースペクトル値Ｐ（ｆ）および位相ス
ペクトル値θ（ｆ）に基いて、その時間の信号に対応し
た周波数毎の信号到来方位Ｂ（ｆ）を算出する。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.
【００１０】前記ベクトル加算部３は、周波数分析部１
からのパワースペクトル値Ｐ（ｆ）と方位算出部２から
の信号到来方位Ｂ（ｆ）とを受けて、これらにより形成
される周波数毎の信号ベクトルＶ（ｆ）を加算し、次の
数１で表される全周波数の和ベクトルＶを求めて出力す
る。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]
【数１】 前記信号ベクトルＶ（ｆ）は、ベクトルの角度がＢ
（ｆ）であり、ベクトルの大きさがパワースペクトル値
Ｐ（ｆ）であるベクトルである。すなわち、この信号ベ
クトルＶ（ｆ）を直交する２つの方向成分ＶN-S （ｆ）
とＶE-W （ｆ）に分解すると、それぞれ次の数２および
数３で表される。(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]
【数２】 (Equation 2)
【００１３】[0013]
【数３】 前記検出部４は、全周波数の和ベクトルＶを受けて、ス
レッショルドレベルを越える信号を検出する。(Equation 3) The detection unit 4 receives the sum vector V of all frequencies and detects a signal exceeding a threshold level.
【００１４】次に、前記ベクトル加算部３の処理を図２
を参照して説明する。Next, the processing of the vector adder 3 will be described with reference to FIG.
This will be described with reference to FIG.
【００１５】前記ベクトル加算部３で受けたパワースペ
クトル値Ｐ（ｆ）と方位算出部２からの信号到来方位Ｂ
（ｆ）とにより形成される周波数毎の信号ベクトルＶ
（ｆ）は、信号成分Ｖs （ｆ）と雑音成分Ｖn （ｆ）の
和になっている。信号成分Ｖs（ｆ）は一定の方位から
到来する成分であり、雑音成分Ｖn （ｆ）は全方位から
到来する等方性の成分であるとする。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.
【００１６】信号成分Ｖs （ｆ）は全周波数で信号成分
方位Ｂs （ｆ）が等しい、すなわちＢs （ｆ1 ）＝Ｂs
（ｆ2 ）＝…＝Ｂs （ｆK ）＝Ｂs となるから、全周波
数のベクトル加算によりレベルの総和ΣＰs （ｆ）が得
られるとともに信号成分方位Ｂs （ｆ）が保存される。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. .
【００１７】雑音成分Ｖn （ｆ）は、各周波数で雑音成
分方位Ｂn （ｆ）がランダムであるため、全周波数のベ
クトル加算によりレベルが小さくなり、ΣＰn （ｆ）＝
０となる。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.
【００１８】前記全周波数の和ベクトルＶは、次の数４
で表される。The sum vector V of all frequencies is given by
It is represented by
【００１９】[0019]
【数４】 この数４におけるＶN-S は、次の数５で表される。(Equation 4) V NS in this equation (4) is expressed by the following equation ( 5).
【００２０】[0020]
【数５】 また、数４におけるＶE-W は、次の数６で表される。(Equation 5) V EW in Equation 4 is expressed by Equation 6 below.
【００２１】[0021]
【数６】 したがって、全周波数の和ベクトルＶは、信号成分のレ
ベル及び方位だけが保存された状態に近くなり、Ｓ／Ｎ
が改善される。(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]
【発明の効果】本発明は、Ｓ／Ｎを改善することができ
るから、広帯域性の水中音響信号の検出が容易である。According to the present invention, since the S / N can be improved, it is easy to detect a wideband underwater acoustic signal.
【図１】本発明の１実施例を示すブロック図である。FIG. 1 is a block diagram showing one embodiment of the present invention.
【図２】本発明の１実施例の作用を説明するための図で
ある。FIG. 2 is a diagram for explaining the operation of one embodiment of the present invention.
【図３】従来の信号検出装置を示すブロック図である。FIG. 3 is a block diagram showing a conventional signal detection device.
１ 周波数分析部 ２ 方位算出部 ３ ベクトル加算部 ４ 検出部 Reference Signs List 1 frequency analysis unit 2 direction calculation unit 3 vector addition unit 4 detection unit
───────────────────────────────────────────────────── フロントページの続き (72)発明者 佐々木 康 千葉県印旛郡白井町大山口２−４−16− 203 (72)発明者 竹井 学 東京都港区芝五丁目７番１号 日本電気 株式会社内 (56)参考文献 特開 昭63−144276（ＪＰ，Ａ) 特開 平６−201811（ＪＰ，Ａ) ──────────────────────────────────────────────────続 き 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 | |
| US8094833B2 (en) | Sound source localization system and sound source localization method | |
| US9479867B2 (en) | Method and circuitry for direction of arrival estimation using microphone array with a sharp null | |
| US9820043B2 (en) | Sound source detection apparatus, method for detecting sound source, and program | |
| JP2648110B2 (en) | Signal detection method and device | |
| JP6862799B2 (en) | Signal processing device, directional calculation method and directional calculation program | |
| EP3232219B1 (en) | Sound source detection apparatus, method for detecting sound source, and 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 | |
| JPS6010190A (en) | Mutual signal processing system for passive sonar | |
| JP2717331B2 (en) | Direction estimation method using narrowband signal | |
| JP2845220B2 (en) | Sound source direction detection device | |
| JP7056739B2 (en) | Wave source direction estimator, wave source direction estimation method, and program | |
| JP3324591B2 (en) | Receiver | |
| 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 | |
| JPS6398577A (en) | Image formation of wave motion | |
| JPH11352204A (en) | Direction finder | |
| JP2007078556A (en) | Signal incoming direction estimation device and method, signal separation device and method, and computer program |
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 |