iirfilter.h

Go to the documentation of this file.

//=============================================================================================================
 
#ifndef IIRFILTER_DSP_H
#define IIRFILTER_DSP_H
 
//=============================================================================================================
// INCLUDES
//=============================================================================================================
 
#include "dsp_global.h"
 
//=============================================================================================================
// EIGEN INCLUDES
//=============================================================================================================
 
#include <Eigen/Core>
 
//=============================================================================================================
// QT INCLUDES
//=============================================================================================================
 
#include <QVector>
 
//=============================================================================================================
// DEFINE NAMESPACE UTILSLIB
//=============================================================================================================
 
namespace UTILSLIB
{
 
//=============================================================================================================
struct DSPSHARED_EXPORT IirBiquad
{
    double b0 = 1.0; 
    double b1 = 0.0; 
    double b2 = 0.0; 
    double a1 = 0.0; 
    double a2 = 0.0; 
};
 
//=============================================================================================================
class DSPSHARED_EXPORT IirFilter
{
public:
    //=========================================================================================================
    enum FilterType
    {
        LowPass = 0,  
        HighPass = 1, 
        BandPass = 2, 
        BandStop = 3  
    };
 
    //=========================================================================================================
    static QVector<IirBiquad> designButterworth(int iOrder,
                                                FilterType type,
                                                double dCutoffLow,
                                                double dCutoffHigh,
                                                double dSFreq);
 
    //=========================================================================================================
    static Eigen::RowVectorXd applySos(const Eigen::RowVectorXd& vecData,
                                       const QVector<IirBiquad>& sos);
 
    //=========================================================================================================
    static Eigen::RowVectorXd applyZeroPhase(const Eigen::RowVectorXd& vecData,
                                             const QVector<IirBiquad>& sos);
 
    //=========================================================================================================
    static Eigen::MatrixXd applyZeroPhaseMatrix(const Eigen::MatrixXd& matData,
                                                const QVector<IirBiquad>& sos);
 
private:
    //=========================================================================================================
    static QVector<std::complex<double>> butterworthPrototypePoles(int n);
 
    //=========================================================================================================
    static IirBiquad poleToDigitalBiquad(std::complex<double> pole,
                                         double dC,
                                         double dGain);
 
    //=========================================================================================================
    static IirBiquad realPoleToDigitalSection(double dPoleReal,
                                              double dC,
                                              double dGain);
};
 
} // namespace UTILSLIB
 
#endif // IIRFILTER_DSP_H