static void MaxAcf (mrs_real& max, mrs_real& mean, const realvec& beatHistogram, realvec& res,mrs_natural startSearchAt, mrs_natural stopSearchAt)
{
  mrs_natural k,len = beatHistogram.getCols ();

  res.setval(0.);

  // compute ACF
  for (k = startSearchAt; k < stopSearchAt; k++) // this can be optimized
  {
    mrs_real val	= 0;

    for (mrs_natural i = k; i < len; i++)
    {
      val += beatHistogram(i) * beatHistogram(i-k);
    }


    res(k)	= val / (len-k);
  }


  //pkr_->process(in, pkres_);

  max = res.maxval ();

  mean = 1e6*res.mean ();
}
Ejemplo n.º 2
0
void
Peak2Rms::myProcess(realvec& in, realvec& out)
{
  mrs_real peak = in.maxval();
  mrs_real rms = 0.0;
  mrs_natural t;

  for (t = 0; t < inSamples_; t++)
  {
    rms += in(0,t) * in(0,t);
  }
  rms /= inSamples_;
  rms = sqrt(rms);
  out(0) = peak/rms;
}