Differential equations and Laplace transforms: Laplace-transformations
Laplace transforms of differential equations
The Laplace transform of the derivative #f'# of a function #f# can be expressed in the Laplace transform of #f# without the use of derivatives:
Derivative in the time domain If #f# is a differentiable function on #\ivco{0}{\infty}#, the Laplace transform #\laplace{f}# exists, and #\lim_{t\to\infty}\ee^{-st} f(t) = 0#, then \[\mathcal{L} f' (s) = s\cdot (\mathcal{ L}f)(s)-f(0)\]
More generally, if #n# is a natural number, #f# is a piecewise #n#-fold differentiable function, and, for all #k# with #0\le k\lt n#, the limit#\lim_{t\to\infty}\ee^{-st} f^{(k)}(t) = 0# exists, then \[ \laplace{\left(f^{(n)}\right)}(s) = s^n\cdot \laplace{f}(s)-s^{n-1}f(0)-s^{n-2}f'(0)-\cdots-s\cdot f^{(n-2)}(0)-f^{(n-1)}(0)\]
Thanks to this property we can convert a linear differential equation with constant coefficients using Laplace transform into an algebraic equation. By calculating the inverse Laplace transform we can solve the differential equation. For a second order ODE are the operations shown below in a diagram
\[\begin{array}{lcr}f''(t)+p\cdot f'(t)+q = r(t)&{\laplace{}\atop\longrightarrow}& P(s)\cdot\laplace{f}(s) = \laplace{r}(s)\\ \text{solution of }\big\uparrow&&\big\downarrow\text{ solve}\\ f(t) =\mathcal{L}^{-1}(F)(t)&{\mathcal{L}^{-1}\atop \longleftarrow}& \laplace{f}(s) = F(s)\end{array}\]
Below are examples of this solution method.
#x(t) =# # 2\cdot \euler^{3 t } + 5 t \cdot \euler^{3 t } + {{t^{14}\cdot \euler^{3 t }}\over{182}} #
In order to find the solution, we write #y=\mathcal{L}(x)#. Then
\[\begin{array}{rcl}
\mathcal{L}(x') (s)&=& s\cdot y(s)-x(0) =s\cdot y(s)-2\\
\mathcal{L}(x'') (s)&=& s^2\cdot y(s)-2s =s^2\cdot y(s)-2\cdot s-11 \\
\mathcal{L}(-t^{12}\cdot \euler^{3\cdot t}) (s)&=&\displaystyle -{{479001600}\over{\left(s-3\right)^{13}}} \\
\end{array}\]
Therefore, after all terms are moved to the left, the Laplace transform applied to the differential equation \( {{d^2}\over{d t^2}} x-6 {{d}\over{d t}} x+9 x=t^{12} \euler^{3 t}\) gives
\[-6\cdot \left(s\cdot y-2\right)+s^2\cdot y+9\cdot y-2\cdot s-{{479001600}\over{\left(s-3\right)^{13}}}-11=0\]
This can be rewritten to
\[\left(s^2-6\cdot s+9\right)\cdot y-2\cdot s-{{479001600}\over{\left(s-3\right)^{13}}}+1=0\]
Solving this equation with unknown #y# gives
\[ y (s)= {{2\cdot s^{14}-79\cdot s^{13}+1443\cdot s^{12}-16146\cdot s^{11}+123552\cdot s^{10}-683397\cdot s^9+2814669\cdot s^8-8756748\cdot s^7+20640906\cdot s^6-36590697\cdot s^5+47849373\cdot s^4-44522946\cdot s^3+27634932\cdot s^2-10097379\cdot s+480595923}\over{s^{15}-45\cdot s^{14}+945\cdot s^{13}-12285\cdot s^{12}+110565\cdot s^{11}-729729\cdot s^{10}+3648645\cdot s^9-14073345\cdot s^8+42220035\cdot s^7-98513415\cdot s^6+177324147\cdot s^5-241805655\cdot s^4+241805655\cdot s^3-167403915\cdot s^2+71744535\cdot s-14348907}} \]
Partial fraction decomposition of the right-hand side leads to
\[ \begin{array}{rcl}y(s)
&=&\displaystyle {{2}\over{s-3}}+{{5}\over{\left(s-3\right)^2}}+{{479001600}\over{\left(s-3\right)^{15}}}\\
\end{array}\]
so linearity of the inverse Laplace transform and determination of the inverse Laplace tansforms of the terms gives:
\[ x (t)= 2\cdot \euler^{3 t } + 5 t \cdot \euler^{3 t } + {{t^{14}\cdot \euler^{3 t }}\over{182}} \]
Or visit omptest.org if jou are taking an OMPT exam.