FreqDomainIdent.h

#ifndef _CREQDOMAIN_IDENT_H_
#define _CREQDOMAIN_IDENT_H_
 
#include <complex>
#include <vector>
 
#include <Eigen/Eigen>
 
using namespace Eigen;
 
namespace Ship{
    class PotentialData;
 
    class FreqDomainIdent
    {
    public:
        enum IDENTIFICATION_TYPE
        {
            INVERSE_FREQUENCY,
            VECTOR_FITTING
        };
    
        struct Polynomial
        {
            Matrix< double, Dynamic, Dynamic > num;
            Matrix< double, Dynamic, Dynamic > den;
        };
    
        struct StateSpace
        {
            Matrix< double, Dynamic, Dynamic > A;
            Matrix< double, Dynamic, Dynamic > B;
            Matrix< double, Dynamic, Dynamic > C;
            Matrix< double, Dynamic, Dynamic > D;   
        };
    
    
        FreqDomainIdent::IDENTIFICATION_TYPE IdentificationMethod() const { return Type; }
        void IdentificationMethod(FreqDomainIdent::IDENTIFICATION_TYPE val) { Type = val; }
    
        double GetAinf( int DOF1, int DOF2, unsigned int SpeedNo=0)
        {
            return Ainf[SpeedNo](DOF1,DOF2);
        }
    
        void SetPotentialData( Ship::PotentialData *Data );
    
        Polynomial IdentifyTF(unsigned int DOF1, unsigned int DOF2, unsigned int SpeedNo=0, bool UseHanWindow=false);
        StateSpace IdentifySS(unsigned int DOF1, unsigned int DOF2, unsigned int SpeedNo=0, bool UseHanWindow=false);
    
    private:
        static const std::complex<double> j;
    
        //static Matrix< std::complex<double>, Dynamic, Dynamic > Poly( Matrix< std::complex<double>,  Dynamic, Dynamic > &coeff );
    
        static Matrix< std::complex<double>, 1, Dynamic > PolyVal(  const Matrix< double, Dynamic, Dynamic > &coeff, const Matrix< std::complex<double>, 1, Dynamic > &val);
    
        //static Matrix< std::complex<double>, Dynamic , Dynamic > Roots( Matrix< double, Dynamic, Dynamic > &coeff);
    
        //static Matrix< std::complex<double>, 1, Dynamic> HannWindow( unsigned int N );
    
        //static Matrix<double, 1, Dynamic> Interp1D(   Matrix<double, 1, Dynamic> &x1, Matrix<double, 1, Dynamic> &val, const Matrix<double, 1, Dynamic> &x2 );
    
        static Matrix< std::complex<double>, 1, Dynamic > FreqResp( const Polynomial &p, const Matrix< std::complex<double>, 1, Dynamic> &freqs );
    
        //static Polynomial InvFreq(    const Matrix< std::complex<double>, 1, Dynamic > &FreqRes, const Matrix<double, 1, Dynamic > &Freq, unsigned int nNum, unsigned int nDen, Matrix< std::complex<double>, 1, Dynamic >   Weight );
    
        static StateSpace Tf2SS( const Polynomial& poly );
        
        //static Polynomial InverseFrequency( const Matrix< std::complex<double>,1,Dynamic> &K, const Matrix< double,1,Dynamic>& w, unsigned int order_num, unsigned int order_den, bool DisciriminateEdges=false);
        static Polynomial VectorFitting( Matrix< std::complex<double>,1, Dynamic> K, Matrix<double,1, Dynamic> w, int i, int dummy, bool UseHanWindow );
        IDENTIFICATION_TYPE Type;
        PotentialData *data = nullptr;
        std::vector<Matrix<double,6, 6> > Ainf;
    };
}
 
#endif