Research Article
Open Access

Source Localization Technique for
Multiple Acoustic Sources by Using
Acoustic Beating Envelopes

Kyoung-Sik Choi

^{*}, Sung-Jae Heo, Tae-Wan Kim
Korea Aerospace Industries Ltd., Sacheon-City, Gyeongnam, 664-710, South Korea

***Corresponding author:**Kyoung-Sik Choi, Korea Aerospace Industries Ltd., Sacheon-City, Gyeongnam, 664-710, South Korea,E-mail:

Received: March 31, 2017; Accepted: April 11, 2017; Published: April 23, 2017

**Citation:**Kyoung-Sik Choi, Sung-Jae Heo, Tae-Wan Kim (2017) Source Localization Technique for Multiple Acoustic Sources by Using Acoustic Beating Envelopes. int J Adv Robot Automn 2(1):1-7. DOI: 10.15226/2473-3032/2/1/00122

Abstract

We propose here new technique to detect multiple acoustic
source positions. This technique is more simple and cost saving method.
We used several acoustic sensors to detect the exact positions of
multiple acoustic sources simultaneously. Especially acoustic waves
higher than audio frequency were used to avoid unnecessary acoustic
noises. In this case, higher frequency means short wave length, which
makes it difficult to extract the phase difference between two sensors
when the distance of the two sensors exceed half of wave length. In
this paper, we applied beating signal which is composed of two slightly
different close frequencies. The wave length of the beating signal’s envelope
which is longer than those of the orignals is controllable. Using
the envelope of beating signal gives us many advantages. For example,
we are free from unnecessary acoustic noises because of using higher
frequency than audio frequency. Additionally, the wave length of beating
signal is under control which means sensor position is not limited
for the signal processing. We applied here a hypothetical simulation to
prove it. by some final considerations.

**Keywords:**Source localization; Beating signal; Envelope; Phase delay; Acoustic wave; TDOA;Introduction

Source localization techniques have been proposed in
many areas for the various applications [1-3]. We proposed here
more simple and cost saving method. The main key idea is using
beating signal which is composed of two slightly different close
frequencies. Especially the frequencies of two waves are higher
than audio frequency. The reason why we use high frequency is
low frequency (audio frequency) yields unnecessary acoustic
noises and most of the mechanical noises are ranging under
audio frequency. Mechanical noises could yield interferences.
Because of these two critical issues, we determined to use high
frequency. But high frequency means short wave length which
makes a difficulty to extract exact phase difference when the
distance of the two sensors is exceeding half of wave length. In
order to solve this problem we used beating signal. Beating signal
is composed of two slightly different close signals [4]. These two
close waves make beating signal as a result of interference. The
wave length of a beating signal’s envelope is longer than those
of the originals. It means we can extract longer wave from high
frequency waves without any transform. This gives us many
advantages. For example, we are free from choosing frequency
range. Even though when we have to deal with high frequency,
we can extract longer wave from original signals without any
transform. For these reasons we used beating signal to estimate
multiple acoustic sources. In order to estimate source location
we have to extract arrival time difference of the received sensor
signals respectively. Time differences could be extracted from
phase differences. In this paper, we extracted beating signal’s
phase difference from original signals. We showed how to find a
beating signal frequencies and phases and how to extract phase
differences among beating signals. Finally showed how these
phase differences are used finding source locations.

Fundamental Theories

Beating Signalm

Beating signal is composed of two slightly different
close waves [4]. The wave length of a beating signal’s envelope
is longer than those of the originals. This allows us to use high
frequency.

When we estimate source location by using phase differences of received signals, we would confront some limitations. The wave length should be two times or longer than the distance between two sensors. For the reason that a wave length of a beating signal’s envelope is much longer than those of the originals, we could put the sensors far away each other. That means sensors’ locations are free for the wave length because we are able to control wave length of a beating signal. We can make beating signal by using frequency mixer or simply summing slightly different two frequencies. In case of frequency mixer (Figure 1), put the high frequency as ${\omega}_{h}$and low frequency as ${\omega}_{l}$ respectively. These two signals satisfy the conditions below. $${\omega}_{h}\gg {\omega}_{l}(abs\left({\omega}_{h}-{\omega}_{l}\right)20KHz)\text{(1)}$$

