- #1
praban
- 13
- 0
Hello,
I am trying to numerically evaluate a convolution integral of two functions (f*g) using Fourier transform (FT) i.e using
FT(f*g) = FT(f) multiplied by FT(g) (1)
I am testing for a known case first. I have taken the gaussian functions (eq. 5, 6 and 7) as given in
http://mathworld.wolfram.com/Convolution.html
I am not using FFT but a naive implementation of FT. The problem I am facing is that the FTs are correct but convolution is not.
When I do inverse FT I do get back the 3 gaussians as given in the wolfram link. However, eq. (1) is not satisfied.
Is there any trick in eq. (1)? I would appreciate any suggestion.
thanks,
Praban
I am trying to numerically evaluate a convolution integral of two functions (f*g) using Fourier transform (FT) i.e using
FT(f*g) = FT(f) multiplied by FT(g) (1)
I am testing for a known case first. I have taken the gaussian functions (eq. 5, 6 and 7) as given in
http://mathworld.wolfram.com/Convolution.html
I am not using FFT but a naive implementation of FT. The problem I am facing is that the FTs are correct but convolution is not.
When I do inverse FT I do get back the 3 gaussians as given in the wolfram link. However, eq. (1) is not satisfied.
Is there any trick in eq. (1)? I would appreciate any suggestion.
thanks,
Praban