Fourier transformation and differential equation 2nd order

152 Views Asked by At

Applying Fourier Transfomation in the following differential:

$$ m\ddot{x}+c\dot{x}+kx =f(t)\tag1 $$ we arrive at the general form for the particular solution:

$$x(t)=f(t)* \int_{-\infty}^\infty{\dfrac{\mathrm{e}^{2\mathrm{i}{\pi}t{\xi}}}{-4{\pi}^2m{\xi}^2+2\mathrm{i}{\pi}c{\xi}+k}}dξ$$

Setting the right term of the convolution in an online integral calculator (https://www.integral-calculator.com/) we get the following solution of ANTIDERIVATIVE. $$A(ξ)=\dfrac{\mathrm{e}^{-\frac{\left(\mathrm{i}\sqrt{4km-c^2}+c\right)t}{2m}}\left(\operatorname{Ei}\left(\frac{t\left(\mathrm{i}\left(4{\pi}m{\xi}+\sqrt{4km-c^2}\right)+c\right)}{2m}\right)-\mathrm{e}^\frac{\mathrm{i}\sqrt{4km-c^2}t}{m}\operatorname{Ei}\left(\frac{t\left(\mathrm{i}\left(4{\pi}m{\xi}-\sqrt{4km-c^2}\right)+c\right)}{2m}\right)\right)}{2{\pi}\sqrt{4km-c^2}}+C$$

where $Ei(ξ)$ the function of exponential integral.

Computing the $A(+\infty)$ and the $A(-\infty)$ the calculator concludes:

$$\int_{-\infty}^\infty{\dfrac{\mathrm{e}^{2\mathrm{i}{\pi}t{\xi}}}{-4{\pi}^2m{\xi}^2+2\mathrm{i}{\pi}c{\xi}+k}}dξ=\dfrac{2\mathrm{e}^{-\frac{ct}{2m}}\sin\left(\frac{\sqrt{4km-c^2}t}{2m}\right)}{\sqrt{4km-c^2}}$$

My questions are: 1. What are the limit values of $Ei(ξ)$ which have set the program and why? 2. Is it possible to define the exponential integral from $-\infty$ to $+\infty$ ?

1

There are 1 best solutions below

3
On

I've another way of solving this DE.

The general form of a second-order inhomogeneous linear differential equation is: \begin{equation} \text{c}_\text{2}\cdot\frac{\text{d}^2\text{f}\left(t\right)}{\text{d}t^2}+\text{c}_\text{1}\cdot\frac{\text{d}\text{f}\left(t\right)}{\text{d}t}+\text{c}_\text{0}\cdot\text{f}\left(t\right)=\text{g}\left(t\right) \tag1\end{equation}

If the Laplace transformation is taken from both sides, one obtains on the left: \begin{align} \mathcal{L}_t\left[\text{c}_\text{2}\cdot\frac{\text{d}^2\text{f}\left(t\right)}{\text{d}t^2}+\text{c}_\text{1}\cdot\frac{\text{d}\text{f}\left(t\right)}{\text{d}t}+\text{c}_\text{0}\cdot\text{f}\left(t\right)\right]_{\left(\text{s}\right)} &= \\ \mathcal{L}_t\left[\text{c}_\text{2}\cdot\frac{\text{d}^2\text{f}\left(t\right)}{\text{d}t^2}\right]_{\left(\text{s}\right)}+\mathcal{L}_t\left[\text{c}_\text{1}\cdot\frac{\text{d}\text{f}\left(t\right)}{\text{d}t}\right]_{\left(\text{s}\right)}+\mathcal{L}_t\left[\text{c}_\text{0}\cdot\text{f}\left(t\right)\right]_{\left(\text{s}\right)} &= \notag\\ \text{c}_\text{2}\cdot\mathcal{L}_t\left[\frac{\text{d}^2\text{f}\left(t\right)}{\text{d}t^2}\right]_{\left(\text{s}\right)}+\text{c}_\text{1}\cdot\mathcal{L}_t\left[\frac{\text{d}\text{f}\left(t\right)}{\text{d}t}\right]_{\left(\text{s}\right)}+\text{c}_\text{0}\cdot\mathcal{L}_t\left[\text{f}\left(t\right)\right]_{\left(\text{s}\right)} &= \notag\\ \text{c}_\text{2}\cdot\left(\text{s}^2\cdot\text{F}\left(\text{s}\right)-\text{s}\cdot\text{f}\left(0\right)-\text{f}'\left(0\right)\right)+\text{c}_\text{1}\cdot\left(\text{s}\cdot\text{F}\left(\text{s}\right)-\text{f}\left(0\right)\right)+\text{c}_\text{0}\cdot\text{F}\left(\text{s}\right) \tag2 \end{align}

And on the right: \begin{equation} \mathcal{L}_t\left[\text{g}\left(t\right)\right]_{\left(\text{s}\right)}=\text{G}\left(\text{s}\right) \tag3\end{equation}

Now the left and right sides can be set equal: \begin{equation} \text{c}_\text{2}\cdot\left(\text{s}^2\cdot\text{F}\left(\text{s}\right)-\text{s}\cdot\text{f}\left(0\right)-\text{f}'\left(0\right)\right)+\text{c}_\text{1}\cdot\left(\text{s}\cdot\text{F}\left(\text{s}\right)-\text{f}\left(0\right)\right)+\text{c}_\text{0}\cdot\text{F}\left(\text{s}\right)=\text{G}\left(\text{s}\right) \tag4\end{equation} From the equivalent equations, $\text{F}\left(\text{s}\right)$ can be solved: \begin{align} \text{F}\left(\text{s}\right) & = \frac{\text{G}\left(\text{s}\right)+\text{f}\left(0\right)\cdot\left(\text{s}\cdot\text{c}_\text{2}+\text{c}_\text{1}\right)+\text{c}_\text{2}\cdot\text{f}'\left(0\right)}{\text{s}\cdot\left(\text{s}\cdot\text{c}_\text{2}+\text{c}_\text{1}\right)+\text{c}_\text{0}} \\ \notag & = \frac{\text{G}\left(\text{s}\right)}{\text{s}\cdot\left(\text{s}\cdot\text{c}_\text{2}+\text{c}_\text{1}\right)+\text{c}_\text{0}}+\frac{\text{f}\left(0\right)\cdot\left(\text{s}\cdot\text{c}_\text{2}+\text{c}_\text{1}\right)}{\text{s}\cdot\left(\text{s}\cdot\text{c}_\text{2}+\text{c}_\text{1}\right)+\text{c}_\text{0}}+\frac{\text{c}_\text{2}\cdot\text{f}'\left(0\right)}{\text{s}\cdot\left(\text{s}\cdot\text{c}_\text{2}+\text{c}_\text{1}\right)+\text{c}_\text{0}} \\ \notag & = \frac{\text{G}\left(\text{s}\right)}{\text{c}_\text{2}\cdot\text{s}^2+\text{c}_\text{1}\cdot\text{s}+\text{c}_\text{0}}+\frac{\text{f}\left(0\right)\cdot\left(\text{s}\cdot\text{c}_\text{2}+\text{c}_\text{1}\right)}{\text{c}_\text{2}\cdot\text{s}^2+\text{c}_\text{1}\cdot\text{s}+\text{c}_\text{0}}+\frac{\text{c}_\text{2}\cdot\text{f}'\left(0\right)}{\text{c}_\text{2}\cdot\text{s}^2+\text{c}_\text{1}\cdot\text{s}+\text{c}_\text{0}} \\ \notag & = \frac{1}{\text{c}_\text{2}}\cdot\frac{\text{G}\left(\text{s}\right)}{\text{s}^2+\frac{\text{c}_\text{1}}{\text{c}_\text{2}}\cdot\text{s}+\frac{\text{c}_\text{0}}{\text{c}_\text{2}}}+\frac{1}{\text{c}_\text{2}}\cdot\frac{\text{f}\left(0\right)\cdot\left(\text{s}\cdot\text{c}_\text{2}+\text{c}_\text{1}\right)}{\text{s}^2+\frac{\text{c}_\text{1}}{\text{c}_\text{2}}\cdot\text{s}+\frac{\text{c}_\text{0}}{\text{c}_\text{2}}}+\frac{1}{\text{c}_\text{2}}\cdot\frac{\text{c}_\text{2}\cdot\text{f}'\left(0\right)}{\text{s}^2+\frac{\text{c}_\text{1}}{\text{c}_\text{2}}\cdot\text{s}+\frac{\text{c}_\text{0}}{\text{c}_\text{2}}} \\ \notag & = \frac{1}{\text{c}_\text{2}}\cdot\frac{\text{G}\left(\text{s}\right)}{\text{s}^2+\frac{\text{c}_\text{1}}{\text{c}_\text{2}}\cdot\text{s}+\frac{\text{c}_\text{0}}{\text{c}_\text{2}}}+\frac{\text{f}\left(0\right)}{\text{c}_\text{2}}\cdot\frac{\text{s}\cdot\text{c}_\text{2}+\text{c}_\text{1}}{\text{s}^2+\frac{\text{c}_\text{1}}{\text{c}_\text{2}}\cdot\text{s}+\frac{\text{c}_\text{0}}{\text{c}_\text{2}}}+\text{f}'\left(0\right)\cdot\frac{1}{\text{s}^2+\frac{\text{c}_\text{1}}{\text{c}_\text{2}}\cdot\text{s}+\frac{\text{c}_\text{0}}{\text{c}_\text{2}}} \tag5 \end{align} The denominator formula above are all the same and can be written via sperating the squares as follows: \begin{equation} \text{s}^2+\frac{\text{c}_\text{1}}{\text{c}_\text{2}}\cdot\text{s}+\frac{\text{c}_\text{0}}{\text{c}_\text{2}}=\left(\text{s}-\left(-\frac{1}{2}\cdot\frac{\text{c}_\text{1}}{\text{c}_\text{2}}\right)\right)^2-\left(\sqrt{\frac{1}{4}\cdot\left(\frac{\text{c}_\text{1}}{\text{c}_\text{2}}\right)^2-\frac{\text{c}_\text{0}}{\text{c}_\text{2}}}\right)^2 \tag6\end{equation}

To write it in a clearer way, the following applies: \begin{equation} \alpha:=-\frac{1}{2}\cdot\frac{\text{c}_\text{1}}{\text{c}_\text{2}}\;\;\;\wedge\;\;\;\beta:=\sqrt{\frac{1}{4}\cdot\left(\frac{\text{c}_\text{1}}{\text{c}_\text{2}}\right)^2-\frac{\text{c}_\text{0}}{\text{c}_\text{2}}} \tag7\end{equation}

This also applies to: \begin{equation} \text{s}^2+\frac{\text{c}_\text{1}}{\text{c}_\text{2}}\cdot\text{s}+\frac{\text{c}_\text{0}}{\text{c}_\text{2}}=\left(\text{s}-\alpha\right)^2-\beta^2 \tag8\end{equation}

It can therefore be rewritten as follows: \begin{align} \text{f}\left(t\right) & = \frac{1}{\text{c}_\text{2}}\cdot\frac{\text{G}\left(\text{s}\right)}{\left(\text{s}-\alpha\right)^2-\beta^2}+\frac{\text{f}\left(0\right)}{\text{c}_\text{2}}\cdot\frac{\text{s}\cdot\text{c}_\text{2}+\text{c}_\text{1}}{\left(\text{s}-\alpha\right)^2-\beta^2}+\frac{\text{f}'\left(0\right)}{\left(\text{s}-\alpha\right)^2-\beta^2} \\ \notag & = \frac{1}{\text{c}_\text{2}\cdot\beta}\cdot\frac{\text{G}\left(\text{s}\right)\cdot\beta}{\left(\text{s}-\alpha\right)^2-\beta^2}+\frac{\text{f}\left(0\right)}{\text{c}_\text{2}}\cdot\frac{\text{s}\cdot\text{c}_\text{2}+\text{c}_\text{1}}{\left(\text{s}-\alpha\right)^2-\beta^2}+\frac{\text{f}'\left(0\right)}{\beta}\cdot\frac{\beta}{\left(\text{s}-\alpha\right)^2-\beta^2} \tag9 \end{align}

In order to finally solve the differential equation, the inverse Laplace transformation must be looked at. Using well-known Laplace transformation tables, it can be said that: \begin{equation} \text{f}_\text{1}\left(t\right):=\mathcal{L}_\text{s}^{-1}\left[\frac{\beta}{\left(\text{s}-\alpha\right)^2-\beta^2}\right]_{\left(t\right)}=\exp\left(\alpha\cdot t\right)\cdot\sinh\left(\beta\cdot t\right) \tag{10}\end{equation}

And with the help of convolution theory also applies: \begin{align} \text{f}_\text{2}\left(t\right) & := \mathcal{L}_\text{s}^{-1}\left[\frac{\text{G}\left(\text{s}\right)\cdot\beta}{\left(\text{s}-\alpha\right)^2-\beta^2}\right]_{\left(t\right)} \\ \notag & = \int_0^t\mathcal{L}_\text{s}^{-1}\left[\text{G}\left(\text{s}\right)\right]_{\left(x\right)}\cdot\mathcal{L}_\text{s}^{-1}\left[\frac{\beta}{\left(\text{s}-\alpha\right)^2-\beta^2}\right]_{\left(t-x\right)}\;\text{d}x \\ \notag & = \int_0^t\text{g}\left(x\right)\cdot\exp\left(\alpha\cdot\left(t-x\right)\right)\cdot\sinh\left(\beta\cdot\left(t-x\right)\right)\;\text{d}x \tag{11} \end{align}

Now the inverse Laplace transformation can be taken, so that one obtains the following: \begin{align} \text{f}\left(t\right) & = \mathcal{L}_\text{s}^{-1}\left[\text{F}\left(\text{s}\right)\right]_{\left(t\right)} \\ \notag & = \frac{\text{f}_\text{2}\left(t\right)}{\text{c}_\text{2}\cdot\beta}+\frac{\text{f}\left(0\right)}{\text{c}_\text{2}}\cdot\left(\text{c}_\text{2}\cdot\frac{\text{s}}{\left(\text{s}-\alpha\right)^2-\beta^2}+\frac{\text{c}_\text{1}}{\beta}\cdot\text{f}_\text{1}\left(t\right)\right)+\frac{\text{f}'\left(0\right)}{\beta}\cdot\text{f}_\text{1}\left(t\right) \\ \notag & = \frac{\text{f}_\text{2}\left(t\right)}{\text{c}_\text{2}\cdot\beta}+\frac{\text{f}\left(0\right)}{\text{c}_\text{2}}\cdot\left(\frac{\text{c}_\text{2}}{\beta}\cdot\frac{\text{s}\cdot\beta}{\left(\text{s}-\alpha\right)^2-\beta^2}+\frac{\text{c}_\text{1}}{\beta}\cdot\text{f}_\text{1}\left(t\right)\right)+\frac{\text{f}'\left(0\right)}{\beta}\cdot\text{f}_\text{1}\left(t\right) \\ \notag & = \frac{\text{f}_\text{2}\left(t\right)}{\text{c}_\text{2}\cdot\beta}+\frac{\text{f}\left(0\right)}{\text{c}_\text{2}}\cdot\left(\frac{\text{c}_\text{2}}{\beta}\cdot\text{f}_\text{1}'\left(t\right)+\frac{\text{c}_\text{1}}{\beta}\cdot\text{f}_\text{1}\left(t\right)\right)+\frac{\text{f}'\left(0\right)}{\beta}\cdot\text{f}_\text{1}\left(t\right) \\ \notag & = \frac{1}{\beta}\cdot\left\{\frac{\text{f}_\text{2}\left(t\right)+\text{f}\left(0\right)\cdot\left(\text{c}_\text{2}\cdot\text{f}_\text{1}'\left(t\right)+\text{c}_\text{1}\cdot\text{f}_\text{1}\left(t\right)\right)}{\text{c}_\text{2}}+\text{f}'\left(0\right)\cdot\text{f}_\text{1}\left(t\right)\right\} \tag{12} \end{align}

This provides the solution of the differential equation: \begin{equation} \text{f}\left(t\right)=\frac{1}{\beta}\cdot\left\{\frac{\text{f}_\text{2}\left(t\right)+\text{f}\left(0\right)\cdot\left(\text{c}_\text{2}\cdot\text{f}_\text{1}'\left(t\right)+\text{c}_\text{1}\cdot\text{f}_\text{1}\left(t\right)\right)}{\text{c}_\text{2}}+\text{f}'\left(0\right)\cdot\text{f}_\text{1}\left(t\right)\right\} \tag{13}\end{equation}