When we estimate source location by using phase differences of received signals, we would confront some limitations. The wave length should be two times or longer than the distance between two sensors. For the reason that a wave length of a beating signal’s envelope is much longer than those of the originals, we could put the sensors far away each other. That means sensors’ locations are free for the wave length because we are able to control wave length of a beating signal. We can make beating signal by using frequency mixer or simply summing slightly different two frequencies. In case of frequency mixer (Figure 1), put the high frequency as ${\omega}_{h}$and low frequency as ${\omega}_{l}$ respectively. These two signals satisfy the conditions below. $${\omega}_{h}\gg {\omega}_{l}(abs\left({\omega}_{h}-{\omega}_{l}\right)20KHz)\text{(1)}$$

**Figure 1:**Frequency mixer can get beating signal

For the conditions above,
${\omega}_{h}$should be much higher than
${\omega}_{l}$Then we can make beating signal with these two frequencies by using frequency mixer. This can be expressed as below

$$\mathrm{cos}\left({\omega}_{h}t\right)\times \mathrm{cos}\left({\omega}_{l}t\right)=\frac{1}{2}\mathrm{cos}(({\omega}_{h}+{\omega}_{l})t)+\frac{1}{2}\mathrm{cos}(({\omega}_{h}-{\omega}_{l})t)\text{}\text{}\text{(2)}$$

In this case, each wave length of ${\omega}_{h}-{\omega}_{l}$and ${\omega}_{h}+{\omega}_{l}$are slightly different. As a result, these two frequencies make beating signal.

The other way to make beating signal is using two frequencies. This is a simpler way to make beating signal. That directly summates slightly different two waves $({\omega}_{1},{\omega}_{2})$ . We can assume that each wave $({\omega}_{1},{\omega}_{2})$ is composed of two frequencies of ${\omega}_{h},{\omega}_{l}$ These two frequencies might satisfy Eq.(1) Then we can express ${\omega}_{1},{\omega}_{2}$as below. $${\omega}_{1}={\omega}_{h}+{\omega}_{l},{\omega}_{2}={\omega}_{h}-{\omega}_{l}\text{(3)}$$

Then we can get slightly different two close waves ${\omega}_{1},{\omega}_{2}$ . Figure 2 shows how to make beating signal.

$$\mathrm{cos}\left({\omega}_{h}t\right)\times \mathrm{cos}\left({\omega}_{l}t\right)=\frac{1}{2}\mathrm{cos}(({\omega}_{h}+{\omega}_{l})t)+\frac{1}{2}\mathrm{cos}(({\omega}_{h}-{\omega}_{l})t)\text{}\text{}\text{(2)}$$

In this case, each wave length of ${\omega}_{h}-{\omega}_{l}$and ${\omega}_{h}+{\omega}_{l}$are slightly different. As a result, these two frequencies make beating signal.

The other way to make beating signal is using two frequencies. This is a simpler way to make beating signal. That directly summates slightly different two waves $({\omega}_{1},{\omega}_{2})$ . We can assume that each wave $({\omega}_{1},{\omega}_{2})$ is composed of two frequencies of ${\omega}_{h},{\omega}_{l}$ These two frequencies might satisfy Eq.(1) Then we can express ${\omega}_{1},{\omega}_{2}$as below. $${\omega}_{1}={\omega}_{h}+{\omega}_{l},{\omega}_{2}={\omega}_{h}-{\omega}_{l}\text{(3)}$$

Then we can get slightly different two close waves ${\omega}_{1},{\omega}_{2}$ . Figure 2 shows how to make beating signal.

**Figure 2:**One beating signal is composed of slightly different two close waves

Envelope of Beating Signal

