annotate dsp/onsets/DetectionFunction.h @ 14:68801ecbab6a

* add adaptive whitening and simple power metric
author cannam
date Thu, 09 Aug 2007 16:30:26 +0000
parents f2b5c4251bf3
children 2e3f5d2d62c1
rev   line source
cannam@0 1 /* -*- c-basic-offset: 4 indent-tabs-mode: nil -*- vi:set ts=8 sts=4 sw=4: */
cannam@0 2
cannam@0 3 /*
cannam@0 4 QM DSP Library
cannam@0 5
cannam@0 6 Centre for Digital Music, Queen Mary, University of London.
cannam@0 7 This file copyright 2005-2006 Christian Landone.
cannam@0 8 All rights reserved.
cannam@0 9 */
cannam@0 10
cannam@0 11 #ifndef DETECTIONFUNCTION_H
cannam@0 12 #define DETECTIONFUNCTION_H
cannam@0 13
cannam@0 14 #include "dsp/maths/MathUtilities.h"
cannam@0 15 #include "dsp/maths/MathAliases.h"
cannam@0 16 #include "dsp/phasevocoder/PhaseVocoder.h"
cannam@0 17 #include "base/Window.h"
cannam@0 18
cannam@0 19 #define DF_HFC (1)
cannam@0 20 #define DF_SPECDIFF (2)
cannam@0 21 #define DF_PHASEDEV (3)
cannam@0 22 #define DF_COMPLEXSD (4)
cannam@12 23 #define DF_BROADBAND (5)
cannam@14 24 #define DF_POWER (6)
cannam@0 25
cannam@0 26 struct DFConfig{
cannam@0 27 double stepSecs; // DF step in seconds
cannam@0 28 unsigned int stepSize; // DF step in samples
cannam@0 29 unsigned int frameLength; // DF analysis window - usually 2*step
cannam@0 30 int DFType; // type of detection function ( see defines )
cannam@12 31 double dbRise; // only used for broadband df (and required for it)
cannam@14 32 bool adaptiveWhitening; // perform adaptive whitening
cannam@14 33 double whiteningRelaxCoeff; // if < 0, a sensible default will be used
cannam@14 34 double whiteningFloor; // if < 0, a sensible default will be used
cannam@0 35 };
cannam@0 36
cannam@0 37 class DetectionFunction
cannam@0 38 {
cannam@0 39 public:
cannam@0 40 double* getSpectrumMagnitude();
cannam@0 41 DetectionFunction( DFConfig Config );
cannam@0 42 virtual ~DetectionFunction();
cannam@0 43 double process( double* TDomain );
cannam@2 44 double process( double* magnitudes, double* phases );
cannam@0 45
cannam@0 46 private:
cannam@14 47 void whiten();
cannam@2 48 double runDF();
cannam@2 49
cannam@0 50 double HFC( unsigned int length, double* src);
cannam@0 51 double specDiff( unsigned int length, double* src);
cannam@14 52 double phaseDev(unsigned int length, double *srcPhase);
cannam@0 53 double complexSD(unsigned int length, double *srcMagnitude, double *srcPhase);
cannam@14 54 double broadband(unsigned int length, double *srcMagnitude);
cannam@14 55 double power(unsigned int length, double *src);
cannam@0 56
cannam@0 57 private:
cannam@0 58 void initialise( DFConfig Config );
cannam@0 59 void deInitialise();
cannam@0 60
cannam@0 61 int m_DFType;
cannam@0 62 unsigned int m_dataLength;
cannam@0 63 unsigned int m_halfLength;
cannam@13 64 double m_stepSecs;
cannam@13 65 unsigned int m_stepSize;
cannam@12 66 double m_dbRise;
cannam@14 67 bool m_whiten;
cannam@14 68 double m_whitenRelaxCoeff;
cannam@14 69 double m_whitenFloor;
cannam@0 70
cannam@2 71 double* m_magHistory;
cannam@2 72 double* m_phaseHistory;
cannam@2 73 double* m_phaseHistoryOld;
cannam@14 74 double* m_magPeaks;
cannam@0 75
cannam@0 76 double* m_DFWindowedFrame; // Array for windowed analysis frame
cannam@0 77 double* m_magnitude; // Magnitude of analysis frame ( frequency domain )
cannam@0 78 double* m_thetaAngle;// Phase of analysis frame ( frequency domain )
cannam@0 79
cannam@0 80 Window<double> *m_window;
cannam@2 81 PhaseVocoder* m_phaseVoc; // Phase Vocoder
cannam@0 82 };
cannam@0 83
cannam@0 84 #endif