Click here to Skip to main content
Click here to Skip to main content

FFT Guitar Tuner

, 10 Aug 2010 MIT
Rate this:
Please Sign up or sign in to vote.
Using a Fast Fourier Transform to calculate the fundamental frequency of the captured audio sound

Introduction

This article shows how to use a Fast Fourier Transform (FFT) algorithm to calculate the fundamental frequency of a captured audio sound. Also, we will see how to apply the algorithm to analyze live sound to build a simple guitar tuner: the code provides a solution to the problem of calculation of the fundamental frequency of the played pitch.

Background

The computer can capture live sound/music using a microphone that is connected to the sound card. Modern sound cards can capture digital signals. A digital signal is a set of quantized sound values that were taken in uniformly spaced times. The digital signal does not provide any information about frequencies that are present in the sound. To determine that, the data need to be analyzed.

The Short-Time Fourier Transform (STFT) makes representation of the phase and magnitude of the signal. The result of the STFT can be used to produce the spectrogram of the signal: the magnitude squared over time and frequencies. We will use a Fast Fourier Transform (FFT) to generate the spectrogram of the signal of short periods of time. After the spectrogram is calculated, the fundamental frequency can be determined by finding the index of the maximum value of the magnitude squared. The improved algorithm finds several such places, candidate frequency bins, with the magnitude squared in the top of the maximum values, and further analyzes them to verify the candidate fundamental frequencies by using the signal data.

When a note is played on a musical instrument, the sound waves are generated by strings, air, or the speaker - an instrument generates a musical note. One of the characteristics of a musical note is a pitch (fundamental frequency). Traditionally musical alphabet frequencies are divided by octaves, and then by semitones. An octave has 12 named pitches: C (prime), C#, D, D#, E, F, F#, G, G#, A, A#, and B. Octaves also have names: great, small, one-lined, two-lined, etc. The "standard pitch" (A one-lined or A4) has a fundamental frequency of its sound waves equals to 440 Hz. The frequencies of two neighboring notes are different by 21/12, and frequencies of the notes with the same name in two neighboring octaves are different by 2.

Table: Notes and Their Fundamental Frequencies
Note Name Traditional Octave Names (Scientific), Hz
Great (2) Small (3) One-lined (4) Two-lined (5)
C 65.4064 130.8128 261.6256 523.2511
C# 69.2957 138.5913 277.1826 554.3653
D 73.4162 146.8324 293.6648 587.3295
D# 77.7817 155.5635 311.1270 622.2540
E 82.4069 164.8138 329.6276 659.2551
F 87.3071 174.6141 349.2282 698.4565
F# 92.4986 184.9972 369.9944 739.9888
G 97.9989 195.9977 391.9954 783.9909
G# 103.8262 207.6523 415.3047 830.6094
A 110.0000 220.0000 440.0000 880.0000
A# 116.5409 233.0819 466.1638 932.3275
B 123.4708 246.9417 493.8833 987.7666

The typical (six string) guitar normally plays pitches of great through two-lined octaves. The pitches of the open strings (E2, A2, D3, G3, B3, and E4) are selected in the table in bold.

Using the Code

The solution contains three projects: the main windows application (FftGuitarTuner), the sound analysis library (SoundAnalysis), and the sound capture library (SoundCapture). The heart of the solution and the SoundAnalysis project is the FFT algorithm (see the Calculate method of the SoundAnalysis.FftAlgorithm class):

// bit reversal
ComplexNumber[] data = new ComplexNumber[length];
for (int i = 0; i < x.Length; i++)
{
    int j = ReverseBits(i, bitsInLength);
    data[j] = new ComplexNumber(x[i]);
}

// Cooley-Tukey 
for (int i = 0; i < bitsInLength; i++)
{
    int m = 1 << i;
    int n = m * 2;
    double alpha = -(2 * Math.PI / n);

    for (int k = 0; k < m; k++)
    {
        // e^(-2*pi/N*k)
        ComplexNumber oddPartMultiplier = 
           new ComplexNumber(0, alpha * k).PoweredE();

        for (int j = k; j < length; j += n)
        {
            ComplexNumber evenPart = data[j];
            ComplexNumber oddPart = oddPartMultiplier * data[j + m];
            data[j] = evenPart + oddPart;
            data[j + m] = evenPart - oddPart;
        }
    }
}

