用C++实现Matlab中一维向量的fftshift函数
最编程
2024-08-14 15:01:39
...
【直播预告】上云 vs 下云:降本增笑?割韭菜?
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;
}
}
是不是很简单~
本文同步分享在 博客“悲恋花丶无心之人”(****)。
如有侵权,请联系 support@oschina.cn 删除。
本文参与“OSC源创计划”,欢迎正在阅读的你也加入,一起分享。
上一篇: 学习Python图像处理的必备技能:傅里叶变换原理及实现详解
下一篇: np.fft