Nothing Special   »   [go: up one dir, main page]

Next Article in Journal
Light Energy Conversion Surface with Gold Dendritic Nanoforests/Si Chip for Plasmonic Polymerase Chain Reaction
Previous Article in Journal
Study of Postural Stability Features by Using Kinect Depth Sensors to Assess Body Joint Coordination Patterns
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Underwater Location Beacon Signal Detection Method Based on Mixing and Normalizing Stochastic Resonance

1
Acoustic Science and Technology Laboratory, Harbin Engineering University, Harbin 150001, China
2
Key Laboratory of Marine Information Acquisition and Security (Harbin Engineering University), Ministry of Industry and Information Technology, Harbin 150001, China
3
College of Underwater Acoustic Engineering, Harbin Engineering University, Harbin 150001, China
4
Qingdao Haina Underwater Information Technology Co., Ltd., Qingdao 266500, China
5
715th. Institute of China Shipbuilding Industry Company, Hangzhou, Zhejiang 310023, China
*
Author to whom correspondence should be addressed.
Sensors 2020, 20(5), 1292; https://doi.org/10.3390/s20051292
Submission received: 16 January 2020 / Revised: 16 February 2020 / Accepted: 25 February 2020 / Published: 27 February 2020
(This article belongs to the Section Remote Sensors)
Figure 1
<p>The variation of the received ULB signal <span class="html-italic">SNR</span> at different distances.</p> ">
Figure 2
<p>The received signal at 2 km distance. (<b>a</b>) waveform of the received signal; (<b>b</b>) frequency spectrum of the received signal.</p> ">
Figure 3
<p>The output SNR versus noise intensity <span class="html-italic">D</span>.</p> ">
Figure 4
<p>A simulation of classical SR. The relative parameters corresponding to Equation (7) configured as: <span class="html-italic">A</span> = 0.1, <span class="html-italic">f</span> = 0.01, <span class="html-italic">a</span> = 1, <span class="html-italic">b</span> = 1, <span class="html-italic">D</span> = 0.31, <span class="html-italic">f<sub>s</sub></span> = 5. (<b>a</b>,<b>b</b>) are the waveform and frequency spectrum of the input; (<b>c</b>,<b>d</b>) are waveform and frequency spectrum of the SR output.</p> ">
Figure 5
<p>The SR of high-frequency signal. The relative parameters: <span class="html-italic">A</span> = 0.1, <span class="html-italic">f</span> = 10, <span class="html-italic">a</span> = 1, <span class="html-italic">b</span> = 1, <span class="html-italic">D</span> = 0.31, <a href="#sensors-20-01292-f005" class="html-fig">Figure 5</a>. (<b>a</b>,<b>b</b>) are the waveform and frequency spectrum of the input; (<b>c</b>,<b>d</b>) are waveform and frequency spectrum of the SR output.</p> ">
Figure 6
<p>The SR of signal with large noise. The relative parameters: <span class="html-italic">A</span> = 0.5, <span class="html-italic">f</span> = 0.01, <span class="html-italic">a</span> = 1, <span class="html-italic">b</span> = 1, <span class="html-italic">D</span> = 7.75, <span class="html-italic">f<sub>s</sub></span> = 5. (<b>a</b>,<b>b</b>) are the waveform and frequency spectrum of the input; (<b>c</b>,<b>d</b>) are waveform and frequency spectrum of the SR output.</p> ">
Figure 7
<p>The frequency spectrum of the SR output when the ratio <span class="html-italic">k</span> is small. The relative parameters: <span class="html-italic">A</span> = 0.1, <span class="html-italic">f</span> = 0.2, <span class="html-italic">a</span> = 1, <span class="html-italic">b</span> = 1, <span class="html-italic">D</span> = 0.31, <span class="html-italic">f<sub>s</sub></span> = 5.</p> ">
Figure 8
<p>The flow of MNSR.</p> ">
Figure 9
<p>An example shown the effectivity of MNSR. The relative parameters: <span class="html-italic">A</span> = 0.5, <span class="html-italic">f</span> = 1000, <span class="html-italic">a</span> = 1000, <span class="html-italic">b</span> = 1.45 × 10<sup>8</sup>, <span class="html-italic">f<sub>c</sub></span> = 990 Hz, <span class="html-italic">D</span> = 7.75, <span class="html-italic">f<sub>s</sub></span> = 5000. (<b>a</b>,<b>b</b>) are the waveform and frequency spectrum of the input. (<b>c</b>,<b>d</b>) are the waveform and frequency spectrum of the low-pass filter output, <span class="html-italic">m</span>(<span class="html-italic">t</span>). (<b>e</b>,<b>f</b>) are the waveform and frequency spectrum of the MNSR output.</p> ">
Figure 9 Cont.
<p>An example shown the effectivity of MNSR. The relative parameters: <span class="html-italic">A</span> = 0.5, <span class="html-italic">f</span> = 1000, <span class="html-italic">a</span> = 1000, <span class="html-italic">b</span> = 1.45 × 10<sup>8</sup>, <span class="html-italic">f<sub>c</sub></span> = 990 Hz, <span class="html-italic">D</span> = 7.75, <span class="html-italic">f<sub>s</sub></span> = 5000. (<b>a</b>,<b>b</b>) are the waveform and frequency spectrum of the input. (<b>c</b>,<b>d</b>) are the waveform and frequency spectrum of the low-pass filter output, <span class="html-italic">m</span>(<span class="html-italic">t</span>). (<b>e</b>,<b>f</b>) are the waveform and frequency spectrum of the MNSR output.</p> ">
Figure 10
<p>The numerical simulation results of MNSR. The relative parameters: <span class="html-italic">A</span> = 0.15, <span class="html-italic">f</span> = 37.5 kHz, <span class="html-italic">a</span> = 3 × 10<sup>4</sup>, <span class="html-italic">b</span> = 4.65 × 10<sup>13</sup>, <span class="html-italic">f<sub>c</sub></span> = 37.8 kHz, <span class="html-italic">D</span> = 0.9, <span class="html-italic">f<sub>s</sub></span> = 150 kHz, <span class="html-italic">N</span> = 1500. (<b>a</b>,<b>b</b>) are the waveform and frequency spectrum of the input; (<b>c</b>,<b>d</b>) are waveform and frequency spectrum of the MNSR output.</p> ">
Figure 11
<p>The result of RFSR applied to the ULB signal detection when the sampling frequency is <span class="html-italic">f<sub>s</sub></span> = 150 kHz, and the sampling data length is <span class="html-italic">N</span> = 1500. The re-scaling ratio of RFSR is <span class="html-italic">R</span> = 3.75 × 10<sup>6</sup>. In addition, the other relative parameters: <span class="html-italic">A</span> = 0.15, <span class="html-italic">f</span> = 37.5 kHz, <span class="html-italic">a</span> = 1, <span class="html-italic">b</span> = 0.344, <span class="html-italic">D</span> = 0.9. (<b>a</b>,<b>b</b>) are the waveform and frequency spectrum of the output.</p> ">
Figure 12
<p>The result of RFSR applied to the ULB signal detection when the sampling frequency is <span class="html-italic">f<sub>s</sub></span> = 18.75 MHz. The sampling data length is still <span class="html-italic">N</span> = 1500 and the other parameters are the same as <a href="#sensors-20-01292-f011" class="html-fig">Figure 11</a>.</p> ">
Figure 13
<p>The result of MSR applied to the ULB signal detection when the sampling frequency is <span class="html-italic">f<sub>s</sub></span> = 150 kHz, and the sampling data length is <span class="html-italic">N</span> = 1500. The carrier frequency is <span class="html-italic">f<sub>c</sub></span> = 37.50001 kHz. In addition, the other relative parameters: <span class="html-italic">A</span> = 0.15, <span class="html-italic">f</span> = 37.5 kHz, <span class="html-italic">a</span> = 1, <span class="html-italic">b</span> = 0.344, <span class="html-italic">D</span> = 0.9. (<b>a</b>,<b>b</b>) are the waveform and frequency spectrum of the output.</p> ">
Figure 14
<p>The schematic diagram of experiment.</p> ">
Figure 15
<p>The input and output of a tank experiment where MNSR used to detect the ULB signal. The relative parameters: <span class="html-italic">f</span> = 37.5 kHz, <span class="html-italic">a</span> = 3 × 10<sup>4</sup>, <span class="html-italic">b</span> = 1.9 × 10<sup>14</sup>, <span class="html-italic">f<sub>c</sub></span> = 37.8 kHz, <span class="html-italic">D</span> = 0.22, <span class="html-italic">f<sub>s</sub></span> = 150 kHz. (<b>a</b>,<b>b</b>) are the waveform and frequency spectrum of the input; (<b>c</b>,<b>d</b>) are waveform and frequency spectrum of the MNSR output.</p> ">
Versions Notes

