| 查看: 4034 | 回复: 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操作呢？ |
回复此楼» 猜你喜欢
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有3人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有4人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有4人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有4人回复
-
不要再数国自然申请书的 filecode 的分隔符个数了
已经有15人回复
-
祈祷青基必中
已经有16人回复
-
生命口会评
已经有12人回复
-
大龄残疾硕士的一点执念
已经有26人回复
-
27届辽宁大学应届毕业生申博
已经有3人回复
-
关于如何从代码看上不上会
已经有23人回复