Re: Fourier Transform challenge
- From: Michael Soyka <msoyka_nospam@xxxxxxxxxxxxxxxxxxxxxxxxxx>
- Date: Fri, 02 Dec 2005 21:02:05 -0500
John wrote:
MATLAB:
fourier(exp(-x-exp(-x)))
ans =
gamma(1+i*w)
This "solution" is actually meaningless because
gamma(1+iw) = iw gamma(iw)
and the gamma function is analytic only in the right half plane (real part of z positive).
That said, the original poster's function exp(-x-exp(-x)) is itself integrable- its integal is
exp(-exp(-x))
as can be verified by differentiation.
Given this result, the Fourier Transform can be integrated "by parts". Start by integrating only over the interval [-T,T]. Then, integrating by parts and taking the limit as T becomes large we get
exp(-iwT-exp(-T)) - exp(iwT-exp(T))
FT{} = lim -------------------------------------
T->oo 1 + iwClearly the second term in the numerator converges to zero and so we are left with
exp(-iwT-exp(-T))
FT{} = lim -------------------
T->oo 1 + iw exp(-iwT)
= lim -----------
T->oo 1 + iwThis limit does not exist in the conventional sense but does exist as a "distribution" or "generalized function" (the Dirac delta function is the usual example- sorry to bring this subject up but I must). In this sense, the limit exists and is zero by the Riemann-Lebesgue Lemma (see Papoulis' book "The Fourier Integral and its Applications").
I suspect the limit being zero explains the excellent roll-off characteristics seen by the original poster.
Mike .
- Follow-Ups:
- Re: Fourier Transform challenge
- From: Michael Soyka
- Re: Fourier Transform challenge
- References:
- Fourier Transform challenge
- From: robert . w . adams
- Re: Fourier Transform challenge
- From: John
- Fourier Transform challenge
- Prev by Date: Re: Filter design.
- Next by Date: A question about bilinear xform and Butterworth LPF
- Previous by thread: Re: Fourier Transform challenge
- Next by thread: Re: Fourier Transform challenge
- Index(es):
Relevant Pages
|