Abstract

:
A flight data recorder (FDR) is an electronic recording device placed in an aircraft for the purpose of facilitating the investigation of aviation accidents. If an aircraft crashes into water, an underwater locator beacon (ULB), which is installed on the FDR, is triggered by water immersion, and emits an ultrasonic 10 ms pulse signal once per second at 37.5 kHz. This pulse signal can be detected by sonar equipment. However, the ULB signal only can be detectable 1–2 kilometers from the surface in normal conditions. Stochastic resonance (SR) is a rising theory in the field of weak signal detection. The classical stochastic resonance limits state that the input must be small-parameter and the sampling frequency must be 50 times higher than the signal frequency. It cannot be applied to the ULB signal detection. To resolve this problem, this paper presents a novel approach named mixing and normalizing stochastic resonance (MNSR). By mixing the ULB signal and normalizing SR system parameters, MNSR provides a new way to detect weak ULB signal. Meanwhile, we propose the parameters adjustment method of MNSR. We prove the effectiveness through numerical simulation. An experiment in a tank is employed to verify the practicability of this method.

1. Introduction

When an aircraft crashes into water, an underwater locator beacon (ULB) [1,2], which is installed on the FDR (commonly known as the black box) for fast searching in the ocean, is triggered by water immersion, and emits an ultrasonic pulse signal. Due to the ULB volume, the size of the transducer and the battery capacity are both limited. The ULB signal frequency is high and the source level (SL) is low. However, the transmission loss level (TL) in water is large. These factors result in the ULB signal being extremely weak at a long distance. In addition, the signal-to-noise ratio (SNR) is very low. In contrast to conventional methods [3,4,5], stochastic resonance (SR) is a popular method in weak signal detection. SR can transfer the energy from noise to the signal, which will amplify the weak signal and improve the output SNR [6,7,8,9].
The classical stochastic resonance theory in bistable systems is studied under the adiabatic approximation conditions. In addition, there is a small-parameter restriction [10,11,12], which means the signal amplitude and frequency and noise intensity should be much less than 1. Therefore, the classical stochastic resonance cannot detect the large-parameter (large-amplitude or high-frequency) signal directly. To make SR more practical, researchers have proposed several large-parameter stochastic resonance (LPSR) methods, such as twice sampling stochastic resonance [11], re-scaling frequency stochastic resonance (RFSR) [13], system parameters normalization transform stochastic resonance, or parameters normalized stochastic resonance (PNSR) [14]. These methods transform the large-parameter signal into a small-parameter signal which satisfies the small-parameter conditions numerically. The disadvantage is the requirement of a high sampling frequency (at least 50 times to the signal frequency, and generally 200 times or more). It is hard to achieve in engineering application scenarios when the weak signal frequency is high. Another LPSR method, modulated stochastic resonance (MSR) [15], modulates the high-frequency weak signal immersed with noise to a low difference frequency signal to implement LPSR. MSR requires a long signal duration to ensure that the duration contains at least two complete periods of the low difference frequency signal. To realize the detection of ULB signal which is high frequency and short pulse by SR method, mixing and normalizing stochastic resonance (MNSR) is proposed in this paper. Through mixing the ULB signal with a carrier signal, the ULB signal detection can be replaced by the “low” frequency difference signal detection. Here, the “low” means the difference signal frequency is lower than the ULB signal frequency and the sampling frequency. Under the restriction that the sampling frequency must be 50 times higher than the signal frequency, MNSR allows the sampling frequency to be much lower than RFSR, which detects the high-frequency signal directly. In addition, MNSR can detect the large-parameter difference signal. By normalizing the barrier parameters, MNSR allows the difference signal frequency to be greater than 1, so that MNSR allows a signal duration shorter than MSR while MSR needs the difference signal frequency to be less than 1. Compared with frequency-shifted and re-scaling stochastic resonance (FRSR) [16] or frequency exchange and re-scaling stochastic resonance (FERSR) [17], they are more focused on how to choose the re-scaling ratio and mention the selection of non-linear system parameters less. We give a parameter adjustment method is in this paper.
The paper is organized as follows. Section 2 introduces the ULB signal and its underwater propagation. Section 3 presents the classical SR theory and discuss its existing restrictions. In Section 4, a method named mixing and normalizing stochastic resonance (MNSR) is proposed and its parameter adjustment method is given as well. In Section 5, both a numerical simulation and an experiment in results verify the effectiveness of the MNSR and the parameter adjustment method. Conclusions are provided in Section 6.