The next step is extracting frequency and phase of a beating signal. One method of extracting the beating signal’s envelope is to use the Hilbert Transform [8,9]. The Hilbert Transform of a real-valued time domain signal
$\text{x}(\text{t})$is another real-valued time domain signal denoted by
$\hat{x}(t)$, such that
$z\left(t\right)=x\left(t\right)+j\hat{x}(t)$
. The Hilbert transform of
$\text{x}(\text{t})$extending over the range
$-\infty <t<\infty $
is a real-valued function
$\hat{x}(t)$defined by
$$\hat{x}\left(t\right)=H\left\{x\left(t\right)\right\}=\underset{-\infty}{\overset{\infty}{{\displaystyle \int}}}\frac{x(\tau )}{\pi (t-\tau )}d\tau \text{(4)}$$
$\hat{x}(t)$is the convolution integral of
$x(t)$and
$\raisebox{1ex}{$1$}\!\left/ \!\raisebox{-1ex}{$\pi t$}\right.$
, written as
$$\hat{x}\left(t\right)=x\left(t\right)*(\raisebox{1ex}{$1$}\!\left/ \!\raisebox{-1ex}{$\pi t$}\right.)\text{(5)}$$
A useful way to compute the Hilbert Transform
$\hat{x}\left(t\right)$of
$x\left(t\right)$ is via the analytic signal
$\text{z}\left(\text{t}\right)=x\left(t\right)+\text{j}\hat{x}(t)$
. One can also write
$\text{z}\left(\text{t}\right)=A\left(t\right){e}^{j\theta (t)}$,where $\text{A}(\text{t})$
denotes the envelope signal of
$x\left(t\right)$ and
$\text{\theta}(\text{t})$
denotes the instantaneous phase signal of
$x\left(t\right)$
. In terms of
$x\left(t\right)$and
$\hat{x}(t)$
, it is clear that
$$\text{A}\left(\text{t}\right)={[{x}^{2}\left(\text{t}\right)+{\hat{x}}^{2}(\text{t})]}^{1/2}\text{(6)}$$
$$\text{\theta}\left(\text{t}\right)=tan{}^{-1}\left[\frac{\hat{x}\left(t\right)}{x\left(t\right)}\right]$$So, we can obtain envelopes by using Eq.(6).

We showed here how we can get beating signal’s envelop by using Hilbert Transform. But in many cases, received signals have lots of noises and it is not always limited to one sources. If there are many sources then it is not easy to find beating signal’s envelope. So we introduced here more general method to find frequency and phase of a beating signal. Let’s put two frequencies as
${\omega}_{1}(={\omega}_{h}+{\omega}_{l}),$
${\omega}_{1}(={\omega}_{h}+{\omega}_{l}),$and its correspoding phases
$\text{\phi}({\omega}_{1})$and its correspoding phases
$\text{\phi}({\omega}_{1})$
$\text{\phi}\left({\omega}_{2}\right)$respectively.These two waves are satifying the conditions of Eq.(1) and Eq.(3). Then beating signal’s frequency
${\omega}_{beat}$And its corresponding phase
$\text{\phi}({\omega}_{beat})$can be expressed as below
$${\text{\omega}}_{\text{beat}}={\text{\omega}}_{1}-{\text{\omega}}_{2}({\text{\omega}}_{1}>{\text{\omega}}_{2})\text{(7)}$$
$${\text{\phi}}_{\text{beat}}=\text{\phi}\left({\text{\omega}}_{1}\right)-\text{\phi}\left({\text{\omega}}_{2}\right)$$We can express a beating sginal from Eq.(2) as below
$$\mathrm{cos}\left({\omega}_{h}(t-{\tau}_{h})\right)\times \mathrm{cos}\left({\omega}_{l}(t-{\tau}_{l})\right)=1/2\mathrm{cos}(({\omega}_{h}+{\omega}_{l})t-({\omega}_{h}{\tau}_{h}+{\omega}_{l}{\tau}_{l}))+1/2\mathrm{cos}(({\omega}_{h}-{\omega}_{l})t-({\omega}_{h}{\tau}_{h}-{\omega}_{l}{\tau}_{l}))\text{(8)}$$Then we can extract beating signal’s frequency and phase using Eq.(7) and Eq.(8).
(Figure 3 shows beating signal and it’s envelope)

**Figure 3:**Envelope of beating signal (a) beating signal, (b) envelope for the beating signaleating signal

$${\omega}_{beating}={\omega}_{1}-{\omega}_{2}=2{\omega}_{l}\text{(9)}$$
$${\phi}_{beating}=\phi \left({\omega}_{1}\right)-\phi \left({\omega}_{2}\right)=2{\omega}_{l}{\tau}_{l}$$

Frequency and Phase of a Beating Signal

In order to estimate an acoustic source locations, the time delays between acoutic sensors signals must be calculated. The time delay is related to the phase delay of the envelopes. But as shown in Figure 4,

