csd.h

Go to the documentation of this file.

//=============================================================================================================
 
#ifndef CSD_H
#define CSD_H
 
//=============================================================================================================
// INCLUDES
//=============================================================================================================
 
#include "dsp_global.h"
 
//=============================================================================================================
// EIGEN INCLUDES
//=============================================================================================================
 
#include <Eigen/Core>
#include <Eigen/Dense>
 
//=============================================================================================================
// QT INCLUDES
//=============================================================================================================
 
#include <QVector>
 
//=============================================================================================================
// DEFINE NAMESPACE UTILSLIB
//=============================================================================================================
 
namespace UTILSLIB
{
 
//=============================================================================================================
struct DSPSHARED_EXPORT CsdResult
{
    Eigen::MatrixXcd matCsd;                
    QVector<Eigen::MatrixXcd> csdByFreq;    
    Eigen::RowVectorXd vecFreqs;            
};
 
//=============================================================================================================
class DSPSHARED_EXPORT Csd
{
public:
    //=========================================================================================================
    static CsdResult computeMultitaper(const Eigen::MatrixXd& matData,
                                       double sfreq,
                                       double fmin = 0.0,
                                       double fmax = -1.0,
                                       double halfBandwidth = 4.0,
                                       int nTapers = -1);
 
    //=========================================================================================================
    static CsdResult computeFourier(const Eigen::MatrixXd& matData,
                                    double sfreq,
                                    double fmin = 0.0,
                                    double fmax = -1.0,
                                    int nFft = 256,
                                    double overlap = 0.5);
 
    //=========================================================================================================
    static CsdResult computeMorlet(const Eigen::MatrixXd& matData,
                                   double sfreq,
                                   const Eigen::RowVectorXd& frequencies,
                                   int nCycles = 7);
};
 
} // namespace UTILSLIB
 
#endif // CSD_H