Laplace Transforms (Part 6. Solving IVP’s with Laplace Transforms)
Posted by Suhaivi Hamdan on 8:18 PM with No comments
Solving IVP’s with Laplace Transforms
It’s now time to get back to differential equations. We’ve spent the last three sections learning
how to take Laplace transforms and how to take inverse Laplace
transforms. These are going to be
invaluable skills for the next couple of sections so don’t forget what we
learned there.
Before proceeding into differential equations we will need
one more formula. We will need to know
how to take the Laplace transform of a
derivative. First recall that f^{(n)} denotes the n^{th} derivative of the
function f. We now have the following fact.
Fact
Suppose that f, f’, f”,…f^{(n1)} are all
continuous functions and f^{(n)}
is a piecewise continuous function.
Then,

Since we are going to be dealing with second order
differential equations it will be convenient to have the Laplace
transform of the first two derivatives.


Notice that the two function evaluations that appear in
these formulas, and ,
are often what we’ve been using for initial condition in out IVP’s. So, this means that if we are to use these
formulas to solve an IVP we will need initial conditions at t = 0.
While Laplace transforms are particularly useful for
nonhomogeneous differential equations which have Heaviside functions in the
forcing function we’ll start off with a couple of fairly simple problems to
illustrate how the process works.
Example 1 Solve
the following IVP.
Solution
The first step in using Laplace
transforms to solve an IVP is to take the transform of every term in the
differential equation.
Using the appropriate formulas from our table of Laplace transforms gives us the
following.
Plug in the initial conditions and collect all the terms
that have a Y(s) in them.
Solve for Y(s).
At this point it’s convenient to recall just what we’re
trying to do. We are trying to find
the solution, y(t), to an IVP. What we’ve managed to find at this point is
not the solution, but its Laplace
transform. So, in order to find the
solution all that we need to do is to take the inverse transform.
Before doing that let’s notice that in its present form we
will have to do partial fractions twice.
However, if we combine the two terms up we will only be doing partial
fractions once. Not only that, but the
denominator for the combined term will be identical to the denominator of the
first term. This means that we are going
to partial fraction up a term with that denominator no matter what so we might
as well make the numerator slightly messier and then just partial fraction
once.
This is one of those things where we are apparently making
the problem messier, but in the process we are going to save ourselves a fair
amount of work!
Combining the two terms gives,
The partial fraction decomposition for this transform is,
Setting numerators equal gives,
Picking appropriate values of s and solving for the constants gives,
Plugging in the constants gives,
Finally taking the inverse transform gives us the solution
to the IVP.

That was a fair amount of work for a problem that probably
could have been solved much quicker using the techniques from the previous
chapter. The point of this problem
however, was to show how we would use Laplace
transforms to solve an IVP.
There are a couple of things to note here about using Laplace transforms to solve an IVP. First, using Laplace
transforms reduces a differential equation down to an algebra problem. In the case of the last example the algebra
was probably more complicated than the straight forward approach from the last
chapter. However, in later problems this
will be reversed. The algebra, while
still very messy, will often be easier than a straight forward approach.
Second, unlike the approach in the last chapter, we did not
need to first find a general solution, differentiate this, plug in the initial
conditions and then solve for the constants to get the solution. With Laplace transforms, the initial
conditions are applied during the first step and at the end we get the actual
solution instead of a general solution.
In many of the later problems Laplace
transforms will make the problems significantly easier to work than if we had
done the straight forward approach of the last chapter. Also, as we will see, there are some differential
equations that simply can’t be done using the techniques from the last chapter
and so, in those cases, Laplace transforms
will be our only solution.
Let’s take a look at another fairly simple problem.
Example 2 Solve
the following IVP.
Solution
As with the first example, let’s first take the Laplace transform of all the terms in the differential
equation. We’ll the plug in the
initial conditions to get,
Now solve for Y(s).
Now, as we did in the last example we’ll go ahead and
combine the two terms together as we will have to partial fraction up the
first denominator anyway, so we may as well make the numerator a little more
complex and just do a single partial fraction. This will give,
The partial fraction decomposition is then,
Setting numerator equal gives,
In this case it’s probably easier to just set coefficients
equal and solve the resulting system of equation rather than pick values of s.
So, here is the system and its solution.
We will get a common denominator of 125 on all these
coefficients and factor that out when we go to plug them back into the transform. Doing this gives,
Notice that we also had to factor a 2 out of the
denominator of the first term and fix up the numerator of the last term in
order to get them to match up to the correct entries in our table of transforms.
Taking the inverse transform then gives,

