HRV1Analyzer.h

/// czestotliwosc probkowania dla fft
#define SAMPLING_FFT 500

/// wlacza development mode: zdefiniowana => development mode, niezdefiniowana => standard mode
//#define DEV

/// wlacza tryb debugowania
#define DEBUG
//#define DEBUG_FFT
//#define DEBUG_SIG
//#define DEBUG_WINDOWS

#pragma once

#include <cmath>
#include <vector>
#include "fft/kiss_fft.h"
#include "fft/kiss_fftr.h"
#include "spline/interpolation.h"

#include "ECGRs.h"
#include "ECGHRV1.h"
#include "ECGChannelInfo.h"
#include "ECGSignal.h"

#ifdef DEV // 4 testing
    #include "ExampleSignal.h"
#else
	#include "ModulesInterfaces.h"
#endif

#define FREQUENCY_FFT (1000/SAMPLING_FFT)

/**
 * @class Class for parameters created in HRV1 module
 */
#ifndef DEV // 4 testing
	class HRV1Analyzer : public HRV1Module {
#else
	class HRV1Analyzer {
#endif

public:
	HRV1Analyzer();
	~HRV1Analyzer();

	void runModule(const ECGInfo &, const ECGRs &, ECGHRV1 &);

	#ifndef DEV
        void setParams(ParametersTypes &);
    #endif

private:
    //input parameters
	ECGRs rPeaksData;
	ECGHRV1* hrv1Data;

	//RR peaks signal parameters
	double *sig;
	double *sigAbsolute;
    unsigned int signalSize;
    unsigned int sizeFftIndex;
    unsigned int signalSampling;

    //methods
    void prepareSignal();
    void prepareSigAbsolute();
	void calculateParameters();

	double* doFFT(double* sigAfterSpline, int size);
	kiss_fft_cpx* copycpx(double *mat, int nframe);
	double mean(double *tab, int start, int end);
	double std(double *tab, int start, int end);
	double* cubicSpline(double* x, double* y, int nframe);
};