{ "cells": [ { "cell_type": "markdown", "id": "32519996", "metadata": {}, "source": [ "# Les verres musicaux\n", "\n", "Un verre d'eau partiellement rempli peut émettre une note distincte lorsque l'on frotte son bord avec un doigt mouillé. C'est l'ingrédient de base du [Verrillon](https://fr.wikipedia.org/wiki/Verrillon). Mais que cache son spectre ?\n", "\n", "La vidéo {download}`data/Singing_Glass.MOV` à été faite avec un verre à vin, de l'eau et un smartphone. Sa bande son à été extraite dans le fichier {download}`data/SingingGlass.wav`.\n", "\n", "1. Télécharger, puis importer la piste sonore SingingGlass.wav :" ] }, { "cell_type": "code", "execution_count": 1, "id": "fca3e2aa", "metadata": {}, "outputs": [], "source": [ "from scipy.io.wavfile import read\n", "samplerate, amplitude = read('data/SingingGlass.wav')" ] }, { "cell_type": "markdown", "id": "7ceae727", "metadata": {}, "source": [ "2. Afficher le nombre de points N contenus dans ce signal, la durée de l'intervalle entre deux points et la durée totale de l'enregistrement.\n", "\n", "3. Afficher le signal. Quelle est approximativement la fréquence du son ?\n", "\n", "4. L'intensité du spectre à une fréquence donnée est proportionnel au carré du module de la transformée de Fourier. Affichez le spectre de SingingGlass.wav.\n", "\n", "5. Zoomer autour de la fréquence du son. Que remarquez-vous ? \n", "\n", "6. On a phénomène de battement que l'on entend à l'oreille. Qu'elle est la fréquence du battement ?\n", "\n", "\n", "Il s'agit d'un phénomène méconnu mais relativement commun avec les verres musicaux, où deux modes de vibration orthogonaux sont excités. Ces deux modes n'ont pas exactement la même énergie due aux imperfections dans la symétrie du verre. La non-dégénérescence de ces modes dépend du niveau d'eau. Voir par exemple:\n", "https://newt.phys.unsw.edu.au/music/people/publications/Jundtetal2006.pdf" ] }, { "cell_type": "markdown", "id": "501f2f54", "metadata": { "tags": [ "sol" ] }, "source": [ "--------------------------------------------------------------------------------------\n", "\n", "## Solution" ] }, { "cell_type": "code", "execution_count": 2, "id": "7d365267", "metadata": {}, "outputs": [], "source": [ "from scipy.io.wavfile import read\n", "from matplotlib.pyplot import figure\n", "\n", "import numpy as np\n", "from numpy.fft import rfft, rfftfreq" ] }, { "cell_type": "code", "execution_count": 3, "id": "f53ef5d7", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "N = 800001\n", "dt = 2.2675736961451248e-05 s\n", "T = 18.14061224489796 s\n" ] } ], "source": [ "Fs , amplitude = read('data/SingingGlass.wav')\n", "amplitude = amplitude/2**16\n", "\n", "N = len(amplitude)\n", "print('N =',N)\n", "dt = 1/Fs\n", "print('dt = ', dt,'s')\n", "T = dt*N\n", "print('T =',T,'s')" ] }, { "cell_type": "code", "execution_count": 4, "id": "f85380d4", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAD4CAYAAADvsV2wAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAAsTAAALEwEAmpwYAABy6klEQVR4nO2dd3xUVfr/P2dqeichJEDAhN4JRUEMTRAUYXUVdVdsi676dXf97Squiqu7rrhf1/W7q1uwgWXFiiC9hiYt9B5aIAkhpGcyyfTz++POnUy5bVomIef9evFi5t577nnmzuQ5z3me5zyHUErBYDAYjOsfVaQFYDAYDEbbwBQ+g8FgdBKYwmcwGIxOAlP4DAaD0UlgCp/BYDA6CZpICyBGWloazcnJCbi90WhEbGxs6AQKE0zO0NIR5OwIMgJMzlDTVnIeOHCgmlLaRfAkpbRd/hs5ciQNhq1btwbVvq1gcoaWjiBnR5CRUiZnqGkrOQEUURG9ylw6DAaD0UlgCp/BYDA6CUzhMxgMRieBKXwGg8EIgE/3XELOgtUwmKyRFkUxTOEzGAxGALz8/XEAwDWDOcKSKIcpfAaDwQiCR5fsj7QIimEKn8FgMIKgpKY50iIohil8RlA0mW2w2ByRFoPBaFNqmjqOG8eddrvSltExGPTKegDA2ddvg1bN7AdG56DWaIm0CAHB/kIZAXOt0eR6/cAHeyMoCYPRtlQ2dkwLnyl8BTRbbJEWoV3y8w/3uV7vu1gbQUkYjLbl4OU60XOT/lqIP6851YbSKIcpfBkKz1zDgIXrMfu9XZEWpd1xptIQaREYjIhgNAsbgUt2XcSFKiMWb7/QxhIpgyl8GR76mEu5OlxaH1lBGAxGu+E/XgqdOvcG/8MPJ13HyuraX/YOU/gMBqNd8frqk+1SWUpRIxDE3XG2OgKSSMMUfgeFUoqcBauRs2A1bHaWFsnoeDRbbLhv8R6P3+/+klq8v+Mixr+5NYKS+U95XYvPsRe+OxYBSaRhCt8P2lO++aJ1p12vc19cG0FJWnE4aKRFYISBhz/ehz+uOil/oZ8MWLgeuy/UePx+f/rv3SHvpy2QCuK2J5jCl8Dbct5f0n4yUf6zrf0FhRyUKfzrka1nqvDhzouRFqNds/7E1UiLoAim8CXYVlzl8f6HI1ciJElkKak24qF1RuQsWC153aZTlSHv+7GlRcj9/Rqf449/WoSZf98R8v4Y4pxt46wss83epv0Fw54L7ccYlIIpfAkeXVrk8X7Z/tIISRJZfvfNEddrk1X8j/DTPZdC3vemU5WwOSjWHKvwOL7+RCVOXGnE6auNIe+T0cr/+6r1u5/6t+1t2veBkvbhJtl8qhI5C1bjwKX2IU8wMIXfiRnxx4249z/yPtMfDxyDzVANSilecpaEFVL8u87VhFQ+91znJz8/KHjN9HeYlR8KTlU04ot9l32Of3uwLALScCze0T7clrzhd9e/foywJMHDFH4H4JEl+/HOpmLJawIp5lRrtGCvghWyDrMRjXu/hWH/9/j0q+W4ePEi9l0IrXIXoqTGKHj8kshxRuDc9n878MJ3x5CzYDXOVzWJXnfmamjcOhUNLZLvAaDwTJXPMX/Zevoachashp0lFABgCr/dU2UwY8vpa3hn01nX4g4htp/174/Dn+mpvltfqKLiEDtoEoguGhs3bsTH//47ms/ugbWmDNQRHl/rzL/vFDz+2g+hzxhhtDL5r9tEz/2z8FxI+vi6yHPmsPJweOJjDztr1a873jGCquEmJAqfEDKdEHKGEHKOELJA4PwEQshBQoiNEHJ3KPrsLIx6fZPr9c5z3EKOagFr/rlvjvp130Cmp+qYRER1H4T58+fDnDsR6pgkmMtPoXHfdzCe3glLVQnM5vAXlbJ5WWssHbTtWBEixfzjec9FSVVh3jXq0z0lAbU7d81zRtPR62oFrfAJIWoA7wG4DcAAAPcRQgZ4XXYZwEMA/htsf50Jb4ueL1Z2+HK9z7VWe3iU3qqjwn/gJ6rt0Gf1Q9yQqYgfcQc0SV1huXYRb7/9Nj7//HMcPHgQRmN4XC/e2VPey9wZ/nHyStsHvhtbPBXnB0Gkfa65YJHNIAs0i6ak2nPFr1RVWCHDo70tigyFhT8awDlK6QVKqQXAMgB3ul9AKS2hlB4F0L4+fTvn7Y3Cfvtj5Q0h7Ucq/W310QrB4w0trRs3q3RR0HfNRdzAiXj22WcxdOhQnD9/Hv/4xz/w8ccfY8+ePWhoCK3M7kQysHg9cI+CwH2oOVnhO8gEauV/Vcz9Fr2VvlhGGb9C/ViZ/G9yrZcr6JCAseXqT+DvyN7O1qaEYgOULADu+YplAMaE4L6dnn9sEfaXKk1//PZAGaJ1aswYnCl53eJtF/A/k/MEz3n/4OXQ6/UYNGgQBg0aBKvVigsXLuD06dPYvn07kpKS0K9fP/Tv3x9dunTx674AF5hOjdP7HD93TTzIyJCnSaTyY7BQStHrhTW4Y2g3/OO+4bLXhyogzLNTppbNHe/uRMmimZLXSMXNfK9VfGnEaFc7XhFC5gOYDwAZGRkoLCwM+F5NTU1BtRdj85atUKsIAOBgpQ1/P8RZJUumxwZ0P3/lLCwsFN1tx/0+D61rdaf8eXw0usWJT+b+urEYg9Xlfsngz7nExETk5+ejsrISP/74I5YtWwatVosePXqgR48eSEtLAyFEtt/1hbtEPwffb7i+91DSEWQsLCxEU1MTAN/vRansF+o5i/eHI1cwPqEOGbHSDoW1uw753ZfZ5qllt2zdCpXzt/T3AyaPc4WFhT6rweXuv/yQr1tSrM26Lb7rFDZu2Y44HSdPe/jeQ6HwywF0d3uf7TzmN5TSxQAWA0B+fj4tKCgIWKjCwkL40768vgWpsTpEadWtB9f5+gVV3QaioF86HA6Kh9xWgD683oiLb0hbC37LKdB/QUGB4HEAGD5mHBKjtT5tf7+zxcOSoZQC6zxXr/ojw9BRNwHrNgpe7n6fFosdJqsdybE6j2sopbhy5QpOnTqFU6dOoayszGX59+jRAyqVSrBvVXouCkb3EH8u8P97jwTtTkaB53nLLbdgxfqtAHzTJZXKvmXFcQDcbPT5HZ6/QaE+7fEZ8HQWyPe16WQlgNYFkoNG3oj0hCgAwMPrPfsoKChAfbMFWL/R45gU1I+/wdju/QF4rhdJyBmECX242Wx7+N5D4cPfDyCPENKLEKIDMBfAyhDct824Ut+CcYu2oN/L62Sv5dO8lu4u8TjeVtM5qen3SoWlH+qbrfIXSVAmUBlQiP4L12H4Hzdiu1eQlRCCrKwsTJkyBU8//TR+9rOfITY2FuvXr8dbb72FFStWYNOeQz7pnpdqQ1cyd8Xhcnx/SNwuqW4yI2fBarRYOs7y/lByvsqIBov8j9potonm7X+y27+V10Ir2eWCnu5FBAHgB7eYk/ffpNlmxwc7PIPD7tt0BouQDthy+lrI7h8Kglb4lFIbgKcBrAdwCsBXlNIThJDXCCGzAIAQMooQUgbgpwD+Qwg5EWy/oeSmRVt8jkn90Mw2e8hWlVYaHX4VZdsgUaTpZecq2AYZhX64rF5xf0Isl1CUl2s4pez+/B78aJ/Y5SCEoEuXLpgwYQIef/xxzJ8/HxkZGVj4n+880j2pzYp/FZ4PSm4eSil+tewwfv3lYfzPF75uBADI/xOXDtt/obwRcD2y9lgFTtYI/w2U17cO+ANfWY/Jf92GnAWrYQ1DRsoqkaQBHu/4jVRVzxNXGvFlkeegcjkAI0Lsc24UqCV14kr4khUCISR5+JTSNZTSPpTSGyilrzuPLaSUrnS+3k8pzaaUxlJKUymlA0PRbzjgSyDXSSjNvi+tEywU9uM5/zc8eH5Hi18lYZ91q20ixqurfMdT9+DTDwK51GuPSf9hufPRLvEUuq8PcH9Qf1od2J6eSUlJGDt2LK5kjEX8iNtd6Z6N+5fDeLIQhw4dgsManFXW64VWd1aoC+IdvFyHzacqw6L8AsVgsroyU5SuWWi22lFhFP4MVxuEn/+vvzwseU+51EkhQll2uKik1icT6L2t/i8kKxYpIvfdQV9DaH87qQfE0+lX2npH4euauYDooQB+aPdL5OgCwIc7L+LF5cKbIngv8BBaXKUUoR/e+251SXYIDEy/FKlV4y/8ZuZLfizxOB7I4iiVLtqV7pkwaja0XXJQXFwMQ9FKNB3bBPOVM3CYg8/1n/RWocd795RTf6CU4if//BGPLi1CXjvZowAABv9hg+t1b6/Ko2JZKN8dLMP2MmH34eurOSvau1yBWApvMPx3r299n0BZts/XZbQ1gPINywX+vjoKnV7hL/jWUwGP+fNmAEB9gH/0YjS0WPHHVSfx+d7L+NWyQzCYrB5T4/mfHvC4fnOISw3/eU2rrzOcqxrFavN4L5Z6bGmRXxYf0eig65KDsVNuR8Lon0DfrR/sTbUwHF4Lw+F1WLtpK6qruYFs40muuiHvXpLjQrXnoDH01Q0e75UOVsHGRtypMpj9mnWJIefeOyGy6KqyUfw3ctCZi/5RAIul5OTxxn1V9b8Kz6OkWn6AF6u15P09e5OzYDXu+fdu2VTMYBaJRZpOr/C9fXo8H+8qCWk/e92Kja04fAWD/7AB49xiBxeqPH+Mz38rPBOQQ6p8sRyhHmTc4YPdPLxLzN+NIxwUIGoNtKnZiOlzIxJG/wRROcNw6NwVfPLJJ1i+fDl+9toHsBlqcPNffGMzgVDb3JoGa7U7RAO5BV4zBX9yuL0Z9fom/PLzg34rSG8e/6xI8vy7Ims9lPD6Gl+3nVyNpt9/fywgd9elGiPeXHfa5xkLccv/yl/jDV9me19JrYfL73qj0yt8ISilOCWwElBpWyG8LfhwEcxeoN71/xtNoZ3lCPH4pwdc/uX/23RW9vptZzyzHghRQZvUFeqckfjNb36D8ePHAwCai3+EoWgFPlr2HUpKSuBwcEpGLKNEqgLnDrfCdHkvrkX/hetwuLTe5zpvV1Ao9gd4/lv/aiR5I1RSwD2gvi7EOzXJ1WhafbQCpQEESt2VuJLCf/66EH/1xWE/JeqYMIUvQFUQ/vO/bw5NNcFAkfL9VypIQXPft/d3X8sHiIWQ+2MTU4R/21Qsm4bn7gZz51+F50EIwerKBETnDEP8iNsRO6AAv19ZjPXr1+Ovf/0rVq5ciXe+2SZY3fOW/y0UHax/8yX3HNzPz35vl6ScALBwhXQy2ldFpchZsFpyVhZqhQyEr+6SUoLtXUnhP/dZmRLOCARiQ12Uz2ilES++1qkVvlgNmbfWnwn4nn+TqVsvxd1h3mCBj09I8cJ3ra6k9ScCc/EcKpW2wPj0USFyX1wrGTR9f4e0/7TQGWgkhEAdm4SoHoPx+OOP47HHHkOXLl2wc9dOt3TPS6C21r7kpvKbT3nOLoIt6cBXOO338joPpf/KCvHnEwqOt0GqoJSyDDQo7o2UcSBXVkEJR0NYs+qawYSnNjdjwML1IbtnIHRqhb9NJEL/VVFkinEVtYMt1L49WIb/bAsu3/2bA8FlMXgHTUNBcnIybrzxRtRlj3dL97zgSve0VJ6Hw2qSDCSv8ErhnPK2eN14ObxnOfyiv5wFq7HUzwVL/uJPGrAYcvGJPRIb5DzvZylvMU5L1N6RSxF1R2jzlZomM3YFkGYtxv/8t3W9h9xmRuGkXdXSCRf8H7F3oaRXw7SRBt/fzucnIjs5BkV+LKxqMtsQp4/s1/LG2tNIEyhSphShrfIALqDsUbqijTDb7NBrWvvl0z31XXNBbRZY667AWl2KlgsHoI5LgTa1O7Sp2VDpPesjSeXs+ztVF5rl/Eli0VCoyf39Gpz784yA2xdXSs9upOrmnw1RsbtQ7TF7sco3fvPcN0eDSo32xj177Z1NZ/HrKX1Cdm9/uO4tfHdLxNt6E/MHh4rxb26F1e7A3X5YVINeWS9aHC0UKA2Y/b8A/fdShHuTC6NI2YmzEsqJT/eM7X+zM92zL5fueYhL9zSVnYC9WXxqz+eif39IWMEpKcHLE+p0P/d4jDc2Bw0qi2jJj+Ky1hototlvAfUlsNDPZLUrLiUix18FypBvPn0NR/z47joK173C/8UnnpknbV1K97cBKM4jQZY+kOKhj8XLHISbUK6aFEJsoLz9H8JbJXrDpXt2d6Z7zkFUzjA4zM0wHt+Mv//jXbRcOgJbU62Horzlf7msqN+LLKhz3xKwvtkSUAng7wKs9y/XVzDph18ILGLiCfXA/ue1p32OXappDpmFH6r7dASue5fOJq9A25S3t8nWwA4lgWwJ9/DH++UvCpDzAtPXtuJvG4tx57CssN1/xt93hOxeRKWGNqkrtEldQXvno9vQ7sDOH9B8eidAKbSp2dCmdkcpTZO8j/t+AsNeE64wKsezXx3BT0Zk+91u9j/lM4nCwccSpTcCQWimclHBAqz2Qrhntv5w3St8RvuhROHK10AxmMKT8kYIwbqLVkT3Gs5Z/c0NsNaUouVCERyWFqxcSWGtvQpNUlcQVdvHKMTwLn3QVghVvQw1T3zWNutaQsHGk77Zbnyszmp3IO/FteidFostvy0IuyydUuHXN1ta68YzGArgVwbz6Z58yqfd1IRmVQzMZSfQXPwjNMndoE3JhjalG4g6tL+xf287j0VO90ZbzlIZwSHk7qsymBGn17hqLsmVfQgV170PX4hhr23sUFPC64lwbaenhCthCNKro+KwYKcZcUNu5dI9E9K5dM99y2E8uQ2Wygtobg7NzGaRmy87mIArI/II7cUQzmQNnk6p8AFg0l8Dz6FmBI6S4lfhoKHZGnZfqkoXDX1mXmt1z7QesNaW4Z133sH7H37srO4ZGuV/Pdd76Qz83+azPjWFRvwxsBiPP1zXLh1mBbU/PvHaKayt+GzvJcnFQKGGaHTQpfeCLr0XRt4+FOWXS2BvWgvT5WNQRcVBm9Yd2tTuUEcnKLqf0G/Z4aBQqeT3Ama0T/aLVJYNJ9e1wo90zRCGL5Faxfy/QZTLCJbH/3sE/TMTENPnJlCHHbaGa7DWlqHp6CaotHponBk/6thk0c3chab7T/33IP71s5GC1zNjp/0jt39GOLiuFb7QrlQMRiTgq68SlRra5ExokzNBe+fDbqiGtabUJ91TnZAGQlo9rq+s9C3E5p7yCQD/2XYeb6w9jYtvzEBNG/iDGR2P61rhPxmiXZwYjHBACIEmoQs0CV0QlTMcjuZ6j3RPXvlrEjMk93Y1We2uWjwA59+/fUhmW3wERoi5ZjAhPT4qbPe/rhU+g9FR4NI9k6GOTUZUjyGwtxhgrSmD6fJxOFp2cemeqd2hTc70SPesM1owXCDYVxjA1n2MyPPLzw7i21/eFLb7M4XPYLRD1NHxUGf3R1R2fzjMzbDWlsNy9Rxazu6BJjGDs/xTsgSVPRDZ9FdG4By4VAdKqWgsJ1iuW4X/yo/hLYzGYLQVKn0M9Jl50GfmwWE1w1ZXDmtNGVouFEEdn8ot9ErtDpU+JtKiMkJArxfW4Bc398KLMweE/N7XpcKnlOJSo//7ZjIY7R2VVg9dem/o0nuD2m2w1VXAWlPKpXtGx7tKOytN92S0T97fcbH9KnxCyHQA/wdADeADSukir/N6AJ8AGAmgBsC9lNKSUPQtRLjLHjMY7QGi1nD5/GndW9M9a0pd6Z6c8u8OVWxS2FwEjI5F0AqfEKIG8B6AqQDKAOwnhKyklLrv5vAogDpKaS4hZC6ANwHcG2zfYkhlNDAY1yMe6Z43jHKlexpPbwcouIyftO5Qx3umezI6F6Gw8EcDOEcpvQAAhJBlAO4E4K7w7wTwB+frbwC8SwghNAyrQ+wO6lFzhMHobPikexrrYa0tRcv5IjgsJmhTs1zpnu2puicj/IRC4WcBcK+HWgZgjNg1lFIbIaQBQCoAj00jCSHzAcwHgIyMDBQWFvotTFUz890zGDyEEKjjkqGOc0/3LHWle2qTs7iVvl7pnozIE4j+k6NdBW0ppYsBLAaA/Px8WlBQENB9frddfCNqBqMzw6V7DkBU9gBnumcZLFfPtqZ7pnWHJjkLKm3gexozQkOg+k+KUCj8cgDd3d5nO48JXVNGCNEASAQXvA0LHz2Uj0eWFMlfyGB0Yrh0zz7QZ/ZpTfes5lw/6vhULuibks3SPa8jQqHw9wPII4T0AqfY5wK43+ualQDmAdgN4G4AW8Lhv+eZ2Dc9XLdmMK5LPNM9rbDVXeVcP5eOQhWd0FrjJzo+0qIygiBohe/0yT8NYD24tMyPKKUnCCGvASiilK4E8CGATwkh5wDUghsUwgYhBPE6wMDqRzEYfkPUWq90z0pYa8rQdHQjVLqo1oVeLN2zwxESHz6ldA2ANV7HFrq9NgH4aSj6Uso/JsXioXVsVysGIxi4dM9u0CZ3A73BWd2z2pnuCbQu9IrvwpR/iOifmYCfjvR/03oltKugLYPB8I89L0zG2Dc2t0lfhKigSUiHJiEdUb1GcOmeNaVoObcfDqvZrbpnOkv3DIK1v7o5bPdmKzAYjA5M18QovPGTwT7HR/ZMDmu/fLpnVM8hiB8xE3FDpkAVFQfT5aNo3LcczWd+hLW6FNRuDasc1xuvzxkU1vtf1wq/fyarJ8K4/rlvdA+seGqc6/0f7hiA+RN6t6kM6ugERGUPQPzQaYgfPgPqhDSYrxa3buZ+7QIc1vDuKXw98MCYnmG9/3Xt0ll4+wDc9/6eSIvBYATNSzP740+rT4meH9o9CWf+NB1VBjOyk2NgMEXOsvZM9zTBVnvFLd0zjXP9sHTPiHBdW/iDsxMjLQKDAQC4c1i3oNoLWX7H/nCrx3u9Ro3sZE6Jxke1j1WzKm0UdBm9ETvgFiSMngNd11zYGqtgOLQahiMbYCo7BXuLIdJiRoQHxvRo8z6va4Ufp7+uJzAdkofH5USk3++fGhe2zAc59r84Ba/dGZxvNlrnGwRtL0pdKUSthS6tB2L7jkPC6J8gqscgOEyNaDq6AYZDa2C6fBR2Y12n2YD95dtDX/5YDqYRGW3KL2+5AR/vKmnzfod1T0KMTo2vD5S1ed9d4vVwOEKrxP5wR9sri1Dime45qjXd8+R2gJDWhV7xaddlumd2cjSitJ6D+Ifz8sPe73Vt4Yvxx9mDcOHPMyItRqckPSF8GzTL0ScjPKtEfzetr+w1KlVoldZD43qF9H6RhE/3jO49EvH5sxDTbzyISoOWc/tg2L8czef3w1pXAeqwR1rUkPHYeN/vb3L/jLD32ykV/s/H9gz5HyDj+iYhSnwy/Iub2yYj5uzrt2FEjyScem16m/QXCQgh0MSluNI9YwdPgUof05ruWfwjrDWloPaOs2fvtIG+ivy2wZkAgKkDwq/k3WEuHQZDAT8ZkY0lP5YIntNp2sZu0qpV+O7JcfIXXkeooxOgzh6IqOyBXHXPmlKYrxSjuXg3NEldnQXeskA0ukiLKsoTt9yA9ScqPY5lOGe67z+Yj3XHK3BzXpc2kaVTWviMyPDQTTlhvf8jYXRz3DUi+IDvht9MAAD8aVx0m/T9qIDboC24N7+7/EUBoNLHQN+tL+IGT8bWj/8CbUo2rNWX0bj/ezQd3wJzxVk4LO1ve9PhPaQXwU0flInYNkowue4V/he/GOvxvmTRzDbtf8v/u8XvNjuemxgGSTi+/eVNYbu3HDOHZIb1/r+anCd4/M27fFei+kteRpzg8dN/5NwrY3qlCJ5399X2yYhHyaKZyI7378/uf+8e4tf1PA/eKL2I59tf3hjQfQHplbzP39Yv4PsK8bd7h/ocG5yTLpDueQ2Gg6vQdNSZ7mlqCqkc1wPXvcK/8YZU1+snC25o8/57d4lDt0TlgcrHxvdC95TwLUgJ95J7KUblCCvFUJEYI5ymOKlf8H5S74wK7+P/9TIseJ69tU/QfQcab+oh8zsa2TPw7+MViSyhlNjQuleEvj+1irhmdN7pnvrsgVy655H1itI9F/98pOBxqbhNR+W6V/gAZ2X/477heG66p+XxkxFZYenvnw+MAAAcf3UaAOAPswYqbvv7Gf3DIpO/XHwj8CymGIGc8UjSJT40uzelSigytYhSjhYZKISUSahnn1LpjEe9Fm35y5DsJMnzL80M3e84MVp4IL+lr6/fm6jU0KZkISZ3DBJGz0F073xQmxXGk9tgOPADWi4egq2x2kP5j8tN87nP89P74elJuSH7DO2FTqHwe3eJwx1DfVc6LhRZ+DAqJzgreMbgTJQsmula+HXrwK6K27aH7KFTr00PKvf5niB9uGJ/4G3BRw+J50LvWjDJ4/0f75QfyMWe44GXp3q8z0jgBqVNz06QvWewFPTtgoQwL9qSCkKKud78ZXCW9Ep6QlTQJPLpnncipt94QKVCy7k9reme9VcRpfH9ju4b3R19u4auFte23xW4Xu9+YZL4hWGmUyh8McR+9IvuCsxnGiwrnw5vBoaSIN6y+WMFV3X6wwIZH+4vZVxrB16aInpOpw7uJyuVBrf45yMl3T/ebp2fBjGwab0+x54XJgMActPjw75G5F8PCLswQonULE8VooVUUoaB98DJp3tG9xyK+BG3I3bQZC7ds+Qw3nrrLZ90z6QYHfqIxG0CoWdqrOt1ZqJ/QftQ0qkVvpg13TstVvC4Ev40O/Al9HLT5GBRspR7bO/WmMe79w8PqB8xfzfP7RLB25JFM6GRUOpPThQeLAZ246yxdwo8/5i8SwfPuzFHtL3YTKxAwHUwIDNB9nPKceq16dj/4hSULJrpMRMI9ywv2AFdCVJxqPvGBDcDvCefy1oSc6MBwA1dpJW1OiYRUdkDET9sOp544gmo41JhvlKMxn3fwXhqB44ePYpolSMoOb1ZMj22zZNGvOnUCl+MYNwZkSiI5I6SVZ9iTOnvuRdwQZj2Bu7rteJ19wuT8MzkPEVxAzGfOF+cLCnK8yd932jP7yM9Qdif/92T4tlLb98zzPV607O3YNX/jMcaBZtU7HxeOtsqWqeWjS+8c+8wyfPBcnOer/86nPTJiEN6vP+rrT9+aJTr9ZsKZuD+/A0nJibi3V/dg7jBkxGfPwvalG44fvw4/vXu39F0YivMV8+1y3TPQGAKX4QhAVbaFPuhHV44VfB4qHk8iDro//m5p/86lMXnxvZuzQhxt+A3/79bkJkYjWen9lH0RzpNxAq/e6Qyq1GsvIJeI271ursOctPjMEjEd/zWTz3TB/nKlYGw47mJmD2sW9BVNv88RzolVex5BgrvmhLjjZ8E5i6d2C8dw7on4fnp/WR/J4G4/W7uww18Km0UXnxwBu6//348++yz0KX3hq3+qivd01zesdM9O73C3/d7zx/o2/dwf7S3hnjJc1JMa4ZHyaKZOPnaNBT+tsB17Nmpnul7634d2DZnUu4QOaSmyMHyyh2eAc7Rzrx1uam3N6lxwpky7qmAF9+YgZVPjwt4+uxdxErpc5klkBgQKN1TYvDO3OFBFw67X2bGOXOw/2sj+Jngvhd9lXtXmRTkYNKCv39qnGz8BwAOvyJsXPXu4uuq5QeH9PgorP3VzTj9x+l4aiKXnaPX66Hr0hOx/ca70j3tzY0wHuXTPY/BbqzvUNU9r79EUz/xLuY1exiXqhlMjrIYhb8tgNHCBYVidBrkpLU+/v/xSgHrkx7aQl8/PD3e9frG3qnYfaEmpPfn4f3+8yf0xuLtF1zHvXcf++rxwBb9KCkJTAjxKx7yo1f2TaCuLPcSC0O7K+8/kiSLpJrOGNwVJy5dw6VGXz/2CKfSTosNTbqrFN1T/A9wxuiE1VpanB4Xqoxex1o/v9QOeXy6pzYlCzMHzcKKXcdgrSmF8WQhoFJDm5KNrF65qKFxHoN0v67hKdgXKJ3ewveGD5gNk/iDPfjyVFcanTvu28wJkZMWi4HdhN0B3pacWOBOLJXUHX4Jvzvum8HkCASljyxUnpc9USCIydPd6ca4d1R4lteHAm8F3y3JU6kEM9M5/uo0rHnmZnwvERNoa0oWzcRTzmC30sJruenx+Fl/4cHgsfGc29D7N/rXn/quiHVHKk4ixoLpocvnF3JfvTjT/zLTr8wa5Jnu2ecmgBCYzu2FYf9ytDjTPanD3u725GAKH8Ir7aRm0imxOsEFUoFYdf5G7qWm6PxUW64MsNCCM7FVqkJIKfM454IiuVWecnw533fl6tAQ7WDmruD50gje/Mf5m/j5WP/2GI3TazCgW0K7q+H+u2n9ULJopuIMnVlDM5EWLfwZ3Gcy7jGGu2Q2mBkhU1NGCLGSFoFw/+geuMlt5T0A3CpQyVIOd68AIQSa+FRE5wzDHXMfRuygySC6aJhKDqNx/3JcO74DZ86cgdXaPjZzD2r4IYSkAPgSQA6AEgD3UErrBK5bB2AsgJ2U0tuD6TMc3DqwKxb/fCTGuv0Y9CIVEJ9xul7uGNINv1p2uC3E80AqFfAvd0tbWDxdg6xJLxWM5H3yWrUK6359M6a/syOggPVogdo0A2UW2viD3CA7bWDXiKfQRZLUWL3ggiRv/m/ucLxz7zAI7e/SNSEKVxtNivskBPB2h8sZL8O6J+Fwab3r/WsSi+GitCo8NTEXP55vdWd6r4cIhon90vFlUSKX8tl9EBxmI34yKha7d+/Gd999B5PJhNTUVOTl5SEqKjL7QgT7aRcA2EwpzQOw2fleiP8F8PMg+wortw7s6rEQS8xC+5mzIJX3dLZnauQ3ZFZawyTTj9o+QvRV6Jfs1zUBJYtmegSslSL0/Ce0UQlZBufbj1ag8AHuuxJyg733gH/rOOYHkGF2v1fa7R1DWmcc3u4UQoikqzZYvI1ElT4WQ0eMwkMPPYRnnnkG2dnZOHbsGP72t7/hs88+w4EDB2A0GkXuFh6CVfh3AljqfL0UwGyhiyilmwFcFzsVp7oFqtxT0Lb9LnwVLgNhnlulxNXPjPc4553J88GD/m2tFkqryB+S/XA7MSKPuwtHSewgkPz8WV5pq+5BaCHlHs4yxKMEZqX8sdjYWOTl5bnSPYcNG4aLFy/iH//4Bz7++GPs3r0b9fX1YZONJ9hPn0EprXC+vgogqFxGQsh8APMBICMjA4WFhQHfq6mpKaj2YuzYvs3j/b+mxKDZSgPuS0rOngkqnyyJwsJCTO2pwcZLvjv+uN9nYiKQPykGBgtFVfEhFBaLy9B4+QQKr50SPKdRATavRA2pzxqOZ87TcPEoCi8LW518v+H63kNJR5CxsLAQTU1NAHyftz+y/2KwDt3iVNj74w7Za+MaLgbVl/f1U7rYsfOc+L3GZqpl7z+6qxr7rnpurSjWZu8u38/4484d0Km5Z+j9vaelpSE5ORlXrlzB5s2bsWTJEsTGxqJnz57o0aMHkpKSQh4LklX4hJBNAIRWZ7zo/oZSSgkhQSWkUkoXA1gMAPn5+bSgoCDgexUWFiKY9gCAdat9DgV9Ty+k5Cy8haLXC2t8+nd0rcTGJUXByeb22X4yXbyY0z9SK/DLzw/69iPwbPyWQQ6vPm6dPFHw+KKfDEaBc2ofku89zLQ3Gd+KL8Nvvz7icaygoMCpnHxdDv7ILnqlwO/nnhkT8fwO39+7LM57DcpKQEFB6/oVSin+sLv1fvy9/i+pHLvOVSuKee03n8a+q+d9ZRKQf/KkicAGkd8s5L93h8OBy5cv49SpUzh16hS0Wi369++P/v37o1u3biFR/rIKn1IqWsmKEFJJCMmklFYQQjIBXAtaIoYL7y+Y30wjkGyHQBkhslAmNVaHGqOlzeRwZ2C3BJy40uh6H+4aRNc7dwzN9FH4kSBQhfbC6CjURnfDc9M8i/aJ3e/OYVm4c5iy0ugT8rrgva2tCn95GNNtVSoVcnJykJOTg+nTp6OiogKnTp3C8uXLYbFYMHHiRAwfHlh9K55gXTorAcwDsMj5/4og78fwYu6o7li2vxQA8KKzxrhQEPSZMNXuzhDJ6Jk7urvHH0Jb0q+rp8Lvn9n2i1vS0tKQk5PT5v0CgNFoRGxs4AX+hPCuqJOf/4rgcfdzgWI0GpEmIH9+/iuicijh25d9j7nfL1C53e9x/wru2Ys9l2Dkl+P1119HTU1wCyaDVfiLAHxFCHkUwCUA9wAAISQfwBOU0sec73cA6AcgjhBSBuBRSun6IPvuFCy6awjqmi0Yn9dF0gLyt1TvX+4egue+OaroWuOZXa6ysTyhKnErxfPT++HNdad9jv96Sh6+PVjmeh+JnPecnBwUFfm61dqC9uYS8peOLH8kZc/P9y+5Qoig0i0opTWU0smU0jxK6RRKaa3zeBGv7J3vb6aUdqGURlNKs5my94///DxfdgGQv9si+rNJiUoXA6JSg1LqSj8Vc/WEErG6KeHcApLBuJ5hK207MSWLZipaXHRLwS1QxybBWn0J63/NlW2YKFBvJicCaxHkqjMyGG2BxeZAWV1zpMWQhSl8CSKx6Xl75P9NG4DovBvRcuEALC3iC0XC6Vop8toFa3SvFKTF6VzVGdetW4e+ffvigQcewKJFi3zab9++HSNGjIBGo8E333zjcW7p0qXIy8tDXl4eli5d6tO2M2OzO2CzB7YRyJHSejQ0+19SwOGgeGv9Gfx4rtrvtsfLGzBu0RYcK2vwu+1HOy8i/08bsfu8/37y331zBOPf3IqaJrPfbZ/54hCe+vwgrAE+Z39gCl8C7+3wlj4yOkKSRJYxvVPxz1ndsfyVB7Fy5UrRcrDf/TL0GQzH/nArdi2YhLQ4z2J1Xz1+I4pe4ko22O12PPXUU1i7di2WLFmCL774AidPnvS4vkePHliyZAnuv/9+j+O1tbV49dVXsXfvXuzbtw+vvvoq6up8qoOEDLPNLn+RAFUGM/YEWOF0z4Ua/GnVSVi8F1QoYN7H+3Df+3v8bldSbcSd7+3C+De3+F0+ePmhcry79Rye+1ZZjMmdZfsvo7y+BX/44URAbaubLNhw8qpf7cw2O1YcvgIArv+VUme0YOWRK1h9rAL7S2r9ahsITOFLMNwr/THYzc1DybQAij4FQ4Ke4Oabb4bRaMTBgwcFrxErtRsM8VFaZCVJl8jdt28fcnNz0bt3b2i1WsydOxcrVngmjOXk5GDIkCFQqTx/8uvXr8fUqVORkpKC5ORkTJ06FevWrQv55wCAHWer0PelddhWXOV326f/exBzF+/Bmav+LVh3OCjmLt6DD3ZexL6L/imU4koDdp2rwf6SOpyqaJRv4AZffttgtqGiQXk9HQDYe5Fre7XBBIPJvxkCb52XVPtXssBkteO8s3TygUv+DfgX3fo6Vu7fzGKv23dyuiL8xQjaV+3Odk6wG2iHksdvuQHrT1QCAD57dEyb9KlWqzFnzhwsWbIEvXsHvrNWqCkvL0f37q1B6OzsbOzduzfgtuXl5YLXLl68GIsXLwYAlJWV+b0K9M97uW3y/rPuIOgV5bXkbQ6KvRc5//Dfvt+N27Mtivu+0tRq1X9deBC2cuWD8rayVmX75aa9KOiuvLTFqqOtro0vN+zC8PRWVSO30nj/We452RwU36zfjl6Jyip82h0UJdXN0KiAGqMFKzdsRYJOmZvxYoMddgdFl2iCk1casHXrVkEXpZDs+yq4DLY4LXDkQgUKC+sV9QkAG85za1mi1MDWQ8XobbukuG0gMIXvB8HsJhVqRvRIxpt3Dcbl2maMb8N9SdPT0zFu3DisWLEClKa1uzLA4WT+/PmYP38+AC5Fzp/0PEop/qdwAwAHSk1av9qerTQAG7YDAEz6JMTFNStu/+2BMgBHoCKAQZuMgoJRsm14dq0+CZ2aU0C61GwUFCivTf/uqR8xJNuBo2UN0Kb1REFBnuucVGqjw0FxdfN61yY9XXr1R8EQZTuJldY2w75hK2YOzsTqoxVIzx2Csb1T5RsCaD5WAew+iEkDs/FlUSkG5d8kuN+wkOxHNp0FIcW4dXAWtpy+5td3u7rqCNLjq9ArLRYmSlFQEN59FNqPBmP4zb2jeuB3XqsL24Ibb7wRDocDlitn2rxvIbKyslBaWup6X1ZWhqwsZSspg2nrD2V1LTCYbOgSr0d5XYtf/vSSGs6675kag7OV/u2neqG6CRoVQUHfdJTX+bcR97lrTejdJRY5aTE+O0XJUVbXgj4Z8UiN1aG8XrlLp8ZoQYvV7jJiLtcqz3y55HxO43O5tlfqlX9e/tp8p9vWn7aldc3IiI9Cv67xqG+2+hWoLq1rRveUGHRN9K+UdKAwhS+D3CbQnRGVSoXZs2fDVHYc9mb/syFCzahRo3D27FlcvHgRVqsVy5Ytw6xZsxS1nTZtGjZs2IC6ujrU1dVhw4YNmDZtWshlPHeNU9ST+6XDQf1TKJdqOGU7pX8GyutbYLYpD4KW17Wga2IUeqTEoLy+xa8A6qXaZk7hp8a6ZFCC2WZHpcGE7slORdag/LNWOpXeDV1ikRqrQ6kfCr/UmRaZ71wjcs2gPGPmSr0JMTo1BnTjtjks9+P7qWw0ISMxCl0TuVhTpUG54i6tbUH35Gh0TYhCZaM57PvjMoUvw/1jemD3C5Nw8Y0ZkRalXZGSkoKYnkPQXLwblIY/nUwKjUaDd999F9OmTcNDDz2Ee+65BwMHDsTChQuxcuVKAMD+/fuRnZ2Nr7/+Go8//jgGDuQ2ykhJScHLL7+MUaNGYdSoUVi4cCFSUkK/nzEfuORdDKV+5GyX1BiRGK3FQKcyqjP7ofDrW5CVFI2spGg0mW1obPGtsipGVaMZ6fFR6JYU7Zf1eaXeBEqB7ORoZCZG+xW0veq8tmtitHOwUN622qnge6bGIk6vcQ0eSqhoaEFmYhSyk2Kcn0G5wr/WaEZ6vB7pThdQlcKBxuGgqGw0oVtSNDISomCxOVAXQBqrPzAfvgIyE/3fSLkzcPJfT+GGn74Ac9lJRHUfFFFZZsyYgRkzZnj4WF977TXX+VGjRqGsrEyw7SOPPIJHHnkkrPJVNLRARYCRTuuztFa5QrnaYEZmYpRrp7I6k38W/tgbUl3bOpbVNyMxRn7nsBaLHQazzeXHNphsMFntkjuu8VxzKtqMhChkJkb5lW7IDyxdE6KQFqf3q0BfdZMZCVEa6DQqpCfoca1RuYVf2WhC18QoJERroFUTVDcp77fSYMKoXsmuZ6VU4TearLA5KFLj9K71JFcbTIo3MgoEZuEzAiYxWoeYvDEwXzmNb+eFbrPp65GKBhPS46Ncf9jX/Jj2VzWZ0SVejwxnW6UWvsNBcbXRhG6J0chI8E8Z8fKlx+tdikypAuWVdKpzYVxDixUmq7L1B5WNJqgIkBanQ1qc3mW1K6G6yYI0p6wZ8VF+Wfi1RgtSYvUghCAlVodao7J+zTY76putyIiP8lvh84NKWpzO9f3487sIBKbwGUHx+9kj8cT9c7CvcD1sNuXugs5GRUMLMpOioFWrkByjRbUfKzKrGrnBgq9cWm9S5kKrb7HCQTnFy1uNtQotZt7/nZ4Q5XJVKFVGLoUfq0Oqn/3WGC1IjtFBo1YhLU6HaqNFsV+7qsnsWqCXFq9HlR/PuNZoccmaEqtXLC+v3LvE6xGv10CvUSnul1+VmxqrR7KzAm59mF06TOEzguKJW27Anx+7A0lJSdi2bZt8g05KtcHSqozi9Kg2KFMolFKXhR+n1yBOr1Fs4fNWakqszrU1p7+KzN3CV2q51jot1+RYnWsxntJ+G1qsSIzm8v3T4vSw2BwwmJUZEjVNZnRxPuPkGC3qFPZptTvQaLK5lK4/ez3wCjo5VgdCiF+zEveZEN+30ucUKEzhM4KGEII77rgDhw4dEvWTd3bqWyyuPXnT4pRbnw0tVljt1KV0U+N0MFiUKnxOGaXG6pEQrYFGRRQrMncrne9b6aykxmhGYrQWWrXK75lFY4sVCU6FnxrHta1R6E+vMVpcbZJidGg02RTVAeKVdkqs1vm/TrG8LoXvVNjJsVrUtyiz0l0WfpwOCdFaqAhQ38wUPqMDEBcXhxkzZmD58uWwWsM7Le2INLRYXRvXpMXrFStP3qpOcyqy5BgdlMYTeQs/OVbb6ptW2LjBqXgSY7RIiub6VppBUuPhHuHb+m/htw4W8s/K4aDcM3a25QfXBgXKl1fuyW4yK31O/OdKcvaXFK1TrLT5QTUlRge1iiAxWotapvAZHYUBAwagW7du2LRpU6RFaVeYrHaYrA43V4VO8bSftxZ55Zcco0WTVZmF32ql61338MdVEaNTQ69RQ6dRIVanVuxfrm2ytCrPGP+sdHeFz1v6jSZ5l47BbAOlQGJM68AIKBukeIWf4qbwDWabouqV/PfDK/zEGOUWfn2zFfF6jWsFf3KsLuxpmUzhM0LKjBkzcOrUKVy8eDHSorQbeCsz0WV96mC02JUpFKcC4K3s5Bg/XDouX3qri0Rp9kldc6u1DHAukvoWhW4ZU2vbRKerIhALPyHKqfAVKFB+dSvfllfASqxt7+8nPorLVm9SMNDwMyH++0mK1ipeadtoanVfAdx3qzTuEChM4TNCSnR0NO644w6sWLECZrP/tcGvR1xKO4ZXZJxCMShQKPXeLoMYHYwKLfyGllYrnetXq6hPrq3FY+/kpBitYgu/wc0Pr3K6KpQofIeDotHDwueekxILX2hQBZRZ+E3OoDA/wMTzA42CSp11zVbE6rhZEOB8Ti1WRZlF7vEKXmZm4TM6HHl5eejduzfWr2c7WQJuSttpBSb64V92KbKYVt+0yQ5FtXiazDaXtQpwCk2JEgOcFn6MpzJS6pt2V9oA55pRMtA0WWxwUARk4fOzj6QYb4UvL3OT85nE6blnFe/XgGz1HBijdbA7qGsQkYKbzbR+P3+aPQhLH1Fe3C4QmMJnhIVp06bhwoULKC4ujrQoEaehxdvC90ORNVuhVhHEO5VRsh9BUIPJ5lJiAKfIlFr49c0Wl9IEnL5pBdanw0FhMNtcsxh/+vV2y0RpOctZySAVjFuGV86xzmeV4IeF32iyegyqra4kZTK7D4xdE6OQHh8l2y4YmMJnhAW9Xo/Zs2fjhx9+QHNz+9/rM5zUeykjfhqv1MJPjNa6ylDz91AyWBjMNsRFeVrazQpjBw0tNi93gzK3DB88dW8br9cq88M7r3FvmxitVVT/x1vhx/lhpRvMNug1Kpdbxh8Lv8nkNYsK4LttS5jCZ4SNnJwcDBw4EGvXro20KBHFZbnymRzRyi3Iei+lwCsXJX7tJpPVNTNwb6vM6rV6WOkJUVo0mW2yvulGAaWdEK3Mwm/0Utpcv5qALHytWoVorVrRjlk+Sts5SCpS+GbfWRR/XInMHUrhE0JSCCEbCSFnnf/77AFICBlGCNlNCDlBCDlKCLk3mD4ZHYvJkyfjypUrOHHC/z1GrxfqWywebhleoSi1Aj3dI7wyUhaMFFJkcgrUanfAZHV4KLK4KA2sdgqzTOzAW/HyMiuRV6htQrSy2UGTyQaNikCvaVVp8VEaRYpXTGkr6dfoNYuK12td8khhsXHPOCGqAyl8AAsAbKaU5gHY7HzvTTOABymlAwFMB/AOISQpyH4ZHQStVos5c+ZgzZo1aGryb/MOf1i3bh369u2LBx54AIsWLfI5bzabce+99yI3NxdjxoxBSUkJAKCkpATR0dEYNmwYhg0bhieeeCLkstU3i7llFFrpUZ4WL6Dc3SCkyOTaGp1KMs59oNErs1xdFn6U56xEyYyEH4gSY9w/r0KFb7YhLkrjsQNbnMLYQZPJ5vlZ/XQHxelbK4i6XElmaZn55+g+ILcFwSr8OwEsdb5eCmC29wWU0mJK6Vnn6ysArgHoEmS/jA5EdnY2RowYgR9++CEsGzzY7XY89dRTWLt2LZYsWYIvvvgCJ0+e9Ljmww8/RHJyMs6dO4ff/OY3eP75513nbrjhBhw+fBiHDx/Gv//975DLV9/imdMepVVBqyaKLHxf61O5u8HgpcgSFPr/+Xt7W/hK+jUIKLJ4pzvI7pD+7kUt/ABcK3y/StxBBq+2GrUKMTrl7iCP56RX5jbjz8fqO5bCz6CUVjhfXwWQIXUxIWQ0AB2A80H2y+hg3HLLLaivr8eRI0dCfu99+/YhNzcXvXv3hlarxdy5c7FixQqPa1asWIF58+YBAO6++25s3rw57LsL8TS2WD2sVkK43HQlykjc+pRu63BQNFlsgj58OQXKW5+eikyZq6JJYLBIUOjXbmjhMpJida0Wc0KURrFLx1vhJ/hj4es9XStKMovsDooWq92jrev7kfmskbLwZXsjhGwC0FXg1IvubyillBAi+hdECMkE8CmAeVRkiyRCyHwA8wEgIyNDcmd7OYR2l2+PdCY5MzIy8K9//Qt33HEH4uLiQiMYgG3btkGj0aCwsBBNTU0wGAw4deqUh7zFxcUoKytzHdPr9Vi5ciVaWlpw/vx55OXlISYmBo8++iiGDBki2M8PP/yAVatWAQAaGhoUP4/SyhYk6IjH9RqHFeculaOwsEaybZ3RhIbqqygsrAPAVc9UgeLYmfMopKWi7VpsFJQCleWXUVjI2WRVzdyf3f7DxxBVfVq0bXEdV7v+/JmTKKzh9i0+V8Md27m3CDXnxDdBOXSJU85HD+xDiZ5zr1wp445t2LoDXWJUor+lU+fMiFFTj6qr9dcsqG+2YuvWrR7uGm9Kr7bA7oDHfVsaTbjW5JD9nqrqm5FIjB7Xqe0WnC+9gsJCz81b3GXnF8BdLStBYWE5AO77URPgZPEFFEK8kOCZWu55njt9EoXVyvaGNhgMwesKSmnA/wCcAZDpfJ0J4IzIdQkADgK4W+m9R44cSYNh69atQbVvKzqbnNu3b6dLly6lDocjJPejlNKvv/6aPvroo5RSTs5PPvmEPvXUUx7XDBw4kJaWlrre9+7dm1ZVVVGTyUSrq6sppZQWFRXR7Oxs2tDQINunP7/P8W9ups98cdDj2J3v7qQ/+2CPZDuHw0FzFqyib60/7XF8wEur6MvfH5NsW1HfQns+v4r+d+8l17F6o4X2fH4V/WDHBcm2W05V0p7Pr6IHL9W6jh0rq6c9n19F1x+vkGz77paztOfzq2iLxeY6tuboFdrz+VX0RDn3XMV+S//z34P0lr9s8Tj2z63naM/nV9Fms02wDc+sf+yg8z7a63Hsua+P0NGvb5RsRymlw1/bQF9a7vk857y3kz7wvu/34y57WV0z7fn8Krps3yWPa4a9ul72+xF6xnIo/c0BKKIiejVYl85KAPOcr+cBWOF9ASFEB2A5gE8opd8E2R+jgzNu3DiYzWYUFRWF7J5ZWVkoLW21dsvKypCVlSV6jc1mQ0NDA1JTU6HX65Gayu0zO3LkSNxwww0hXyzWZLL5ZGMo8U03W+ygFD6uimgNkfele60eBVr98LI+fEE/vDK3TJNZKFtGWWYRH3h1p7W8grzM3s8p0KAtL7OsvC73led3q6Rfg4DbrC0IVuEvAjCVEHIWwBTnexBC8gkhHzivuQfABAAPEUIOO/8NC7JfRgdFpVJhzpw52Lp1K2prle91KsWoUaNw9uxZXLx4EVarFcuWLcOsWbM8rpk1axaWLuXyC7755htMmjQJhBBUVVXBbuem1xcuXMDZs2fRu3fvkMgFcDNoIUWWqCDdUChbBuAVvjKl7d5WrSKI08srI1e/boosTmGWjlEgW0ZpTZwmsw2xOt/nBMgPUkI+/PgoDZotdsma+GabHRa7Q7Ct3HNqEvl+4vTypSTEvttwE1RvlNIaAJMFjhcBeMz5+jMAnwXTD+P6Ii0tDRMmTMDy5cvx8MMPQ6UKzu7QaDR49913MW3aNBiNRjz55JMYOHAgFi5ciPz8fMyaNQuPPvoofv7znyM3NxcpKSlYtmwZAGD79u1YuHAhtFotVCoV/v3vfyMlJSUUHxMAYLY5YLVTwYCiUkvbu22MRoHydJ6PF+pXqeUa5Ts7kFWCJl+lrdjCN9nQLcmztIDStQNGkSwd7pwdiTHCvzHXcxKw8JUHtz1jGvF6DZrk0jIjlKXTtr0xGE7GjBmD06dPY/fu3Rg3blzQ95sxYwZmzJiBwsJCFBQUAABee+011/moqCh8/fXXPu3uuusu3HXXXUH3L4ZQxgvAWa4NzqqKYsFIMWWkxKXTmgXi60pSMjsgBIjRtioyvUYNnVqlyFXhqzwV5v9bbD4KkG8rlcJqd1AYLXYBt0yrO8g9S8odse8nIdqPgVHApSO3/6+rfo+ubVUwK63AiAiEENx5553YtWsXrl27FmlxwoZRVKFoYXNQNFvsom2bBFwrABCtVWYtA74ug/gojeyCryaTDXE6DVQqz4EoLkrechW2tJWlkkpZ6U1m8edktIg8YwUDjdCaA66tFhabA2abRL+iLh2NfPqq2YYYnRpqlXjmUThgCp8RMZKTkzFlyhQsX77c5Ue/3jBIKF5A2icupoxiFFj4jQJBW8BZE192FahV0LesVJF5W+l6jRp6jUrWReJd3RNQNliIDar8QCnVVswPr2RWIuZyU+L/Fxrc2gKm8BkRZfjw4YiLi8OOHTsiLUpYEHMZKPFriy3OidYQ2UJm4v0qsPBFlJGS2jRCAWqurbQryWp3wGxz+AwWSlauSs1meJnk2sZ7u2Wc/RoVtHVfKMbLIbfwSiirqC1gCp8RUQghmDVrFvbv348rV65EWpyQI7TyFGgNpkpZgt4bc/DEaDi/taQ7yGRDrIDLIC5KI6nEeJnELHwlQds4Ab90gsxAI2alx+jUUBGZmZBXPXseJVa6eKaNfFujxYZordq1J62rX71GkTuorTN0AKbwGe2A+Ph4TJ8+HcuXL4fNpmyDjo6Cy78cgMvAe2MOnhgtUdRWWGnLpwyGx8KXbis2IyFEPpWUHyy8M5KUzKLE3DJxCl1uYgMjJ5fcgMwUPqOTMmjQIHTp0gVbtmyJtCghxWASVkaKFIrXxhw80Rpe4UsoMgF/OMApXotd2voUymkHnD58CXn5WYdgW6UKX9QdFLhLRyp2IJqWqaB2kNjAyJdLlm3LLHxGZ4UQgpkzZ+Lo0aO4fPlypMUJGWJWuiIfvsk3xREA+G1QpRSZ925XPIp84qKKTMbSFsmW4Y9J9WkUeU4APzuQt9K9LWa9RgWNikj74c1Wn5XBnBxq53lpl5vYZwWk1w6IPeNwwxQ+o90QGxuL22+/HcuXL4fFomzD7PZOk8mZ0+4d2FPiwxdRCryFL6eMvGcVgPJAppg7KJDgqattABlJ/DElFr734EgIkY1Z8J/Vey1Ea117qUFKeDajJB2UZekwGAD69euHnj17YuPGjZEWJSTwSttHoSgK2gor3hhe4csMFkKzA7l+KfUtq8wj5w4SC7zybaVTK+2ibeXcQVKzg1iddMaMWLYM79KRGiwMAimo7nJIzyyE24YbpvAZ7Y7p06ejuLgY5893/G0TmszCylNJXRsxZcTrcTl3kJjy5OUSwlWwTSTwyt9bTF5AwqUjkUrKu2zEfPhyg5teo4JW7avO4pVY+ELPWKuCWkVk+rUKD6rOY7yLyxuzzQ6rnbZ5LXyAKXxGOyQqKgqzZs3CypUrYTJJL1Fv74hZ6QCvBOWUtq8fns/SCSSDRC4YKba6l5dXql+pgmBxURo4KNBiFZ4d8CtphVI64/TSWyQKlXPgiZUJNIvNhPjsIOmZhUiAWmYWJZa/3xYwhc9ol9xwww3o06cP1q5dG2lRgkJq6i63IlNMGfG1usSUoNBuVzxyFr7YymBAuSKTUoKiA42rmJivEpQL2kr5w+WCxVLBUyWxA+m0TLGB0Tm4tfEG5gBT+Ix2zNSpU1FaWorTp8V3Z2rvSCoUBamKQm1VvPUpooyarZxbxrtwGuCmtEX6da3ulRgsxJSglEtHbus/o8WGKK3KZxETfz+T1QGrSJnjJpP4oCr7jE3C2Ux8v2IDjVhZZYAL0BOJxWJ8aQvvKpttAVP4jHaLTqfD7NmzsWrVKhiNRslr161bh759++KBBx7AokWLfM6bzWbce++9yM3NxZgxY1BSUuI698YbbyA3Nxd9+/bF+vXrQ/oZxKx0QL4Er5Q7SMrqlcqWkatNI9mWdwfJuXQCsPDF1g24yywVOxAdVHXSCl+yrcRgIRVkJoRI9itWZbMtYAqf0a7p0aMHhgwZgtWrV4sG/Ox2O5566imsXbsWS5YswRdffIGTJ096XPPhhx8iOTkZ586dw29+8xs8//zzAICTJ09i2bJlOHHiBNatW4cnn3wypIXcxIKCgLNuuojilbIgAWl3g9BuVzx8brq4D1+8bas7SHqwELK2lfj/pT6rXFuxQZVLy5ReZCbaVq8RrdIp5b4CnLEDkWfMB3OF3Ffhhil8Rrtn0qRJqKqqwvHjxwXP79u3D7m5uejduze0Wi3mzp2LFSs8d9tcsWIF5s3jduO8++67sXnzZlBKsWLFCsydOxd6vR69evVCbm4u9u3bFzLZA/Xhi+WWu7cVdxmIW+l8brqsD18qh18i4KsTWBnsLovo55V5ToG25YO2DoevsWCzO9BiFQ688jKLDci8W0bKlSSWpWOQ+W7DSdv3yGD4iUajwZw5c/D555+jZ8+eSEhI8DhfXl6O7t27u95nZ2dj7969otdoNBokJiaipqYG5eXlGDt2rEfb8vJyQTkWL16MxYsXA+D2zS0sLJSU2+Hc3rDmajkKC6t8ztdXmdHQbBO8z7Vmzl9deuEsCs0lHueamppgadagvoEKtj1ezSmUcyePgVT4WpFaasO5S+UoLKz2OXe4hFNkR4r24rzOc+2Axc4pzSOnitHdSyYAKL5ohl7lkPw8RUeOYVii2eea8soWABBse6GGs7J37NmPaym+n6fW0IzGGotg28oy7vOs31LoWrDGY7Ryn+dqaQkKC32/c0OtGTWNdo/7NjU1obCwEGdqndtinjmBwmrfGJPD3ILLV1oEZTpUysl09OB+lEW1rc3NFD6jQ9CtWzeMGjUKK1euxAMPPCC6S1Q4mT9/PubPnw8AyM/Pd+2sJUaT2QasX49BfW9AwYQbfM4fsZ3FupJijL95gk+w8nh5A7B9J0YNG4yCgV09zhUWFqJnZgJOX20UlKH5WAVQdBATbhyNvl3jfc6nHdmB2KRoFBTk+5w7tvkscLoYt066xcdSp5RCu2Ut0rv1QEFBP5+2y68eQrKxXlCmWqMF2L4R2Tm5iLNe8rnmL0d2IDMxCgUFo3zaJpfWA/t3Ibf/IBT0z/A5b9m0Fn1690BBQX+fcxUxl/HlmWMYPupGdE303D6xrK4Z2LwVwwb1Q0F+d5+2O5tO4sC1yx6y8juqOU5XAvuKMG70SAztnuTT9oNze9FssaGgwHc3t+Lt54ETpzG14GbBwHo4YS4dRofh5ptvhtFoxMGDBz2OZ2VlobS01PW+rKwMWVlZotfYbDY0NDQgNTVVUdtAkQvOSZU5kComxreVTY8UayvhX+YXMQm5ZeQqV0r54Xl/tWhbi8R6BYnnxJUhdgjm73P9iscdpDKS+LZGix12AXeQQSJewR1XiwdtnXEBVi2TwZBArVZjzpw52Lx5M+rq6lzHR40ahbNnz+LixYuwWq1YtmwZZs2a5dF21qxZWLp0KQDgm2++waRJk1y1+JctWwaz2YyLFy/i7NmzGD16dEjkbXL5eYWDc1J+bbGNOVxtJRYFie125d6vuDISD2LybcXyy8UWewHOPXE1KsnAq5wPXyijSWqxF9CqzIWCr7IDo8SKWbHNaXik6g7xexV4byHZFjCXDqNDkZ6ejnHjxrmCsIQQaDQavPvuu5g2bRqMRiOefPJJDBw4EAsXLkR+fj5mzZqFRx99FD//+c+Rm5uLlJQULFu2DAAwcOBA3HPPPRgwYAA0Gg3ee+89qNWhyZ6QC85JFdmSt/C1aHZan96bnIjVlueJ02twsVo4zVWuimOcXiuZw5+RECV4DuCUr8FsA2J8z0mmZUqsDpb7rC4LX6Ct1LoB9+NGsw0JXq4XqRRU7ri4hS81uIUbpvAZHY4bb7wRZ86cwd69e10B1xkzZmDGjBkuHysAvPbaa642UVFR+PrrrwXv9+KLL+LFF18MuZytudpiVrp4XnvrLk7Ss4Mmkw2JMZ73F9vtyr2tlDtIqk67lDtIrgIkl/Xi25bf3lCsrauujZRbRiK1krtOoK3MgOz+jJHo21aoCqp7W752kHe8KVK18IEgXTqEkBRCyEZCyFnn/8kC1/QkhBwkhBwmhJwghDwRTJ8MhkqlwuzZs7F9+3ZUV/tmmrQXpHLaAelFUHIundaVq8JKUE5piy28kttrVc4dJGW5irmhpKpdAlzsQCxmIbbfAE9rnETApSNRN4iXFxBeHWwwc1s5iiUPxOm1cFDAZPVdHRypWvhA8D78BQA2U0rzAGx2vvemAsCNlNJhAMYAWEAI6RZkv4xOTkpKCiZOnIjly5fD4RBech9p5Fw6UsHIJrMVahVBlFb4T1RqT1wp9wjAKTKzzQGLTUAZiRRsc28rlcMv6f8XmR3IBU+VtJV36YgPqlLF7dyvc0duT1q+bILogNxBFf6dAJY6Xy8FMNv7AkqphVJqdr7Vh6BPBgMAlxqp1+uxc+fOSIsiiFLrUygYya/QFbUgJQYLsd2uvNsKBV+VBG2FBhk5twzgzCySyEiSmx0ItpUZVHmXmFFgw3eD2emW0cq4zURklpsJAcL72nZkH34GpbTC+foqAN8kWQCEkO4AVgPIBfA7SukVkevmA5gPABkZGbILW6TgF0i0d5icwZGSkoLPP/8c1dXVSElJaVdytroqhBVKgsQ2h3KuFaktEptMVldAWKptk9mG5FidZ1u5fkUKism5Zfh+m8wGeNt8cpk2XFthN5TcYKHXqKFTq0QzoeJ0GtFsGSkL3yBRsA1oTbkUays1mxHDYDAE/duW7ZUQsglAV4FTHlEuSiklhAgWO6GUlgIY4nTlfE8I+YZSWilw3WIAiwEgPz+fyi1skcI9eNeeYXIGT3Z2Nnbv3o1Zs2Zh586d7UZOg4krNaDXCCt8qbo2UjVeAOlSxXLZMpJtZYK27pUr3Tcc8c8t4znItJZzEM+Oio/S4prBd28EuWwZgI87CA0WVplYh3hQXclMSKyt0RKYhR8fHx/0b1u2V0rpFLFzhJBKQkgmpbSCEJIJ4JrMva4QQo4DuBnAN35Ly2AIMHToUJw6dQrbtm0LWUplKGg02SQtbblgpJQSS5ByN8j48MUWfMkVbAM83UFJMa2KWy6NlD/HZa54uptaYx3SsYMLVcLWMiC9iClOL1xATe4ZS21kbjTbkBEvP6h6t6WUyg6q4SRYf/pKAPOcr+cBWOF9ASEkmxAS7XydDGA8gDNB9stguCCE4I477sChQ4dQVeVbsyZSGExW2aXzYlkvcpk2Uou2pBZAAeKpinIVIN3PefcrVSnTva3VTuGduCKXWgmIxw6azPKLmGJFVgfLPSeNWoVorXA+vZIANeAbJzHbHLA5aIcN2i4CMJUQchbAFOd7EELyCSEfOK/pD2AvIeQIgG0A3qKUHguyXwbDg7i4OMyYMQM7duyA1Sq+O1JboiQbI16vFU3LlFKe0Vq14J6rUrtd8YgNFkqqOIrNDuSyZdzbeuteqXLO7m2FgrZy2TIA52YSC1DLfT9igwXncpOekQC+KZ1KnlM4CapXSmkNgMkCx4sAPOZ8vRHAkGD6YTCUMGDAAKSmpmLz5s2YPn16pMWRtQIBccuVcweJK5TWujaeg4XUblc88SLuBmWuFWG/thJFxp9rsXmG+gzORUxSbpl4vcZZN8fuERORS0EFONdMdZPF53iTyYbMRHG3DCBchpofVJXMwLwHGiXB7XDCUiQZ1xVjx47FyZMncfHixUiL4nTpSP9hJ4gofIPJioRomcFCIFVRLrfc/Zyvhc8NHoqCkSIuHSWuJCGFH6cXz5Zxb+vdb6Mit5k2oNRKvl9vpd1ksYFSSMZnorVqqIivvEoC1OGEKXzGdYVer8cdd9yBFStWwGw2yzcII3LTfoCzpr0X55isdphtDkkLn2vruxhJiXtETBk1KnDpyLoqpIKnUbzC9zyuJE3RPZXUW+aEaBmFL1LXRm6RGdfW9xkriXUQQlybr7jDK3y57zZcMIXPuO7Iy8tD7969Q74/rb8ocTcIKZRWpSCnBH1nB1K7XfHw7iBfZWR19ivhDhKx8BudbhlJ/79TuZrs3ha+suA2d62vzLJuM4Fn7HLLyFjacQKxAyWuL1e/Is+4revg8zCFz7gumTZtGi5cuIDi4uKI9M8rFKVK232/XpfilbVcfRVKY4u80ub61QYUtBXL8Glsscq7ZUQsfLmcdkC8lIRBJvWVk1mLFqtnXXuXW0bmGQstNHPVSApgoHEN5jLuunDBFD6jQ1NbW4upU6ciLy8PU6dOhcFgAMC5dmbPno0ffvgBzc3NWLp0KfLy8pCXl+eqiw8ABQUF6Nu3L4YNG4Zhw4bh2jXJpSSKMToVihLL1eagMLvVtVHiWuHO+2b48G0TZRRKvMBiJKn9bHlidGoQQXeQVXaQkfLhKwluA0KBZvnZgSuf3k1mpQOj0GbkSr8fbgMV3+fEtWUWPoPhN4sWLcLkyZNx9uxZTJ48Gf/9739d53JycjBw4EB89dVXePXVV7F3717s27cPr776qscGKp9//jkOHz6Mw4cPIz09PSRyKd2omv/Db3RT3EqVkVAOf0NL4LMDg8mKGJ3aZ7tFd1zZQT4zC3lfenyUmMK3Stb+4dryPvzW52SxOWCyOmT9/7xcns9YmaUt9IybgnG5KfxdhAum8BkdGn4jFACYN28edu3a5XF+8uTJKCwsRH5+PlJSUpCcnIypU6di3bp1YZXLZS0rdFW4W5Gt0355BepdeM2vwUJAGSnJD48XyE1vlKnfA7SWkvBOSlLi0hFa8KXU9ZXoPM8Phry8gILn5FwsZra1rtRtnQnJzA50vhk+jS1WRGvVHmUp2pLIDDMMRoiorKxEZmYmAKBr166ora31OK/VapGdnY2ioiI0NTUhLi4O2dnZKC8vd13z8MMPQ61W46677sJLL70kWqFy8eLFWLx4MQBu71upQlZn6zgFcfHMSRTWiscRSq5xCqHwx724nMi5HvaXcsro+KH9uBLlqxj4AnFVVyyw2BzYuGUrtE7f+fFiC7QqYM+uHaJ9AkBzgwmVBofHZzhfaoLa4ZAt0KWym3GhtAKFha2zpIrqFqRFE9m2ejVFY4vF47r6Zgvqr1WgsLBGtJ3FGeg9erIYheYSAEClkXODlV08i0JLiWjbi7Xcd7F9TxGqU7lnfLCSe+7FJ47AUiYeuK24zH0X67dsR4KOoKmpCYcvngYAHN6/G1Ea8ZhFU70Z1Q12j89aXGKGXiX/jMMFU/iMds+UKVNw9epVn+Ovv/66x3tCiKCyTkpKQmZmJn744QfMnTvX49znn3+OrKwsGAwG3HXXXfj000/x4IMPCsoxf/58zJ8/HwBXmlmqkBU9cw3Yux/jxozEiB4++wK5iLlYi3cO7kafAUMxPi8NAHBm23ngxGlMmzhBcIEOX8jukq4E3509gZFjxiHFWfVyfe1RJFVfky2ytb72KEpOeV734fm9yNDZUFAwTrJttzO7QQAUFNzoOubYswW9slNQUDBMsm3y3i2wEaurX7PNDtu6dRiQ1wsFBXmi7Sil0G5Ziy5ZPVBQ0A8AcKysAdixE6OHD0HBAMFCvQCAjIpGLNq3A736DEDBYM44qD5QBhw6gknjb0SPVIE9F53UHizDZ6eOYOjI0eiZGovCwkJk6DKhKj6HaZMLRI0DANjRdBIHqy57POMvyw4gzdKEgoJbRNuFE6bwGe2eTZs2iZ7LyMhARUUFMjMzUVFRgeRkX+WalZWFM2fOoL6+HkeOHEFZWZnrjzArKwsAV4nw/vvvx759+0QVvj8oTa0UynoxmGxQq4jo9nnebQ0mq0vhN7RYXS4Muba+bhn5jBcASIrW4nJts1db+aAt32+LrXXVa2s9e+m2QiuLlSwUA1pdPh4uHVesIxBXkvReBTxJ0dy+w+6rg5UEqMMJ8+EzOjSzZs1yZd0sXboUN910k88106ZNw6ZNmzBx4kSsWLEC69atw7Rp02Cz2VxbJFqtVqxatQqDBg0KiVytC6CUBTLdffGNztxyOYUilJve2KJMafOpijZ7a3ZQk0KlnRitRX1zq/J0OCiazPJBW4D7vCabewqq8iCmd9yhUaHCl/Lhy66TEHjGBgUL6gAgKca3X4PCZxwumMJndGgWLFiAjRs3Ii8vD5s2bcL9998PACgqKsJjjz0GgNsk5eWXX8btt9+O999/H1OnTkVycjLMZjOmTZuGIUOGYNiwYcjKysIvfvGLkMgltxMTj9BCpsYWZUohXkAZKbXw+bTNRh9FpsDCj9F6KjGzfKkBnji9Bs1uEwslqaA88XrPEgmNrlmUXPCUKzTnaeFzVrpURpL7vRu9ZhZKnpMrO8hjoImshc9cOowODV8sjYcPhuXn5+ODDz5wHX/kkUfwyCOPwOFw4KOPPkJRURFGjRqFAwcOhEUu/90ynopXycIcoQ06Gk1W9O4SK9s20c365N1ByhW+Di1WO0xWO6K0ajf3iPxAkxCtRbPVzcI3K89Lj/PKSlJapoAQgsRorYClLf9ZhWYHSrKKALj2C3CfDXE1kpiFz2C0CSqVCrNnz8bWrVt9MnpCicFkRaxOLeuWaa257uluiJdxBQHuFn5rW+UWvqcis9odaLHaFSleb8tVaYoj4PRrB+jSifdaBOVymylU3N4uHSWKlx8YG1s8YyxKZiSCrqQW5sNnMNqUtLQ0TJgwAcuXL4fD4ZBvEAD1LVaPHaGk8F6gwy1iUu4y4BUKpVSxO4hXRvXNXABVqQsK4JS2e79KFzHx/TZbOb8/11b5YOG9CKqxxeZy18iREK31jJO0SJefdvWp00BFfP3winz4rmfMtTVZuR3FmA+fwWhjxowZA7Vajd27d4fl/nXNViTHKPvD9q6Jr6TkL8ApT0JaFUqT2QYHVa54gVZFprQgmHvb+gAs/MQYHahbf7zsSbHKZjTeWTpK3SPCFr78c1KpCBK8gtRyO2W59wn4PmMlrqRwwRQ+o1NCCMHs2bOxa9eukNXPcae+2eKHha/1KFWg1EpXqwgSorQuK721jo4Spa1z9QUAdc57JClo68o+aeYtfKsf/fKDhcX1v0ZFZMsjAFxmkXuhOX9SHBOjtV7BU+XZMu6DBaUU9S3KBvMEkYExUnV0AKbwGZ2YpKQkTJkyBcuXL4fd7rvJdTDUNVsUW/iJ0Vo0OBWu2WaH0WJHigKLF+CUb51T8fIK2D+XDtem1tl/cqz8IJXkHCxaFZnyGu/e7qC6ZiuSYrSysQ6A+6w2B0Wzxe66h3KlrQnYl+6u8JttgN1BkaxgMFeruE3qG70t/AhVygSYwmd0coYPH464uDjs2CFdisBf6o3KffipsTrUGDmFW+v8PyVWr6htUozOZZ3zFqQSS1un4YLFvNKuc/UrL7O3/59XaIrcHDGeA01Ds7Igs7tstW7PSom8vMwNLVZQSuFwUL/dQfxzarJwswslCh/gBinv58QsfAYjQhBCMGvWLOzfvx9XrlwJyT2tdgcMZptipZASq3MpsZom5YoXAJLdcuKVVsp077fOe6BRIDO3KMzTHZQYrVUUPE3ycnNwMyGFzynGU+HXGC1IjVOu8O0OCqPFDqPFGetQqHjdn1OTM6U0WeEMzH12wA/MSr/bcMAUPqPTEx8fj+nTp2P58uWw2Xy3wvMX3npVqhRSYnVotnB57byln6ZQkSW7WfjVTWZnW2Wzg7S41plFXbOFiwkoDWRGtVq91U1mxfK6rHSnrPVOl44SeHdTbbMFDgdFXbN/Fj7ADYr8gKG039RYvauNwV8LP1rn9pyc363C2Vs4YAqfwQAwaNAgdOnSBVu2bAn6XrziTVX4h53qVFo1RgtqjVxb5Ra+DrVORVJt4P5XavWmxOpQ4+yv1mhFcoxOkS8d4F0VVle/qQoHmeQYHVSkVfnV+RHc5p9JndGCRpMVdgf1y/XFt60ycJ85PSFKUdvUOB2azDaYrHYYrf4p/ES3GVhNkxkahYNquGAKn8EA59q5/fbbcezYMVy+fDmoe1U2mgAAXROVKaNWq9ficukoHSzSE/Scm8JsQ1WTCSmxOsW11lPj9K7+ao1mxUFmgHPNuCxXoxldFCp8lYogXkdQ3WSGw0FRZTAjPd7P52S0uGYmqQoHRr6PawaTS+Erldl9QHaOqcoVfrTWFUyvaeJcUEoH1XAQlMInhKQQQjYSQs46/xetA0sISSCElBFC3g2mTwYjXMTExOD222/H999/D4vFIt9ABF7hZyi2IDnFU200o8pghk6tUmwFtioyrq1S1wrXrw41TRZQSnG10YyuicrkBYAu8Xpcc37OaoNZ8awCABKcCr/GaIHNQRX3mxClgUZFnDMh//zhvDV/rdGMKucMrIvCgYb/fmqbLKg3O6DXKP9+UmI4l47dQVFjNCseyMNFsBb+AgCbKaV5ADY734vxRwDbg+yPwQgrffv2RY8ePbBx48aA71HZ6HQZxCtTZJlOhVdRb0J5fQu6JUUptgL5Pq41cparUiUGcL5kizPAfLWhBV0VDlAA0DUxClcbTbDYHGg02RTHDQAgUUdQ1WRxDYxKnxMhBBkJUahsMOFqg3+DKm/N8wOjiigfLPjrqo1m1Jm4AUrp95ORGAW7g6K6yYzqJuVB5nARrMK/EwC/I/RSALOFLiKEjASQAWBDkP0xGGFn+vTpKC4uxvnz5wNqf7XRhNRYHXQaZX9eGQlR0KgIyuubUV7fgqzkaMV9pSdwiqzSYEZlo3LXinvb0tpmXDOYkZmkvN/MxGjUN1txsdoIAH4NFol6gmqD2c31pbxtZmIUrjS04Ep9CwCgW5KytjqNCskxWlwzmHCt0YzUOL2irCIAyOCfcYMJdWbq12fNdF5b4RyklA5u4SLY6EEGpbTC+foqOKXuASFEBeCvAH4GYIrUzQgh8wHMB7iNLYLZBozfBq69w+QMjsbGRrz22mu4evUqunbtit/+9reCcj733HM4efIkBg8ejDfeeMN1vKKiAq+99hoaGxvRp08f/P73v4dWq0VGRgbefvtt3HnnndDp/LPKyupakKlQEQHcAp3MpCiU1XGKbEJeF8VteZdOWV0zKhpa0CM1W3HbHincTk9FJXWgtHWmoQTest5XwhWgk9o1ypvUaILdFS0oq2tx3kv5INUtKRqHS+txpb4FCVEav3La0+OjcLXBDKPZhmw/BtWuCVFQqwhK65pRZ6IY4Mdz4gezy7XNqDSY0D1Feb/eGAyGoP8GZRU+IWQTgK4Cp150f0MppYQQKnDdkwDWUErL5KZBlNLFABYDQH5+PpXbpk0Kfhu49g6TMziee+453H333ViwYAEWLVqEFStW4IsvvvC57o033kBzczP+85//eHyOe+65B6+88grmzp2LJ554AufOncMvf/lLAEBiYiKMRiNuvfVWv2Q6f60Jo3LEtzUUIjspBuermnDNYEY3PyztxGgtEqI02F5cBQcFeqUpV7w9U7kyyrvPc3vJ+qPw+Wv3XuDa8oOHEjJiCBwU2HG2ClFalV9Wb2ZSFNYdN6GsrsWv5wQAPVO5Z9xisWNs71TF7TRqFbolReFSDafwMwJ4TgcvcYNqdrLy5+RNfHx80H+DsnNOSukUSukggX8rAFQSQjIBwPm/UFGSGwE8TQgpAfAWgAcJIYuCkprBcLJixQrMmzcPADBv3jzs2rVL8LrJkycjPj7e4xilFFu2bMHdd9/tav/999+7zk+dOhWlpaU4ffq0YnmMZhvK61uQlxEvf7EbN6TH4nh5IygF+mcmKG5HCEGfjHjsucBZ2rwSV0JyjBZxeg02nOT2C+7jh8y56XEAgM2nrkGnVin2pQNAegyndjaduoYbusQpdq0AQPfkGFjsDuy5UOO38szLiMP5KiOuNJj8ek4ANyDvOlcNGwVu6BKnuF1KrA4xOjV2nK3i7uPHzCIcBOvSWQlgHoBFzv9XeF9AKX2Af00IeQhAPqVUKrjLYCimsrISmZncxtRdu3b1q8Z9TU0NkpKSoNFwfwbZ2dkoLy93ndfpdJg9eza++uordO/eHZ9//jkWL14MACgrKxOcXh+t4hZuma+VoLCwTLEscW7bQDWXnURhtfgg4+1ei3OYXa8rzhxG4QXlCrRbjAPFZiBWC5w5tAfFCoORlFLE6wCDxY7u8Srs2L5NcZ9xtAUA108CNfrlprA0cHV0jBY74q21frW11rQ+45aqSygsLJe42hO91Yw653fUXF6Mwibl8Z3sWIriKi7WUVF8BIWXI5cNH6zCXwTgK0LIowAuAbgHAAgh+QCeoJQ+FuT9GQxMmTIFV69e9Tn++uuve7wnhIQ8x7lHjx4YOnQoVq9ejV/84heYP38+AG5HLe/pNaUU7/57N7omEDx1V4Fr42olDDCY8NGJrchOicacabdIfg5v95oprQKFnx1EerweM6dO9OvzHbWfxdsbi9G/WzImTvTdD1iKgWd3Y8+FWkwd0hMFBQMUt+OUNKcAZ4zpj4KxPRW3vdFmx2u71wEA5kwYjgl9lMc7cuuasfjoVgDAE3dOULzgCwAMyVdQ+MUhAMDcGQWI0ir/brcZTqB4VwlSY3WYM32S4nbhICiFTymtATBZ4HgRAB9lTyldAmBJMH0yOh+bNm0SPZeRkYGKigpkZmaioqICycnKfeepqamor6+HzWaDRqNBWVkZsrKyfK6bOHEiFi9ejOPHj2Pw4MGi9yOE4Pnb+qGh2eqXsge4gOKhhVOh16j8HrSmDeyKX03Ow7jcNL/aAcBP87NRUm3E05Ny/W770swBeGXlCdw9UnmgmOfPcwZjyY8XcfcI/9rqNWosuK0fTlc0YnSvFL/aZifH4P9N7YNmq90vZQ8A43PTkBanw4hUh1/KHgB+MjwbH+8qCeg5hRq2py2jQzNr1iwsXboUCxYswNKlS3HTTcqtVEIIJk6ciG+++QZz587F0qVLceedd/pcp9FoMHv2bHz++efIycnxiQW4MyrHPyXkTqyCmvBCEELwm6l9AmqbmRiNt+8dFlDbQVmJ+PaX/s0KeO4f0wP3j+kRUNsnbrkhoHYA8D+T8wJqlxyrw97fT/HLdcUzODsRp/843a9YRbhgpRUYHZoFCxZg48aNyMvLw6ZNm3D//fcDAIqKivDYY62TzJtvvhk//elPsXnzZmRnZ2P9+vUAgDfffBNvv/02cnNzUVNTg0cffVSwn27dumHUqFFYuXKlawMORuciGIUdpVUrLnkRTpiFz+jQpKamYvPmza73fBAvPz8fH3zwgeu4WL373r17Y9++fYr6uvnmm/Hhhx/i4MGDgQvMYESQyA85DEYHQa1WY/bs2R4DDIPRkWAWPoPhB+np6bjzzjvx5ptvIj8/P6B7GAwGyTiAHFVVVejSRXl2SqAEK6cYoZY/XHIKEYzswcpZUlIScFsepvAZDD/p27cvampqAm4f7Krl/Px8FBUVBdxeKeFaXR1q+dtyFXgwsreH1erMpcNgMBidBKbwGQwGo5PAFD6D0cHgV/t2VDqy/B1ZdoApfAajw9HRlU5Hlr8jyw4whc9gMBidBqbwGQwGo5PAFD6DEWHsdjuGDx+O22+/3efcpUuXMHnyZAwZMgQFBQUoK2stufz8889j0KBBGDRoEL788kvX8c2bN2PEiBEYNmwYxo8fj3PnzoVN9pycHAwePBjDhg0TXJdAKcUzzzyD3NxcDBkyxGOV8tKlS5GXl4e8vDwsXbrUdfzAgQMYPHgwcnNz8cwzz4S1lEWo5W9ubsbMmTPRr18/DBw4EAsWtLNK8JTSdvlv5MiRNBi2bt0aVPu2gskZWjqCnN4y/vWvf6X33XcfnTlzps+1d999N12yZAmllNLNmzfTn/3sZ5RSSletWkWnTJlCrVYrbWpqovn5+bShoYFSSmleXh49efIkpZTS9957j86bNy8kcgrRs2dPWlVVJXp+9erVdPr06dThcNDdu3fT0aNHU0oprampob169aI1NTW0traW9urVi9bW1lJKKR01ahTdvXs3dTgcdPr06XTNmjVBy9lW8huNRrplyxZKKaVms5mOHz/eJX9b/TYBFFERvcosfAYjgpSVlWH16tUehd7cOXnyJCZN4mqoT5w4EStWrHAdnzBhAjQaDWJjYzFkyBCsW8fViSeEoLGxEQDQ0NCAbt26tcEnEWbFihV48MEHQQjB2LFjUV9fj4qKCqxfvx5Tp05FSkoKkpOTMXXqVKxbtw4VFRVobGzE2LFjQQjBgw8+6LELWXuXPyYmBhMncnsS6HQ6jBgxwmNWFmmYwmcwIsivf/1r/OUvf4FKJfynOHToUHz33XcAgOXLl8NgMKCmpgZDhw7FunXr0NzcjOrqamzduhWlpaUAgA8++AAzZsxAdnY2Pv3007C6FQghuPXWWzFy5EjXbmDulJeXo3v37q73/K5iUsezs7N9jncU+d2pr6/HDz/8gMmTfbYMiRhM4TMYEWLVqlVIT0/HyJEjRa956623sG3bNgwfPhzbtm1DVlYW1Go1br31VsyYMQM33XQT7rvvPtx4441Qq7mNOf72t79hzZo1KCsrw8MPP4xnn302bJ9h586dOHjwINauXYv33nsP27dvD1tf4SBc8ttsNtx333145pln0Lt375DcMxQwhc9gRIhdu3Zh5cqVyMnJwdy5c7Flyxb87Gc/87imW7du+O6773Do0CHXlo5JSUkAgBdffBGHDx/Gxo0bQSlFnz59UFVVhSNHjmDMmDEAgHvvvRc//vhj2D4Dv0NYeno65syZ41NqOisryzXzAODaVUzquLsLRGwXsvYqP8/8+fORl5eHX//612GTPRCYwmcwIsQbb7yBsrIylJSUYNmyZZg0aRI+++wzj2uqq6vhcDhc1z/yyCMAuMwevoDb0aNHcfToUdx6661ITk5GQ0MDiouLAQAbN25E//79wyK/0WiEwWBwvd6wYQMGDRrkcc2sWbPwySefgFKKPXv2IDExEZmZmZg2bRo2bNiAuro61NXVYcOGDZg2bRoyMzORkJCAPXv2gFKKTz75RHAXsvYqPwC89NJLaGhowDvvvBMWuYOBVctkMNoZCxcuRH5+PmbNmoXCwkK88MILIIRgwoQJeO+99wAAVqsVN998MwAgISEBn332GTQa7s/5/fffx1133QWVSoXk5GR89NFHYZGzsrISc+bMAcC5MO6//35Mnz4d//73vwEATzzxBGbMmIE1a9YgNzcXMTEx+PjjjwEAKSkpePnllzFq1CjXZ05J4baH/Oc//4mHHnoILS0tuO2223Dbbbd1GPnLysrw+uuvo1+/fhgxYgQA4OmnnxYNyrc5Yuk7kf7H0jLbF0zO0NERZKSUyRlqWFomg8FgMNoMpvAZDAajkxCUwieEpBBCNhJCzjr/Txa5zk4IOez8tzKYPhkMBoMRGMFa+AsAbKaU5gHY7HwvRAuldJjz36wg+2QwGAxGAASr8O8EwFc9WgpgdpD3YzAYDEaYCDYtM4NSWuF8fRVAhsh1UYSQIgA2AIsopd8LXUQImQ+A32GgiRByJgjZ0gBUB9G+rWByhpaOIGdHkBFgcoaatpKzp9gJQmVKjxJCNgHoKnDqRQBLKaVJbtfWUUp9/PiEkCxKaTkhpDeALQAmU0rPKxQ+IAghRZRS33qn7QwmZ2jpCHJ2BBkBJmeoaQ9yylr4lNIpYucIIZWEkExKaQUhJBPANZF7lDv/v0AIKQQwHEBYFT6DwWAwPAnWh78SwDzn63kAVnhfQAhJJoTona/TAIwDcDLIfhkMBoPhJ8Eq/EUAphJCzgKY4nwPQkg+IeQD5zX9ARQRQo4A2ArOh98WCt+31mn7hMkZWjqCnB1BRoDJGWoiLqesD5/BYDAY1wdspS2DwWB0EpjCZzAYjE5Ch1b4hJDphJAzhJBzhBCfVb6EED0h5Evn+b2EkJwIyNidELKVEHKSEHKCEPIrgWsKCCENbuUnFra1nE45Sgghx5wyFAmcJ4SQvzuf51FCyIgIyNjX7TkdJoQ0EkJ+7XVNxJ4nIeQjQsg1Qshxt2NKS5DMc15zlhAyT+iaMMr4v4SQ087vdTkhJEmkreRvpA3k/AMhpNztu50h0lZSN7SBnF+6yVhCCDks0rbNnieA9lseWe4fADW41M7eAHQAjgAY4HXNkwD+7Xw9F8CXEZAzE8AI5+t4AMUCchYAWNUOnmkJgDSJ8zMArAVAAIwFsLcd/AauAujZXp4ngAkARgA47nbsLwAWOF8vAPCmQLsUABec/yc7Xye3oYy3AtA4X78pJKOS30gbyPkHAL9V8LuQ1A3hltPr/F8BLIz086S0Y5dHHg3gHKX0AqXUAmAZuFIP7riXfvgGwGRCCGlDGUEpraCUHnS+NgA4BSB8e7aFlzsBfEI59gBIcq6/iBSTAZynlF6KoAweUEq3A6j1OqykBMk0ABsppbWU0joAGwFMbysZKaUbKKU259s9ALJ9GrYxIs9SCUp0Q8iQktOpb+4B8EW4+veHjqzwswCUur0vg68idV3j/DE3AEhtE+kEcLqUhgPYK3D6RkLIEULIWkLIwLaVzAUFsIEQcoBwZS68UfLM25K5EP9Dag/Pk0dJCZL29GwfATeTE0LuN9IWPO10PX0k4h5rT8/yZgCVlNKzIufb9Hl2ZIXfoSCExAH4FsCvKaWNXqcPgnNLDAXwDwDft7F4POMppSMA3AbgKULIhAjJIQshRAdgFoCvBU63l+fpA+Xm8e02F5oQ8iK4mlefi1wS6d/IvwDcAGAYgApw7pL2zH2Qtu7b9Hl2ZIVfDqC72/ts5zHBawghGgCJAGraRDo3CCFacMr+c0rpd97nKaWNlNIm5+s1ALTOVcltCm0tgXENwHJwU2N3lDzztuI2AAcppZXeJ9rL83Sjknd9EfESJBF/toSQhwDcDuAB58Dkg4LfSFihlFZSSu2UUgeA90X6j/izBFw65ycAvhS7pq2fZ0dW+PsB5BFCejmtvbngSj2441764W4AW8R+yOHC6cP7EMApSunbItd05WMLhJDR4L6XNh2YCCGxhJB4/jW4IN5xr8tWAnjQma0zFkCDm6uirRG1nNrD8/RCtgQJgPUAbiVcKZJkcM9/fRvJB0LIdADPAZhFKW0WuUbJbySseMWM5oj0r0Q3tAVTAJymlJYJnYzI82yr6HA4/oHLGikGF5F/0XnsNXA/WgCIAjflPwdgH4DeEZBxPLgp/FEAh53/ZgB4AsATzmueBnACXDbBHgA3RUDO3s7+jzhl4Z+nu5wEwHvO530MQH6EvvdYcAo80e1Yu3ie4AahCgBWcL7jR8HFjTYDOAtgE4AU57X5AD5wa/uI87d6DsDDbSzjOXB+b/43yme3dQOwRuo30sZyfur87R0Fp8QzveV0vvfRDW0pp/P4Ev436XZtxJ4npZSVVmAwGIzOQkd26TAYDAbDD5jCZzAYjE4CU/gMBoPRSWAKn8FgMDoJTOEzGAxGJ4EpfAaDwegkMIXPYDAYnYT/D+yRrp+LGONRAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "time = np.arange(N)*dt\n", "\n", "mask = np.abs(time-5)<3E-3\n", "\n", "\n", "fig = figure()\n", "ax = fig.subplots(1, 1)\n", "\n", "ax.plot(time, amplitude)\n", "ax.set_ylim(-.5, None)\n", "\n", "axins = ax.inset_axes([0.5, 0.10, 0.47, 0.47])\n", "axins.plot(time[mask], amplitude[mask])\n", "ax.indicate_inset_zoom(axins, edgecolor=\"black\")\n", "\n", "for a in [ax, axins]:\n", " a.grid()\n" ] }, { "cell_type": "markdown", "id": "f6579265", "metadata": { "tags": [ "sol" ] }, "source": [ "Il y a a peu près 3 oscillations en 2 ms." ] }, { "cell_type": "code", "execution_count": 5, "id": "4550a298", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "Text(0.5, 0, 'Frequence [Hz]')" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWoAAAERCAYAAABSPe3hAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAAsTAAALEwEAmpwYAAAQzUlEQVR4nO3da5Ac1XnG8edB4mYQIKINRXFbTGESQmIkLxgHimAcy4Ah2FWkEBXHYEgUuzAFqRAXLr7gL7nZIU7wJVGwHEgwF2NToYjNJTEEY4PESixYQtwRQQpGSzAXEYJAvPnQZ6XRana3V5qefVfz/xVb29N95ujtQ8+zPT3dPY4IAQDy2mmqCwAAjI+gBoDkCGoASI6gBoDkCGoASI6gBoDkGgtq24ttr7O9okbbv7E9VH6esP1KU3UBwHTjps6jtn2ipPWSro2IoybxvIskzY2I8xspDACmmcb2qCPiXkkvt86zfZjt220vs/1j27/S5qnnSLq+qboAYLqZ2eV/b5Gkz0bEk7Y/KOkbkk4eWWj7EEmHSvpRl+sCgLS6FtS295T0m5K+a3tk9q6jmi2QdHNEbOxWXQCQXTf3qHeS9EpEHD1OmwWSLuxOOQAwPXTt9LyIeE3Ss7Z/V5Jcef/I8nK8erak+7tVEwBMB02enne9qtA9wvYa2xdI+j1JF9h+WNJKSWe2PGWBpBuC2/kBwBYaOz0PANAZtfaobe9j+2bbj9leZftDTRcGAKjU/TDxbyXdHhFn2d5F0nvGazxnzpzo7+/f3toAoGcsW7bspYjoa7dswqC2vbekEyWdJ0kRsUHShvGe09/fr8HBwclXCgA9yvZzYy2rc+jjUEnDkr5t+yHbV9veo2PVAQDGVSeoZ0qaJ+mbETFX0huSLhvdyPZC24O2B4eHhztcJgD0rjpBvUbSmohYUh7frCq4txARiyJiICIG+vraHmYBAGyDCYM6In4u6XnbR5RZH5H0aKNVAQA2qXvWx0WSritnfDwj6TPNlQQAaFUrqCNiSNJAs6UAANrhq7gAILmeCOoVa1/V0POvTHUZALBNuv3FAVPi9KvukySt/ouPT3ElADB5PbFHDQDTGUENAMkR1ACQHEENAMkR1ACQHEENAMkR1ACQHEENAMkR1ACQHEENAMkR1ACQHEENAMkR1ACQHEENAMkR1ACQHEENAMkR1ACQHEENAMkR1ACQHEENAMkR1ACQHEENAMkR1ACQ3Mw6jWyvlvS6pI2S3omIgSaLAgBsViuoiw9HxEuNVQIAaItDHwCQXN2gDkl32l5me2G7BrYX2h60PTg8PNy5CgGgx9UN6hMiYp6kUyVdaPvE0Q0iYlFEDETEQF9fX0eLBIBeViuoI2Jt+b1O0i2Sjm2yKADAZhMGte09bM8amZY0X9KKpgsDAFTqnPWxn6RbbI+0/05E3N5oVQCATSYM6oh4RtL7u1ALAKANTs8DgOQIagBIjqAGgOQIagBIjqAGgOQIagBIjqAGgOQIagBIjqAGgOQIagBIjqAGgOQIagBIjqAGgOQIagBIjqAGgOQIagBIjqAGgOQIagBIjqAGgOQIagBIjqAGgOQIagBIjqAGgOQIagBIjqAGgORqB7XtGbYfsn1bkwUBALY0mT3qiyWtaqoQAEB7tYLa9oGSPi7p6mbLAQCMVneP+quSviDp3bEa2F5oe9D24PDwcCdqAwCoRlDbPl3SuohYNl67iFgUEQMRMdDX19exAgGg19XZoz5e0u/YXi3pBkkn2/6XRqsCAGwyYVBHxBcj4sCI6Je0QNKPIuJTjVcGAJDEedQAkN7MyTSOiHsk3dNIJQCAttijBoDkCGoASI6gBoDkCGoASI6gBoDkCGoASI6gBoDkCGoASI6gBoDkCGoASI6gBoDkCGoASI6gBoDkCGoASI6gBoDkCGoASI6gBoDkCGoASI6gBoDkCGoASI6gBoDkCGoASI6gBoDkCGoASI6gBoDkJgxq27vZXmr7YdsrbX+pG4UBACoza7R5S9LJEbHe9s6S7rP9w4h4oOHaAACqEdQREZLWl4c7l59osigAwGa1jlHbnmF7SNI6SXdFxJI2bRbaHrQ9ODw83OEyAaB31QrqiNgYEUdLOlDSsbaPatNmUUQMRMRAX19fh8sEgN41qbM+IuIVSXdLOqWRagAAW6lz1kef7X3K9O6SPirpsYbrAgAUdc762F/SNbZnqAr2myLitmbLAgCMqHPWxyOS5nahFgBAG1yZCADJEdQAkBxBDQDJEdQAkBxBDQDJEdQAkBxBDQDJEdQAkBxBDQDJEdQAkBxBDQDJEdQAkBxBDQDJEdQAkBxBDQDJEdQAkBxBDQDJEdQAkBxBDQDJEdQAkBxBDQDJEdQAkBxBDQDJEdQAkBxBDQDJTRjUtg+yfbftR22vtH1xNwoDAFRm1mjzjqQ/iYjltmdJWmb7roh4tOHaAACqsUcdES9ExPIy/bqkVZIOaLowAEBlUseobfdLmitpSZtlC20P2h4cHh7uUHkAgNpBbXtPSd+TdElEvDZ6eUQsioiBiBjo6+vrZI0A0NNqBbXtnVWF9HUR8f1mSwIAtKpz1oclfUvSqoi4svmSAACt6uxRHy/p9yWdbHuo/JzWcF0AgGLC0/Mi4j5J7kItAIA2uDIRAJIjqAEgOYIaAJIjqAEgOYIaAJIjqAEgOYIaAJIjqAEgOYIaAJIjqAEgOYIaAJIjqAEgOYIaAJIjqAEgOYIaAJIjqAEgOYIaAJIjqAEgOYIaAJIjqAEgOYIaAJIjqAEgOYIaAJIjqAEgOYIaAJIjqAEguQmD2vZi2+tsr+hGQQCALdXZo/4nSac0XAcAYAwTBnVE3Cvp5S7UAgBoo2PHqG0vtD1oe3B4eLhT3QJAz+tYUEfEoogYiIiBvr6+TnULAD2Psz4AIDmCGgCSq3N63vWS7pd0hO01ti9oviwAwIiZEzWIiHO6UQgAoD0OfQBAcgQ1ACRHUANAcgQ1ACRHUANAcgQ1ACRHUANAcgQ1ACRHUANAcgQ1ACRHUANAcgQ1ACRHUANAcgQ1ACRHUANAcgQ1ACRHUANAcgQ1ACRHUANAcgQ1ACRHUANAcgQ1ACRHUANAcgQ1ACRHUANAcrWC2vYpth+3/ZTty5ouCgCw2YRBbXuGpK9LOlXSkZLOsX1k04UBACoza7Q5VtJTEfGMJNm+QdKZkh7tdDFnXHWf/u/tjZseh6SIUJQHIendCEVIofI7Nrdpnf9ueVLE5v4/euV/drpkANhk9nt20U2f/VDH+60T1AdIer7l8RpJHxzdyPZCSQsl6eCDD96mYg7r20MbNr67Zb+yyn+yrZ1apl01kGXZKsuqaZeFtvSdJf8lSTp8vz23qS4AqGOv3XZupN86QV1LRCyStEiSBgYGYoLmbX11wdxOlbOFP/vkrzfSLwB0Q50PE9dKOqjl8YFlHgCgC+oE9YOSDrd9qO1dJC2QdGuzZQEARkx46CMi3rH9eUl3SJohaXFErGy8MgCApJrHqCPiB5J+0HAtAIA2uDIRAJIjqAEgOYIaAJIjqAEgOUds07Up43dqD0t6bhufPkfSSx0sZ0fAmGyNMdkaY7K16TQmh0REX7sFjQT19rA9GBEDU11HJozJ1hiTrTEmW9tRxoRDHwCQHEENAMllDOpFU11AQozJ1hiTrTEmW9shxiTdMWoAwJYy7lEDAFoQ1ACQXJqg7rUv0LW92vbPbA/ZHizz9rV9l+0ny+/ZZb5t/10Zm0dsz2vp59zS/knb507V+mwL24ttr7O9omVex8bA9gfKGD9VnuvuruHkjTEmV9heW7aVIduntSz7Ylm/x21/rGV+29dTuV3xkjL/xnLr4tRsH2T7btuP2l5p++Iyv3e2lYiY8h9Vt099WtJ7Je0i6WFJR051XQ2v82pJc0bN+ytJl5XpyyT9ZZk+TdIPVX0L2XGSlpT5+0p6pvyeXaZnT/W6TWIMTpQ0T9KKJsZA0tLS1uW5p071Om/jmFwh6dI2bY8sr5VdJR1aXkMzxns9SbpJ0oIy/feSPjfV61xjTPaXNK9Mz5L0RFn3ntlWsuxRb/oC3YjYIGnkC3R7zZmSrinT10j6RMv8a6PygKR9bO8v6WOS7oqIlyPiF5LuknRKl2veZhFxr6SXR83uyBiUZXtFxANRvRKvbekrrTHGZCxnSrohIt6KiGclPaXqtdT29VT2Ek+WdHN5fuv4phURL0TE8jL9uqRVqr7LtWe2lSxB3e4LdA+Yolq6JSTdaXtZ+WJgSdovIl4o0z+XtF+ZHmt8dsRx69QYHFCmR8+frj5f3sYvHnmLr8mPyS9JeiUi3hk1f9qw3S9prqQl6qFtJUtQ96ITImKepFMlXWj7xNaF5S97T587yRhs8k1Jh0k6WtILkv56SquZIrb3lPQ9SZdExGuty3b0bSVLUPfcF+hGxNrye52kW1S9XX2xvA1T+b2uNB9rfHbEcevUGKwt06PnTzsR8WJEbIyIdyX9o6ptRZr8mPyPqsMAM0fNT8/2zqpC+rqI+H6Z3TPbSpag7qkv0LW9h+1ZI9OS5ktaoWqdRz6JPlfSv5bpWyV9unyafZykV8tbvjskzbc9u7wdnl/mTWcdGYOy7DXbx5Vjs59u6WtaGQmj4pOqthWpGpMFtne1faikw1V9KNb29VT2Ou+WdFZ5fuv4plX+/31L0qqIuLJlUe9sK1P9aebIj6pPap9Q9Wn15VNdT8Pr+l5Vn8Q/LGnlyPqqOob4H5KelPTvkvYt8y3p62VsfiZpoKWv81V9iPSUpM9M9bpNchyuV/VW/m1VxwUv6OQYSBpQFWpPS/qaypW4mX/GGJN/Luv8iKoQ2r+l/eVl/R5Xy5kKY72eyra3tIzVdyXtOtXrXGNMTlB1WOMRSUPl57Re2la4hBwAksty6AMAMAaCGgCSI6gBIDmCGgCSI6gBIDmCGo2xvbHljm9D5fLfac92v+03bQ+1zFs/qs15tr82Th9nlzu13dZgqdhBzJy4CbDN3oyIo9stKBcWOKqr7aajp8datzoi4kbbL0q6tHMlYUfFHjW6puyJPm77WlUXFxxk+09tP1huOPSllraX237C9n22r7d9aZl/j+2BMj3H9uoyPcP2l1v6+qMy/6TynJttP2b7upF7Dds+xvZPbT9se6ntWWP1s53r3fqu4k3bv7W9faK3sEeNJu3ecnjgWUl/rOoy53Mj4gHb88vjY1VdTXZruTnVG6ouez5a1Ta6XNKyCf6tC1RdKnyM7V0l/cT2nWXZXEm/Jum/Jf1E0vG2l0q6UdLZEfGg7b0kvTlWP1HdRrTuukrVPY9vlaSRPW/bZ0j6gqSfTtAXsAWCGk3a4tBHOUb9XFT3CJaqey3Ml/RQebynquCeJemWiPjf8rw6932ZL+k3bI/cx2Lv0tcGSUsjYk3pa0hSv6RXJb0QEQ9KUpS7sZU/Hu36mSioR6/reaouSx55fLikL0v6cES8XWN9gE0IanTbGy3TlvTnEfEPrQ1sXzLO89/R5kN2u43q66KI2OKmVLZPkvRWy6yNGn+7b9vP9nB1e86bJP1hbL5/MlAbx6gxle6QdH4JMtk+wPYvS7pX0ids717uMnhGy3NWS/pAmT5rVF+fc3U7TNl+n6s7E47lcUn72z6mtJ/l6vafk+2njsWSvh0RP97OftCj2KPGlImIO23/qqT7y+d76yV9KiKW275R1d0F16m6beeIr0i6ydW34vxby/yrVR3SWF4+LBzWOF+nFBEbbJ8t6Srbu6s6Pv3bk+1nIrYPUfUH5X22zy+z/yAiBre1T/Qe7p6H9GxfIWl9RHxlqmuRNh1rvy0ijtrOfk5S9aW1p3egLOzAOPQBTN5GSXuPOstjUsre/Dck/aJTRWHHxR41ACTHHjUAJEdQA0ByBDUAJEdQA0ByBDUAJPf/nJYYgQNObvQAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "frequency = rfft(amplitude)\n", "freqs = rfftfreq(N, dt)\n", "\n", "fig = figure()\n", "ax = fig.subplots(1, 1)\n", "\n", "ax.plot(freqs, np.abs(frequency)**2, '')\n", "\n", "ax.set_xlabel('Frequence [Hz]')" ] }, { "cell_type": "code", "execution_count": 6, "id": "5b9478d5", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "Text(0.5, 0, 'Frequence [Hz]')" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYAAAAERCAYAAABy/XBZAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAAsTAAALEwEAmpwYAAA0IklEQVR4nO2deZhcdZnvv2+t3enurN0JWYAk7ItCoEEFQQaEQRSRUUe8ozOMC6N3xnG5juLFbbw+o47OuI3jiNuIKwqKiMrIFiUQErKTkJA9pLN1d5Leu6u66rz3j7PUqVPn9FJdVaer6vt5nn6q6tSpc94+Xf37nnf5vT9RVRBCCKk/ImEbQAghJBwoAIQQUqdQAAghpE6hABBCSJ1CASCEkDqFAkAIIXVKxQRARL4nIp0isnUC+35ZRDZZPztFpKcCJhJCSF0hlZoHICJXAxgAcI+qXjiJz70PwApVfUfZjCOEkDqkYh6Aqv4JwAn3NhE5Q0QeFpH1IvKkiJzr89G3AvhpRYwkhJA6Ihby+e8G8B5V3SUiLwPwnwCutd8UkdMBLAPweEj2EUJIzRKaAIhIM4ArAPxCROzNSc9utwG4T1WzlbSNEELqgTA9gAiAHlW9eIx9bgPw95UxhxBC6ovQykBVtQ/APhF5MwCIyUX2+1Y+YA6A1SGZSAghNU0ly0B/CnMwP0dEOkTknQD+CsA7RWQzgG0AbnF95DYAP1O2KyWEkLJQsTJQQggh04uiPAARmS0i94nIDhHZLiKvKLVhhBBCykuxSeCvAnhYVd8kIgkAM8baubW1VZcuXVrkqQghpP5Yv359t6q2lfMckxYAEZkF4GoAtwOAqqYBpMf6zNKlS7Fu3bpi7COEkLpERA6U+xzFhICWAegC8H0R2Sgi3xGRphLbRQghpMwUIwAxAJcA+KaqrgAwCOBO704icoeIrBORdV1dXVM0kxBCSKkpRgA6AHSo6hrr9X0wBSEPVb1bVdtVtb2traxhLEIIIUUwaQFQ1aMADorIOdam6wA8X1KrCCGElJ1iq4DeB+DHVgXQXgB/WzqTCCGEVIKiBEBVNwFoL60phBBCKgmXhCSEkDqFAkDIGKza1Y393YNhm0FIWaAAEDIGH/7FZnxn1d6wzSCkLFAACBmDjGEgk2XDRFKbUAAIGQNDATbMJbUKBYCQMVBVKKgApDahABAyBgrTCyCkFqEAEDIGyhAQqWEoAISMgaEKrppHahUKACFjoSjIAGzp6MGff/lPGExlQjGJkFJBASBkDBQo8ABeONqPF471o3sgFY5RhJQICgAhY6CqBUlg+2WW2WFS5VAACBkDwycEZHsEHP9JtUMBIGQMFArDEwKyXzI5TKodCgAhY6CKAhfAvvOnB0CqHQoAIWNgTgTzeACWIjAHQKodCgAhY6CqBRPBch4ABYBUNxQAQsZAFYW9gKyBn+M/qXYoAISMgV8vIHoApFagABAyBn4hILv6J0sBIFUOBYCQMfCbCWywDJTUCBQAQsZA/SaCWY8sAiLVDgWAkADUSfZ6J4JZM4GpAKTKoQAQEoAGTPiytzMHQKodCgAhAajn0cZgGSipESgAhARgBIWAPO8TUq3EivmQiOwH0A8gCyCjqu2lNIqQ6UCu6Vv+dnvgZysIUu0UJQAWf6aq3SWzhJBphj0D2DsTOEgYCKk2GAIiJAAnCWx4t9vrAVABSHVTrAAogD+IyHoRucNvBxG5Q0TWici6rq6u4i0kJCScO/0AD4ARIFLtFCsAr1TVSwC8BsDfi8jV3h1U9W5VbVfV9ra2tikZSUgY2AN/UC8g5gBItVOUAKjqIeuxE8CvAFxeSqMImQ5oQB2okxtgCIhUOZMWABFpEpEW+zmAGwBsLbVhhIRNULknVwQjtUIxVUALAPxKROzP/0RVHy6pVYRMA5x5AN43mAQmNcKkBUBV9wK4qAy2EDKtCFr8nesBkFqBZaCEBBHUCwj0AEhtQAEgJIDcRLB8HA/AMz+AkGqDAkBIAEEhIGUIiNQIFABS9xzpHcYnf70VmWz+LX1Q10/OBCa1AgWAhM7+7kFs6egJ7fwfvf853LP6AJ7eczxve24aQFA30PLbRkg5oQCQ0PnKozvxkfu2hHZ+DSj3DOoFZK8ERg+AVDsUABI6qYyBdCa8jKo1p8Wn73+AMFiPXBKSVDsUABI6hmqoyyuK9RjkARTOA/DvEURItUEBIKFjaLjhFAlQgKC+/6wCIrUCBYCEjqqGWlNvj//eAT14QRh6AKQ2oACQ0DE03M6auRxA/vagvv/MAZBagQJAQsdQDfVuOigHELQovMF5AKRGoACQ0DEU4SaBLQUImvE7Uc+AkGqDAkBCR1WnRwgo4P3AXkD0AEiVQwEgoTNtQkAT7vljhYDoApAqhwJAQscwpkkZqIegXkB2xRLHf1LtUABI6BiqoS6wLgioAnIe/ctDGQIi1Q4FgISOauHgW0lsD6Cg3NOu9vH2AgqYIUxItUEBIKFj5gDCG0wjThLYv+unF9vUMCuXCCkFFAASOmELAJwy0PzNQUlgzgQmtQIFgISO2QsobCv8msEFLAhjPTIHQKodCgAJHbMX0DQIARX0ArIfA2YCTwfVImQKUABI6ITeDdR6nHAvIM4EJjUCBYCETugTwewcQEC5Z8E8APYCIjUCBYCETthllUEegF3+qaoYTmfxl/+1GtuP9OVCQxz/SZVTtACISFRENorIQ6U0iNQf9sAf1mSwwHbQriUhN754Emv3n8BnfvN86PYSUiqm4gG8H8D2UhlC6pewl1gMXBDG5ZnEY+a/ymjW4IpgpGYoSgBEZAmA1wL4TmnNIfVI6N01nZnA/uc3FEhEcwIQtmARUiqK9QC+AuAjAAIX8hORO0RknYis6+rqKvI0pB4IarpWKewy0GxBy4fcgjCxqLlPOquBi8UTUm1MWgBE5HUAOlV1/Vj7qerdqtququ1tbW1FG0hqn7BDKnYIyNvawb0gjC0Spgdg7U8XgFQ5xXgAVwJ4vYjsB/AzANeKyI9KahWpK+yBP6zeOoErgrke7bfSGcN5h+M/qXYmLQCq+jFVXaKqSwHcBuBxVX1byS0jdYMTagkMKJYXux20944+1/Mn16vI7QEwBESqHc4DIKGTW2BF8eSuLuztGqjo+W0PwCsAhisElBMAzRMGQqqZ2FQ+rKorAawsiSWkbnEPqG//7loAwP7Pv7Zi57fnAQQt/ajIzVROZ7K5HADHf1Ll0AMgoZMbUMPNAYzV8yfPA4C9nQpAqhsKAAmdsMtAnSogbw7A9cTu/GlOBMuVhxJSzVAASOiEPRHM8QAKksDWdlezuoyRmwfAMlBS7VAASOjYd9KZkILqThWQR4AczwT54sSZwKRWoACQ0LEH1JHRbMh25L92z/h1CwBnApNagQJAQsceeEdGw5kIYKd1C0JArglfhss0egCkVqAAkNCxB9ThkDyAwCok10v3e/Yz5gBItUMBIKGjjgcQjgA4yd6AiWDe9zgRjNQKFAASOmHnAIIGdPcSkVnDLwdQftsIKScUABI6jgBkwskBOM3oPKd3D/AZlwBwTWBSK1AASOg4SeB0yCGggG6gQL4HYFerMgdAqh0KAAkdOwSTyoScBC7IAbgH/dzzlBWqogNAqh0KAAmdaVMGOlYVkKsO1M5VMAREqh0KAAmd8JPA+XY4210K4J6lbJerhtW8jpBSQQEgoaKa660zEloIyE4C+/cC8r43nLY9gPLbRkg5oQCQUHEPssPpkEJArrbPbtyv3Xf7dqiKrSBItUMBIKHiDruE7QEUdgP1nweQtupFmQMg1Q4FgISKe8y1cwCxiATsXR6c1g5jlIH6dSr1zhsgpNqgAJBQMfLKK80RNVJpAZhADsDvbp8hIFLtUABIqOTlACwPICqVFQC7wrNwPHdVAflkfBkCItUOBYCESl4OILQQkL8HkJcE9hWAsppFSNmhAJBQ8ROACjsAge2gg8pAnc9RAUiVQwEgoZKfBLbLKytrQ1A7aGUIiNQ4FAASKupTBuo32FbCBkMV2w734sHNh63tuX387vbpAJBqhwJAQiXPA0iH02IhVwYK/OiZF/GZ3zxv2UYPgNQ2kxYAEWkQkbUisllEtonIP5fDMFIf5E8EM0NAlW6z7J4INpo18hq/2fhtYw6AVDuxIj6TAnCtqg6ISBzAKhH5vao+U2LbSB3glwTOGgpVhVQoG+xuB23/AMELwng/R0i1MmkPQE0GrJdx64f/CqQoNC8JnGsFUcnB1Z0DGM0azrndSWD/HAC/9qS6KSoHICJREdkEoBPAI6q6xmefO0RknYis6+rqmqKZpFZxD6JG3h135fosuNtBuz0AtwnMAZBapCgBUNWsql4MYAmAy0XkQp997lbVdlVtb2trm6KZpFYJutOvZB7APRFsNKtOEjpoSUgbhoBItTOlKiBV7QHwBIAbS2INqTuCEqmVLAW17/SzaiZ77ZBQUDdQ53P0AEiVU0wVUJuIzLaeNwK4HsCOEttF6oSgMbSSFTaGa8DPuJPArn04E5jUIsVUAS0E8AMRicIUkJ+r6kOlNYvUC0F30ZX0AJx5AIYik1UYaq9UNt48gAoZSEiZmLQAqOoWACvKYAupQ4IEoKI5AFc7aDv5bIqAyx4fOxkCItUOZwKTUJkWSWBXFZB9p581ND8E5LMgDAWAVDsUABIqQYuqVFIAnJnAmlv5y3AtVg8wBERqEwoACRXvIBqPmrN/K1oF5OoGmu8BuKuADMc253NUAFLlUABIqHjDKA3xKIBKzwMwyaoiYy30m1XNXxBGgVgk/9+FISBS7VAASKhMCwHISwLnGsO5Y0CGoYhGJG+xGjoApNqhAJBQ8d5EJ2PmVzKMVhCqufN6k8AZw4BI4XrFXBieVDMUABIqXg+g0fIAKjj+OzbY8wAAMwTkXRIyGhFEPOsVV7p1NSGlpJiJYISUDO/4aYeAKukB5K0J7CSECxeEiYgUeAAc/0k1QwEgoeL1AOwQUBg5AFVXJ1CPB2AYiogA0YhXAKgApHphCIiEijeG3pgIIwkM55xBE8EyhrlAjT3+244ABYBUMxQAEirecT4ZC6MMtDAHYHh6AWUNRVTE8QBi1iNzAKSaoQCQUPFOpkrG7SqgECaCaW5AN5elzO2T9YSA7EeO/6SaoQCQUPEOoHYVkF/ztXKRtySk0wzOOxPYDgGZA3/cmhTG2cCkmmESmISKNwfQYHkAfs3XymeD+ZjJqisfULgofCQCRNUKAYXQsoKQUkMPgIRKUA6gsiEg81zpbK70tKAbqJUDsD2AqO0BMAlMqhgKAAmVwIlglQwB+WwzVPNsyFrzAOx2QGE0rSOk1FAASKgU9gIKIwlceC5vEtgsA4XjAdghIOYASDVDASChUtgLyC4DrXwvIDfeJHTWMBARgT0NLB6t/IQ1QkoNBYCESoEH4EwEq5wNfgJgGD7zACKFVUAMAZFqhgJAQqWgF5DTCqKSvYD8V/vyzgMQEdguQG4eAAWAVC8UABIqQesBVPLO2jcEZHgXhDEngjkegJ0ErmC5KiGlhgJAQqVwHoDdDjrcJLDfRDB3DiAWZRkoqX4oACRUCttBVz627ncmvyqgiCsHwF5ApBagAJBQmS5LQiai+f8KWc2fCKYKRCTXBdSuAmISmFQzFAASKoVJ4MoLgKFAUzKav81TBQSY8X/xzgNgCIhUMZMWABE5VUSeEJHnRWSbiLy/HIaR+iCoF1Blk8CKpmR+WyxvCAgw1wN2cgB2GSiTwKSKKaYZXAbA/1HVDSLSAmC9iDyiqs+X2DZSB0yHEJChQLNHALxJYMAM/3hbQdADINXMpD0AVT2iqhus5/0AtgNYXGrDSH3gLfe31wOodA6g0AMoLA81q4DsEBBnApPqZ0o5ABFZCmAFgDU+790hIutEZF1XV9dUTkNqGO8ddDwSQUQqvySk1wPIqhbkJyIROEtCxlkFRGqAogVARJoB3A/gA6ra531fVe9W1XZVbW9ra5uKjaSG8d5lx2MRxKORvNbMhqE4eGKobDYYqoUhIKMwBBSRXBlQlAJAaoCiBEBE4jAH/x+r6i9LaxKpJ7weQCwimD8ziWN9I862P+7swjVfWonO/hHvx0uCorAKKGtowQSBiGtR+BjLQEkNUEwVkAD4LoDtqvrvpTeJ1BPe8TMejWDRrEYc7hl2tnUNpJA1FH3Do2WyoTAHYKj3/t+aB+DYySQwqX6K8QCuBPB2ANeKyCbr56YS20XqBO8AGo0IFs9uxOGe3N2+HWYp1922KpCI5f8rGKoF7Sgi4p4JzCQwqX4mXQaqqquQuxEiZEp45wEAwKLZjTjaN4JM1kAsGkHGygeUq+benOWb/5XOGoUtIiIRcc0EZg6AVD+cCUxCxW/8XDS7EVlD0dmfAgCMZsvrARiqBXc0WS2cCGa2gsifCUwBINUMBYCEil8MfdHsBgBw8gD2IFuuNQIUhR5AUBWQdyYwBYBUMxQAEipBHgAAHLIEYNQobwjIUIVn/PdtBeHOATghICaBSRVDASCh4pcDWDDT9AA6+8wQkD3wl+tuWzUX2rEx1KcZnCsHwJnApBagAJBQ8QsBNVr9gEZGswBysf9y5ADsQb4gB2D4l4FyPQBSS1AASKj4jZ/xqHmnnXaqf6zHMuQAbP0pqALyTQK7PAAKAKkBKAAkVPw8ABFBIhpBKmMP/JYHUIYcgH3+iMcF8EsCS14VEENApPqhAJBQCcqhJmMRpDP5yd9yDLb2IQuTwIXeiXs9ACaBSS1AASChErT4eyIWdXkA+Z5AKbHv8kUErzo717TQCAgBOb2AWAZKagAKAAkVe/xsjEedDpuAxwNwksDlywGIAN96+6VY/bFrERE7NOStAuJEMFJbFLMiGCElw47Bb/zk9Xnbk7EIUhmrCqiMrSDcSeCGeBQLZzUiGpEx5gGYz+PMAZAagB4ACRW1JmE1xKPOcpCA2ZytMjmAwiRwRMRaEEbzvBKzUkic/SMCfHPlHrzlW6tLbhchlYACQELF8GnEBlgCkPWGgMonAOKaCRCNiFkFpF5hyE8WRyOCdNbAzmP9JbeLkEpAASChYqgWlGACVghoND8JXA4PwD6ieDyArv4Ujg+m88RJXFVAqrlVwUbL1KKCkHJDASCh8NCWwzj7rt9jKJ0taMMA5HsA9gA7mi1DEtg6pNuGiAAPbDqMx3d05glANJLrBaQwy0IB5C1fSUg1QQEgofC53+1AOmugs3/E1wNIRHM5gFw30PKVgUY8oR2/5+4QkNsDSGcM355GhEx3KAAkFOzBXFwdNt0kY1GnCmg0W755AIarCsjGO+jnnrs9gPwEMdcGJtUIBYCEgj1gjmaM4CRwBauAvDkA53kkPwdgJwEMzReKcoSnCCk3FAASCvbiLqmMUdCGAcgXAGdN4HLkAJyJYEEegHt77rV6SkRHM/QASPVBASChYHsAqUw2IASUawY3Ws5WED7toPM8AM/zvCog13upbLbkthFSbigAJBTsHkDpjOGfBK5QCMg+Yn7YJ/d+1PVcXDOBFYpo1B0CogdAqg8KAAkF+25+ZDQ4B5Cq4EQwtwi5W05EPeWhdqjIMPLfG80wB0CqDwoACQV74B0ezTp9ddwkY1GnvNKO/VeqHXTKNZjn5QZcC8Io8hPETAKTaoQCQELBvpsfSGWQiPkJgLktnTVyFUNlSQLnylFtUqO5eH5eQjgiTssIQ9VZFQzIFw1CqgUKAAkFu/pmcDwByBjlbQVhewCube6ZvdG8MlBXqMjTw4geAKlGKAAkVIbS/iGghFsAsuVcFN58zB/Mc+fxzg+wXxuqzpoA3s8QUi1MWgBE5Hsi0ikiW8thEKk//DyAhCUKqYzhDK5lbQcd9J/gOmVU/HsBAfQASHVSjAfw3wBuLLEdpI5J+iWB4zkPwJ40Vo5B1q8dtBv3mr+RSL4H4JajNHMApAqZtACo6p8AnCiDLaRO8fcAzMVh0lmjIvMA/GYje8+5rHWGkyxWzS8XZUdQUo2ULQcgIneIyDoRWdfV1VWu05AawFcArG2pUaMyM4EDFMAtABcunuWaCayO9wAwBESqk7IJgKrerartqtre1tZWrtOQKsTb0yceLRx8c2Wg2fK2g3aSwP7vu0VnfksDrjijFQBw/qJZefYwBESqES4KTyrOiGewTMSiBfvkeQBlrAJyJoIF5QA853ztSxfiijOux5ymRF5+gB4AqUZYBkoqznA6v3FaYowy0GHXpKyydAP1WRDGjT2w337FUmfbnKYEgHxxSLMMlFQhxZSB/hTAagDniEiHiLyz9GaRWmZk1CMAY0wEG3SJRVk8AJ8lId1kDcW+z92ET918fsF77iTweL2Ath/pwyd/vZUrh5FpxaRDQKr61nIYQuoHrwAkxxCAoVTG2VapBWHcZLIaKA7uJPB4VUB/+a3V6B/J4P3XnYV5zcnijCWkxDAERCrOyKg3B+DfDA4wewXZlMIDyGQNpxW1G7+OpEBuLQI/3II0ngfQP2L+HiwXJdMJCgCpOCOZfA/ArwqopcF0Tk8OpZ1tpcgBnHnX7/Hue9Y5rw2fBWHcjBWxKaYM1Ct+hIQJBYBUnMIkcGEV0MyGOESA7v6cAJQqBPTYjk7nuVMG6vpP+Mdrz0Rrc2Lc42SKSAJ7f3dCwoQCQCrORJLAkYhgVmMcXQMpZ1s5F4Rxx/k/dMM5+Nbb28f9rN88gJHRbIGn0js06jwfHqUAkOkDBYBUnL6RTN5rPwEAgDkzEuh2CcBUPQC/EFJuHkA+fqWpBZ81CkNA537iYbzLFWICgAMnBp3nXvEjJEwoAKTi/HFnF+Y2JbB4diOAYAGY1RhHd78pAIloxFkXoFjcE9B6h+27cnseQL4ExGNBWYEcGR8BAICVL3Rh6Z2/xX3rOwAAL54YytlAASDTCAoAqSgjo1k8tv0YbrzwFDQmzNh/wicJDABzZsTRPWDmAJKxSF7dfbHntjncMwzAf0lIAIgF9ofOYYePohFB2qe66Lur9gEATjIERKYpFABSUdbsO4GhdBY3XnCKsxBMkAcwe0bCKZtMxqNTzgG4E7C2APgtCANMLARk29OUMNcv9g7udiJ5wBXyYhKYTCcoAKSi7OkcAACcv2imU/7pVwUEALNnxJ3nyVhkyjkAfw/AfyLYREJAtng0J2MYzRoYTOfnNuZZLSMGUjkPgCEgMp2gAJCKsv/4IFoaYpjXlBjfA2jMlWI2JqJTngfgvkN/aMsRZA0NXBBmIiEgm6ZkDKNZxWAqf3C35WpgJOMsIM8QEJlOUAAIPvHAVvxx59TWbNhxtA8fuW/zuHfp+7oHsay1CSKS8wCCqoCach5AQzwy5RCQPQnrDRcvwpp9J/CrjYecUdrbDG4iISCbJtsDSOV7APbrgVQW86xwECeCkekEBaDOyRqKHz5zAE+4JkcVw8oXuvDzdR041jcy5n77jw9i6bwmAHA8AL+ZwIBZBWSTjJUgB2Ddfb/18tMAAIdODruSwJOvAnrZsrmWbRF09adwtDf/d7fbPwykRjG7MYFENEIPgEwrKAB1jl0O2e+pzZ8sfdZxcuWVhaQyWRw6OYylrfkC4NcMDjCTwDZL5zVNOQdgJ2CbG2JojEfRcXIIb/vuGgCFHsBEQkDfu/0yPPHhaxCLCnYc7Xfq/6PWwQYcDyCD5oYYkvEIk8BkWkEBqHN6rF477kRlMdgC0jNUeBxVxRce3oHHt3fCUHNtXQDjJoFPnWPOE3jftWdi8ZxGZA2dUjtlOwHbGI+ipSGGVbu7nfcKPIAAr8RNUzKGZa1NeGr38bztv3//VXjDxYtyAjCSQXPSFB0mgcl0giuC1Tl2jfpAaooewEiwB3B8MI1vrtzjvF7W2gwA4yaBl7c1Y93HX43W5iS+/tguAGbIKjaBwdkPRwASUTQ3xNBxcth5L+VpUBfUAtqPmy9ahN9sPuy8nmEd3y7/7E9lsGTuDDQmKABkekEPoM5xPIAShYD6fATA6xUs8+QAggQAAFqt3vlRa9CfSh5g2O0BJGN56/ie2dZc9HG/+paL8e6rljmvmxIxNCfj6Hd7AAnTA5hoDsAwtCCpTEipoQDUOfbg3D9lD8AKAQ2nC97rcbV0njMjjllWff94VUBu7DLKwVQmrz+Qm87+EVz1r49j57F+3/ftwbchbt6h22z+1A2YP7NhXBuCiEQEp8xqdF7PSEbRnDQnh73/ZxvR2Z+ycgBRDE+wCuj+DR14xeceY86AlBUKQJ1zskgPIJM18IWHdzgTqoKSwId7hvN64dgJYGD8KiA3jXEzT/CqL65E+2cf9d3niR2dOHhiOC/c5MYuwUzGImhOmgIgArQkpx4JnduUX7FkH//Xm8zQkJkDiGBkggP6poM96BvJ5F07m1Qmizd+82ms3nPc55OETBwKQJ1jD9iTzQHsONqPb67cg4e3HgWQywF4wz3/69vP4FMPbnNe2+EfwFUFFJAEdnPdeQsQkZydfpPC7OP1j/gntEdGs2iMRyEiaE6aA3ZzMoZI0Irwk2BuU/4yj80N8bzXLVblkXcxnCD2Hzc7iB70EYBDJ4ex/sBJrN7TXfAeIZOBAlDn2B7AUDo7qTLLPV1mS4ejVt1/37A5MLs9gEzWwIsnhvJKTJfleQATDwEtmt2Ia89d4Lw+MVgYarK3BZW0DqezTgM6e8Ux91yDqWC3fbBp9ngVzckYGhPRCYd09nebA/+P1hzAN57YnffesT4zBHa4d+w5F4SMBwWgznHfsY/lBZwcTGProV4ncbqny7xDPdI7gtFsrhGaWwC6B9Jwa8pLFs/CK89qdV5PJAns5uOvPQ+vPs8UgTX7TmC31VfIfT7ArDpy223vNzyaRYN1LnuAntlQGgGY4xGAlgaPADTE0BCbWBJ4ZDSLw71maG3lC1348iM786qHOvvNgf9o7whUFdsO907VfFKnUABqmI6TQ2j/7CPY0tETuM9EBeCDP9+E1319Fe68fwsAlwfQO5x3x907PIrRrIGR0azjHQDAgplJ/OZ9r8SK0+Y4285dOBMvXTLLmTg1Hktbm/Dea5ab9ty7Cf903+a89+3k8IsnhrD1UC8+/sBz+MLDO3Db3c8AMAfWBssDsJPAMxuD4/+LZjXg1efNn5Btc2fkC0CTxwNIxqJoSETRcXIY91vrBLjpGxnFG77xFLYe6sWLJ4by1iLOGIqth3KDvD3b+nDvMB55/hhe+7VV2HSwZ0J2EuKG8wBqmMd3dKJ7II01e0/gpUtm++7jrtoZKxG87XAfAGCLNRDZXT2P9I44CeCImAJw293P4GjvCD558/nO5+fMKFxj9/UXLcLrL1o0qd/JLgvNGIrdxwagqk7N/nFLANIZA6/7+ioAZkfRnqFRnBxMOzkAYGIewNMfu27CdtmhJRtvYntkNOtc3//zi824+aJFeZ7Pcx292HSwBytf6MSZ8wtLUjcd7EH7UrP1xNHelPU4gmf2ngAAPLvvBC4+dfaE7SUEoAdQM9x5/xZ89L4teduesma6BpVFGobiaO8I2lrMQdWeDdxxcggfvW+Lk0ztHxlFl7Uy14HjgxgZzWJvtxkCOtY3gh5LABbOasSB40NYf+AkDvUMY29XbilEPwEoBlsAALN01b1mcPdAGqd4yjltD2ff8UEMuwSgxfEAShMC8nLBoln42ltXYO1d1+G915yBGy5YgKtc4a/tR/rQcXIIn35wG1bt6sYu62+042g/vv/UfsxrSmB5m5kvmdeUwEbXHf4xKwQ0lM5i5U6zh9PGgyex6WBPQb6AkLGgB1DFHOkdxoYDPbj4tNn42bMHAQCfvfVCxKNm73y7THCXK1b+g6f349Htx/CDv70cz+w9ju6BNG6/Yin+++n9eP5IP/7uhxuwvLUJa/efQGMiik+//gJnIL/+/AV45Plj+MX6DqQzBl6xfB5W7z2OVbvMTqJnzG/GoZ7c7NrHdxxznrs7e06FpmQMMxJRDFnJ1D2dg/jQvZvxhhWLcXwghSvPbMUrz5qHb/1xL3YczQnfgeODGE5nMSNhfuVtD6BUSWAA+PTN5+OUWTkBsr2bj954LgDgze2n4sozW3HF5x/HhhdP4sdrXsTuzgHsONqHM6yJaA9tOQIA+MwtF+D68xdgx9F+/HLDIazZexyqin/53Xb81toHgPO32fhiD7795F78dssR/MUli7HhQA+iEeDGCxfm2djZb3oNN114CmKT6HhKapOiBEBEbgTwVQBRAN9R1c+X1CoyLpmsgff8cD02d/TihvNz1TFbOnpx6elz8LvnjqBvJIPFsxuxu3MAmayBjKH46mO7cGIwjT/u6sIDGw+hJRnDrSsW47+f3o9frDuI7oGUE0v/yZoXcfsVS7FmnykkN1gC8NVHd6EpEcVbX3YaVu89ji/9YScA4IOvPgtvvGQx0hkD/3TfFjy7/yRakjH0pzIl8wAAoK0liQPHzSqZX286hFW7u7Grsx/H+lJobU7g1hVLcOuKJXjt1550QlebXuzB0d4RnL9oFoDSJ4EB4PYrl427z8JZDVgwM4kHNh7C7s4BRATYfLAXKdes5KZEFG+57FQkY1EsnNWIrr4UfrP5MB7cfBjfftJcZrK1Oen8na46qxVP7up2hOHR7Z344sM7EIkIrjlnPpKxCH6z5QhmNsTwnh+tx8iogdWXn4Z/ufXCSbW8ILXHpAVARKIAvgHgegAdAJ4VkQdV9flSG1eLdJwcQlMihqHRLGIRwbymBI71p9CcjKGrP4UzLLffULOr5MhoFj965gCakjEsmJnEJ3+9DbdcvAgrX+jCtsN9EAH+8PwxZxD44er9WDy7EV9+ZCfOXtCMv3rZ6fjUg9tw5l2/x+LZjTgxmEZDPIL/99Dz2Nc9iHdeuQytVghoS0cu0fjnFyzAY9s7ccNX/uRU/lx7rpkQ7R5I4dYVi3HeKS0AgItPnY2XL5+Hly6Z7SR5v/LoLhzqGcaytiYkYxFc4kr+TpXW5iQ6Tg4jGhHH87FLI+2++wBw3sKZ2Ha4D8lYBD9YfQDxqODN7UsATCwJXA5EBC9bNg8PWr2D/u9N5+Gzv92OjS/2IBYRZAzFBYtmIRnL5RSuOtsMHb3/Z5sQEfO7sbytCafPm4Grz2rDWy47FS//3GPO/p94YKvz/NebDqExEcM//nQjAGBuUwJ/2b4Q96w+gOZk1PRIzmjFyaE0BlMZzGyMY+2+EzjH+tt+d9U+pDMGZjfGccvFi3HuwhaMZg3MSJitNEYyWbQkYxSSKqWYb//lAHar6l4AEJGfAbgFQKAA7DzWj+v+bWVRBk6pAbDCubNKxiMF677anSXVtb/7tfd9dd7X/Nc+Rno/GxFBJAIcPDHs/KP7Ma8pgawqhlJZNMQjGPTU58ejgm88sQenzZ2Bf3vzRXh0+zE8uv0Y/uXWl+Dd96zDA5sO44FNhxGPCr53+2WY39KARCyC689fgEzWwK0rFuOsBc248/7nsGROIz54/dl5tixvbcLe7kG88ZIlmNuUxL3PvggAWDpvBuY1J5GIRZDOGPjQ9Wfj1Lkz8OiHXoXlrU0Fk6nesGIRvvHEHpy/cCY+/8aX+v6uxbJkTiP6hkcxe0Ycz+4/ifbT5+AlS2bh+0/tx3kLZzr7vfXyUzGvKYFfbjyErv4UvvX2S525BAtmNiAeFZw2d0ZJbZsIH77hHEcA3nTpEnz2t9sBAO+8ahl+teEQ7nrteXn7L5zViHNPacGOo/34i0uW4LKlc/Dy5fNwumtS3dJ5M7D/+BAuXzoXa/efwKlzG9GUiOHOXz6HeCSChbMa0NmfwodvOAe3XXYqeoZG8e0n9+HbT+5DPGp+H/2+xzMSUcxujKN7MI3vrNqHZMxcmGfJnEYc6RlBOmtgzow4WhriSGcMpLOG+ZgxIGKG2LKGImOo44XObUrkkuauc7pP7+76OtYY4P7W2SIkQTuUiFqSOplse10ReROAG1X1XdbrtwN4mar+g2e/OwDcAQAzFy2/9G1f/EXxRhb9SbP8TsTqA+P3q0r+ObxfIhnn/dznxXd/e1vWAEazBs6c34z+kVG0NicxIxFFZ38K81uS6BvJYFZjHFs6ehARwczGOEZGs2hKxnDVma2Y2RjHyhc6cfNFi7C7cwBXntmKhngU3QMpHOkZwUuWzMKhnmE819GL5w714Kqz2vDy5fMAIK9Sxub4QAoigrlNCaiaoaGhdBbvuHIZ1uw7jte9dBGyhuJI7zBam5MYHs2itTmJPV0DmDMjgblN44d0htIZJGPRCZd5TpSu/hSG01k0JCL4yZoX8cozW9G+dC4yWcM3rr1u/wkc6hnGLRcvLjhOa3MilLvXtftOIBYVXHLaHDyxoxMKxdVntQXG5fd2DeBo7wguWzbXmT/h5mjvCP60qwuvufAUbDrYg6XzmjCzIY57Vu/H8cE0/uaKpZjXnHBCXqqKPV0DOHB8CGv3n0BLMoYFMxswkMrgnAUt2Hq4F4YCf7FiMebPbEDv8Ch+taED+7oHMSMZw6GTw1gwM4n5LQ3Y3TmAVCaLRCxi/kSjSMQiMFTROzSKWFQQiwhi0QiiEUH3QCqvEZ/7+nv/b/y22+QLht+2qa0f4Ufpjxh8ov9826XrVbW9nKcpmwC4aW9v13Xr1k3JUEIIqSdEpOwCUEwZwCEAp7peL7G2EUIIqSKKEYBnAZwlIstEJAHgNgAPltYsQggh5WbSSWBVzYjIPwD4H5hloN9T1W3jfIwQQsg0o6gaOFX9HYDfldgWQgghFYRTAQkhpE6hABBCSJ1CASCEkDqFAkAIIXXKpCeCFXUSkX4AL5T9RFOjFUA1LLJKO0sL7SwttLN0nKOqLeU8QaU6Yb1Q7hltU0VE1k13GwHaWWpoZ2mhnaVDRMrePoEhIEIIqVMoAIQQUqdUSgDurtB5pkI12AjQzlJDO0sL7SwdZbexIklgQggh0w+GgAghpE6hABBCSJ0yrgCIyPdEpFNEtnq2v09EdojINhH5V9f2j4nIbhF5QUT+3LX9RmvbbhG5M+BcSRG519pnjYgsnegvUgo7ReRUEXlCRJ639n9/wLmuEZFeEdlk/XyyknZa2/eLyHPW+X3LxcTka9bnt4jIJZW0U0TOcV2jTSLSJyIf8DlXRa6niMyz/r4DIvIfnv0vta7nbuuaFSxCVez1LIWNIjJDRH7r2v/zAedaKiLDrmv5XxOxsVR2Wu+ttL4Htg3zA87n+92uhJ0i0uL5bnaLyFd8zlWp63m9iKy3voPrReRa1/5l+25CVcf8AXA1gEsAbHVt+zMAjwJIWq/nW4/nA9gMIAlgGYA9MFtGR63nywEkrH3O9znX/wbwX9bz2wDcO559JbZzIYBLrH1aAOwMsPMaAA9N1LZS22m9tx9A6zjnugnA72GuqPdyAGsqbafrs1EARwGcHuL1bALwSgDvAfAfnuOsta6RWNfsNaW6nqWwEcAMAH9mPU8AeDLAxqXu84R0LVcCaB/nXON+Z8ptp+eY6wFcHeL1XAFgkfX8QgCHKvHdHNcDUNU/ATjh2fxeAJ9X1ZS1T6e1/RYAP1PVlKruA7Ab5iLyzkLyqpoGYC8k7+UWAD+wnt8H4Do/tSuXnap6RFU3WPv2A9gOYDFKSImu50S5BcA9avIMgNkisjAkO68DsEdVD0zC/pLaqaqDqroKwIh7Z+uazFTVZ9T8b7oHwBt8TlfU9SyFjao6pKpPWM/TADbAXI2vZJTCzklQ9He71HaKyNkA5sMU1ZIxSTs3qupha59tABrFjIiU9btZbA7gbABXiRmm+aOIXGZtXwzgoGu/Dmtb0HYvzn6qmgHQC2BekTYWY6eDmOGnFQDWBBz7FSKyWUR+LyIXTMHGYu1UAH+w3MU7Ao470eteTjttbgPw0zGOXYnrGcRimDbbjPv9HGe/ctjoICKzAdwM4LGAXZaJyEbruFcVad9U7fy+FTL5RMBNXKW+mxPBjjYElURW+nq+EcAGSyTK+t0sthVEDMBcmK7GZQB+LiLLizxWOSnKThFpBnA/gA+oap/PLhtghjIGROQmAA8AOKvCdr5SVQ+JGV99RER2WHcc5aTY65kA8HoAHwvYpSLXc4x/8DAoykYRicEU0q+p6l6fXY4AOE1Vj4vIpQAeEJELAr7H5bLzr6zvZgvM/6O3w7xzLSdT+ZvfBtNGPyp6Pa2bny8AuKHI40+KYj2ADgC/tNyNtQAMmM2VghaMn+hC8s5+1hd9FoDjRdpYjJ0QkTjML+2PVfWXfgdV1T5VHbCe/w5AXERaK2mnqtqPnQB+BX/3eaLXvWx2WrwG5h3NMb+DVvB6BnEI+eGUcb+f4+xXDhtt7gawS1W/4vemFVI5bj1fDzO2fnaRNhZlp+u72Q/gJwj3uzkmInIRgJh1rQqo5PUUkSUw/5f/WlX3WPuX9btZrAA8ADOZYcfPEjA76z0I4DYrdrUM5l3cWkx8IfkHAfyN9fxNAB6f4l3bpOy0XNXvAtiuqv8edFAROcV2a0XkcpjXcSpCNVk7m6y7K4hIE8y7ha0+x30QwF9bFQIvB9CrqkcqZafrc2/FGOGfCl5PX6xr0iciL7fs+GsAv/bZtZTXc1I2Wvt9FuZN0QfG2KdNRKLW8+Uw/xZ+nkJZ7BSRmC3e1s3U6xD83RzrO1NWO12M992syPW0wnq/BXCnqj5l71z276aOn8n+KUw3aBSmer3TMvpHMP+wGwBc69r/Lpgq+QJc2WqYWeqd1nt3ubZ/BsDrrecNAH4BMyG0FsDy8ewrpZ0wqwUUwBYAm6yfm6z33gPgPdbzf4CZqNkM4BkAV1TYzuXWuTdbdrivp9tOAfAN6/PPYZzKjDL93ZtgDuazPMcP63ruh5mYG7D2P9/a3m7tvwfAfyA3S37K17MUNsK8o1OYhQn2d/Nd1v6vB/AZ6/kbrWu5yTruzZW8ltbfez3M/6FtAL6KXOWaY+dY35lK/c2t9/YCONdz/IpfTwAfBzDo+ttuQq5CqGzfTbaCIISQOoUzgQkhpE6hABBCSJ1CASCEkDqFAkAIIXUKBYAQQuoUCgCpGCKSlfwOjEvDtqkUiKtjpGvbgGef28XTNdPz/lvE7OT4UBlNJSSPYltBEFIMw6p6sd8b1iQXUVWjsiaVjD1Bv9tEUNV7ReQYgA+XziRCxoYeAAkN6875BRG5B+ZEl1NF5J9E5Fkxe5r/s2vfu0Rkp4isEpGfisiHre0rRaTdet4qIvut51ER+aLrWH9nbb/G+sx9YvZk/7FrFvJlIvK0mE3p1orZM973OFP8vd1e0LCIvGqqxySkGOgBkErS6AqT7APwQZhT6/9GVZ8RkRus15fDnNn4oIhcDXOG5G0ALob5nd0Ac7bpWLwT5nT4y0QkCeApEfmD9d4KABcAOAzgKQBXishaAPcCeIuqPisiMwEMBx1HzRbGE/1dAbMB2IMAYHsKInIzgI8AeHqcYxFSFigApJLkhYCsHMABNfuXA2ZPoxsAbLReN8MUhBYAv1LVIetzfn2kvNwA4KUi8ibr9SzrWGkAa1W1wzrWJpiLfvQCOKKqzwJmgzrr/aDjjCcA3t/1dphT+u3XZwH4IsxFXkYn8PsQUnIoACRsBl3PBcDnVPVb7h3EZxlJFxnkQpkNnmO9T1X/x3OsawCkXJuyGPv/wPc4U0HMduM/B/BunVpzPkKmBHMAZDrxPwDeYQ2QEJHFYq538CcAbxCRRjG7oN7s+sx+AJdaz9/kOdZ7rY6UEJGzxeycGsQLABaKtUCHFf+PFXGcifA9AN9X1ZKuQEXIZKEHQKYNqvoHETkPwGorLzsA4G2qukFE7oXZLbQTZntxmy/BXFTjDpjtdG2+AzO0s8FK8nbBfyk9+9xpEXkLgK+LSCPM+P+rJ3uc8RCR02EK1dki8g5r87tUdV2xxySkWNgNlFQdIvJpAAOq+qWwbQGcXMZDqnrhFI9zDYAPq+rrSmAWIePCEBAhUycLYJan6mdSWN7HfwI4WSqjCBkPegCEEFKn0AMghJA6hQJACCF1CgWAEELqFAoAIYTUKRQAQgipU/4/fWVxy3cHCHgAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "fig = figure()\n", "ax = fig.subplots(1, 1)\n", "\n", "ax.plot(freqs, np.abs(frequency)**2, '')\n", "ax.set_xlim(1600, 1620)\n", "\n", "ax.set_xlabel('Frequence [Hz]')" ] }, { "cell_type": "markdown", "id": "7b61c9d0", "metadata": {}, "source": [ "On observe bien deux fréquences séparée d'environ 2 Hz." ] }, { "cell_type": "code", "execution_count": null, "id": "f3d56d21", "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "celltoolbar": "Tags", "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.8.8" }, "varInspector": { "cols": { "lenName": 16, "lenType": 16, "lenVar": 40 }, "kernels_config": { "python": { "delete_cmd_postfix": "", "delete_cmd_prefix": "del ", "library": "var_list.py", "varRefreshCmd": "print(var_dic_list())" }, "r": { "delete_cmd_postfix": ") ", "delete_cmd_prefix": "rm(", "library": "var_list.r", "varRefreshCmd": "cat(var_dic_list()) " } }, "types_to_exclude": [ "module", "function", "builtin_function_or_method", "instance", "_Feature" ], "window_display": false } }, "nbformat": 4, "nbformat_minor": 5 }