| 查看: 4075 | 回复: 0 | |||
xiashaoyan新虫 (正式写手)
|
[求助]
fftshift和ifftshift的使用场合,求助大神，这个问题已经困扰很久了~~
|
|
我知道，fftshift和ifftshift对于偶数长度的数组，两者无区别；对于奇数长度数组，fftshift相当于把数组的第一个放到数组中间，而ifftshift则是把最后一个放到中间。但是在用fft进行滤波等处理时，什么时候用fftshift或ifftshift感到很糊涂了，因为这个问题已经遇到很多次了. 现在想实现对一组数据进行低通滤波，代码如下： Nf = 2^12; xAxisf = linspace(-10,10,Nf); % plot(xaxis, data); % FFT xsize = numel(data); Xf = zeros([1 Nf]); indices = Nf/2-floor(xsize/2):Nf/2-floor(xsize/2)+xsize - 1; Xf(indices) = data; % Xf = fftshift(Xf); Xf = fft(Xf); Xf = fftshift(Xf); % plot Xfa = abs(Xf); plot(xAxisf, Xfa); % generate super-gaussian filter function Nf = numel(Xf); widthfilter = 0.12; filterpower = 2; filter = exp(-(xAxisf.^2./widthfilter^2).^filterpower); % filter filtertimes = 20; Xf = Xf .* filter.^filtertimes; % plot Xfa = abs(Xf); plot(Xfa); % iFFt Xfs = ifftshift(Xf); Xif = ifft(Xfs); % Xif = ifftshift(Xif); result = abs(Xif); plot(result(indices)) 现在针对这个程序，我有一些疑问，希望大家能帮忙看看，谢谢大神门了！！ 问题1： 在fft用之前，不需要用fftshift对data进行操作吗？就像如下操作： Xf = fftshift(data); % NOT NEEDED？？？ Xf = fft(Xf); 问题2： 程序中是用ffitshift对Xf操作的(Xf = fftshift(Xf)),我想为什么不能用ifftshift操作呢？ 问题3： 在使用ifft之前，为什么用ifftshift操作呢？(Xfs = ifftshift(Xf)  而不是用ffishift呢？问题4; Xif = ifft(Xfs)反傅里叶变换之后，要用iffshift操作吗？如果需要，为什么不是用ffishift呢？ 其实，问题总的来说，就是fftshift和ifftshift在什么场合下使用，为什么有时候用fftshift操作，而不是用ifftshift操作呢？或者为什么有时候用ifftshift操作，而不是用fftshift操作呢？ |
回复此楼» 猜你喜欢
面上合作单位盖章
已经有9人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
-
为什么国自然不能直接公布
已经有5人回复
-
梦想
已经有7人回复
-
学科评审组评审是指会评吗？
已经有4人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
-
怎么看青基中了没有啊
已经有6人回复
-
国社科又开始会评了，不知道这次命运如何
已经有11人回复
投票:
30