Laplace Transform of f(t) Related to smoothed f(t)?
When reading (Comment#18839) I started to wonder if there was a relationship between the Fourier Transform of a smoothed signal and the Laplace transform. I assumed there was a relationship (Comment#18854). After further derivation, I recommenced that if the goal is to derive the Laplace transform from the Fourier transom of the filtered signal:
1) The signal be properly windowed.
2) The FFT of the windowed Fourier Transform, needs to be compensated for the frequency effects that resulted from the low pass filter.
Weather it is a good idea to compute the Laplace transform from a windowed FFT of a filtered signal is outside of the scope of this thread (but feel free to comment bellow) .
The Laplace transform is given by:
1)
The Fourier transform is given by:
2)
The Two Sided Laplace Transform is given by:
3)
Therefor the Fourier transform is the two sided Laplace transform evaluated at
Returning to the one sided Laplace transform:
4)
5)
Let:
6)
7)
8)
where the low pass filtered version of f(t):
9)
and is the convolution of f(t) and the impulse response of a filter (or atleast aproximatly so) with bandwidth .
Plugging this result into integration by parts gives:
10)
or equivalently:
11)
The first two terms show how the endpoints chosen effect the transfrom. These two terms will cancel for a given frequency if the distance between the endpoints is some multiple of the period. The last term is the Fouier transform of the smoothed function with the frequencies weighted by and using a windowing function
(note the multiple is there because the Fourier transform variable is the Laplace transform variable but rotated by 90 degrees.)
The effect of the windowing function is to smooth the frequency response. This is because multiplication in the time domain is equivalent to convolution in the frequency domain. The following Fourier transform relationship is useful (relationship 205):
12)
Note, that if a non causal filter was used for the smoothing the relationship is much simpler.
13)
In both cases to properly deal with the end points the time shifting property of the Fourier transform is needed (relationship 102):
14)
Applying this property to the last two relationships gives:
15)
16)
Strictly dealing with the case where a causal filter is used and applying the rule for the Fourier transform of a convolution (Rule 109) we obtain:
17)
of equivalently:
18)
Some Comments:
1) If is negative the system is causal, and the filtered version of the signal will be causal.
2) Computing the smoothed signal does not save any computations with regards to computing the Laplace transform.
3) The derivation seems to show that their is a relationship between Laplace trancform and a windowed Fouier transform of the filtered signal.
4) To compute the Laplace transform based on the orginal signal use equation (5). To compute it based on the filtered signal use equation (11).
8 Comments »
Leave a reply to s243a Cancel reply
-
Recent
- Laplace Transform Via Limits
- log(CO2) and Scary Graphs
- Numeric Solutions to The Heat Equation
- Coriolis Forces in Hopkins and Simmons Vorticity Equation
- The Cross Product in Non Orthogonal Coordinate Systems
- Lagrangian Mechanics and The Heat Equation
- Laplace Transform of f(t) Related to smoothed f(t)?
- Coriolis Forces
- Vector Operations in Hoskins and Simmons Coordinates
- API/Object Viewers/Memory Mapping/
- Defining a Microsoft access Datasource
- Fractal Modeling of Turbulence
-
Links
-
Archives
- July 2012 (1)
- September 2009 (5)
- August 2009 (19)
- March 2009 (2)
-
Categories
-
RSS
Entries RSS
Comments RSS
Hmmmmmm………..I forgot about pingbacks, maybe I should only link to other blogs after most of the editing is done.
Comment by s243a | August 30, 2009 |
As a matter of contour integration, if you take the real axis for a Fourier transform, and move it so that the ends fold up along the imaginary axis, but the contour still passes near or below the origin, the result is the difference between two Laplace transforms. Often, because of the effect of a sqrt sign, the two parts ends up being equal and adding.
But I don’t think you’ll find any magic bullet in the FFT, unless you can reverse that contour manoeuvre to get an integral suitable for the purpose.
Comment by Nick Stokes | August 31, 2009 |
It wasn’t my intent to use the FFT to compute the Laplace transform along the real axis. (However, your idea sounds interesting). Rather, I wish to use the computational efficiency of the FFT to compute the Laplace transform along lines which are parallel to the imaginary axis in Laplace space. Out of the above equations equation (11) is the most suitable for computing this based on the filtered signal and equation (5) is most suitable to compute this based on the original signal.
There is of course a lot of algebra to check here and I was still making corrections after you posted (or at least up to when you posted). The thing that bugs me the most about my derivation is that factor I have as a coefficient for the last term in equation (11). I think this should go away as approach zero but it doesn’t.
Comment by s243a | August 31, 2009 |
I think my next thread related to this topic will be about compensating for the effects of a finite data set when attempting to evaluate the Laplace transform. The important relationship is the Fourier transform of the unit step function:
This relationship may be why I’m getting that factor of in the last term of equation (11)
(Relationship 13)
Comment by s243a | August 31, 2009 |
I know why I didn’t study numerical methods of computing the Laplace transform in school. I didn’t study it because engineers don’t study measure theory :
http://en.wikipedia.org/wiki/Laplace_transform#Fourier_transform
So I guess I have to learn how Lebesque Integration, relates to the Laplace transform.
Comment by s243a | August 31, 2009 |
With regards to my last comment here are some relevent links:
This first link looks pretty complete. I’ll tell you later how easy (or hard) it is to follow for me:
Does there exist the Lebesgue measure in the infinite-dimensional Space
(Last two links are books, so you can’t view the whole book online)
Borel-Laplace Transform and Asymptotic Theory: Introduction to
http://www.worldscibooks.com/etextbook/p245/p245_chap1.pdf
Comment by s243a | August 31, 2009 |
I found a good link for learning mathematics. It is from the university of Colorado.
http://www.uccs.edu/~math/vidarchive.html
The site contains free video’s of the lectures. Apparently the course on real analysis is mostly about measure theory. I found out about it when reading this thread from physics forums.
http://www.physicsforums.com/showthread.php?p=2342489#post2342489
Comment by s243a | September 10, 2009 |
Just posting some old links when I was investigating this stuff further. Haven’t looked at it in a while:
http://www.physicsforums.com/showthread.php?t=334191
http://www.physicsforums.com/showthread.php?t=338132
Comment by s243a | August 19, 2010 |