Matlab 一维向量fftshift函数 C++实现

fftshift网上有很多二维的实现，但是对于提高精度的另一种办法就是把二维向量看成是一维一维的向量，也可以这么做fftshift。

博主这里提供一个简单的算法实现，希望大家喜欢！

1.针对于向量元素个数为奇数个，这里以25为例：

  

   

    

     

    
    

     

      void FFTshift1(double (&xin)[1][25], int h, int l)
     
    
   

    

     

    
    

     

      {
     
    
   

    

     

    
    

     
	double xinf1[1][13];//前半段
     
    
   

    

     

    
    

     
	double xinf2[1][12];//后半段
     
    
   

    

     

    
    

     
	//for (int i = 0; i < h; i++) {
     
    
   

    

     

    
    

     
		for (int j = 0; j < l/2+1; j++) {
     
    
   

    

     

    
    

     

      			xinf1[0][j] = xin[0][j];
     
    
   

    

     

    
    

     

      		}
     
    
   

    

     

    
    

     
		for (int j = l/2+1; j < l; j++) {
     
    
   

    

     

    
    

     

      			xinf2[0][j-l/2-1] = xin[0][j];
     
    
   

    

     

    
    

     

      		}
     
    
   

    

     

    
    

     
		for (int j = 0; j < l / 2; j++) {
     
    
   

    

     

    
    

     

      			xin[0][j] = xinf2[0][j];
     
    
   

    

     

    
    

     

      		}
     
    
   

    

     

    
    

     
		for (int j = l / 2; j < l; j++) {
     
    
   

    

     

    
    

     

      			xin[0][j] = xinf1[0][j-l/2];
     
    
   

    

     

    
    

     

      		}
     
    
   

    

     

    
    

     
	//}
     
    
   

    

     

    
    

     

      }
     
    
  
 

因为这里h代表只有1行，所以可以注释掉；

说白了这个函数就是把前半段和后半段换一个位置，但是奇数的时候前半段长度是总长度除以2向上取整，后一半是向下取整，偶数的直接前后交换即可，这里以20为例：

  

   

    

     

    
    

     

      void FFTshift2(double (&xin)[1][20], int h, int l)
     
    
   

    

     

    
    

     

      {
     
    
   

    

     

    
    

     
	double temp;
     
    
   

    

     

    
    

     
	for (int j = 0;j<l/2;j++) {
     
    
   

    

     

    
    

     

      		temp = xin[0][j+l/2];
     
    
   

    

     

    
    

     

      		xin[0][j+l/2] = xin[0][j];
     
    
   

    

     

    
    

     

      		xin[0][j] = temp;
     
    
   

    

     

    
    

     

      	}
     
    
   

    

     

    
    

     

      }
     
    
  
 

是不是很简单~

文章来源: nickhuang1996.blog.csdn.net，作者：悲恋花丶无心之人，版权归原作者所有，如需转载，请联系作者。

原文链接：nickhuang1996.blog.csdn.net/article/details/79578806