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;
}