**Figure 4:**Time delay of a beating signal’s envelope

it is difficult to find the arrival time difference using the envelope in the time domain. So, we tried to find the time delays using the envelope’s phase delay in the frequency domain [7-9]. We are dealing with two beating signals
$p(t)$and
$q(t)$.We could think
$p(t)$and
$q(t)$are received signals at a different two sensors, and
${\tau}_{d}$means the time delay between
$p(t)$and
$q(t)$.Then we can express it as
$$p\left(t\right)=A\{\mathrm{cos}\left({\omega}_{1}t-{\phi}_{1}\right)+cos\left({\omega}_{2}t-{\text{\phi}}_{2}\right)\}\text{(10)}$$

In this case, the frequency and phase of beating signal $p\left(t\right)$are ${\text{\omega}}_{1}-{\text{\omega}}_{2}\left({\text{\omega}}_{1}>{\text{\omega}}_{2}\right),$ ${\text{\phi}}_{1}-{\text{\phi}}_{2}$respectively. Then, $$q\left(t\right)=\alpha \times p(t-{\tau}_{d})=\alpha \times A\{\mathrm{cos}\left({\omega}_{1}\left(t-{\tau}_{d}\right)-{\phi}_{1}\right)+cos\left({\omega}_{2}\left(t-{\tau}_{d}\right)-{\text{\phi}}_{2}\right)\}\text{(11)}$$ where ${\text{\omega}}_{1,}{\text{\omega}}_{2}$are the frequency components composing beating signal and α is the magnitude ratio between $p\left(t\right)$and $q\left(t\right)$The cross-correlation of Eqs. (10) and (11) is derived as the following equation. $${R}_{p,q}\left(\tau \right)=\alpha E[p\left(t\right)q(t+\tau )]=\alpha {R}_{p,p}\left(\tau -{\tau}_{d}\right)\text{(12)}$$where E[a] is the expected value of a. This equation shows that the cross-correlation of two signals can be expressed as the autocorrelation that has a time delay τ

In this case, the frequency and phase of beating signal $p\left(t\right)$are ${\text{\omega}}_{1}-{\text{\omega}}_{2}\left({\text{\omega}}_{1}>{\text{\omega}}_{2}\right),$ ${\text{\phi}}_{1}-{\text{\phi}}_{2}$respectively. Then, $$q\left(t\right)=\alpha \times p(t-{\tau}_{d})=\alpha \times A\{\mathrm{cos}\left({\omega}_{1}\left(t-{\tau}_{d}\right)-{\phi}_{1}\right)+cos\left({\omega}_{2}\left(t-{\tau}_{d}\right)-{\text{\phi}}_{2}\right)\}\text{(11)}$$ where ${\text{\omega}}_{1,}{\text{\omega}}_{2}$are the frequency components composing beating signal and α is the magnitude ratio between $p\left(t\right)$and $q\left(t\right)$The cross-correlation of Eqs. (10) and (11) is derived as the following equation. $${R}_{p,q}\left(\tau \right)=\alpha E[p\left(t\right)q(t+\tau )]=\alpha {R}_{p,p}\left(\tau -{\tau}_{d}\right)\text{(12)}$$where E[a] is the expected value of a. This equation shows that the cross-correlation of two signals can be expressed as the autocorrelation that has a time delay τ

_{d}. If Eq. (12) is transformed to the frequency domain, we can obtain the cross-spectrum. Let $P(\omega )$and $\text{Q}(\omega )$denote the components of the Fourier transform of $p\left(t\right)$ and $q\left(t\right)$;the cross-spectrum is $${S}_{p,q}\left(\omega \right)={P}^{*}\left(\omega \right)\times Q\left(\omega \right)=\left|P\left(\omega \right)\right||Q(\omega )|{e}^{j\{{\text{\phi}}_{p}(\omega )-{\text{\phi}}_{q}\left(\omega \right)\}}\text{(13)}$$where ${\text{\phi}}_{p}(\omega )$and ${\text{\phi}}_{q}(\omega )$ mean the phase of $p\left(t\right)$and $q\left(t\right)$ respectively. Figure 4 shows two beating signals having some phase differences.TDOA(Time Difference Of Arrival)

