Edge Detection 1.2 2006/12/18 22:46:31 US/Central 2006/12/23 13:55:11.814 US/Central Barron D. Stone barron@rice.edu Rob Smith rob@rice.edu Neil K. Narayan narayann@rice.edu Brent E. Stephens brents@rice.edu Barron D. Stone barron@rice.edu Rob Smith rob@rice.edu Neil K. Narayan narayann@rice.edu Brent E. Stephens brents@rice.edu Richard G. Baraniuk richb@rice.edu Mark A. Davenport md@rice.edu beat beat detection detection edge edge detection gaussian music sound This module investigates the use of Matlab filters to determine when notes are played within a recording from an electric guitar. In order to determine the speed at which the user is playing it is first necessary to determine where in the recording individual notes are plucked. When a note is played on the guitar, a sharp peak occurs which is then followed by a decay as the vibrations in the string dampen. To determine where these edges occur the signal is first low pass filtered to remove some of the noise, then it is convolved with an edge detector function which produces a signal with peaks located at large edges (where notes are played).Due to the non-uniform qualities of guitar waveforms, an edge detector function with rounded edges produced the best results. A good function to use for an edge detector function is the derivative of a Gaussian. This produces a filter that can be tuned to work well with certain frequencies specific to the instrument being tested.

Derivative of a Gaussian This derivative of a Gaussian was produced in Matlab with f(x) = y(x)*exp(-y(x)^2/100000);

After the original signal is convolved with the Gaussian edge detector the result is a signal with peaks located where edges occurred in the original signal. A simple algorithm is used to create a vector of zeros with ones located at the locations of the peaks. Due to the nature of the guitar signal, it is not uncommon for the convolved signal to contain "double peaks" where two or more closely spaced peaks are created to represent what should be considered a single edge. To accommodate for this error, the algorithm is designed that once it finds a peak, it ignores any other peaks that occur within a given threshold range.

Orignal Signal and Convolved Signal The top plot represents the original signal from the guitar. The bottom plot represents product of the convolution of the orignal signal with the derivative of a Gaussian.

Location of Edges The vector of zeros and ones representing edge locations overlayed on top of the original signal.

To use the information produced by the edge detection to determine how fast the user is playing, the vector of ones and zeros representing peak locations is windowed into equal length chunks. The values within each of these chunks are summed together to determine the number of notes played within that window of time. From this, the program can determine whether the user was playing fast or slow during that period of time based on whether the number of notes played is above or below a previously decided threshold.

This is the Matlab code for the function EdgeDetect.

This is the Matlab code for the function GaussianFilter which works with EdgeDetect.