linear_systems_I-solutions.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Linear Systems I - Solutions\n",
    "\n",
    "[back to exercise](linear_systems_I.ipynb)\n",
    "\n",
    "## Preparations"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2020-11-26T13:08:45.486618Z",
     "iopub.status.busy": "2020-11-26T13:08:45.486182Z",
     "iopub.status.idle": "2020-11-26T13:08:45.849175Z",
     "shell.execute_reply": "2020-11-26T13:08:45.849580Z"
    }
   },
   "outputs": [],
   "source": [
    "import tools\n",
    "import sounddevice as sd  # for playback\n",
    "import soundfile as sf  # for reading a soundfile"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "And some other stuff:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2020-11-26T13:08:45.852847Z",
     "iopub.status.busy": "2020-11-26T13:08:45.852392Z",
     "iopub.status.idle": "2020-11-26T13:08:46.000164Z",
     "shell.execute_reply": "2020-11-26T13:08:45.999709Z"
    }
   },
   "outputs": [],
   "source": [
    "# remove \"inline\" to get a separate plotting window:\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "%matplotlib inline"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## One-dimensional time-continuous Systems\n",
    "\n",
    "First, we will have a quick review on linear systems in one dimension. In our case, we will use time signals depending on $t \\in \\mathbb{R}$. Generally, the input signal $x(t) \\in \\mathbb{C}$ and the corresponding output signal $y(t) \\in \\mathbb{C}$ of a system $\\mathcal H$ are related via:\n",
    "\n",
    "$$y(t) = \\mathcal{H}\\{x(t)\\}\\,.$$\n",
    "\n",
    "### Linear Time-Invariant (LTI) Systems\n",
    "\n",
    "As simple as it sounds, LTI-system are linear and time-invariant\n",
    "\n",
    "#### Linearity\n",
    "\n",
    "*Exercise*: Explain the term \"linear\" in your own words.  \n",
    "  \n",
    "<font color='blue'>\n",
    "*\n",
    "The system response (output signal) to a linear combination of different input signals is equal to the linear combination of the system responses to each invididual input signal\n",
    "*\n",
    "</font>  \n",
    "\n",
    "*Exercise*: What does this mean mathematically?\n",
    "\n",
    "<font color='blue'>\n",
    "$$\\mathcal{H}\\{A \\cdot x_1(t) + B \\cdot x_2(t)\\} = A \\cdot \\mathcal{H}\\{x_1(t)\\} + B \\cdot \\mathcal{H}\\{x_2(t)\\}\\, \\text{ for all } A,B \\in \\mathbb{C}$$\n",
    "</font>\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Time-Invariance\n",
    "\n",
    "*Exercise*: Explain the term \"time-invariance\" in your own words.\n",
    "<font color='blue'>\n",
    "*\n",
    "The properties of the system to not change over time: If a distinct input signal causes corresponding system response, a delayed version of that input signal results a delayed version of the system responses.\n",
    "*\n",
    "</font>\n",
    "\n",
    "*Exercise*: What does this mean mathematically?\n",
    "\n",
    "<font color='blue'>\n",
    "$$\\mathcal{H}\\{x(t-\\tau)\\} = y(t-\\tau)\\, \\text{ for all } \\tau \\in \\mathbb{R}$$\n",
    "</font>\n",
    "\n",
    "if $y(t) = \\mathcal{H} \\{ x(t) \\}$ is known.\n",
    "\n",
    "#### Are these systems LTI?\n",
    "\n",
    "*Exercise*: Vote for your LTI-system.\n",
    "\n",
    "1. $\\displaystyle y(t) = a \\cdot x(t) $ with $a \\in \\mathbb{C}$  \n",
    "    Yes:  \n",
    "    No:  \n",
    "    Result: <font color='blue'>Yes</font>\n",
    "2. $\\displaystyle y(t) = a \\cdot x(t) + b $ with $a,b \\in \\mathbb{C}$  \n",
    "    Yes:  \n",
    "    No:  \n",
    "    Result: <font color='blue'>No: time-invariant but non-linear</font>  \n",
    "3. $\\displaystyle y(t) = a \\cdot x(t-t_0) $ with $a \\in \\mathbb{C}$ and $t_0 \\in \\mathbb{R}$  \n",
    "    Yes:  \n",
    "    No:  \n",
    "    Result: <font color='blue'>Yes</font>  \n",
    "4. $\\displaystyle y(t) = a \\cdot x(t-b \\cdot t) $ with $a \\in \\mathbb{C}$ and $b \\in \\mathbb{R}$   \n",
    "    Yes:  \n",
    "    No:  \n",
    "    Result: <font color='blue'>No: linear, but time-variant</font>  \n",
    "5. $\\displaystyle y(t) = \\frac{\\mathrm d x(t)}{\\mathrm d t}$  \n",
    "    Yes:  \n",
    "    No:  \n",
    "    Result: <font color='blue'>Yes</font>  \n",
    "6. $\\displaystyle y(t) = \\int x(t)\\,\\mathrm d t $  \n",
    "    Yes:  \n",
    "    No:   \n",
    "    Result: <font color='blue'>Yes</font>  \n",
    "7. $\\displaystyle y(t) = \\int_{-\\infty}^{\\infty} h(t_0) \\cdot x(t - t_0)\\,\\mathrm d t_0 $  \n",
    "    Yes:  \n",
    "    No:  \n",
    "    Result: <font color='blue'>Yes</font>\n",
    "    \n",
    "#### Listen to a linear and a non-linear system\n",
    "\n",
    "We will investigate two unknown systems. The only information we have about these systems is that the first is LTI (linear and time invariant) and the second is non-linear. They are defined by the functions `tools.blackbox()` and `tools.blackbox_nonlinear()`. Have a quick look at the documentation:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2020-11-26T13:08:46.019826Z",
     "iopub.status.busy": "2020-11-26T13:08:46.019428Z",
     "iopub.status.idle": "2020-11-26T13:08:46.041818Z",
     "shell.execute_reply": "2020-11-26T13:08:46.042290Z"
    }
   },
   "outputs": [],
   "source": [
    "tools.blackbox?"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2020-11-26T13:08:46.046371Z",
     "iopub.status.busy": "2020-11-26T13:08:46.045941Z",
     "iopub.status.idle": "2020-11-26T13:08:46.048446Z",
     "shell.execute_reply": "2020-11-26T13:08:46.048030Z"
    }
   },
   "outputs": [],
   "source": [
    "tools.blackbox_nonlinear?"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "*Exercise:* Load the audio file [data/xmas.wav](data/xmas.wav) and apply both functions to it."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2020-11-26T13:08:46.053254Z",
     "iopub.status.busy": "2020-11-26T13:08:46.052857Z",
     "iopub.status.idle": "2020-11-26T13:08:46.055205Z",
     "shell.execute_reply": "2020-11-26T13:08:46.054777Z"
    }
   },
   "outputs": [],
   "source": [
    "# how to read an audio file\n",
    "sf.read?"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2020-11-26T13:08:46.058067Z",
     "iopub.status.busy": "2020-11-26T13:08:46.057618Z",
     "iopub.status.idle": "2020-11-26T13:08:46.073575Z",
     "shell.execute_reply": "2020-11-26T13:08:46.073977Z"
    }
   },
   "outputs": [],
   "source": [
    "[x, fs] = sf.read('data/xmas.wav')\n",
    "\n",
    "ylin = tools.blackbox(x, fs)\n",
    "ynonlin = tools.blackbox_nonlinear(x, fs)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "*Exercise*: Listen to the input signal and both output signals.\n",
    "\n",
    "**CAUTION: CHECK THE PLAYBACk LEVEL**"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2020-11-26T13:08:46.078109Z",
     "iopub.status.busy": "2020-11-26T13:08:46.077670Z",
     "iopub.status.idle": "2020-11-26T13:08:46.079231Z",
     "shell.execute_reply": "2020-11-26T13:08:46.079653Z"
    }
   },
   "outputs": [],
   "source": [
    "# how to play back the signal\n",
    "sd.play?"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2020-11-26T13:08:46.082755Z",
     "iopub.status.busy": "2020-11-26T13:08:46.082325Z",
     "iopub.status.idle": "2020-11-26T13:08:54.287531Z",
     "shell.execute_reply": "2020-11-26T13:08:54.287945Z"
    }
   },
   "outputs": [],
   "source": [
    "sd.play(x, blocking=True)\n",
    "sd.play(ylin, blocking=True)\n",
    "sd.play(ynonlin, blocking=True)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### The Impulse Response\n",
    "\n",
    "The impulse response $h(t)$ of an LTI system characterises it completely. It is the LTI system's response to a Dirac impulse $\\delta(t)$\n",
    "\n",
    "$$h(t) = \\mathcal{H}\\{\\delta(t)\\}\\,.$$\n",
    "\n",
    "So why is the impulse response sufficient to describe the whole LTI system? The input signal $x(t)$ can be described as a sequence of Dirac impulses \n",
    "\n",
    "$$ x(t) = \\int_{-\\infty}^{\\infty} x(t_0) \\cdot \\delta(t-t_0)\\,\\mathrm d t_0\\,,$$ \n",
    "\n",
    "where the Dirac impulse at $t_0$ is weighted by the value $x(t_0)$ of the signal at $t_0$. Applying the system onto $x(t)$ yields\n",
    "\n",
    "$$ y(t) = \\mathcal{H}\\{x(t)\\} = \\mathcal{H}\\left\\{\\int_{-\\infty}^{\\infty} x(t_0) \\cdot \\delta(t-t_0)\\,\\mathrm d t_0\\right\\}$$\n",
    "\n",
    "As a next step, we can exchange the integral operator $\\int$ and the system operator $\\mathcal{H}$:\n",
    "\n",
    "$$ y(t) = \\mathcal{H}\\{x(t)\\} = \\int_{-\\infty}^{\\infty} x(t_0) \\cdot \\mathcal{H}\\{\\delta(t-t_0)\\}\\,\\mathrm d t_0$$\n",
    "\n",
    "*Exercise*: What property has to be fulfilled by $\\mathcal{H}$ in order to be able to exchange integral and system operator?"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<font color='blue'>\n",
    "*\n",
    "Linearity  \n",
    "*\n",
    "</font>"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "As the last step, we re-write the system response of the Dirac $\\delta$ shifted about $t_0$ as the shifted impulse response $h(t-t_0)$:\n",
    "\n",
    "$$ y(t) = \\mathcal{H}\\{x(t)\\} = \\int_{-\\infty}^{\\infty} x(t_0) \\cdot h(t-t_0)\\,\\mathrm d t_0$$\n",
    "\n",
    "*Exercise*: What property has to be fulfilled by $\\mathcal{H}$ in order to replace $\\mathcal{H}\\{\\delta(t-t_0)\\}$ by $h(t-t_0)$?"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<font color='blue'>\n",
    "*\n",
    "Time-Invariance\n",
    "*\n",
    "</font>"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Hence, we can describe the output signal $y(t)$ by the so-called **linear convolution** integral of the corresponding input signal $x(t)$ and the impulse response $h(t)$. Its short version reads\n",
    "\n",
    "$$ y(t) = x(t) * h(t) $$\n",
    "\n",
    "where $*$ is a common notation of the convolution."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### A Naive Implementation of the Linear Convolution\n",
    "\n",
    "Time-continuous signals cannot be handled by computers. Signals must be sampled in time with the sample period $T_s$.\n",
    "This yields discrete-time signals. The discrete-time counterpart of a linear convolution is given as\n",
    "\n",
    "$$ y[n] = x[n] \\ast h[n] = \\sum_{k = -\\infty}^{\\infty} x[k] \\cdot h[n-k] $$\n",
    "\n",
    "where $y[n] = y(n T_s)$, $x[n] = x(n T_s)$, and $h[n] = h(n T_s)$ denote the discrete-time versions of the involved entities.\n",
    "\n",
    "*Exercise:* Write a function called `naive_convolution()` that computes the convolution of two one-dimensional arrays (discrete-time signals of finite length) by means of two nested loops according to the equation above, where $x$ and $h$ are one-dimensional arrays of finite lengths. The infinite sum can be changed to a finite sum by assuming that all values before index 0 and all values after the last array element are equal to zero.\n",
    "\n",
    "Following this assumption, at which indices $n$ does $y[n]$ have its first and last non-zero value?"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2020-11-26T13:08:54.308282Z",
     "iopub.status.busy": "2020-11-26T13:08:54.307748Z",
     "iopub.status.idle": "2020-11-26T13:08:54.424036Z",
     "shell.execute_reply": "2020-11-26T13:08:54.424503Z"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Text(0.5, 0, 'n')"
      ]
     },
     "execution_count": 1,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAEGCAYAAABo25JHAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAURklEQVR4nO3df4xd5X3n8fc3g9l680NeyaMCYxdnK8sShG0djVxcqgi1zRpYFLxsVkPU/KIbGSLoBrV1FEfbsKnUotYS2qQkeFzwFtIEaBLH6yQGJ82WJZbGMWMbx9hmFocEPB4nHkBj4zDEP/j2j3sZDcN4PDOeM3fmPu+XdDX3Oee553wfe+79zPl1T2QmkqRyva3RBUiSGssgkKTCGQSSVDiDQJIKZxBIUuEuaHQBEzV//vxctGhRo8uQpFll586dL2Zm62jzZl0QLFq0iO7u7kaXIUmzSkQ8f7Z57hqSpMIZBJJUOINAkgpnEEhS4QwCSSpcZWcNRcSvAU8A/6a+nm9k5p0j+gTwBeA64FXg45m5q6qaNLNs2n2YtVt76BsY5JJ5c1m9Ygkrl7Y1uqxJaZaxNMs4NDFVnj76K+D3M/NERMwBtkXEo5m5fVifa4HF9cfvAPfWf6rJbdp9mDUb9zJ46gwAhwcGWbNxL8Cs++BplrE0yzg0cZXtGsqaE/XmnPpj5Hde3wA8WO+7HZgXERdXVZNmjrVbe4Y+cN4weOoMa7f2NKiiyWuWsTTLODRxlR4jiIiWiHgKOAp8PzN/NKJLG3BoWLu3Pm3kclZFRHdEdPf391dXsKZN38DghKbPZM0ylmYZhyau0iDIzDOZ+dvAAmBZRLxnRJcY7WWjLGd9ZrZnZntr66hXSGuWuWTe3AlNn8maZSzNMg5N3LScNZSZA8DjwDUjZvUCC4e1FwB901GTGmv1iiXMndPypmlz57SwesWSBlU0ec0ylmYZhyausiCIiNaImFd/Phf4Q+CZEd02Ax+NmiuBY5l5pKqaNHOsXNrGXTdewYUttV/BtnlzuevGK2blQclmGUuzjEMTV+VZQxcDD0REC7XA+afM/E5E3AqQmeuALdROHT1I7fTRmyusRzPMyqVtPLTjBQAeuWV5g6s5P80ylmYZhyamsiDIzB8DS0eZvm7Y8wRuq6oGSdK5eWWxJBXOIJCkwhkEklQ4g0CSCmcQSFLhDAJJKpxBIEmFMwgkqXAGgSQVziCQpMIZBJJUOINAkgpnEEhS4QwCSSqcQSBJhTMIJKlwBoEkFc4gkKTCGQSSVDiDQJIKZxBIUuEMAkkqnEEgSYUzCCSpcJUFQUQsjIh/iYgDEbEvIj41Sp+rI+JYRDxVf3yuqnokSaO7oMJlnwb+LDN3RcQ7gZ0R8f3M3D+i3w8z8/oK65AkjaGyLYLMPJKZu+rPXwEOAG1VrU+SNDnTcowgIhYBS4EfjTJ7eUTsiYhHI+Lys7x+VUR0R0R3f39/hZVKUnkqD4KIeAfwTeCOzDw+YvYu4NLM/C3g74BNoy0jM9dnZntmtre2tlZbsCQVptIgiIg51ELgq5m5ceT8zDyemSfqz7cAcyJifpU1SZLerMqzhgK4HziQmXefpc9F9X5ExLJ6PS9VVZMk6a2qPGvoKuAjwN6IeKo+7bPAbwBk5jrgg8AnI+I0MAjclJlZYU2SpBEqC4LM3AbEOfrcA9xTVQ2SpHPzymJJKpxBIEmFMwgkqXAGgSQVziCQpMIZBJJUOINAkgpnEEhS4QwCSSqcQSBJhTMIJKlwBoEkFc4gkKTCGQSSVDiDQJIKZxBIUuEMAkkqnEEgSYUzCCSpcAaBJBXOIJCkwhkEklQ4g0CSCmcQSFLhLqhqwRGxEHgQuAh4HVifmV8Y0SeALwDXAa8CH8/MXVXV1Aw27T7M2q099A0Mcsm8uaxesYSVS9saXZY04/heGb/KggA4DfxZZu6KiHcCOyPi+5m5f1ifa4HF9cfvAPfWf2oUm3YfZs3GvQyeOgPA4YFB1mzcC+AvuDSM75WJqWzXUGYeeeOv+8x8BTgAjPwfuAF4MGu2A/Mi4uKqaprt1m7tGfrFfsPgqTOs3drToIqkmcn3ysRMyzGCiFgELAV+NGJWG3BoWLuXt4YFEbEqIrojoru/v7+qMme8voHBCU2XSuV7ZWIqD4KIeAfwTeCOzDw+cvYoL8m3TMhcn5ntmdne2tpaRZmzwiXz5k5oulQq3ysTU2kQRMQcaiHw1czcOEqXXmDhsPYCoK/Kmmaz1SuWMHdOy5umzZ3TwuoVSxpUkTQz+V6ZmMqCoH5G0P3Agcy8+yzdNgMfjZorgWOZeaSqmma7lUvbuOvGK7iwpfbf1jZvLnfdeIUHv6QRfK9MTJVnDV0FfATYGxFP1ad9FvgNgMxcB2yhduroQWqnj95cYT1NYeXSNh7a8QIAj9yyvMHVSDOX75XxqywIMnMbox8DGN4ngduqqkGSdG5eWSxJhTMIJKlwBoEkFc4gkKTCGQSSVDiDQJIKZxBIUuEMAkkqnEEgSYUzCCSpcAaBJBXOIJCkwhkEklQ4g0CSCmcQSFLhDAJJKpxBIEmFMwgkqXAGgSQVziCQpMKd9eb1EXHjOF7/WmZumcJ6JEnT7KxBAPw98H+AGKPP+wCDQJJmsbGC4NHM/OOxXhwR/zjF9UiSptlZjxFk5ofP9eLx9JEkzWxjbREMiYjfBRYN75+ZD57jNRuA64GjmfmeUeZfTW3X00/rkzZm5l+Oq2pJ0pQ5ZxBExFeA3wSeAs7UJycwZhAA/wDcc45+P8zM689dpiSpKuPZImgHLsvMnMiCM/OJiFg0maIkSdNnPNcRPA1cVNH6l0fEnoh4NCIuP1uniFgVEd0R0d3f319RKZJUpvFsEcwH9kfEDuBXb0zMzA+c57p3AZdm5omIuA7YBCwerWNmrgfWA7S3t09oy0SSNLbxBMH/rGLFmXl82PMtEfHliJifmS9WsT5J0ujOGQSZ+f+qWHFEXAT8IjMzIpZR2031UhXrkiSd3VhfMfGdc53RM1afiHgIuBqYHxG9wJ3AHIDMXAd8EPhkRJwGBoGbJnpAWpJ0/sbaIvi9iNg8xvwALjvbzMz80Fgrzsx7qJ1eKklqoLGC4Ib6z/8A9AIvj9Ln5JRXJEmaVmcNgjeODUTE+4FPUTvLZwOw1V04ktQ8znkdQWb+D2qndd4PfBx4NiL+OiJ+s+LaJEnTYFw3pqlvAfy8/jgN/DvgGxHxtxXWJkmaBuP5rqH/DnwMeBG4D1idmaci4m3As8Cnqy1RklSl8V5ZfGNmPj98Yma+HhF+YZwkzXLjuaDsc2PMOzC15UiSpps3r5ekwhkEklQ4g0CSCmcQSFLhDAJJKpxBIEmFMwgkqXAGgSQVziCQpMIZBJJUOINAkgpnEEhS4QwCSSqcQSBJhTMIJKlwBoEkFc4gkKTCjedWlZMSERuA64GjmfmeUeYH8AXgOuBV4OOZuauKWjbtPszarT30DQxyyby5rF6xhJVL26pYlSRNuao/w6rcIvgH4Jox5l8LLK4/VgH3VlHEpt2HWbNxL4cHBkng8MAgazbuZdPuw1WsTpKm1HR8hlUWBJn5BPDyGF1uAB7Mmu3AvIi4eKrrWLu1h8FTZ940bfDUGdZu7ZnqVUnSlJuOz7BGHiNoAw4Na/fWp71FRKyKiO6I6O7v75/QSvoGBic0XZJmkun4DGtkEMQo03K0jpm5PjPbM7O9tbV1Qiu5ZN7cCU2XpJlkOj7DGhkEvcDCYe0FQN9Ur2T1iiXMndPypmlz57SwesWSqV6VJE256fgMa2QQbAY+GjVXAscy88hUr2Tl0jbuuvEKLmypDbVt3lzuuvEKzxqSNCtMx2dYlaePPgRcDcyPiF7gTmAOQGauA7ZQO3X0ILXTR2+uqpaVS9t4aMcLADxyy/KqViNJlaj6M6yyIMjMD51jfgK3VbV+SdL4eGWxJBXOIJCkwhkEklQ4g0CSCmcQSFLhDAJJKpxBIEmFMwgkqXAGgSQVziCQpMIZBJJUOINAkgpnEEhS4QwCSSqcQSBJhTMIJKlwBoEkFc4gkKTCGQSSVDiDQJIKZxBIUuEMAkkqnEEgSYUzCCSpcJUGQURcExE9EXEwIj4zyvyrI+JYRDxVf3yuynokSW91QVULjogW4EvA+4Fe4MmI2JyZ+0d0/WFmXl9VHZKksVW5RbAMOJiZz2XmSeBh4IYK1ydJmoQqg6ANODSs3VufNtLyiNgTEY9GxOWjLSgiVkVEd0R09/f3V1GrJBWryiCIUabliPYu4NLM/C3g74BNoy0oM9dnZntmtre2tk5xmZJUtiqDoBdYOKy9AOgb3iEzj2fmifrzLcCciJhfYU2SpBGqDIIngcUR8e6IuBC4Cdg8vENEXBQRUX++rF7PSxXWJEkaobKzhjLzdETcDmwFWoANmbkvIm6tz18HfBD4ZEScBgaBmzJz5O4jSVKFKgsCGNrds2XEtHXDnt8D3FNlDZKksXllsSQVziCQpMIZBJJUOINAkgpnEEhS4QwCSSqcQSBJhTMIJKlwBoEkFc4gkKTCGQSSVDiDQJIKZxBIUuEMAkkqnEEgSYUzCCSpcAaBJBXOIJCkwhkEklQ4g0CSCmcQSFLhDAJJKpxBIEmFMwgkqXCVBkFEXBMRPRFxMCI+M8r8iIgv1uf/OCLeW2U9kqS3qiwIIqIF+BJwLXAZ8KGIuGxEt2uBxfXHKuDequqRJI2uyi2CZcDBzHwuM08CDwM3jOhzA/Bg1mwH5kXExVUVtP/Icb7efQiAU2dep6Ozi2/t7gVg8OQZOjq7+PaePgCOv3aKjs4uHnv6CAAv//IkHZ1d/PP+XwBw9JXX6Ojs4vGeowD0DQzS0dnFtmdfBOCFl16lo7OL7c+9BMBP+k/Q0dnFzudfBqDn56/Q0dnFnkMDAOzrO0ZHZxf7+o4BsOfQAB2dXfT8/BUAdj7/Mh2dXfyk/0StvsFafS+89CoA2559kY7OLvoGBgF4vOcoHZ1dHH3lNQD+ef8v6Ojs4uVfngTgsaeP0NHZxfHXTgHw7T19dHR2MXjyDADf2t1LR2cXp868DsDXuw/R0dk19G/50I4X+KP7tg+1v9L1Mz62YcdQe8O2n/KJB54caq9/4ifc+pWdQ+0vP36Q27+2a6j9xR88yx0P7x5q3/29Hv7863uG2n/z2DOs2fjjofZffXc/f7Hp6aH257+9j89/e99Q+y82Pc1ffXf/UHvNxh/zN489M9T+86/v4e7v9Qy173h4N1/8wbND7du/tosvP35wqH3rV3ay/omfDLU/8cCTbNj206H2xzbs4BfHXxtq/9F923loxwtD7Y7Orlnzu3fiV6fZf+T4WX/3tj/30qz63TtybLCpfveqUGUQtAGHhrV769Mm2oeIWBUR3RHR3d/fP6liHrllOZdd/K5JvXameeSW5dz5gcsbXcaUeOSW5Txyy/JGlzEl/tvvvbspxvLX//mKpnqvfPjKSxtdxpSo8r0SmVnNgiP+K7AiMz9Rb38EWJaZfzKsz3eBuzJzW739A+DTmblztGUCtLe3Z3d3dyU1S1Kzioidmdk+2rwqtwh6gYXD2guAvkn0kSRVqMogeBJYHBHvjogLgZuAzSP6bAY+Wj976ErgWGYeqbAmSdIIF1S14Mw8HRG3A1uBFmBDZu6LiFvr89cBW4DrgIPAq8DNVdUjSRpdZUEAkJlbqH3YD5+2btjzBG6rsgZJ0ti8sliSCmcQSFLhDAJJKpxBIEmFq+yCsqpERD/w/CRfPh94cQrLaSTHMjM1y1iaZRzgWN5waWa2jjZj1gXB+YiI7rNdWTfbOJaZqVnG0izjAMcyHu4akqTCGQSSVLjSgmB9owuYQo5lZmqWsTTLOMCxnFNRxwgkSW9V2haBJGkEg0CSCldMEETENRHRExEHI+Izja5nsiJiQ0QcjYinz9175oqIhRHxLxFxICL2RcSnGl3TZEXEr0XEjojYUx/L5xtd0/mKiJaI2B0R32l0LecjIn4WEXsj4qmImLV3tIqIeRHxjYh4pv6emdJblRVxjCAiWoD/D7yf2s1wngQ+lJn7x3zhDBQR7wNOULvX83saXc9k1e9NfXFm7oqIdwI7gZWz9P8kgLdn5omImANsAz5Vvw/3rBQRfwq0A+/KzOsbXc9kRcTPgPbMnNUXlEXEA8APM/O++v1d/m1mDkzV8kvZIlgGHMzM5zLzJPAwcEODa5qUzHwCeLnRdZyvzDySmbvqz18BDjDK/apng6w5UW/OqT9m7V9YEbEA+E/AfY2uRRAR7wLeB9wPkJknpzIEoJwgaAMODWv3Mks/dJpRRCwClgI/amwlk1fflfIUcBT4fmbO2rEA/wv4NPB6owuZAgl8LyJ2RsSqRhczSf8e6Af+d3133X0R8fapXEEpQRCjTJu1f7E1k4h4B/BN4I7MPN7oeiYrM89k5m9Tu+/2soiYlbvtIuJ64Ghm7mx0LVPkqsx8L3AtcFt91+pscwHwXuDezFwK/BKY0uOcpQRBL7BwWHsB0NegWlRX35/+TeCrmbmx0fVMhfom++PANQ0uZbKuAj5Q37f+MPD7EfGPjS1p8jKzr/7zKPAtaruJZ5teoHfYVuY3qAXDlCklCJ4EFkfEu+sHWm4CNje4pqLVD7DeDxzIzLsbXc/5iIjWiJhXfz4X+EPgmcZWNTmZuSYzF2TmImrvk/+bmR9ucFmTEhFvr5+IQH1Xyn8EZt3Zdpn5c+BQRCypT/oDYEpPqqj0nsUzRWaejojbga1AC7AhM/c1uKxJiYiHgKuB+RHRC9yZmfc3tqpJuQr4CLC3vm8d4LP1+1zPNhcDD9TPTnsb8E+ZOatPu2wSvw58q/Y3BxcAX8vMxxpb0qT9CfDV+h+yzwE3T+XCizh9VJJ0dqXsGpIknYVBIEmFMwgkqXAGgSQVziCQpMIZBJJUOINAkgpnEEjnKSIW1b8j/u/r9yP4Xv0KY2lWMAikqbEY+FJmXg4MAP+lwfVI42YQSFPjp5n5xldl7AQWNbAWaUIMAmlq/GrY8zMU8j1eag4GgSQVziCQpML57aOSVDi3CCSpcAaBJBXOIJCkwhkEklQ4g0CSCmcQSFLhDAJJKty/AqcseoH//wYJAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "def naive_convolution(x, h):\n",
    "    # in python, you have to indent everything inside your function\n",
    "\n",
    "    Nx = len(x)  # length of x\n",
    "    Nh = len(h)  # length of h\n",
    "    Ny = Nx + Nh - 1  # resulting length of y\n",
    "\n",
    "    y = np.zeros(Ny)  # initialise output array\n",
    "    for k in np.arange(0, Nx):  # help: for which indices k is x non-zero?\n",
    "        for n in np.arange(k, Nh+k):  # help: for which indices (n-k) is h non-zero?\n",
    "            y[n] = y[n] + x[k] * h[n-k]\n",
    "\n",
    "    return y\n",
    "    # end of function\n",
    "\n",
    "\n",
    "# try out your function\n",
    "x = np.array([1, 1, 1, 1, 1])\n",
    "h = np.array([1, 1, 1])\n",
    "\n",
    "y = naive_convolution(x, h)\n",
    "plt.stem(y, use_line_collection=True, basefmt='C0:')\n",
    "plt.ylabel('y[n]')\n",
    "plt.xlabel('n');"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### The Transfer Function\n",
    "\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The transfer function of an LTI-system is the temporal Fourier transform its impulse response"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "$$ H(\\omega) = \\int_{-\\infty}^{\\infty} h(t) e^{-j\\omega t} \\mathrm d t$$"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "*Exercise*: $y(t)$, $x(t)$ and $h(t)$ are related to each other by the convolution. How are the respective Fourier spectra $Y(\\omega)$, $X(\\omega)$, $H(\\omega)$ related?"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    " "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### A naive implementation of the Fourier Transform\n",
    "\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Time-continuous signals cannot be handled by computers. They are sampled in time with the sample period $T_s$. In order to compute the Fourier transform of a signal, one has to discretize the integral\n",
    "\n",
    "$$ H(\\omega) = \\sum_{n=-\\infty}^{\\infty} h[n] e^{-j\\omega n T_s} $$\n",
    "\n",
    "where $h[n] = h(n T_s)$ denotes the discrete version of the involved $h(t)$. This is also known as DTFT, discrete-time Fourier transform.\n",
    "\n",
    "*Exercise:* Write a function called `naive_ft()` that computes the Fourier transform of a given signal $x$ of finite length for different frequencies. Use two nested loops according for this, again. The infinite sum can be changed to a finite sum by assuming that all values before index 0 and all values after the last array element are equal to zero."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2020-11-26T13:08:54.429356Z",
     "iopub.status.busy": "2020-11-26T13:08:54.428861Z",
     "iopub.status.idle": "2020-11-26T13:08:54.431154Z",
     "shell.execute_reply": "2020-11-26T13:08:54.430662Z"
    }
   },
   "outputs": [],
   "source": [
    "def naive_ft(x, f, fs):\n",
    "    # inputs:\n",
    "    # x - signal vector\n",
    "    # f - time-frequencies\n",
    "    # fs - sample rate\n",
    "\n",
    "    # outputs:\n",
    "    # X - frequency spectrum\n",
    "\n",
    "    Nsig = len(x)  # length of signal\n",
    "    Nspec = len(f)  # length of spectrum\n",
    "\n",
    "    omega = 2*np.pi*f  # angular frequency\n",
    "\n",
    "    X = np.zeros(Nspec, dtype=complex)  # initialise output spectrum\n",
    "    for k in np.arange(0, Nspec):  # angular frequency\n",
    "        for n in np.arange(0, Nsig):  # help: for which indices (n-k) is h non-zero?\n",
    "            X[k] = X[k] + x[n]*np.exp(-1j*omega[k]*n/fs)\n",
    "\n",
    "    return X\n",
    "    # end of function"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "*Exercise:* Compute the transfer function of the blackbox system `tools.blackbox`."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {
    "execution": {
     "iopub.execute_input": "2020-11-26T13:08:54.436783Z",
     "iopub.status.busy": "2020-11-26T13:08:54.436143Z",
     "iopub.status.idle": "2020-11-26T13:08:55.498636Z",
     "shell.execute_reply": "2020-11-26T13:08:55.498221Z"
    }
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAcEAAAEaCAYAAABpQuwEAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nO3deZwcdZ3/8ddn7nsmmUkm5CAJJOFWjgCKIAOigD+VVVdlf6xr3F1YXNDVPWQRd5ddl1Xcn7qHJyrLCgqiK4quCkYZIEAICaKYiwy5CTnmSDKTZM7+/P6o6tAZ5u6eqe6u9/PxmMd0V1VXffrb1fXp7/db9S1zd0REROKoIOoAREREoqIkKCIisaUkKCIisaUkKCIisaUkKCIisaUkKCIisZU3SdDM3mBmm8ysy8x+bwzLLzAzN7OiDGy7ycx2jjD/LjP758naxmjbl5GZ2TVm9nBE23YzWxTFtkVkDEnQzLaa2REz6zSz/Wb2pJldb2YF4fyfhYmny8z6zKw35flXwwN0ImVal5n9OJyXfN4bvjb5/GcTeC//BHzR3avc/YfDvI/LJrBemSRmdquZ3RN1HO7+bXd/S9RxyNQY7pgUdVxRMLNlZjYQlsFBM/uNmb0t6rim0lhrQW939+VmVgtcDPw7cD7wQXe/MrmQmd0F7HT3T6ZMawJ2ufvcIdZ7fbjMrcAid//DibyJ0HxgbRqvF8k7Zlbk7v1Rx5GFhjsmHSMm5feUu18YVmyuBe4zs7nuvj/qwKbCuJpD3f2Auz8IvA/4gJmdPjlhDc3MrjWzFjNrN7MHzWx2OP1F4ATgx+EvmtJBr7sbOD5l/sdTZl9jZtvNrNXMbkl5TYGZ/a2ZvWhmbWZ2v5lNHyW+T4Tr2Wpm1wyzzDQz+4mZ7TOzjvDx3JT5083sv8xsVzj/VbXacLmPmNm6Qa8dcvtmVmtm3wq3uc3MPplSk/+KmX0/ZdnbzeyXZmZDbHORmT1qZgfC7Xw3nP4lM/vcoGV/bGYfDR/fZGYvha0JG83sTWZ2BfAJ4H3hZ/KblFi/aWYvh6/5ZzMrDOctM7MnzOwLFrRKbDazC8LpO8xsr5l9ICWGt4Zl1Bmu66+HKctlZrYi5bmb2Z9b0LzeaWafMrMTzeyp8Nfy/WZWEi7bZGY7Ryj7ZjP70+G2NSiOYeM1s7eZ2XP2SmvMa4ZaR0r8N5jZJmBTOO1kM/uFBd+djWb23tG2O4b3NtJ+tczMVpjZ/wv34y1mlvqDeVn4+XWG81LX+8dmtj583UNmNn+495ppFrROfN/M7jGzg8AyG+VYYGbvD99/m5ndYimtTjaoK8QGdV2Y2Wwz+5+wDLeY2UcGxXJ/WMadZrbWzJamzJ9nZj8IX9tmZl80s9LwMz4jZbmZFrTmzRjpvbt7ArgbqAQWp7z+deE+t9+CmmJTyrwhP0d75bv6nxYcLzaY2ZsGve8Hw1hbzOzacbzvVx1PwunjPmYn3/iIf8BW4LIhpm8HPjRo2l3APw+a1kRQOxxpG7cC94yyzKVAK3A2UAr8J/DYaHEONx9YADjwdaAceC3QA5wSzv8osBKYG27va8C9w6y7CegHPh8uezFwCDhpcLkA9cC7gQqgGvge8MOUdf0v8F1gGlAMXDy4HIG/A54FZoxx+98CfhRubwHwAvAn4byK8Pky4KKwjOcO8z7vBW4h+PFUBlwYTj8P2AUUhM8bgMNAI3ASsAOYnVLuJw73uQM/DMu6EpgJrAL+LJy3LHyfHwQKgX8m2A+/FL7vtwCdQFW4/MvAReHjacDZw7yvZcCKlOcOPAjUAKeF+8UvCX5o1QLrgA+MseybgT8dZVuLRoqXYJ/fS9D6Ugh8gGB/Lh3m/TjwC2A6wb5dGX4GHyRo/Tk7/JxPG2W7o723kfarZUAfQc2iEPgQwT5iYTwHU9ZzXEosvwe0AKeEsX4SeHKE7/X+Ef7+doTv65DHJIJ9si+MoyAsv2GPBcCpQBfwxnDe58Myu2yoYyLHfo8LgDXA3wMlBPvXZuDylFi6gbeGZfhpYGU4rxD4DfCFsDxTv49fBm5P2eZfAD8ebd8P13kD0AvMDKfNAdrCGAqAN4fPZ4zyOS4Ly+FjBMex9wEHgOnh/EfDOMuAM4F9wJvG8L5HOp6M+Zh9TBmMusDwSXAlcMugacd84CkfeoJjd873DrHjjZYEvwl8NuV5FcHOumCkOId7H7ySBOemTFsFXB0+Xp/8UFI+4D6gaJgvVT9QmTLtfuDvhiuXlOXOBDpStpEApg2zjZcIvmQrgNqxbD/ciXqAU1Pm/RnQnPL8PKAd2Ab8wQhl+C3gDoZIkmF5vTl8fCPw0/DxIoID+GVA8UifO0HS7AHKU6b9AfBIyhdrU8q8M8LPsDFlWhtwZvh4e/hea0bZt5bx6sT0hpTna4CbUp5/Dvi3MX72zYw9CQ4ZL/AV4FODpm0k/IE0xPtx4NKU5+8DHh+0zNeAfxhluxPer8L32ZIyryKMaxbBwXM/wY/B8kHb/BlhIg2fFxD8oJo/0mc4nj9GOCaF++Rjg5Yf9lhAkMDuS5lXSZBExpIEzwe2D9rWzcB/pcSyPGXeqcCR8PHrCRLHUMej8wkSRfJH6WoGHXMH7Y/9YRn0AUdSlwVuAu4e9JqHCH6IjfQ5LiP80ZMybRXwfmAeMABUp8z7NHDXGN73SMeTMR+zU//SOTt0DsGBcyx2uXtdyt/9E9jebIKDNADu3kVwwJszgXWl2p3y+DBBcoWgj/GBsAlgP0EBDxAcqIfS4e6HUp5vC2M+hplVmNnXwuaTg8BjQJ0FTX7zgHZ37xhmG3XAdcCn3f3AGLffQPArc9ugeUfLzd1XEfwCNYKD3HA+Hi6zKmyi+OOUef8NJPt0/5CgWQV3byH4hXYrsNfM7rOwGXsI8wl+Nb6cUu5fI6gRJu1JeXwk3MbgacnP8N0Evya3WdCM+/oR3ttgg9c53DZgjJ/9GAwX73zgr5JlEpbLvFG2sSPl8Xzg/EGvv4YgIY203ZHe26j7FSnfLXc/HD6sCtf3PoJzAl42s/81s5NTYv33lDjbCfa5dL/ng410TNoxaNmRjgWzU5cP31vbGGOYD8we9Ll8gmOPMYOPT2UWnNE+D9jmQ/RXuvvTBDX2i8NyXUTQsjGcle5eR9AK8CBBi1BqjO8ZFOOFwHGjfI4AL3mYjULJfWc2wXGuc9C8Ifed1Pc9yvFkvMdsYIKXSJjZuWHAQ/ZtTJJdBG8yGUMlQdPiS2N8vY++yDF2AFcO+qKUuftw25sWxpR0fBjzYH9FUKU/391rCJpRIPii7wCmm1ndMNvoAN4G/JeZvWGM228l+DU0f9C8o+/DzG4gaD7YRZDohuTuu939WnefTfCr/8v2yun99wBXmdlrCZqyfpjyuu+4+4VhDA7cnpw1aBM7CGoXDSllXuPupw0X00jc/Rl3v4ogif6QkRN8Okb67A8R1IKSZjGMEeLdAdw2aF+scPd7R4gptWx3AI8Oen2Vu39olO2O9N5G3a9G4u4PufubCX6tbyDolkjG+meDYi139yeHWo8de4bn4L9PjCWWocIb9HykY8HLBAkpGU8FwXEpaaTPfwewZdB6q939rWOIcQdwvA1/iVfyR+n7ge+7e/doKwwrFn8OvN/MzkrZzt2DYqx098+ErxnucwSYY3bMuQXJfWcXwXGuetC8se47wx1PxnvMBsaZBM2sxoLTZ+8jaMZ6fjyvT9N3gA+a2ZkWnPjyL8DT7r51jK/fQ9DmPlZfBW6zsFPezGaY2VWjvOYfzazEzC4iSFbfG2KZaoKaxP6w0/YfkjPc/WWC5qAvW3ACTbGZvTH1xe7eTPAr/gEzO3+07bv7AMFB7TYzqw7fz18SJC3MbAlB31ryC/NxMztzqDdnZu+xV07E6SDYAQfCuHYCzxDUAP/H3Y+ErznJzC4NP7Pu8L0PhOvYAyyw8GSK8P0/DHwu3NcKLDgh5eKh4hlJWA7XmFmtu/cR9F0MjPa6NAz32T8HvCtsAVgE/MkE4v06cL2ZnW+BSjP7P4MOIiP5CbDEghM4isO/c83slDGW07j3q5GYWaOZvSNMrj0EfWrJbX4VuNnMTguXrTWz9wy3rjCZD/f3L2Msn9GMdCz4PvA2M7vQgpOl/oljj6vPAW+14IS3WQS1mKRVwEELTvQoN7NCMzs9rGSMZhVBAv5MuD+UDfphfDfwToLv9bfG+kbdvQ34BkEzLwSf59vN7PIwvjILTu6ZO8rnCMGPqo+E+9t7CH4c/9TddwBPAp8O1/cagu/Ft0eLb5TjyUSO2WNOgj82s06CTHsLQb/UB8f42oxw918S9EX8D8GHfyJw9ThW8WngkxZUlYc8S3CQfydoGng4fO8rCdrah7ObIDHsIvgwr3f3DUMs928Ene2t4Tp/Pmj++wl+YW8gaPv+6KD5uPsvCMr/QTM7Zwzb/zDBL9LNBLX37wB3hr8i7yHoRP+Nu28iaI652wadYRs6F3jazLrCsvkLd9+SMv+/Cfrp7k6ZVgp8Jny/uwm+GMlf6MlE0WZmz4aP/4igmW1d+H6+T/ArcyLeD2y1oNn5el5prs20kcr+CwR9RHsIymekL/qQ8br7aoITTL4YbqeFoM9lTMJmp7cQfF92hfHeTvDZDLvdMby3IferMYRUQNAisougufNighoI7v5AGNt9YTy/A64cZj1TZdhjgbuvJTiZ5DsEx6UOIHXgirsJTmDZSvAD77vJGeEPibcTnBewheA78g2Ck69GlPLaRQR9ujsJmiaT83cSnDznwOPjfL//RpC4XxMmrKsIvrP7CHLA3xB8hsN+jqGnCc4ybQVuA34/TLIQ9PUvCF/7AEH/9C/GENtIx5PxHrOBsNNSJBPCWus9BCcrJaKOZypYcLr4PT6Ga85yTT6/t8lkZlsJToZaHnEcdxL0fX5y1IUzv+1lBGVw4VRve7zSHjJMBMDMiglOxf5GXBKgSLYyswXAu4CzRl5S8mbsUImOmZ1CcKr0cQRNKSISETP7FEEz8r8O6q6QIag5VEREYks1QRERiS0lQRERia1YnBjT0NDgCxYsiDqMSB06dIjKysrRFxRA5TURKrPxyYXyWrNmTau7jzjwdq6LRRJcsGABq1evjjqMSDU3N9PU1BR1GDlD5TV+KrPxyYXyMrNtoy+V29QcKiIisaUkKCIisaUkKCIisaUkKCIisaUkKCIisaUkKCIisaUkKCIisaUkKCIisaUkKCIisaUkKCIisaUkKCIisZWTSdDMrjCzjWbWYmZ/G3U8IiKSm3IuCZpZIfAl4ErgVOAPzOzUaKMSEZFclIt3kTgPaHH3zQBmdh9wFbAu0qhEJllrVw/7D/dxsLuP7t4ByksKqSwtoqKkkLLiQkqKCigpDH7XusOAO339CXr6E/T2J+juH+BI7wDdfQP0hNO7+wboG0jQN5Cgd8AZGEgw4ODuuL+ybTMoLLCjfyWFBUe3V1pcQGlRIZs6BqjfeYCy4oKj8ZQWBcsVFhiFFrzWHRLuDLjTG8aWGk9Pf4KBROLo9kuKCqgsLaKypIjqsuD9mtkxZTOQcDq7+9h/uI+EO7Nqy6goycXDm0y1XNxL5gA7Up7vBM4fvJCZXQdcB9DY2Ehzc/OUBJeturq6Yl8G45FN5bX7UIJvr+/l+daBqEMZ3dMrJn0TBQYVRUFS7k84AwnoGQAftFxlMTSUF3B8dfB3cn0h86qzp/Erm/axOMvFJGhDTBu8/+PudwB3ACxdutSz/b5dky0X7l2WTbKhvI70DvAfv9rEN57cTFlRIR+7bAkLGiqoKSumrLiQ7v4BDvX0c7hngJ7+V2p3ZlBgRoER1tiCWllZcQHlxUGtsSysvSVrasWFBRQVGsUFBcFrCzimtjWQcBKJoPbWP+BHa489KTW5VWueZckppx+t0fUOJI7W9IJkFfwlYytIrVGG8ZUVFVJaXEBhQcHRL3pvf4JDvf109fTT1d3Pwe4+DhzpYyABJYVGcWEBFSWF1FWUUFdRDMDug93sPtDNltZDrNt1kMdf6qWkqIDnb30LpUWFEXyar5YN+5jkZhLcCcxLeT4X2BVRLCKT5p9+so57V23n3WfP5aYrT2JmdVnUIY3o8LYimk6bFXUYr+Lu3PXkVv7xx+toP9TLcbXlUYckWSR72gbG7hlgsZktNLMS4GrgwYhjEsmoF/d1cf/qHSy7YAGfe+9rsz4BZjMzY05dkPhaO3sjjkayTc7VBN2938xuBB4CCoE73X1txGGJZNTnHt5IaVEBN166KOpQ8kJ9VSkQnFwkkirnkiCAu/8U+GnUcYhMht/s2M9Pn9/NR960mIbw4C3pmaEkKMPIxeZQkbzl7tz+8w1Mryzh2osWRh1O3qivKgGgtUvNoXIsJUGRLLKipZUnX2zjxksWUV1WHHU4eaOytIjy4kLaVBOUQZQERbLIt1duZ0Z1Kde87vioQ8k7DdUlag6VV1ESFMkSPf0DPL5pH28+tTFrrmXLJ/WVpbQdUnOoHEtJUCRLPL25nUO9A1x2ysyoQ8lLDVWl7OtUTVCOpSQokiV+tWEvZcUFXHBiQ9Sh5KUZ1SU6MUZeRUlQJAu4O8vX7+HCRQ2UFaspdDLUV5bSfqiHROJVoyxKjCkJimSBTXu72NlxhEtPbow6lLzVUFVCwqHjsGqD8golQZEssHz9HgAuPVn9gZMlOWqMTo6RVEqCIlngV+v3cvqcGmbVaozQyZIcfadVJ8dICiVBkYi1H+rl2e0dagqdZDOqg1Fj9ulaQUmhJCgSseaNe0k4ujRiktVXhs2hOkNUUigJikSseeM+GqpKOX12bdSh5LXa8mKKCkyjxsgxlARFIuTurNrSzutPrKegwEZ/gUxYQYExvbJENUE5hpKgSIR2tB9h98FuzlswLepQYqGhqlQ1QTmGkqBIhFZtbQfg3IXTI44kHhqqlQTlWEqCIhF6Zks7teXFLJlZHXUosdBQqaHT5FhKgiIRemZrO0vnT1N/4BRJ1gTdNXSaBJQERSKyr7OHza2H1BQ6hRqqSujpT3CodyDqUCRLKAmKRGR1sj9wgZLgVEleK6hRYyRJSVAkIqu2tlNWXMAZc3R94FRpqE6OH6okKAElQZGIPLO1nbPmTaOkSF/DqVJfGQ6d1qmTYySgb59IBDq7+1i366D6A6fYjLAmqMskJElJUCQCa7Z1kHA4T/2BU2p6WBPUqDGSlJVJ0MzeY2ZrzSxhZksHzbvZzFrMbKOZXR5VjCLpeGZrO4UFxlnH10UdSqwUFxZQV1GsmqAcVRR1AMP4HfAu4GupE83sVOBq4DRgNrDczJa4u853lpzyzNYOTp9dQ2Vptn4F81dDValOjJGjsrIm6O7r3X3jELOuAu5z9x533wK0AOdNbXQi6entT/CbHfs5Z76aQqPQUFVCq06MkVBWJsERzAF2pDzfGU4TyRlrdx2gpz/BUg2aHYl6DaItKSJrizGz5cCsIWbd4u4/Gu5lQ0wbcvwjM7sOuA6gsbGR5ubmiYSZN7q6umJfBuMxmeX10NY+AHpeWk9z21ANHrkpV/axnv097N7fH3msuVJe+S6yJOjul03gZTuBeSnP5wK7hln/HcAdAEuXLvWmpqYJbC5/NDc3E/cyGI/JLK/v3rOGudMO8M4rLp2U9UclV/ax5wc2sXz7C7z+wosoLSqMLI5cKa98l2vNoQ8CV5tZqZktBBYDqyKOSWTM3J3V2zpYOl9NoVGprwquFWw/pH5BydIkaGbvNLOdwOuB/zWzhwDcfS1wP7AO+Dlwg84MlVyys+MI+zp7OEfXB0ZmWkUxAAeO9EUciWSDrDw/290fAB4YZt5twG1TG5FIZqzeFgyarZpgdGrDJLj/sJKgZGlNUCRfrd7aQXVpEUsadRPdqNSVB6PGKAkKKAmKTKk12zo48/g6CnUT3cjUHa0Jqk9QlARFpszB7j427ulkqS6Sj9TRJKg+QUFJUGTK/Hr7ftzhHPUHRqq8uJCSogI1hwqgJCgyZdZsbafA4EwNmh0pM6OuvJgDR9QcKkqCIlNmzfYOTjmuhioNmh25uopi1QQFUBIUmRL9Awl+vX2/Lo3IEnXlJUqCAigJikyJ9S93crh3QBfJZ4naimKdGCOAkqDIlEheJH+u7hyRFerKi3WJhABKgiJTYvXWDubUlXNcbXnUoQjqE5RXKAmKTLJg0Ox23T8wi9RVlHCkb4DuPg09HHdKgiKTbGfHEfYc7NFJMVkkecH8QfULxp6SoMgkS/YHnqORYrLG0fFDlQRjL60kaGZFZmbh43lm9vtmdlZmQhPJD8lBs0+apUGzs0Wd7iQhoQknQTO7FtgLbAsf/xL4feA+M7spQ/GJ5Lw12zo4a/40DZqdRWrLNYi2BNIZuuKjwIlANbAemO/urWZWATwD3J6B+ERy2oEjwaDZ/+eM46IORVJoEG1JSicJ9rp7B9BhZi3u3grg7ofNTD+vRIBnt3cEg2brzNCsUleRvKegDlVxl04SLA/7/wqAkvCxhX9lmQhOJNet3tpOYYFx5jwNmp1NKksKKSow9QlKWklwN/D5IR4nn4vE3uqtHZw+u4aKEg2anU3MjLqKEjWHysSToLs3ZTAOkbzT25/guR37ueb8+VGHIkOoqyjmgGqCsTfhJGhm7xppvrv/YKLrFskHz7+0n57+BOct1PWB2aiuvJj9uqdg7KXTRvP28P9M4ALgV+HzS4BmQElQYm3Vlg5Ag2Znq7qKYl4+0B11GBKxdJpDPwhgZj8BTnX3l8PnxwFfykx4Irnrma3tnDijkvqq0qhDkSHUlpew/uXOqMOQiGVi2LQFyQQY2gMsycB6RXJWIuGs3tquptAsFtxJQs2hcZeJU9aazewh4F7AgauBRzKwXpGctXFPJwe7+zlXN9HNWnXlxRzqHaC3P0FJkYZRjqu0P3l3vxH4KvBa4EzgDnf/cDrrNLN/NbMNZvZbM3vAzOpS5t1sZi1mttHMLk8vepHJsWpL8ia6SoLZqq4yuGD+gC6TiLWM/Pxx9wfc/WPh3wMZWOUvgNPd/TXAC8DNAGZ2KkFN8zTgCuDLZlaYge2JZNSqre3Mri1j7jTdRDdb1YXjhx7QGaKxlpVtAO7+sLv3h09XAnPDx1cB97l7j7tvAVqA86KIUWQ47s4zW9o5d+F0wpusSBbSnSQEMtMnONn+GPhu+HgOQVJM2hlOexUzuw64DqCxsZHm5uZJDDH7dXV1xb4MxiOd8tp7OMHezh5qe1tjVea5to9tPRDcVf7xVc/StXXqD4W5Vl75Kp2L5e8AfgYsd/dxn2dsZsuBWUPMusXdfxQucwvQD3w7+bIhlveh1u/udwB3ACxdutSbmprGG2JeaW5uJu5lMB7plNf3Vu8AfssfXv46ljTG5x6CubaP7Wg/zK1PPcK8E0+m6Zy5o78gw3KtvPJVOj9/7iTol/vL8K4RDwM/d/ffjOXF7n7ZSPPN7APA24A3uXsy0e0E5qUsNhfYNd7ARSbTqi3t1FUUs2hGVdShyAhqK3RPQUmjT9DdV7r7re5+EfBeYDvwV2b2azO708zeO9F1m9kVwE3AO9z9cMqsB4GrzazUzBYCi4FVE92OyGR4Zms7S+dPp0A30c1q1aVFFOpOErGXkYZwd28juE7wXgAzO4egljhRXwRKgV+EJxasdPfr3X2tmd0PrCNoJr3B3QfSCl4kg/Yc7GZr22ENmp0DzEzjh8rknBjj7muANWm8ftEI824DbpvoukUm08rNbQC87oT6iCORsaitKFZNMOay8hIJkVz19JZ2qkuLOHV2TdShyBjUlRfrYvmYUxIUyaCVm9tYumAaheoPzAl1FSWqCcackqBIhuzt7GbzvkOcr6bQnKE+QUnnOsEtDHONXnKRcP6/uft/THQ7IrkiOV6o+gNzh/oEJZ37CS7MZCAiue7pze1UlhRyuvoDc0ZdeQmd3f30DyQoKlTDWBzpUxfJkJWb2zhnwXQdTHNIcvxQnRwTX+k0h3YydHOoAe7u+jkssdHW1cOmvV288+whh7KVLJVMgh2H+6ivKo04GolCOs2hRwdFNLNfu/tZmQlJJPck+wPPX6j+wFwyPbynYIeGToutTLXbjHSCjEjeW7m5jfLiQl4ztzbqUGQc6iuD2l9bV0/EkUhU1HkhkgFPb2nnnPnTKFZ/YE6prwpqgm2HVBOMq3T6BN+V8rRu0HPc/QcTjkokh7Qf6mXD7k7+6s3HRR2KjNO0ijAJdikJxlU6Y4e+PeXxo4OeO6AkKLGwakswXujrT1R/YK4pKSqgpqyIdtUEYyudE2M+mMlARHLVUy8m+wProg5FJqC+qpRW9QnGVjrNocePcdH97n5wotsRyXYrN7ezdME0SorUH5iL6itLVBOMsXSaQ/97DMs4cBfwrTS2I5K12rp62Link3ecOTvqUGSCpleWsK3t8OgLSl5Kpzn0kkwGIpKLntZ4oTmvvqqUZ7fvjzoMiYjab0TSsHJzGxUluj4wl9VXltBxuJdEQpc7x5GSoEgannqxjaULpuv6wBxWX1XCQMI1fmhM6ZsrMkGt4XihrzthetShSBqSQ6fpgvl4ykgSNLPpZjYtE+sSyRUrN4fXB6o/MKc1VGnotDibcBI0s+PN7D4z2wc8DTxjZnvDaQsyFaBItlq5uS24f+Ac9QfmsmRNUJdJxFM6NcHvAg8As9x9sbsvAo4Dfgjcl4ngRLKZ+gPzQ3L80FYlwVhK59vb4O7fdfeB5AR3H3D3+wC1D0le29vZzYv7DunSiDyQHD+0XeOHxlI6F8uvMbMvE1w0vyOcNg/4APDrdAMTyWZPvRj0B16g8UJzXnFhAXUVxbQdUp9gHKVTE/wj4HngH4GHgIeBW4HfAe9PJygz+5SZ/dbMnjOzh81sdsq8m82sxcw2mtnl6WxHZKKeerGN6rIiTptdE3UokgHTK0t0dmhMpTNiTC/wlfAv0/7V3f8OwMw+Avw9cL2ZnQpcDZwGzAaWm9mS1CZZkanw1OY2zl9YT5H6A/NCQ2Wpzg6NqUn5BpvZ36fz+kEDblfyyp3rrwLuc/ced98CtADnpbMtkfF6af8RtrUdVlNoHpmuQbRja7J+xv5puisws9vMbAdwDUFNEGAOr/Q/AuwMp4lMmaP9gaFX4HUAABPnSURBVIuUBPNFfVWJbqwbU+ncSmm42yMZUD6G1y8HZg0x6xZ3/5G73wLcYmY3AzcC/xCue7AhB/wzs+uA6wAaGxtpbm4eLaS81tXVFfsyGI+RyuuB3/ZQXQy71q9h94ahdsl4yuV9rLO1l/ZDffzqkUcosKn5THO5vPJJOmeH7gfOdfc9g2eENbgRuftlY9zOd4D/JUiCOwnOQE2aC+waZv13AHcALF261Juamsa4ufzU3NxM3MtgPIYrL3fn5qd+xRtPnsGll5w99YFlsVzex7aVbOXBF9fy2nMvoD4cQWay5XJ55ZN0mkO/BcwfZt530lgvZrY45ek7gA3h4weBq82s1MwWAouBVelsS2Q8trUd5uUD3bxe/YF5RaPGxFc6Z4d+coR5N010vaHPmNlJQALYBlwfrnetmd0PrAP6gRt0ZqhMpSd1fWBeOjpqTFcvixsjDkamVDrNoQCY2VBtQgeAbe7eP5F1uvu7R5h3G3DbRNYrkq4nX2ylsaaUhQ2VUYciGVRfGTSBqiYYP2knQeDLwNnAbwlOXDk9fFxvZte7+8MZ2IZI5NydlZvbuGjxDGyKTp6QqZGsCWrUmPjJxCUSW4Gz3H2pu58DnEUwasxlwGczsH6RrLBpbxetXb26dVIemlZRghm6TCKGMpEET3b3tckn7r6OICluzsC6RbLGik2tALxhcUPEkUimFRYY0ypKVBOMoUw0h240s6/wyu2T3ge8YGalQF8G1i+SFZ58sZWFDZXMqRv1MljJQRo1Jp4yURNcRjB82UeBjwGbw2l9wCUZWL9I5PoGEqzc3K6zQvNYfWUJrWoOjZ20a4LufsTM/pPgLhIObHT3ZA2wK931i2SD3+7cT1dPPxcuUlNovqqvKuGFPTpkxU0mLpFoIrin4FaCs0PnmdkH3P2xdNctki2eaGnDDF0kn8emV5boThIxlIk+wc8Bb3H3jQBmtgS4FzgnA+sWyQorWlo5fXYtdeFdyCX/1FeWsv9IH/0DCd0iK0Yy8UkXJxMggLu/ABRnYL0iWeFwbz+/3t6hu0bkufqqEtyh47DO54uTTNQEV5vZN4G7w+fXAGsysF6RrLBqSzt9A67+wDyXOmrMjOqpGURbopeJmuCHgLXAR4C/IBjX8/oMrFckKzzR0kpJYQFL50+POhSZRK+MH6p+wTjJxNmhPcDnwz+RvPNESxvnzJ9GeUlh1KHIJErW/pQE4yWdm+o+zzA3tAVw99dMdN0i2aKtq4d1Lx/kby4/KepQZJIlk+C+TiXBOEmnJvi2jEUhkqWe0K2TYqO6tIjSogIlwZhJ536C2zIZiEg2emJTKzVlRbxmbl3UocgkMzNmVJcqCcaMLoYRGYa7s6KllQtObKCwQLdOioMZ1aXsU59grCgJigxjS+shXtp/hAt114jYmFGlmmDcKAmKDGNFS3DrpIuUBGNDzaHxM+EkaGa1ZvYZM9tgZm3h3/pwmjpQJOc9vqmVudPKOX56RdShyBSZUV1K++Fe+gYSUYciUySdmuD9QAfQ5O717l5PcOukDuB7mQhOJCoDCWfli21ctLgBM/UHxkVDVSnu6L6CMZJOElzg7re7++7kBHff7e63A8enH5pIdDYfSNDZ08+Fi2ZEHYpMIV0rGD/pJMFtZvZxM2tMTjCzRjO7CdiRfmgi0VnbNoCZrg+MGyXB+EknCb4PqAceNbN2M2sHmoHpwHszEJtIZNa2DnDGnFqmVerWSXEyo0pJMG7SuVi+A7gp/BPJG53dfbx4IMH1Z+qs0Lg5WhPUtYKxkdYlEmZ2spm9ycwqB02/Ir2wRKKzcnM7CUfXB8ZQWXEh1WVFqgnGSDqXSHwE+BHwYWCtmV2VMvtf0g0s3MZfm5mbWUPKtJvNrMXMNprZ5ZnYjkiqJ1paKSmAc+ZPizoUiYCuFYyXdAbQvhY4x927zGwB8H0zW+Du/w6kfU65mc0D3gxsT5l2KnA1cBowG1huZkvcfSDd7YkkPb5pH0umF1JapFsnxdGMKg2dFifpNIcWunsXgLtvBZqAK83s82QgCQJfAD7Osbdrugq4z9173H0L0AKcl4FtiQDw8oEjvLjvEKfVKwHG1YzqUlpVE4yNdGqCu83sTHd/DiCsEb4NuBM4I52gzOwdwEvu/ptBFyrPAVamPN8ZThtqHdcB1wE0NjbS3NycTkg5r6urK/ZlMBaP7+wD4ISKHpXXOOXLPtZzoIeX9/dP+nvJl/LKdekkwT8C+lMnuHs/8Edm9rXRXmxmy4FZQ8y6BfgE8JahXjbEtCFv7OvudwB3ACxdutSbmppGCymvNTc3E/cyGIsf3PtrGqraWDKzUOU1Tvmyj62jhV9s28j5F1xEecnktQjkS3nlunQukdg5wrwnxvD6y4aabmZnAAuBZC1wLvCsmZ1HUPObl7L4XGDXOMIWGVYi4TzR0hoOlXYg6nAkIslrBVu7epincWPzXtbdRcLdn3f3me6+wN0XECS+s8Ph2R4ErjazUjNbCCwGVkUYruSRDbs7aTvUy4WLNVRanCWvFdyrfsFYyLokOBJ3X0swcPc64OfADTozVDJlRcs+AC5cpOsD40xDp8VLOn2CUyKsDaY+vw24LZpoJJ89vqmVxTOrmFVbxoaog5HIaNSYeMmpmqDIZOnuG2DVlnaNEiPUV5ZSYKoJxoWSoAiwZlsHPf0JNYUKhQXG9EqNGhMXSoIiwIqWVooKjPNP0K2TREOnxYmSoAiwYlMrZx8/jarSrO8mlykwo1pDp8WFkqDEXsehXn636wBvUFOohGZUaei0uFASlNh78sU2XLdOkhQN1SXs6+zBfcgBqSSPKAlK7K1oaaW6tIjXzq2NOhTJEjOqSukdSHDwSP/oC0tOUxKU2FvRso/XnVhPUaG+DhJ45VrB7ogjkcmmb73E2va2w+xoP6JLI+QYGjotPpQEJdYeTw6Vpv5ASTFTQ6fFhpKgxNoTLa0cV1vGCQ2VUYciWWRmTRkAew8qCeY7JUGJrYGE80RLGxcuamDQzZsl5qpLiygvLmTPQfUJ5jslQYmttbsOcOBIn5pC5VXMjJk1peoTjAElQYmtxze1AnDBiUqC8mqN1WWqCcaAkqDE1opNrZw8q/romYAiqWbUaPzQOFASlFg60jvAmm0dXKSmUBmGaoLxoCQosfT0ljZ6BxJcuHhG1KFIlppZU8qh3gG6ejRqTD5TEpRYWrGplZLCAs5bMD3qUCRLNdaEF8yrNpjXlAQllla0tLJ0wTTKSwqjDkWyVGN1cK3gHl0rmNeUBCV29nZ2s2F3py6NkBHNTNYEO1UTzGdKghI7T7QEl0ZctEj9gTI8jRoTD0qCEjuPb2plWkUxp82uiToUyWLVpUWUFRfoDNE8pyQoseLurNjUyhsWNVBQoKHSZHhmRmNNmUaNyXNKghIrm/Z2sbezR9cHypjMrC5VTTDPZWUSNLNbzewlM3su/HtryrybzazFzDaa2eVRxim5JzlUmq4PlLGYWVOmUWPyXFHUAYzgC+7+/1InmNmpwNXAacBsYLmZLXH3gSgClNzz+KZ9nNBQyZy68qhDkRzQWF1G88G9UYchkygra4IjuAq4z9173H0L0AKcF3FMkiO6+wZYubmNNy5RLVDGRqPG5L9sToI3mtlvzexOM5sWTpsD7EhZZmc4TWRUq7d20N2X4I1L1B8oY6NRY/JfZM2hZrYcmDXErFuArwCfAjz8/zngj4GhTufzYdZ/HXAdQGNjI83NzekHncO6urpiXwbf3dhLoUHvznU0714/4rIqr/HLxzJ7uS3oaXnosac5pT6zowvlY3nlosiSoLtfNpblzOzrwE/CpzuBeSmz5wK7hln/HcAdAEuXLvWmpqYJx5oPmpubiXsZfOa5xzhvYQ1XXPa6UZdVeY1fPpbZ3L2dfPaZx5h94sk0nZnZRqd8LK9clJXNoWZ2XMrTdwK/Cx8/CFxtZqVmthBYDKya6vgk9+w9GAyVpv5AGY8Z1Ro1Jt9l69mhnzWzMwmaOrcCfwbg7mvN7H5gHdAP3KAzQ2UskpdG6PpAGY+asmDUGI0fmr+yMgm6+/tHmHcbcNsUhiN54LFN+2ioKuHU4zRUmoxdctQY3Ukif2Vlc6hIJiUSzuObWrlo8QwNlSbjNrO6VDXBPKYkKHlv7a6DtB/qVVOoTMjMmjL1CeYxJUHJe49t2gfARRoqTSZA44fmNyVByXuPvrCPU4+rYUZ1adShSA5qrCnTqDF5TElQ8trB7j7WbOug6STVAmViZlZr1Jh8piQoee2JTa0MJJymk2ZGHYrkqMbkHeZ1N4m8pCQoea154z6qy4o4+/i6qEORHJUcP1T9gvlJSVDylrvz6Av7uGhxA0WF2tVlYmbWaNSYfKYjg+StDbs72X2wm4s1VJqkobq0iPLiQnarJpiXlAQlbzVvDC6NuHiJ+gNl4syMWbVlag7NU0qCkrcefWEvJ8+qZlZtWdShSI6bWV2q5tA8pSQoeamzu4/VWzt0VqhkRGNNmZpD85SSoOSlJ1ra6E+4rg+UjEg2h7oPeQ9vyWFKgpKXHn1hL9WlRZwzf1rUoUgemFldSk9/ggNH+qIORTJMSVDyjrvzyIZ9XLi4gWJdGiEZkOxX1i2V8o+OEJJ31r8cXBpxycnqD5TMSI4ao37B/KMkKHnnkY17AdQfKBkzqyZZE1QSzDdKgpJ3HtmwlzPm1DKzWpdGSGYk70Cy54CSYL5REpS80nGol2e3d6gpVDKqrLiQaRXF7NEd5vOOkqDklcc27SPhcImaQiXDGmvK2H1AJ8bkGyVBySuPbNhLfWUJr52ru0ZIZjXWlLFXNcG8oyQoeWMgEdw14uKTZlBQYFGHI3mmsaaU3eoTzDtKgpI3ntvRQcfhPi7RUGkyCWbVlNHa1UP/QCLqUCSDlAQlbzyyYR+FBcYbdeskmQQza8pIOLR29UYdimSQkqDkjV9t2Ms586dRW14cdSiSh3StYH7K2iRoZh82s41mttbMPpsy/WYzawnnXR5ljJI9dh/oZt3LB7lUl0bIJNGoMfmpKOoAhmJmlwBXAa9x9x4zmxlOPxW4GjgNmA0sN7Ml7j4QXbSSDZKjxKg/UCZLY21wwfxeJcG8kq01wQ8Bn3H3HgB33xtOvwq4z9173H0L0AKcF1GMkkUe2bCXOXXlLGmsijoUyVP1laUUFphqgnnGsvH+WGb2HPAj4AqgG/hrd3/GzL4IrHT3e8Llvgn8zN2/P8Q6rgOuC5+eBGwcRwi1wIEJhD6W1420zHDzxjp9pOcNQOsosU1UVOU10vyhpqu8Jq+8YPLKTOU1fhMps6FeM9/d8/tMM3eP5A9YDvxuiL+rwv//ARhBTW9L+PhLwB+mrOObwLsnIbY7Jut1Iy0z3LyxTh/pObB6Ej/LSMprvGWm8pq88prMMlN5TU2ZTbScc/0vsj5Bd79suHlm9iHgBx58MqvMLEHwq2knMC9l0bnArkkI78eT+LqRlhlu3linj/Z8skRVXiPNH2q6ykvlNd5lcrG8JrqtqYwva2Rrc+j1wGx3/3szWwL8EjgeOBX4DkHtcHY4fbHrxJhRmdlqd18adRy5QuU1fiqz8VF5ZYesPDsUuBO408x+B/QCHwhrhWvN7H5gHdAP3KAEOGZ3RB1AjlF5jZ/KbHxUXlkgK2uCIiIiUyFbL5EQERGZdEqCIiISW0qCIiISW0qCMWVmJ5jZN83sVQMNyKuZ2e+Z2dfN7Edm9pao48l2ZnaKmX3VzL4fXvIkozCzSjNbY2ZvizqWOFESzCNmdqeZ7Q3Pqk2dfkU44HiLmf0tgLtvdvc/iSbS7DDO8vqhu18LLAPeF0G4kRtnea139+uB9wKxvAxgPOUVugm4f2qjFCXB/HIXwVBzR5lZIcFIO1cSXGf5B+FA5DKx8vpkOD+O7mIc5WVm7wBWEFzPG0d3McbyMrPLCC792jPVQcadkmAecffHgPZBk88DWsKaXy9wH8HQdLE3nvKywO0EY9U+O9WxZoPx7l/u/qC7XwBcM7WRZodxltclwOuA/wtca2Y6Nk+RbL1YXjJnDrAj5flO4HwzqwduA84ys5vd/dORRJd9hiwv4MPAZUCtmS1y969GEVwWGm7/agLeBZQCP40grmw1ZHm5+40AZrYMaHX3RASxxZKSYP6zIaa5u7cB1091MDlguPL6D4JB3eVYw5VXM9A8taHkhCHL6+gD97umLhQBNYfGwVQNOp4vVF7jo/IaH5VXllESzH/PAIvNbKGZlQBXAw9GHFM2U3mNj8prfFReWUZJMI+Y2b3AU8BJZrbTzP7E3fuBG4GHgPXA/e6+Nso4s4XKa3xUXuOj8soNGkBbRERiSzVBERGJLSVBERGJLSVBERGJLSVBERGJLSVBERGJLSVBERGJLSVBkQwzs4+Y2Xoz+/Yw89eEF0qnTms2s6UpzxcMvgWPiGSexg4Vybw/B6509y2DZ5jZAuCl8A4CIhIx1QRFMsjMvgqcADxoZh8bYpErgZ+Pc53fMLPnwr99ZvYPmYhVRDRijEjGmdlWYKm7tw4x70fAx9x986DpzcBxwJFwUgmQcPfTU5aZTzDc1uXuvm1yoheJFzWHikyRsB9w7uAEmOIad18dLrsA+EnKa8uA7wE3KgGKZI6aQ0WmzkXAigm+9qvAD9x9eQbjEYk9JUGRqXMF8LPxvsjMbgCq3f0zmQ9JJN6UBEWmThPw6ARe99fAGSknx1yf2bBE4ksnxohMATObC3zd3a+MOhYReYWSoIiIxJaaQ0VEJLaUBEVEJLaUBEVEJLaUBEVEJLaUBEVEJLaUBEVEJLaUBEVEJLb+P0EmGBB5p5YyAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "fs = 44100  # sample rate\n",
    "# frequency axis\n",
    "f = np.logspace(-4, 0, 100)*fs/2\n",
    "# dirac as input signal\n",
    "x = np.zeros(2**10)\n",
    "x[0] = 1\n",
    "\n",
    "h = tools.blackbox(x, fs)  # impulse response of black system\n",
    "H = naive_ft(h, f, fs)  # transfer function of blackboard system\n",
    "\n",
    "plt.semilogx(f, 20*np.log10(abs(H)))\n",
    "plt.ylim([-60, 5])\n",
    "plt.xlabel('f / Hz')\n",
    "plt.ylabel('20 log10 |H| / dB')\n",
    "plt.title('DTFT of the blackbox system''s impulse response = Frequency Response')\n",
    "plt.grid(True)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Solutions\n",
    "\n",
    "If you had problems solving some of the exercises, don't despair!\n",
    "Have a look at the [example solutions](linear_systems-solutions.ipynb)."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<p xmlns:dct=\"http://purl.org/dc/terms/\">\n",
    "  <a rel=\"license\"\n",
    "     href=\"http://creativecommons.org/publicdomain/zero/1.0/\">\n",
    "    <img src=\"http://i.creativecommons.org/p/zero/1.0/88x31.png\" style=\"border-style: none;\" alt=\"CC0\" />\n",
    "  </a>\n",
    "  <br />\n",
    "  To the extent possible under law,\n",
    "  <span rel=\"dct:publisher\" resource=\"[_:publisher]\">the person who associated CC0</span>\n",
    "  with this work has waived all copyright and related or neighboring\n",
    "  rights to this work.\n",
    "</p>"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "mystiasp",
   "language": "python",
   "name": "mystiasp"
  },
  "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.8"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 1
}