houdini fft模型

2023-12-29 16:48:29

按照博友给的代码,进行可视化,发现绘出的图形没有错误,但计算的结果16,32,64.。。居然是错的;哪位知道正确的,敬请指正,或是后续啥时候复习完后再改正吧

主要函数代码



i@c_lay = detail(1,'iteration',0);
//printf("current _lay: %d ",@c_lay);

/*
float wn_r[@len];
float wn_i[@len];

//init wn
for(int k=0;k<@len;k++)
{
    wn_r[k] = cos(2*3.14/@len*k);
    wn_i[k] = -1*sin(2*3.14/@len*k);
    
}
*/



int len = @numpt;
i@t_lay = @lays;  //log2 8
//printf("%d",@numpt);



i@dist = shl(1,@c_lay);        //计算每一级的间隔距离
 
for( int j=0;j<@dist;j++ )      
{

    for(int k=j;k<len;k+=(shl(@dist,1)))   //k计算的序号1   
    {
        int q = k+@dist; //q计算的序号2
        int z = shl(k, (@t_lay - @c_lay -1));    //确定旋转因子的指数
        
        //w = rotation_factor(len,1,z);

        f@w_real = cos(2* PI * z / len);
        f@w_img = -sin(2* PI * z / len);    
        
        f@q_real = point(0,'real',q);
        f@q_img = point(0,'img',q);
        
        f@k_real = point(0,'real',k);
        f@k_img = point(0,'img',k);
     //   printf("k_real for %d ",@c_lay);
     //   printf("%f  ",@k_real);
                        
        float tmp_real =@k_real;
        float tmp_img = @k_img;
        
        //complex_multi x
        f@wxq_real = @w_real*@q_real -@w_img*@q_img;
        f@wxq_img = @w_real*@q_img + @w_img*@q_real;
        
        
        f@wxq_real = @wxq_real * @q_real - @wxq_img * @q_img;
        f@wxq_img =  @wxq_real * @q_img + @wxq_img * @q_real;
      
        
        
        
        //complex_add  x[k]
        tmp_real = @k_real + @wxq_real;
        tmp_img  = @k_img +  @wxq_img;
        
        setpointattrib(0,'real',k,tmp_real);
        setpointattrib(0,'img',k,tmp_img);
        
        //printf("%f",tmp_img);
        
        //complex sub  x[q]
        tmp_real =@k_real-@wxq_real;
        tmp_img = @k_img -@wxq_img;
        
        setpointattrib(0,'real',q,tmp_real);
        setpointattrib(0,'img',q,tmp_img);
        
        //printf("%f",tmp_img);
         //add prim
         
         int plen = detail(0,"len",0);
              v@pk1 = point(0,'P',plen+k);
              v@pq1 = point(0,'P',plen+q);
              
              v@pk2 = point(0,'P',k);
              v@pq2 = point(0,'P',q);  
              
          //    i@p1 = addpoint(0,@pk1);
          //    i@p2 = addpoint(0,@pq1);
              
              addprim(0,'polyline',k,plen+q);
              addprim(0,'polyline',q,plen+k);
        
       

    }
}

每次循环中进行两次计算(一次算两)

8点 = log2(8) = 3?分层? ;32点? 5层,64点 6层

每层点的跨度?dist << n(n代表层)

文章来源:https://blog.csdn.net/qq_39239990/article/details/135290459
本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。