Energy Separation Algorithms Applied to Sonar Data 

R.K. Lennartsson  , S. McLaughlin , M.J. Levonen , J.W.C. Robinson and L. Persson



Swedish Defence Research Agency, SE-17290 Stockholm, Sweden Institute for Digital Communication, School of Engineering and Electronics University of Edinburgh,Edinburgh, EH9 3JL, UK.

Abstract

hull and the engine. The radiated sound in the water was recorded by bottom-mounted wide-band hydrophones. The sampling rate was 25 kHz for all recording channels. Ambient sea noise was recorded before and after the sea trial. The weather conditions were good throughout the whole experiment with a wind speed below 5 m/s.

Energy separation algorithms (ESA) using Teagers energy operator have been used to analyze modulations in speech resonances, energy separation in bird song and in the resynthesis of musical instrument sounds. The data set used in the analysis here is from a sea trial in the Baltic sea. A 23ft fiberglass boat was equipped with accelerometers for source signal recordings from the engine and the hull. The sound propagated through the water was picked up with a passive hydrophone array. The analysis of the accelerometer data was used to identify frequencies associated with engine and propeller rotation. Frequency components not harmonically related to either the engine crank shaft or drive were identified in the hydrophone data. The ESA was then applied to the none related frequencies. It is shown that ESA can be used as a tool to identify stability and time dependency of frequencies.

Previously [1] bispectral analysis was used to assess the level of quadratic coupling present in the dataset discussed above. However, it was difficult to discern precisely the level of coupling because much of the coupling lay between harmonics of the generating signals which were themselves phase coupled by the stiff mechanism in the engine/drive gear assembly. This paper reports an analysis of the frequency content of the data, in particular frequency components not harmonically related to either the rotation of the engine crank shaft or drive (gear box and propeller) were identified. In the work presented here we use the ESA to assess the level of stationarity and stability of the unrelated frequency components in the original data set. The aim is to use this assessment at a later date to select particular segments for an analysis by bicoherence.

I. I NTRODUCTION In sonar signal processing identification and classification of different underwater sources are crucial. In passive sonar analysis the estimation problems are often focused on spectral separation and characterization of the narrow-band signals emitted from different parts of the target, e.g. engine, drive and hull. In many times the power spectrum is able to separate these sources. If a more precise description of targets is required the relation and coupling between frequencies may need to be uncovered. One possible way to do this is by using higher-order spectra, e.g. [1]. From the bispectrum the phase coupling patterns can be obtained and thereby be used as a tool for association of different physical mechanisms such as engine and drive. In all spectral analysis it is important to use as significant spectral values as possible for source parameter estimation. In spectral analysis second-order stationarity is normally assumed, at least within the observation window. The estimated spectra may be erroneous if the data contains disturbances and transients. In this paper use is made of Teagers Energy operator as proposed by Maragos and his co-workers [2, 3, 4] to determine the stability of the power spectral peaks. The data is from a sea trial conducted in a shallow bay in the Baltic Sea off the east coast of Sweden. The water depth at the site was approximately 30 m. The trials were performed using a 23-ft. fiberglass motor boat with a 4-cyl. 4-stroke, turbo-charged diesel engine as sound source. The boat was equipped with accelerometers for source signal recordings. We allowed the boat run on straight tracks from far distance with a constant speed and approaching and passing an array of bottom mounted hydrophones. Engine and hull vibrations were recorded with accelerometers fitted directly on the

II. E NERGY S EPARATION The use of the Teager Energy operator to separate the amplitude and frequency modulated (AM/FM) components of mixed signals with AM and FM components. Much of the early work to elucidate the properties of this operator was performed by Maragos and his co-workers [2, 3]. For example, for discrete-time signals the energy operator is given by,



 

       

(1)

which can be used to find the amplitude and frequency of signals of the form,

  !#"%$'&)(* ,+.-/ 021'3

(2)

The energy operator is applied to both the signal and the first difference, given by

4  56

7!

 83

(3)

The instantaneous amplitude and frequency can be obtained from, A

 and

1

