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.
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

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 . We can assume that each wave is composed of two frequencies of These two frequencies might satisfy Eq.(1) Then we can express as below.
Then we can get slightly different two close waves . 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}\left(\text{t}\right)$is another real-valued time domain signal denoted by $\stackrel{̂}{x}\left(t\right)$, such that $z\left(t\right)=x\left(t\right)+j\stackrel{̂}{x}\left(t\right)$ . The Hilbert transform of $\text{x}\left(\text{t}\right)$extending over the range $-\infty is a real-valued function $\stackrel{̂}{x}\left(t\right)$defined by $\stackrel{̂}{x}\left(t\right)$is the convolution integral of $x\left(t\right)$and $1}{\pi t}$ , written as A useful way to compute the Hilbert Transform $\stackrel{̂}{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}\stackrel{̂}{x}\left(t\right)$ . One can also write $\text{z}\left(\text{t}\right)=A\left(t\right){e}^{j\theta \left(t\right)}$,where $\text{A}\left(\text{t}\right)$ denotes the envelope signal of $x\left(t\right)$ and $\text{θ}\left(\text{t}\right)$ denotes the instantaneous phase signal of $x\left(t\right)$ . In terms of $x\left(t\right)$and $\stackrel{̂}{x}\left(t\right)$ , it is clear that 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}\left(={\omega }_{h}+{\omega }_{l}\right),$ ${\omega }_{1}\left(={\omega }_{h}+{\omega }_{l}\right),$and its correspoding phases $\text{φ}\left({\omega }_{1}\right)$and its correspoding phases $\text{φ}\left({\omega }_{1}\right)$ $\text{φ}\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{φ}\left({\omega }_{beat}\right)$can be expressed as below We can express a beating sginal from Eq.(2) as below 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
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\left(t\right)$and $q\left(t\right)$.We could think $p\left(t\right)$and $q\left(t\right)$are received signals at a different two sensors, and ${\tau }_{d}$means the time delay between $p\left(t\right)$and $q\left(t\right)$.Then we can express it as
In this case, the frequency and phase of beating signal $p\left(t\right)$are ${\text{ω}}_{1}-{\text{ω}}_{2}\left({\text{ω}}_{1}>{\text{ω}}_{2}\right),$ respectively. Then, where 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. 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\left(\omega \right)$and $\text{Q}\left(\omega \right)$denote the components of the Fourier transform of $p\left(t\right)$ and $q\left(t\right)$;the cross-spectrum is where ${\text{φ}}_{p}\left(\omega \right)$and ${\text{φ}}_{q}\left(\omega \right)$ 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 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}=1}{f}$The time delay of acoustic beating envelope $\Delta {t}_{p,q}\left(f\right)$can be obtained by multiplying the period by the ratio of phase delay. Therefore, the time delay can be expressed as following equation. The distance difference transferred from the impact source to the two sensors is derived by multiplying Eq.(11) by the sound velocity.
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.
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
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
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
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.
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

Listing : ICMJE