How to properly use the ifft function ?
20 次查看(过去 30 天)
显示 更早的评论
Hello, I'm trying to apply the ifft function to the Fourier transform of a real signal S I know, so as to get S in the end. (basically ifft(ffft(S)) = S).
My fft was made with a number N of points ( fft(S,N) ) : N = 2^nextpow2(Ns), where Ns is the number of points in S ( Ns = length(S) ).
I don't know what operations (symmetry, normalization, scaling ...) I have to do before applying ifft but for the moment I don't get S. And the signal I get has a non-neglegtable imaginary part.
Thank you for your help, Léo.
0 个评论
回答(4 个)
Wayne King
2011-10-20
x = randn(1000,1);
xdft = fft(x);
y = ifft(xdft,'symmetric');
max(abs(x-y))
The 'symmetric' option enforces that the imaginary part should be zero in the inverse DFT.
0 个评论
Daniel Shub
2011-10-20
Working from Wayne's answer:
x = randn(1000,1);
xdft = fft(x, pow2(nextpow2(length(x))));
y = ifft(xdft);
isreal(y)
I don't need to use the symmetric option and I get a real y. I cannot do x-y, because y is not the same length as x.
Can you show your code Leonard?
0 个评论
Walter Roberson
2011-10-20
It is not uncommon for ifft(fft(S)) to have an imaginary part, due to round-off error.
1 个评论
Daniel Shub
2011-10-20
Really? Is it because the fft of a real signal results in a non-symmetric complex signal, the ifft of a symmetric signal results in a non-real signal, or both? I understand there is round-off error, I am surprised it is asymmetric and that the FFT algorithm doesn't prevent this. Can you provide an "x" that will demonstrate this?
Temur Dzigrashvili
2022-3-11
Is it possible to approximate a complex-valued function of real variable in Matlab.
Values of function (complex numbers) and arguments (real numbers) are given.
0 个评论
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Digital Filter Analysis 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!