{ "cells": [ { "cell_type": "markdown", "id": "6c3e1560", "metadata": {}, "source": [ "# Machine à laver le linge\n", "\n", "Le fichier {download}`data/machine_a_laver.wav` contient le bruit d'une machine à laver le linge lors de l'essorage. Quelle est la fréquence de rotation du tambour ? " ] }, { "cell_type": "markdown", "id": "501f2f54", "metadata": { "tags": [ "sol" ] }, "source": [ "--------------------------------------------------------------------------------------\n", "\n", "## Solution" ] }, { "cell_type": "code", "execution_count": 6, "id": "75d2d9ad", "metadata": {}, "outputs": [], "source": [ "from scipy.io.wavfile import read\n", "from numpy.fft import rfft, rfftfreq\n", "from matplotlib.pyplot import figure\n", "import numpy as np" ] }, { "cell_type": "code", "execution_count": 7, "id": "0282d641", "metadata": {}, "outputs": [], "source": [ "samplerate, amplitude = read('data/machine_a_laver.wav')" ] }, { "cell_type": "code", "execution_count": 8, "id": "396e39d6", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "44100\n", "414380\n", "9.396371882086168\n" ] } ], "source": [ "print(samplerate)\n", "delta_t = 1/samplerate\n", "N = len(amplitude)\n", "print(N)\n", "T = N*delta_t\n", "print(T)" ] }, { "cell_type": "code", "execution_count": 9, "id": "456bf33f", "metadata": {}, "outputs": [], "source": [ "amplitude_tilde = rfft(amplitude)\n", "freq = rfftfreq(N, delta_t)" ] }, { "cell_type": "code", "execution_count": 10, "id": "fd722ed6", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYoAAAEKCAYAAAAMzhLIAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAAsTAAALEwEAmpwYAAA/70lEQVR4nO3dd5hU5fnw8e892wuwdKVIBykiVQVRVzSKGqIxvlGTmGgsMVETTdRoTCKpmp9JTIzEaJCgibHENCMKKgIigjRBuiC9Ln3ZXuZ5/5iZ3bOzM2fO9LL357q4mFPmnOfZnT33PF2MMSillFLBuJKdAKWUUqlNA4VSSilbGiiUUkrZ0kChlFLKlgYKpZRStjRQKKWUspWd7ATEQ5cuXUzXrl0pKipq2ldZWdm0Hei1dV80nF7H7rxgx/z32207eR0pzaOz85zm0X+fkzwnMo+hzg2Vn2D79G8y9T6vmzZtOmyM6drqoDEm4/6NHTvWzJ8/31hZtwO99j8/Uk6vY3desGN2efLfdvI6UppHZ+c5zaP/Pid5TmQeQ50bKj/B9unfZHjpsxOrPAIrTIBnqlY9KaWUsqWBQimllC0NFEoppWxpoFBKKWVLA4VSSilbGiiUUkrZ0kChlFLKlgYKpZRStjRQKKWUsqWBQimllC0NFEoppWxpoFBKKWVLA4VSSilbGiiUUkrZ0kChlFLKlgYKpZRStjRQKKWUsqWBQimllC0NFEoppWxlJzsBoYjIecCX8aR1mDFmYpKTpJRSbUpSShQiMlNEykRknd/+KSKyWUS2isgDAMaYRcaY24HXgeeSkV6llGrLklX1NAuYYt0hIlnAdOAyYBhwvYgMs5zyJeDviUqgUkopDzHGJOfGIn2B140xI7zbE4BpxphLvdsPAhhjHhGR04AfGWNutbnebcBtAN27dx87Y8YMiouLm45XVFQ0bQd6bd0XDafXsTsv2DH//XbbTl5HSvPo7DynefTf5yTPicxjqHND5SfYPv2bTL3P69SpU1caY8a1OmiMSco/oC+wzrJ9DTDDsn0D8KT39U+AiU6vPXbsWDN//nxjZd0O9Nr//Eg5vY7decGO2eXJf9vJ60hpHp2d5zSP/vuc5DmReQx1bqj8BNunf5Phpc9OrPIIrDABnqkp35gNYIx5ONlpUEqptiqVusfuBXpbtnt59ymllEqiVAoUy4FBItJPRHKB64DXkpwmpZRq85LVPfZFYAkwRET2iMjNxpgG4E5gLrAReMUYsz4Z6VMqUzS6DdPnb6WitiHZSVFpLCltFMaY64PsfwN4I8HJUSpjvbF2P4/N3cyBEzX87KoRyU6OSlOpVPWklIqx2gY3AJV1WqJQkdNAoVRbkJzhUipDaKBQSillSwOFUm2BJDsBKp1poFCqLQij6sntNry94aBvhoRW5qw7QFl5TYwSptKBBgql/Ax+6E3+uqE22clImlkf7ODW51fw2pp9rY7VNjRy+99Wcv2flyYhZSpZNFAo5aeu0c28XRnWS8iv6mnXkSq2Hm8MeOre49UAHDrZOlj6Chl7jlXHNHkqtWmgUCrD1DY0cjBE1dD5j83n50sDn5OkCaVVCtNAoVSa2X20iuUHgpd4vv3iR5z9y3kRX9/tjRQiwVvANZa0LRoolEozl/9+EdNXB29Dmbv+YEzuox2llI8GCqXSzMkEzdtkU6BQbYwGCqUyVLDurU7fZxsntO6pTdFAoVRbEMaDXWOA8qeBQinVgq8gYteYrQ0YbYsGCqUyVIuapyAP9lW7jtH3gdms23ui+X34ej3ZX9/tNjpCu43QQKFUWxCkPuntDZ4eUgs/OdR8qpMShYHfvfMJZ/1yHkdr3LFKZZuy60hVxO1IiaaBQinVwpo9x4HAhRBr7FjgDS7Ha5sfdit2HGX74co4pi4z7DjRyPmPzee5D3YkOymOJGWFO6VU9IwxzgfFOWxT2HG4knV7yyNO0zV/WgLArClFEV+jLThQ5fntrNh5jL49kpwYB7REoVSaCqvWwuG5x6rqml7b1zylR5WJig0NFEqlqXg/qiVEMSRNqtdTWrr8CDVQKGWRLo2LEDqtLY47rHqyXtFl2ztWLK9VuNLtZ6aBQimLNIoT4X0bjSBfvh5RgS+XRj8oFTUNFEpZuNMoUsQiqTX1LdeksH7T3XqoIvobKHtp8nHTXk9KWST673bPsSq2HYqsO2mob/VO8vKXxTuCvidU9YiWKiLn+9mmy89QSxRKWSS6RHHxbxfy1ZnLAh5bufMoX3x6CQ3uwGlyktRo6sL9u97uOFzJb976JKb3ALjjhVW262sky90vfcT5/zc/PhdPs0YKLVEoZRGrOFF2soY9J0OPWK6pD37O/a9+zKeHKrmqZ0HA407SGk12/J9lX/vLMnYeqXJ8b6dmr90PwH3Xxe6asfCf1a3XDG+rNFAoFQeTHp1PXaObr0yN3z1CVj1F8DC3Bgf/cRT1Dc6n6qiqT48qlWRLlyYxrXpSyiJWVU91jbGb/yhYioIl1fqAD7eGo0Ubhc2Iu1ATBv50SXWYd05tx2vdPPDPj6kLI1jaSbOap9QPFCJSKiKLRORPIlKa7PSozJYu3/Cg+aF+pKKWaa+tpz5AcIpl1VOL64a4sG+KimTZeaSS45ZR5tF6YWMdLy3fzVsbDsTsmukkKYFCRGaKSJmIrPPbP0VENovIVhF5wLvbABVAPrAn0WlVbUt6dY/1pPXnszcy64MdvLmu5UMs2h41/qWG9PnJwAWPLWDybxbG7HpNM+rGuCyQLh+3ZJUoZgFTrDtEJAuYDlwGDAOuF5FhwCJjzGXA94GfJDidqo0J0sEoKWyn+ab5we3rFRVopHZUvZ4s795b4Wb/ifRae+JoZexKFLG2s9xT+lv86eEkp8SZpAQKY8x7wFG/3WcBW40x24wxdcBLwJXGGF95+hiQl8BkqjaoIYZtCzETJHj54oLdVB5RVT1ZoszivS27rxoc9royhgWby9JqapRQjDFRt1VsOuoZ6HiyJvW6BQciyfoFikhf4HVjzAjv9jXAFGPMLd7tG4CzgXeBS4ES4CljzIIg17sNuA2ge/fuY2fMmEFxcXHT8YqKiqbtQK+t+6Lh9Dp25wU75r/fbtvJ60hlch6P1bi5Z4GnIdbJVNnB7nXjnMqg17C+x3fek5NMq/T/YFEV+yoND402DOrefOzO9z1P8AfOyuf0Tln8cXUNyw40cvvIPM7pkc1NcyoxwDOfKWTZgQZmrK1jYo9sbhvZ/D3Ld1+fLwzKYeqAXLYeb+TnSz0lh97tXPzsXE/X3L+treCdvc2RwyWe4zvL3Xx/lGHoKcWtrtu1QJg6IIeZ6+q4eUQuo0tqW/0uq12F3Pee5+f95ORCnlpTw5cHNNKjU3R/k9aff6SfV+s1fre8gtVHhG+NyuNoteGlzXU8MbkQV11lRJ/Xny6uYNtJz8/T+ruPVKz+JqdOnbrSGDPO/1jKd481xvwL+JeD854BngEYN26cKS4uprS0tOn4ggULmrYDvbbui4bT69idF+yY/367bSevI5XJedx9tAoWeAZZRZXHObODXqPFe7znWT+vvuOFqxZCZQWFhYUtjoHnAfboshp2PHoFr+5bBQf2M3TYUEpH9UTmzsYYOP/88zn58X5Yu4ZTunentHRUq/T59O/fn9LSgbTbeQyWfmBJ03kAvLL5LaC+6XwR8TxwysspKCholR+AQ9WG47ldgb2YDj0pLi5r9bt0dR4CrAJgV25v1h/ZzPvtc/nj1VH+TVp+/hF/Xi3XePKjOUAjI4YP58l3twJ19B8+hsNbPoro8/qLpW8CnlKJ/7MqErH6mwwmlXo97QV6W7Z7efcplTCBeg4lS7jtC44XMXLIOnus//udpu1fqzx/wjMXb2dneWOr43YVGrWNhh/+Zy0na+qDn5SmQnUvTjWpFCiWA4NEpJ+I5ALXAa8lOU2qjalvTL269ESmqMUYDPsls1s95GsbWgcCKycj1a3m7arnb0t38ccFn7bY/+rKPTHt+qpCS1b32BeBJcAQEdkjIjcbYxqAO4G5wEbgFWPM+mSkT7VdqVSiiAXfwz6Stkhfr6ea+kbbb/6+e/zgX+uCnxSEXRde3z2tXZY/PVTBvf9Yw10vfhT2vVJJmhUoktNGYYy5Psj+N4A3EpwcpZqkY6CIZYmjxVpHAmv3nGDqk+/TvTD0o23ZjiNR3c+JWu/cWIdO1oZ9r2ikXjkzsVKp6kmppEulqqdo67FbPvQ9F3t7w0H6PjA7yDta+njPCf672tPGcNBvpHWsekvGomSwevfxhHW/TbeSQKxooFDKItklCk81j99DOcxrBGrU/vdHe5n5/nZufX5F0Pc1NLp5/O2W04jPWZ+4KSv8Swm+4PT0wm38+b1tAYPBgs1lXDV9Mc8v2RnXtO3yDpBbsi38UlMg2pitVBoLN1D85INqnvJrbI3U4YpaTv/RHObsiGwQVtNiOJYHqvXZ+oZ3Ou9gXluzj/e3Rj5SONrpLfwXUfI9nAF+8cZGlnx6pKlNY9OBk4C3OzOwpexkVPcO5VC15752y8OGI83ihH0bhYg46XV01BhzY2ySo1RyhVv1tL3cza/mbOKbpQOivvf+456Bbkv3RxYo/FNe1+hmz7HqoMf9HauKrhvqLu9DO5hQU5KE8vPZG9mwv7zFvtSpKMxsoRqzhwK32BwXPPMzKZURkjmFh/9zNNJv6CICxvDj/67jv5bFd+zq8V9ctqtFUPEJtK/pehGlzjn/3PsHCaApf7GerE+1FCpQPGSMsZ2CUUR0oj6VMazrSBhjHH8L3n20it6dCmOShkjbZf1T+t4nh1oet8mLXUAIxG1g26GKsN7j1O6Tbtxu46h+ZuXOY0Di6vxjdZuMaqMwxrwS6gJOzlEqXVirnmrDmPht2Xb/OS7D51tmNGzBJg30345xz6Bwfj7g/CG79nAjT87fGvJ8t2Wq30Q9dw2BSzbhSrM44awxW0S6isivReQNEXnX9y/eiVMq0axVT7U261n7i8Uj+I6/e+Y8inZNjGAPoXSqz/94z4mQ57y7qazpdWWd/ahwJ45Uu/nnyvCWvGlwm7CqK8tr6tNyJl2nvZ5ewDNauh+eNSF24JlyQ6mMYu31VBNiSoon5m1peh3sjz+Sh8KeCs97IqmeqK5rbFqfItW+tcY6PU+82/zzfzXMB3wgv/ywhu/9Y43tFOLH/Rr8b32rikm/mu/o+sdr3Yyc9hZ/XPBp2rWpOA0UnY0xzwL1xpiFxpivA5PjmC6lkqLOUvVUU28fKH5rGXPgdF3rhZ8cYnVZ6F5NEU25IfCVZz8MmqZ4fZE1Bi553Nlqci8t28WGfdFX3YCzUkc4jteG/gFV+30mDHCg3NmCTsdrPNd/c519N+VU5DRQ+MLofhG5QkRGA53ilCalksZajVATRtWT03aCr81cxu9WhZ5+oi7C3le+xt1AVu8+HtE1Q6lphE8OOmjYFnjgX2u5/IlFjq7rStKX7uNVdfz+nS2hT7R4f4vz8SfGkHrFvRCcBoqfi0gH4HvAvcAM4J64pUqpJLFWPQWaDbWytoGa+sZW1RPBJreLtL3BlW7dYhyw5qjvA7NpDLHu7Nbjyemq/MP/rOPxd5pLi2uCBNjth5sXafrKsx+GLIFapdtv11GgMMa8bow5YYxZZ4y50Bgz1hijU4CrjNOy6qn1g2r4w3O5+LcLOehX3RDsmReL6p5Q1VDBgpR/fXqqsRvbeKQy8kn/5jqcduRwRW3Akfj+1UtXTl8c8P3PvLetxfbpP5pDRV3oX7hx1vM3pUQ8hYeIfDaWCVEqFbSsegr8DTHQmINgz/JoezCFY+HmQ6FPioN4ZHFrWWRjNB7+7zq+8deVIc+rb3Qz7ufv8P1/fhzRfSBwgD5h087RopCYZpEimrmexscsFUqliHoHgQJaPxwbY/y0NKZ5PiOn/hGDnj+RmL3NWcklEc/G5xxODtjgLc6Emv/KTll5ZKWe9OscG0WgMMY8HMuEKJUKrAPuasIYUFZdF7gnU6QlilhNPpcI645ENobhufXxX6XOrqurv3B/U2VhronhW+Fv4/5yTjroYZVKnA64+2qgf/FOnFKJ5rREUVHbMjD8dWngb7KRFjR8g+/A0/4RquE3HfhnYfG+yCY/DKYqQLB+9M1NYV8n2skLgzlgWdNju2Vm3BUHGrjn5dVB37doyyHKk7xuuNMSxXjLv/OAacDn4pQmpZKmoraBotwsAGptAoX/H/buo4HnSnJaorBrsJ62pIYBP0j/hR//vNb5N/BIHtXDfjy31b4tZScTGmTt7lQZpKH7ydW1/PujvQGPHamo5YZnl3HHC6sCHk8UR0uhGmPusm6LSAnwUjwSpFQyHa+qp3v7fLYdrrQdR7H5oLP2g/QvB8ROODO4x+rntv1wJTdvqeb3nfZy5aiebNhXzsGTNZzTr3Pre3pvGovyxP4T1ZzaoaDl9SO4jm8+rVgNUoxUpG0UlXim81Aqo5yorqdb+zwg9MhsJ4zDKnINKC2drIlNtZSvh5qvy+zlTyzipr8Enn3I9zuIZvEm8Ky6N+GRdx1307Xjm2zySGX823PsOCpRiMj/aP45uoBhgM4aqzJOeU09PUvaI4Se68nq7H6BJypoXpGtnNW7jkeVtnScTC5V+C+zaifaqqr13m//a3Yf56z85v2RXPVokgOEj6NAAfza8roB2GmMSU5fPKXiqLy6gfYFOeRmhTeFx/mDuwbc73vmTPmds2kr7Hx6qDLgfo0foS3fEXxqk1gyQV6nO6dtFM5m/FIqzZXX1NOhIIdcV+uqp0Vbgg9o2xhkjYJYDrj716o9Lb6hqtjYsK+cwxWRjwSPp1SZycW2jUJEXg91ASfnKJUOahs8czi1y88mJ0talShueHZZ0Pe+/vF+Nu4vZ7ffutENQVpw3W7T4lwn8WTm4u0B94czVkC1dvkTi/jqzOC/23BsO9HIY3M3Bz6YxkWMUCWKSSJiN6eT4GmvUCrtVXsXvynIyaIg2zMBYDgu+72nemnHo1c07WtwB36IPzl/a4tpyv+5JXQ/+Zp6NxsDDG6bZ1nARwWXiEGMf1nX3KbgH/zTOE6EDBRXOrhGarS2KOXQ1rKTPPrmJr44rjcLPjnEL64agYhQ5Q0UhblZFGZLxIOc5lse3MFKFP7VWG9sd3avXy13tvaBau3W51c0vfYF8LCmkm/DbAOFtk2oTPTcBzt5Z2MZ72z0PNAfuOx02ufnNC1AU5CbRWGOcLSyjka3ISvMhRHmb7YEiiAlCq0uSq6Xl++O+z38Jw1MkeaGiEQzKaBSaanBr/ujb+2Hq//4AQCFudkUZnsm5bv+z0vDvr41CPjfq+mccEafqZgLd8LFaG3YV05ZtfMvB/uOV8dkHE+sOO0eq1QGsf+mV5ibRXGuZ69vwFM4rMEhWNVTY5CShkqM/JwEfEf2/uq3HDzpeFU/n4mPvst5g7ow+fRucUhY+Bz/tESkQESGxDMxQe47VET+JCKvisg3E31/lT6ufPJ9HvxX+OsL+D/Kc7JcLaabCLbCWTBbLGspBFoYBxwuHaoywtVPfRDR+xaFsbxqvDmdPXYqsBqY490eFaI3VKjrzRSRMhFZ57d/iohsFpGtIvIAgDFmozHmduCLwLmR3lNlvjV7TvDistB1z616oxiD220Y1bsEgJG9OjCwJKvp+Dsbw+stYw0swaqeVHKFM1I7UgZ4d1d9VNORpEq7htMSxTTgLOA4gDFmNdHN9TQLmGLdISJZwHTgMjxdbq8XkWHeY58DZgPpP4WmSjluA7+as4nV3gd8fk4W9ZYHvNsYhvzwzYiuHazqKVbeSaN1K1LJ3PXx/7kt2nKY5zeE3ynU100bPF9+UoE4mT9GRJYaY84RkY+MMaO9+z42xoyM+MYifYHXjTEjvNsTgGnGmEu92w8CGGMesbxntjHmigCXQ0RuA24D6N69+9gZM2ZQXFzcdLyioqJpO9Br675oOL2O3XnBjvnvt9t28jpSqZrHO9/3fP96cpKxTd8d8yqptPRGnX5RIfe/V9W0b9aUImZ/UsE/tnmud3m/HMfdV/3dOy6f/h1cfGteVeiTVZs2a0oRANtPNPKTJa27QfuOBxKrv8mpU6euNMaM8z/mtDF7vYh8CcgSkUHAt4HIKt6C6wlY6w32AGeLSClwNZCHTYnCGPMM8AzAuHHjTHFxMaWlpU3HFyxY0LQd6LV1XzScXsfuvGDH/PfbbTt5HalUzaNnUmPw/91bbT5wkso577XYN3HiuVTOe7tpu7S0lDe3v41viFBFTgkQ2XrUw0ecwU2zAs9WqpTVjXMqefue8+lU3whLFrc63nv4OAZ0DfyQj9XfZDBOq57uAoYDtcDfgRPA3Q7fGxVjzAJjzLeNMd8wxkxPxD1V5lm39wQHy2ta9W2H5qmorSac2txG8d4nkQUJ0DYKFZ4Ptx8NOp3LRb9ZyMqd4ffCiwWnkwJWAQ95/8XLXqC3ZbuXd59SIYWafO+zf3if3GwXb919fqtjU598v9W+kvzYdJ9sCNLrSalAQk0CuONwFWP7BJ7SPp6c9np627uqnW+7o4i0XncwOsuBQSLST0RygeuAiHtWqbbFyUDnugY3jQmek7teSxQqDJIy/Zxacvq1qYsx5rhvwxhzDIh4JIiIvAgsAYaIyB4RudkY0wDcCcwFNgKvGGPWR3oP1baUVbXspWTtOWL10L/XOr7m+L4do05XIrphqswhkpqTBzoNFG4ROc23ISJ9iCI/xpjrjTGnGmNyjDG9jDHPeve/YYwZbIwZYIz5RaTXV+lr45FGDpwIf+K7VzY3d0N8eXMdQ388h9oAK9Qt3WZfx3vjxL5Nr8Od4ymQn72+IeprqLYjNcsTzgPFQ8D7IvJXEfkb8B7wYPySpdqCmvpGlm470mLfr5bXcM4j83h+yY6wrnVa++aP8oLdngFOkUy8N+1zw5te//yqEbbnpsr0CipzhGqjCLYmSbw5ChTGmDnAGOBl4CVgrDEm1m0Uqo354/ytXPfMUlbsaP0t/8f/XR90PYj6RjeTf7OA1WUN1DY08umhiha9i2q9BQnfnkjXmh7YrZ3t8SkjTonoukoFI4jt53X9vvKIp7+PRjiTAuYBR73vGSYiGGPeC/EepYLa561iWrD5EOP6dmrVQyhYT6Z5G8vYdqiSFyqEv/xiHieqA//huL3Bo97h6Oj/3hHeDDGuUF//lArTqyv3cPdnBtmek4w10h0FChH5FXAtsB7w/TUbPFVQSkXEV2LwTZ1R61dV5Jtg9XhVHfk5WeTnZHG0xs13/7YSgPa5wqcngn+7anAbth2q4E8LP3WUnjO9cz05pWFCxdqyHUf56f9Sr13LaYniKmCIMUa7cKiY2eVdM/p4tach2r9Nod4bKa5+6gMuHX4K359yOuW1zV+nskJUnJaV1zqe3vknlrYJqz9cP5q7Xvwo4DGXruai4iCctTJOVNcnZMp0p3fYBuTEMyGq7WkKFFWeUoF/iaLRW3W0/3gNW7zTctdYOjJlhfhKH84aACWFgT/eednB/0S6FOc5vr5S8XDmT95qWnArnpwGiipgtYg8LSJP+P7FM2Eqs7mNaZp+ec+xasb9/O1WJYqdR6pocBuq6xs5UO6ZZuNkXWwqaNvntyxMB6v3DdZO8o0L+nPeoK7cPUaDhUqu9fvK434Pp1VPr6GjpFUM+WJCXraL2gY3hyvqWo17+OLTSyjt7fmIHjjhqfWssASKjUcjnx5DHDZEB5uBwzeCdlS3bDxToCmVuZzO9fRcvBMSSydrGlhzqIHGjQdpdBvcBtYeaKDi4324Dazf18DRVXtwG9i4p54Dy3axaXc9e5bupF1+Nv27FNOvaxHFebpSbLzUex/A3drnsfuop7TgX/UEsPKAp9RxpLKWugY3ZdWtv+G3y8vmZG0DHQpygvaA8ud2OLVGsCk/Ak0uqFSmctrraRDwCJ4FhfJ9+40x/eOUrqjsOFLJ4ytrYeWKlgdWWxolP17T/Hqdd1qH9S0W3KN7+zwGdC2mf9ci+nfx/D+gazE9SgpiMmq3LfN9U+9SbB8ofAHFGCg7WcOu8tajrdsX5HCytoGeJQWOA8XJIGM0/DkNKEplMqdfmf8CPAw8DlwI3EQY620n2oCuxdx7Tj7jxo3FJYJLhJUrV3DOWeMREVYsX8aEc87GJcKHHy7l3IkTWLpkCeeeO5HjVfVsO1TBp4cq+fRQBdsOVfLa6n2UW5YzzM120b9LUYsA0t8bUJQzDd5v6tbG4kAjqessuw6W17Cz3M2143rz8ormpUsKcz1Tgl80tBsb9se2vrYxSKBI1cnbVBuQquMogAJjzDwREWPMTmCaiKwEfhzHtEWsMDeL/iVZjOxV0rSv7BMXg7p7RtruKXLRp7Pnof5pgYtTOxTQMd9Ft3b5dGuXz+DuLUfkGmM4UlnHtqbg4QkgG/efZO76gy0eJh3yhNM3L6F/12IGdC1qKoX06liopRALX0woKcht2lcXoEHA+pz+56q9VNTD8J7twVJYzPN2D4zm5xusyaJ0SNeIr6lUpnAaKGpFxAVsEZE78awTEf26oWlCROhSnEeX4jzO6tdyLvi6Bje7jlZ5gsfhShZ/vJUqt2HOuv0cq2quBsnNctGnc2FT6aPucD3d9pUzrEf7RGcnJfgCxWVnnMLavSfYe7ya2vrAM776/P3DXQD0LClg+pfGcMffVwGeny1AVhQjpYP1euqsXWBVGiiriu+6J04DxXeAQjxLoP4MmAx8LV6JSie52S4GditmYDdP3Dzd7Ka0dCIAxyrr2HbYU43lK41sLavg3U1l1Dcanl23iCtGnsqPPzuM7u3z7W6TUcpO1vDUGk9PobzsLL4wpid/mL+Vvy7d6ej97fJzGNit+QGe4w0UriSV2OwG5SmVCPe/V80XL4/f9Z32evIt+lshIjcDxcaY+HfeTXMdi3IZW9Sp1YpUDY1uXp2zgIN5vZm+YCsLNx/ie5cM5oZz+pAdarhxBvjuy2vYfdLzDSg3W8jOcmEMLNpy2NH72+VnU5DTvFRpThJ+ZtbCy9Qze2igUBnN6Qp3fxeR9iJSBKwDNojIffFNWubKznJxSpGL71w8iLfuPp8xfTryk/9t4Mrpi5vmPcpk729tDgg5Wa6w2xaK87LJswSKbO8Q7UhniY1EMiZmU8oqkZ93p1/FhnlLEFcBbwL9gBvilai2pG+XIp67aTzTvzSGwxW1fP6Pi3lufS0nquppaHSz43Al8zeVMfP97fx1Qy2zFm8P2mWzrtGwdNuRhH6AwuWfttwsFzmh5uIA/nLj+KbX7fNz6FDQPOVGCmdXqYzgtI0iR0Ry8ASKJ40x9SKif54xIiJcMfJUzh/chd++/QmzFu9gwqPzqGtwt1hnITcL5u3awIfbj/KbL57Z4hrHq+r4v+U1bH17KX/88hgKE50Jh/zHSuRku8hyMLvehZZFgoryPKWJPu1d7Cx3B51mIxx2A+h6lhSw93g1Q7q3Y/NBz4Rtlw7v3uKcAV2L+PRQZdTp8NepKJejlXWhT1RtRiw+7+FyGiieBnYAa4D3vEuhahtFjLXLz+HhqcPpaw6wqaErHQtz6duliP5diujbpYi1yz/g0+w+/PKNjex8qoqbB3seuvuOV/PVmcvYccLNKe3z+fVbm3lodOQfpvpGN9kucTzNRSCvrNjNwUMNlPrtL/cbEOe0RAEwoX9nlmw70tSOU5CggfOLH5gMwDf+uoLNB0/y1JfHMPq0lutp/+2Ws5nwyLt0b5/HwfLapvRVOxvXF1Qqlw5Vcjz+zif89MoRCS1JO23MfgKwTgK4U0QujE+SVJ/2WXytdGSr/SLCLef1Z2C3Yu568SOmLWnAdNvDr+duprK2gXvH59N38HBu/9tKFu/L5aII7n2wvIaLfrOQru3yuGZsL75VOiDsgGGM4f5XPwZg4llHWzTmWwcugqfXWLalRDH5tGze3RX46Trr6+OZO695CZRvjMxju6snDW7Doi2HMcazPOmq7Yc4Xmt48kujufPvsW9kDvTj8P3RCsLDU4dRUdPAMws+ifpeGiaUv7c3HOSnV46goi7KbyFhcNqY/R1vY7aIyLMisgpPF1mVBKVDuvGfO86lMFu49x9rcBvDK7dP4PROWVw6vDtn9i7hP1vrqfGOSzDGMH9zGWUna0Je+9WVe6iobaCkMIfH5m7m4z0nHKerqq6B9ftOsP1wcxXM/9bsb3GO/zKOOVkusi2N2fl+pYv/O7+ARfd7vpPkZWfRPq/5eMd8F9+9ZEiLB/fMG8fzuwsL2fSzKVw8tLl6KNuvwbxdBPN4OfkGJwI3nduPuy4aFNEAwFF+iycZA8M7Oy91qcy337sy5DMLtyXsnk4bs7/ubcy+BOiIpyH70bilSoU0oGsxP55QwHc/M5h/fWsiQ0/1DNwTEe6/dAhHawwvfLiLitoGpq+u5aa/LGfSr+bz/IZa9hyrCnhNt9vwyordnN2vE099eSwAK3YeC5mWytoGquoauPX5FVzxxPtM/s3CpmO7j7a8V6uqp2xXU68lgLysFofpVuiid6fwW1zyc7JaBAf/9SYe/txwNv98CoO7e8a/OJmSozlOtD7XBDhyx6jwB+sV5mbxu9IC3rrnfM91jeG+8QV85Zw+YV9LZa4XPtzZov0y3pwGCt/n/3Lgr8aY9ehKkElXlCN8+6JB9OrY8kF67sAuDOvsYvr8rXx++mJWHmzknosH84UxPVm4u4HJv1nIi8t2tar//nD7UXYeqeLa8b05pUM+PUsKWLnzKBC8rvx4VR1XPLGISb+az+KtR+jbuTktfdq72H2siroGNwe834L8q55ysqTFN++8CL45B3vIW6/b0+9nJHhKKL4A60RT9VLAqifjPdZ8cFDHrNYnhiACJfkuurfzDMAM9ii4eURukCOqLXjo3+scL/EbC04DxUoReQtPoJgrIu1oXjtbpaBrBnl6yxyprOO+8fl85+JBPHL1SB67oICz+3XiwX+t5Z6XV1Nd1zxtxisrdtMuL5vLRpwKwPi+HVmx4xgrdx5l5E/e4umFn7YKGN//58fsO16D2xi6t8/jsf/X3BtrSEcXO45Ucc4j8zjnkXlAsMbs5o9hfhQN1P4PVetD+5sX9Od/d06yeW/ob2d3TR5I13Z5nO03jQtE1kX3qxNalxK+ecFAoHlsSI8OBQHfO6mnToGvEsfpp+1mYBSwzRhTJSKd8cwgq1JU/5IsXrjlbAZ0LWbTR0ub9nfKd/HcTWcxff5WfvvOJ3QpzuOHnx1GWXkNs9fu54vjelHgnY11bN9O/Gf1Pn70n/VU1jbwyJubqG1wM9L7RdntNiz85BDXn9WbOy4cSL3bcKplKpJuhS7qGho42lDXdL5/G0VutivqEoUTU7zBz184dzuzdwnLH7o4ovuLtA4mXxzXm+eXNE9bsv2RyxERFuyForxs/nD9aM7u14kNq5a2KjX5dzB4/a5JfPYP70eUNpUZGt0mbhOP2pYoROR078tR3v/7i8gYoA/Og4xKknMHduGUDq3nkHK5hLsuGsQXxvTi+SU72Xu8mj8t3Eaj23Drec1LjIzv6+kCumF/OXdeOJAze5ew2DKqeu/xamrq3Zx+anu6tfdUVblcwuPXnskfvzyGLgUtP7SVdQ2UVzc0TeIHnsbsHJs2ili7enTPuFzX9wdqHQhodemwU1rtc4n9w3/qmT3oFmAOsG7tWrd9jOjZIWjaZt00PugxlTmCTYkfC6Ee9t8FbgN+E+CYQXs+pbV7PjOY19bs4+H/rmPRlsNcNapn0/TrAIO7taNdfja19W6+MqEP+07U8P6Ww4Dnab71UAVA04SIPp8f3QuAit0bAc+DrexkLRW1DWw5eJJ2+dkc8Q4iy3ZJi+6xkZQo+nbxpPm0MBq9Y/0n1aOkgGlTh3HpiNYBAYK0a0SYin9+cyKffrzM8fmlQ7qFPkkpG7aBwhhzm/d/HTORgXqWFPC1CX3486LtuATunDywxXGXS/jqhD7kZHnW6ujTqZBXy2uobfQ8kD8t8wSKAV0DzzjfrdDF+9+/kJU7j/Gdl1bzyvI9zNtUxjdLB/DUAk9DnIi06J2UG0GJ4gtjenJk52auHNUj9Mlx7IJx47n9gh675bx+vLnuQIt9kQ6Y6t2pkMQ1YyrlfBzFVwP9i3fivPfu7x278Woi7tfWfKt0IB0Kcrh6TC/6dWm9Qt99l57O3RcPBqCP9/ihKs8TbmtZBZ2KculUFLwHTq+OhRTler6PvLx8F93b53HfJUNanGOdMdc6EexLt53jKA8iwpBOWY4GBo70VtH06dyy9BHvUa5n9iqhd6fADdNhXadr8Eh6hk31k8p8UUykEJLTdgZrJWc+cBGwCng+kpuKyEzgs0CZMWaEZf8U4Pd46jZmGGMeNcZsA27WQBEfHYtymfe9C2jnoLuRr+vrQe8iKVvLKhgYpDRhVeQd3LbvRA1fm9Cn1boR1ga4bIE/fWUMQ09tT5/ORSzY5TgrjnxtYl8mDuzStIphNNOUhMMlEpPlU4d2Ch4oiiMYRKgyx6/f2syDlw2Ny7UdlSiMMXdZ/t0KjCG6Fe5mAVOsO0QkC5gOXAYMA64XkWFR3EM51KU4j7zs0HU+fTp5ShRlVYbK2ga2lFUwoFvoj4H1ATb4lHatjlsbs7NcwpQRp7ZoK4klEWm11G08+dYEt5v2w4lo49myhyKZ0EWlk6cXbmPJp0ficu1IV3ypxDPVeESMMe8BR/12nwVsNcZsM8bUAS8BV0Z6DxV7HQpzKCnMYdPRRi55/D3Ka+odrSldaOnK1NlbTWWtQvEvUSRSvCfd+++d53L/lCGICN+7ZHCLY75ZQEf0bM/MS+M732+3dm1nBcW27PG3o59fLBBx8ociIv+juaOIC883/leMMQ9EfGORvsDrvqonEbkGmGKMucW7fQNwNvAw8AvgM3iqox4Jcr3b8PTQonv37mNnzJhBcXHzt92Kioqm7UCvrfui4fQ6ducFO+a/327byetI/HRJNdtOuMkWuP+sfAbbjD723etYjZt7FlQD8OBZ+QzplNUiHXtPunlosef4T8cZTusSuzzeOMcz79SsKYFLKE+vqWHJ/kZuG5nHxB7NJZ94/B63HW/kp0ub59v69QUF3Luwmkv7ZjO1V33Q31FFRQX/253D3J0NfL6f4cohns/rne97omq2C2ZcUsSvllWz8WjrcbC+vPt+FlaTe2fz7u7ETS6n4i/YZx1Cf16nTp260hgzzv+Y00rNX1teNwA7jTF7HL43KsaYI8DtDs57BngGYNy4caa4uJjS0tKm4wsWLGjaDvTaui8aTq9jd16wY/777badvI7Evw98xLbV+7hhYl9umzrc9lzfvSpqG2DBXAAumnQWA7u1a5GObYcqYLFnfqj2xYUxzeP1Rz9mYLd2lE4KXAD+z4GPYP8+hg49nVJvt95A9w2VRyfHSnYfh6WLm7avuWwyo8dV0KdTIe8vei/o72jBggX06t0Ndm4nLzev6fPqKdiDy+WitLSUZ7YshaNHmHnjOJZtP9Y0xUPTtebMbpXGmXdcynXPLGHptqNcOrw7c9cfDJlnldr6n3EWp3UOXEIN9XkNxjZQiEg+nof0QGAt8KwxJl5fP/YCvS3bvbz7VAoZ0aMDb63bx7dKB4Y+2avQsmxpp6LWg8Ws4yhiPbL0katbT9ceSLKWfQjWtTgYazKnTR3GtP9taHVOXnYWD1x2OpefcQrvbHD+4LeuQ67SV1V97B/RodoongPG4QkSlxF44F2sLAcGiUg/EckFrgNei+P9VAS+Pqkfv7mgkK4BRgcHY+3lFGjksnX22OxIW80ilKheT/Fw3VmneV4ECXIje5XwXb+uyIHEojeWymyhqp6GGWPOABCRZwHnw0FtiMiLQCnQRUT2AA8bY54VkTuBuXi6x870zlKrUkiWSyjMifzBEqjEYB1wp8su2LP+eGId49I5aKpm8SgdhwoUTTO4GWMaYvVBMsZcH2T/G8AbMbmJShstBtwl+FlVOqQr//5ob1jTjSeT9Rngq7K7eoxn/ipdNVXFS6hAcaaI+NbGFqDAuy2AMcakx1+XSmnW6qg4TX4Z1JWjenLx0O5NgwJTVaAfS5ZLWDvtEgpzs0Oea3ttLUioEELN9aStWyom/vSVMUEH9WW5hNsvGMDstfuSUv2RqCARjzEb7fIDz1arVCwluOlQtVVTRpzKhacHn8X0gctOZ9H9OhlxLEQajrRgkRniUQWpgUKpDBFpYUyrnlQoqV0xq5SKi5dvO4f1a1e32Ne1vfMuz6pt0RKFUm3Q2f0707+Dd8nbPp41wK8Z08vuLSpNRLoglh0NFEolSDR/vvGsHvrORYN457sXMCjErLq/+HzTigCcN6hL/BKkUo4GCqXSQDzHSGS5pNVytv5Gdsniy2f3adqeeeN41k67JH6JUilF2yiUUiHdMdrTfvH4tWdSVl5LTpaLnCz9ntlWaKBQKg04qXryrVKYHYdRi3neIfOfH63tGKkuGVN4KKVSiF1D5aNXj2TMabs5q1+nBKZItQUaKJTKEB2LcvnGBQOSnQyVgbSSUak0kowpwe+80PnaI+GaOKBz3K6tYkcDhVJpJJo+8ucODP+hfNHp3bj30tBrWljdffEgx+f+/dZzwk2SSgINFEolwRPXjw7r/FhMlvjCLeE/lAeE6DYbyN0XDw64//HSAm4PUDVWlOsZ+KdjM1KXBgqlEqx/BxefO7NHspPhyH1hliZ6dSwIeqxjvov7Lh3Cqh99psX+Rd+fzK/OK+CvN59Nz5Lg71fO1DW6Y35NDRRKJUg03RZ7ex/AnfIT+ycb7liJ3BDnZ7mETkW5LfZ1Ksqle5E+imJl55HKmF9Tez0plQa+fHYf+nQuonHvumQnxVawWDj59G5A7B9gKjE0jCuVBlwu4fzBXdN2XesZXx2X7CSoKGigUErFzOVnnBJwvyvRa9y2YfHoQq1VT0qpmPh42iUU58b/kVJSmMPxqvq43yddxWOacQ0USiVYMr9bz7ppfMgG50i1T9D63Wf07MCiLYcTci/loYFCqQQZ3qM9w05tz9V96pKWhtIhwdctT0cdCnI4Ua2lC6t4VD1pG4VSCZKfk8Ub3zmPgSVZyU5KSLO/PYnX7jw34vdH2iRhwuxDHO7AxbZAV7hTSiXE8B4dGNmrJOL3P3Z+Af+4fUJM0tK7U/BBePnZ+ghLBP0pK6Widsukfi22Oxe4GN83NtOdL7p/ctBj2puqNa16UkqlpL5dimJynXArTYoS0MtKaaBQSqWxYT3aJzsJbYKGY6VU1C4Y3DXg/j9/dRz9u8amtOHja+8e37cjAFeMPJXZH++P6T0A8nNc1NTHfoK9dJTyJQoR6S8iz4rIq8lOi1IqsN6dCgPu/8yw7gzo6nyq8n4OqrD8e/X8/tpRjq8fju9+ZjA9OuTH5drpJq6BQkRmikiZiKzz2z9FRDaLyFYRecDuGsaYbcaYm+OZTqVUanjqy2NDnuPfgzY7TgMIBYlDR9P0FO8SxSxginWHiGQB04HLgGHA9SIyTETOEJHX/f5l1uggpZStDoU5TL+okIX3lYb1voKcLL49eSB/+kroQBPM6Z1SvoIlaSTcAS5h30CkL/C6MWaEd3sCMM0Yc6l3+0EAY8wjIa7zqjHmGpvjtwG3AXTv3n3sjBkzKC5uLvJWVFQ0bQd6bd0XDafXsTsv2DH//XbbTl5HSvPo7DynefTf5yTP8crjz5ZUU9to+PmkwpDnAtw4xzN1+KwpRSHPt+7zf1+gPFrP8b0GGNrJxcajbgZ3dPGDs1uOsdhz0s0PF1c7yH1rgzoYtpxo7lp67ZBc3tpRz7Ha9CpX3DYyj4k9Ajc/h/q8Tp06daUxptVUv8lozO4J7LZs7wHODnayiHQGfgGMFpEHgwUUY8wzwDMA48aNM8XFxZSWljYdX7BgQdN2oNfWfdFweh2784Id899vt+3kdaQ0j87Oc5pH/31O8hyvPAa7ZND7zZntfV9pyPNb7PN7X6A8/uDYPNyd+lJ6wYCm8wFKSjrC0SN06NCB0tKJrZL0w8WzW+1zIisrC2huvB4woD95+3dAbU1E10uW4cOGUjqqZ8BjoT6vwaR8rydjzBHg9mSnQynV2rzvXcCq5cvicu3BHbM8QSJMhblZVNU1hv2+7oUuNh3VXk6BJKNSbi/Q27Ldy7tPKZVmBnQtpmthYh8j8ZjLaGSvDpzbs/X35mmfG06uThOSlECxHBgkIv1EJBe4DngtCelQSmUgXxfbcKdT95/4QhCmjDiFpQ9eFKOUpa94d499EVgCDBGRPSJyszGmAbgTmAtsBF4xxqyPZzqUUpkjVP+bHG+AGNOnxPE17WZHytMSRXzbKIwx1wfZ/wbwRjzvrZRq23yT4+VmQagmi2E9OgBlAY8V5WVz1+g8/vBRbYxTmD40VCqlMsrEAZ0BOLN3CQA3DM0Neu7o0zznXDq8u+01exS17Udlyvd6Ukopq1BN2d+7ZAjXn3UaPUsKGHNaCTllG3l2XeBVBe2qnOLRaJ6u2naYVEplnCyX0LtTIS6XcMnwUxBpDgd2iyCp4DRQKKUSbu7d53PnqLyI3jtpYBcA2uXnxDJJrcRjAaB0pVVPSqmEG3JKO/afEv7jZ91PLqUoN4vC3CyuHtMr6nRYSxv9Ouj35mA0UCil0kZxnueRdct5/aO6znNfP4tPDpxk0dbDTfuydVnVoDSEKqXaDF91Ut/Ohdx6vvNgU5IffRA5pSh9A5EGCqVUm1WclwVAtsv+UViQLVw3vneLfUNPDW8Z1kElWeElLoVooFBKtVm//PwZ3HvJ4KaxF3Z84zJ8rhkbfRtJPFjbXWJFA4VSqs0qKczlzsmDcPm1T5zaIZ+pZ/Zose+68b357x3nNm1H8zjuXig8c0PkiyzZ6VIUfIBhpLQxWymV8v5zx7k0umM3BXio+aKWBJgIUETo17V5gabi/MCPz4cuH8qOI5W88OGuoNe/vH8Olww/xVliwxSPbsNaolBKpbxRvUsY26dT1NeJZa3MF8b04isBpgeJQ81PWOJxfw0USqmM9+v/d2bT9ON2fvvFM/mcX5WTlfUZnOUSLu7j/Nv70M7amK2UUinrmrG9mH9vacjzrh7TiyeuHx2XNEzskc3FQ+0nH0xVGiiUUm3GVyf0BaBzcewafIf3cN5Ntms7z7Ql6TaiQgOFUqrNuHlSP3Y8ekVMG3z/fss5rbrOBvP9KUO4cWJfJvRIr35EGiiUUioKHQpzKB3ctcW+yad3C3huSWEu0z43nBxvd1zfehipTgOFUko5lBNkHe5vXTig6bWIcNHQ7txz8eCQ1/vrzWfHLG3N94/5JTVQKKWUU/k5gXsu5WU37/c9p3OyQz+xi/Oy+eWkAv79rYmxSF7cpFdFmVJKJdk/bp/Aa6v3xex6PYpdjD6tY8yuFw8aKJRSKgzj+3ZifN/oB/+lE616UkopZUsDhVJKxVCkjcmv3zUpNvePwygNDRRKKZUCRvTskOwkBKWBQimllC0NFEopFUMTB3SJ+L3+a2BEQsdRKKVUihtySjsArhnTix7F6TarU2AaKJRSKg66tc/nl5MKk52MmNBAoZRSGaRNVj2JyFAR+ZOIvCoi30x2epRSKl6sz/iB3YqTlg5/cQ0UIjJTRMpEZJ3f/ikisllEtorIA3bXMMZsNMbcDnwRONfuXKWUautccShSxLtEMQuYYt0hIlnAdOAyYBhwvYgME5EzROR1v3/dvO/5HDAbeCPO6VVKqaSJxTM+Hs3nYoyJw2UtNxDpC7xujBnh3Z4ATDPGXOrdfhDAGPOIg2vNNsZcEeTYbcBt3s0hwAHghOWUDpbtQK+7AIcdZyw467UjPS/YMf/9dtvBXscin5pHZ+c5zaP/Pid5TmQeQ50bKj/B9unfpEcqfV5LjDFdWx0xxsT1H9AXWGfZvgaYYdm+AXjS5v2lwBPA08AdYdz3mWDbgV4DK2KU32eiPS/YMbs8OcljrPKpeYxtHu3yEuxYIvMYST6d7NO/yfT4vBpjUn/2WGPMAmBBBG/9n812sNex4PR6ducFO2aXJ/9tzWN0EplH/31O8xytcK4Vbj6d7MvE32Um5jG9qp7iSURWGGPGJTMNidAW8ql5zAxtIY+QHvlMRvfY5cAgEeknIrnAdcBrSUiHv2eSnYAEaQv51DxmhraQR0iDfMa1RCEiL+JpY+gCHAQeNsY8KyKXA78DsoCZxphfxC0RSimlohL3qiellFLpLeVHZiullEouDRRKKaVsaaBwQET6i8izIvJqstMSSyJSJCLPicifReTLyU5PvGTq789KRK7y/h5fFpFLkp2eeGgr8755/y5XiMhnk50Wn4wPFDGab2qbMebm+KY0NsLM79XAq8aYW4HPJTyxUQgnn+n0+7MKM4//8f4ebweuTUZ6IxFmHtNy3rcInkHfB15JbCrtZXygIEbzTaWRWTjML9AL2O09rTGBaYyFWTjPZ7qaRfh5/KH3eLqYRRh5TNN532bh/Bn0GWADUJboRNpJ+ZHZ0TLGvOcd9Gd1FrDVGLMNQEReAq70DvpLmeJeJMLJL7AHT7BYTZp9aQgznxsSnLyYCCePIrIReBR40xizKrEpjVy4v0djzGvAayIyG/h7QhMboTDzWAwU4Qke1SLyhjHGncj0BpJWD4cY6knzN2nwPDB7BjtZRDqLyJ+A0b6R5GkmWH7/BXxBRJ4i9lMKJEPAfGbA788q2O/yLuBi4BoRuT0ZCYuhYL/HUhF5QkSeJr1KFIEEzKMx5iFjzN14guCfUyFIQBsoUcSCMeYInrrfjGKMqQRuSnY64i1Tf39Wxpgn8EyembGimPct7RhjZiU7DVZttUSxF+ht2e7l3Zep2kp+20I+NY+ZIa3y2FYDRarONxUvbSW/bSGfmsfMkFZ5zPhAIZ75ppYAQ0Rkj4jcbIxpAO4E5gIbgVeMMeuTmc5YaSv5bQv51DxqHlOFzvWklFLKVsaXKJRSSkVHA4VSSilbGiiUUkrZ0kChlFLKlgYKpZRStjRQKKWUsqWBQmUsEWkUkdWWf32TnaZYEJG+IlItIqst+yr8zrlRRJ60uca13umtX49jUlWG0LmeVCarNsaMCnRARATPOKKUmHQtAp8Gy5sTxpiXReQgcG/skqQylZYoVJvh/Sa+WUSeB9YBvUXkPhFZLiIfi8hPLOc+JCKfiMj7IvKiiNzr3b9ARMZ5X3cRkR3e11ki8pjlWt/w7i/1vudVEdkkIi94gxQiMl5EPhCRNSKyTETaBbtOlPm2lqqqReSCaK+p2hYtUahMVmCpntkO3AMMAr5mjFkqniVDB+FZG0DwrHNwPlCJZ+6dUXj+RlYBK0Pc62bghDFmvIjkAYtF5C3vsdHAcGAfsBg4V0SWAS8D1xpjlotIe6A62HWMMdvDyCtAJ7xzB/lKHiIyFbgf+CDEtZRqQQOFymQtqp68bRQ7jTFLvbsu8f77yLtdjCdwtAP+bYyp8r7PyWRtlwAjReQa73YH77XqgGXGmD3ea60G+gIngP3GmOUAxphy7/Fg1wkVKPzzeiMwzrI9CHgMuNAYU+8gP0o10UCh2ppKy2sBHjHGPG09QUTutnl/A81Vtvl+17rLGDPX71qlQK1lVyP2f3cBrxMNESnGswbzrcaY/bG6rmo7tI1CtWVzga97H6SISE/xrJH+HnCViBSISDtgquU9O4Cx3tfX+F3rmyKS473WYBEpsrn3ZuBUERnvPb+diGRHcB0nZgJ/McYsivI6qo3SEoVqs4wxb4nIUGCJt325AviKMWaViLwMrMGzyP1yy9t+DbwiIrcBsy37Z+CpUlrlbaw+BFxlc+86EbkW+IOIFOBpn7g43OuEIiJ98AS0wSLyde/uW4wxKyK9pmp7dJpxpUIQkWlAhTHm18lOCzS1tbxujBkR5XVKgXuNMZ+NQbJUBtOqJ6XSTyPQwa+XU1i8pZk/AsdilSiVubREoZRSypaWKJRSStnSQKGUUsqWBgqllFK2NFAopZSypYFCKaWULQ0USimlbP1/lGziLu2J/6UAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "fig = figure()\n", "ax = fig.subplots(1, 1)\n", "\n", "ax.loglog(freq, np.abs(amplitude_tilde)**2)\n", "ax.set_ylim(1E-3, None)\n", "\n", "ax.set_xlabel('Frequence [Hz]')\n", "ax.set_ylabel('Puissance [u.a]')\n", "\n", "ax.grid(which='both')" ] }, { "cell_type": "markdown", "id": "b3221149", "metadata": { "tags": [ "sol" ] }, "source": [ "La fréquence typique de rotation est de 1200 tours par minutes, soit environ 20 Hz. Il s'agir du pic de gauche. On zoom dessus." ] }, { "cell_type": "code", "execution_count": 11, "id": "ea109613", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "Text(0, 0.5, 'Puissance [u.a]')" ] }, "execution_count": 11, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAEGCAYAAACKB4k+AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAAsTAAALEwEAmpwYAABIH0lEQVR4nO2dd3zb9Zn4348sec/EiROP2NmDQJYJIxCcQlllttcWSnvQArkOuOv1Omi543rXo+39SnsHhes1ZV8ZHdA2QCgzhhCy995xvGLHM/GULX1+f2hYtmVbsiXLsp/36+VXrK++34+eTyzp+T5bjDEoiqIoCoAl0gIoiqIoIwdVCoqiKIoXVQqKoiiKF1UKiqIoihdVCoqiKIoXa6QFGAqZmZmmoKAg0mIETXNzM0lJSZEWY1gZi3uGsbnvsbhniK59b9u2rcYYM8Hfc1GtFAoKCti6dWukxQia4uJiioqKIi3GsDIW9wxjc99jcc8QXfsWkZK+notK95GI3CgiqxobGyMtiqIoyqgiKpWCMeY1Y8zKtLS0SIuiKL14/PHHmTFjBiKC741LfX09t956KxdccAFLly5l7969fq83xvDggw8ya9Ys5s6dy2OPPRbU9QUFBdTU1HgfFxcXc8MNN4Rwh8poJiqVgqKMZJYtW8a7775Lfn5+t+M//vGPWbhwIbt37+b555/nH/7hH/xe/+yzz1JaWsrBgwc5cOAAt912W1DXK8pQUKWgKCFm0aJF+EuA2L9/P5/4xCcAmDNnDidPnqSqqqrXeb/61a946KGHsFhcH8+JEycGdX1/XH/99SxcuJCFCxeSlpbGc889F9T1yugnKpWCxhSUaGTBggW8+uqrAGzevJmSkhLKysp6nXfs2DF+97vfUVhYyHXXXceRI0eCuh5gxYoV3i//e+65x3t8zZo17Ny5k6eeeor8/HxuueWWEO9SiXaiUiloTEGJFC+88AIFBQVYLBYKCgp44YUXAr72gQceoKGhgYULF/LLX/6SRYsWERMT0+u89vZ24uPj2bp1K/feey9f+cpXgroeYO3atezcuZOdO3fy5JNPdnuupqaGL33pS7z44ovoZ0jpSVSnpCrKcPLCCy+wcuVKWlpaACgpKWHlypUA3HHHHQNen5qayjPPPAO4gslTp05l2rRpvc7Lzc3l05/+NAC33norX/7yl4O6vj8cDge33XYbDz30EPPnzw/qWmVsEJWWgqJEggcffNCrEDy0tLTw4IMPBnR9Q0MDdrsdgCeffJLly5eTmpra67xbbrmFtWvXAvDBBx8wa9asoK7vjwceeIALLrjAG7xWlJ6oUlCUADl16lRAxx977DFyc3MpKyvj7rvv9vr0Dxw4wPz585k9ezZvvvkmjz76qPea66+/noqKCsD1xf3KK69w/vnn8/3vf9/r/unv+kB55JFHePvtt73xhtWrVwe9hjK6kWgeslNYWGi0ojk6GA17LigooKSkdyFofn4+J0+e9HvNaNh3sIzFPUN07VtEthljCv09F5WWgmYfKZHg4YcfJjExsduxxMREHn744QhJpCihJyqVgmYfKZHgjjvuYNWqVSSOmwQI47OyWbVqVUBBZkWJFjT7SFGC4I477uDF2insqzjL3ZdN5Y4b5kVaJEUJKVFpKShKJGlo6QCgrL5lgDMVJfpQpaAoQdLQ4koLLa1rjbAkihJ6VCkoShDYO5002x2AWgrK6ESVgqIEQUOry0qYMi6Rs22dNLZ2RFgiRQktqhQUJQg88YTzc1yZb2otKKONEaMURMQiIg+LyC9F5M5Iy6Mo/qhvdlkK5+W42kuU1WtcQRldhFUpiMjTIlItInt7HL9WRA6JyFERecB9+GYgF+gA/PcDVpQIU++2FC7ISQdUKSijj3BbCs8C1/oeEJEY4AngOmAecLuIzANmAx8bY74FfC3McinKoGh0xxSmTkgiKTZG3UfKqCOsxWvGmA9FpKDH4aXAUWPMcQAReRmXlVAK2N3nOPpaU0RWAisBsrKyKC4uDrHU4aepqSkq5R4Ko2XP24673qJ7t24kI9bJziOlFBef6fP80bLvYBiLe4bRs+9IVDTn4FIAHsqAi4BHgV+KyOXAh31dbIxZJSKVwI0pKSlLoqUBlS/R1DgrVIyWPW9oPUDssZNcfWURvy/bSnlDK0VFy/s8f7TsOxjG4p5h9Ox7xASajTEtxpi7jTH3G2OeGOBc7X2kRISG5g7SE22ICHnjEimvbyWaOw0rSk8ioRTKgTyfx7nuYwGjXVKVSNHQaicjMRaA3IwEzrVrrYIyuoiEUtgCzBSRqSISC9wGBDXpQy0FJVLUt3SQlmgDXEoBNANJGV2EOyX1JWADMFtEykTkbmNMJ3Af8BZwAPi9MWZfkOuqpaBEhIYWOxlepeCaraAZSMpoItzZR7f3cXwNsGYI674GvFZYWHjvYNdQlMHQ0NLhdR/leZWCWgrK6GHEBJoVZaRjjKGhpYN0t1JITbCSEmeltE4tBWX0EJVKQd1HSiRosTuwO5yku91HIkLuuES1FJRRRVQqBQ00K5Gg3j1HwRNTAFewWZWCMpqISqWgloISCTwdUj3uI3AphdL6Fq1VUEYNUakU1FJQIoFHKWT4KIW8jERa7A5vozxFiXaiUikoSiTwuI/Se7iPQNNSldGDKgVFCZCGVo/7yFcpaFqqMrqISqWgMQUlEjS4B+ykJ3S5j3LSXZZCRYMqBWV0EJVKQWMKSiSob+kgKTaGWGvXxyY1weqeq6BKQRkdRKVSUJRI0NBi75Z5BK5ahZyMBLUUlFGDKgVFCZCG1g4ykmy9jmenJ1DRqEpBGR1EpVLQmIISCepb7N3SUT3kpCdQru4jZZQQlUpBYwpKJGho6SAtwb+lUN/SQYu9MwJSKUpoiUqloCiRoD9LATQDSRkdqFJQlABwOg2NrR3d+h55yHEXsJU3tA23WIoSclQpKEoAnG3rwBhI82MpZKuloIwiVCkoSgDUe/se9bYUslLiiLGIBpuVUUFUKgXNPlKGmwZv2+zeloI1xsKk1Hi1FJRRQVQqBc0+UoabrrbZvS0FgOz0eMpUKSijgKhUCooy3HR1SO1tKYC7gE2VgjIKUKWgKAHQX0wBXGmppxvbcDh12I4S3ahSUJQAaGyxIwKp8X25jxLodBqqz2laqhLdqFJQlACod1czWyzi93lPrYK6kJRoR5WCogRAX9XMHjxVzdpCW4l2VCkoSgA0tHT0mXkEvgVs6j5SopsRoxREpEhE1onI/4pIUaTlURRfTtQ0e60BfyTHWUlLsKn7SIl6wqoURORpEakWkb09jl8rIodE5KiIPOA+bIAmIB4oC6dcihIMdc12yhtaOT+n/7qY7PQEylUpKFFOuC2FZ4FrfQ+ISAzwBHAdMA+4XUTmAeuMMdcB3wP+LcxyKUrA7C13Vc4PpBRy0rWqWYl+rOFc3BjzoYgU9Di8FDhqjDkOICIvAzcbY/a7n68H4vpaU0RWAisBsrKyKC4uDrXYYaepqSkq5R4K0bzn14+5CtfqT+yhuMx/9hGAaW6npKaz2z6jed+DZSzuGUbPvsOqFPogByj1eVwGXCQinwauAdKBx/u62BizSkQqgRtTUlKWFBUVhVHU8FBcXEw0yj0UonnPvyvbRv74s3zqkyv6Pe+QHOO9UwdZfPEybz1DNO97sIzFPcPo2feICTQbY141xvydMebzxpjiAc7V3kfKsLGnvJH52QO/17RWQRkNREIplAN5Po9z3ccCRrukKsNFfbOdsvpW5g8QT4CutFRtoa1EM5FQCluAmSIyVURigduA1cEsoJaCMlzsrQgsyAw6llMZHYQ7JfUlYAMwW0TKRORuY0wncB/wFnAA+L0xZl+Q66qloAwLe8vPAjA/J3XAcyckxxEbY6FULQUligl39tHtfRxfA6wZwrqvAa8VFhbeO9g1FCUQ9pY3kjcuoc+W2b5YLMKcySnsLmsIv2CKEiZGTKBZUUYie8obA3IdeVg8JYNdpY10OJxhlEpRwkdUKgV1HynDQWNLB6fqWjgvgMwjD0vyM2jtcHCw8lwYJVOU8BGVSkEDzcpwEEyQ2cOS/AwAtpbUhUUmRQk3UakU1FJQhoM9Aba38CU7PYHJafFsK6kPl1iKElaiUimopaAMB3vLG8lJTyAjaeAgsy9L8jPYrkpBiVKiUikoynCwN8ggs4cl+RlUNLZpvYISlahSUBQ/1DS1c7K2hfNzB6cUALafUmtBiT6iUiloTEEJN2sPVgNwxawJQV87d3IqCbYYjSsoUUm/xWsiEkj7iTpjzF2hEScwtHhNCTfvHagmKzWO87IHrmTuiS3GwoK8NLaX1HPF/DAIpyhhZKCK5rnAPf08L7gG5ijKqKG908G6I2e4aWEOIn3PT+iPJfkZ/PqD47TP7XuEp6KMRAZSCg8aYz7o7wQR0Slpyqhi4/E6mu0Orpo7cdBrLMnPoNNpONGolc1KdNFvTMEY8/uBFgjknFCjMQUlnLx/oIp4m4VlMzIHvcaiPFew+Wi9I1RiKcqwEFCgWUQmiMgjIrJGRN73/IRbuL7QOgUlXBhjePdANZfNyCTeFjPodTKSYpk+IYkjDWopKNFFoNlHL+Bqcz0V+DfgJK65CIoyqjhUdY7yhlaunJs15LWWzcjkQK2DxtaOEEimKMNDoEphvDHmKaDDGPOBMeYrwCfCKJeiRIT3DrhSUa+cM/h4gofPFeZhd8JfdgY1WFBRIkqgSsFzq1MpIp8SkUXAuDDJpCgR490DVVyQm8bE1PghrzU/J438VAsvbjqFMSYE0ilK+AlUKfyHiKQB/wR8G3gS+MewSaUoEaCmqZ2dpQ1cOWforiMPRblWDp4+x64yTYpQooOAlIIx5nVjTKMxZq8xZoUxZokxJqi5yqFEs4+UcPDG7kqMgavmDd115OHibCsJthhe2nQqZGsqSjgZdJsLEbkhlIIEg2YfKaHGGMNzG06yIC89qKE6A5FgFW5akM1ruys416YBZ2XkM5TeRxeGTApFiTDrj9Zy/Ewzd16SH/K1b79oCi12B6t3VYR8bUUJNYNWCsaYfw2lIIoSSZ7bcJLxSbFcf/7kkK+9IDeNOZNSeHlzacjXVpRQM1CbCwBE5G/9HTfGPB9acRRl+Cmta+G9A1V8rWj6kArW+kJE+MJFU3joL/vYfqqexVMyQv4aihIqArUULvT5uRz4IXBTmGRSlGHlBXcQ+I6LQu868vCZxbmkJdhY9cHxsL2GooSCgCwFY8z9vo9FJB14ORwCKcpw0tbh4OUtp7h63iSy08PX0TQpzsoXL57C/xQf40RNM1Mzk8L2WooyFAYbU2jG1fJCUaKa1bsqaGjp4M5LC8L+WndeWoDNYuHJdWotKCOXQBvivSYiq90/rwOHgD+FWhgRSRKRrZFMd1XGFq9uL2PahCQunhb+Av2JKfF8enEOf9hWRk1Te9hfT1EGQ6CWwiPAz90/PwGWG2MeGOgiEXlaRKpFZG+P49eKyCEROSoivut8Dxj2VtzK2KT6bBubTtRx4wXZgx6mEyz3XD4Ne6eT5z8+OSyvpyjBEmhF8wc+P+uNMWUBrv8scK3vARGJwTWt7TpgHnC7iMwTkU8C+4HqgKVXlCGwZo+rgvnGBaFPQ+2LGROT+eS8LJ7fWEKLvXPYXldRAkX6a9QlIq8bY/p15Qx0jogUAK8bY+a7H18C/NAYc4378ffdpyYDSbgURStwqzGmVzN6EVkJrATIyspa8vLL0RfvbmpqIjk5OdJiDCsjcc8Pb2yltdPwH5clhu01/O37SL2Dhze18cW5sVyVbwvba0eKkfi3Hg6iad8rVqzYZowp9PfcQNlHl4lIfz2OBNeXeDDkAL5VPGXARcaY+wBE5C6gxp9CADDGrAJWARQWFpqioqIgXz7yFBcXE41yD4WRtufKxlaO/PV9/umTsygqmhm21/G37yuM4Y2Kj/mo2s6/f+kKLJbhcV0NFyPtbz1cjJZ9D6QUbg5gDXsoBPFgjHl2oHNE5EbgxhkzZoTypZUxxBu7KwG4YUH2sL+2iHD3ZVO5/6UdvHewmk/OC11XVkUZKv0qBWPMB2F4zXIgz+dxrvuYogwbr++u5Lzs1IjVC1w3fxLZafE89dFxVQrKiGIoDfEGyxZgpohMFZFY4DYgqDbc2iVVGQqldS3sLG3ghguG30rwYI2xcNeyAjYer2NvubaAV0YOYVUKIvISsAGYLSJlInK3MaYTuA94C9fc598bY/YFua7OU1AGzZo9btfRBcOXdeSPz184hcTYGJ7+6ERE5VAUXwJWCiKSICKzg1ncGHO7MWayMcZmjMl1z3nGGLPGGDPLGDPdGPNwsEKrpaAMlrNtHfxuSykLctPIGxe+rKNASEuw8bnCPFbvqqDqbFtEZVEUD4FWNN8I7AT+6n68cICsJEUZcbR1OLjnua2cqmvh29cEdX8TNr68rACHMdr6QhkxBGop/BBYCjQAGGN2EsHeR+o+UoKlw+HkGy9sZ8vJOn7x+YVcPnNCpEUCIH98En+zOJenPjrB5hN1kRZHUQJWCh3GmJ7fwH1XvYUZdR8pwdDW4eA7f9jFewer+dHN87kpAmmo/fGvN51H3rhEvvnyDhpbdGSnElkCVQr7ROQLQIyIzBSRXwIfh1GuflFLQQkEYwzv7K/i6v/6kD/vrOA718zmixeHb2bCYEmOs/LYbYuoPtfO9/+0m/66DChKuAlUKdwPnAe0Ay8CjcA3wyTTgKiloAxETVM7dz2zhXuf30qs1cJv776Ib6wYucWOC/LS+fY1s1mz5zS/26JjO5XIEeiQnRbgQfePoox4Hn7jABuO1/LPn5rrmmMQE4mSnOBYefk01h05w49e38+18yeRnhgbaZGUMUig2UfvuKeteR5niMhbYZNKUYbAgcqz/HlnOV9eVsA9l0+LCoUAYLEI/3LDPJrtDp7V1tpKhAj005JpjGnwPDDG1AMTwyJRAGhMQemP//fXg6TEWfn6FSPXXdQXcyal8sl5WTyz/iRN7dpaWxl+AlUKThGZ4nkgIvlo9pEyAtl0vJa1h87wtaIZpCVGZ1vq+1bMoLG1g99uLIm0KMoYJFCl8CDwkYj8n4j8FvgQ+P4A1yjKsGKM4ad/PUhWahx3DcPM5XCxIC+d5bMm8OS647TaHZEWRxljBDp57a/AYuB3wMvAEmOMxhSUEcVf955mx6kGvnnVLBJiYyItzpC4b8UMaprsvLzlVKRFUcYYwUTg4oA64CwwT0SWh0ekgdGYguKLvdPJf71zmPtf2sGsrGQ+uyQ30iINmaVTx7F06jhWfXic9k61FpThI9Dso/8E1uNyI33H/fPtMMrVLxpTUDzsLW/kpsc/4tH3jnDDBZN5eeUlWKMk22gg7v/EDCob2/jjtkBHoivK0AmoTgG4BZhtjGkPoyyKEhRtHQ7ueHITcVYLv/nbwlE3rOayGZksnpLOE+8f5bNL8oi1jg5lp4xsAn2XHQeiM5VDGbV8cPgMja0d/PxzC0adQgDX2M5vXjWLisY2/rBNq5yV4SFQS6EF2Cki7+FqdQGAMebvwyJVFHO0+hxHqpq47vzIDnAZC6zZU0lGoo1Lpo2PtChh4/KZmSxSa0EZRgJ9h60GfoSrCd42nx+lB0+sPcZ3X9kdaTFGPW0dDt47UM01500aNTEEf6i1oAw3gfY+ei7cggSDe+jPjTNmjLyK1WNnmjjX1onDaYixSKTFGbV8ePgMTe2dXD8GLLLlai0ow0ig2UczReSPIrJfRI57fsItXF+M1OwjYwzHzzQDaIuCIfKTNw/w/Vf7biO9Zk8l6Yk2Lpk+el1HHnythc/+78f897uH2VZSR6fDGZL1Ox1Onlx3nNomzSNRAncfPQP8CugEVgDPA78Nl1DRSvW5dq8yONuqw1KGwuqdFby0uZQn1/Ueat/W4eDdA9VcM29S1DS7GyrLZ2by4PVzMcCj7x3hM7/awA2//IhTtS1DXvuV7WX8xxsHeHGTFsopgSuFBGPMe4AYY0qMMT8EPhU+saKTY2eavL+fbVOlMFjOtnVQ2dhGcpyVn/71IJuO13Z7ft2RGpfr6ILR7zryICLcu3waq++7jB3/8kl+8bkFVDa2cePjH7HuyJlBr9ve6eDRd48A8PGx2gHOVsYCgWYftYuIBTgiIvcB5UBy+MSKTo65XUcA59rUfTRYjlS5lOuPbjmPX753lPte2sEb91/GxNR4oMt1dOkYcB35Iz0xlk8vzqUwfxz3Pr+VO5/ezN2XTWV8chztHU5E4O7LppIUN/DH+8VNp6hobGNBbhrbTtXT1uEg3hbdLUKUoRGoUvgHIBH4e1xZSJ8A7gyXUOFm1YfH+J/iY+z4l08iErpg8HFfS0HdR4PmSNU5ABZPyeBXX1zCLU+s54tPbeKSaeOZkBLHu/uruO78seM66osp4xN59euX8t0/7uY3PdxsSXFW7r5sar/Xt9g7eWLtUS6ZNp6Vy6fx5We3sPVkPZfNzAxalmfXn0BEuDOKGxEqLgLNPtri/rVJRO4Gko0xZ8MnVvhpaOmg2e4gOYC7qUA5dqaZ1HgrZ9s61VIYAoermoi3WcjLSMRiER69bSE/e+sQf9pRzln3/+uti6K/v1EoSIqz8sQdi/lJWwc2i4VYq4VP/896XtlWNqBSeGb9SWqa7Pz6S7OZMykFq0VYf6xmUErhN+tO0NTeyRcumjLwycqIJqBvRBF5Efgq4AC2AKki8qgx5mfhFC5cpMa7irPPtnaEVCkcP9PEgrx01h2p0ZjCEDhSfY6ZE1OwuFN6rz5vElefNwlwBZlb7Q4yknRUpS+e9zTAZ5bk8tBf9rG/4izzslP9nt/Y0sGvPzjGlXMmsiQ/A4CFeemDiivUN9spb2gFYPOJukFIr4wkArW/57ktg1uAN4GpwJdCKYiIzBWR/3Wnvn4tlGv3JMX9AQrl3Xxbh4PyhlYW5qWHfO1ooL3TwYHK0BiPh6vOMTPLf8gq3hajCmEAbrwgG1uM8Mp2/4303jtQxQ2Pr6OpvZNvXT3Le/zSGZnsKWugMUjX536fv/vb+04PTmhlxBCoUrCJiA2XUlhtjOkggMlrIvK0iFSLyN4ex68VkUMiclREHgAwxhwwxnwV+BywLKhdBElqgss6COXd/ImaZoyBWVkpJMbGjKmYgtNpuO/FHVz36Dq2ngz8TrG2qZ0frznQbZBMY2sHVWfbmZWVEg5RxwQZSbFcNTeLP+8op8OnlqGioZV7ntvK3c9tJc4awwv3XMx52V21PpdOH4/T0CvbayD2lrta2F9YkMHb+6v6rC0Z7ZxodDD/X9+iwm01RSuBKoVfAyeBJOBD9zjOQG4LnwWu9T0gIjHAE8B1wDzgdhGZ537uJuANYE2Acg2KLkshdF/cnqK16ROSSYm3jin30X++dZB39lcRG2Ph8bVHA77utxtPserD47xzoMp7zBNkntWHpaAExmcW51LbbOeDQ6501YqGVv7mVx+z/mgND1w3hzV/f3mvwr9FU9KJt1mCdiHtqzhLdlo8n79wCpWNbZw4G5qiumjjaL2TpvZODp4OzmIuqW0e+KRhJNBA82PAYz6HSkRkRQDXfSgiBT0OLwWOGmOOA4jIy8DNwH5jzGpgtYi8Abzob00RWQmsBMjKyqK4uDiQLXTjdLPrTbtx+x4spw8Efb0/3j1qB6D0wDZinHaOl1ZSXFzv99ympqZByT0SWVfWwVN77Xwiz0pGvPDKoTM8+5f3KEjrSmt0GkNTU3O3PRtjeGmD647q5Q/2kFp/GIDiUpcyrT2+j+IQ/W0iScT+1k5DSiz86q0dtJXF8eNNrdS3Gx5YGk++KeXjj/z3UZqRKry9u4Si1MBrH7YcbWFSkoX42iNYBDaWtjJtlLy/g6GssR0QijfvxnI6sKbSJxsd/HBDG19fGMfSSaGLbw6FQAPN/4Crqvkc8CSwCHgAeHsQr5kD+L4jy4CLRKQI+DSuCW99WgrGmFXAKoDCwkJTVFQUtABnzrXDunfJmzqDoksKgr7eH38+vYOc9HquuXIFvz60nrjYGIqKLvZ7bnFxMYORe6Sx+UQdz7+zkctnZvLruy6krcPBOz99nw1n07jr5kIAKhtbuX3VRqYk2Hj+viLvtXvLGzn91kekxFs52Gjh8uVXEGMRilfvIzG2lM9cu8IbaI5mIvm3/lzLfp7fcJJfH7JR097G81+5iIsH6Ch7UI7x0zcPMm/JxUxMiR/wNVrsnZx+6y0+f8k0brhqFi+WbGTP6bpR8f4Oll9s/SvgIGliHkVFcwK65vkNJ4F9bKxL4Lu3BeY1t3c6eWNPBRcWjCM3I3HQ8vZFoO6jr7gDzVcDGbiCzD8NpSDGmGJjzN8bY/7OGPNEf+cOdRxnSrwnphC6YPCxM81Mm5AEQGqCbUwEmv/rncNMTInn8S8sxhZjISXexl3LpvLWvioOnT5HXbOdLz21mZO1Lawr6+xWx/HargqsFuE718ymrtnOrrIGwJN5lDwqFEKk+cziXDochv2VZ3niC4sHVAgAy6a70lE3BOhCOlB5FmPwxiauOW8Slc2Go9VNA1wZOI2tHbTYR/7nqarF5YGoDCKmsL/C5WrafqqBPWWBfZ+VN7Tyj7/bxcbj4cn0ClQpeD6h1wP/Z4zZ53MsWMqBPJ/Hue5jw0a8LYZYqyVkfn9XI7wmpk9w+cFT421RF2g2xuBwBh4gbOtwsK2knuvPn0RaQpep/OVLC0iKjeGRtw9x1zObKa1r4YkvLMZqgV8VHwNcgenVuyq4YtYEbl6Qg0Vg7cFqwFWjMFODzCFhXnYqXy+azuO3Lwp4CNG87FTSEmw8v6GkWxaSMYbnN5zkyp8XU1rX1W9pn/tL7Tx36qvndd7ef7rbtUPh87/ewCU/eZ8n1h6leYQ2mux0OKlpde2zoqEt4Ov2VZxlQW4aibExPPvxyYCu8fz/52UkBC1nIASqFLaJyNu4lMJbIpICDDaatAWYKSJTRSQWuA3XvIaACUWXVNcXd2jeYFVn22m2O5juthRS3AVs/fHx0ZqgszzCycr/28Zdz2zu9QHucDh5ZVtZr+HxW0/WY3c4uXR690KnjKRYvnhJPu/sr2JfxVn+547FfOqCyVyRa+VPO8opq29ha0k9lY1t3LQwm7REG0vyM3j/YDX1zXbOnGvXIHMI+e61c4Ia+BRjEf7tpvPYXdbArf+znhM1zdg7nfzgT3t46C/7OHammd9uKvGev6/8LBmJNianuVxN2ekJTE218NLmU3zjhe2seKSY2f/yV7aV+I+vDURNUzsHT58jOc7Kz946xOX/by2vjMCZ1RUNbTgMWC1CRWNglkKHw8mh0+e4aNp4Pr04h9d2VwTUqfaUWylMGR961xEErhTuxhVDuNAY0wLEAl8e6CIReQnYAMwWkTIRudsY0wncB7wFHAB+77Y8Amao7iOA1HhryLKPPG6RaR5LIcHGubaOfu+Qvvfqbh55+1BIXn+onGvrYO3BatYdqeHNvd3zzJ9Ye5R/+sOuXsPj1x+rwWoRlk4d12u9ey6bxsK8dH7xuQVcOdd153jdVBsisOrD4/xlZzkJthiucj+3Ys5E9lWc5aOjNQBqKUSYWxbl8Nu7L6K+2c4tT6zns7/ewEubS/la0XQ+OS+LP24tw97puifcW9HI/Jy0bu1iLs2xUlrXyp7yRuZMSiEuxsJvN5b09XL9suNUAwD/fdtC/vT1S8lOj+ef/7x3xKW9nnRnEC2aks7pxraArO5jZ5qwO5ycl53KnZcUYO908vKWgQcplda3EBtjISuAmM9g6FcpiIgnWrLQ/e80EVkM5BNAkNoYc7sxZrIxxmaMyTXGPOU+vsYYM8sYM90Y83CwQofCUkhJsIUspuDpjupxH6XEW+lwGNo6/BtTZ1qclNa1UtdsD8nrD5X1R2vpdBrSEmz8eM0B2jpcVsGRqnM84U4xfX1XZbdrPj5Wy8K8dL9N1yakxPHnbyzj5oU53mPjEyx8ZnEuL28p5fXdlVw1L8t77SfmTATgN+tcIzq0RiHyXDRtPH/5xmVMTInjYOVZHr1tId+7dg53XDSF2mY7b+8/jb3TyeGqc72qpj+Zb+PQf1zLh99dwa++uISbFmazZk/loNy1O07VY7UI5+eksWhKBp9ZnEtrh2PEfHY8eNJKL5k2nk6ncSWzDMC+8i7X28ysFJbNGM8LG0sGnJNRWtdCbkZC2OJuA1kK33L/+3M/P4+ERaJhIpSWwrEzzSTFxpCVGudeu/86iP21ri/dhpaREXf44HA1KXFWHrt9EWX1rTz10QmcTsP3XtlNcpyVL148hU0naqk+5/KVNrZ2sKesIegupV+9YjqdDieNrR3cvCDbe3x2VgrZafHsLmskOc5Kdlp47oCU4JgyPpHX7r+Mdd9b4VXwy2dOICc9gZc2n+JI9Tk6HKZbAZyHOGtXSvJnC/No73T2urEIhO2n6pmXnert3JqT7vKjl9WPrAKxk7UtxMbAAndHg0BcSPsrzxJvszA103UzeeclBVQ0tvHO/qp+ryutayV3XHhcRzCAUjDGrHT/u8LPzyfCJtUAhMZ9FLpg8PGaZqZNSPaa0KnuwGtfd0ZepdDav4tpODDGUHzoDMtmZHLFrAlcPS+LJ9Ye5RfvHGb7qQb+5YZ53HlJAU4Db+5xuZY2Ha/FaVxtEYKhIDOJWxblkJkcx/JZE7zHRYQVbmthxsTkkHauVYZGvC2mW2qqxSLcvjSP9UdrWbPH9SU/v4/+Sh4W5KYxKys56BnTnQ4nu0obWTwlw3ssxx1cLR9hVcMltc1kJVq88gVS1byvopE5k1K9Y3uvnJtFbkYCT6/vPVjKl1N1LUwZF54gMwQ+jvNv/f2ETaoBCIn7KN4asrTRktpmCjKTuq0N/lNejTEcqHMgAg6n4VyEsykOVzVR2dhG0WzXl/QPrp9Lh8PJ42uPcvnMTG5dlMPMrBRmZ6Xw+u4KwOU6irdZWDQlPejX+/Gt5/PWNy/vNWfY40LSIPPI57OFecRYhN98eIKk2BgKxif1e76I8Nkleew41cDR6nMBv86hqnO0dji6vc9y0113yOUDWAp7yxu96Z7DwcnaFiYmCpPTAlMKxpheDQtjLMKXl01ly8l6dpU2+L2usbWDxtYO8sJQn+Ah0EDzhT4/lwM/BG4Kk0wDEhJLIcEWspTU2iY7E5Ljutb26cLak0NV5zhrh6UFrgBtQ3NkXUjFh1ypoFe4lUJBZhJfu2K6K75w6/neu/YbLpjMlpP1VDa28vGxGi4sGNfNRRAo8bYYxvv8X3m4dHomOekJLAvS+lCGn6zUeK6cMxG7w8ncyakB+bZvWZSD1SL8YWvgmUPb3UFmX0shNcFKcpy1X0vBGMPXX9jO91/dHfBrDQWH03CqtoWsRAup8S75BkpLLatv5WxbpzeV18PnCnNJibPy1Ef+rQVPOuqUSLmPPBhj7vf5uRdYTAQnr4XEUoiz0tbh9GZRDJa2DgdN7Z2MT+7q3Jnaj6Ww/qgrDfVT7lGS9S2DC5jtLW/krmc20zjEuETxoTPMmZTivcMB+NbVs9n0gyvJ83nj3eCOATz78UkOVzX1SkUdKgmxMax/4BPdgtPKyOV299yE+TmBfQYnpMSxYs5EXtnevUlff+woqSczOY5cn3x8ESEnPaHfmMKhqnOcqmvh4OlzQdXeDJbKxlbsDicTEwURITs9fkBLwVPfMW9yd6WQEm/jtqV5vLGn0u8a3hqFSCsFPzTjap8dtXj8/kMNNte6syDG+7Rz7m/tDcdqyEoU7x3CYJXC1pN1FB86w3+/d3hQ1wM0tXeytaTOayX40nMk49TMJObnpPK0+w5mrI7CVFwsnzmBuy4t4G+WBD7s6HOFedQ0tXsLFQdiR2kDi6ak94ox5WQk9GspvL3PFaht73Ryoib8zeZKal1f1FmJrq/TyWkJVDb2bynsrzyLRWDOpN7xGM/0uuf8FLOV1o8QpSAir4nIavfP68Ah4E9hk2oY8Pj9hxpXqGtyKwW/7qPua3c6nGw6Xsfc8TGkJ7qUyGAzkJrd7aaf31Di7SwaLOuP1tDhMBTNmhjQ+TdckE2Hw5ASbw34DlEZncRYhB/edF5Q74Oi2a7Mpe+/usfbbtvDm3squf+lHd4q6rpmOydqmru5jjzkpCdQXt/S67iHt/efZpz7Ji3YjqUeGlrsvLz5lN9EkLf2ne5WjOepUZiY6FJe2ekJA1oK+ysamT4hmYTY3i7Y3IxErps/iRc3n6KpR8zxVF0LaQm2bl0EQk2glsIjdKWi/gRYbox5IGxSDUCoso9g6DMVappd+ci+7qN4mwWrRXpZCrvLGznX3sm8cTFkuJXCQJbCnrJGv1WOze2dxFiEpNgY/v31/YPKYio+dIbkOCuFBb0/eP74lLsy9uJp470ZE4oSKLYYC/9391LirBZu/81Gtp6so63DwYN/2sPXXtjOa7sq+O4fd2GMYccp15fuYj/JDDkZCZxt6/T72S1vaGVv+VnuurSAGIsMevDTH7aW8cCre7zVw77861/28Z0/7MLpdk2V1LYQa7WQEe9WCmnx1DbbvfU+/uhvKh7APZdP41xbJ3/Y2j1jq7SulbwwZh7BwMVr8SLyTeCzwBxgvTFmvTEmonXmoco+gqFbCrVNvd1HIuI3kP2xu2J3zvgY0hJcFb71/VgKDqfhtlUbWPXh8V7PtbjnS3/zqlmsO1LDewf6N8kbWuz81zuHueoXH3DT4x9x59ObWbOnkmUzxmOLCezeIG9cIg/dMI9vrJgR0PmK0pNpE5L5w9cuJTM5ji89tZkbfvkRL2w6xd8tn8YD183hrX1VPLP+JNtP1RNjEc7P7f0Z99Qq+MtAesc9+e3GBdlMn5DEwcrBWdFH3FlSPd1UHQ4nVefaOF7T7K3AP1nTTP64RCzSZSkAfbqQ6pvtVDS29Qoy+7IwL50l+Rm8sOlUt+OldS1hDTLDwJbCc0AhsAfXUJyfh1WaYcRbSzDEWoU6r6XQPaMmJd7ay320/mgtcyenkhorxFiE1HgbDf1YCuX1rTTb/VdvNrV3khQbw5cuyWf6hCR+9Mb+Xs3CPI36fvLmAZb99H0efe8Ik1LjGZcUS32LnYxEG5+/MK/X2v3xlcumekeOKspgyElP4Pd/dwkFmUnUNdt55ssX8v3r5/J3y6dx1dwsfvLmAV7bVcncySkkxvaumPcEnv0phbf3VzFzYjJTM5OYMyl10JaCp8vr6R5f7NXn2vEY5a621y5LwTcl3aMU+nIhdQWZ+7+pvf78yRytbqLM7SpzOg1l9a1hTUeFgVtVzDPGnA8gIk8Bm8MqzTASSkshzmohqYdvMDXe1s19ZO90su1UPV+6OB9w3dVnJNr6tRQ87TOa/bQNbrF3khhnxRZj4aEbz+POpzez4N/e5rzsVBbnZ9DY2sGGY7VUNrZhEVc84BsrZjB7kraQUCLPhJQ4Vt+3jE6H8frVRYRHPnsBn3rsI07VtfC3l+T7vbavAraGFjubTtTx1SumATB3ciqrd1XQ2NJBWmLgPnhjulp/97zb97TFXpCbxnsHqympbaakrpnlszIB1zXZ6a5iv76Uwv5Kl9u7P/cRwBWzMvkRsO5IDbcvnULVuTbsDmdYg8wwsKXg/cZyN7IbEYSqTgFCEFNosjM+KbZXhkRqQvdOqeUNrdg7XXndHtITY/u1FDxKoam9t2+yud3h7R10xawJvLzyYlYun0acLYYXN52i+NAZFk/J4D9umc8H31nBY7cvUoWgjChsMZZegdb0xFge/8IiEmwxLJ/ZOysOIDMpjlirpZdSeP9gNQ6n4ep5kwCYM9n1fg822FzTZPd+dnt+sVe4lcQ/XT0biwg/f/swbR1O8n0K+CaleZSCf/fRodNNTEyJ8wbD+2L6hGQmp8Xz4WHXFLzSOpcs4VYKA1kKC0TE8z8qQIL7sQDGGNO/qgsTxpjXgNcKCwvvHewaybFWRIY+aKeuud1vMVZKnI0z57oGjXgaZuWPT6TZ7eZMT7R5YxL+8FoKfqqem93uIw8XTxvvHaLicBoEdFCNEpUsmpLB7h9e3Wesy2IRdwZS9y/st/dVMSk1nvPdGVGeGoADlWe5KIABQx58BwT1ZSksmpLONedlsXqXq8q/YHwSne6pMHHWGCakxPVpKRypPhdQ00cRYfnMCby5t5JOh7OrZXYkLQVjTIwxJtX9k2KMsfr8HhGFECosFiE5zjrkmEJts71b5pGH1ITuMQV/f9CMxNh+s4+OVbsUiV+lYHf49beCK11QFYISzQyU/JCTnkCZz5duW4eDDw6f4ap5E73v/YkpcWQk2jh4Orhg81H3zdh52am9vtgrG9tIibOSEm/jTp9Rvvk9Zhtkp8X7bYrndBqOVDUxM8B2LstnTeBsWye7yhoprWtBpMs9FS4GW7w2KnD5/YceU/BnBqb0iCmcqm0h3mZhYkqXVZGeaOu3TqHLfeQ/ppAcF3ybCUUZDfS0FDadqKO1w+Gd3wGuO+25k4MPNh+rbiIpNoZFU9J7WQoVDa1Mdn8pL506jjmTUrDFiDe47KGvWoXyhlZaOxwBt4dfNmM8FoEPD5+htK6Fyanxg2ovEwxjWim4JqQN3lIwxlDb3E6mH/dRaryNZrvD2xu9xJ1K5ht7yEiMpam902+rjYYWO7XNdizSt/so0c8sA0UZC+RkJFDT1O6tBVh7sJo4q4VLeriJ5kxK5VBVcO0ujp1pYvrEZLLTE3rNh65sbPO2hBER/v3m+fzg+rm96nay011VzT3rhw67C00DbfyYnhjLBbnprDtyhtL6lrC2zPYwppVCzwyhYGmxO2jrcPZhKXTPbjpV2zu/OMOdEdHQ2tuFdOyMy3U0Y2IyzX0Fmv1UQyrKWMBbq+C+G//g8BkumT6+V3uWuZNTaOtwequOA+FotWveera342mXtVDZ2OodPQoua+HLy3p3/JmcFk+L3dFtzjW4uhIDzJgYeNLH8lkT2FnawKHT58IeT4AoVQqhyD6C3n7/YPFXuNa1tqf/USfGGHcP9O4thvtrdeFxHV2Qm47d0b1xn8NpaO1w+J16pihjgRyfWoWTNc2cqGlmxeze7Vo82X6BFrE1tXdS2djGjInJ3iwiT61Ce6eDmiZ7t+aRfcqX7j9t9kjVOSalxgfVpmL5zEycxpUUE+4aBYhSpRCKimZw+/3bA7cUWnrUC9S6C9f8u488nVI7ONPUTmuHo1cwytvqwk9x2rEzTcTGWJjjTiP1dSF55EjqI9CsKKMd3y9dT/v3Ij+NHWdMTCbGIgGnpR73Ga3rtRTcAeOqRtfnfXIAgV5vVXOPtNTD1ecCDjJ7WJiX7vU8TBkf3hYXEKVKIVSk+qk67oufvXWQRf/+TrcvcI+l0FegGVxK4VSt/1SydLf7yF8B27HqZgoyE70WR1M3peByJyVqoFkZo0xOiyfGIpTXt7L20BmmZSZ1qxXwEG+LYVpmUsDBZk866oyJSWSluW72PF/sHuWQHYCl4FEcZT6N+5xOV1FcsDPIrTEWlrlb1aulEGY8GUIDNZN7YVMJT6w9Rnunk0M+HUlr/TTD85Ca4LYUWju9rXWn9LQU3Mqk0U9M4fgZl18z2e0i8q1q9lgNyeo+UsYo1hgLk1LjOXamiY3Ha/22f/cwd3Iqe8obu31B98XR6iasFiF/fBJx1hgyk+OodCsDz7+BWAoTkuPISo1j4/E677HS+hbaOpyDmi547fxJxNssTJ8Q/jE2Y1oppCZYcZquNtT+WHuomof+ss/b7+f4ma6AVdcsBf/ZR+C2FNz5xb7DQqAr0NzTUrB3Oimpa2H6hGRv3MDXfeQJPPdVp6AoY4Gc9ATeO1hNe6fTbzzBw0XTxlF1tp3L/nMtV/68mJ+8eaDP4VpHq5vIH5/orZPITo/3VjF7As6BWAoiwjXnTaL4cDWt7u8XT5B5ZpCWAsDNC7PZ9IOrvDeS4WRMKwWPi6evDKR9FY1844XtzJmU4m356/E5gst9lBgb47cnemp8V6D5VB/5xQm2GGKtll4FbKfqmnE4DdMnJnlrEXxbXTR7YwrqPlLGLjkZCdg7nSTYYlg6dVyf591xUT7v/ONy/vlTc5mUFs+vPzjOa+5K5J4cc1voHianxXurmCsbW0lPtPn9vPvjmvMm0dbh5MMjrjYVnnTUmRODv9sXkbDOUPBlTCuFvobhgEtRfPW320hLsPH0XReSEm9jamYSx30mOdU2tft1HQEkewLNrR2U1Db3ch2B6w+dnmDrNaf5qLuSuS9LwRtoVveRMobxBJsv9ZOK2pOZWSncc/k0/u8rF5ESb2Wrz5AcDx0OJyW1LcyY6KsUuqaoVTa0MSk18GripVPHkZ5o4629rnbeR6rOkZ0W770ZHamMaaXQVUvQ21L44er9lNe38svbF5HlfiNMm5DU3VJotvt1HYGr1URynNVtKbSSP653EAz8t7rwpKNOm5DszTBq8uM+StJAszKG8aSl+ss66guLRVg8JYNtJXW9niupbabTabophez0eJraXQN9KhrbelUu94ctxsKVc7J490AVHQ4nh6uaBuU6Gm5GlFIQkVtE5Dci8jsRuTrcr9dXp9Q3dlfyyvYyvrFiBoUFXWbptMxkSutbvf7IWneH1D7Xj7dSdbaNmqZ2v5YC+G91cexME5NS40mOs3YFmtt7B5o1pqCMZS6eNp6FeelcM39SUNcV5mdwuKqJxpaeFnpXOqoHT01CZUMbp3sUrgXCtfMncbatk/VHazh2pmlQQebhJuxKQUSeFpFqEdnb4/i1InJIRI6KyAMAxpg/G2PuBb4KfD7csvmbqVDZ2MoP/rSHBXnp/P2VM7udP31iEg6n4VSdy71T29y3+8i1vo29Fa4Cu74qEf1bCs1Mn+iyLPwGmt2BK61TUMYyUzOT+PM3ljExJbgv6iXu8bPbT3V3IXm6CEzvYSkAnKhpor6lIyhLAeDymZkkxsbw5LoTtHc61VJw8yxwre8BEYkBnsA1zW0ecLuIzPM55Z/dz4eVFB+/v4fvvbKHDoeT//78wl6dGqdlut4sR6ubMcZQ12z32zbbQ2qC1ZuO2rNwzUNGUvdBO8YYjld3BbtirRZiYyzdAs0tHktB3UeKEjQL89KJsQhbe7iQDp4+x+S0+G6p3pPclsKOUw0AQVsK8bYYimZP8I7uDLZGIRKE/VbTGPOhiBT0OLwUOGqMOQ4gIi8DN4vIAeCnwJvGmO3+1hORlcBKgKysLIqLiwctm93hqk/Ytf8wee0nOWs3fHi4hVtn2CjZu4WSHue3drrOf3fzbhyVNjochvrKUoqLT/tfv7mrmrH0wA7qjrqaZjU1NXnlPnvGTn1zB2vXrkVEaGhzcq69E2dDJcXFrjdSrMXJ4RMl3tc5cNSO1QLr13046L0PN757HkuMxX1Hw56nJAvv7jzBhXGuz1Rbp+HtvS0snWTtJrtnNsn7u08CUHXiEMVnj/pds69950mXlX/60A6Kj43stvaR8j/kAKU+j8uAi4D7gauANBGZYYz5354XGmNWAasACgsLTVFR0ZAEiX3/TTKzp1BUNIe3950GtnHHVYXdYgm+TNj0LqRMYO7C6fDeB1y0cC5Fi3L9nvvn0zvYdaaC1Hgrn/rkCu/x4uJiPHIfthzjjRMHKbzkMlLibaw9VA3FW7h+2SIuneGqYkzf9D7p48dRVLQQgPca9pJSVcFQ9z6c+O55LDEW9x0Ne/7g3D5e2nyKZZcvxxZj4Q9bS2l37Obvb7yQC3t89rM2vkeJu3vB9UUX+62chr73vbitg6f3vcPElHiuvWpF7wtHGCMq0GyMecwYs8QY81V/CsFDqBrigWdspst9s+VkHbFWC+fn9t1TaVqmKwOprp/Cta61XYHsvt5E0Lsp3qbjddhihIVT0r3nJMdZu2cf2Ts1HVVRhkBh/jjaOpzsq3C1v/jDtjKmZiZRmJ/R69zJ6fHe5JKsIFJSPaTG2/jM4lw+OS9r4JNHAJH6ZikH8nwe57qPDTspPoN2tpysZ2Fuer9DLKZNSObNvZXU9NP3qGttTxOrvvuVZPgohbxxsPF4LRfkpnfLLEqMjenV5kKDzIoyeArdweatJ+vISLSx+UQd37lmdq9Z6+CqYN5BA+OTYgesh+iLn37mgiHJO5xEylLYAswUkakiEgvcBqwO9OJQdUkFT1O8DlrtDvaWN3rfLH0xfUISDS0d3loCfx1Su9Z2WQr99UDvanVhp7m9kz3ljVw8rbv5mhRn7R5otjs0yKwoQyArNZ68cQlsPVnPH7eVYRH49OIcv+d6gsuB9DwaDQxHSupLwAZgtoiUicjdxphO4D7gLeAA8HtjzL4g1gyZ+8jTFG9naQOdTtPLn9iTaRNcrqDNJ1yZC/1ZCl73UT9KweM+qm+xs7WkHofTcNHU7tOjkuOsveoUtBmeogyNwvxxbC2p55VtZVw2c0KfcxImu9NQA5mjMBoIu1IwxtxujJlsjLEZY3KNMU+5j68xxswyxkw3xjwc5JqhsxQSrJxt62TryTpEYPGU/i0FT1rq9pJ6UuKtxFr7/i/0Wgr9uI887bMbWjrYdLwWq0VY0sOvmdRLKThI1L5HijIkluRnUNPUTkVjG59d4j9ZBCDbbSlkB5mOGq2MqEBzJEiJc1kKW0rqmZ2VQlpi/31JcjMSsMUI59o7+3UdAVw2M5OvXjG915e8L+kJXe6jjcdrOT83rVcQ2W+gWWMKijIkPK7i1Hhrv0Fgr6UQZOFatBKVSiHU2UcNLR1sL6kfMJ4Arj7unmyi/lpcAKQl2Hjgujn9Bq6tMRZS4q1UNLSyu6yRi3sMHgdXj6Pm9k7v3IcWu47iVJShMmtiCpNS4/lsYV6/AeTpE5KYMTF5QNfyaCEqv1mMMa8BrxUWFt471LVS4m20dzpp73QG/EeflpnE0eqmfuMJwZCRGMv7B8/Q6TRc5KcFcFKca+5DW4eThNgYmto7NdCsKEPEYhHe/tZyEgbIKEqJt/Hut64YJqkiT1RaCqHEM0sZ6LNgrSfT3C0o+mtxEQwZiTZqmtqJsYhfGTxB5ab2TjocTuydTnUfKUoISI239WpnM9aJyv+NUGcfgSuIlBOgz9CTgZTZTzO8YPBkIM3PSfObVeRRAM3tnV3zmTXQrChKGIhKpRDa7COXUgjUSoCu1rqhcx+5ZOhZn+AhycdS0PnMiqKEkzH/zeKpOr4wgCCzh/OyU7lyzkQunZ4ZEhk8lsLFU3sHmYFuMxXibS49nqhKQVGUMBCV3ywiciNw44wZM4a81gW5ady+dArXnz854GvibTE8ddeFQ35tD7kZCSTYYvrMfvJMWGu2dxLf7vpd5zMrihIOxrz7KDHWyk8+fX7IgsaD4YsX5/POt5b3Obu1y1JweN1HmpKqKEo40G+WEUC8LYbcjL6rnpO6uY88loL+6RRFCT36zRIF+AaaE9xuI61TUBQlHESl+yiUKanRgCd+4HIfuVJSNftIUZRwEJVKIZQxhWjAGmMh3mah2d6Vkqp1CoqihAO93YwSPE3xPMogUWMKiqKEgai0FMYinvbZLXYHCbYYYiwje/i3oijRiSqFKCEp1qUUmto7vXULiqIooUaVQpTgcR+1tHeq60hRlLARlUphrGUfgWemgoNmnaWgKEoYiUqlMNayj6ArptDc3qktLhRFCRtRqRTGIh73UbPdoc3wFEUJG6oUogRv9pFaCoqihBFVClFCUpyVZrvDnX2kloKiKOFBlUKUkOxOQ61paldLQVGUsKFKIUrwWAcdDqMxBUVRwoYqhSjBtwGeNsNTFCVcjBilICLTROQpEfljpGUZifjOT9BmeIqihIuwKgUReVpEqkVkb4/j14rIIRE5KiIPABhjjhtj7g6nPNGM7/wEHbCjKEq4CLel8Cxwre8BEYkBngCuA+YBt4vIvDDLEfX4uow0+0hRlHAR1m8XY8yHIlLQ4/BS4Kgx5jiAiLwM3AzsD2RNEVkJrATIysqiuLg4ZPIOF01NTUHLXdnk9P5+7NA+iusOhViq8DKYPY8GxuK+x+KeYfTsOxK3nDlAqc/jMuAiERkPPAwsEpHvG2N+4u9iY8wqYBVAYWGhKSoqCrO4oae4uJhg5a462wYfvQfAxYWLWTp1XBgkCx+D2fNoYCzueyzuGUbPvkeMH8IYUwt8NZBzReRG4MYZM2aEV6gRRFI395EGmhVFCQ+RyD4qB/J8Hue6jyn9kGjTQLOiKOEnEkphCzBTRKaKSCxwG7A6mAXGYpdUi0W8lcyJaikoihImwp2S+hKwAZgtImUicrcxphO4D3gLOAD83hizL8h1x9w8BehyIWnxmqIo4SLc2Ue393F8DbBmCOu+BrxWWFh472DXiEaS46ycaWon3qqWgqIo4WHEVDQrA5MUZyXRFoPFIpEWRVGUUUpUKoWx6z6K0WZ4iqKElahUCmMx0Awu95G2zVYUJZzobWcUcdelU6lpao+0GIqijGKiUimMxeI1gMtmZkZaBEVRRjnqPlIURVG8RKVSUBRFUcJDVCqFsZp9pCiKEm6iUimo+0hRFCU8RKVSUBRFUcKDKgVFURTFS1QqBY0pKIqihIeoVAoaU1AURQkPYoyJtAyDRkTOACWRlmMQZAI1kRZimBmLe4axue+xuGeIrn3nG2Mm+HsiqpVCtCIiW40xhZGWYzgZi3uGsbnvsbhnGD37jkr3kaIoihIeVCkoiqIoXlQpRIZVkRYgAozFPcPY3PdY3DOMkn1rTEFRFEXxopaCoiiK4kWVgqIoiuJFlUKYEZGnRaRaRPb2OH6/iBwUkX0i8v8iJV848LdnEVkoIhtFZKeIbBWRpZGUMdSISJ6IrBWR/e6/6T+4j48TkXdE5Ij734xIyxpK+tn3z9zv790i8icRSY+wqCGjrz37PP9PImJEJCqnYmlMIcyIyHKgCXjeGDPffWwF8CDwKWNMu4hMNMZUR1LOUNLHnt8G/ssY86aIXA981xhTFEExQ4qITAYmG2O2i0gKsA24BbgLqDPG/FREHgAyjDHfi5ykoaWffecC7xtjOkXkPwFGy7772rMxZr+I5AFPAnOAJcaYaClm86KWQpgxxnwI1PU4/DXgp8aYdvc5o0YhQJ97NkCq+/c0oGJYhQozxphKY8x29+/ngANADnAz8Jz7tOdwfWGOGvratzHmbWNMp/u0jbiUxKign781wH8B38X1fo9KVClEhlnA5SKySUQ+EJELIy3QMPBN4GciUgo8Anw/suKEDxEpABYBm4AsY0yl+6nTQFak5Ao3Pfbty1eAN4ddoGHAd88icjNQbozZFVmphoYqhchgBcYBFwPfAX4vIhJZkcLO14B/NMbkAf8IPBVhecKCiCQDrwDfNMac9X3OuHy1UXsH2R997VtEHgQ6gRciJVu48N0zrj3+AHgokjKFAlUKkaEMeNW42Aw4cTXTGs3cCbzq/v0PwKgKNAOIiA3Xl8QLxhjPXqvcPmiPL3pUuQqhz30jIncBNwB3mFEWvPSz5+nAVGCXiJzE5S7bLiKTIifl4FClEBn+DKwAEJFZQCzR011xsFQAV7h//wRwJIKyhBy3pfcUcMAY8wufp1bjUoi4//3LcMsWTvrat4hci8u3fpMxpiVS8oUDf3s2xuwxxkw0xhQYYwpw3fgtNsacjqCog0Kzj8KMiLwEFOGyBKqAfwX+D3gaWAjYgW8bY96PkIghp489HwIexeU6awO+bozZFikZQ42IXAasA/bgsvzA5U7YBPwemIKrzfvnjDE9g/BRSz/7fgyIA2rdxzYaY746/BKGnr72bIxZ43POSaAwGrOPVCkoiqIoXtR9pCiKonhRpaAoiqJ4UaWgKIqieFGloCiKonhRpaAoiqJ4UaWgjApExOHuwOr5KYi0TKFARApEpFVEdvoca+pxzl0i8ng/a3xeRI6KyOthFFUZJVgjLYCihIhWY8xCf0+4i43EGOP093wUcKyvvQWCMeZ3IlIFfDt0IimjFbUUlFGJ+w77kIg8D+wF8kTkOyKyxd3j/998zn1QRA6LyEci8pKIfNt9vFhECt2/Z7oLkhCRGPe8AM9af+c+XuS+5o/uWQIveHpaiciFIvKxiOwSkc0iktLXOkPct6+11CoiVwx8laJ0oZaCMlpI8HGxnMDVdG8mcKcxZqOIXO1+vBQQYLV77kMzcBuu6nIrsB1Xf/z+uBtoNMZcKCJxwHr3vAhwdcw8D1dbj/XAMhHZDPwO+LwxZouIpAKtfa1jjDkRxF7B1VxxNYDHohCRG3G1mfh4gLUUpRuqFJTRQjf3kTumUGKM2eg+dLX7Z4f7cTIuJZEC/MnTn0dEVgfwWlcDF4jI37gfp7nXsgObjTFl7rV2AgVAI1BpjNkC4Oki6lZU/tYZSCn03OtdQKHP45nAz4AVxpiOAPajKF5UKSijmWaf3wX4iTHm174niMg3+7m+ky4Xa3yPte43xrzVY60ioN3nkIP+P2N+1xkK7nbOvwfu9ZnjoCgBozEFZazwFvAV95cmIpIjIhOBD4FbRCTBPVrxRp9rTgJL3L//TY+1vuZun4yIzBKRpH5e+xAw2TNMyR1PsA5inUB4GnjGGLNuiOsoYxS1FJQxgTHmbRGZC2xwx36bgC+65+z+DtiFa9bBFp/LHsE1AGkl8IbP8SdxuYW2uwPJZ+hnzKYxxi4inwd+KSIJuOIJVwW7zkCISD4u5TVLRL7iPnyPMWbrYNdUxh7aJVVRfBCRHwJNxphHIi0LeGMjrxtj5g9xnSJcLdpvCIFYyihG3UeKMrJxAGk9so2Cwm2l/A9QHyqhlNGLWgqKoiiKF7UUFEVRFC+qFBRFURQvqhQURVEUL6oUFEVRFC+qFBRFURQv/x+eVRMW73EQWwAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "mask = np.abs(freq-20)<5\n", "\n", "fig = figure()\n", "ax = fig.subplots(1, 1)\n", "\n", "ax.semilogy(freq[mask], np.abs(amplitude_tilde[mask])**2)\n", "ax.grid()\n", "\n", "i = np.argmax(np.abs(amplitude_tilde[mask])**2)\n", "freq_max = freq[mask][i]\n", "amp_max = np.abs(amplitude_tilde[mask][i])**2\n", "\n", "ax.semilogy(freq_max, amp_max, 'ko')\n", "ax.text(freq_max, amp_max, f' {freq_max:.2f} Hz')\n", "\n", "ax.set_xlabel('Frequence [Hz]')\n", "ax.set_ylabel('Puissance [u.a]')" ] } ], "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 }