example_predictor_kernel_fourier_series.cpp

Programa para el testeo de las funciones.Para compilar o código example_predictor_kernel_fourier_series.cpp:

g++ -static -o example_predictor_kernel_fourier_series example_predictor_kernel_fourier_series.cpp -lpdsmlmm 

Para executar o programa:

./example_predictor_kernel_fourier_series

Retornando por consola:

W:
0.00033443141643536     -0.0013788569187263     -0.002951637702057      -0.0017915203641805     -0.0018105859720968     0.0025881890045695      -0.0004863605736278     0.00018092888756878     -0.0016466673924326     -0.0049004521365083     -0.0014126034500854     -0.0010221018171096     0.39753973539033        -0.0012625166580139     0.15153048098468        -0.005054085903324      0.253549559735180.00149669658349        0.15155719808668        -0.0017099977179967     0.10083484544658
Training :: R2 : 0.98082060837963
Testing  :: R2 : 0.98231113063073

Ajuste de curva

Código example_predictor_kernel_fourier_series.cpp:

 
#include <iostream>
#include <Pds/Ra>
#include <Pds/Ml>
 
    
int main(void)
{
    Pds::Vector Yp;
    Pds::IterationConf Conf;
    
    // Generating data
    unsigned int L=1000;
    unsigned int K=10;
    Pds::Matrix X;
    Pds::Matrix F;
    Pds::Vector Y;
    Pds::DataSet::Regression::LoadDataFourierSeries(L,X,Y);
    
    // Create Predictor
    F=Pds::Kernel::FourierSeries(X,Pds::Vector("2.0"),K);
    Pds::Predictor Data(Conf,F,Y);
    Data.Print("\nW:\n");
    
    // Evaluate training data
    Yp=Data.Evaluate(F);
    
    // Testing W with training data set
    std::cout<<"Training :: R2 : "<<Y.R2(Yp)<<std::endl;
    
    // Evaluate testing data
    Pds::DataSet::Regression::LoadDataFourierSeries(L,X,Y);
    
    
    F=Pds::Kernel::FourierSeries(X,Pds::Vector("2.0"),K);
    Yp=Data.Evaluate(F);
 
    // Testing W with training data set
    std::cout<<"Testing  :: R2 : "<<Y.R2(Yp)<<std::endl;
    
    Pds::Octave::XLabel="x1";
    Pds::Octave::YLabel="x2";
    Pds::Octave::Plot::PointsXYXY(X,Y,X,Yp,"testando.m","example_predictor_kernel_fourier_series.png");
    
    return 0;
}