test_fft.cc File Reference

#include <Nsound/NsoundAll.h>
#include "Test.h"
#include <stdlib.h>
#include <cmath>
#include <iostream>
#include <sstream>

Include dependency graph for test_fft.cc:

Go to the source code of this file.

Functions
int	main (int argc, char **argv)

Variables
static const char *	THIS_FILE = "test_fft.cc"

Function Documentation

int main	(	int	argc,
		char **	argv
	)

Definition at line 47 of file test_fft.cc.

References Nsound::Buffer::abs(), Nsound::FFTransform::fft(), Nsound::Generator::generate(), Nsound::Buffer::getMax(), Nsound::FFTransform::ifft(), Nsound::FilterHighPassIIR::plot(), Nsound::FilterLowPassIIR::plot(), Nsound::FilterLowPassFIR::plot(), Nsound::FilterHighPassFIR::plot(), Nsound::FilterBandPassFIR::plot(), Nsound::FilterBandRejectFIR::plot(), Nsound::FilterBandPassIIR::plot(), Nsound::FilterBandRejectIIR::plot(), SAMPLE_RATE, Nsound::Plotter::show(), TEST_ERROR_HEADER, TEST_HEADER, Nsound::Tic(), Nsound::Toc(), Nsound::Plotter::xlim(), and Nsound::Plotter::ylim().

{
    uint32 SAMPLE_RATE = 44100;

    Sine sine(SAMPLE_RATE);

    Buffer input = sine.generate(2.0,  4000.0)
                 + sine.generate(2.0,  8000.0)
                 + sine.generate(2.0, 12000.0);

    FFTransform fft(SAMPLE_RATE);

    Tic();
    cout << TEST_HEADER
         << "FFTransform::fft(), FFTransform::ifft() ... " << flush;

    FFTChunkVector vec = fft.fft(input, 2048);

    // Call the inverse fft.
    Buffer ifft = fft.ifft(vec);

    Buffer diff = input - ifft;

    diff.abs();

    if(diff.getMax() > 0.00005)
    {
        cerr << TEST_ERROR_HEADER
             << "fft() - ifft() > 0.00005, max(diff) = "
             << diff.getMax()
             << endl;
        exit(1);
    }

    cout << Toc() << " seconds: SUCCESS" << endl;

    float64 fl = 440;
    float64 fh = fl * 2;

    FilterLowPassFIR    f1(SAMPLE_RATE, 512, fl);
    FilterHighPassFIR   f2(SAMPLE_RATE, 512, fl);
    FilterBandPassFIR   f3(SAMPLE_RATE, 512, fl, fh);
    FilterBandRejectFIR f4(SAMPLE_RATE, 512, fl, fh);

    FilterLowPassIIR    f11(SAMPLE_RATE, 6, fl, 0.01);
    FilterHighPassIIR   f22(SAMPLE_RATE, 6, fl, 0.01);
    FilterBandPassIIR   f33(SAMPLE_RATE, 6, fl, fh, 0.01);
    FilterBandRejectIIR f44(SAMPLE_RATE, 6, fl, fh, 0.01);

    Tic();
    cout << TEST_HEADER
         << "Plotting Frequecy Responses of Filters ... " << flush;

    Plotter pylab;

    f4.plot();
    pylab.xlim(0, fh * 2);
    pylab.ylim(-40, 5);

    f44.plot();
    pylab.xlim(0, fh * 2);
    pylab.ylim(-40, 5);

    f3.plot();
    pylab.xlim(0,fh * 2);
    pylab.ylim(-40,5);

    f33.plot();
    pylab.xlim(0,fh * 2);
    pylab.ylim(-40,5);

    f2.plot();
    pylab.xlim(0,fh * 2);
    pylab.ylim(-40,5);

    f22.plot();
    pylab.xlim(0,fh * 2);
    pylab.ylim(-40,5);

    f1.plot();
    pylab.xlim(0,fh * 2);
    pylab.ylim(-40,5);

    f11.plot();
    pylab.xlim(0,fh * 2);
    pylab.ylim(-40,5);

    cout << Toc() << " seconds" << endl << flush;

    Plotter::show();

    return 0;
}

Variable Documentation

const char* THIS_FILE = "test_fft.cc"

static

Definition at line 45 of file test_fft.cc.