13.7. Solving Differential Equations Using Convolution

Consider a second-order linear ordinary differential equation with initial conditions:

y′′+ay′+by =f(t)

This initial setup presents a second-order linear ordinary differential equation (ODE) where 'y' is the unknown function, 'a' and 'b' are constants, and 'f(t)' is a function representing external inputs (a forcing function). The prime notation indicates derivatives of 'y' with respect to time. This equation models many physical systems and phenomena, including mechanical vibrations and electrical circuit responses.

Think of a swing in a playground. The forces acting on it (like someone pushing the swing) can be modeled as 'f(t)', the initial position and speed of the swing represent the initial conditions, and the formulas help predict how high and how far the swing will go over time.

Taking Laplace Transform and applying initial conditions:

s²Y(s)+asY(s)+bY(s)=F(s)+initial terms

The Laplace Transform converts the differential equation into an algebraic equation in the 's' domain, which is often easier to work with. Here, 'Y(s)' is the Laplace Transform of 'y(t)', and 'F(s)' represents the Laplace Transform of the forcing function 'f(t)'. The initial conditions provide the necessary values to solve for 'Y(s)' completely.

Imagine translating a recipe (the differential equation) into a simpler form to understand the steps (algebraic equation) without having to worry about timing and cooking methods. The Laplace Transform simplifies the process, making it easier to arrive at a solution.

Solving for Y(s):

Y(s)=F(s)+terms from initial conditions

(s²+as+b)

In rearranging the transformed equation, Y(s) is isolated on one side. This equation expresses Y(s) as a combination of the transformed input function F(s) and terms derived from initial conditions, divided by a polynomial related to the system characteristics (represented as (s² + as + b)). This step is crucial for finding the solution in the time domain later.

Consider a problem where you need to calculate the final height of a structure but first need to consider its foundation and materials. Here, the denominator (s² + as + b) represents the foundation, while F(s) includes the external forces acting upon it, showing how all factors influence the final outcome.

Now, using convolution:

y(t)=(f ∗h)(t) Where h(t) is the inverse Laplace Transform of 1/(s²+as+b) which acts as the impulse response.

The solution y(t) is expressed as the convolution of the input function f(t) and the impulse response h(t). The impulse response is derived from the inverse Laplace Transform of the polynomial in the denominator, allowing for the effect of the initial conditions and the external forcing function to be combined to provide the overall response of the system over time. This convolution encapsulates how the system reacts to the inputs throughout its time evolution.

Imagine playing with a trampoline. When you jump (f(t)), the trampoline's response (h(t)) depends on how springy and bouncy it is. The overall bounce you feel as you jump at different times is like the convolution, combining your action and the trampoline's reaction to give a complete view of your bouncing experience.

--