Group Abstract Group Abstract

Message Boards Message Boards

5
|
4.4K Views
|
19 Replies
|
21 Total Likes
View groups...
Share
Share this post:
GROUPS:

[WSG25] Daily Study Group: Introduction to Laplace Transforms

Posted 6 months ago
POSTED BY: Juan Ortiz
19 Replies

Hi;

In trying to get a better understanding of complex numbers in the Laplace Transform, I have formed the following understanding from the webinar. Please correct me if my understanding is incorrect.

It seems that there is quite a bit of terminology associated with the complex numbers that were created during the Laplace Transform, with the first being “analytical”. My understanding of the meaning of analytical is that the derivative can be taken of the function. In other words, the function has no features which would inhibit taking a derivative, such as: cusps, jumps, esc. (just a smooth continuous function). Now sometimes in complex numbers, you encounter features such as poles (where the function runs off to infinity) or zeros (where you have a divide by zero – undefined), which these features make the function non-analytic. However, with these non-analytic functions, you can incorporate techniques such a pruning (which is piecewise definitions) to eliminate the areas of the function that are non-analytic (poles & zeros), thus making the function differentiable. If this is not correct, please tell me what I am misinterpreting.

Additionally, what causes the complex numbers of the Laplace Transform?

Thanks,

Mitch Sandlin

POSTED BY: Mitchell Sandlin

Dear Mitch:

We used 'analytical' essentially meaning that there was a closed formula for it, particularly in the numerical evaluation lesson to distinguish between the two.

The parameter s, in the integral definition of the Laplace transform is complex, and this is the variable in the Laplace transform. That is where complex numbers come from.

POSTED BY: Juan Ortiz
Posted 5 months ago
POSTED BY: Phil Earnhardt

Thank you Phil. This is a very nice real LIVE example.

POSTED BY: Juan Ortiz
Posted 5 months ago
POSTED BY: Phil Earnhardt

I asked in today's lecture but I stepped away for a moment and when I came back the lecture was over and never saw if I got a reply:

Why do you choose to use the 2D version of ComplexPlot? I find the rainbow colors to be difficult to interpret. The 3D version of the plot more clearly shows poles and zeroes and seems to provide a better understanding of the behavior of the functions. Specially with the ability to rotate it around and look at if from all angles.

Do you have a special reason for that? Or is it just a question of taste?

POSTED BY: Carl Hahn
POSTED BY: Juan Ortiz
POSTED BY: Carl Hahn
POSTED BY: Juan Ortiz
POSTED BY: Michael Ulrey
POSTED BY: Juan Ortiz

Thanks!

POSTED BY: Michael Ulrey
Posted 5 months ago
POSTED BY: Graham Gyatt

Dear Graham:

Both seem fine to me. For problem 6, I will suggest comparing your answer with the answers provided considering that s>0 in this case. For example, you can try the command FullSimplify[yourAnswer==answerCompared, s>0] or you can cheat a bit to very using InverseLaplaceTransform on each of the posible aswers provided. For 10, make sure to set the WorkingPrecision to 10 as an option in your call to LaplaceTransform.

POSTED BY: Juan Ortiz
POSTED BY: Carl Hahn

Dear Carl:

I believe that in signal analysis the Fourier transform is the way to go. Note that in the Laplace transform s is complex and in fact if h(t) is the Heaviside function and s=c+i y, then the Laplace transform of f(t) is equal to the Fourier transform of exp(-c t) h(t)f(t) depending on y. This is actually used to get the inversion formula for the Laplace transform using the Fourier transform. For more on this I can refer you to the documentation page of the FourierTransform function in the Wolfram language, where I included a subsection on signals and systems in the Applications section.

POSTED BY: Juan Ortiz

I noticed in the exercises for Part 6 that when WorkingPrecision (WP) was increased, the resulting numerical answer differed before the last digit. For example, in exercise 2, when WP->10, the result was 0.1979119699 and when WP was increased to 20, the answer was 0.19791196966502245964 yielding a difference of 3x10^10. Granted, not much of a difference, but I expected a difference of no greater than 1x10^-10. I suspect the difference is due to the accumulation of rounding "errors". There is discussion in the WP documentation that final results from internal calculations done to n-digit precision may have much lower precision. So I wonder if there are examples where a long series of calculations resulted in significant differences. If so, it would be useful to know what to watch out for. I need to review the Precision & Accuracy Control Guide in the WL documentation to see if any guidelines are suggested.

POSTED BY: James Kralik

Dear James:

WorkingPrecision->n causes all internal computations to be done to at most n-digit precision. This does not imply that the answer will be correct to that number of digits. The final results you get may have much lower precision.

Instead of comparing WorkingPrecision->10 with WorkingPrecision->20, we should compare each to infinite precision. This will require to have an analytic representation for the transform. Consider exercise 1 in the same lesson:

exact = LaplaceTransform[
   Cos[Sqrt[t^2 - 4]]/Sqrt[t^2 - 4] UnitStep[t - 2], t, s] /. 
  s -> 17/10

approximate = 
 LaplaceTransform[Cos[Sqrt[t^2 - 4]]/Sqrt[t^2 - 4] UnitStep[t - 2], t,
   1.7]

approximateWP10 = 
 LaplaceTransform[Cos[Sqrt[t^2 - 4]]/Sqrt[t^2 - 4] UnitStep[t - 2], t,
   1.7, WorkingPrecision -> 10]

approximateWP20 = 
 LaplaceTransform[Cos[Sqrt[t^2 - 4]]/Sqrt[t^2 - 4] UnitStep[t - 2], t,
   1.7, WorkingPrecision -> 20]

Now:

exact - approximate=4.50937*10^-10

exact - approximateWP10=-1.14043*10^-12

exact - approximateWP20=-3.46945*10^-18

POSTED BY: Juan Ortiz

Juan is a highly-experienced instructor and has worked hard to create this wonderful introduction to Laplace transforms and their many applications.

I strongly recommend this study group to everyone!

POSTED BY: Devendra Kapadia
Reply to this discussion
Community posts can be styled and formatted using the Markdown syntax.
Reply Preview
Attachments
Remove
or Discard