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;
}
}
是不是很简单~