{ "cells": [ { "cell_type": "markdown", "metadata": { "nbsphinx": "hidden" }, "source": [ "# The Laplace Transform\n", "\n", "*This Jupyter notebook is part of a [collection of notebooks](../index.ipynb) in the bachelors module Signals and Systems, Communications Engineering, Universität Rostock. Please direct questions and suggestions to [Sascha.Spors@uni-rostock.de](mailto:Sascha.Spors@uni-rostock.de).*" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Properties\n", "\n", "The Laplace transform has a number of specific properties which can be concluded directly from its definition. The most important ones in the context of signals and systems are reviewed in the following." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Linearity\n", "\n", "The Laplace transform is a linear operation, as can be concluded straightforwardly from its definition, . For two signals $x_1(t)$ and $x_2(t)$ with Laplace transforms $X_1(s) = \\mathcal{L} \\{ x_1(t) \\}$ and $X_2(s) = \\mathcal{L} \\{ x_2(t) \\}$ the following holds\n", "\n", "\\begin{equation}\n", "\\mathcal{L} \\{ A \\cdot x_1(t) + B \\cdot x_2(t) \\} = A \\cdot X_1(s) + B \\cdot X_2(s)\n", "\\end{equation}\n", "\n", "with $A, B \\in \\mathbb{C}$. The region of convergence (ROC) of the superposition $A x_1(t) + B x_2(t)$ is as least as large as the intersection of the ROCs of $x_1(t)$ and $x_2(t)$\n", "\n", "\\begin{equation}\n", "\\text{ROC} \\{ A \\cdot x_1(t) + B \\cdot x_2(t) \\} \\supseteq \\text{ROC} \\{ x_1(t) \\} \\cap \\text{ROC} \\{ x_2(t) \\}\n", "\\end{equation}\n", "\n", "The Laplace transform of a weighted superposition of signals is equal to the weighted superposition of the individual Laplace transforms. This property is useful to derive the Laplace transform of signals that can be expressed as superposition of other signals for which the Laplace transforms are known or can be calculated easier." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "#### Transformation of the cosine and sine signal\n", "\n", "The Laplace transform of the causal cosine $\\epsilon(t) \\cdot \\cos(\\omega_0 t)$ and sine $\\epsilon(t) \\cdot \\sin(\\omega_0 t)$ signal is derived by expressing the cosine/sine as harmonic exponential signals using [Euler's formula](https://en.wikipedia.org/wiki/Euler's_formula)\n", "\n", "\\begin{align}\n", "\\cos(\\omega_0 t) &= \\frac{1}{2} \\left( e^{j \\omega_0 t} + e^{-j \\omega_0 t} \\right) \\\\\n", "\\sin(\\omega_0 t) &= \\frac{1}{2j} \\left( e^{j \\omega_0 t} - e^{-j \\omega_0 t} \\right)\n", "\\end{align}\n", "\n", "The Laplace transform of the [causal complex exponential signal](definition.ipynb#Transformation-of-the-Causal-Complex-Exponential-Signal) reads\n", "\n", "\\begin{equation}\n", "\\mathcal{L} \\{ \\epsilon(t) \\cdot e^{-s_0 t} \\} = \\frac{1}{s + s_0} \\qquad \\text{for } \\Re \\{ s \\} > \\Re \\{ - s_0 \\}\n", "\\end{equation}\n", "\n", "Using the linearity of the Laplace transform yields\n", "\n", "\\begin{equation}\n", "\\mathcal{L} \\{ \\epsilon(t) \\cdot \\cos(\\omega_0 t) \\} = \n", "\\frac{1}{2} \\left( \\frac{1}{s - j \\omega_0} + \\frac{1}{s + j \\omega_0} \\right) = \\frac{s}{s^2 + \\omega_0^2}\n", "\\qquad \\text{for } \\Re \\{ s \\} > 0\n", "\\end{equation}\n", "\n", "and\n", "\n", "\\begin{equation}\n", "\\mathcal{L} \\{ \\epsilon(t) \\cdot \\sin(\\omega_0 t) \\} = \n", "\\frac{1}{2j} \\left( \\frac{1}{s - j \\omega_0} - \\frac{1}{s + j \\omega_0} \\right) = \\frac{\\omega_0}{s^2 + \\omega_0^2}\n", "\\qquad \\text{for } \\Re \\{ s \\} > 0\n", "\\end{equation}" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "**Example - Computational transformation of the cosine signal**\n", "\n", "The Laplace transform of the causal cosine signal is calculated in `Sympy`" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAQ8AAAAyCAYAAABCvF75AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAOs0lEQVR4Ae2d73HcNhPGzxoVoNgdxB3YcQWxO3D8VmCng2T8zd80SQdJKnCSDhJXkNgdxG8FVtSB8vwgLAySIAkeeRRP2p2hQAILYLHEPlj84ene1dXVbozevHnzh3i+V/hhjNfTXQOugePVgGz8TNL/pusb3V8OteRkKJE0FUBBvyl04BhTlqe7Bo5cAxEwvlcz3useIOmlQfCIwPFR4c+9JXiCa8A1cKs0IHvHUQgAMtSwe33TFhXwShm/Vfh4qABPcw24Bm6nBmT7P6hlXyr8ptTCouch5i/FTMZiplJBHucacA3cLg0IB/A+HinEkehQETzExQLpuTJ97OTwCNeAa+AuaQAA+UFY0Fn/6IBHRJn7Cn+8SxrytroGXANdDQgHflcsTgQzkQZ1wEOpMN34AqmE/qkhaetB6czF/tH1qJXkj64B18CyGsD7eIXN5cU2wEOJ3ykR9+Q8Z1r7XnI8V51MncaIxjwdY/J014BrYH8NyB7/VG52YBrex2mryG/1/LuYBw+HtPIc4vGZZECWXlI6W8hMrZKserbGfVL8Az2DmE6uAdfAfA3gUHDe60xXsLnkeSgC95+R/Ka9DjwJDqYligKXPJF/xAQq7sTDivAnhZyEBVT+UFjKA7uTa8A1ME0Dwc6UJe28JPBQ5GtdlzK4mz5JitdhglrzvkI2e8jCh+K1HSG8jpQvlvFUIYDo5BpwDczQgOwI+8O+0owgBw9G/GR8M+qZm7Vk7HhFb/OCc1CI92dKNyAxVhrsC6qmDQ9dA/M0wIyAjYpgowE89ICBYXxbcPOZUyWDj/d4Fc90BVIcstIQ25HhuUQXirxfSvA414BrYLIGfo05wiaFeR7hQQlb8DyYNr0TQHAwBYB4p4sj8mwVsY4BYPxf10+6b3saim4QwNEHLA1Gf5imAem+5CFOK2Ql7jVkramjhmcllexVjeS/VEauMJAbePDAeseYMe5V6cRMLNj+rYttYzroY8kFoDDXYu3jhS5Ov+ZnUWhQiQCOLbSpJNuxx/Hdkw06m22LZKQfJU/2gILizlPXEB2FzoYaoDRsM+jzNDJilETeOOkFJGTLhVE8YJEDRkpWGtu25ANs2gu+7eeUb+s3alO+9fxQ8uKNLQqGKg+d2ZY2/eCCZ8Ufrd4kfyC14blu2LI/+Glp1fGnrvAdiMJiP41i3Xgg+eb0K/rF9bRFBTE634YRGo8ljYRqF/ecWVnU2FTmKiS536uivxTa1jMGzrQNY1+EYllM/xgRuZge0jn4LYeky0UqW7kQyU+ffq3QgPHgEqguQAo9UvcmSbLN7Vccj9ipnKcnCq0zhshNtrhCKDWGF8cow+iM+8gXwS8rsm6ORfKzl87CMd8VBNL9pW54tkXi64R5fxmB0tYbRakejI26GmdtSDsyom1L6qq2+dRpI3ttnlX49G6X6Fc2Q3mUg8dRjtC51un48fpRISMARnCMBPCVpg1/KZ6zK0uNbHgXfB/ULo+Fc8DLBhY9Hh29kPyrTx9indTd1ukWFLhEvzKceAh42FamRW6hkXddBoz6oqAEe0dLTSkACVsvKlS325wB1ACaeFjrMF2V2nXoOOpmYX9rNLtfSbeXsVH3T3VjHcQit9bgOyVP5YhlgD9LN6qr78eewmq60kvez2CdyhNc48j0IIbnig/9SyEdGLcer4Z1KgztiS7i+dHdYPSRjykVz/ByBom+Ctix9sMCch+xe9h77EB5wzpPX2alW33IM1kHsW5kWN3zGWiT2XkfC/FT+tUZ4GEvuDTSDVXkaYfRgL3AITCv6Qh7SSdjATgwnskLjcrLiA8wsL1uIMAagJ3V2Skeoyb9X4UAAHzwsE4FgPysNGQALL7QfdBD5GdaWvPhJrtGxfUO5UfGmsOQ6AB59gEP1g/JuxdJRmSfmv+D8vUNBsixZL/inTQ8DypwOg4NGOAfQloWSjHQpbY3Ke+VyuMcRACUKDSDFes3YcFWYQIKxb/W1Z5OAToAGuDS61UoHQJc+wbDNb7Ypm7AZy8yneyVeV6m2n5F+4LnERBJAgeEr61b/OP/8KW2sFvAJ33cW6gZfZ2e4m304CcHFie1gREPox0awXrrVT52g8IOke4xYDwA3PcS0d8SmIh/Uv8rFZjFoadOeaqD0byxixTl5FPztpx4DwGklMZUjC+2A6AqBPSYOrXziC0Q7aL9W6Kl+1XwPPZqoBS3lLHsVf9tzSS9ctKX5pU6n8Ulo1tKD6oTA+HnJ/sMoqoq5UdGpi4Y8Ftd7BD1UV+HBsQwTnZ8LmPmRwopdw7hdbSnYwCc1ZGXzRfbtmZBvV9bouI5EBbO3CgsvYsieFn+mwgl5+L96lQNCS9Qhecv6iba53V+1gAjXsntpVNCY277NVflX7171gEwluRx6D7Ur7BkHMWSYzm/KJGFRhu1Mfo+KhktvJwlwHB/ifUjC+d3zJj1OEj0aUCsTSWdIh8gl0j1JL54T1ltPSA7edvxlAN/HzCSPkiq8xBrHtS5ZL+6ADxQAkTHtPsQcYg/UgyKZU4L2Ut6qfiD131d5Tp/aeeMNuFal0bZx4pnYSzpamY9O+XHAJ4obI/IAEqtsYo1EHK31yToV4FiXYyAJYMzNsIwvRDfviBJ+da38nJ5Jxwh/0Ak9wrQM+3Mp1u0w4CU/loiwCG1rcVA/FgbW1k+P0quxsG9zymz76r71UhNoX2Ah82f+5Q0Us7kZEaQpBzdg7IcmR3aeptcyZIZJCO6qQa8yP+vQgwdg59EyoMBsrPwXFfeqTk7kLvPyDWnHgyMDoUbznvIiXn9EoumGCiErNQXDDc+9/W5S6XTT14q5B660LPdX8f0/6WOJ4Vk4vliG7BAFkCK98NxfJ4xeHSM7seMHwPqkz+AvNI3RbRb12i/qhT644kYa19IZZmjbKy889KMQH5W4q2TWfyWQjoySudiRGLUAfCKJB50SudjpOvrYMW8WSQdkDk6dbONyXTga92b8e10P7cetiwxGtY72tc+/QKZWTf5ThfAh9wYKiAESFEX79ruX3APr+IT6dk8DnTMwiUXIMlpWMocI6Yhpf50rnimRJSBLGwZo08GM9Y+XujiTErucfXpgffaBzD075rtYLGtTqP9qkIi2h6mLRgCBJKuQbwoXuAxEYBHJ7dODeARl1zgdmOUxhpCwyjaPEPPykunTR5aH++cesjbV+4+8SoPQzR3Py+CKVE+LTI95jzpXuVgeBhx7nXRPwEEwHSnq9crUhoeHzyN96NndNpZEFY8YJEDhh6vSWm2ZQzYJOCOye3nnfjhIxxsY8y/eiC5qvpVn2DKD3BAl6f6Y+gZGh2iD/hHlbdfEgbCC+q8iAOKMbXofQGPtYRgAFMrnMi/Vj0TxZrOLn0xau9yvemeDs9FPwFEAKhe8FAaBMDz3kYBGOYROlc6coU+Khm47/tiG5Ck7ttKhhOfTqQIM1rcmVVJdTOSMDrvXTcvUhcu98FI5TNXpPMa0SEHAU/8ILStJ1m+xcO16llc8P4C8UqHprHovr0+0ymNdxbLsc7e4amNUFkA1egX2+KjLmRvD5C1VR0Dn+nzI54HhPfBnG81iooGoZl35oY5VQaMlGsVkqy1gMe0Zmx0XELmtepZQtbRMugLuhhM+C2O/yk0AH6ge94zO3M24OlxkPBQWF/pTFUGcxUSVWc+7SpwhChAbQlPp6/8LcQbeHww8OBl4IqtQnoRCMDiY3ip8Xmn0KZQq8gxtZIoZxXgiXcN4EBnq9RT0BVGfVGInx2lNjGY1BjrYF2Uo4ufZmAB96B6og4JQ11DffhgOhtUxLKJTyiOdt67urriBrcf1My/L4BnlJR30heK4gc4qCvvHKA1YEKnmUTKw7QHV3FW51B+QMFGORYSG/Io3QAvjCzxeadwqLNMaoszuwa2rgH1dz5o/FvhM/M8bPcD76N6gU8FYLg1W1IYHmXj4bzXhftJmEhl3Zi7p7ppA2AR3GGFyAqY5EARAE9pTFugAHjXt/7XNXD7NaC+jx1ju8HmA3hgNLoY9ZljVoOHeDmHMGj0SmdhEa8geBW6/0L3myHJE4BNYQCOKBheEXNlo80BngnmoWtgRQ0wqEIBIwJ4XD/vflXIIZkqkrFRUG5gO8WBSoNfKFYVXmBS2Yz8JnzOwdYddZdADFAsnTvI84d78YGmtIfTlo0FNj1vCvA6wnuEa2AdDWAXOANhqp6DB4bzSglsfdYccMHryNctEJ8dm+Bh8JBR/oViFl1/q7pK4LBT/Kw1D+U3z4gFrwBOiuM3UNttqxfWOV0Dt1MD2FqyixNro4wFwABRqkZq8eHut4n1gLd5pMrt8CmOU4Lh2DVhzn8T95KB9Q5+YgBkRQ+szuNFObkGXAPSgOzBzlKlMywJPKKGMGS8jxrDCV8ommaVB+Agf3L5Yzl4NEw5AkUh+GEVDJa1EH4XoWbRNZawXKB6ATC+mwik+86UxdI8dA3ccQ3gcTR+ArIBHjIeUIVpx+sKRbHAuM//lAVg0rQIg9UzU6WOh1Ihw1wWpikJ2ChMcuCatU+UkuTkGriTGpBNYCfYZ5qyoIgGeBAhgqHGbT8XH1u8rBVQ8OgXihEg8GraZyMALNsC1e1q9JKaJBfeFm2mLXwVWlxfWU0qr8g1sC0NmNfRsNt8wTSIK8Nh1MV48D4aSJO3RzwYfJqiWBr5dZ/mRRYfQ4CjRBeKDLsmpcRDxUlWvKd8i/ZQVXm5roGj1IBsBK+Dq7PjWPI8aCQjMqMxHsUaBHD0ActY/YAYl5NrwDWwvAaY1vOzkh0b63ge1C1GzkcwLSFjx7uAZ0/qCBDLOVPYcIlqy5ecaf2kNo/zuQZcA+MakG2xPgkWFA+O9nkeO2UIp0JjAeM1VXCoLAACACl5ND59qNChs7gG1tCAbDVMVxT2Ht3oBQ8EjBmfKmQHYinCo0GwQCqb+74fVolcHrgGXANraUA2accu0u/lluoOX9WWEvI4Fdb4cCxP2+de5eEOQXzF2vmCNaT4H9eAa2B1Dcg2WUJ4p4vfy+1bZghy/QetAYeLcO5geAAAAABJRU5ErkJggg==\n", "text/latex": [ "$\\displaystyle \\left( \\frac{\\omega_{0}}{\\omega_{0}^{2} + s^{2}}, \\ 0, \\ 2 \\left|{\\arg{\\left(\\omega_{0} \\right)}}\\right| = 0\\right)$" ], "text/plain": [ "⎛ ω₀ ⎞\n", "⎜────────, 0, 2⋅│arg(ω₀)│ = 0⎟\n", "⎜ 2 2 ⎟\n", "⎝ω₀ + s ⎠" ] }, "execution_count": 1, "metadata": {}, "output_type": "execute_result" } ], "source": [ "import sympy as sym\n", "sym.init_printing()\n", "\n", "t, w0 = sym.symbols('t omega0', real=True)\n", "s = sym.symbols('s', complex=True)\n", "\n", "x = sym.sin(w0*t)\n", "X, a, cond = sym.laplace_transform(x, t, s)\n", "X, a, cond" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "**Exercise**\n", "\n", "* Determine the zeros and poles of the Laplace transform $\\mathcal{L} \\{ \\epsilon(t) \\cdot \\sin(\\omega_0 t) \\}$ by manual calculus and by extending above example. What specific symmetries do the zeros and poles show in the $s$-plane?" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Symmetry for Real-Valued Signals\n", "\n", "The Laplace transform of the complex conjugate $x^*(t)$ of a signal $x(t)$ is given as\n", "\n", "\\begin{equation}\n", "\\mathcal{L} \\{ x^*(t) \\} = \\int_{-\\infty}^{\\infty} x^*(t) \\, e^{-s t} dt \n", "= \\left( \\int_{-\\infty}^{\\infty} x(t) \\, e^{-s^* t} dt \\right)^* \n", "= X^*(s^*)\n", "\\end{equation}\n", "\n", "where $X(s) = \\mathcal{L} \\{ x(t) \\}$. The ROC for $x^*(t)$ is equal to the ROC for $x(t)$. Above result is used to derive the Laplace transform of the real part of the signal $x(t)$\n", "\n", "\\begin{equation}\n", "\\mathcal{L} \\left\\{ \\Re \\{ x(t) \\} \\right\\} = \\mathcal{L} \\left\\{ \\frac{1}{2} [ x(t) + x^*(t) ] \\right\\} = \\frac{1}{2} \\left[ X(s) + X^*(s^*) \\right]\n", "\\end{equation}\n", "\n", "and imaginary part\n", "\n", "\\begin{equation}\n", "\\mathcal{L} \\left\\{ \\Im \\{ x(t) \\} \\right\\} = \\mathcal{L} \\left\\{ \\frac{1}{2j} [ x(t) - x^*(t) ] \\right\\} = \\frac{1}{2j} \\left[ X(s) - X^*(s^*) \\right]\n", "\\end{equation}\n", "\n", "For a real-valued signal $x(t) \\in \\mathbb{R}$ with $\\Im \\{ x(t) \\} = 0$ it follows that\n", "\n", "\\begin{equation}\n", "X(s) = X^*(s^*)\n", "\\end{equation}\n", "\n", "It can be concluded, that the Laplace transform of a real valued signal shows complex conjugate symmetry. The real/imaginary part of the Laplace transform show even/odd symmetry\n", "\n", "\\begin{align}\n", "\\Re \\{ X(s) \\} &= \\Re \\{ X(-s) \\} \\\\\n", "\\Im \\{ X(s) \\} &= - \\Im \\{ X(-s) \\}\n", "\\end{align}\n", "\n", "respectively and accordingly its magnitude/phase show also even/odd symmetry\n", "\n", "\\begin{align}\n", "|X(s)| &= |X(-s)| \\\\\n", "\\varphi(s) &= - \\varphi(-s)\n", "\\end{align}\n", "\n", "where $\\varphi(s) = \\arg \\{ H(s) \\}$.\n", "\n", "The poles and zeros of a [Laplace transform which can be expressed as a rational function](definition.ipynb#Representation) in $s$ \n", "\n", "\\begin{equation}\n", "X(s) = \\frac{\\sum_{m=0}^{M} \\beta_m s^m}{\\sum_{n=0}^{N} \\alpha_n s^n} = K \\cdot \\frac{\\prod_{\\mu=0}^{Q} (s - s_{0 \\mu})}{\\prod_{\\nu=0}^{P} (s - s_{\\infty \\nu})}\n", "\\end{equation}\n", "\n", "with $Q=M-1$ and $P=N-1$ show also specific symmetries. From the specific symmetry $X(s) = X^*(s^*)$ of the Laplace transform of a real-valued signal $x(t) \\in \\mathbb{R}$, it can be concluded that the coefficients $\\alpha_n, \\beta_m$ of the numerator and denominator are real-valued. The [complex conjugate root theorem](https://en.wikipedia.org/wiki/Complex_conjugate_root_theorem) states that the roots of a polynomial with real valued coefficients are either real-valued or complex conjugate pairs. Hence, the zeros and poles of the Laplace transform $X(s)$ of a real-valued signal $x(t)$ are either real-valued or constitute complex conjugate pairs." ] }, { "cell_type": "markdown", "metadata": { "nbsphinx": "hidden" }, "source": [ "**Copyright**\n", "\n", "This notebook is provided as [Open Educational Resource](https://en.wikipedia.org/wiki/Open_educational_resources). Feel free to use the notebook for your own purposes. The text is licensed under [Creative Commons Attribution 4.0](https://creativecommons.org/licenses/by/4.0/), the code of the IPython examples under the [MIT license](https://opensource.org/licenses/MIT). Please attribute the work as follows: *Sascha Spors, Continuous- and Discrete-Time Signals and Systems - Theory and Computational Examples*." ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.7.10" } }, "nbformat": 4, "nbformat_minor": 1 }