example_goertzel_create.cpp

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

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

Para executar o programa:

./example_goertzel_create

Retornando por consola:

FT{X}[0]: 3.63798e-12
FT{X}[1]: 0.266025
FT{X}[2]: 0.664383
FT{X}[3]: 1.70215
FT{X}[4]: 255.724
FT{X}[5]: 2.24209
FT{X}[6]: 1.20616
FT{X}[7]: 0.851853

Código example_goertzel_create.cpp:

 
#include <Pds/Ca>
#include <Pds/Ra>
    
int main(void)
{
    unsigned int N=512;
    unsigned int K=8;
    unsigned int k;
    double fc=4;
    
    Pds::Vector T=Pds::LinSpace(0,1,N);
    Pds::Vector X=Pds::Sin(2*Pds::Ra::Pi*fc*T);
    
    Pds::Goertzel FT[K];
    Pds::Complex Out[K];
    
    for(k=0;k<K;k++)
    FT[k].Configure(k,N);
    
    for(k=0;k<K;k++)
    Out[k]=FT[k].Evaluate(X);
    
    for(k=0;k<K;k++)
    std::cout<<"FT{X}["<<k<<"]: "<<Out[k].Module()<<std::endl;
    
}