In order to extract Time-Difference-Of-Arrival (TDOA) between two received beating signals. We should compare the arrived time difference of the two received beating signals. Put received beating signal’s phases as
${\phi}_{p},{\phi}_{q}$respectively then phase difference between two beating signals could be expressed as
$${\phi}_{p}-{\phi}_{q}={\phi}_{p}\left({\omega}_{1}\right)-{\phi}_{p}\left({\omega}_{2}\right)-({\phi}_{q}\left({\omega}_{1}\right)-{\phi}_{q}\left({\omega}_{2}\right))={\phi}_{p}\left({\omega}_{1}\right)-{\phi}_{q}\left({\omega}_{1}\right)-({\phi}_{p}\left({\omega}_{2}\right)-{\phi}_{q}({\omega}_{2}))\text{(14)}$$Here, we could get some conclusion. Phase difference between two beating signals is a subtraction of a same frequency. This means we can get beating signal’s phase differences by crossing spectrum. It is not necessary to find beating signal envelopes. The relationship between period
$\text{'T'}$and frequency
$\text{'}f\text{'}$is
$\text{T}=\raisebox{1ex}{$1$}\!\left/ \!\raisebox{-1ex}{$f$}\right.$The time delay of acoustic beating envelope
$\Delta {t}_{p,q}(f)$can be obtained by multiplying the period by the ratio of phase delay. Therefore, the time delay can be expressed as following equation.
$$\Delta {t}_{p,q}(f)=\frac{1}{f}\times \frac{{\phi}_{p}\left(f\right)-{\phi}_{q}(f)}{2\pi}\text{(15)}$$The distance difference transferred from the impact source to the two sensors is derived by multiplying Eq.(11) by the sound velocity.
$(v=340\text{m}/\text{s,}@15\mathbb{C})$

$$\Delta {r}_{p,q}={r}_{p}-{r}_{q}=v\times \Delta {t}_{p,q}\left(f\right)\text{(16)}$$

$$\Delta {r}_{p,q}={r}_{p}-{r}_{q}=v\times \Delta {t}_{p,q}\left(f\right)\text{(16)}$$

Triangular Method

Triangulation method was applied to point out an exact position. Figure 5 shows the triangular method for only two sensors. And figure 6 shows for more than three sensors.

In figure 7 Triangulation method looks for the hyperbola intersection point based on time-of-arrival differences (TOADs) of a burst between three different sensors [3,5,6]. Assume sensor number is N, then we can get hyperbola lines number of
$$M=\frac{N!}{2!\left(N-2\right)!}=\frac{N(N-1)}{2}\text{(17)}$$

**Figure 5:**Hyperbola line when only two sensors are used

**Figure 6:**Hyperbola line when more than 3 sensors are used

**Figure 7:**Sample of software simulation for triangulation method (a) only two sensors are used, (b) more than 3 sensors are used

Simulation and Verification

In order to verify our proposed theory by experiment, we put two acoustic emission sources and three sensors. Figure 8 shows the conditions of a hypothetical simulation.

The coordinates of two acoustic emission sources are E1(250,350) and E2(350,250). The coordinates of three acoustic sensors are S1(0,0), S2(600,0) and S3(300,600) respectively. We can derive two beating signals from the frequencies of the two acoustic emission sources. The first beating signal is composed of 25KHz, 25.5KHz, and the second beating signal is composed of 27KHz and 27.5KHz. Figure 9 shows received signals from S1, S2 and S3. Figure 10 shows the result of cross power spectrum of the received signals. In order to find phase differences between beating signals depending on Eq.(14) we should find phase differences of the same frequencies first. The phase differences of the same frequencies for the first beating signal (E1) are as below

The coordinates of two acoustic emission sources are E1(250,350) and E2(350,250). The coordinates of three acoustic sensors are S1(0,0), S2(600,0) and S3(300,600) respectively. We can derive two beating signals from the frequencies of the two acoustic emission sources. The first beating signal is composed of 25KHz, 25.5KHz, and the second beating signal is composed of 27KHz and 27.5KHz. Figure 9 shows received signals from S1, S2 and S3. Figure 10 shows the result of cross power spectrum of the received signals. In order to find phase differences between beating signals depending on Eq.(14) we should find phase differences of the same frequencies first. The phase differences of the same frequencies for the first beating signal (E1) are as below

