example_regression_fitting.cpp

Programa para el testeo de las funciones.

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

g++ -static -o example_regression_fitting example_regression_fitting.cpp  -lpdsmlmm -lpdsramm -lpdsspmm

Para executar o programa:

./example_regression_fitting

Retornando por consola:

C.T():
-0.3 0.1  
C.T():
0.99314856871846 2.0109980511082 

Código example_regression_fitting.cpp:

 
#include <Pds/Ra>
 
double Func(const Pds::Vector &x,const Pds::Vector &c)
{
    if(c.IsEmpty()) return 0;
    return c[0]+c[1]*x[0];
}
 
Pds::Vector DcFunc(const Pds::Vector &x,const Pds::Vector &c)
{
    if(c.IsEmpty()) return Pds::Vector();
    
    Pds::Vector Dc(c.Nel());
 
    Dc.SetRaw(0,1.0);
    Dc.SetRaw(1,x[0]);
    return Dc;
}
 
 
int main(void)
{
    unsigned int L=1000;
    Pds::Matrix X=Pds::LinSpace(0,1,L);
    Pds::Vector Y=2*X+1+0.1*Pds::RandN(L,1);
    
    Pds::Vector Yres;
    
    Pds::Vector C="-0.3 0.1";
    
    C.T().Print("C.T():\n");
    C=Pds::Regression::FitModel(X,Y,Func,DcFunc,C);
    C.T().Print("C.T():\n");
 
    Yres=X.OperateRows(Func,C);
    
    Pds::Octave::MarkerSize=1;
    Pds::Octave::Plot::CurveXYXY(X,Y,X,Yres,"testando.m","testando.png");
 
}