5. Solutions by transform¶
Faced with an initial-value problem, we can take the Laplace transform of the ODE and solve for \(X(s)\), the transform of the solution. This part is fairly routine. Most of the effort comes in finding the inverse transform \(x(t)\) corresponding to \(X(s)\), for which there is no easy formula.
The table below summarizes the transforms we have learned so far. (We will add to it in the next chapter.) When faced with finding an inverse transform, the game is to manipulate the expression algebraically until it can be matched with one or more lines in the right column.
Function |
Transform |
---|---|
\(x(t)\) |
\(X(s)\) |
\(x'(t)\) |
\(sX(s)-x(0)\) |
1 |
\(\dfrac{1}{s}\) |
\(e^{at}\) |
\(\dfrac{1}{s-a}\) |
\(H(t-T)\) |
\(\dfrac{e^{-sT}}{s}\) |
\(H(t-T)f(t-T)\) |
\(e^{-sT}F(s)\) (shift theorem) |
\(\delta(t-T)\) |
\(e^{-sT}\) |
One major tool at our disposal is the partial fraction decomposition you probably learned for integrating rational functions.
Example
Solve \(x'-ax=e^{c t}\), where \(a\neq c\) are constants.
Solution
Transforming both sides, we get
and thus
We can immediately recognize the first term as the transform of \(e^{at}x_0\). To invert the transform of the second term, we need to express it using partial fractions. There is a general procedure for this, but it’s worth knowing a simple mnemonic for this particular form:
Thus the inverse transform of this fraction is a sum of two exponentials,
This is a particular solution that has \(x_p(0)=0\). The general solution is \(e^{at}x_0 + x_p(t)\).
We saw before that both steps and impulses in the forcing term lead to step functions in the solution. The next example shows how that fact fits with the shift theorem.
Example
Solve \(x'=ax + \delta(t-T)\), with \(x(0)=0\).
Solution
Transform the ODE to get
thus
In a previous example we showed that the inverse transform of this is
5.1. Transfer function¶
Perhaps it’s clear by now that
always leads to
as a particular solution. The term \(1/(s-a)\) has a special name.
For the first-order operator \(\opA[x]=x'-ax\) with constant \(a\), the function \(1/(s-a)\) is called the transfer function of the operator. The value \(s=a\) is called a pole of the transfer function.
The transfer function is a multiplier that “transfers” \(F(s)\), the transform of the forcing function, to \(X(s)\), the transform of a particular solution. The pole of the transfer function corresponds to an exponential solution in the original \(t\) value.
The impulse response is the case when \(f(t)=\delta(t)\), so that \(F(s)=1\). This implies the following.
The transfer function of a linear operator is the Laplace transform of the impulse response of the operator.