2. Characteristics of the ULB Signal and Its Underwater Propagation

The ULB is a cylinder with a size of φ33 × 100 mm, which will transmit an acoustic CW pulse underwater and its frequency is 37.5 kHz, pulse width is 10 ms, repetition interval is 1 s, and SL is 160.5 dB. According to the passive sonar equation [18], the SNR level (SNR) of the signal received at a long distance is the difference between SL, TL and the environmental noise level (NL), i.e., SNR = SLTLNL. TL and NL. It can be calculated through the empirical formula below:
T L = 20 lg r + α r × 10 3
N L = 101 g   f 1.7 + 6 S + 55
where r (m) is the signal propagation distance; α (dB/km) is the oceanic sound absorption coefficient, which is related to the signal frequency f (kHz); and S is the level of sea state. Equation (3) is used to calculate α.
α = 0.109 f 2 1 + f 2 + 40.7 f 2 4100 + f 2 + 3.01 × 10 - 4 f 2
We can get the curve of the received signal SNR from 100 m to 5000 m in Figure 1 using Equations (1)–(3), where the level of the sea state S is set to 3 and the bandwidth of noise is set to 5 kHz.
Under the assumption that the sea noise has Gaussian distribution, Figure 2 is a received signal numerical simulation at 2 km. The sampling frequency fs is 150 kHz. It shows that the received signal totally immersed in noise in the time waveform and the spectral spike can hardly be found at 37.5 kHz.

3. Classical SR

The overdamped motion of a Brownian particle in a bistable system driven by noise and external periodic force can be described by Langevin equation (LE) as
d x d t = U ( x ) + A cos ( 2 π f t ) + n ( t )
where x(t) is the particle trajectory; Acos(2πft) is a periodic signal with amplitude A and frequency f. n(t) is a zero-mean, Gaussian white noise with noise intensity D, i.e.,
n ( t ) = 0 n ( t ) n ( t + τ ) = 2 D δ ( τ )
is the statistical mean value calculation. U(x) is the classical bistable symmetric potential considered to be written as below:
U ( x ) = a 2 x 2 + b 4 x 4
in which a and b are barrier parameters with positive real values. Then Equation (4) can be written as:
d x d t = a x b x 3 + A cos ( 2 π f t ) + n ( t )
For small parameters under the adiabatic approximation conditions, the response of the system to the periodic input signal can be written as
x ( t ) = x ¯ cos ( 2 π f t + ϕ ¯ )
with amplitude x ¯ and a phase lag ϕ ¯ , and the approximate expressions of the amplitude x ¯ is
x ¯ = A x 2 0 D r k r k 2 + π 2 f
Here, r k is Kramers rate and can be written as
r k = a 2 π exp ( a 2 4 b D )
and x 2 0 is the variance of the stationary unperturbed system (A = 0). In addition, then the output SNR of the classical SR can be deduced approximately as follows:
S N R = 2 a 2 A 2 4 b D 2 exp ( a 2 4 b D )
According to Equation (11), the SNR curves when a = 1, b = 1, A = 0.3 is shown in Figure 3. The SNR first increases to maximum and then decreases as the noise increases continually. This phenomenon is called stochastic resonance (SR).
A simulation is performed in Figure 4 with the parameters of Equation (7) is configured as: A = 0.1, f = 0.01, a = 1, b = 1, D = 0.31, respectively. The sampling frequency is fs = 5. The data length is 2000 points; and all the SR output points are used to calculate the frequency spectrum. The output is numerically solved with a fourth-order Runge–Kutta method (this way is used with all the following numerical calculations in this paper, its calculating time step is h = 1/fs). Comparing the spectral spikes at f = 0.01 in Figure 4b,d, the merit of SR in the weak signal detection can be observed: there is an obvious sharp spectral spike at f = 0.01 in Figure 4d, but the spike can hardly be found in Figure 4b. SR can enhance the weak periodic signal immersed in noise, and the signal has been successfully detected by SR.
A limitation of the classical SR theory is the small-parameter restriction: the signal amplitude A and frequency f and the noise intensity D should be far less than 1. An example is shown in Figure 5; the frequency f of the periodic signal is changed to 10 Hz. This is a large parameter greater than 1. In addition, the sampling frequency fs = 5000, but the other parameters are the same as Figure 4. The spike at f = 10 in Figure 5d nearly disappears. This appearance indicates that classical SR cannot deal with large-parameter problems.
A similar conclusion can be obtained in Figure 6, when the amplitude A of the periodic signal is changed to 0.5, and the noise intensity D set to 7.75 (greater than 1). From the contrast result between Figure 6b,d, we know that the weak periodic signal cannot be enhanced by using classical SR theory when the noise intensity D is greater than 1. More information on why SR can only handle small-parameter problems can be obtained from the works of adiabatic approximation theory [19,20] and linear response theory [21].
Another limitation of the classical SR theory is that the sampling frequency fs should be greater than 50 times of the periodic signal frequency f, i.e.,
k = f s / f 50
where k is the ratio of the sampling frequency fs to the input frequency f.
A simulation is performed and the frequency spectrum of the SR output is shown in Figure 7, only increasing the frequency f to 0.2 and holding the other parameters in Figure 4. The spectral spike at f = 0.2 cannot be found completely. The requirement for high sampling frequency is also verified in [13].
Because of the above limitations, a novel method of LPSR named mixing and normalizing stochastic resonance (MNSR) is proposed in the next section. The proposed method makes it favorable to overcome the small parameters and high sampling frequency restriction. Therefore, the weak ULB signal detection can be dealt with by SR theory in practical application.

