refman/html/example__regression__fitting_8cpp_source.html

#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");


}

Pds::Matrix
La clase tipo Pds::Matrix . Esta clase genera una matriz de Nlin lineas y Ncol columnas....
Definition: Matrix.hpp:96

Pds::Vector
La clase tipo Pds::Vector . Esta clase genera una matriz de Nlin lineas y 1 columna....
Definition: Vector.hpp:80

Pds::RandN
Pds::Matrix RandN(unsigned int N)
Crea una matriz con datos aleatórios con distribución Normal de media cero 0.0 y desvío padrón 1....

Pds::LinSpace
Vector LinSpace(double start, double end, unsigned int N)
Retorna un vector con elementos desde start hasta end con N elementos.

Pds::Matrix::T
Pds::Matrix T(void) const
Transpuesta de si mismo (A), el resultado es cargado en B.

Pds::Matrix::OperateRows
Pds::Matrix OperateRows(double(*func)(const Pds::Matrix &Row)) const
Opera la función func usando como entrada cada fila de la matriz.

Pds::Matrix::IsEmpty
bool IsEmpty(void) const
Verifica si la matriz es nula es decir con lineas o columnas cero o arreglo NULL.

Pds::Matrix::Nel
unsigned int Nel(void) const
Retorna el número de elementos de la matriz (Nlin x Ncol).

Pds::Matrix::Print
void Print(const std::string &str, unsigned int precision) const
Imprime en pantalla el contenido de la matriz después del texto indicado en str.

Pds::Octave::Plot::CurveXYXY
bool CurveXYXY(const Pds::Vector &X, const Pds::Vector &Y, const Pds::Vector &Xs, const Pds::Vector &Ys, const std::string &octfilename, const std::string &imgfilename, bool EnableLine=false)
Plot de X vs. Y junto con Xs vs. Ys.

Pds::Octave::MarkerSize
unsigned int MarkerSize
Line width usado en el código octave. Por defecto:

Pds::Regression::FitModel
Pds::Vector FitModel(const Pds::Matrix &X, const Pds::Vector &Y, double(*Func)(const Pds::Vector &, const Pds::Vector &), Pds::Vector(*DcFunc)(const Pds::Vector &, const Pds::Vector &), Pds::Vector c0, double Lambda=0.01, double MinError=1.0e-5, unsigned int MaxIter=2048)
Esta función encaja un modelo multidimensional  com parametros  en un conjunto de muestras .