**Figure 8:**Simulation conditions which motion box area is 600mm x 600mm, red circle denotes the emission sources

**Figure 9:**Measured signals for each acoustic sensors, (a) S1, (b) S2, (c) S3

**Figure 10:**The result of cross power spectrum (a) S1&S2 (b) S1&S3, (c) S2&S3. Blue line denote magnitude, red line denote phase

${\phi}_{1,2}\left(25KHz\right)=0.594,$

${\phi}_{1,3}\left(25KHz\right)=0.0,$

${\phi}_{2,3}\left(25KHz\right)=-0.594,$

${\phi}_{1,2}\left(25.5KHz\right)=0.0,$

${\phi}_{1,3}\left(25.5KHz\right)=0.0,$

${\phi}_{2,3}\left(25.5KHz\right)=1.615$

and then we can extract phase difference for the two beating signals by using Eq.(14)

${\phi}_{1,2}\left(25.5KHz\right)-{\phi}_{1,2}\left(25KHz\right)=-0.594$

${\phi}_{1,3}\left(25.5KHz\right)-{\phi}_{1,3}\left(25KHz\right)=1.615$

${\phi}_{2,3}\left(25.5Hz\right)-{\phi}_{2,3}\left(25KHz\right)=2.209$

Time delays for the first beating signal is as below

$\Delta {t}_{1,2}\left(0.5KHz\right)=-0.000189$

$\Delta {t}_{1,3}\left(0.5KHz\right)=0.000514$

$\Delta {t}_{2,3}\left(0.5KHz\right)=0.000703$

Then the delay distance can be obtained by multiplying sound velocity.

$\Delta {r}_{1,2}=-0.06426\text{m}$

$\Delta {r}_{1,3}=0.17476\text{m}$

$\Delta {r}_{2,3}=0.23902\text{m}$

The second beating signals are as below(E2)

${\phi}_{1,2}\left(27KHz\right)=0.0$

${\phi}_{1,3}\left(27KHz\right)=0.0$

${\phi}_{2,3}\left(27KHz\right)=0.0$

${\phi}_{1,2}\left(27.5KHz\right)=0.701$

${\phi}_{1,3}\left(27.5KHz\right)=0.713$

${\phi}_{2,3}\left(27.5KHz\right)=0.021$

For the second beating signal

${\phi}_{1,2}\left(27,500Hz\right)-{\phi}_{1,2}\left(27,000Hz\right)=0.701$

${\phi}_{1,3}\left(27,500Hz\right)-{\phi}_{1,3}\left(27,000Hz\right)=0.713$

${\phi}_{2,3}\left(27,500Hz\right)-{\phi}_{2,3}\left(27,000Hz\right)=0.012$

time delays for the second beating signal is

$\Delta {t}_{1,2}\left(0.5KHz\right)=0.000223$

$\Delta {t}_{1,3}\left(0.5KHz\right)=0.000227$

$\Delta {t}_{2,3}\left(0.5KHz\right)=0.000004$

Then the delay distance can be obtained by multiplying sound velocity.

$\Delta {r}_{1,2}=0.07582\text{m}$

$\Delta {r}_{1,3}=0.07718\text{m}$

$\Delta {r}_{2,3}=0.00136\text{m}$

Figure 11 shows the hypothetical simulation result. We were able to verify our proposal. This is one of the best solutions to find exact source location for the two acoustic beating signals or more.

${\phi}_{1,3}\left(25KHz\right)=0.0,$

${\phi}_{2,3}\left(25KHz\right)=-0.594,$

${\phi}_{1,2}\left(25.5KHz\right)=0.0,$

${\phi}_{1,3}\left(25.5KHz\right)=0.0,$

${\phi}_{2,3}\left(25.5KHz\right)=1.615$

and then we can extract phase difference for the two beating signals by using Eq.(14)

${\phi}_{1,2}\left(25.5KHz\right)-{\phi}_{1,2}\left(25KHz\right)=-0.594$

${\phi}_{1,3}\left(25.5KHz\right)-{\phi}_{1,3}\left(25KHz\right)=1.615$

${\phi}_{2,3}\left(25.5Hz\right)-{\phi}_{2,3}\left(25KHz\right)=2.209$

Time delays for the first beating signal is as below

$\Delta {t}_{1,2}\left(0.5KHz\right)=-0.000189$