4. MNSR and Its Parameter Adjustment Method

4.1. Mixing and Normalizing Stochastic Ronance

This section presents a novel method of LPSR, mixing and normalizing stochastic resonance (MNSR). The method consists of three stages, as shown in Figure 8. First, the periodic signal immersed in noise is mixed with a carrier signal (the frequency of the carrier signal is fc). The output of the mixer is sm(t), which consists of the difference frequency signal component, the sum frequency periodic signal component, and the modulated noise. Secondly, a low-pass filter is used to filter sm(t). The pass-band cutoff frequency of the filter is fcut. The filtered output, m(t), only contains the difference frequency component and the noise. Finally, PNSR is used to detect the difference frequency component in m(t).
The carrier signal is cos(2πfct), the mixing output sm(t) is
s m ( t ) = A cos ( 2 π f t ) cos ( 2 π f c t ) + n ( t ) cos ( 2 π f c t ) = A d i f cos ( 2 π ( f c f ) t ) + A s u m cos ( 2 π ( f c + f ) t ) + n ( t ) cos ( 2 π f c t )
which consists of 3 components: the difference frequency signal component, the sum frequency periodic signal component, and the modulated noise. With a proper fc value, the difference frequency Δf = |fcf| can be changed much lower than the raw signal frequency f, i.e.,
Δ f = | f c f | < < f
To meet the limitation in Equation (12), MNSR requires the sampling frequency fs to be higher than 50·Δf, i.e., Equation (15). This requirement is much less than RFSR, which need fs to be higher than 50·f.
f s > 50 Δ f
Using a low-pass filter to remove the sum frequency periodic signal component, the filtered output can be written as
m ( t ) = A d i f cos ( 2 π Δ f t ) + n ( t )
where n′(t) is the filtered output of the modulated noise. From [15], the modulated noise n(t)·cos(2πfct) is a zero-mean, Gaussian white noise with noise intensity D/2. If the cutoff frequency fcut of the low-pass filter is much higher than the difference frequency Δf, i.e.,
Δ f < < f c u t < ( f c + f )
n′(t) can be regarded as a zero-mean, Gaussian white noise with noise intensity D′ in the low-frequency band. Because the modulated noise n(t)·cos(2πfct) is a white noise, the ratio of D/2 to D′ is equal to the ratio of fs/2 to fcut approximately, i.e.,
D / 2 D = f s / 2 f c u t
The intensity of n′(t), D′, can be calculated by Equation (18).
Then we can detect the weak difference frequency signal m(t) by PNSR method. With the input is m(t), the LE of MNSR is
d x d t = a x b x 3 + A d i f cos ( 2 π Δ f t ) + n ( t )
Placing a numerical transform to Equation (19), i.e.,
z = x b / a ,   τ = a t
Equation (7) can be rewritten as
a a b d z d τ = a a b z a a b z 3 + A d i f cos ( 2 π Δ f a τ ) + n ( τ a )
And Equation (21) can be further simplified as
d z d τ = z z 3 + b a 3 [ A d i f cos ( 2 π Δ f a τ ) + n ( τ a ) ]
n′(τ/a) can be regarded as the compressed form of n′(τ) in frequency domain, and n′(τ) is a zero-mean, Gaussian white noise with noise intensity D′ in the low-frequency band. Because n′(τ) is a white noise, the noise energy evenly distributed in frequency domain, so that the power of n(τ/a) and n(τ) is equal, i.e., the power is (σ′)2 = 2D′. Therefore, n′(τ/a) is still a zero-mean, Gaussian white noise with intensity D′,
n ( τ a ) = σ ξ ( τ ) = 2 D ξ ( τ )
where ξ(τ) denotes a zero-mean and unit-variance Gaussian white noise. Equation (22) can be rewritten as
d z d τ = z z 3 + A p cos ( 2 π f p τ ) + 2 D p ξ ( τ )
Equation (24) is another LE where the driving force is a periodic signal with amplitude Ap and frequency fp, and the noise is a zero-mean, Gaussian white noise with noise intensity Dp. There is
A P = A d i f b / a 3
f p = Δ f / a
D p = ( b / a 3 ) D
Compared with Equation (19), the barrier parameters of Equation (24) are changed to be constants, both are 1. The potential U(z) in Equation (24) is normalized as:
U ( z ) = 1 2 z 2 + 1 4 z 4
To make the transformed periodic signal frequency fp in Equation (24) less than 1, a large a is selected according to Equation (26). MNSR can handle the high-frequency signal detection problem when the difference signal frequency Δf in Equation (19) is higher than 1. To ensure there are at least two complete periods of the difference component in signal m(t), the signal duration Tsd should satisfy the following equation:
Δ f > 2 T s d
MNSR allows the signal duration Tsd to be shorter than MSR. Because the difference frequency Δf in MSR should be set to much lower than 1, the lower difference frequency Δf, the longer signal duration is required. On the other hand, the other barrier parameter b can be selected to ensure Dp is less than 1 according to Equation (27). No matter if the input frequency or noise strength or both are greater than 1, MNSR still works.

