tomwalters@10: // Copyright 2008-2010, Thomas Walters
tomwalters@0: //
tomwalters@0: // AIM-C: A C++ implementation of the Auditory Image Model
tomwalters@0: // http://www.acousticscale.org/AIMC
tomwalters@0: //
tomwalters@0: // This program is free software: you can redistribute it and/or modify
tomwalters@0: // it under the terms of the GNU General Public License as published by
tomwalters@0: // the Free Software Foundation, either version 3 of the License, or
tomwalters@0: // (at your option) any later version.
tomwalters@0: //
tomwalters@0: // This program is distributed in the hope that it will be useful,
tomwalters@0: // but WITHOUT ANY WARRANTY; without even the implied warranty of
tomwalters@0: // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
tomwalters@0: // GNU General Public License for more details.
tomwalters@0: //
tomwalters@0: // You should have received a copy of the GNU General Public License
tomwalters@0: // along with this program.  If not, see <http://www.gnu.org/licenses/>.
tomwalters@0: 
tomwalters@8: /*! \file
tomwalters@8:  *  \brief Dick Lyon's Pole-Zero Filter Cascade - implemented in C++ by Tom
tomwalters@8:  *  Walters from the AIM-MAT module based on Dick Lyon's code.
tomwalters@8:  *
tomwalters@8:  *  \author Thomas Walters <tom@acousticscale.org>
tomwalters@8:  *  \date created 2008/02/05
tomwalters@8:  * \version \$Id: ModulePZFC.h 2 2010-02-02 12:59:50Z tcw $
tomwalters@8:  */
tomwalters@8: 
tomwalters@0: #ifndef _AIMC_MODULES_BMM_PZFC_H_
tomwalters@0: #define _AIMC_MODULES_BMM_PZFC_H_
tomwalters@0: 
tomwalters@0: #include <vector>
tomwalters@0: 
tomwalters@0: #include "Support/Module.h"
tomwalters@0: #include "Support/Parameters.h"
tomwalters@0: #include "Support/SignalBank.h"
tomwalters@0: 
tomwalters@0: namespace aimc {
tomwalters@0: using std::vector;
tomwalters@0: class ModulePZFC : public Module {
tomwalters@0:  public:
tomwalters@8:   explicit ModulePZFC(Parameters *pParam);
tomwalters@0:   virtual ~ModulePZFC();
tomwalters@0: 
tomwalters@8:   /*! \brief Process a buffer
tomwalters@8:    */
tomwalters@0:   virtual void Process(const SignalBank &input);
tomwalters@0: 
tomwalters@3:  private:
tomwalters@8:   /*! \brief Reset all internal state variables to their initial values
tomwalters@8:    */
tomwalters@8:   virtual void ResetInternal();
tomwalters@0: 
tomwalters@8:   /*! \brief Prepare the module
tomwalters@8:    *  \param input Input SignalBank
tomwalters@8:    *  \param output true on success false on failure
tomwalters@8:    */
tomwalters@0:   virtual bool InitializeInternal(const SignalBank &input);
tomwalters@0: 
tomwalters@8:   /*! \brief Set the filterbank parameters according to a fit matrix from Unoki
tomwalters@8:    *  and Lyon's fitting routine
tomwalters@8:    */
tomwalters@0:   bool SetPZBankCoeffsERBFitted();
tomwalters@0: 
tomwalters@8:   /*! \brief Sets the general filterbank coefficients
tomwalters@8:    */
tomwalters@0:   bool SetPZBankCoeffs();
tomwalters@0: 
tomwalters@8:   /*! \brief Automatic Gain Control
tomwalters@8:    */
tomwalters@0:   void AGCDampStep();
tomwalters@0: 
tomwalters@8:   /*! \brief Detector function - halfwave rectification etc. Used internally,
tomwalters@8:    *  but not applied to the output.
tomwalters@8:    */
tomwalters@0:   float DetectFun(float fIN);
tomwalters@0: 
tomwalters@8:   /*! \brief Minimum
tomwalters@8:    */
tomwalters@0:   inline float Minimum(float a, float b);
tomwalters@0: 
tomwalters@0:   int channel_count_;
tomwalters@0:   int buffer_length_;
tomwalters@0:   int agc_stage_count_;
tomwalters@0:   float sample_rate_;
tomwalters@0:   float last_input_;
tomwalters@0: 
tomwalters@0:   // Parameters
tomwalters@0:   // User-settable scalars
tomwalters@0:   float pole_damping_;
tomwalters@0:   float zero_damping_;
tomwalters@0:   float zero_factor_;
tomwalters@0:   float step_factor_;
tomwalters@0:   float bandwidth_over_cf_;
tomwalters@0:   float min_bandwidth_hz_;
tomwalters@0:   float agc_factor_;
tomwalters@0:   float cf_max_;
tomwalters@0:   float cf_min_;
tomwalters@0:   float mindamp_;
tomwalters@0:   float maxdamp_;
tomwalters@0:   bool do_agc_step_;
tomwalters@0: 
tomwalters@0:   // Internal Buffers
tomwalters@0:   // Initialised once
tomwalters@0:   vector<float> pole_dampings_;
tomwalters@0:   vector<float> agc_epsilons_;
tomwalters@0:   vector<float> agc_gains_;
tomwalters@0:   vector<float> pole_frequencies_;
tomwalters@0:   vector<float> za0_;
tomwalters@0:   vector<float> za1_;
tomwalters@0:   vector<float> za2_;
tomwalters@0:   vector<float> rmin_;
tomwalters@0:   vector<float> rmax_;
tomwalters@0:   vector<float> xmin_;
tomwalters@0:   vector<float> xmax_;
tomwalters@0: 
tomwalters@0:   // Modified by algorithm at each time step
tomwalters@0:   vector<float> detect_;
tomwalters@0:   vector<vector<float> > agc_state_;
tomwalters@0:   vector<float> state_1_;
tomwalters@0:   vector<float> state_2_;
tomwalters@0:   vector<float> previous_out_;
tomwalters@0:   vector<float> pole_damps_mod_;
tomwalters@0:   vector<float> inputs_;
tomwalters@0: };
tomwalters@0: }
tomwalters@0: 
tomwalters@0: #endif  // _AIMC_MODULES_BMM_PZFC_H_