// calculate spectrogram
double[] spectrogram = new double[length];
for (int i = 0; i < spectrogram.Length; i++)
{
    spectrogram[i] = data[i].AbsPower2();
}

The data for the algorithm is provided from the sound card capture buffer. The abstract SoundCapture.SoundCaptureBase utility class is an adapter for DirectSound's Capture and CaptureBuffer classes, that helps to encapsulate buffering and setting up the audio format parameters. The application requires Microsoft DirectX 9 runtime components for the live sound capture from the microphone.

Main Application Form

Figure: Main Application Form

After the application is started, select the sound device and play a note. The application will capture the live sound and will calculate the current fundamental frequency of the signal. The information can be used to tune the guitar.

Points of Interest

To calculate the Fast Fourier Transform, the Cooley-Tukey algorithm was used. It gives good performance for the required task. To challenge the algorithm, the application analyses about 22,000 sample blocks in real time: the sound is captured at a 44,100 Hz rate and a 16 bits sample size, and the analysis is performed twice a second.

The sound analysis library can be used for tone, background noise, sound, or speech detection. Series of the spectrogram of the continued sound can be displayed as a 2D (or 3D) image to present it visually.

References

  1. "Musical Note", Wikipedia
  2. "Short-Time Fourier Transform", Wikipedia
  3. "Fast Fourier Transform", Wikipedia
  4. "Cooley-Tukey FFT Algorithm", Wikipedia

History

  • 1st January, 2009: Initial version
  • 2nd January, 2009: Added algorithm code snippet
  • 3rd August, 2010: Corrected article typos; new frequency detection algorithm

License

This article, along with any associated source code and files, is licensed under The MIT License

Share

About the Author

notmasteryet
Software Developer
United States United States
No Biography provided
Follow on   Twitter

Comments and Discussions

 
QuestionFundamental frequency algoritm PinmemberKvanttt1-Jul-11 9:46 
BugI'm having problems using the VB .Net version, running on Windows 7 64-bit PinmemberTheJediMaster30-Jun-11 15:24 
GeneralRe: I'm having problems using the VB .Net version, running on Windows 7 64-bit Pinmembernotmasteryet30-Jun-11 15:26 
GeneralRe: I'm having problems using the VB .Net version, running on Windows 7 64-bit PinmemberTheJediMaster5-Jul-11 12:23 
GeneralRe: I'm having problems using the VB .Net version, running on Windows 7 64-bit PinmemberMember 826810726-Sep-11 11:41 
GeneralRe: I'm having problems using the VB .Net version, running on Windows 7 64-bit PinmemberMember 766644016-Dec-11 7:26 
GeneralRe: I'm having problems using the VB .Net version, running on Windows 7 64-bit Pinmemberbmac20-Jul-13 7:05 
QuestionIs it easy to convert to VB. Net? PinmemberTheJediMaster30-Jun-11 6:23 
AnswerRe: Is it easy to convert to VB. Net? Pinmembernotmasteryet30-Jun-11 15:25 
GeneralRe: Is it easy to convert to VB. Net? PinmemberTheJediMaster30-Jun-11 15:26 
GeneralReverseBits function optimization PinmemberKvanttt15-May-11 4:37 
QuestionSystem sound PinmemberPareen Vatani10-May-11 21:53 
GeneralMy vote of 5 PinmemberP1l19r1m8-May-11 1:11 
QuestionNice solution thanks for providing it and a ? PinmemberAlasdair Macleod22-Jan-11 12:24 
Hey
 
Will add my thanks as your code is proving very useful with my dissertation project.
 
I have been trying to add a routine that accesses the volume of the capturebuffer but although this is available in other buffers it is not readily available as a property in a capturebuffer. However I have also looked at another project which is a staright forward volume meter which manages to do so. I am new to C# but chose this route as your library is the simplest most complete for the frequency to musical note conversion, but find I need to add a noise gate to achieve my goals. Pretty sure this can be done if The sample amplitude can be calculated to create a threshhold feature.
 
Thanks
 
