Wednesday 20 June 2012

Laplace Transforms


In this chapter we will be looking at how to use Laplace transforms to solve differential equations.  There are many kinds of transforms out there in the world.  Laplace transforms and Fourier transforms are probably the main two kinds of transforms that are used.  As we will see in later sections we can use Laplace transforms to reduce a differential equation to an algebra problem.  The algebra can be messy on occasion, but it will be simpler than actually solving the differential equation directly in many cases.  Laplace transforms can also be used to solve IVP’s that we can’t use any previous method on.

For “simple” differential equations such as those in the first few sections of the last chapter Laplace transforms will be more complicated than we need.  In fact, for most homogeneous differential equations such as those in the last chapter Laplacetransforms is significantly longer and not so useful.  Also, many of the “simple” nonhomogeneous differential equations that we saw in the Undetermined Coefficients and Variation of Parameters are still simpler (or at the least no more difficult than Laplace transforms) to do as we did them there.  However, at this point, the amount of work required for Laplacetransforms is starting to equal the amount of work we did in those sections.

Laplace transforms comes into its own when the forcing function in the differential equation starts getting more complicated.  In the previous chapter we looked only at nonhomogeneous differential equations in which g(t) was a fairly simple continuous function.  In this chapter we will start looking at g(t)’s that are not continuous.  It is these problems where the reasons for using Laplace transforms start to become clear.

We will also see that, for some of the more complicated nonhomogeneous differential equations from the last chapter,Laplace transforms are actually easier on those problems as well.

Here is a brief rundown of the sections in this chapter.

The Definition  The definition of the Laplace transform.  We will also compute a couple Laplace transforms using the definition.

Laplace Transforms  As the previous section will demonstrate, computing Laplace transforms directly from the definition can be a fairly painful process.  In this section we introduce the way we usually compute Laplacetransforms.

Inverse Laplace Transforms  In this section we ask the opposite question.  Here’s a Laplace transform, what function did we originally have?

Step Functions  This is one of the more important functions in the use of Laplace transforms.  With the introduction of this function the reason for doing Laplace transforms starts to become apparent.

Solving IVP’s with Laplace Transforms  Here’s how we used Laplace transforms to solve IVP’s.

Nonconstant Coefficient IVP’s  We will see how Laplace transforms can be used to solve some nonconstant coefficient IVP’s

IVP’s with Step Functions  Solving IVP’s that contain step functions.  This is the section where the reason for using Laplace transforms really becomes apparent.

Dirac Delta Function  One last function that often shows up in Laplace transform problems.

Convolution Integral  A brief introduction to the convolution integral and an application for Laplace transforms.

Table of Laplace Transforms  This is a small table of Laplace Transforms that we’ll be using here.

The Definition

You know, it’s always a little scary when we devote a whole section just to the definition of something.  Laplace transforms (or just transforms) can seem scary when we first start looking at them.  However, as we will see, they aren’t as bad as they may appear at first.

Before we start with the definition of the Laplace transform we need to get another definition out of the way.

A function is called piecewise continuous on an interval if the interval can be broken into a finite number of subintervals on which the function is continuous on each open subinterval (i.e. the subinterval without its endpoints) and has a finite limit at the endpoints of each subinterval.  Below is a sketch of a piecewise continuous function.
pw_cts
In other words, a piecewise continuous function is a function that has a finite number of breaks in it and doesn’t blow up to infinity anywhere.

Now, let’s take a look at the definition of the Laplace transform.

Definition
Suppose that f(t) is a piecewise continuous function.  The Laplace transform of f(t) is denoted  and defined as
                                                                                        (1)

There is an alternate notation for Laplace transforms.  For the sake of convenience we will often denote Laplace transforms as,


With this alternate notation, note that the transform is really a function of a new variable, s, and that all the t’s will drop out in the integration process.

Now, the integral in the definition of the transform is called an improper integral and it would probably be best to recall how these kinds of integrals work before we actually jump into computing some transforms.

Example 1  If , evaluate the following integral.
                                                                  
Solution
Remember that you need to convert improper integrals to limits as follows,
                                                      

Now, do the integral, then evaluate the limit.
                                                    

Now, at this point, we’ve got to be careful.  The value of c will affect our answer.  We’ve already assumed that c was non-zero, now we need to worry about the sign of c.  If c is positive the exponential will go to infinity.  On the other hand, if c is negative the exponential will go to zero.

So, the integral is only convergent (i.e. the limit exists and is finite) provided c<0.  In this case we get,
                                                                                (2)

Now that we remember how to do these, let’s compute some Laplace transforms.  We’ll start off with probably the simplest Laplace transform to compute.

Example 2  Compute L{1}.

Solution
There’s not really a whole lot do here other than plug the function f(t) = 1 into (1)
                                                           