4.2. The Parament Adjustment Method

Based on the above discussion, we can conclude the parameters adjustment method of MNSR.
First, choosing the carrier signal frequency fc according to Equations (15) and (29), i.e.,
2 T s d < Δ f = | f c f | < f s / 50
As demonstrated in Equation (9), the amplitude of the response decreases monotonically with the input signal frequency. For this reason, the difference frequency Δf should be set as small as possible under the condition in Equation (30). In addition, then we can determine the pass-band cutoff frequency fcut by Equation (17).
Secondly, we determine a set of effective small parameters of MNSR in Equation (24). In this paper, we select the parameters: fp = 0.01, Dp = 0.31 in Figure 4.
Thirdly, we calculate the barrier parameters a, b from Equations (26) and (27). The computing way is:
a = Δ f / f p
b = a 3 ( D p / D )
where the value of D′ can be calculated by Equation (18) or by the mixer output variance. D′ is approximately equal to 1/2 of the variance of m(t) when the input SNR is low. The weak difference frequency signal can be found from the output of Equation (19), which is numerically solved with the fourth-order Runge–Kutta method.

4.3. An Example of MNSR

An example is cited in Figure 9. The parameters of input are: A = 0.5, f = 1000, D = 7.75. The sampling frequency fs = 5000, and the data length is 2000 points, i.e., the signal duration Tsd is 0.4 s. (a) and (b) are the waveform and frequency spectrum of the input. (c) and (d) are the waveform and frequency spectrum of the low-pass filter output, m(t). (e) and (f) are the waveform and frequency spectrum of the MNSR output. Based on the parameter adjustment method, we determine the carrier signal frequency fc at first. From Equation (30), the range of Δf is (5, 100). Δf is set to be 10, and fc is set to be 990. By Equation (15), the pass-band cutoff frequency of the filter fcut is 1500. Then, we choose fp = 0.01 and Dp = 0.31 in Figure 4 as the effective small parameters of MNSR in Equation (24). Therefore, the barrier parameters are calculated by Equations (31) and (32): a = Δf/fp = 10/0.01 = 1000; b = a3(Dp/D′) = a3[Dp/(D·fcut/fs)] = 1.45 × 108.
We can find that the signal frequency and noise intensity are greater than 1, and the ratio of fs to f is k = 5. The input cannot satisfy the two restrictions of classical SR, i.e., the input must be a small-parameter signal and k must be greater than 50. In Figure 9f, the frequency spectrum spike at Δf = 10 is evidently found. However, the height of the spike at f = 1000 in Figure 9b and the spike at Δf = 10 in Figure 9d are not prominent. Through applying MNSR, the difference frequency component in m(t) can be enhanced and detected. Therefore, the weak periodic signal in the input is found. This result shows the effectiveness of MNSR, even if the input is large-parameter signal and the sampling frequency is not high enough.

5. The Detection of the weak ULB Signal via MNSR

5.1. Numerical Simulation

Based on Section 2, the received signal of the ULB at a distance can be written as
s r e ( t ) = A cos ( 2 π f t ) + n ( t )
where f = 37.5 kHz, A = 0.15, the intensity D of the Gaussian white noise n(t) D = 0.9. The sampling frequency fs = 150 kHz, and the signal duration is the pulse width, Tsd = 0.01 s. The sampling data length is N = fs·Tsd = 1500. Obviously, the frequency is a large parameter and the ratio k (= fs/f = 4) is less than 50. In terms of the parameter adjustment method in Section 4.2, the difference frequency Δf should be set as small as possible in the range (2/Tpd, fs/50), i.e., (200, 3000). Because the frequency resolution is 1/Tpd = 100 Hz, Δf is set to 300 (200+100) Hz and the carrier frequency is fc = 37.8 kHz. Moreover, the cutoff frequency of the low-pass filter is set to fcut = 30 kHz. The filtered noise intensity D′ changed to 0.18, it is calculated by Equation (18). The barrier parameters: a = 300/0.01 = 3 × 104, b = a3(0.31/0.18) = 4.65 × 1013. They are calculated by Equations (31) and (32), respectively. The output of MNSR is numerically solved with a fourth-order Runge–Kutta method, shown in Figure 10.
There is a very distinct spectrum peak at Δf = 300 Hz in Figure 10d, but the spectrum spike at f = 37.5 kHz in Figure 10b are hardly distinguished. This appearance indicates that MNSR can be applied to detect and enhance the weak ULB signal. Even if the ULB signal frequency is much greater than 1 and the ratio k is far less than 50, MNSR is still effective. Those parameters of MNSR are determined by the parameter adjustment method proposed in Section 4.
Figure 11 shows the result of RFSR applied to the ULB signal detection. To meet the small-parameter limitation of SR, the re-scaling ratio of RFSR is set to R = 3.75 × 106. The re-scaling signal frequency is fr = f/R = 0.01, and the re-scaling sampling frequency is fsr = fs/R = 0.04. The other relative parameters are: A = 0.15, f = 37.5 kHz, a = 1, b = 0.344, D = 0.9. We can find that the numerical simulation result of RFSR is divergent. The ratio of the sampling frequency to the signal frequency k (= fsr/fr = 4) is less than 50. This low-frequency sampling cannot satisfy another limitation of SR, i.e., k > 50.
In Figure 12, the sampling frequency is changed to be fs = 500·f = 18.75 MHz. The sampling data length is still N = 1500, and the other parameters are the same as Figure 11. After the ratio is changed to be k = 500, we can find a distinct spectral spike at f = 37.5 kHz in Figure 12b. We know that RFSR can handle the large parameter SR problem only when the sampling frequency is more than 50 times to the signal frequency. The large ratio k requirement of RFSR is not exercisable: the quite huge data amount will produce challenges to the data storage, transmission, and processing. However, MNSR allows the sampling frequency to be much lower than RFSR.
Figure 13 is the output of MSR applied to the ULB signal detection. To meet the small-parameter limitation of SR, the carrier frequency of MSR is fc = 37.50001 kHz. The difference signal frequency is Δf = |fc − f| = 0.01 Hz, which is a small parameter. Meanwhile, the ratio k (= fsf = 150 × 103/0.01 = 1.5 × 105) is much greater than 50. The other relative parameters are: A = 0.15, f = 37.5 kHz, a = 1, b = 0.344, D = 0.9. Although Δf and k satisfy the two limitations of classical SR, we cannot find the spectral spike at Δf = 0.01 Hz in Figure 13b. This result is caused by the short duration of the ULB signal. In 10ms ULB pulse, a complete period of the difference signal cannot be contained, which is 100 s. MSR is not effective in ULB signal detection. However, MNSR allows a greater Δf and a shorter signal duration.