Alasdair
AnswerRe: Nice solution thanks for providing it and a ? PinmemberAlasdair Macleod25-Jan-11 1:14 
GeneralA Bug?. Given a 44100HZ wav, the FFT result repeat and shrink many times in above 11025HZ. Pinmemberjasonhh14-Aug-10 0:48 
GeneralRe: A Bug?. Given a 44100HZ wav, the FFT result repeat and shrink many times in above 11025HZ. Pinmembernotmasteryet14-Aug-10 3:51 
GeneralRe: A Bug?. Given a 44100HZ wav, the FFT result repeat and shrink many times in above 11025HZ. Pinmemberjasonhh17-Aug-10 1:16 
GeneralRe: A Bug?. Given a 44100HZ wav, the FFT result repeat and shrink many times in above 11025HZ. Pinmembernotmasteryet17-Aug-10 15:00 
GeneralRe: A Bug?. Given a 44100HZ wav, the FFT result repeat and shrink many times in above 11025HZ. Pinmemberjasonhh17-Aug-10 20:34 
GeneralMy vote of 5 PinmemberJoel Ivory Johnson10-Aug-10 18:04 
QuestionHow do you compare fundamental frequencies? Pinmembermarluv23-Jul-10 5:04 
AnswerRe: How do you compare fundamental frequencies? Pinmembernotmasteryet8-Aug-10 15:23 
GeneralOctave Problem PinmemberJohn Parse13-Jul-10 21:55 
GeneralRe: Octave Problem Pinmembernotmasteryet8-Aug-10 15:37 
GeneralCooley-Tukey PinmemberPKultimPK28-Jun-10 5:32 
GeneralRe: Cooley-Tukey Pinmembernotmasteryet8-Aug-10 17:26 
GeneralC# to VB Pinmemberpeter879-Jun-10 23:12 
GeneralRe: C# to VB Pinmembernotmasteryet9-Aug-10 14:14 
GeneralViolin PinmemberZee-J20-Apr-10 19:37 
GeneralRe: Violin Pinmembernotmasteryet21-Apr-10 18:23 
GeneralRe: Violin PinmemberZee-J21-Apr-10 20:51 
GeneralRe: Violin Pinmember5thWheel7-Sep-10 23:45 
GeneralVb.net PinmemberKrio889-Apr-10 4:15 
GeneralRe: Vb.net Pinmemberpeter879-Jun-10 23:16 
QuestionHow do I FFT only part of spectrum? PinmemberBengt Thulin23-Mar-10 7:49 
QuestionHow to determine how long a note is held for? Pinmemberme8meme9me13-Mar-10 5:43 
Questionhow about in VB.net? Pinmemberpeter8712-Mar-10 0:57 
AnswerRe: how about in VB.net? Pinmemberme8meme9me13-Mar-10 5:59 
Generallanguage Pinmemberbintangtan7-Mar-10 19:01 
QuestionUSB Input ? Pinmemberwijays@gmail.com4-Feb-10 10:52 
GeneralProblems running FFT Guitar Tuner on Windows 7... Hmmmmm PinmemberFlemming Jensen217-Jan-10 2:30 
Questionmsg problem: is not a win32 valid application Pinmemberangelo carlt6-Jan-10 7:35 
AnswerRe: msg problem: is not a win32 valid application Pinmembernotmasteryet6-Jan-10 17:25 
QuestionRe: msg problem: is not a win32 valid application Pinmemberangelo carlt7-Jan-10 13:03 
AnswerRe: msg problem: is not a win32 valid application Pinmemberdaylightdj17-Mar-10 2:26 
GeneralRe: msg problem: is not a win32 valid application Pinmemberangelo carlt17-Mar-10 18:20 
GeneralRe: msg problem: is not a win32 valid application Pinmemberdaylightdj24-Mar-10 2:35 
Generalpls help me. . . Pinmemberclareese6-Jan-10 2:34 
GeneralRe: pls help me. . . Pinmembernotmasteryet6-Jan-10 17:20 

General General    News News    Suggestion Suggestion    Question Question    Bug Bug    Answer Answer    Joke Joke    Rant Rant    Admin Admin   

Use Ctrl+Left/Right to switch messages, Ctrl+Up/Down to switch threads, Ctrl+Shift+Left/Right to switch pages.

| Advertise | Privacy | Terms of Use | Mobile
Web01 | 2.8.141220.1 | Last Updated 10 Aug 2010
Article Copyright 2009 by notmasteryet
Everything else Copyright © CodeProject, 1999-2014
Layout: fixed | fluid