+ - :9<;=$>$>&(.?@. B A4   C=D E E " GF A   ,  +.C-J  3 H*I

(4)

(5)

145

s1: lopa_la004 Frequency [Hz]

300

200

130 0

100

0 0

20

40

60 Time [sec]

III. R ESULTS

80

60 Time [s]

80

100

20

40

60 Time [s]

80

100

0.2

100

Figure 2: Results of ESA analysis at frequency component at 136.7 Hz in hydrophone data (in Block 1).

AND

includes hydrophone data and data from an engine mounted accelerometer.

A NALYSIS

An extensive analysis of the data set was carried out considering the power spectra, short-time fourier transform (STFT) and finally ESA was applied to the frequency components which appeared to be not harmonically related to any rotating elements. In Fig. 1 the STFT is displayed. In this case the boat is moving with the engine running at approximately 3000rpm (actually 2712rpm) with the drive connected. The frequencies are varying a lot during the first 20 seconds. The reason is that the recording is started when the boat was standing still and that it takes some time for the boat to reach a stable speed with a fixed rpm. For the rest of the run the frequencies are quite stable except for the 90 to 95 [sec] section (the rpm was shortly lowered by the driver).

A. Block 1 In this sub-section the data in block 1 is analyzed with power spectrum and ESA. In block 1 the boat is far away from the hydrophone and slowly approaching it. The spectra were generated using Matlab with a 65536 point FFT resulting in a frequency resolution of approximately 0.4 Hz. The power spectrum of both the hydrophone (see Fig. 3 and 4) and accelerometer (see Fig. 5 and 6) data contains strong peaks for frequencies up to 600 Hz. All strong peaks have been identified and these are indicated in the Figures. s1: lopa_004 (25 to 40[sec])

3

In Fig.2 the output of ESA applied to the strongest component identified in the STFT of the hydrophone data (136.7 Hz later found to be the rotation frequency of the propeller shaft modulated by the number of propeller blades) is displayed. The ESA shows the same over all picture as the STFT but in greater detail, especially for the amplitude. The frequency is very stable (except for the time segments mentioned above) but the amplitude shows some variance.

9

2

10 7

1

23 1 4

Power

0

5

6

8

11 13 12

16

17

18

14 15

−1 −2 −3 −4

For the rest of the analysis the data was considered as three distinct data sets divided as follows:

−5 −6 0

Block 1: The engine running at 3000rpm with the boat moving (25-40 seconds)

K

40

0.4

0 0

Figure 1: Spectrum of frequency range 0 to 300 Hz for engine accelerometer with the boat moving

K

20

150

50

K

135

Amplitude

Frequency [Hz]

250

140

50

100

150 200 Frequency [Hz]

250

300

Figure 3: Spectrum of frequency range 0 to 300 Hz for the hydrophone data in Block 1.

Block 2: The engine running at 3000rpm with the boat moving (65-80 seconds)

In addition, all identified peaks were tabulated in three tables (not shown here) which indicate the following:

Block 3: The engine running at 3000rpm with the boat moving (100-115 seconds)

K

The distinction between the three blocks relates to the position of the boat with relation to the hydrophone array. In block 1 the boat is far away and approaching, block 2 contains the closest point of approach and in block 3 the boat is moving away from the array. The data analyzed here

K

The frequency peaks that are related to the rotation of the engine axes

K 2

The frequency peaks related to the propeller shaft The frequency peaks that are not related to anyone of the previously mentioned frequencies in the other

s1: lopa_004 (25 to 40[sec])

3 2

2 2 1 34

1

5

6

7

9

8

10

11

13 14 12

1 15 16

−1 −2

−4

−4

−5

−5 400

450 500 Frequency [Hz]

550

Figure 4: Spectrum of frequency range 300 to 600 Hz for the hydrophone data in Block 1.

7

0 −1

1

−2

3

8 4 5

2

−3

9

14

12

350

400

450 500 Frequency [Hz]

550

600

120

15 16 Frequency [Hz]

13 11 12

11

6

and the sound source remains at fairly constant frequencies so the influence from earlier and later arrivals (surface and sea-floor reflections) interferes with main propagation path giving the varying amplitude that can be seen in the figures mentioned above.

2 1

10

Figure 6: Spectrum of frequency range 300 to 600 Hz for the engine accelerometer data in Block 1.

ch1: lopa_004 (25 to 40[sec])

3

4 5

−6 300

600

7

8

−2 −3

350

2

3

−1

−3

−6 300

9 1

0 Power

Power

0

Power

ch1: lopa_004 (25 to 40[sec])

3

10

6

115 110 105 25

−4

30

Time [s]

35

40

35

40

0.2

−6 0

50

100

150 200 Frequency [Hz]

250

Amplitude

−5 300

Figure 5: Spectrum of frequency range 0 to 300 Hz for the engine accelerometer data in Block 1.

0.15 0.1 0.05 0 25

tables

30

Time [s]

Figure 7: Results of ESA analysis at frequency component at 109.9 Hz in hydrophone data (in Block 1).

From the analysis of these spectra and others not illustrated here it is apparent that many of the peaks are multiples of 22.6 Hz, which corresponds to one half of the rotation frequency of the engine axes. In addition there are multiples of 19.5 Hz which corresponds to the frequency of rotation of the propeller shaft. The gear ratio between the propeller shaft and the engine axes is 2.31. If these frequencies are discounted the frequencies which are not related to either of these are only found in the hydrophone data. Discounting very low frequencies where the SNR was low, the following frequency peaks were found 109.9 Hz, 333 Hz, 392.1 Hz, 431 Hz, 489.6 Hz and 508.9 Hz.

Frequency [Hz]

340 335 330 325 25

30

Time [s]

35

40

35

40

Amplitude

0.2

The ESA was then applied to each of the unrelated frequencies. In Fig. 7 to 12 the output of the ESA is displayed. In the upper panel the ESA estimated frequency is displayed and in the lower panel the amplitude is shown. The most stable components are 333 Hz, 489.6 Hz, 508.9 Hz, while the 109.9 Hz peak shows the largest variance. The variations in amplitude can be explained by multipath propagation of the sound waves. The site is very shallow

0.15 0.1 0.05 0 25

30

Time [s]

Figure 8: Results of ESA analysis at frequency component at 333 Hz in hydrophone data (in Block 1).

3

395 390 385 25

30

Time [s]

35

505

0.2

0.2

0.15

0.15

0.1 0.05 0 25

30

Time [s]

35

30

Time [s]

35

40

80

30

Time [s]

35

Frequency [Hz]

Frequency [Hz]

430

40

0.15

70 60 50 65

0.1 0.05

70

Time [s]

75

80

75

80

30

Time [s]

35

Amplitude

0.2 40

Figure 10: Results of ESA analysis at frequency component at 431 Hz in hydrophone data (in Block 1).

0.15 0.1 0.05 0 65

500 Frequency [Hz]

40

identifying the remaining frequencies. Interestingly only one frequency was identified in the hydrophone data at 64.5 Hz. In Fig. 13 the output of the ESA is displayed. The amplitude is quite stable while the frequency shows some variance.

0.2

70

Time [s]

Figure 13: Results of ESA analysis at frequency component at 64.5 Hz in hydrophone data (in Block 2)

495 490

C. Block 3

485 480 25

30

Time [s]

35

In this sub-section the data in Block 3 is analyzed, here the boat is moving away from the hydrophone.

40

Following the same procedure as for Blocks 1 and 2 and discounting those frequencies associated with the mechanical vibrations results in identifying the remaining frequencies. In this case two frequencies were identified in the hydrophone data at 469.6 Hz and 538.6 Hz. In Fig. 14 to 15 the output of the ESA is displayed. The 469.6 Hz is quite stable in both amplitude and frequency, while the 538.6 Hz peak shows quite large variance in frequency.

0.2 Amplitude

35

Figure 12: Results of ESA analysis at frequency component at 508.9 Hz in hydrophone data (in Block 1).

435

0 25

Time [s]

0.05 0 25

440

425 25

30

0.1

40

Figure 9: Results of ESA analysis at frequency component at 392.1 Hz in hydrophone data (in Block 1).

Amplitude

510

500 25

40

Amplitude

Amplitude

515 Frequency [Hz]

Frequency [Hz]

400

0.15 0.1 0.05 0 25

30

Time [s]

35

40

Figure 11: Results of ESA analysis at frequency component at 489.6 Hz in hydrophone data (in Block 1).

IV. C ONCLUSIONS

B. Block 2

This paper has presented some preliminary analysis of a sonar data set which has been analyzed previously using normalized bispectrum. The ESA gives a more detailed picture of the frequency and amplitude variations than the spectrum or the STFT (for one frequency component at a time). The ESA results shown indicate some consistency

In this sub-section data from Block 2 is analyzed. The boat moving towards and passing over the hydrophone. Following the same procedure as for Block 1 (the spectra are not shown here) and discounting those frequencies associated with the mechanical vibrations results in 4

analysis,” IEEE Trans. Signal Processing, pp. 3024– 3051, October 1993. [4] P. Maragos, J. Kaiser and T. Quateri, “On amplitude and frequency demodulation using energy operators,” IEEE Trans. Signal Processing, pp. 1532–1550, April 1993.

Frequency [Hz]

480 475 470 465 460 100

105

Time [s]

110

115

110

115

Amplitude

0.2 0.15 0.1 0.05 0 100

105

Time [s]

Figure 14: Results of ESA analysis at frequency component at 469.6 Hz in hydrophone data (in Block 3).

Frequency [Hz]

550 545 540 535 530 100

105

Time [s]

110

115

110

115

Amplitude

0.2 0.15 0.1 0.05 0 100

105

Time [s]

Figure 15: Results of ESA analysis at frequency component at 538.6 Hz in hydrophone data (in Block 3).

in the unrelated frequency components and help identify areas where the components are likely to be strongest. The source of the unrelated frequencies are not known to us, but most likely they are generated by some rotating mechanism (other than engine crank shaft or drive) on the boat. The analysis presented here has indicated that in order to perform the data analysis using normalized bispectrum it is necessary to identify those segments where stationarity is strongest and few transients exist. The ESA has enabled this and can hence be a useful tool in the analysis of sonar data.

V. R EFERENCES [1] R.K. Lennartsson, L. Persson, J.W.C. Robinson and S. McLaughlin, “Passive signature characterization by bispectral analysis,” Proc. IEEE Nordic Signal Processing Symposium, NORSIG 2000, Kolmarden, Sweden, June 2000, pp. 445–448. [2] B. Santhanam and P. Maragos, “Energy demodulation of two component am-fm signal mixtures,” IEEE Signal Processing Letters, vol. 3, pp. 294–298, November 1996. [3] P. Maragos, J. Kaiser and T. Quateri, “Energy separation in signal modulations with application to speech 5

Energy Separation Algorithms Applied to Sonar Data - Semantic Scholar

The analysis of the accelerometer data was used to identify frequencies associated with engine and propeller rotation. Frequency components not harmonically ...

180KB Sizes 12 Downloads 287 Views

Recommend Documents

Energy Separation Algorithms Applied to Sonar Data - Semantic Scholar
The analysis of the accelerometer data was used to identify frequencies associated with engine and propeller rotation. Frequency components not harmonically ...

Weighted Automata Algorithms - Semantic Scholar
The mirror image of a string x = x1 ···xn is the string xR = xnxn−1 ··· x1. Finite-state transducers are finite automata in which each transition is augmented with an ...

Broadcast Gossip Algorithms - Semantic Scholar
Email:{tca27,mey7,as337}@cornell.edu. Abstract—Motivated by applications to wireless sensor, peer- to-peer, and ad hoc networks, we study distributed ...

Weighted Automata Algorithms - Semantic Scholar
A finite-state architecture for tokenization and grapheme-to- phoneme conversion in multilingual text analysis. In Proceedings of the ACL. SIGDAT Workshop, Dublin, Ireland. ACL, 1995. 57. Stephen Warshall. A theorem on Boolean matrices. Journal of th

Renewable Energy and Energy Efficiency ... - Semantic Scholar
Oct 18, 2013 - Energy Efficient Appliance Tax Credit for Manufacturers . ..... FY2011; $0 for FY2012; $0 for FY2013; data for FY2014 is currently unavailable.

Renewable Energy and Energy Efficiency ... - Semantic Scholar
Oct 18, 2013 - Bioenergy Technologies Office – Financial Opportunities online ... contract for cost-shared research, development, and demonstration projects. ... production, delivery, and storage technologies; overcome technical barriers to.

On Approximation Algorithms for Data Mining ... - Semantic Scholar
Jun 3, 2004 - Since the amount of data far exceeds the amount of workspace available to the algorithm, it is not possible for the algorithm to “remember” large.

On Approximation Algorithms for Data Mining ... - Semantic Scholar
Jun 3, 2004 - The data stream model appears to be related to other work e.g., on competitive analysis [69], or I/O efficient algorithms [98]. However, it is more ...

On Approximation Algorithms for Data Mining ... - Semantic Scholar
Jun 3, 2004 - problems where distance computations and comparisons are needed. In high ..... Discover the geographic distribution of cell phone traffic at.

Fast exact string matching algorithms - Semantic Scholar
LITIS, Faculté des Sciences et des Techniques, Université de Rouen, 76821 Mont-Saint-Aignan Cedex, France ... Available online 26 January 2007 ... the Karp–Rabin algorithm consists in computing h(x). ..... programs have been compiled with gcc wit

The WebTP Architecture and Algorithms - Semantic Scholar
satisfaction. In this paper, we present the transport support required by such a feature. ... Multiple network applications run simultaneously on a host computer, and each applica- tion may open ...... 4, pages 365–386, August. 1995. [12] Jim ...

MATRIX DECOMPOSITION ALGORITHMS A ... - Semantic Scholar
solving some of the most astounding problems in Mathematics leading to .... Householder reflections to further reduce the matrix to bi-diagonal form and this can.

Adaptive Algorithms Versus Higher Order ... - Semantic Scholar
sponse of these channels blindly except that the input exci- tation is non-Gaussian, with the low calculation cost, com- pared with the adaptive algorithms exploiting the informa- tion of input and output for the impulse response channel estimation.

all pairs shortest paths algorithms - Semantic Scholar
Given a communication network or a road network one of the most natural ... ranging from routing in communication networks to robot motion planning, .... [3] Ming-Yang Kao, Encyclopedia of Algorithms, SpringerLink (Online service).

all pairs shortest paths algorithms - Semantic Scholar
In this paper we deal with one of the most fundamental problems of Graph Theory, the All Pairs Shortest. Path (APSP) problem. We study three algorithms namely - The Floyd- Warshall algorithm, APSP via Matrix Multiplication and the. Johnson's algorith

Minimax Optimal Algorithms for Unconstrained ... - Semantic Scholar
Jacob Abernethy∗. Computer Science and Engineering .... template for (and strongly motivated by) several online learning settings, and the results we develop ...... Online convex programming and generalized infinitesimal gradient ascent. In.

Non-Negative Matrix Factorization Algorithms ... - Semantic Scholar
Keywords—matrix factorization, blind source separation, multiplicative update rule, signal dependent noise, EMG, ... parameters defining the distribution, e.g., one related to. E(Dij), to be W C, and let the rest of the parameters in the .... contr

MATRIX DECOMPOSITION ALGORITHMS A ... - Semantic Scholar
... of A is a unique one if we want that the diagonal elements of R are positive. ... and then use Householder reflections to further reduce the matrix to bi-diagonal form and this can ... http://mathworld.wolfram.com/MatrixDecomposition.html ...

The WebTP Architecture and Algorithms - Semantic Scholar
bandwidth-guaranteed service, delay-guaranteed service and best-effort service ..... as one of the benefits of this partition, network functions can be integrated ...

Fast data extrapolating - Semantic Scholar
near the given implicit surface, where image data extrapolating is needed. ... If the data are extrapolated to the whole space, the algorithm complexity is O(N 3. √.

Reactive Data Visualizations - Semantic Scholar
of the commercial visualization package Tableau [4]. Interactions within data visualization environments have been well studied. Becker et al. investigated brushing in scatter plots [5]. Shneiderman et al. explored dynamic queries in general and how

CASA Based Speech Separation for Robust ... - Semantic Scholar
and the dominant pitch is used as a main cue to find the target speech. Next, the time frequency (TF) units are merged into many segments. These segments are ...

Multiple-Layer Neural Network Applied to Phase ... - Semantic Scholar
collimated fringes (usually with sinusoidal intensity pattern) will be projected onto the target. 20 surface. Cameras would be placed from a different view angle. Deformation of fringe pattern. 21 will appear in the captured images and surface shape