FundamentalFromPartials createEstimator( void )
{
FundamentalFromPartials eparts;
eparts.setAmpFloor( Floor );
eparts.setAmpRange( Range );
eparts.setFreqCeiling( Ceiling );
eparts.setPrecision( Precision );
return eparts;
}
int main( int argc, char * argv[] )
{
std::cout << "Unit test for fundamental estimation functions." << endl;
std::cout << "Tests FundamentalFromPartials and FundamentalFromSamples." << endl << endl;
std::cout << "Relies on AiffFile, Analyzer, Partial, PartialList, and LinearEnvelope." << endl << endl;
std::cout << "Built: " << __DATE__ << endl << endl;
string path("");
if ( std::getenv("srcdir") )
{
path = std::getenv("srcdir");
path = path + "/";
}
// --- parameters that are sample-specific ---
string fname = path + "clarinet.aiff";
double fmin = 200;
double fmax = 500;
double res = 415*.8;
double win = 415*1.6;
double interval = 0.05; // 50 ms
const double approx = 414;
double x;
try
{
// import (mono) samples
std::vector< double > buf;
double rate;
importSamples( fname, buf, &rate );
// step 1. analyze the samples
Analyzer anal( res, win );
anal.buildFundamentalEnv( fmin, fmax );
anal.setHopTime( interval );
cout << "--- step 1 analyzer ---" << endl;
cout << "analysis resolution is " << anal.freqResolution() << endl;
cout << "window width is " << anal.windowWidth() << endl;
cout << "amplitude threshold is " << anal.ampFloor() << endl;
cout << "lower bound is " << fmin << endl;
cout << "upper bound is " << fmax << endl;
PartialList plist = anal.analyze( buf, rate );
LinearEnvelope est1 = anal.fundamentalEnv();
x = dumpEnvelope( est1 );
if ( (approx-1) > x || (approx+1) < x )
{
throw std::runtime_error( "that isn't right" );
}
double tbeg = est1.begin()->first;
double tend = (--est1.end())->first;
// step 2. estimate fundamental from Partials
FundamentalFromPartials eparts;
eparts.setAmpFloor( -65 );
eparts.setAmpRange( 40 );
eparts.setFreqCeiling( 5000 );
cout << "--- step 2 fundamental estimator from partials ---" << endl;
cout << "amplitude threshold is " << eparts.ampFloor() << endl;
cout << "amplitude range is " << eparts.ampRange() << endl;
cout << "frequency ceiling is " << eparts.freqCeiling() << endl;
cout << "precision level is " << eparts.precision() << " Hz" << endl;
LinearEnvelope est2 = eparts.buildEnvelope( plist, tbeg, tend,
interval, fmin, fmax, 0.95 );
x = dumpEnvelope( est2 );
if ( (approx-1) > x || (approx+1) < x )
{
throw std::runtime_error( "that isn't right" );
}
// step 3. estimate fundamental from the samples
FundamentalFromSamples esamps( win );
esamps.setAmpFloor( -65 );
esamps.setAmpRange( 40 );
esamps.setFreqCeiling( 5000 );
cout << "--- step 3 fundamental estimator from samples ---" << endl;
cout << "window width is " << esamps.windowWidth() << endl;
cout << "amplitude threshold is " << esamps.ampFloor() << endl;
cout << "amplitude range is " << esamps.ampRange() << endl;
cout << "frequency ceiling is " << esamps.freqCeiling() << endl;
cout << "precision level is " << esamps.precision() << " Hz" << endl;
LinearEnvelope est3 = esamps.buildEnvelope( buf, rate, tbeg, tend,
interval, fmin, fmax, 0.95 );
x = dumpEnvelope( est3 );
if ( (approx-1) > x || (approx+1) < x )
{
throw std::runtime_error( "that isn't right" );
}
}
catch( Exception & ex )
{
cout << "Caught Loris exception: " << ex.what() << endl;
return 1;
}
catch( std::exception & ex )
{
cout << "Caught std C++ exception: " << ex.what() << endl;
return 1;
}
// return successfully
cout << "Fundamental estimation passed all tests." << endl;
return 0;
return 0;
}