Differential Equations
Also, depending upon the E-Book, it will be possible to download the complete E-Book, the chapter containing the current section and/or the current section. You can do this be clicking on the E-Book, Chapter, and/or the Section link provided below.
For those pages with mathematics on them you can, in most cases, enlarge the mathematics portion by clicking on the equation. Click the enlarged version to hide it.
Step Functions
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
Before proceeding into differential equations we will need one more formula. We will need to know how to take the
Fact
| Suppose that f, f’, f”,…f(n-1) 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
| | |
|
Notice that the two function evaluations that appear in these formulas,
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
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
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 for the previous chapter. The point of this problem however, was to show how we would use
There are a couple of things to note here about using
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
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
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
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
| Example 4 Solve the following IVP. Solution The first thing that we will need to do here it to take care of the fact that initial conditions are not at t = 0. The only way that we can take the
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 that 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
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 denominator where 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 just as easy, if not easier, if we had used techniques from the previous chapter. They were worked here using
No comments:
Post a Comment