test_iir_ga.cc File Reference

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

Include dependency graph for test_iir_ga.cc:

Go to the source code of this file.

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

Function Documentation

int main	(	int	argc,
		char **	argv
	)

Definition at line 42 of file test_iir_ga.cc.

References Nsound::FilterIIR::designFrequencyResponse(), Nsound::Generator::drawLine(), Nsound::Plotter::figure(), Nsound::Filter::getFrequencyAxis(), Nsound::Filter::getFrequencyResponse(), Nsound::Buffer::getLength(), Nsound::Buffer::getMax(), Nsound::Filter::plot(), Nsound::Plotter::plot(), Nsound::Plotter::show(), sr, Nsound::Plotter::text(), Nsound::Plotter::title(), Nsound::Plotter::xlabel(), Nsound::Plotter::xlim(), Nsound::Plotter::ylabel(), and Nsound::Plotter::ylim().

{
    Plotter pylab;

    float64 sr = 44100.0;

    FilterIIR f(sr, 15);

    Generator g(256);

    Buffer ref_response;

    // custom response
    ref_response << g.drawLine(0.25, 0.0, 0.8)
                 << g.drawLine(0.25, 0.8, 0.8)
                 << g.drawLine(0.15, 0.8, 1.4)
                 << g.drawLine(0.15, 1.4, 1.4)
                 << g.drawLine(0.20, 1.4, 0.0);

    // Some round off error will cause f_response
    while( ref_response.getLength() < 256 ) ref_response << 0.0;

    Buffer ref_freq_axis = g.drawLine(1.0, 0.0, sr/2.0);

//~    // Band Pass
//~    f_response << g.drawLine(0.2, 0.0, 0.0)
//~               << g.drawLine(0.2, 1.0, 1.0)
//~               << g.drawLine(0.6, 0.0, 0.0);

    pylab.plot(ref_freq_axis, ref_response);
    pylab.title("desired frequency response");

    Buffer evolution = f.designFrequencyResponse(ref_response, 0.001, 10000);

    cout << "f = " << endl << f << endl;

    pylab.figure();
    pylab.plot(evolution);
//~     pylab.ylim(0.0, 1.0);

    pylab.xlabel("Generation");
    pylab.ylabel("Error");

    int32 n_generations = evolution.getLength();
    float64 error = evolution[n_generations - 1];

    char buffer[1024];
    sprintf(buffer, "%5.2f RMS error in %d generations",
        error,
        n_generations);

    pylab.title(std::string(buffer));

    Buffer fr = f.getFrequencyResponse();
    Buffer faxis = f.getFrequencyAxis();

    pylab.figure();
    pylab.plot(ref_freq_axis, ref_response);
    pylab.plot(faxis, fr);
    pylab.xlabel("Frequency Hz");
    pylab.ylabel("Frequency Response");

    pylab.ylim(0.0, ref_response.getMax() * 1.10);

    sprintf(buffer, "Generation %06d", n_generations);

    pylab.title(std::string(buffer));

    pylab.xlim(0.0, faxis[fr.getLength() - 1]);

    sprintf(buffer, "Error = %.3f RMS", error);

    pylab.text(
        faxis[static_cast<uint32>(faxis.getLength() * 0.5)],
        0.25,
        std::string(buffer));

    f.plot();

    Plotter::show();


    return 0;
}