Differential Equations Project

Computer Lab Modules

 

Laplace Transforms Solution of D.E.'s Involving Square Wave Pulses

 

 

Background. Laplace transforms are useful in solving initial value problems for ordinary differential equations and systems of ordinary differential equations. Step functions are used in various applications involving electrical circuits and mechanical vibrations.
In this exercise we explore how to make pulses from step functions and the use the method of using Laplace transforms to find the solution to the D. E. First load Mathematica's built in "LaplaceTransform" subroutine package.

[Graphics:e19.txtgr2.gif][Graphics:e19.txtgr1.gif]
 

Computer Lab Work.

Exercise 1. Use Laplace transforms to solve the initial value problem
[Graphics:e19.txtgr3.gif]
where f[t] is consists of 6 pulses given by the following construction:

[Graphics:e19.txtgr2.gif][Graphics:e19.txtgr4.gif]

First set up the initial conditions and find the Laplace transform of f(t).

[Graphics:e19.txtgr2.gif][Graphics:e19.txtgr5.gif]

Second, set up the D.E. in Laplace transformation format and solve for Y[s].

[Graphics:e19.txtgr2.gif][Graphics:e19.txtgr6.gif]

Since Y[s] is in its expanded form, we can find the inverse transform of each term in the above sum is done in the following loop. In order to rapidly find the inverse Laplace transform of Y[s] it speeds things up to find the inverse of each term in the sum and add them up as we go along. The answer is stored in the variable named "sum".

[Graphics:e19.txtgr2.gif][Graphics:e19.txtgr7.gif]

Now plot the solution over the interval 0 <= t <= 21.

[Graphics:e19.txtgr2.gif][Graphics:e19.txtgr8.gif]


Exercise 2. Solve the initial value problem
[Graphics:e19.txtgr9.gif]
where f[t] is consists of 2 pulses given by the following construction:

[Graphics:e19.txtgr2.gif][Graphics:e19.txtgr10.gif]

First set up the initial conditions and find the Laplace transform of f(t).

[Graphics:e19.txtgr2.gif][Graphics:e19.txtgr11.gif]

Second, set up the D.E. in Laplace transformation format and solve for Y[s].

[Graphics:e19.txtgr2.gif][Graphics:e19.txtgr12.gif]

Since Y[s] is in its expanded form, we can find the inverse transform of each term in the above sum is done in the following loop. In order to rapidly find the inverse Laplace transform of Y[s] it speeds things up to find the inverse of each term in the sum and add them up as we go along. The answer is stored in the variable named "sum".

[Graphics:e19.txtgr2.gif][Graphics:e19.txtgr13.gif]

Now plot the solution over the interval 0 <= t <= 12.

[Graphics:e19.txtgr2.gif][Graphics:e19.txtgr14.gif]
 

Solutions.

 

 

Return to the Differential Equations Project

Return to the Numerical Analysis Project

Return to the Complex Analysis Project

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews, 1998