| 查看: 3942 | 回复: 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操作呢？ |
回复此楼» 猜你喜欢
遇见不省心的家人很难过
已经有24人回复
-
依托企业入选了国家启明计划青年人才。有无高校可以引进的。
已经有6人回复
-
博士延得我，科研能力直往上蹿
已经有9人回复
-
天津大学招2026.09的博士生，欢迎大家推荐交流（博导是本人）
已经有6人回复
-
AI 太可怕了，写基金时，提出想法，直接生成的文字比自己想得深远，还有科学性
已经有6人回复
-
有院领导为了换新车，用横向课题经费买了俩车
已经有9人回复
-
酰胺脱乙酰基
已经有13人回复
-
同年申请2项不同项目，第1个项目里不写第2个项目的信息，可以吗
已经有4人回复
-
有时候真觉得大城市人没有县城人甚至个体户幸福
已经有10人回复