fftwf频域相乘

最新推荐文章于 2025-05-22 08:00:00 发布

原创最新推荐文章于 2025-05-22 08:00:00 发布 · 2.4k 阅读

0 ·

CC 4.0 BY-SA版权

c++优化专栏收录该内容

45 篇文章

订阅专栏

 bool HOG = true;
    bool FIXEDWINDOW = false;
    bool MULTISCALE = true;
    int w=10;
    int h=10;
    fftwf_complex * r   =
            (fftwf_complex*) fftwf_malloc(sizeof (fftwf_complex) *w*h );
    for (size_t i = 0; i < h; i++) {
        for (size_t j = 0; j < w; j++) {
                r[i*w+j][0]= i*j*0.1+1;
                r[i*w+j][1]= i*j*0.2+1;
        }
    }
    fftwf_complex * r2   =
            (fftwf_complex*) fftwf_malloc(sizeof (fftwf_complex) *w*h );
    for (size_t i = 0; i < h; i++) {
        for (size_t j = 0; j < w; j++) {
                r2[i*w+j][0]= i*j*0.5+2;
                r2[i*w+j][1]= i*j*0.9+2;
        }
    }
      SimpleTracker  tracker(HOG, FIXEDWINDOW, MULTISCALE );
fftwf_complex * res=tracker.mulSpectrums(r,r2,w,h,0);
for (size_t i = 0; i < h; i++) {
    for (size_t j = 0; j < w; j++) {
           printf("%f+%fi ", r2[i*w+j][0],r2[i*w+j][1]);
    }
     printf("\n");
}


fftwf_complex * SimpleTracker::mulSpectrums(fftwf_complex * x1, fftwf_complex *  x2,
                                            int w,int h,bool conj )
{
    fftwf_complex * r   =
            (fftwf_complex*) fftwf_malloc(sizeof (fftwf_complex) *w*h );
    for (size_t i = 0; i < h; i++) {
        for (size_t j = 0; j < w; j++) {
            if(conj){
                r[i*w+j][0]= x1[i*w+j][0]* x2[i*w+j][0] - x1[i*w+j][1]* (-x2[i*w+j][1]);
                r[i*w+j][1]= x1[i*w+j][0]* (-x2[i*w+j][1] )+    x1[i*w+j][1]* x2[i*w+j][0];
            }
            else{
                r[i*w+j][0]= x1[i*w+j][0]* x2[i*w+j][0] -     x1[i*w+j][1]* x2[i*w+j][1] ;
                r[i*w+j][1]= x1[i*w+j][0]* x2[i*w+j][1] +    x1[i*w+j][1]* x2[i*w+j][0];
            }
        }
    }
    return r;
}