5.2. Experiment

The schematic diagram of the experiment in the anechoic tank is shown in Figure 14. We use a 37.5kHz acoustic source instead of a ULB. The reason for this is the SL of a ULB is fixed; it is hard to reduce the received SNR in anechoic tank. We can change the received SNR by changing the acoustic source SL. The sampling frequency of the data acquisition system is 150 kHz.
The source transmits the identical signal as a ULB signal. A piece of data in Figure 15a have the ULB signal known in advance, but it is hard to find the signal from Figure 15b.
Based on MNSR, the input noise intensity is estimated by variance, D = 0.22. The cutoff frequency fcut and the carrier frequency fc are 30 kHz and 37.8 kHz respectively, because input frequency and data length are the same as Figure 10.
The filtered noise intensity D′ changed to 0.044, which is calculated by Equation (18). Therefore, the barrier parameters are a = 3×104, and b = a3(Dp/D′) = a3(0.31/0.044) = 1.9 × 1014. A sharp peak in Figure 15d indicates that it has the ULB signal in this data piece. With an effective set of parameters (fc, a, b), this method MNSR can be applied in practical scenarios.

6. Conclusions

In view of the characteristics of high frequency and short pulse width of the ULB signal, a method named mixing and normalizing stochastic resonance (MNSR) is proposed in this paper. In addition, its parameter adjustment method is given as well. MNSR is a method of large-parameter stochastic resonance (LPSR) which can overcome the small-parameter limitation of classical SR. It has several advantages over conventional LPSR, such as RFSR and MSR. As shown in the simulations, RFSR requires higher sampling frequency and MSR requires longer signal duration time. Because we can adjust the carrier frequency fc and the barrier parameters a and b, MNSR is more exercisable, even if the signal sampling frequency is less than 50 times the signal frequency and the signal duration time is short. For a given received ULB signal, the valid parameters of MNSR can be calculated rapidly and easily by the proposed parameter adjustment method. Its effectiveness is verified by the experiment results. MNSR can be an active approach in the practical weak signal detection scenes, especially for ULB signal detection.

Author Contributions

