| 查看: 3969 | 回复: 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操作呢？ |
回复此楼» 猜你喜欢
300求调剂
已经有9人回复
-
291求调剂
已经有7人回复
-
求调剂推荐
已经有7人回复
-
289 分105500药学专硕求调剂(找B区学校)
已经有4人回复
-
0854求调剂
已经有13人回复
-
初试324 中药学 一志愿天中医 求调剂
已经有3人回复
-
药学求调剂
已经有14人回复
-
327求调剂
已经有27人回复
-
急需调剂
已经有5人回复
-
271求调剂
已经有33人回复
