<< Chapter < Page Chapter >> Page >
In this section, we will go over the Matlab code we used to simulate our project, the various algorithms we tried, how we simulated "real-time", and how the matlab simulation dealt with real signals.

t 2 2 t 1 1 2 t 1 2 1 2 2 1 2 1 2 2
shows that the sum of two sinusoids that are out of phase is just another sinusoid with an amplitudedirectly related to the phase difference. Our goal is to adjust the phase difference by adding time delays to theincoming signals so as to maximize the amplitude of the output. The maximum occurs when the phase difference is zero,because the signals will add constructively. Once the maximum is found, the time delays used to achieve it tell us fromwhich direction the signals originated.

Since we are working with discrete-time signals, the time delays we tried were limited by the sampling frequency of theDSP boards, which is 48 kHz. By dividing the desired time delay by that sampling period and rounding to the nearestinteger, we converted our trial time delays into indices that could be used to select the correct sample out of the buffer.

So, the algorithm for delay-and-sum beamforming is straightforward, but there is room for a little bit ofcreativity in finding the amplitude of the summed sinusoids. We experimented with two methods to accomplish that task. Wewill call the first method "amplitude extraction," and the second "signal integration."

Amplitude extraction flow diagram

A flow diagram for the amplitude extraction method is shown in . We split the signal into two parts and multiply one part by t and the other by t . By low-pass filtering the results, we obtain the DC part of the signal which contains the amplitudeinformation. This algorithm has the obvious disadvantage that it is dependent on knowing the frequency of the incomingsignal so that the correct w is used in the multiplication step. In spite of that, we were originally selected it forimplementation on the DSP board because it showed extremely robust performance in the presence of loud noise. Addinggaussian white noise with a variance of 1 to a signal in the range -1 1 had no affect on the performance of the beamformer in our Matlab simulation. Unfortunently, the algorithm is tooslow to be used in real time. Evaluating two low-pass filters for every time delay combination tried was simply notpractical.

The signal integration method is much simpler computationally, which made it a better choice for our final implementation.We only had to square the beamformer output to make all the numbers positive, and sum the results over approximately onecycle of the incoming signal. The sum is similar to an integral over one period of the signal, except that thesamples aren't multiplied by the sampling period to make an "area." Our matlab simulation showed that the algorithm shouldwork, but that it is somewhat more sensitive to noise than amplitude extraction.

We were unable to try either of our Matlab simulations with real signals recorded from our microphone array because we had difficulty making stereo recordings. The computers we useddefaulted to recording from the microphone input (which is mono) instead of the line-in input, and we didn't haveadministrator access to change the settings.

Get Jobilize Job Search Mobile App in your pocket Now!

Get it on Google Play Download on the App Store Now




Source:  OpenStax, Audio localization. OpenStax CNX. Jun 06, 2011 Download for free at http://cnx.org/content/col10250/1.2
Google Play and the Google Play logo are trademarks of Google Inc.

Notification Switch

Would you like to follow the 'Audio localization' conversation and receive update notifications?

Ask