Now, at this point notice that this is nothing more than the integral in the previous example with .  Therefore, all we need to do is reuse (2) with the appropriate substitution.  Doing this gives,
                             

Or, with some simplification we have,
                                                 

Notice that we had to put a restriction on s in order to actually compute the transform.  All Laplace transforms will have restrictions on s.  At this stage of the game, this restriction is something that we tend to ignore, but we really shouldn’t ever forget that it’s there.

Let’s do another example.

Example 3  Compute 

Solution
Plug the function into the definition of the transform and do a little simplification.
                                            

Once again, notice that we can use (2) provided .  So let’s do this.
                                  

Let’s do one more example that doesn’t come down to an application of (2).

Example 4  Compute L{sin(at)}.

Solution
Note that we’re going to leave it to you to check most of the integration here.  Plug the function into the definition.  This time let’s also use the alternate notation.
                                             

Now, if we integrate by parts we will arrive at,
                             

Now, evaluate the first term to simplify it a little and integrate by parts again on the integral.  Doing this arrives at,
           

Now, evaluate the second term, take the limit and simplify.
             

Now, notice that in the limits we had to assume that s>0 in order to do the following two limits.
                                                         

Without this assumption, we get a divergent integral again.  Also, note that when we got back to the integral we just converted the upper limit back to infinity.  The reason for this is that, if you think about it, this integral is nothing more than the integral that we started with.  Therefore, we now get,
                                                        

Now, simply solve for F(s) to get,
                               

As this example shows, computing Laplace transforms is often messy.

Before moving on to the next section, we need to do a little side note.  On occasion you will see the following as the definition of the Laplace transform.


Note the change in the lower limit from zero to negative infinity.  In these cases there is almost always the assumption that the function f(t) is in fact defined as follows,


In other words, it is assumed that the function is zero if t<0.  In this case the first half of the integral will drop out since the function is zero and we will get back to the definition given in (1).  A Heaviside function is usually used to make the function zero for t<0.  We will be looking at these in a later section.
As we saw in the last section computing Laplace transforms directly can be fairly complicated.  Usually we just use a table of transforms when actually computing Laplace transforms.  The table that is provided here is not an inclusive table, but does include most of the commonly used Laplace transforms and most of the commonly needed formulas pertaining toLaplace transforms.

Before doing a couple of examples to illustrate the use of the table let’s get a quick fact out of the way.

Fact
Given f(t) and g(t) then,
                                                  
for any constants a and b.

In other words, we don’t worry about constants and we don’t worry about sums or differences of functions in takingLaplace transforms.  All that we need to do is take the transform of the individual functions, then put any constants back in and add or subtract the results back up.

So, let’s do a couple of quick examples.

Example 1  Find the Laplace transforms of the given functions.
(a)    [Solution]
(b)    [Solution]
(c)    [Solution]
(d)    [Solution]

Solution
Okay, there’s not really a whole lot to do here other than go to the table, transform the individual functions up, put any constants back in and then add or subtract the results.

We’ll do these examples in a little more detail than is typically used since this is the first time we’re using the tables.

(a) 





(b) 

(c) 

(d) 

Make sure that you pay attention to the difference between a “normal” trig function and hyperbolic functions.  The only difference between them is the “+ a2” for the “normal” trig functions becomes a “- a2” in the hyperbolic function!  It’s very easy to get in a hurry and not pay attention and grab the wrong formula.  If you don’t recall the definition of the hyperbolic functions see the notes for the table.

Let’s do one final set of examples.

Example 2  Find the transform of each of the following functions.
(a)    [Solution]
(b)    [Solution]
(c)    [Solution]
(d)    [Solution]
(e)    [Solution]
Solution
(a)  

This function is not in the table of Laplace transforms.  However we can use #30 in the table to compute its transform.  This will correspond to #30 if we take n=1.

So, we then have,
                                       
Using #30 we then have,
                                                                 
(b) 

This part will also use #30 in the table.  In fact we could use #30 in one of two ways.  We could use it with .
                              

Or we could use it with .
                               

Since it’s less work to do one derivative, let’s do it the first way.  So using #9 we have,
                                       
The transform is then,
                                                                
(c) 

This part can be done using either #6 (with  ) or #32 (along with #5).  We will use #32 so we can see an example of this.  In order to use #32 we’ll need to notice that
                                    
Now, using #5,
                                                   
we get the following.
                                                        

This is what we would have gotten had we used #6.

(d) 

For this part we will use #24 along with the answer from the previous part.  To see this note that if
                                                                       
then
                                                                   

Therefore, the transform is.
                                                             

(e)  

This final part will again use #30 from the table as well as #35.
                                                   

Remember that g(0) is just a constant so when we differentiate it we will get zero!

As this set of examples has shown us we can’t forget to use some of the general formulas in the table to derive new Laplace transforms for functions that aren’t explicitly listed in the table!



No comments:

Post a Comment