example_vector_poly.cpp

Programa para el testeo de las funciones.

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

g++ -static -o example_vector_poly example_vector_poly.cpp -lpdsramm 

Para executar o programa:

./example_vector_poly

Retornando por consola:

X.Size(): 1000  1
Y.Size(): 1000  1

P=Pds::PolyFit(X,Y,2);  P.T():
-0.0053770932721563 0.017219023067994   0.98828000248775    
X.T():  1  -1 2  -2 
P.T():  1 0 1 
Pds::PolyVal(P,X).T():  2 2 5 5 

U:
0 1 2 
3 4 5 
6 7 8 

Pds::PolyVal(P,U):
1  2  5  
10 17 26 
37 50 65 
M:
1  1  1  1  
1  2  4  8  
1  -1 1  -1 
1  3  9  27 
Pd.T(): 0 2 
Pm.T(): 1 3 3 1 
Ps.T(): 2 3 1 
Acc.T():        1 2 1 
Acc.T():        2 1 1 
Acc.T():        3 2 3 1 
Pp.T(): 3  -6 4

Código example_vector_poly.cpp:

 
#include <cmath>    
#include <Pds/Ra>
    
int main(void)
{
    unsigned int L=1000;
    
    // Creating Vectors X in [0,1> and Y: Y=X^2+noise
    Pds::Vector X=Pds::LinSpace(0,1,L);
    Pds::Vector Y=Pds::Pow(X,2)+0.05*Pds::RandN(L,1);
    
    // Print size of vectors
    X.Size().Print("X.Size(): ");
    Y.Size().Print("Y.Size(): ");
    
    // Getting the polynomial P of order 2 that fit (X,Y)
    Pds::Vector P=Pds::PolyFit(X,Y,2);
    P.T().Print("\nP=Pds::PolyFit(X,Y,2);\tP.T():\n");
    
    // Creating vector X
    X="1 -1 2 -2";
    X.T().Print("X.T():\t");
    // Creating vector X
    P="1 0 1";
    P.T().Print("P.T():\t");
    
    // Getting Y=p(X)
    Y=Pds::PolyVal(P,X);
    Y.T().Print("Pds::PolyVal(P,X).T():\t");
    
    // Creating vector Xa
    Pds::Matrix U("0 1 2\n3 4 5\n6 7 8");
    U.Print("\nU:\n");
    
    // Getting Ya=p(U)
    Pds::Matrix Ya=Pds::PolyVal(P,U);
    Ya.Print("\nPds::PolyVal(P,U):\n");
    
    X="1\n 2\n -1\n 3";
    Pds::Matrix M=Pds::PolyMat(X,3);
    M.Print("M:\n");
    
    Pds::Vector Pd=Pds::PolyDer(P,1);
    Pd.T().Print("Pd.T():\t");
    
    Pds::Vector P1="1 1";
    Pds::Vector P2="1 2 1";
    Pds::Vector Pm=Pds::PolyMul(P1,P2);
    Pm.T().Print("Pm.T():\t");
    
    Pds::Vector Ps=Pds::PolyAdd(P1,P2);
    Ps.T().Print("Ps.T():\t");
    
    Pds::Vector Acc;
    
    Pds::PolyAcc(Acc,"1 2 1");
    Acc.T().Print("Acc.T():\t");
    
    Pds::PolyAcc(Acc,"1 -1");
    Acc.T().Print("Acc.T():\t");
    
    Pds::PolyAcc(Acc,"1 1 2 1");
    Acc.T().Print("Acc.T():\t");
    
    Pds::Vector Pp=Pds::PolyPoly("1 1 1","1 -2");
    Pp.T().Print("Pp.T():\t");
    return 0;
}