$\Delta {t}_{1,3}\left(0.5KHz\right)=0.000514$

$\Delta {t}_{2,3}\left(0.5KHz\right)=0.000703$

Then the delay distance can be obtained by multiplying sound velocity.

$\Delta {r}_{1,2}=-0.06426\text{m}$

$\Delta {r}_{1,3}=0.17476\text{m}$

$\Delta {r}_{2,3}=0.23902\text{m}$

The second beating signals are as below(E2)

${\phi}_{1,2}\left(27KHz\right)=0.0$

${\phi}_{1,3}\left(27KHz\right)=0.0$

${\phi}_{2,3}\left(27KHz\right)=0.0$

${\phi}_{1,2}\left(27.5KHz\right)=0.701$

${\phi}_{1,3}\left(27.5KHz\right)=0.713$

${\phi}_{2,3}\left(27.5KHz\right)=0.021$

For the second beating signal

${\phi}_{1,2}\left(27,500Hz\right)-{\phi}_{1,2}\left(27,000Hz\right)=0.701$

${\phi}_{1,3}\left(27,500Hz\right)-{\phi}_{1,3}\left(27,000Hz\right)=0.713$

${\phi}_{2,3}\left(27,500Hz\right)-{\phi}_{2,3}\left(27,000Hz\right)=0.012$

time delays for the second beating signal is

$\Delta {t}_{1,2}\left(0.5KHz\right)=0.000223$

$\Delta {t}_{1,3}\left(0.5KHz\right)=0.000227$

$\Delta {t}_{2,3}\left(0.5KHz\right)=0.000004$

Then the delay distance can be obtained by multiplying sound velocity.

$\Delta {r}_{1,2}=0.07582\text{m}$

$\Delta {r}_{1,3}=0.07718\text{m}$

$\Delta {r}_{2,3}=0.00136\text{m}$

Figure 11 shows the hypothetical simulation result. We were able to verify our proposal. This is one of the best solutions to find exact source location for the two acoustic beating signals or more.

**Figure 11:**Source localization test result

Conclusion

We proposed here a localization technique for multiple acoustic sources using acoustic beating envelopes of which wave length are longer than audio frequency. Then we used a hypothetical simulation to prove it. The reason why we used high frequencies is not to make unnecessary acoustic noises. Additional purpose for using high frequency is to avoid unnecessary mechanical noise interferences. Because most mechanical noises are ranging under 20 KHz. We also newly proposed here how to use beating signal’s envelope. Wave length of a beating signal is much longer than those of the originals. For this reason we are free from sensor locations. We were able to get several advantages by using acoustic beating signal.

ReferencesTop

- Lubin BT, Longo R, Hammel T. Analysis of internals vibration monitoring systems data related to the St. Lucie 1 thermal shield failure. Progress in Nuclear Energy. 1988;21:117-126.
- Olma BJ. Source location and mass estimation in loose parts monitoring of LWRs. Progress in Nuclear Energy. 1985;15:583-594.
- Choi KS, Choi YC, Park JH, Kim HW. Source Localization Technique for Metallic Impact Source by Using Phase Delay between Different Type Sensors. KSNVE. 2008;18(11):1143-1149. DOI : 10.5050/KSNVN.2008.18.11.1143
- Lee HC. Robust Damage Diagnostic Method Using Short Time Fourier Transform and Beating. KSNVE. 2005;15(9):1108-1117. DOI : 10.5050/KSNVN.2005.15.9.1108
- Kryter RC, Shahrokhi F. Summary of studies on methods for detecting, locating, and characterizing metallic loose parts in nuclear reactor coolant system. US nuclear regulatory commission report NUREG/CR-2344. 1981.
- Choi KS, Kim JP, Won HK, Park JH, Kim IG. Source Localization Technique for Radar Pulse Emission by Using Scanning Method of Interest Area. KSAS. 2011;39(9):889-895.
- Cohen L. Time-frequency Distribution. Prentice Hall PTR, Engle-wood Cliffs. New Jersey. 1995.
- Randall RB. Frequency Analysis. 3rd Edition. 1987:19-31:58-62.
- Julius SB, Allan GP. RANDOM DATA, Analysis and Measurement Procedures. 2nd Edition. 1986;109-159.

FacebookTwitterLinkedInYouTubeRSS