annotate dsp/onsets/DetectionFunction.h @ 16:2e3f5d2d62c1

* Move dsp/maths to maths ; bring PCA and HMM across from Soundbite
author cannam
date Wed, 09 Jan 2008 10:31:29 +0000
parents 68801ecbab6a
children 5bec06ecc88a
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@16 14 #include "maths/MathUtilities.h"
cannam@16 15 #include "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@0 24
cannam@0 25 struct DFConfig{
cannam@0 26 double stepSecs; // DF step in seconds
cannam@0 27 unsigned int stepSize; // DF step in samples
cannam@0 28 unsigned int frameLength; // DF analysis window - usually 2*step
cannam@0 29 int DFType; // type of detection function ( see defines )
cannam@12 30 double dbRise; // only used for broadband df (and required for it)
cannam@14 31 bool adaptiveWhitening; // perform adaptive whitening
cannam@14 32 double whiteningRelaxCoeff; // if < 0, a sensible default will be used
cannam@14 33 double whiteningFloor; // if < 0, a sensible default will be used
cannam@0 34 };
cannam@0 35
cannam@0 36 class DetectionFunction
cannam@0 37 {
cannam@0 38 public:
cannam@0 39 double* getSpectrumMagnitude();
cannam@0 40 DetectionFunction( DFConfig Config );
cannam@0 41 virtual ~DetectionFunction();
cannam@0 42 double process( double* TDomain );
cannam@2 43 double process( double* magnitudes, double* phases );
cannam@0 44
cannam@0 45 private:
cannam@14 46 void whiten();
cannam@2 47 double runDF();
cannam@2 48
cannam@0 49 double HFC( unsigned int length, double* src);
cannam@0 50 double specDiff( unsigned int length, double* src);
cannam@14 51 double phaseDev(unsigned int length, double *srcPhase);
cannam@0 52 double complexSD(unsigned int length, double *srcMagnitude, double *srcPhase);
cannam@14 53 double broadband(unsigned int length, double *srcMagnitude);
cannam@0 54
cannam@0 55 private:
cannam@0 56 void initialise( DFConfig Config );
cannam@0 57 void deInitialise();
cannam@0 58
cannam@0 59 int m_DFType;
cannam@0 60 unsigned int m_dataLength;
cannam@0 61 unsigned int m_halfLength;
cannam@13 62 double m_stepSecs;
cannam@13 63 unsigned int m_stepSize;
cannam@12 64 double m_dbRise;
cannam@14 65 bool m_whiten;
cannam@14 66 double m_whitenRelaxCoeff;
cannam@14 67 double m_whitenFloor;
cannam@0 68
cannam@2 69 double* m_magHistory;
cannam@2 70 double* m_phaseHistory;
cannam@2 71 double* m_phaseHistoryOld;
cannam@14 72 double* m_magPeaks;
cannam@0 73
cannam@0 74 double* m_DFWindowedFrame; // Array for windowed analysis frame
cannam@0 75 double* m_magnitude; // Magnitude of analysis frame ( frequency domain )
cannam@0 76 double* m_thetaAngle;// Phase of analysis frame ( frequency domain )
cannam@0 77
cannam@0 78 Window<double> *m_window;
cannam@2 79 PhaseVocoder* m_phaseVoc; // Phase Vocoder
cannam@0 80 };
cannam@0 81
cannam@0 82 #endif