Conceptualization, G.L. and G.W.; Formal analysis, G.W.; Funding acquisition, G.L.; Methodology, G.W.; Project administration, J.W.; Resources, J.W. and X.W.; Supervision, J.W. and X.W.; Writing—original draft, G.W.; Writing—review & editing, G.L., G.W., and J.W. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported by the National Key Research and Development Plan (2017YFC0306900), the National Natural Science Foundation of China (61801138), Qingdao National Laboratory for Marine Science and Technology Open Found (QNLM2016ORP0102), and the Acoustic Science and Technology Laboratory Foundation.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Stone, L.D.; Keller, C.M.; Kratzke, T.M.; Strumpfer, J.P. Search analysis for the underwater wreckage of Air France Flight 447. In Proceedings of the 14th International Conference on Information Fusion, Chicago, IL, USA, 5–8 July 2011; pp. 1–8. [Google Scholar]
  2. Pranjoto, H. Innovative Tester for Underwater Locator Beacon Used in Flight/Voyage Recorder (Black Box). In Proceedings of the Second International Conference on Electrical Systems, Technology and Information 2015 (ICESTI 2015), Singapore, 9–12 September 2016; pp. 299–307. [Google Scholar]
  3. Xiao, Y.; Tadokoro, Y. LMS-based notch filter for the estimation of sinusoidal signals in noise. Signal Process. 1995, 46, 223–231. [Google Scholar] [CrossRef]
  4. Amarnath, M.; Praveen Krishna, I.R. Local fault detection in helical gears via vibration and acoustic signals using EMD based statistical parameter analysis. Measurement 2014, 58, 154–164. [Google Scholar] [CrossRef]
  5. Ioana, C.; Quinquis, A. Transient signal detection using overcomplete wavelet transform and high-order statistics. In Proceedings of the 2003 IEEE International Conference on Acoustics, Speech, and Signal Processing, Hong Kong, China, 6–10 April 2003. [Google Scholar]
  6. Benzi, R.; Sutera, A.; Vulpiani, A. The mechanism of stochastic resonance. J. Phys. A Math. Gen. 1981, 14, L453–L457. [Google Scholar] [CrossRef]
  7. Collins, J.J.; Chow, C.C.; Imhoff, T.T. Aperiodic stochastic resonance in excitable systems. Phys. Rev. E 1995, 52, R3321. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Wiesenfeld, K.; Moss, F. Stochastic resonance and the benefits of noise: From ice ages to crayfish and SQUIDs. Nature 1995, 373, 33–36. [Google Scholar] [CrossRef] [PubMed]
  9. Liu, Q.; Wang, M. Application of Stochastic Resonance in Directions of Arrival. In Proceedings of the 2010 2nd International Conference on Information Engineering and Computer Science, Wuhan, China, 25–26 December 2010; pp. 1–3. [Google Scholar]
  10. Fauve, S.; Heslot, F. Stochastic resonance in a bistable system. Phys. Lett. A 1983, 97, 5–7. [Google Scholar] [CrossRef]
  11. Leng, Y.G.; Wang, T.Y. Numerical research of twice sampling stochastic resonance for the detection of a weak signal submerged in a heavy Noise. Acta Phys. Sin. 2003, 52, 2432. [Google Scholar] [CrossRef]
  12. Lu, S.; He, Q.; Kong, F. Stochastic resonance with Woods–Saxon potential for rolling element bearing fault diagnosis. Mech. Syst. Signal Process. 2014, 45, 488–503. [Google Scholar] [CrossRef]
  13. Leng, Y.G.; Leng, Y.S.; Wang, T.Y.; Guo, Y. Numerical analysis and engineering application of large parameter stochastic resonance. J. Sound Vib. 2006, 292, 788–801. [Google Scholar] [CrossRef]
  14. Yang, D.X.; Hu, Z.; Yang, Y.M. The analysis of stochastic resonance of periodic signal with large parameters. Acta Phys. Sin. 2012, 61, 80501. [Google Scholar] [CrossRef]
  15. Min, L.; Yong-Mei, H. Modulation and demodulation for detecting weak periodic signal of stochastic resonance. Acta Phys. Sin. 2006, 55, 3277. [Google Scholar] [CrossRef]
  16. Tan, J.; Chen, X.; Wang, J.; Chen, H.; Cao, H.; Zi, Y.; He, Z. Study of frequency-shifted and re-scaling stochastic resonance and its application to fault diagnosis. Mech. Syst. Signal Process. 2009, 23, 811–822. [Google Scholar] [CrossRef]
  17. Liu, J.; Leng, Y.; Lai, Z.; Fan, S. Multi-Frequency Signal Detection Based on Frequency Exchange and Re-Scaling Stochastic Resonance and Its Application to Weak Fault Diagnosis. Sensors 2018, 18, 1325. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Urick, R.J. Principles of Underwater Sound, 3rd ed.; Peninsula Publising: Los Atlos, CA, USA, 1983. [Google Scholar]
  19. Nicolis, C.; Nicolis, G. Stochastic aspects of climatic transitions–additive fluctuations. Tellus 1981, 33, 225–234. [Google Scholar] [CrossRef] [Green Version]
  20. McNamara, B.; Wiesenfeld, K. Theory of stochastic resonance. Phys. Rev. A 1989, 39, 4854–4869. [Google Scholar] [CrossRef] [PubMed]
  21. Dykman, M.; Luchinsky, D.; Mannella, R.; McClintock, P.V.; Stein, N.; Stocks, N. Stochastic resonance: Linear response and giant nonlinearity. J. Stat. Phys. 1993, 70, 463–478. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The variation of the received ULB signal SNR at different distances.