Example 3 Solve
the following IVP.
Solution
Take the Laplace
transform of everything and plug in the initial conditions.
Now solve for Y(s)
and combine into a single term as we did in the previous two examples.
Now, do the partial fractions on this. First let’s get the partial fraction
decomposition.
Now, setting numerators equal gives,
Setting coefficients equal and solving for the constants
gives,
Now, plug these into the decomposition, complete the
square on the denominator of the second term and then fix up the numerators
for the inverse transform process.
Finally, take the inverse transform.

To this point we’ve only looked at IVP’s in which the
initial values were at t = 0. This is because we need the initial values to
be at this point in order to take the Laplace
transform of the derivatives. The
problem with all of this is that there are IVP’s out there in the world that
have initial values at places other than t
= 0. Laplace
transforms would not be as useful as it is if we couldn’t use it on these types
of IVP’s. So, we need to take a look at
an example in which the initial conditions are not at t = 0 in order to see how to handle these kinds of problems.
Example 4 Solve
the following IVP.
Solution
The first thing that we will need to do here is to take
care of the fact that initial conditions are not at t = 0. The only way that
we can take the Laplace transform of the
derivatives is to have the initial conditions at t = 0.
This means that we will need to formulate the IVP in such
a way that the initial conditions are at t
= 0. This is actually fairly
simple to do, however we will need to do a change of variable to make it
work. We are going to define
Let’s start with the original differential equation.
Notice that we put in the (t) part on the derivatives to make sure that we get things
correct here. We will next substitute
in for t.
Now, to simplify life a little let’s define,
Then, by the chain rule, we get the following for the
first derivative.
By a similar argument we get the following for the second
derivative.
The initial conditions for u(η) are,
The IVP under these new variables is then,
This is an IVP that we can use Laplace
transforms on provided we replace all the t’s
in our table with η’s. So, taking the Laplace
transform of this new differential equation and plugging in the new initial
conditions gives,
Solving for U(s)
gives,
Note that unlike the previous examples we did not
completely combine all the terms this time.
In all the previous examples we did this because the denominator of
one of the terms was the common denominator for all the terms. Therefore, upon combining, all we did was
make the numerator a little messier, and reduced the number of partial
fractions required down from two to one.
Note that all the terms in this transform that had only powers of s in the denominator were combined for
exactly this reason.
In this transform however, if we combined both of the
remaining terms into a single term we would be left with a fairly involved
partial fraction problem. Therefore,
in this case, it would probably be easier to just do partial fractions
twice. We’ve done several partial
fractions problems in this section and many partial fraction problems in the
previous couple of sections so we’re going to leave the details of the
partial fractioning to you to check.
Partial fractioning each of the terms in our transform gives us the following.
Plugging these into our transform and combining like terms
gives us
Now, taking the inverse transform will give the solution
to our new IVP. Don’t forget to use η’s
instead of t’s!
This is not the solution that we are after of course. We are after y(t). However, we can get
this by noticing that
So the solution to the original IVP is,

So, we can now do IVP’s that don’t have initial conditions
that are at t = 0. We also saw in the last example that it isn’t
always the best to combine all the terms into a single partial fraction problem
as we have been doing prior to this example.
The examples worked in this section would have been just as easy,
if not easier, if we had used techniques from the previous chapter. They were worked here using Laplace
transforms to illustrate the technique and method.
0 comments:
Post a Comment