Figure 1. The variation of the received ULB signal SNR at different distances.
Sensors 20 01292 g001
Figure 2. The received signal at 2 km distance. (a) waveform of the received signal; (b) frequency spectrum of the received signal.
Figure 2. The received signal at 2 km distance. (a) waveform of the received signal; (b) frequency spectrum of the received signal.
Sensors 20 01292 g002
Figure 3. The output SNR versus noise intensity D.
Figure 3. The output SNR versus noise intensity D.
Sensors 20 01292 g003
Figure 4. A simulation of classical SR. The relative parameters corresponding to Equation (7) configured as: A = 0.1, f = 0.01, a = 1, b = 1, D = 0.31, fs = 5. (a,b) are the waveform and frequency spectrum of the input; (c,d) are waveform and frequency spectrum of the SR output.
Figure 4. A simulation of classical SR. The relative parameters corresponding to Equation (7) configured as: A = 0.1, f = 0.01, a = 1, b = 1, D = 0.31, fs = 5. (a,b) are the waveform and frequency spectrum of the input; (c,d) are waveform and frequency spectrum of the SR output.
Sensors 20 01292 g004
Figure 5. The SR of high-frequency signal. The relative parameters: A = 0.1, f = 10, a = 1, b = 1, D = 0.31, Figure 5. (a,b) are the waveform and frequency spectrum of the input; (c,d) are waveform and frequency spectrum of the SR output.
Figure 5. The SR of high-frequency signal. The relative parameters: A = 0.1, f = 10, a = 1, b = 1, D = 0.31, Figure 5. (a,b) are the waveform and frequency spectrum of the input; (c,d) are waveform and frequency spectrum of the SR output.
Sensors 20 01292 g005
Figure 6. The SR of signal with large noise. The relative parameters: A = 0.5, f = 0.01, a = 1, b = 1, D = 7.75, fs = 5. (a,b) are the waveform and frequency spectrum of the input; (c,d) are waveform and frequency spectrum of the SR output.
Figure 6. The SR of signal with large noise. The relative parameters: A = 0.5, f = 0.01, a = 1, b = 1, D = 7.75, fs = 5. (a,b) are the waveform and frequency spectrum of the input; (c,d) are waveform and frequency spectrum of the SR output.
Sensors 20 01292 g006
Figure 7. The frequency spectrum of the SR output when the ratio k is small. The relative parameters: A = 0.1, f = 0.2, a = 1, b = 1, D = 0.31, fs = 5.
Figure 7. The frequency spectrum of the SR output when the ratio k is small. The relative parameters: A = 0.1, f = 0.2, a = 1, b = 1, D = 0.31, fs = 5.
Sensors 20 01292 g007
Figure 8. The flow of MNSR.
Figure 8. The flow of MNSR.
Sensors 20 01292 g008
Figure 9. An example shown the effectivity of MNSR. The relative parameters: A = 0.5, f = 1000, a = 1000, b = 1.45 × 108, fc = 990 Hz, D = 7.75, fs = 5000. (a,b) are the waveform and frequency spectrum of the input. (c,d) are the waveform and frequency spectrum of the low-pass filter output, m(t). (e,f) are the waveform and frequency spectrum of the MNSR output.
Figure 9. An example shown the effectivity of MNSR. The relative parameters: A = 0.5, f = 1000, a = 1000, b = 1.45 × 108, fc = 990 Hz, D = 7.75, fs = 5000. (a,b) are the waveform and frequency spectrum of the input. (c,d) are the waveform and frequency spectrum of the low-pass filter output, m(t). (e,f) are the waveform and frequency spectrum of the MNSR output.
Sensors 20 01292 g009aSensors 20 01292 g009b
Figure 10. The numerical simulation results of MNSR. The relative parameters: A = 0.15, f = 37.5 kHz, a = 3 × 104, b = 4.65 × 1013, fc = 37.8 kHz, D = 0.9, fs = 150 kHz, N = 1500. (a,b) are the waveform and frequency spectrum of the input; (c,d) are waveform and frequency spectrum of the MNSR output.
Figure 10. The numerical simulation results of MNSR. The relative parameters: A = 0.15, f = 37.5 kHz, a = 3 × 104, b = 4.65 × 1013, fc = 37.8 kHz, D = 0.9, fs = 150 kHz, N = 1500. (a,b) are the waveform and frequency spectrum of the input; (c,d) are waveform and frequency spectrum of the MNSR output.
Sensors 20 01292 g010
Figure 11. The result of RFSR applied to the ULB signal detection when the sampling frequency is fs = 150 kHz, and the sampling data length is N = 1500. The re-scaling ratio of RFSR is R = 3.75 × 106. In addition, the other relative parameters: A = 0.15, f = 37.5 kHz, a = 1, b = 0.344, D = 0.9. (a,b) are the waveform and frequency spectrum of the output.
Figure 11. The result of RFSR applied to the ULB signal detection when the sampling frequency is fs = 150 kHz, and the sampling data length is N = 1500. The re-scaling ratio of RFSR is R = 3.75 × 106. In addition, the other relative parameters: A = 0.15, f = 37.5 kHz, a = 1, b = 0.344, D = 0.9. (a,b) are the waveform and frequency spectrum of the output.
Sensors 20 01292 g011
Figure 12. The result of RFSR applied to the ULB signal detection when the sampling frequency is fs = 18.75 MHz. The sampling data length is still N = 1500 and the other parameters are the same as Figure 11.
Figure 12. The result of RFSR applied to the ULB signal detection when the sampling frequency is fs = 18.75 MHz. The sampling data length is still N = 1500 and the other parameters are the same as Figure 11.
Sensors 20 01292 g012
Figure 13. The result of MSR applied to the ULB signal detection when the sampling frequency is fs = 150 kHz, and the sampling data length is N = 1500. The carrier frequency is fc = 37.50001 kHz. In addition, the other relative parameters: A = 0.15, f = 37.5 kHz, a = 1, b = 0.344, D = 0.9. (a,b) are the waveform and frequency spectrum of the output.
Figure 13. The result of MSR applied to the ULB signal detection when the sampling frequency is fs = 150 kHz, and the sampling data length is N = 1500. The carrier frequency is fc = 37.50001 kHz. In addition, the other relative parameters: A = 0.15, f = 37.5 kHz, a = 1, b = 0.344, D = 0.9. (a,b) are the waveform and frequency spectrum of the output.
Sensors 20 01292 g013
Figure 14. The schematic diagram of experiment.
Figure 14. The schematic diagram of experiment.
Sensors 20 01292 g014
Figure 15. The input and output of a tank experiment where MNSR used to detect the ULB signal. The relative parameters: f = 37.5 kHz, a = 3 × 104, b = 1.9 × 1014, fc = 37.8 kHz, D = 0.22, fs = 150 kHz. (a,b) are the waveform and frequency spectrum of the input; (c,d) are waveform and frequency spectrum of the MNSR output.
Figure 15. The input and output of a tank experiment where MNSR used to detect the ULB signal. The relative parameters: f = 37.5 kHz, a = 3 × 104, b = 1.9 × 1014, fc = 37.8 kHz, D = 0.22, fs = 150 kHz. (a,b) are the waveform and frequency spectrum of the input; (c,d) are waveform and frequency spectrum of the MNSR output.
Sensors 20 01292 g015

Share and Cite

MDPI and ACS Style

Liang, G.; Wan, G.; Wang, J.; Wang, X. A Novel Underwater Location Beacon Signal Detection Method Based on Mixing and Normalizing Stochastic Resonance. Sensors 2020, 20, 1292. https://doi.org/10.3390/s20051292

AMA Style

Liang G, Wan G, Wang J, Wang X. A Novel Underwater Location Beacon Signal Detection Method Based on Mixing and Normalizing Stochastic Resonance. Sensors. 2020; 20(5):1292. https://doi.org/10.3390/s20051292

Chicago/Turabian Style

Liang, Guolong, Guangming Wan, Jinjin Wang, and Xue Wang. 2020. "A Novel Underwater Location Beacon Signal Detection Method Based on Mixing and Normalizing Stochastic Resonance" Sensors 20, no. 5: 1292. https://doi.org/10.3390/s20051292

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop