diff --git a/whiten_signal.ipynb b/whiten_signal.ipynb index 4e5cc20c35c312843aeed9436cb57443c766a8e6..256190686fbd1e0c3cb7fceb0c7c1dfec4e11225 100644 --- a/whiten_signal.ipynb +++ b/whiten_signal.ipynb @@ -46,7 +46,7 @@ { "data": { "text/plain": [ - "[]" + "[]" ] }, "execution_count": 3, @@ -55,7 +55,7 @@ }, { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYAAAAEWCAYAAABv+EDhAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3de5xcdX3/8ddnZnZ2N3vLbXMnd+4h3CKgVMWKwk+5qHjBolYt+LPWUov+bKmtrba/1ksVtK211FrrpVChUsU7QS4KBki4hXCJJEBIyGWTkNtuMrMz8+kf58zs7P3MZndnzu77+XjsY+dyds6HYXI+8/1+vhdzd0REZPJJVDsAERGpDiUAEZFJSglARGSSUgIQEZmklABERCYpJQARkUlKCUBEZJJSAhCpgJk9Z2bn93ksbWa3hM+5mZ1XpfBEKqIEIDI6fgW8C9hR7UBEokpVOwCRuHP3LHA9gJnlqxyOSGRqAYiITFKxSwBm9nUz22Vmj4/Ca51mZr82sw1m9piZvWOAY75sZoeO9lwiIrUmdgkA+AZw4Si9VhfwHnc/OXzN681savFJM1sFTBulc4mI1JTYJQB3vwfYW/6YmS0zs5+a2Toz+6WZnRDxtTa6+2/C2y8Cu4D28DWTwOeBj4/qf4CISI2YKEXgG4APuvtvzOxs4CvAb1fyAmZ2FpAGNoUPfRj4gbtvN7NRDVZir87MGsru54AkUPygpMPnM6711qWGxT4BmFkz8Arg5rILdX343FuATw/wZ9vc/YKy15gLfAv4XXcvmNk84G3AeWMYusTXj/vc//8EQ0AXhfd/Fv5eAjw3TjGJVMzi+AXFzBYDP3T3FWbWCjzt7nNH+FqtwF3A37r7LeFjbwT+DTgSHrYQ2Ozuy48ydBGRmhG7GkBf7n4AeNbM3gZggVOj/K2ZpYFbgW8WL/7ha/7I3ee4+2J3Xwx06eIvIhNN7BKAmd0I/Bo43sy2mtnvAVcAv2dmjwIbgEsjvtzbgVcB7zWzR8Kf08YkcBGRGhPLLiARETl6sWsBiIjI6IjVKKCZM2f64sWLqx2GiEisrFu3bre7t/d9PFYJYPHixaxdu7baYYiIxIqZPT/Q4+oCEhGZpJQAREQmKSUAEZFJSglARGSSUgIQEZmklABERCYpJQARkUlKCUBkAti27zB3Pr2r2mFIzCgBiEwA3/z1c/zhfz5c7TAkZpQARCaAbK5Ad75Q7TAkZpQARCaAQsEpaGVfqZASgMgEkCs4BV3/pUJKACITQMHVApDKKQGITAD5guMO2uBJKqEEIDIB5ML+H3UDSSWUAEQmgEIpASgDSHRKACITQD687isBSCWUAEQmgGILQNd/qYQSgMgEkCsEk8DyKgJIBZQARCaA4iRgdQFJJZQARCaA4oVfDQCphBKAyASQK9UAlAEkOiUAkQmgWARWDUAqoQQgMgHkNRFMRkAJQGQCyLu6gKRySgAiE4BaADISSgAiE0AxAeTVApAKKAGITAClYaBqAkgFlABEJoBcXktBSOWUAEQmgJ6JYMoAEp0SgMgEoBqAjIQSgMgEkNdMYBmBqicAM0ua2cNm9sNqxyISV3mtBSQjUPUEAPwR8GS1gxCJs7x2BJMRqGoCMLMFwBuBr1UzDpG4K20JWahyIBIr1W4BXA98HBj0Y2tmHzCztWa2tqOjY/wiE4mRnFoAMgJVSwBmdhGwy93XDXWcu9/g7qvcfVV7e/s4RScSLxoGKiNRzRbAucAlZvYccBPw22b27SrGIxJbWgtIRqJqCcDdr3X3Be6+GLgc+IW7v6ta8YjEmbqAZCSqXQMQkVHQUwRWApDoUtUOAMDd7wLuqnIYIrGleQAyEmoBiEwAmgcgI6EEIDIBKAHISCgBiMScu5e6fjQRTCqhBCASc+X9/moBSCWUAERiLlf2tV8JQCqhBCASc+XdPrr+SyWUAERirnwTmLzGgUoFlABEYi6f77noqwtIKqEEIBJz5S0ANQCkEkoAIjFX3u2jLSGlEkoAIjFXngC0KbxUQglAJObUBSQjpQQgEnMFdQHJCCkBiMRceReQRgFJJZQARGIuV14D0FpAUgElAJGYK7haADIySgAiMadhoDJSSgAiMde7BlDFQCR2lABEYq7XPABlAKmAEoBIzJXPA1AXkFRCCUAk5grqApIRUgIQibmc5gHICCkBiMScWgAyUkoAIjHXay0gZQCpgBKASMypC0hGSglAJObUBSQjpQQgEnNaDE5GSglAJOZ6JQA1AaQCSgAiMacNYWSklABEYk5dQDJSSgAiMVfQUhAyQkoAIjGXy2tTeBkZJQCRmCuoBiAjpAQgEnPl20CqBiCVUAIQibl8oScD6PovlVACEIk5bQgjI1W1BGBmx5jZnWb2hJltMLM/qlYsInFWVgNWF5BUJFXFc+eAj7r7Q2bWAqwzs9vd/YkqxiQSO8UuoFTC1AUkFalaC8Ddt7v7Q+Htg8CTwPxqxSMSV8UicCppagFIRWqiBmBmi4HTgfsHeO4DZrbWzNZ2dHSMd2giNa940a9LJlQDkIpUPQGYWTPw38BH3P1A3+fd/QZ3X+Xuq9rb28c/QJEaV7zo1yUTmgcgFalqAjCzOoKL/3fc/XvVjEUkroobwgQ1AGUAia6ao4AM+DfgSXf/YrXiEIm7QsFJGCQTqgFIZarZAjgXeDfw22b2SPjzhirGIxJLeXdSiQQJs16zgkWGE2kYqJlNBd4DLC7/G3e/eqQndvdfATbSvxeRQL7gJBJgFqwG6u5cd/tG3v6yY1gwbUq1w5MaFnUewI+BNcB6QN8xRGpIvuAkzUhY0AW0+1CWL//iGaY1pXnfuUuqHZ7UsKgJoMHdrxnTSERkRPIFJ5mwsAbQMyook9N3NRla1BrAt8zsKjOba2bTiz9jGpmIRFJMAGZBPSAXzgzOdCsByNCitgCywOeBTwDFYQYOLB2LoEQkurwHCSBhwTDQ4uKgR3L56gYmNS9qAvgosNzdd49lMCJSuUKhmACgUEAtAIksahfQM0DXWAYiIiPTtwjcUwNQC0CGFrUF0Ak8YmZ3Apnig0czDFRERkcwDLSYAHr2BT6iFoAMI2oC+J/wR0RqTKkGkAgWhituEq8WgAxn2ARgZkng9e5+xTjEIyIVyhd6isC9u4DUApChDVsDcPc8sMjM0uMQj4hUqFgDsLALqLg43JFutQBkaFG7gDYD95rZDwjqAQBoETeR6itNBAuXgiguCKcWgAwnagLYFP4kgJaxC0dEKlUomweQL5TVANQCkGFESgDu/ikobd6Cux8ay6BEJDrVAGSkIs0DMLMVZvYwsAHYYGbrzOzksQ1NRKLIFZyEBUtBlA8DVQKQ4USdCHYDcI27L3L3RQQzg/917MISkaiKXUDJcEewfDgTWEVgGU7UBNDk7ncW77j7XUDTmEQkIhUp7wLqVQNQC0CGEXkUkJn9BfCt8P67CEYGiUiV5QvBjmDFLqDSKCC1AGQYUVsA7wfage8RbOI+M3xMRKqsvAXg7j3zANQCkGEM2QIws2+5+7uB92jdH5HalHcG3BAm6A4qkEpWc+tvqWXDfTLONLN5wPvNbFr5ZjDaEEakNuQLhdJy0OU1AFAdQIY2XA3gq8AdBBu/PNTnOW0II1ID8gXKhoF6aRgoBCOBmuqjlvpkshmyBeDuX3b3E4Gvu/uSPj+6+IvUgELBSYUtAC/rAgK1AGRokToH3f33zey3zOx9AGY208yWjG1oIhJF3ss3he8pAoPmAsjQos4E/kvgT4Brw4fSwLfHKigRia64IYyZkXenoBaARBR1eMCbgUsIVwJ19xfRonAiNSFYDppwGCi9WgBKADKUqAkg6+5OUPjFzDQLWKRGBPMAEsGm8GVLQYC6gGRoURPAd83sX4CpZnYVsBqtBSRSE4IEAMnSaqA9z6kFIEOJuhz035vZ64ADwPHAJ9399jGNTEQiKRaB8wUoFOjVAtByEDKUqEXga4An3P3/ufvHdPEXGV55MfbGB7Zw4wNbAPjlbzq44Lp7eHrHQQC27z/M/q7uozpPcSJYv1FAagHIEKLOEGkBfm5me4H/Am52951jF5ZI/HQczNCVzbFoRhM/Wb+dj938KB989TLmT2vk2u+tB2DXgQz/ft+z7Ovq5trvPcZnL1vJm79yH6ce08Z3rjxnROfNhXsC990QBtQCkKFVsiPYp8xsJfAO4G4z2+ru549pdCI1rCub49aHt3HK/DYSZrzvGw+ytzPLxSvn8qP122lpqOMLt28E4Jyl06lPJblu9UZaGlJc/dpj+fIdv+EtX7mPQ5kc9z6zh80dh1ja3lxxHIVwGGiiz1pAoBqADK3SOeK7gB3AHmDW6IcjUnuOdOdJJxNk8wX++a5NPLH9AK86dibfWvM8G3cGu6PWJY325nouO2M+N6/bysnzWvnOledw89oXuHtjB1++/HTSqQSfum0Dl542n1csm8H9m/ew7vmX+MoVZ3D1jQ9z4wNbuOZ1x7N59yFOntcWOb68F1sAQTLIayKYRBQpAZjZh4C3EywJfTNwlbs/MZaBiVRDvuDc+8xuntpxgIXTm7j3md3c9OAWGuuSTEmn2HHgCLNa6rn9iZ1Mb0rzL+8+k+f3dPLwln188uKTmNvWyO+ft5zZrfVMSae48pVLufKVPaumfO6tp5Zuf+13V7Fj/xGOnd3CDx97kZvXbeWejbvZuOsgd1zz6sitgXzBSSaNRK5nJnBd0ujOu1oAMqSoLYBjgI+4+yNjGYzIWHJ3DmZy7DpwhJ0HMuwMf6/fto/7Nu0hnw+6UvYf7inIphLGZWcsIJGAbfuO8IW3n8orls3g6Z0Hmd3SwLSmdL/zLJkZbZpMS0MdLQ11APzOWYv48fod5PKOOzz43N7KEkCxBRB2ATXUJckVcmS68+ztzHLgcDeLI8Ylk0fUGsC1AGY2C2goe3zLGMUlUpHOTK50Qd918Ai7ihf4g8Hv4kX/8ABdInPbGnjdibNpbkjRlcnz6uPbOWfpDLbs7WJ2az1z2xr7/c0Jc1pHNf5zl8/gS5efxqrF07noy79k7XMv8Y6XLYz0t8VhoIlETxG4LpmgPpUgkyvwuZ8+xU837GDNta+loS45qnFLvEXtAroY+CIwj6AOsAh4Ejh57EKTya5QcPYd7mZvZ4bdh7LsOpgJL+RH2FW6sAe/O7P9L+yNdUlmtdYzu7WBUxZM5fyW4HbxsdmtDcxqqR90ueTpA3y7HytmxqWnzQfgzEXTWLflpUh/VygELYaeHcHCUUEJoz6VJJMrsPWlw+zr6uZnG3aUzuEedA8pIUxuUbuA/gY4B1jt7qeb2WsI9gU+KmZ2IfAlIAl8zd0/c7SvKbXJ3TncnWf/4W4OHM6xryvL3s4sezqz7DmUDS7ynVn2HsqypzPD3s7g+bJ6ZklDXaJ08T5xXivnHT+L2a31wYW9pYFZ4UW+pT6FmY3/f+xROmPRNFY/uYu9ndlhk1Bx7f9iF1C+ECwGlzSjoS7Bke48HQczQDAXoZgAvn3/Fj73k6e4++Ovobk+xQXX38N7Xr6I9527hGd3d/LQ8y9x2ZkLAMjmCqRT2lVsIoqaALrdfY+ZJcws4e53mtn1R3NiM0sC/wS8DtgKPGhmP1BxufZkcwUOZ/N0ZnN0ZXN0ZsLbxd/ZPJ2ZHPsPdw/4cyC86GfzgxckWxtSzGyuZ3pTmiUzmzhz0XRmNqeZ3pRmRnM9M5rSzGqpZ1ZrA60N8bywR7VqUbDZ3kPPv8T5J80e8tjiiJ9EomceQLEFkEwEXUAdhzKkUwnWbN7Ls7s7WTh9Cv96z2YOZnL8aP12ZrfU8+zuTv793uf43Zcv5tO3beDOpzs4fk4LU6fUcdE//IoPv2Y5V75yKbc+vJWHnt/HX11yMrlCgW+v2cJFK+cyu7WBF/Z20ZXNc/ycYJ3IbfsOM6+tIViltOBYGKfUjqgJYJ+ZNQO/BL5jZrsIVwY9CmcBz7j7ZgAzuwm4FBj1BHDHkzt55IV9QLBhRpHTc8f7fNMsvzvY3wxyE+/zYj7ocdFi6f1aZcdV+FoO5PNONl8gmyuQyRXC23m68042FzyezRfIdOfpzObpyubozg8RUBkzaG2oo62x52deWyOtjb0fK/5Mb0ozsznNtKY0ddq3tmTlgjZSCWPdluETQKHYAghrAMGGMAVSSSOdTHAok2NvZ5Z3nrWQ7659ga/c+QwXrpjDlr1dpFMJvv/wNma11gOwZW8Xt6zbyl0bOwD457s30ViXZF9XN1/4+UZOmtfKn/73ejK5AvOmNrJlbyc3PvACN699gb950wqu/OZaujJ5rr/8NNZs3sM3f/08F586j7evWsBf/M/jAFz7hhNZ+9xefrZhJ29cOZeT5rby3bUv0FiX5A2nzOWpHQd5escBzlw0jbYpaR59YR9N6SRL25vZ05llX1eWaVPS1CWNzmyegnuwGxrB5694eyJ68xkLIg8wiGq4TeE/AtwHvAnoAj4CXAG0AZ8+ynPPB14ou78VOHuAGD4AfABg4cJoRbG+7t7YwbfXPF/+mj23e52rz7nLnx34Zq+/KT++/2uV/83A5x/8HL1fzCIcN1iMyYSRTiVIJxOkU8nwdtBd0NqQCu6nktSnEjSlk0ypTwW/0yma6vv8TqeYUp8s/W5Op/QNbxQ01CU5eX4bv/rNbj72+uDb/GCKyz6k+iwFEXQBJdn20mEAVsxvpa1xKV+9exN3PLWLWS31XHH2Iq5bvZG6pHH5y47hh49t58+//zhJMy49fT63PrwVgDeunMvtT+zk3f/2AE3pJC9fNoPP/ewp3OGCk2dzx5O7eOtXf83s1noWTZ/Ch74T7B573vHt/OixF7nt0RdZMK2RVML4v99ahxmcuXAaX717E+4wf2oj2XyBnz+xk1TCWDhjCnc+HSShGU1pjoRfRhIWjJw6cKQb9/C/OWHgwX+305MQJ6IzFk0b3wQALACuB04A1gP3EiSE29x976hGMgh3vwG4AWDVqlUj+r/76UtX8OlLV4xqXCJj6R2rjuHPbl3Px295jM+/deWgibW43lAiXAoi707Bi0XgBBt3BusNtTfX886XLWTXwSN876Ft/PH5x/GWM+Zz3eqNdOedd54VfLm66cEXuPjUefzJ/zme2x57kXQywV9fuoKlM5v4h188w59fdBIXnDSHi//xVxw7q5mvXHEmtz36Il+/91mue8dpzGlt4K9/+ARnLprG21Ydw/2b93Dfpj1c9aqlpBLGTx7fzglzWjlxbiubOg6xfd8RXr5sBu7OY9v2s3RmE1OnpNlzKENXNs+CaY24Q8ehDNOmpEmnEuTyBfLu1KdUwD5aQyYAd/8YgJmlgVXAK4D3ATeY2T53P+kozr2NYH5B0YLwMZFJ73fOXkjHwQzXrd7ICXNauOpVA2/BXawBJMMdwdyhOx8mgLoEB47kAGhvqSeRMD572UrOO34Wrz9pNg11Sc5aMp29nVlWLmgjmVjEDx/bzlWvXMKslgb+9s2nMCWdZHpTmo+cfxznnziblQvaMDNWX/Nq6pLBOd90+nzedPr8UkyfuWxl6fbZS2dw9tIZpftvPn1B6fay9maWleY6GGcsnFZ6bkZzPcW/MoPZraXR56SSiYqXMJCBRX0fG4FWgq6fNuBFghbB0XgQODbcW3gbcDnwO0f5miITxtWvXc6tD28t1a8Gki+rASTDvr5cvkAyYTSUfUNubwn6+euSCS45dV7p8X++4oygQGvGivltPP6pC0rPvfXMnot1MmGceszU0n2NCpoYhqsB3EAw1v8gcD9B988X3T3aIOUhuHvOzD4M/IxgGOjX3X3D0b6uyERhZsyb2sj2/YcHPaa8BVDsJcoVnFTYAiia2Vw/4N/PGORxmRyGawEsBOqB3xB8S98KDP51pELu/mPgx6P1eiITzZy2BtZs2jPo86UEYFaqE2RyhdJEMICWhpQmfMmAhqsBXGjB0JKTCfr/PwqsCPcF+LW7/+U4xCgyac1ta2DnwUy47WP/QnBx86+gBhDc7s4XSCWCkV3Q0/0j0tewHXkeeJzgm/pPCEYCLQP+aIxjE5n05rQ1ki84uw9lBnw+F2aA3jUA79UCaFc3jwxiyARgZleb2U1mtgW4G7gIeAp4CzB9HOITmdTmtQWjX7bvPzLg88Vx78WZwFBsAQSLwYFaADK44WoAiwnW//9jd98+9uGISLk5YQLYsf8wlI3CKSqurpG0ni6gbLEGEPb7KwHIYIarAVwzXoGISH/FpagHawH0HgUUZIBsvlCaCAZKADI4DeYVqWHTptSRTiXYESEBFIvEPTWA4J/3YENARZQARGqYmTG3rWHwFkBpIhileQA9o4DUBSRDUwIQqXFzWhsitAASpcUAs/kCifIuILUAZBBKACI1bm5bA9sPDDwbuNdEsF6jgIxzl8/knWct5LjZLeMWq8SLEoBIjZvT1sjO/ZnSyp/lejaECbqBoGcxuHlTG/m7t5yidXtkUPpkiNS4uW0NZPMF9nZl+z1XKNsSstgFlA/3AxAZjhKASI3rmQvQvw6QG2AYKEAqqQQgw1MCEKlx88K5AFtf6l8HyOaCmWD1qSSJPju/iQxHCUCkxi2cMQWA5/f034a7lADqEr1aAOoCkiiUAERqXFtjHTOa0jy7u38CyOTyANSnEr22jUwm9E9bhqdPiUgMLJnZxOYBE0DQAkinEr26gFQDkCiUAERiYMnMJp4bIAH0rgGUtwCUAGR4SgAiMbCkvYldBzMcyuR6Pd6rC6i8CKwagESgBCASA0tnNgH0awVkusu7gNQCkMooAYjEwJKZzQD96gCZXIGEQarvPAAlAIlACUAkBhbNmIIZPNvROwFk8wXqU0nMjPKBPwklAIlACUAkBhrqksxra+TZ3Yd6PZ7pzpfW+jG1AKRCSgAiMbG0vanfXIBMrlBa9jmpGoBUSAlAJCaKcwHce1YFzeYK1NcF/4xVBJZKKQGIxMSCaY0cPJLrNRQ0aAEEO3/1mgimBCARKAGIxERzfR0AnZl86bFMLk862b8GoKUgJAp9SkRioqk++KbfrwUQdgGVd/uoBSBRKAGIxERzfQqAzn5dQMUaQM+xGgYqUSgBiMRE0yAJIB3WADQMVCqlBCASE8UWQK8uoO78gC0AjQKSKJQARGKiaYAEEMwE7l8DUAKQKJQARGKiWATu1QXUXT4MVAlAKqMEIBITPV1A5cNAC2VLQfQcqxqARKEEIBITjXXBxu/lLYBsrrwG0HPR1yggiUIJQCQmzIym+pTmAcioUQIQiZHm+lSpBeDuQQJIahSQjExVEoCZfd7MnjKzx8zsVjObWo04ROKmqT5FZzZIAN35YFG4+rr+8wC0JaREUa0WwO3ACndfCWwErq1SHCKxEnQBBUXg8v2AoXcNIJVUApDhVSUBuPvP3b3YkbkGWFCNOETiprk+WeoCyuR69gOGvvsBqHdXhlcLn5L3Az8Z7Ekz+4CZrTWztR0dHeMYlkjtaUqnOHSkdwKo1zBQGaHUWL2wma0G5gzw1Cfc/fvhMZ8AcsB3Bnsdd78BuAFg1apVPthxIpNBc9kooGwpAYQTwcou+gnVACSCMUsA7n7+UM+b2XuBi4DXevkWRyIyqPIicP8aQM9xqgFIFGOWAIZiZhcCHwde7e5d1YhBJI6ayoaBZrqHqgEoAcjwqlUD+EegBbjdzB4xs69WKQ6RWGlpSNGddzK5PNl87y4gDQOVSlWlBeDuy6txXpG4a0oXF4TLl1oAPZvC9xynFoBEUQujgEQkovJNYYo1gHRS8wBkZJQARGKkfFOY0jDQYgsgoS4gqYwSgEiMlLcA+g0DVReQVEgJQCRGmnq1AMIuoIGWgtBMYIlAnxKRGBmwC2jA/QDGPzaJH31MRGKkfFvIbN8EUPavWS0AiUKfEpEYKd8WMtOvBqCJYFKZqswDEJGRKS8C58KJYHXhkM/eNQAlABmeEoBIjNQlE9SnEsFyEBZ0/xRnAJdf87UnsEShBCASM8UVQYvJoMjMMNMcAIlONQCRmCkuCJfJFUiH/f9FCTP1/0tkagGIxExxW8hEwnq1ACDoBlICkKjUAhCJmeb6JIcy3WRzhdIyEEWmFoBUQAlAJGZmtTSwY/8RMrlCaQhoUdJMI4AkMiUAkZhZ1t7Elr1dHDzSXVoGokhdQFIJJQCRmFk2q5mCw9M7Dg5QA1AXkESnBCASM8vamwF4qau7XwIw0zIQEp0+KSIxs7S9qXS7XwsgYVoITiLTR0UkZqakU8yf2ggwSBFY/6wlGn1SRGJo2aygG6h/F5BqABKdEoBIDC0P6wADjQLSMFCJSglAJIaWzQrqAAONAkpoLSCJSAlAJIaKI4Hq6/rUABJGKqkEINEoAYjE0PKwBpBO9h8GqhqARKUEIBJDM5rSfPg1y7lwxZxejyfMtBy0RKbVQEViyMz42AXH93tcS0FIJdQCEJlAEqoBSAWUAEQmEI0CkkooAYhMIJoHIJVQAhCZQILVQPXPWqJREVhkAvnQa5YztbGu2mFITCgBiEwgl5w6r9ohSIyorSgiMkkpAYiITFJKACIik5QSgIjIJFXVBGBmHzUzN7OZ1YxDRGQyqloCMLNjgNcDW6oVg4jIZFbNFsB1wMcBr2IMIiKTVlUSgJldCmxz90cjHPsBM1trZms7OjrGIToRkcnB3MfmC7iZrQbmDPDUJ4A/A17v7vvN7DlglbvvjvCaHcDzIwxpJjDsOWqI4h1bindsKd6xVWm8i9y9ve+DY5YABmNmpwB3AF3hQwuAF4Gz3H3HGJ53rbuvGqvXH22Kd2wp3rGleMfWaMU77ktBuPt6YFbxfiUtABERGT2aByAiMklVfTE4d188Tqe6YZzOM1oU79hSvGNL8Y6tUYl33GsAIiJSG9QFJCIySSkBiIhMUhMuAZjZhWb2tJk9Y2Z/OsDz9Wb2X+Hz95vZ4vGPslc8w8X7KjN7yMxyZvbWasTYJ57h4r3GzJ4ws8fM7A4zW1SNOMviGS7eD5rZejN7xMx+ZWYnVSPOsniGjLfsuMvCdbSqOnQxwvv7XjPrCN/fR8zsymrEWRbPsO+vmb09/AxvMLP/HO8Y+8Qy3Pt7Xdl7u9HM9lV0AnefMD9AEtgELAXSwKPASX2O+f/9KkUAAATGSURBVBDw1fD25cB/1Xi8i4GVwDeBt8bg/X0NMCW8/fsxeH9by25fAvy0luMNj2sB7gHWEAyhrtl4gfcC/1itGEcQ77HAw8C08P6sWo63z/F/CHy9knNMtBbAWcAz7r7Z3bPATcClfY65FPiP8PYtwGvNzMYxxnLDxuvuz7n7Y0ChGgH2ESXeO929OMlvDcFEv2qJEu+BsrtNVHdtqiifX4C/Bj4LHBnP4AYQNd5aESXeq4B/cveXANx91zjHWK7S9/edwI2VnGCiJYD5wAtl97eGjw14jLvngP3AjHGJrr8o8daSSuP9PeAnYxrR0CLFa2Z/YGabgM8BV49TbAMZNl4zOwM4xt1/NJ6BDSLq5+GysEvwlnAV4GqJEu9xwHFmdq+ZrTGzC8ctuv4i/3sLu1qXAL+o5AQTLQFIjTCzdwGrgM9XO5bhuPs/ufsy4E+AP692PIMxswTwReCj1Y6lArcBi919JXA7Pa3vWpUi6AY6j+Ab9b+a2dSqRhTN5cAt7p6v5I8mWgLYBpR/w1gQPjbgMWaWAtqAPeMSXX9R4q0lkeI1s/MJFv27xN0z4xTbQCp9f28C3jSmEQ1tuHhbgBXAXeESKucAP6hiIXjY99fd95R9Br4GnDlOsQ0kyudhK/ADd+9292eBjQQJoRoq+fxeToXdP8CEKwKngM0ETaFi0eTkPsf8Ab2LwN+t5XjLjv0G1S8CR3l/TycoXB0bk8/DsWW3LwbW1nK8fY6/i+oWgaO8v3PLbr8ZWFPj8V4I/Ed4eyZBF8yMWo03PO4E4DnCib0VnaNa/zPG8E17A0HW3gR8Inzs0wTfRgEagJuBZ4AHgKU1Hu/LCL6VdBK0VDbUeLyrgZ3AI+HPD2o83i8BG8JY7xzqglsL8fY5tqoJIOL7+3fh+/to+P6eUOPxGkE32xPAeuDyWo43vP9XwGdG8vpaCkJEZJKaaDUAERGJSAlARGSSUgIQEZmklABERCYpJQARkUlKCUAmLTObUbaS4g4z21Z2/75RPM+bzOyTQzx/ipl9Y7TOJxKVhoGKAGb2V8Ahd//7MXjt+wjGbe8e4pjVwPvdfcton19kMGoBiAzAzA6Fv88zs7vN7PtmttnMPmNmV5jZA+E+AsvC49rN7L/N7MHw59zw8eOATPHib2ZvM7PHzexRM7un7JS3EcxMFxk3SgAiwzsV+CBwIvBu4Dh3P4tgbZs/DI/5EnCdu78MuCx8DuBc4KGy1/okcIG7n0qw/0DRWuCVY/ZfIDKAVLUDEImBB919O0C4bPTPw8fXE2yAA3A+cFLZ1hKtZtYMzAU6yl7rXuAbZvZd4Htlj+8C5o1N+CIDUwIQGV75iqaFsvsFev4NJYBz3L3XJi1mdphgxVkA3P2DZnY28EZgnZmd6e57CNaoOjxG8YsMSF1AIqPj5/R0B2Fmp4U3nwSWlz2+zN3vd/dPErQMisv9Hgc8Pk6xigBKACKj5WpgVbjz1RMENQMI9u49vWzb0c+HxePHgfsIVsmEoCupFnb5kklEw0BFxpiZfQm4zd1XD/J8PXA38FsebFMqMi7UAhAZe38LTBni+YXAn+riL+NNLQARkUlKLQARkUlKCUBEZJJSAhARmaSUAEREJiklABGRSep/AXP7pYHR7mr1AAAAAElFTkSuQmCC\n", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYAAAAEWCAYAAABv+EDhAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3deZwcdZ3/8denu+fKXLlvkpBwE8IVAWUVXBBQuRREFHXFBVddZV3g5y66y6q7D1fRFUUXXTzWVVhQEFYODw455AiQcIVwRBIgJCRkkpBrJpmZ7v78/qjqmZq7ejIz3TXzfj4e80h3VXXVZzo99envbe6OiIiMPalSByAiIqWhBCAiMkYpAYiIjFFKACIiY5QSgIjIGKUEICIyRikByKhkZl80sx8P9bExzuVmtk+RrznPzO4ciusPcJ3jzWztcF9HksM0DkDKnZl9HLgEWABsB24BLnP3raWMqzdm5sC+7v5SqWPpzsyOB65199mljkXKg0oAUtbM7BLgG8D/AxqBY4C5wF1mVtnHazIjF6FIcikBSNkyswbgK8Dn3P337t7u7q8A5wDzgI+Ex33ZzG4ys2vNbDvw8XDbtZFzfczMXjWzzWb2z2b2ipmdGHn9teHjeWE1zl+Z2Roz22RmX4qc5ygze8TMtprZejP7fl+JqJff5+NmttrMdpjZy2Z2XmT7g5HjTjKzF81sm5ldbWb3m9kF0WPN7Ftm9mZ4nndHXnu+mT0fXmO1mf3N4N59GQuUAKScvQ2oBm6ObnT3ncBvgXdFNp8B3ASMB66LHm9mBwFXA+cBMwhKErMGuPZfAPsDJwCXm9mB4fYc8PfAZOCt4f7PDPSLmFktcBXwbnevD3+3p3o5bnL4e1wGTAJeDI+NOjrcPhm4AviJmVm4byNwKtAAnA9caWZHDBSfjE2JSwBm9lMz22hmzw7BuQ4Lv82tMLNnzOyDvRxzlZnt3NNryaBMBja5e7aXfevD/QWPuPv/uXve3Xd1O/Zs4DZ3f9Dd24DLgYEav77i7rvc/WngaeBQAHdf5u5L3D0blkb+Czgu5u+TBxaaWY27r3f3Fb0c8x5ghbvfHP7eVwEbuh3zqrv/yN1zwP8QJLVpYXx3uPsqD9wP3Am8PWZ8MsYkLgEAPwNOGaJztQAfc/eDw3N+x8zGF3aa2WJgwhBdS4q3CZjcR53+jHB/wWv9nGdmdL+7twCbB7h29KbbAtQBmNl+Zna7mW0Iq5u+RtdE1Ct3bwY+CHwKWG9md5jZATFidaB7z50Nkf0t4cNCfO82syVmtsXMthIklAHjk7EpcQnA3R8AtkS3mdkCM/u9mS0zsz/18YfV27lWuvufw8evExSfp4TnTAPfBL4wpL+AFOMRoBV4f3SjmdUB7wbuiWzu7xv9eqCj54uZ1RBUrwzGD4AXCHr6NABfBKz/l4QBuv/B3d9FkLxeAH4UI1aLPu+PmVUBvwa+BUxz9/EEVWWx4pOxJ3EJoA/XEDQUHglcSlDfWxQzOwqoBFaFmz4L3Oru64csSimKu28jaAT+npmdYmYVZjYP+BXBt+JfxDzVTcBpZva2sMH2ywz+plhP0BV1Z/hF49NxXmRm08zsjLAtoBXYSVAl1N0dwCFmdmZY8vlbYHrM2CqBKqAJyIaNwyfFfK2MQYlPAOG3wbcBN5rZUwR1sjPCfe83s2d7+flDt3PMILiZnO/ueTObCXwA+N7I/jbSnbtfQfAt+1sEN95HCapITnD31pjnWAF8DriB4Bv2ToLSXqzXd3Mp8GFgB8E3+F/GfF0KuBh4naAEexy9JA9330Tw2buCoJrqIGBpnFjdfQdwEUGCfDOM89aY8ckYlMiBYOG3wNvdfWHYVfBFd58xyHM1APcBX3P3m8Jt7wV+AuwOD5sDrHb3okZ4SnkKvzRsJajGebnU8fTHzFIEpZ3z3P3eUscjo0viSwDuvh142cw+AEGdqZkdGue1YXXALcDPCzf/8Jx3uPt0d5/n7vOAFt38k83MTjOzcWEVzLeA5cArpY2qd2Z2spmND+v0C20MS0ocloxCiUsAZnY9QePg/ma21sz+mqB/91+b2dPACoI+4XGcA7yDYODQU+HPYcMSuJTaGQTVL68D+wLnevkWf99K0Ba1CTgNOLOXrq0ieyyRVUAiIrLnElcCEBGRoZGoSbMmT57s8+bNK3UYIiKJsmzZsk3uPqX79kQlgHnz5rF06dJShyEikihm9mpv21UFJCIyRikBiIiMUUoAIiJjlBKAiMgYpQQgIjJGKQGIiIxRSgAiImOUEoDIKLBu6y7ufXFjqcOQhFECEBkFfv7IK3zuf58sdRiSMEoAIqNAWzZPe663BcZE+qYEIDIK5PNOXjP7SpGUAERGgWzeyev+L0VSAhAZBfKuEoAUTwlAZBTI5R130AJPUgwlAJFRIBvW/6gaSIqhBCAyCuQ7EoAygMSnBCAyCuTC+74SgBRDCUBkFCiUAHT/l2IoAYiMAtl8MAgsp0YAKYISgMgoUBgErCogKYYSgMgoULjxqwAgxVACEBkFsh1tAMoAEp8SgMgoUGgEVhuAFEMJQGQUyGkgmAyCEoDIKJBzVQFJ8ZQAREYBlQBkMJQAREaBQgLIqQQgRVACEBkFOrqBqgggRVACEBkFsjlNBSHFUwIQGQU6B4IpA0h8SgAio4DaAGQwlABERoGcRgLLIJQ8AZhZ2syeNLPbSx2LSFLlNBeQDELJEwDwd8DzpQ5CJMlyWhFMBqGkCcDMZgPvBX5cyjhEkq5jSch8iQORRCl1CeA7wBeAPj+2ZvZJM1tqZkubmppGLjKRBMmqBCCDULIEYGanAhvdfVl/x7n7Ne6+2N0XT5kyZYSiE0kWdQOVwShlCeBY4HQzewW4AfhLM7u2hPGIJJbmApLBKFkCcPfL3H22u88DzgX+6O4fKVU8IkmmKiAZjFK3AYjIEOhsBFYCkPgypQ4AwN3vA+4rcRgiiaVxADIYKgGIjAIaByCDoQQgMgooAchgKAGIJJy7d1T9aCCYFEMJQCThovX+KgFIMZQARBIuG/narwQgxVACEEm4aLWP7v9SDCUAkYSLLgKTUz9QKYISgEjC5XKdN31VAUkxlABEEi5aAlABQIqhBCCScNFqHy0JKcVQAhBJuGgC0KLwUgwlAJGEUxWQDJYSgEjC5VUFJIOkBCCScNEqIPUCkmIoAYgkXDbaBqC5gKQISgAiCZd3lQBkcJQARBJO3UBlsJQARBKuaxtACQORxFECEEm4LuMAlAGkCEoAIgkXHQegKiAphhKASMLlVQUkg6QEIJJwWY0DkEFSAhBJOJUAZLCUAEQSrstcQMoAUgQlAJGEUxWQDJYSgEjCqQpIBksJQCThNBmcDJYSgEjCdUkAKgJIEZQARBJOC8LIYCkBiCScqoBksJQARBIur6kgZJCUAEQSLpvTovAyOEoAIgmXVxuADJISgEjCRZeBVBuAFEMJQCThcvnODKD7vxRDCUAk4bQgjAxWyRKAme1lZvea2XNmtsLM/q5UsYgkWaQNWFVAUpRMCa+dBS5x9yfMrB5YZmZ3uftzJYxJJHEKVUCZlKkKSIpSshKAu6939yfCxzuA54FZpYpHJKkKjcCZtKkEIEUpizYAM5sHHA482su+T5rZUjNb2tTUNNKhiZS9wk2/Ip1SG4AUpeQJwMzqgF8Dn3f37d33u/s17r7Y3RdPmTJl5AMUKXOFm35FOqVxAFKUkiYAM6sguPlf5+43lzIWkaQqLAgTtAEoA0h8pewFZMBPgOfd/dulikMk6fJ5J2WQTqkNQIpTyhLAscBHgb80s6fCn/eUMB6RRMq5k0mlSJl1GRUsMpBY3UDNbDzwMWBe9DXuftFgL+zuDwI22NeLSCCXd1IpMAtmA3V3rrxrJee8ZS9mTxhX6vCkjMUdB/BbYAmwHNB3DJEykss7aTNSFlQBbdrZxlV/fIkJtZWcf+zepQ5PyljcBFDt7hcPayQiMii5vJNOWdgG0NkrqDWr72rSv7htAL8wswvNbIaZTSz8DGtkIhJLIQGYBe0B2XBkcGu7EoD0L24JoA34JvAloNDNwIH5wxGUiMSX8yABpCzoBlqYHHR3NlfawKTsxU0AlwD7uPum4QxGRIqXzxcSAOTzqAQgscWtAnoJaBnOQERkcLo3Ane2AagEIP2LWwJoBp4ys3uB1sLGPekGKiJDI+gGWkgAnesC71YJQAYQNwH8X/gjImWmow0gFUwMV1gkXiUAGciACcDM0sBJ7n7eCMQjIkXK5TsbgbtWAakEIP0bsA3A3XPAXDOrHIF4RKRIhTYAC6uACpPD7W5XCUD6F7cKaDXwkJndStAeAIAmcRMpvY6BYOFUEIUJ4VQCkIHETQCrwp8UUD984YhIsfKRcQC5fKQNQCUAGUCsBODuX4GOxVtw953DGZSIxKc2ABmsWOMAzGyhmT0JrABWmNkyMzt4eEMTkTiyeSdlwVQQ0W6gSgAykLgDwa4BLnb3ue4+l2Bk8I+GLywRiatQBZQOVwTLhSOB1QgsA4mbAGrd/d7CE3e/D6gdlohEpCjRKqAubQAqAcgAYvcCMrN/Bn4RPv8IQc8gESmxXD5YEaxQBdTRC0glABlA3BLAJ4ApwM0Ei7hPDreJSIlFSwDu3jkOQCUAGUC/JQAz+4W7fxT4mOb9ESlPOafXBWGC6qA8mXQpl/6WcjbQJ+NIM5sJfMLMJkQXg9GCMCLlIZfPd0wHHW0DALUDSP8GagP4IXAPwcIvT3TbpwVhRMpALk+kG6h3dAOFoCdQbVXcpj4Za/otAbj7Ve5+IPBTd9+7249u/iJlIJ93MmEJwCNVQKASgPQvVuWgu3/azP7CzM4HMLPJZrb38IYmInHkPLoofGcjMGgsgPQv7kjgfwH+Abgs3FQJXDtcQYlIfIUFYcyMnDt5lQAkprjdA94HnE44E6i7v44mhRMpC8F00ITdQOlSAlACkP7ETQBt7u4EDb+YmUYBi5SJYBxAKlgUPjIVBKgKSPoXNwH8ysz+CxhvZhcCd6O5gETKQpAAIN0xG2jnPpUApD9xp4P+lpm9C9gO7A9c7u53DWtkIhJLoRE4l4d8ni4lAE0HIf2J2wh8MfCcu/8/d79UN3+RgUUbY69/bA3XP7YGgD/9uYmTr3yAFzfsAGD9tl1sa2nfo+sUBoL16AWkEoD0I+4IkXrgTjPbAvwSuNHd3xi+sESSp2lHKy1tWeZOquV3y9dz6Y1P86njFjBrQg2X3bwcgI3bW/nvh19ma0s7l938DN84axHvu/phDt2rkesuOGZQ182GawJ3XxAGVAKQ/hWzIthXzGwR8EHgfjNb6+4nDmt0ImWspS3LLU+u45BZjaTMOP9nj7OluY3TFs3gjuXrqa+u4D/uWgnAMfMnUpVJc+XdK6mvznDRCfty1T1/5v1XP8zO1iwPvbSZ1U07mT+lrug48mE30FS3uYBAbQDSv2LHiG8ENgCbgalDH45I+dndnqMynaItl+cH963iufXbece+k/nFkldZ+UawOmpF2phSV8VZR8zixmVrOXhmA9ddcAw3Ln2N+1c2cdW5h1OZSfGV21ZwxmGzeNuCSTy6ejPLXn2Tq887gouuf5LrH1vDxe/an9WbdnLwzMbY8eW8UAIIkkFOA8EkplgJwMw+A5xDMCX0jcCF7v7ccAYmUgq5vPPQS5t4YcN25kys5aGXNnHD42uoqUgzrjLDhu27mVpfxV3PvcHE2kr+66NH8urmZp5cs5XLTzuIGY01fPr4fZjWUMW4ygwXvH0+F7y9c9aUK84+tOPxj/9qMRu27WbfafXc/szr3LhsLQ+s3MTKjTu45+LjYpcGcnknnTZS2c6RwBVpoz3nKgFIv+KWAPYCPu/uTw1nMCLDyd3Z0Zpl4/bdvLG9lTfCf5ev28rDqzaTywVVKdt2dTbIZlLGWUfMJpWCdVt38x/nHMrbFkzixTd2MK2+mgm1lT2us/fkeMNk6qsrqK+uAODDR83lt8s3kM057vD4K1uKSwCFEkBYBVRdkSabz9LanmNLcxvbd7UzL2ZcMnbEbQO4DMDMpgLVke1rhikukaI0t2Y7bugbd+xmY+EGvyP4t3DT39VLlciMxmredeA06qoztLTmOG7/KRwzfxJrtrQwraGKGY01PV5zwPSGIY3/2H0m8d1zD2PxvImcetWfWPrKm3zwLXNivbbQDTSV6mwErkinqMqkaM3mueL3L/D7FRtYctkJVFekhzRuSba4VUCnAd8GZhK0A8wFngcOHr7QZKzL552tu9rZ0tzKpp1tbNzRGt7Id7Ox48Ye/Nvc1vPGXlORZmpDFdMaqjlk9nhOrA8eF7ZNa6hman1Vn9MlT+zl2/1wMTPOOGwWAEfOncCyNW/Gel0+H5QYOlcEC3sFpYyqTJrWbJ61b+5ia0s7f1ixoeMa7kH1kBLC2Ba3CujfgGOAu939cDN7J8G6wHvEzE4BvgukgR+7+9f39JxSntydXe05tu1qZ/uuLFtb2tjS3Mbm5jY272wLbvLNbWzZ2cbm5la2NAf7I+2ZHaorUh037wNnNnD8/lOZ1lAV3Njrq5ka3uTrqzKY2cj/snvoiLkTuPv5jWxpbhswCRXm/i9UAeXywWRwaTOqK1Lsbs/RtKMVCMYiFBLAtY+u4YrfvcD9X3gndVUZTv7OA3zsrXM5/9i9eXlTM0+8+iZnHTkbgLZsnsqMVhUbjeImgHZ332xmKTNLufu9ZvadPbmwmaWB/wTeBawFHjezW9W4XH7asnl2teVobsvS0paluTV8XPi3LUdza5Ztu9p7/dke3vTbcn03SDZUZ5hcV8XE2kr2nlzLkXMnMrmukom1lUyqq2JSbSVT66uY2lBNQ3Uyb+xxLZ4bLLb3xKtvcuJB0/o9ttDjJ5XqHAdQKAGkU0EVUNPOViozKZas3sLLm5qZM3EcP3pgNTtas9yxfD3T6qt4eVMz//3QK/zVW+fx1dtWcO+LTew/vZ7x4yo49XsP8tl37sMFb5/PLU+u5YlXt/Ll0w8mm89z7ZI1nLpoBtMaqnltSwstbTn2nx7ME7lu6y5mNlYHs5TmHQvjlPIRNwFsNbM64E/AdWa2kXBm0D1wFPCSu68GMLMbgDOAIU8A9zz/Bk+9thUIFswocDqfeLdvmtGnfb2mj4d4t5N5n8fFi6XruSLHFXkuB3I5py2Xpy2bpzWbDx/naM85bdlge1suT2t7jua2HC1tWdpz/QQUYQYN1RU01nT+zGysoaGm67bCz8TaSibXVTKhtpIKrVvbYdHsRjIpY9magRNAvlACCNsAggVh8mTSRmU6xc7WLFua2/jQUXP41dLXuPrelzhl4XTWbGmhMpPiN0+uY2pDFQBrtrRw07K13LeyCYAf3L+Kmoo0W1va+Y87V3LQzAb+8dfLac3mmTm+hjVbmrn+sde4celr/NuZC7ng50tpac3xnXMPY8nqzfz8kVc57dCZnLN4Nv/8f88CcNl7DmTpK1v4w4o3eO+iGRw0o4FfLX2Nmoo07zlkBi9s2MGLG7Zz5NwJNI6r5OnXtlJbmWb+lDo2N7extaWNCeMqqUgbzW058u7BamgEn7/C49HofUfMjt3BIK6BFoX/PPAwcCbQAnweOA9oBL66h9eeBbwWeb4WOLqXGD4JfBJgzpx4jWLd3b+yiWuXvBo9Z+fjLtfqdu3o3t4fdnlN9Pie54q+pvfr932NriezGMf1FWM6ZVRmUlSmU1Rm0uHjoLqgoToTPM+kqcqkqK1MM64qE/xbmaG2qtu/lRnGVaU7/q2rzOgb3hCorkhz8KxGHvzzJi49Kfg235fCtA+ZblNBBFVAada9uQuAhbMaaKyZzw/vX8U9L2xkan0V5x09lyvvXklF2jj3LXtx+zPr+affPEvajDMOn8UtT64F4L2LZnDXc2/w0Z88Rm1lmrcumMQVf3gBdzj54Gnc8/xGzv7hI0xrqGLuxHF85rpg9djj95/CHc+8zm1Pv87sCTVkUsbf/GIZZnDknAn88P5VuMOs8TW05fLc+dwbZFLGnEnjuPfFIAlNqq1kd/hlJGVBz6ntu9txD3/nlIEHv7fTmRBHoyPmThjZBADMBr4DHAAsBx4iSAi3ufuWIY2kD+5+DXANwOLFiwf1v/vVMxby1TMWDmlcIsPpg4v34ou3LOcLNz3DN89e1GdiLcw3lAqngsi5k/dCI3CKlW8E8w1NqaviQ2+Zw8Ydu7n5iXX8/Yn78f4jZnHl3StpzzkfOir4cnXD469x2qEz+Yd3789tz7xOZTrFv56xkPmTa/neH1/in049iJMPms5p33+QfafWcfV5R3Lb06/z04de5soPHsb0hmr+9fbnOHLuBD6weC8eXb2Zh1dt5sJ3zCeTMn737HoOmN7AgTMaWNW0k/Vbd/PWBZNwd55Zt435k2sZP66SzTtbaWnLMXtCDe7QtLOVCeMqqcykyOby5NypyqgBe0/1mwDc/VIAM6sEFgNvA84HrjGzre5+0B5cex3B+IKC2eE2kTHvw0fPoWlHK1fevZIDptdz4Tt6X4K70AaQDlcEc4f2XJgAKlJs350FYEp9FamU8Y2zFnH8/lM56aBpVFekOWrviWxpbmPR7EbSqbnc/sx6Lnz73kytr+Zr7zuEcZVpJtZW8vkT9+PEA6exaHYjZsbdFx9HRTq45pmHz+LMw2d1xPT1sxZ1PD56/iSOnj+p4/n7Dp/d8XjBlDoWdIx1MI6YM6Fj36S6KgqvMoNpDR29z8mkU0VPYSC9i/s+1gANBFU/jcDrBCWCPfE4sG+4tvA64Fzgw3t4TpFR46IT9uGWJ9d2tF/1JhdpA0iHdX3ZXJ50yqiOfEOeUh/U81ekU5x+6MyO7T8474iggdaMhbMaefYrJ3fsO/vIzpt1OmUcutf4jufqFTQ6DNQGcA1BX/8dwKME1T/fdvd4nZT74e5ZM/ss8AeCbqA/dfcVe3pekdHCzJg5vob123b1eUy0BFCoJcrmnUxYAiiYXFfV6+sn9bFdxoaBSgBzgCrgzwTf0tcCfX8dKZK7/xb47VCdT2S0md5YzZJVm/vc35EAzDraCVqz+Y6BYAD11RkN+JJeDdQGcIoFXUsOJqj/vwRYGK4L8Ii7/8sIxCgyZs1orOaNHa3hso89G4ILi38FbQDB4/Zcnkwq6NkFndU/It0NWJHngWcJvqn/jqAn0ALg74Y5NpExb3pjDbm8s2lna6/7s2EG6NoG4F1KAFNUzSN96DcBmNlFZnaDma0B7gdOBV4A3g9MHIH4RMa0mY1B75f123b3ur/Q770wEhgKJYBgMjhQCUD6NlAbwDyC+f//3t3XD384IhI1PUwAG7btgkgvnILC7Bpp66wCaiu0AYT1/koA0peB2gAuHqlARKSnwlTUfZUAuvYCCjJAWy7fMRAMlACkb+rMK1LGJoyroDKTYkOMBFBoJO5sAwj+vPvqAiqiBCBSxsyMGY3VfZcAOgaC0TEOoLMXkKqApH9KACJlbnpDdYwSQKpjMsC2XJ5UtApIJQDpgxKASJmb0VjN+u29jwbuMhCsSy8g49h9JvOho+aw37T6EYtVkkUJQKTMTW+s4Y1trR0zf0Z1LggTVANB52RwM8fX8O/vP0Tz9kif9MkQKXMzGqtpy+XZ0tLWY18+siRkoQooF64HIDIQJQCRMtc5FqBnO0C2l26gAJm0EoAMTAlApMzNDMcCrH2zZztAWzYYCVaVSZPqtvKbyECUAETK3JxJ4wB4dXPPZbg7EkBFqksJQFVAEocSgEiZa6ypYFJtJS9v6pkAWrM5AKoyqS7LRqZT+tOWgelTIpIAe0+uZXWvCSAoAVRmUl2qgNQGIHEoAYgkwN6Ta3mllwTQtQ0gWgJQApCBKQGIJMDeU2rZuKOVna3ZLtu7VAFFG4HVBiAxKAGIJMD8ybUAPUoBre3RKiCVAKQ4SgAiCbD35DqAHu0Ardk8KYNM93EASgASgxKASALMnTQOM3i5qWsCaMvlqcqkMTOiHX9SSgASgxKASAJUV6SZ2VjDy5t2dtne2p7rmOvHVAKQIikBiCTE/Cm1PcYCtGbzHdM+p9UGIEVSAhBJiMJYAPfOWUHbsnmqKoI/YzUCS7GUAEQSYvaEGnbsznbpChqUAIKVv7oMBFMCkBiUAEQSoq6qAoDm1lzHttZsjsp0zzYATQUhcehTIpIQtVXBN/0eJYCwCiha7aMSgMShBCCSEHVVGQCae1QBFdoAOo9VN1CJQwlAJCFq+0gAlWEbgLqBSrGUAEQSolAC6FIF1J7rtQSgXkAShxKASELU9pIAgpHAPdsAlAAkDiUAkYQoNAJ3qQJqj3YDVQKQ4igBiCREZxVQtBtoPjIVROexagOQOJQARBKipiJY+D1aAmjLRtsAOm/66gUkcSgBiCSEmVFbldE4ABkySgAiCVJXlekoAbh7kADS6gUkg1OSBGBm3zSzF8zsGTO7xczGlyIOkaSprcrQ3BYkgPZcMClcVUXPcQBaElLiKFUJ4C5gobsvAlYCl5UoDpFECaqAgkbg6HrA0LUNIJNWApCBlSQBuPud7l6oyFwCzC5FHCJJU1eV7qgCas12rgcM3dcDUO2uDKwcPiWfAH7X104z+6SZLTWzpU1NTSMYlkj5qa3MsHN31wRQpW6gMkiZ4Tqxmd0NTO9l15fc/TfhMV8CssB1fZ3H3a8BrgFYvHix93WcyFhQF+kF1NaRAMKBYJGbfkptABLDsCUAdz+xv/1m9nHgVOAEjy5xJCJ9ijYC92wD6DxObQASx7AlgP6Y2SnAF4Dj3L2lFDGIJFFtpBtoa3t/bQBKADKwUrUBfB+oB+4ys6fM7IclikMkUeqrM7TnnNZsjrZc1yogdQOVYpWkBODu+5TiuiJJV1tZmBAu11EC6FwUvvM4lQAkjnLoBSQiMUUXhSm0AVSmNQ5ABkcJQCRBoovCdHQDLZQAUqoCkuIoAYgkSLQE0KMbqKqApEhKACIJUtulBBBWAfU2FYRGAksM+pSIJEivVUC9rgcw8rFJ8uhjIpIg0WUh27ongMhfs0oAEoc+JSIJEl0WsrVHG4AGgklxSjIOQEQGJ9oInA0HglWEXT67tgEoAcjAlABEEqQinaIqkwqmg7Cg+qcwAjh6zzwetZoAAAgYSURBVNeawBKHEoBIwhRmBC0kgwIzw0xjACQ+tQGIJExhQrjWbJ7KsP6/IGWm+n+JTSUAkYQpLAuZSlmXEgAE1UBKABKXSgAiCVNXlWZnaztt2XzHNBAFphKAFEEJQCRhptZXs2Hbblqz+Y4uoAVpM/UAktiUAEQSZsGUWtZsaWHH7vaOaSAKVAUkxVACEEmYBVPryDu8uGFHL20AqgKS+JQARBJmwZQ6AN5sae+RAMw0DYTEp0+KSMLMn1Lb8bhHCSBlmghOYtNHRSRhxlVmmDW+BqCPRmD9WUs8+qSIJNCCqUE1UM8qILUBSHxKACIJtE/YDtBbLyB1A5W4lABEEmjB1KAdoLdeQCnNBSQxKQGIJFChJ1BVRbc2gJSRSSsBSDxKACIJtE/YBlCZ7tkNVG0AEpcSgEgCTaqt5LPv3IdTFk7vsj1lpumgJTbNBiqSQGbGpSfv32O7poKQYqgEIDKKpNQGIEVQAhAZRdQLSIqhBCAyimgcgBRDCUBkFAlmA9WftcSjRmCRUeQz79yH8TUVpQ5DEkIJQGQUOf3QmaUOQRJEZUURkTFKCUBEZIxSAhARGaOUAERExqiSJgAzu8TM3MwmlzIOEZGxqGQJwMz2Ak4C1pQqBhGRsayUJYArgS8AXsIYRETGrJIkADM7A1jn7k/HOPaTZrbUzJY2NTWNQHQiImODuQ/PF3AzuxuY3suuLwFfBE5y921m9gqw2N03xThnE/DqIEOaDAx4jTKieIeX4h1eind4FRvvXHef0n3jsCWAvpjZIcA9QEu4aTbwOnCUu28YxusudffFw3X+oaZ4h5fiHV6Kd3gNVbwjPhWEuy8HphaeF1MCEBGRoaNxACIiY1TJJ4Nz93kjdKlrRug6Q0XxDi/FO7wU7/AaknhHvA1ARETKg6qARETGKCUAEZExatQlADM7xcxeNLOXzOwfe9lfZWa/DPc/ambzRj7KLvEMFO87zOwJM8ua2dmliLFbPAPFe7GZPWdmz5jZPWY2txRxRuIZKN5PmdlyM3vKzB40s4NKEWcknn7jjRx3VjiPVkm7LsZ4fz9uZk3h+/uUmV1Qijgj8Qz4/prZOeFneIWZ/e9Ix9gtloHe3ysj7+1KM9ta1AXcfdT8AGlgFTAfqASeBg7qdsxngB+Gj88Fflnm8c4DFgE/B85OwPv7TmBc+PjTCXh/GyKPTwd+X87xhsfVAw8ASwi6UJdtvMDHge+XKsZBxLsv8CQwIXw+tZzj7Xb854CfFnON0VYCOAp4yd1Xu3sbcANwRrdjzgD+J3x8E3CCmdkIxhg1YLzu/oq7PwPkSxFgN3HivdfdC4P8lhAM9CuVOPFujzytpbRzU8X5/AL8K/ANYPdIBteLuPGWizjxXgj8p7u/CeDuG0c4xqhi398PAdcXc4HRlgBmAa9Fnq8Nt/V6jLtngW3ApBGJrqc48ZaTYuP9a+B3wxpR/2LFa2Z/a2argCuAi0Yott4MGK+ZHQHs5e53jGRgfYj7eTgrrBK8KZwFuFTixLsfsJ+ZPWRmS8zslBGLrqfYf29hVevewB+LucBoSwBSJszsI8Bi4JuljmUg7v6f7r4A+Afgn0odT1/MLAV8G7ik1LEU4TZgnrsvAu6is/RdrjIE1UDHE3yj/pGZjS9pRPGcC9zk7rliXjTaEsA6IPoNY3a4rddjzCwDNAKbRyS6nuLEW05ixWtmJxJM+ne6u7eOUGy9Kfb9vQE4c1gj6t9A8dYDC4H7wilUjgFuLWFD8IDvr7tvjnwGfgwcOUKx9SbO52EtcKu7t7v7y8BKgoRQCsV8fs+lyOofYNQ1AmeA1QRFoUKjycHdjvlbujYC/6qc440c+zNK3wgc5/09nKDhat+EfB72jTw+DVhazvF2O/4+StsIHOf9nRF5/D5gSZnHewrwP+HjyQRVMJPKNd7wuAOAVwgH9hZ1jVL9Zwzjm/Yegqy9CvhSuO2rBN9GAaqBG4GXgMeA+WUe71sIvpU0E5RUVpR5vHcDbwBPhT+3lnm83wVWhLHe298Ntxzi7XZsSRNAzPf338P39+nw/T2gzOM1gmq254DlwLnlHG/4/MvA1wdzfk0FISIyRo22NgAREYlJCUBEZIxSAhARGaOUAERExiglABGRMUoJQMYsM5sUmUlxg5mtizx/eAivc6aZXd7P/kPM7GdDdT2RuNQNVAQwsy8DO939W8Nw7ocJ+m1v6ueYu4FPuPuaob6+SF9UAhDphZntDP893szuN7PfmNlqM/u6mZ1nZo+F6wgsCI+bYma/NrPHw59jw+37Aa2Fm7+ZfcDMnjWzp83sgcglbyMYmS4yYpQARAZ2KPAp4EDgo8B+7n4Uwdw2nwuP+S5wpbu/BTgr3AdwLPBE5FyXAye7+6EE6w8ULAXePmy/gUgvMqUOQCQBHnf39QDhtNF3htuXEyyAA3AicFBkaYkGM6sDZgBNkXM9BPzMzH4F3BzZvhGYOTzhi/ROCUBkYNEZTfOR53k6/4ZSwDHu3mWRFjPbRTDjLADu/ikzOxp4L7DMzI50980Ec1TtGqb4RXqlKiCRoXEnndVBmNlh4cPngX0i2xe4+6PufjlByaAw3e9+wLMjFKsIoAQgMlQuAhaHK189R9BmAMHavYdHlh39Zth4/CzwMMEsmRBUJZXDKl8yhqgbqMgwM7PvAre5+9197K8C7gf+woNlSkVGhEoAIsPva8C4fvbPAf5RN38ZaSoBiIiMUSoBiIiMUUoAIiJjlBKAiMgYpQQgIjJGKQGIiIxR/x8v/+oIbQFBPQAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] @@ -69,7 +69,7 @@ "source": [ "\n", "# Plot the waveform of L1.\n", - "plt.title('L1')\n", + "plt.title('Original signal')\n", "plt.xlabel('Time(s)')\n", "plt.ylabel('Waveform')\n", "plt.plot(times_L1,data_L1)" @@ -112,9 +112,9 @@ "def whiten_signal(signal, asd):\n", " \n", " signal_rfft = np.fft.rfft(signal)\n", - " whiten_signal_rfft = signal_rfft / asd\n", - " whiten_signal = np.fft.irfft(whiten_signal_rfft, n=Nt)\n", - " return whiten_signal" + " whitened_signal_rfft = signal_rfft / asd\n", + " whitened_signal = np.fft.irfft(whitened_signal_rfft, n=Nt)\n", + " return whitened_signal" ] }, { @@ -159,24 +159,63 @@ "f_min = 20.\n", "f_max = 128. \n", "plt.axis([f_min, f_max, 1e-24, 1e-19])\n", - "plt.show()\n", - "\n", - "# A small problem here: freq_rfft contanins some values which is out of freq.\n", - "# I am not sure if this has a big impact on the calculation of the whitening process." + "plt.show()\n" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 7, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZAAAAEKCAYAAAA8QgPpAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3deXjU5bn/8fc9SVhktS5UQQWKS2UpwSjyEy2IC2pRcKdabfGg1FLUnqNCpVVPVVB6EKla1KOi1SJqNUJRaVVyAIsVBQRBqcqigBVF2UGy3L8/ZjIkYZKZzJKZzHxe15WLmee7zJ3HmDvP+jV3R0REpL4C6Q5AREQaJyUQERGJixKIiIjERQlERETiogQiIiJxUQIREZG45Kc7gFQys0HAoFatWg0/6qijANixYwctWrRIb2CNnOowcarDxKkOExetDt99992v3P2g2o5bLqwDKSoq8nfeeQeAkpIS+vXrl96AGjnVYeJUh4lTHSYuWh2a2bvuXlTbcXVhiYhIXJRAREQkLlmdQMxskJk9vGXLlnSHIiKSdbJ6EN3dZwIzi4qKhqc7Fkm+0tJS1q1bx+7du9MdSr21adOGDz74IN1hJE2zZs3o0KEDBQUF6Q5FGlBWJxBJTPHi9UyYvZINm3dxaNvm3Hjm0QwubJ/usMLWrVtHq1at6NixI2aW7nDqZdu2bbRq1SrdYSSFu7Np0ybWrVtHp06d0h2ONCAlEImoePF6xrywjF2l5QCs37yLMS8sA6BtOgOrYvfu3Y0yeWQbM+OAAw7gyy+/THco0sA0BiIRTZi9Mpw8Ku0qLWfC7JVpiigyJY/MoP8OuSmrE4i7z3T3q9u0aZPuUBqdDZt31as8V7Vs2XKfsrlz59KrVy/y8/N5/vnn67y+uLgYM+PDDz8Ml1VUVDBq1Ci6detG9+7dOf7441m9ejUAHTt2pHv37nTv3p1jjz2WsWPHNsoxIMkOWZ1AJH6Htm1er/LGoHjxek4a/wadRs/ipPFvULx4fUo+5/DDD2fq1Kn8+Mc/jnrutGnT6Nu3L9OmTQuXTZ8+nQ0bNrB06VKWLVvGiy++SNu2ezsO58yZw7Jly3j77bdZtWoV11xzTUq+D5FolEAkohvPPJrmBXnVypoX5HHjmUenKaLEVI7prN+8C2fvmE4qkkjHjh3p0aMHgUDd/3tt376d+fPn8+ijj/LMM8+Eyz///HMOOeSQ8PUdOnRg//333+f6li1bMmXKFIqLi/n666+T+02IxECD6BJR5WyrG59/j9Jy56CWTbnlnO8zuLA9JSUfpTm6fd0+czkrNmyt9fjiTzezp7yiWtmu0nJuen4p097+NOI1xx7amlsHdU1qnFW99NJLDBw4kKOOOooDDjiAd999l+OOO46LL76Yvn37Mm/ePAYMGMDll19OYWFhxHu0bt2aTp068dFHH9G7d++UxSoSiVogUqvBhe055rutAXj0p0UZNYW3vmomj2jlDWHatGlceumlAFx66aXhbqwOHTqwcuVKxo0bRyAQYMCAAbz++uu13icX9rOTzJTVLZDK3Xi7dOmS7lAavUz/HRWtpXDS+DdYH2ECQPu2zZl+TZ9UhVWrr7/+mjfeeINly5ZhZpSXl2NmTJgwATOjadOmnHXWWZx11lm0a9eO4uJiBgwYsM99tm3bxpo1a6jcbVqkIWV1C0SzsKRSpo3pvPTSS/zkJz9h7dq1rFmzhs8++4xOnToxb948Fi1axIYNG4DgjKylS5dyxBFH7HOP7du3c+211zJ48OCIYyQiqZbVCUQSly3T+wcXtmfc+d1p37Y5RrDlMe787gl3y+3cuZMOHTqEvyZOnMjChQvp0KEDzz33HNdccw1du+7bOnr++ecZMmRItbILLriAadOmsXHjRgYNGkS3bt3o0aMH+fn5jBw5Mnxe//796datGyeccAKHH344Dz30UELfg0i8sroLSxKX6V1X9TG4sH3Sx3EqKiKPoaxbt67O62bNmrXPViajRo0Kvx44cGDE69asWVO/AEVSSC0QiUm2tEREJHmUQEREJC5KICIiEhclEBERiYsSiNTJyaJRdBFJqqxOINrOPXkMjaKLSHVZnUC0kDB51BKJLNJ27jVNmjSJnTt3pjyWqVOnVlsvEklJSQn/+Mc/wu+nTJnCk08+merQJEtldQKRxGVNy2P+JFg9t3rZ6rnB8hSLJ4GUl5dHPykONRPIiBEjuOKKK1LyWZL9lEAkN7TvBc/9dG8SWT03+L59r6TcvqSkhH79+nHhhRdyzDHHcNlll+HuTJ48mQ0bNtC/f3/69+8PwN/+9jcGDBhAr169uOiii9i+fTsQ3Ab+5ptvplevXjz33HP069eP6667jp49e9KtWzfefvttILiP1uDBg+nRowcnnngiS5cu3SeemTNn0rt3bwoLCznttNP44osvWLNmDVOmTOHee++lZ8+ezJs3j9tuu43f//73ACxZsoQTTzyRHj16MGTIEL755hsA+vXrx80338wJJ5zAUUcdxbx585JSZ9L4KYFInbKm66rTKXDR1GDSeOPO4L8XTQ2WJ8nixYuZNGkSK1asYNWqVbz55puMGjWKQw89lDlz5jBnzhy++uor7rjjDmbMmMGiRYsoKipi4sSJ4XsccMABLFq0KLxL786dO1myZAkPPvggw4YNA+DWW2+lsLCQpUuXctddd0VsQfTt25e33nqLxYsXc+mll3LPPffQsWNHRowYwQ033MCSJUs4+eSTq11zxRVXcPfdd7N06VK6d+/O7bffHj5WVlbG22+/zaRJk6qVS27TViYSk6zoyup0ChRdBXPvgVNuSmryADjhhBPo0KEDAD179mTNmjX07du32jlvvfUWK1as4IwzziAQCLBnzx769Nm7G/All1xS7fyhQ4cCcMopp7B161Y2b97M/Pnz+ctf/gLAqaeeyqZNm9i6tfqzUNatW8cll1zC559/zp49e+jUqVOdsW/ZsoXNmzfzwx/+EIArr7ySiy66KHz8/PPPB+C4447TdioSphaI5I7Vc+GdR4PJ451H9x0TSVDTpk3Dr/Py8igrK9vnHHfn9NNP580332TJkiWsWLGCRx99NHy8RYsW1c63GnvI1Hxfm1/+8peMHDmSZcuW8dBDDyX83PTK762270tykxKI5IbKMY+LpsKpt+ztzkpyEomkVatWbNu2DYATTzyRN998k08++QSAHTt28K9//avWa6dPnw7A/PnzadOmDW3atOHkk0/m6aefBoJjLwceeCCtW7eudt2WLVto3z64ceQTTzwRMZaq2rRpw/777x8e3/jTn/4Ubo2I1EZdWJIb1i+qPuZROSayflHSu7Jquvrqqxk4cGB4LGTq1KkMGzYs/Jf8HXfcUesDoZo1a0ZhYSGlpaU89thjANx2220MGzaMHj16sN9++1VLEJVuu+02LrroIvbff39OPfVUVq9eDcCgQYO48MILeemll/jDH/5Q7ZonnniCESNGsHPnTjp37szjjz+ezGqQbOTuWf913HHHeaU5c+a4xO7s++b6ETf/1Zd+tjlclil1uGLFinSHELetW7dGPeeHP/yhL1y4sAGiSY6G/u+RKT+HjVm0OgTe8Tp+t6oLS0RE4qIuLJEMVVJSku4QROqU8S0QMxtsZo+Y2XQzOyNU1sLMngiVX5buGLNV8eL1/OuL4IDrsCcWUrx4fZojEpFMktIEYmaPmdlGM3u/RvlAM1tpZh+b2ei67uHuxe4+HBgBVE6SPx94PlR+bkqCz3HFi9cz5oVllJYHFxJ+ue1bxrywLOOSiGfTM3cbMf13yE2pboFMBao93NnM8oAHgLOAY4GhZnasmXU3s7/W+Dq4yqVjQ9cBdAA+C71OzaZBOW7C7JXsKq1etbtKy5kwe2WaItpXs2bN2LRpk355pZm7s2nTJpo1a5buUKSBpXQMxN3nmlnHGsUnAB+7+yoAM3sGOM/dxwE/qnkPC66cGg+84u6LQsXrCCaRJdSSBM3sauBqgHbt2oX7k7dv366+5Ris37yr1vLt2z0j6tDMaNGiBZ999ln0kzOMu8e8KLAxKC8vZ8eOHaxdu7bBPlP/Lycu0TpMxyB6e/a2HiCYDHrXcf4vgdOANmbWxd2nAC8A95vZOcDMSBe5+8PAwwBFRUXer18/YO+md1K39m+9ETGJtG/bnJYtA6rDBOnnMHGqw8QlWocZPwvL3ScDk2uU7QB+Fu1aMxsEDOrSpUuKosteN555NGNeWFatG6t5QR43nnk0bPkojZGJSKZIxyys9cBhVd53CJUlneuBUnEbXNieced3pyAv2M1yUMumjDu/O4ML26c5MhHJFOlIIAuBI82sk5k1AS4FZqQhDolicGF7jmrXCoBHrixS8hCRalI9jXcasAA42szWmdlV7l4GjARmAx8Az7r78hR9vp6JniBNcBKR2qQ0gbj7UHc/xN0L3L2Duz8aKn/Z3Y9y9++5+50p/Hx1YSWg6kLC4U+8k3FrQEQkvTJ+EF3So3IhYVlFaCHh9uBCQoC26QxMRDJGxm9lkgh1YcWvMSwkFJH0yuoEoi6s+G2oZSFhbeUiknuyOoFI/A5t27xe5SKSe7I6gagLK343nnk0zQvyqpWFFxKKiJDlCURdWPGrXEiYHwguJDywZRMtJBSRarI6gUhiqi4kfPCy45Q8RKQaJRCpVdV1INc+/a7WgYhINVm9DkSbKcav5jqQr7bv0ToQEakmq1sgGgOJn9aBiEg0WZ1AJH5aByIi0SiBSERaByIi0SiBSESR1oEY0P+Yg9ITkIhknKxOIFpIGL/Bhe254Ljq03Yd+Mu76/nHhtL0BCUiGSWrE4gG0RMz58Mv9ynbVVrOX/6lBCIiWZ5AJDG1DZhv2q2nTImIEojUoe1+BRHLD2hmDRyJiGQiJRCJqHjxerbvLtunvCDPuOCoyIlFRHKLEohENGH2Skor9u2qatEkn/93qBKIiGR5AtEsrPjVNv6xZZcG0EUkKKsTiGZhxa+2BYNtmqv1ISJBWZ1AJH43nnk0BYF9B8t37CnTOhARAZRApBaDC9vTstm+mzWXlrvWgYgIEON27mZWBJwMHArsAt4H/u7u36QwNkmzb3ZGThRaByIiEKUFYmY/M7NFwBigObAS2Aj0BV4zsyfM7PDUhykNrXjxempb7aF1ICIC0Vsg+wEnuXvEKTlm1hM4Evg02YFJek2YvZJI7QwDrQMRESBKAnH3BwDM7DB3/6zqMTP7rrsvSWVwkj61TeN10DoQEQFiH0RfbWbTzGy/KmUvpyKgZEpoHcj8SbB6bvWy1XOD5ckQ7/1jvS7aeXUdnz+JH7X6uNqhPoHlXJM3EwPNwhIRIPYEsgyYB8w3s++FyjK+IzyhdSDte8FzP937S3b13OD79r2SE1y894/1umjn1XW8fS/GV0ykT2A5EEwe9xdMZql3Dm7prllYIkKMs7AAd/cHzew9YKaZ3QwRu8izR6dT4KKpwV+qRVfBO48G33c6Jb33j/W6aOdFOT468CvuL5jAU+WncXnea4wsHcWCiq6AZmGJSFCsLRADcPc3gQHATcAxqQoqY3Q6JfjLde49wX+TlTwSvX+s10U7r47jf93WhafKT+O6/Bd5qvy0cPIAzcISkaBYE8jZlS/c/XOgPzAwJRFlktVzg3+Zn3JT8N+aYwbpun+s10U7r47jP2r1MZfnvcZ9ZUO4PO+1cHcWaBaWiATV2YVlZr+q8jrSKUn+jZpBKscEKrt1Op1c/X267h/rddHOq+s48Hu7l6vKrmN++bG8VXEs9xdMZmTpKN6mW+Lfu4hkhWgtkFZVvv6rxvtWqQ0tzdYvijxmsH5Reu8f63XRzqvr+PpFNB36JMub/gCABRVdGVk6ih62ivIKbWUiIkHR1oHcXvnazAZXfZ/1+l6/b1mnU5I3DhLv/WO9Ltp5Mdxn885Z4dcLKrqyAA2ii8he9dlMUb81ckxtW7prEF1EQLvxSh36H3NQxPIeB+nHRkSiD6IvI9jyMOB7Zra08hDBtSE9UhyfpNGcD7+MWL70y4oGjkREMlG0hYQ/apAoJCPVth+WxkBEBKJ0Ybn7WndfC1xb+bpqWUMEaGaDzewRM5tuZmeEyjqb2aNm9nxDxJCrant8bYtY9y8QkawWa2f26RHKzop2kZk9ZmYbzez9GuUDzWylmX1sZqPruoe7F7v7cGAEcEmobJW7XxVj7BKnyEt/ai8XkdwSbQzk5wRbGlXHPyC4BuTNGO4/FbgfeLLKPfOABwgmpXXAQjObAeQB42pcP8zdN4Zejw1dJw1kcy1PJNyuZSAiQvQxkD8DrxD8xV61pbDN3b+OdnN3n2tmHWsUnwB87O6rAMzsGeA8dx9HhDEXCy6BHw+84u4xr+Izs6uBqwHatWtHSUkJANu3bw+/lrrtlw87yiKU57nqMEH6OUyc6jBxidZhtIWEW8xsO1AYGvdIhvZA1YdTrQN613H+L4HTgDZm1sXdp5jZAcCdQKGZjQkln5qxPww8DFBUVOT9+vUDoKSkhMrXUrcmc//GjrJ9mxuBgKkOE6Sfw8SpDhOXaB1GHQ519/LQeMXh7t7gj65198nA5BplmwiOiUgKqQtLROoS63ya/YHlZvY2sKOy0N3PjeMz1wOHVXnfIVSWdGY2CBjUpUuXVNw+67VpXsDmXftmC83CEhGIPYH8JomfuRA40sw6EUwclwI/TuL9w9x9JjCzqKhoeCrun+00C0tE6hJtFtZs4FWCA9gf1vfmZjYN6AccaGbrgFvd/VEzGwnMJjjz6jF3X17HbeKmFkhivlEXlojUIVoL5EqCD466zcyOAv5JMKG85u476rwScPehtZS/DLxcz1jrTS2QxOSZUe77rjrXTlgiAtFnYf2b4FqOqWYWIDhb6izgJjPbBfzN3e9JeZSSFpGSB4B2whIRqMcfk+5e4e4L3P237n4SwbGLlAx+J4uZDTKzh7ds2ZLuUBqlvFoGO9QCERGIcRDdzA4ChgMdq17j7sNSE1ZyqAsrMWqBiEhdYp2F9RIwD3gNKE9dOJJJNAYiInWJNYHs5+43pzSSFNAsrMSoBSIidYn1j8m/mtnZKY0kBdx9prtf3aZNm3SH0ihpDERE6hLr74LrCCaRXWa21cy2mdnWVAYm6acWiIjUJaYE4u6t3D3g7s3dvXXofetUByfpVVsLBKB4cUZPwBORBlBnAjGzY0L/9or01TAhxk/TeBNTWwsEYMLslQ0YiYhkomiD6L8i+EyN/4lwzIFTkx5REmkab2Lat23O+lqei15buYjkjmgr0a8O/du/YcKRTHLjmUdz/fQlEY9pP0URiXlCjZl1M7OLzeyKyq9UBibpN7iwfa3Hau/cEpFcEetK9FsJ7qp7LMFNEM8C5lPlWeciIpJbYm2BXAgMAP7t7j8DfgBk/OIKDaKLiKROrAlkl7tXAGVm1hrYSPWnCmYkLSQUEUmdWLcyecfM2gKPAO8C24EFKYtKREQyXtQEYmYGjHP3zcAUM3sVaO3uS1MenYiIZKyoCcTd3cxeBrqH3q9JdVAiIpL5Yh0DWWRmx6c0EhERaVRiHQPpDVxmZmuBHQTXkbm790hZZEmg7dxFRFIn1gRyZkqjSBFtZSIikjqxdmHd4e5rq34Bd6QyMMl82pFXJLfFmkC6Vn1jZnnAcckPRxqTMS9oIp5ILou2nfsYM9sG9Ag9SGpr6P1Ggs9Jlyy3/34FtR7bVapHS4nksjoTiLuPc/dWwITQg6QqHyZ1gLuPaaAYJY1uHdQ1+kkikpPq80z0FgBmdrmZTTSzI1IYl2SIunbkFZHcFmsC+SOw08x+APwn8AnaiVdEJKfFmkDK3N2B84D73f0BoFXqwhIRkUwXawLZZmZjgMuBWWYWAGofXc0Q2s5dRCR1Yk0glwDfAle5+7+BDsCElEWVJNrOPfW0FkQkd0WbxmsA7v5vd5/o7vNC7z919yerniO5SWtBRHJXtBbIHDP7pZkdXrXQzJqY2alm9gRwZerCk0xxTd5M+gSWVyvrE1jOFRVaDiSSq6IlkIFAOTDNzDaY2QozWwV8BAwFJrn71BTHKGnWokkeS70z9xdMDieRPoHl3F8wmaXeOc3RiUi61LmZorvvBh4EHjSzAuBAgo+33dwQwUlmuHNId66fXs7I0lHcXzCZp8pP4/K81xhZOooFFVpoKJKrYh1Ex91L3f1zJY/cU7mYcEFFV54qP43r8l/kqfLTlDxEclzMCUSkT2A5l+e9xn1lQ7g877Vwd9ZljyxIc2Qikg6xPg9EclzlmEdlt9VbFceG37/5iVoiIrkophaImbU1s+NDX1pUkYN62KpqYx4LKroysnQUPWxVmiMTkXSpswViZk2Bh4DBwGqCj7I9wsxeBEa4+57UhyiZYMfxv2DBW59WK1tQ0ZUFqPUhkquitUBuIbhlyWHuXujuPYHDCSae36Q6uEpmNtjMHjGz6WZ2Rm1lkjp3DO5e53GNg4jknmgJ5HxguLtvqywIvb4WGBLLB5jZY2a20czer1E+0MxWmtnHZja6rnu4e7G7DwdGENxWJWKZpM+bn3yd7hBEpIFFG0SvcPedNQvdfbuZeYyfMRW4nyrbv4ceifsAcDqwDlhoZjOAPGBcjeuHufvG0OuxoeuqilQmIiIpFq0F4ma2v5l9p+YXENPzTN19LlDzz9MTgI/dfVVoHOUZ4Dx3X+buP6rxtdGC7gZecfdFENyDq2aZpNaRB7eo8/jpE0saJhARyQjRWiBtgHcJDp7XFGsLJJL2wGdV3q8Detdx/i+B04A2ZtbF3afUUhZmZlcDVwO0a9eOkpISALZv3x5+LfVzSy/46au1H/9o4w7VbYz0c5g41WHiEq3DaFuZdIz7zknk7pOBydHKahx/GHgYoKioyPv16wdASUkJla8lDq/OqvOw6jY2+jlMnOowcYnWYbTt3I+ouu7DzPqb2X1mdoOZNYn7U2E9cFiV9x1CZUmlB0olX36g7t37e9/59waKRETSLdoYyLNACwAz6wk8B3wK9CS4yWK8FgJHmlmnUCK6FJiRwP0i0gOlku/3F/0g/DrSFu+ddyyG+ZMaOiwRSYNoCaS5u28Ivb4ceMzd/wf4GcGB8KjMbBqwADjazNaZ2VXuXgaMBGYDHwDPuvvyuu4TD7VAkq9yY0Wg1i3er52TyPCYiDQW0QbRq/ZXnAqMAXD3ilgfROjuQ2spfxl4OaabxMndZwIzi4qKhqfyc3JN2yawec/e7Uz22eL926PSHaKINIBoLZA3zOxZM7sP2B94A8DMDgG0jUmOmnTq3um8tW3xrim9ItkvWgK5HngBWAP0dffSUPl3gV+nMK6kUBdW6lS2P2vb4v2jjTvSF5yINIg6E4gHPePu97p71VlSLYHzUhta4jSInjr3XtKz2hbv95ZdFO7OqkwimpElkt1ifqCUmRWa2QQzWwP8N8HBb8lRgwvbR93i/Ytt6uUUyWbR1oEcZWa3mtmHwB8ITuE1d+/v7vc3SISSsXYc/4t9Hmu7oKIrD5UPCr/vMqbuhYci0nhFa4F8SHD21Y/cva+7/wEoT31YyaExkNSqucV7pHUhx9tyXn2ozs2WRaSRimU798+BOaFnbwwg8r5YGUljIKl3+YmHh1/Xti7kibXfSVd4IpJC0QbRi939UuAYYA7BWVkHm9kf9RAngeqtkKrrQm7If67aM9Q1rVck+8Q0iO7uO9z9z+4+iOC+VYuBm1MaWRKoC6thTLqkZ/h1betCNK1XJPvEPAurkrt/4+4Pu/uAVASUTOrCahiDC9vTLC/Ys1nbuhCAjqM1oC6STeqdQEQi+fDOs6OuCwE9O10kmyiBSNKc2mpdtXUhPWwVD5adG14XAlCxep526xXJEkogkjTDx0yuti5kqXfm2vwZLPXOwN5ZWUNf1gJDkWwQbTfeRs3MBgGDunTpku5Qcsaa8eeExzpq3a23oiu97/w7/7zl9DRHK5KdjrnlZXaXV3+sQuumeSy9fWBSPyerWyAaRE+Pdq32PqyytllZX2zbQ/HipD+EUiSnXfbIAjqOnrVP8gDY+m05PW59Namfl9UJRNKjasuirllZ109fko7wRLJSx9GzePOTr+s8Z+u3yd1IRAlEUmJSjd16d3ozHiw7d5+V6uNvGZHmSEUat2NueTltU+SVQCQlBhe254SC1eExj8oB9cpZWZXJ5T3vrKm9InGoq7uqoSiBSMrccPsf99nq/dr8Gexnu6ttcxKt2S0i1cXSXRVJ66Z5SY0jqxOItjJJvzXjzwm/rm1AHbRKXSQWiXZXaRZWPWgWVmY46XvB3XgrB9Tnl3flp3mvVhtQ13iISO0S7a5q16pJtT/mkiWrE4hkhqeH96k2oP5A+WAMY0rBRPoElms8RKQO8XZXVVoz/pyUrbnK6oWEkjmmnd2EoS/v3ebkmtIbeKjgXn6RV8z3A5/u3QJF4yEiQPBpnmUJjI9PuqQngwvbJy+gCNQCkYbR93ouueiy8NsFFV15vPxM+uYtp6TiB9XGQ4b++h7tlyU56/SJJXQcHX/yaN00jzXjz0l58gAlEGlAgwvbh1epV46H/KW8L0MCb3JV3qxwufbLklzVcfSshJ6ds2b8OUkfKK+LurCkQf3zltMZ+ut7qk3jXZF3BLfk/5ljA2vpF3gvXH7MLS/z4Z1npztkkZRrDN1VkSiBSIOrOR7yaPk5HBtYywV587mvbEi4fHe5M7Z4WbXH5opkk9MnliTU4mjXqklaNyVVApGG1/d6pvXdu/ajT2A5/QLvhaf3vlVxbDiJfPL2K3Dg69D3+nRGLJJ0ia59SsW03PrK6jEQLSTMbCd97zsxTe/VeIhkk06jZyWUPCZd0jMjkgdkeQLRQsLM9vTwPvQMrAqPeSyo6Mo1pTdgGL/IK642TqKV6tLY9bj1VTqOnkW8Qx2ViwHTMdZRm6xOIJL5br5jSrUpvAsqurK0ohN985ZX2+5EK9WlsRpbvIyOo2fFvZW6kdrFgIlQApG0q9oc7xNYTvfAKnZ6k/B2J1VXqve+8+9pjFSkfjqOnsVTb30a9/Vrxp/D6gzpropECUQywprx54QTxYjSX3FV6Y00pYzHC+7moYJ7w11ZnXcs5tWHRqc7XJE6Jbrp4ZEHt7WWVHsAAA/kSURBVMiYcY66aBaWZIxzD/w3I7/aO713VkVvLsibz0flHVhQ0XXvgPvaUTTcUimR2CU6LTff4ONxmZ84KqkFIhlj6K8mstD3jnn0C7zHX8r70jWwlv8peFCD6pKxihevT8oq8saUPEAJRDLMx+POqTa19z9Lr2W5H8EFefOr7ZnVJ7CccRpUlwzQcfQsrp++JO7rG0t3VSTqwpKMU3Wlep/Acg63L/jW8zkn8E+eD/wQIJxgOo6e1Wj/55PGLdHtR5rlWaPfqkcJRDJPaKV65Z5ZI0p/BcDjBffweMHd7KEJ15TeEE4w426ZyZg7p6Q5aMkFxYvXJ9TaqJQtf/SoC0sy1i09d1dbZDirojfNrIy1FQdXG1Rf6p01JiIpd/rEkoSTRyatIk+GjE8gZjbYzB4xs+lmdkao7PtmNsXMnjezn6c7RkmNbhf/lkCnk4Gag+preLLgrmqD6lpoKKmSjAHyTFxFngwpTSBm9piZbTSz92uUDzSzlWb2sZnVOanf3YvdfTgwArgkVPaBu48ALgZOSlX8kn5PD+/DgKYfVBtUn1/RjVPy3md5xRHVWiLveWe6jFFLRJIn0QHyfMvcVeTJkOoxkKnA/cCTlQVmlgc8AJwOrAMWmtkMIA8YV+P6Ye6+MfR6bOi6yvucC/wc+FOqgpfM8Ojp+dUG1bsG1jK3vBsnB97nyYK76BpYW60lMuE3M7nxdxoTkfh1SmDPqkrZ1FVVm5S2QNx9LlDzIdcnAB+7+yp33wM8A5zn7svc/Uc1vjZa0N3AK+6+qMq9Z7j7WcBlSHbrez3T7rqp2vTeK0p/zbxaWiKLyjtzzC0vpztqaYS6jJmV0IaHAJefeHhOJA8Ac080z0b5ALOOwF/dvVvo/YXAQHf/j9D7nwC93X1kLdePAq4EFgJL3H2KmfUDzgeaAkvd/YEI110NXA3Qrl2745555hkAtm/fTsuWLZP5LeacdNXhYZ++wNgPq69KX15xBCcH3mdeRbd9WiI/sFWcePqlDR5nLPRzmLhk1uHw2TsoTfBXYdsmMOnUFkmJp6FEq8P+/fu/6+5FtR3P+ASSDEVFRf7OO+8AUFJSQr9+/VL1UTkh3XVY85G4TxbcxSl57/NeeUfOK70rnFweLDuXAioYnYFTfNNdh9kgGXXY49ZX494lt1JjXs8RrQ7NrM4Eko51IOuBw6q87xAqSzozGwQM6tKlSypuL2lSc6Fh18Ba3ivvSI/AGt5ocgMH2FYml53PtfkzGFk6iqG/vodpZzfRUw0lLBmJA3JjnKMu6ZjGuxA40sw6mVkT4FJgRio+SA+UylIRxkTOK72LeRXd6Bz4ghZ8y6j8FxhZOgpATzWUsMqHOiWaPNaMPyfnkwekfhrvNGABcLSZrTOzq9y9DBgJzAY+AJ519+WpjEOy07Szm1Qb8+gaWMsHFYeRbxU051tODKwIJxhAe2flsMrt1RNNHNm2EDBRKe3CcvehtZS/DKR8moy6sLJcjS1PHiw7l2vzZ4Sn+F6X/yJ/Ke8L7N07a+iv7+HKI75m4DXj0xy8NIRkTMcFOOl73+Hp4X2ScKfskvEr0ROhLqzcMO3sJuHkMbJ0FH8sP4/dFOAO5wfm81jBPdW6sxZ+ulWr1rNcx9GJT8eFvTvlKnlEps0UpfHrez2/YRJDX96bJIaV3sSFef/HBXnzaeal3Jn/vxEH12/puZtuF/82zd+AJMNljyzgzU9qLjuLT7tWTbJ29XgyZXULxMwGmdnDW7ZsSXcokmqhgfUetirc2ugXeI/7yoZQSoDOgS9oyW5uyH+uWmvkxff+zd1j1RppzCoX/yUjebRumpfVW48kW1YnEHVh5Z7Kbd0rxzzeqjiW3TSj1APkmdOCPdyZ/79MKZgY7vba4wGeHjuER8aNSnP0jVvx4vWcNP4NOo2exUnj36B4cUpm5wMwtngZP311Bx1HJ/ZMjkqViWPp7XpYcn2oC0uyzrSzm3DtnP+E0lLuL5jMfWVDuC7/RT6sOJhutobOgS8o8wD/lT+dFRUduSH/OSrI55ptfXjj1/fwH20XMeD/nah1I/VQvHg9Y15Yxq7S4Cyn9Zt3MeaFZQBJ24E2Wc/iqKp10zwljQRkdQLRLKwc1fd6HuwL7z/73zz4XrCVcU3pDQD8b8Hvaebfkm8VmDu9Ah8DsKjicI61NYzKfwHbYfzule/Q9fULOH/AyUokMZgwe2U4eVTaVVrOhNkrE04giT75LxKNcSRHVicQd58JzCwqKhqe7lik4XW7+Ld0O7T64PrEsgu5Lv9FdnkB7Sw4NrbHA/QKfEyvwMeUEWBaef9gIqkwfvdKOzq/OoQf9CjUYHsdNmzeVa/yaJI1/bYmTcdNrqxOICKVa0XuvfXn4TGP+8qGMCr/BfZ4HgWUA4ZZ8PR8r+Anea9TRoAvfH/+K386pTThvvcO4b2lQ2hn33BY26Yc/Z+z0/ptZZpD2zZnfYRkcWjb5jFdn4ruqaq0+C81lEAkJ9xw+x9h/iR+98q5wdYFxpWlozkr8BY/yXsddyjHyLfg3735XkF72wTAV96Gm/KfIY8KAjh/31zIO2OHhO/93dZNc37MpP8xB/HUW59GLK9NQzyGWIkjtbI6gWgMRKoJrRd58fWTePbb3gCcl/cPdtKUDysOo3tgNe7BfvzKFkm5Q4dQIgGowDgjsJhyDAMc8B3Gu7M/4sNX/4929g3fYSsf+BFcNvCHOZNU5nz4ZZ3lydq8MBYa32g4Kd/OPRNoO/fkypY6vHvsCA7jCz7xQ8KtkuLyPgzNm0M+FcDeRALBZFLXvPcKCCcVgF3ehF00pTnfst4PpPJW22nO17Sm2yEt+e61f03+N5YGqRqzqA+1NuqvMW7nLpIRbr4juGZk/C0jmFXeJ5xIdtFsbyLxYCJxIM/2vUfVv78qk0vlaS1sDy0I7gJ8lG3Ye03o39J/57Hnt/vjGHlUUIGxw5tRZnl8461ozrc0sz0s884MK70ZSP8gcEN0O9XH5Scezh2Du6c7jJylBCI5r/KBUxN+Uz2RlJFHPhV8VnEAhwU2UbOxbla9hVIpWqO+8pImVr1LJw+nre0E4EDbFi5/s7zb3teffJ1xv8QbWmUXVfCvZyWPdFICEQm58XehJxfOn8SfX+3D920ty7wj5+X9g90UkE85AZxAqA0RKVHUllSqqk+v8R1ll/Foubpm0t3yksiyOoFoEF3i0vd6fhwa/C6cP4k/v7r3N/7B9g0nBd4nn3IMJ69Gz38sySFagqn0z4qjczp5aEwj82V1AtFCQklYlWQChFon+wPwfVvLJlqxP9voEVgdTijO3m6qSGJtgZxgK7kqb1bOJJFJl/RM2rYn0jCyOoGIJF3NhAIwfxJvL3qXgi+Xs4lWHMA2DrLNtGQXpeTRnG/Z4c04OLC13h83Nv9pgKxLItqDKjsogYgkqu/1nNC37lNaPnUhH/3r/Wplh7CJJpRFnYV1UuD9RptADFitrqispQQi0hAuf54jq7wtKSnhyAjz75uE/q1cv53pM67U7ZTblEBEMpgGkiWTZfUDpfREQhGR1MnqBKInEoqIpE5WJxAREUkdJRAREYmLEoiIiMQlJ7ZzN7MvgbWhtwcCX6Xw49oA8Yzax3pdtPPqOh7pWCxlNd+rDutflm11GO3cXK7D+lwb789iQ9XhEe5e+1PB3D2nvoB3Unz/h1N5XbTz6joe6VgsZRHeqw5zvA7jqadcqcP6XBvvz2Km1KG6sJJvZoqvi3ZeXccjHYulLN7vKV6qw8Slug6jnZvLdVifa+P9WcyIOsyJLqyqzOwdr+MJWxKd6jBxqsPEqQ4Tl2gd5mIL5OF0B5AFVIeJUx0mTnWYuITqMOdaICIikhy52AIREZEkUAIREZG4KIGIiEhccjqBmFkLM3vCzB4xs8vSHU9jZGadzexRM3s+3bE0ZmY2OPRzON3Mzkh3PI2NmX3fzKaY2fNm9vN0x9OYhX4vvmNmP4p2btYlEDN7zMw2mtn7NcoHmtlKM/vYzEaHis8Hnnf34cC5DR5shqpPHbr7Kne/Kj2RZrZ61mNx6OdwBHBJOuLNNPWsvw/cfQRwMXBSOuLNVPX8nQhwM/BsLPfOugQCTAWqPWzZzPKAB4CzgGOBoWZ2LNAB+Cx0WnkDxpjpphJ7HUrtplL/ehwbOi71rD8zOxeYBbzcsGFmvKnEWI9mdjqwAtgYy42zLoG4+1zg6xrFJwAfh/5a3gM8A5wHrCOYRCAL6yJe9axDqUV96tGC7gZecfdFDR1rJqrvz6G7z3D3swB1R1dRz3rsB5wI/BgYbmZ1/l7MlUfatmdvSwOCiaM3MBm438zOoeG3SWhsItahmR0A3AkUmtkYdx+Xlugaj9p+Fn8JnAa0MbMu7j4lHcE1ArX9HPYj2CXdFLVAYhGxHt19JICZ/RT4yt0r6rpJriSQiNx9B/CzdMfRmLn7JoL99pIAd59M8A8aiYO7lwAlaQ4ja7j71FjOy5Vum/XAYVXedwiVSexUh8mhekyM6i85klKPuZJAFgJHmlknM2sCXArMSHNMjY3qMDlUj4lR/SVHUuox6xKImU0DFgBHm9k6M7vK3cuAkcBs4APgWXdfns44M5nqMDlUj4lR/SVHKutRmymKiEhcsq4FIiIiDUMJRERE4qIEIiIicVECERGRuCiBiIhIXJRAREQkLkogIjEys3IzW1Llq2OKPud6M7si9HqqmV1Y4/j2Oq5tYmZzzSyntymShqEfMpHY7XL3nrUdNLP80AKtuIV+8Q8DesVzvbvvMbPXCT5T5OlEYhGJRi0QkQSY2U/NbIaZvQG8Hiq70cwWmtlSM7u9yrm3mNm/zGy+mU0zs/+KcMtTgUWxJCIz++8qraH1ZvZ46FAx2tJcGoBaICKxa25mS0KvV7v7kNDrXkAPd//ago+jPZLg8xYMmGFmpwA7CO431JPg/3eLgHcjfMZJEconmNnYmie6+2+B35pZW2AecH/o0PvA8XF+jyIxUwIRiV1tXVh/d/fKB/acEfpaHHrfkmBCaQW86O47Acysto3rDiG4N1FVN7p7+JnzVcdAzMyAp4CJ7v4ugLuXm9keM2vl7tvq9R2K1IMSiEjidlR5bcA4d3+o6glmdn2M99oFNKvHZ98GrHP3x2uUNwV21+M+IvWmMRCR5JoNDDOzlgBm1t7MDgbmAoPNrLmZtQIG1XL9B0CXWD7IzAYRfIrhqBrlBxB8mlxpnN+DSEzUAhFJInf/m5l9H1gQ7F1iO3C5uy8ys+nAe8BGgs9jiOQV4E8xftyvCD6a9O3QZ80IjYv0B2bF/12IxEbbuYukgZndBmx3999HOPYicJO7fxTnvV8ARrv7vxKLUqRu6sISyTyjCQ6m11vo6XLFSh7SENQCERGRuKgFIiIicVECERGRuCiBiIhIXJRAREQkLkogIiISFyUQERGJy/8HTPN2wJgua2gAAAAASUVORK5CYII=\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], "source": [ - "whiten_signal_L1 = whiten_signal(data_L1,asd_interp)" + "# This step is not necessary. Just to show a potential issue: interpolation of ASD has some points out of the curve in left side,\n", + "# I am not sure if this has a big impact on the calculation of the whitening process.\n", + "\n", + "plt.loglog(freq, asd,'-o',label='L1 ASD')\n", + "plt.loglog(freq_rfft, asd_interp, 'x',label='Interpolation')\n", + "\n", + "plt.grid('on')\n", + "plt.ylabel('ASD (strain/rtHz)')\n", + "plt.xlabel('Freq (Hz)')\n", + "plt.legend(loc='upper center')\n", + "\n" ] }, { "cell_type": "code", "execution_count": 8, + "metadata": {}, + "outputs": [], + "source": [ + "whitened_signal_L1 = whiten_signal(data_L1,asd_interp)" + ] + }, + { + "cell_type": "code", + "execution_count": 9, "metadata": { "scrolled": true }, @@ -184,16 +223,16 @@ { "data": { "text/plain": [ - "[]" + "Text(0, 0.5, 'Waveform')" ] }, - "execution_count": 8, + "execution_count": 9, "metadata": {}, "output_type": "execute_result" }, { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAD4CAYAAADvsV2wAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAcdUlEQVR4nO3de3xc5X3n8c9vZnSzLfluY8BYdjAEwsVgBTAJDVtIAmyLyYYEyCYLKZeEwIbXNtsWmt2EkG5Ddl8taRcW6kIaSBsu4dKY1oRwTwKYWIAvwdxsY7AM2PJFtmxLI83Ms3/MmYtkSTNnRtbM4/m+Xy+/mMuZOT8Po68f/c5znmPOOURE5OAXqXQBIiIyNhT4IiI1QoEvIlIjFPgiIjVCgS8iUiNilS5gONOmTXOtra2VLkNExCsvv/zyNufc9KGeq9rAb21tpb29vdJliIh4xczeHe45tXRERGqEAl9EpEYo8EVEaoQCX0SkRijwRURqhAJfRKRGKPBFRGqEAl/EQ4+t+YDte+KVLkM8o8AX8UxPX5Kr/+UVHnl1c6VLEc8o8EU805dMAdCf1MWLJBwFvohnkql00Kd0tToJSYEv4plEKj3CT6UU+BKOAl/EM0Heo7yXsBT4Ip7JjvDV0pGQFPginlEPX0qlwBfxTEKBLyVS4It4JpUN/AoXIt5R4It4RiN8KZUCX8Qz2R6+hvgSkgJfxDNJtXSkRAp8Ec+opSOlUuCLeCYzwlfeS1gKfBHPZE68SqqnIyEp8EU8k1taQYEv4SjwRTyTW1qhwoWIdxT4Ip7J9fCV+BKOAl/EM1pLR0qlwBfxTCbwgwtfiRRtVALfzH5sZlvN7PfDPG9m9vdmts7MVpvZyaOxX5FalFBLR0o0WiP8nwDnjPD8ucD84M9VwO2jtF+RmqOWjpRqVALfOfdrYMcImywG7nFpy4FJZjZrNPYtUmu0tIKUaqx6+IcBm/LudwSPDWBmV5lZu5m1d3Z2jlFpIn7J9vA1wpeQquqgrXNuiXOuzTnXNn369EqXI1KV1MOXUo1V4G8GZufdPzx4TERCyozsU5qlIyGNVeAvBf5LMFvnNGCXc+6DMdq3yEElmdRFzKU0sdF4EzO7FzgTmGZmHcB3gToA59wdwDLgPGAdsA/46mjsV6QWaXlkKdWoBL5z7pICzzvgmtHYl0itywS9ZulIWFV10FZECtMIX0qlwBfxTDKpEb6URoEv4pmELmIuJVLgi3gm18NX4Es4CnwRz6iHL6VS4It4RmvpSKkU+CKeSaqHLyVS4It4RssjS6kU+CKe0UXMpVQKfBHP6CLmUioFvohndNBWSqXAF/FMInsRcyW+hKPAF/GMDtpKqRT4Ip7J9fArXIh4R4Ev4hmN8KVUCnwRzyR0EXMpkQJfxDNq6UipFPginlFLR0qlwBfxjAJfSqXAF/FMdmmFVIULEe8o8EU8kwl6jfAlLAW+iGdyi6cp8CUcBb6IZ7SWjpRKgS/imcz8e10ARcJS4It4JpHULB0pjQJfxDNq6UipFPgintE8fCmVAl/EM+rhS6kU+CKeyfXwK1yIeEeBL+IZtXSkVAp8Ec9kWjrKewlLgS/imaTWw5cSKfBFPJNIamkFKY0CX8Qz+RdAcQp9CUGBL+KZ/FaO8l7CUOCLeCaZNx9TbR0JY1QC38zOMbM3zWydmV0/xPOXmVmnma0M/lwxGvsVqUWJlCMaMUAHbiWcWLlvYGZR4Dbg00AHsMLMljrn1g7a9H7n3LXl7k+klqVSDuegLmYkg9sixRqNEf4pwDrn3AbnXB9wH7B4FN5XRAbJjOjroukfXbV0JIzRCPzDgE159zuCxwb7vJmtNrMHzWz2UG9kZleZWbuZtXd2do5CaSIHl0z/vj4b+JWsRnwzVgdtHwVanXMnAE8Adw+1kXNuiXOuzTnXNn369DEqTcQfidTAEX5SiS8hjEbgbwbyR+yHB49lOee2O+fiwd07gYWjsF+RmpMJ+Fg0fdBW8/AljNEI/BXAfDOba2b1wMXA0vwNzGxW3t3zgddHYb8iNUctHSlH2bN0nHMJM7sWeByIAj92zr1mZjcB7c65pcA3zex8IAHsAC4rd78itSiRSi+rUB/TQVsJr+zAB3DOLQOWDXrsO3m3bwBuGI19idSyIO9zs3RSjv5kiv5kinH1o/LjLAcxnWkr4pHMCD/Tw085+H/PrOdzt71QybLEEwp8EY8kU/vPw/9wdy/vd/VUsizxhAJfxCP7H7R1JFMp4olUJcsSTyjwRTySG+EHLZ0UJFPQl0xpTr4UpMAX8cjgE68yI3yAPo3ypQAFvohHsiP8vGmZyWBgH08kK1WWeEKBL+KRoU68yozw1ceXQhT4Ih5JDO7hO0ciGOLH+xX4MjIFvohHhpqWmTnbtlctHSlAgS/ikf0CP5Ub9WuEL4Uo8EU8ku3h5x+0zQS+RvhSgAJfxCOZpRXqsssjkxf4GuHLyBT4Ih4Z3NJJOpdr6WiELwUo8EU8MuRBW/XwpUgKfBGPDJ6H7/JG+JqlI4Uo8EU8sv88/NxFUDTCl0IU+CIeyYR7ZmmFZCrvxCsdtJUCFPgiHsmEe11E0zIlPAW+iEdyi6flTctUS0eKpMAX8Ugie9A2Cgwc4eugrRSiwBfxSGY0nzlom0zltXQ0wpcCFPgiHkkmM2faZqZl6kxbKZ4CX8QjQ13xKpFdD18tHRmZAl/EI5lpmbnF09LXtAWN8KUwBb6IRzIj/NiAHn4wwlcPXwpQ4It4JJnU0gpSOgW+iEcys3Ri+UsraJaOFEmBL+KRZMoRMYhF8q5pqzNtpUgKfBGPJFKOaMSIWC7ws4un6aCtFKDAF6kSvf1JVnd0jbhNcojAT2gevhRJgS9SJe55cSPn3/o8b2/pHnabZMoRi0SygZ9IOoIBPr39aunIyBT4IlXiN29vA+ChVzYPu02mhx/kfXZ0DxrhS2EKfJEq0JdI0b5xJwCPvNqRXS5hsEQqRSwaIRoctO1P5kI+rhG+FKDAF6kCqzq66OlPsnjBoWzZHefZN7fyxNotdHbHB2yXTDGgh98XjOojlh7hJ5Ipvv9va+nYuW/M/w5S/WKVLkCkVjnnuOyfVnD0Ic2Mr49hBt8+7xieeWMrV9zTjnPwtU/N44Zzj8m+JplKEYsYwQCf/uBErPH1MbrjCTZs28tdv32H/mSKmxYfV4m/llSxURnhm9k5Zvamma0zs+uHeL7BzO4Pnn/JzFpHY78iPurpS7deVm7q4rm3Olny6w3c8+JGPnZoCzNaGvnmWfM566MzaWmMsWNP34DXJlKOiBlmA1s6TfXp9fE/2NULwKOr3s+O/gGeXLuFd7fvBeDDXb38YmXuOEEimcK5oVtIcnApO/DNLArcBpwLHAtcYmbHDtrscmCnc+5I4Bbgh+XuV6Sadff2Z2+/sG4bKzelp1s+9HIHJ970K17asJ2fv9xBY12E4w+byPa9fSyaNxWAK86Yx52XtjFrYhO7894Hglk6Ucv28BNB4I9vSP+y/n5XDwA79/Xz3FudAGzasY+rftrOTY+uBeBvfvUm1923kjc/7GZXTz+n/vVT3PXbd9K1rt/G3/zqTZxzOOd4bM0H7NqXrmF3b3/2/QH29SWytzPbS3UbjZbOKcA659wGADO7D1gMrM3bZjFwY3D7QeBWMzN3AL4hfYkUyzdsD/26UgoppfyS/sIlvMiVtidK+T9S0mvCv6TkQOlNpOja18f4+hgzWhr4oKuX3b39HHfYRJxL98/nThvPwjmT+fVbnWzbE+eCBYexq6efpaveZ9G8qZwydwr3/u49tnbH+caZR/LOtr3836ff5tLTW1kwexI3PLyGnr4kt1y0gPtXvMf3/m0t/+uC4zlx9kQu+6cVNNRFePjq07n5l2/Ql0hxwyNr6OyOc+5xs7jurPl8419e4Y9OOHRA3S1NMXb17B/40byWTl/Q0hkXjPAzgdxUF+WRVzv49LEzufuFjaQcPPdWJ5u7enjs9x8C8PP2TRw6qYnte/u4/dn1fKFtNn/289Vs7urh6EOa2RdP8ucPreaM+dO4/csLufD2F9i8s4cHvr6I1z/o5i8eWs2VZ8zjyjPmcs3PXqGnL8nfX3ISL72zg7t+8w6Xnt7KJ4+cxm3PrKMuZnz1E3NZ+V4XL72znT/86AwOnzyOp9/YyoSGGKfMnULHzn107OzhqJnNjG+I0rGzh4ZYhMnj6ulLpuhLpILfbiBiFvwp6StR4vev1H2Ff2FzYx0LZk8qbYcjsHIz18wuBM5xzl0R3P8KcKpz7tq8bX4fbNMR3F8fbLNt0HtdBVwFcMQRRyx89913Q9ezfU+chX/1ZKl/HREgvXRB/pTHaRPq2Ra0Vw6f3MTW7jh9iXQ//cgZE3gzmDt/9Mxm3tzSzfj6GL39SWa2NNKXTLEvnsDM2BNP8N/OPopbnnwLgJ9dcSqnHzltyBquuHsF73f1suy6M7KPXf3PL7Nu6x4evPp0Tvzer7j8k3O567fv8PHWyazYuJPPn3w4D7/awZdPncP9KzbxD19ZyDfvfZXWaeNZs3kXp82bwvINOzh8chM9fUlamurYG0+wtTvOibMnsWpTF9MmNFAXNfb1JamLRti2J07r1HG8u2MfU8fXk3LQta+P6c0NbNkdp7khRl8yRUMsQk9/kv6kY9qEBrbtSR9wbqyLkEpBX/DbSGNdhF6t+zOiBbMn8a/XfKKk15rZy865tqGeq6qDts65JcASgLa2tpL+JWppquOhqxeVWEH44YKVMMIoZVBipeyo5H2Vsp+x+exK0RCLMGlcPd29/WztjjNrYiPjG2Ks2tSFGZx4+CRe/6CbV9/byelHTmXK+AbuX7GJCQ1RLlw4m1+s3Mxv3t7GVxbNobkxxl8+vIY/mD+d//lHx/Ddpa/x/Lpt/OiiBaSc41sPrOKk2ZO44ysL+c//+BJvb93DP19+Ku/v6uHPH1zNfzxhFtedPZ/3u3p45b2dnBa0cYbS0ljHG70DT8JKDBrhZ3r44+pzLZ1JTXV8/cyP8Py6bXz1JysA+P4Fx3H9Q6tZvmEHsyY28t0//hhX3tPO9r193HLRifz0xXd55b0uTp07hT/99FFctGQ59dEID113Oj/85Rs8sXYL//0zR3H2sTP5wh0vctIRk7nnT07hjufW868rN/Oji05iRnMD33t0LafOncKffHIuS1dt5q0te7h0USsAD7/awfGHTWTRvKn8Zt02tu/p48yjp7MvnuTl93YwZ+p4jpgyjre2dNPbn2T25HH0Jx079/VRH4tQH43gSJ9d7JwjmUrfLvVrVMrPVKnf2bAvy/z/HG2jMcJfBNzonPtscP8GAOfcD/K2eTzY5kUziwEfAtNHaum0tbW59vb2smoTORCcc9mwcM6xuyfBxHF1AKzb2s2hk5oYVx9j+544G7bt5eOtU3DO8cTaLZw6dyoTx9WRSqV/0Y+O0JO4celrPPRKB2tu/Gz2sct/soIt3b088LVFHPudx7nklNnc+7tNnHf8ISxb8yGtU8cRjRhPfetMdvf285cPrwHg1i+dzD88t54fPPYGX/vUPP7sM0ez6OanSaUcL9zwh7y4fjtX3tPOfVedxsI5U7j92fUcOqmRxQsOo7u3n+fXbeczx84kEjG69vUxoSFGLKpZ3dXoQI/wVwDzzWwusBm4GPjSoG2WApcCLwIXAk8fiP69yFjIHxmaWTbsAY6c0Zy9PXVCA1MnNGS3+8zHDsk+Fymi+dzSVMeeeIJUymW3T4/wI3nz8DM9/GCEv6s32/ttaazj1i+dnH2/CxcezoqNO/jyqXOIRSP86KIFADTEopx59AxWf/ez2dk+V5/5kezrmhvrOOe4XO2TxtUXrF2qU9mB75xLmNm1wONAFPixc+41M7sJaHfOLQXuAn5qZuuAHaT/URCREbQ0xnAOuntzv0EkU45o3tIKuZZOOqj7Eimmjh86kKdOaODOSz+evf+JQccOMmEvB69RaRQ555YBywY99p28273AF0ZjXyK1YmJTOuR39/ZnAz+RSg1cPC01sIcPMGWYwBdRE06kSrUEgZ8/NTMVLK0QHdTSGZ83Os+0kUQGU+CLVKnsCD8v8PtTKaIR27+l05Ab4Q/X0hFR4ItUqZbGXEsno6cvSVN9FAtOQMoEfv4IXy0dGY4CX6RKZfr2+S2dPfEEzcFoPmK231o6AFMnKPBlaAp8kSrV0pgO9t09uTVr9sYT2XVzIjZwtcyMqePVw5ehKfBFqtT4+hgRG9jS2TMg8C27Iua4BrV0pLCqWlpBRHIiEaOlqS7b0okn0uvUNDcO0dKpi2YP5E7OOxFMJJ8CX6SKtTTWZWfp7I2n19HPHKCNWO6atrFIhIZYhKa6qJY8kGHpmyFSxSbmjfD39KZ7+UO1dCKR9BIJaufISBT4IlWspSnG7iDo98TT/822dCKWPdM2M8LXSVcyEgW+SBWb2JTX0ukbPMLPXcQ8GkmvpzNdgS8jUA9fpIq1NOa1dOL7t3Ti2cCP8NefO55pzQp8GZ4CX6SKtTTVZadlZnr4mROvLG+WTixiw145SyRDLR2RKjaxqY7e/hTxRJK9g0b40UhuaYVi1tcXUeCLVLH8s22HaulkLrsbU+BLERT4IlWsJW9N/EzgT8gL/IzIWF0gWLymwBepYvlr4u+NJ2iqi2avg5uf8RrhSzEU+CJVbHJw/diufX3siSez7RwYeAH0aFSBL4Up8EWq2PRgmuXW3fH00siNucDPb+NE1dKRIijwRapY5kSqrd3xYGnk3KqY+RkfVUtHiqDAF6li9bEIk8fVsbW7N700cv0wI3wFvhRBgS9S5WY0N6ZbOr0DWzpRtXQkJAW+SJWb0dLAlu44e/sSAw7aZjLeTCdeSXEU+CJVbkZzI527ewdc3hByLR1NyZRiaS0dkSo3o6WBzj1xzCy7jg6k18AH9e+leAp8kSo3o7khuFi5GzgPPxjhq38vxVJLR6TKzWhuzN4e2MMPAl8jfCmSAl+kys1oya1xP6ClE+S8Al+KpcAXqXIz8i5qMtRB22hEP8ZSHH1TRKpcfktnwhBLK2iWjhRLgS9S5Zrqo9kTribkLa2gWToSlgJfxAOZts7QLR0FvhRHgS/igUxbZ4ICX8qgwBfxQGamzoQhllZQ4EuxFPgiHhiqpZMJep14JcXSmbYiHvhi22xmtjRSF82N0dTSkbDKCnwzmwLcD7QCG4EvOud2DrFdElgT3H3POXd+OfsVqTXzZzYzf2bzgMcyOR/T5Q2lSOW2dK4HnnLOzQeeCu4Ppcc5tyD4o7AXGQWZpRUiaulIkcoN/MXA3cHtu4ELynw/ESlSVCdeSUjlBv5M59wHwe0PgZnDbNdoZu1mttzMhv1HwcyuCrZr7+zsLLM0kYNb5sQrXfxEilWwh29mTwKHDPHUt/PvOOecmblh3maOc26zmc0DnjazNc659YM3cs4tAZYAtLW1DfdeIkKupaMRvhSrYOA7584e7jkz22Jms5xzH5jZLGDrMO+xOfjvBjN7FjgJ2C/wRaR4mqUjYZXb0lkKXBrcvhT4xeANzGyymTUEt6cBnwDWlrlfkZoX1YlXElK5gX8z8Gkzexs4O7iPmbWZ2Z3BNscA7Wa2CngGuNk5p8AXKZNWy5SwypqH75zbDpw1xOPtwBXB7ReA48vZj4jsT9MyJSwtrSDiKZ14JWEp8EU8lenda4QvxVLgi3hK0zIlLAW+iKdyFzHXj7EUR98UEU/l5uFXuBDxhr4qIp7SCF/C0jdFxFOZNXQ0wpdi6asi4qnciVf6MZbi6Jsi4qmIllaQkBT4Ip7S4mkSlgJfxFO5Hr4CX4qjwBfxVLalozNtpUgKfBFPqaUjYSnwRTylpRUkLAW+iKcyrRxd01aKpcAX8VR2eWQFvhRJgS/iKc3SkbAU+CKeMp14JSEp8EU8FdUsHQlJgS/iKU3LlLAU+CKeMh20lZAU+CKeyozwdU1bKZYCX8RTmVZOLKrAl+Io8EU8lenkaIQvxVLgi3jKdAEUCUnfFBFPaZaOhKXAF/FU5lq2CnwplgJfxFNaLVPCUuCLeCqi1TIlJAW+iKe0WqaEpcAX8ZROvJKwFPginoroxCsJSYEv4qmIlkeWkBT4Ip7KzsNXS0eKpMAX8ZROvJKwygp8M/uCmb1mZikzaxthu3PM7E0zW2dm15ezTxFJW/SRqXztU/M4+pDmSpcinih3hP974D8Bvx5uAzOLArcB5wLHApeY2bFl7lek5k1squOGc4+hLqpf1KU4sXJe7Jx7HXJn/A3jFGCdc25DsO19wGJgbTn7FhGRcMZiaHAYsCnvfkfw2H7M7Cozazez9s7OzjEoTUSkdhQc4ZvZk8AhQzz1befcL0azGOfcEmAJQFtbmxvN9xYRqXUFA985d3aZ+9gMzM67f3jwmIiIjKGxaOmsAOab2VwzqwcuBpaOwX5FRCRPudMyP2dmHcAi4N/N7PHg8UPNbBmAcy4BXAs8DrwOPOCce628skVEJKxyZ+k8AjwyxOPvA+fl3V8GLCtnXyIiUh5N4BURqRHmXHVOhjGzTuDdMt5iGrBtlMoZC6r3wFK9B5bqPbDC1DvHOTd9qCeqNvDLZWbtzrlhl3uoNqr3wFK9B5bqPbBGq161dEREaoQCX0SkRhzMgb+k0gWEpHoPLNV7YKneA2tU6j1oe/giIjLQwTzCFxGRPAp8EZEa4XXgF7qSlpk1mNn9wfMvmVnr2Fc5oJ5C9f6Bmb1iZgkzu7ASNQ6qp1C9f2pma81stZk9ZWZzKlFnXj2F6v26ma0xs5Vm9ttquBBPsVeDM7PPm5kb6cpyB1oRn+9lZtYZfL4rzeyKStSZV0/Bz9bMvhh8h18zs5+NdY1D1FPoM74l7/N9y8y6Qu3AOeflHyAKrAfmAfXAKuDYQdt8A7gjuH0xcH+V19sKnADcA1zowef7H4Bxwe2rPfh8W/Junw/8sto/42C7ZtJXlVsOtFVrrcBlwK2V/ExD1jsfeBWYHNyfUe01D9r+vwI/DrMPn0f42StpOef6gMyVtPItBu4Obj8InGUFLs91ABWs1zm30Tm3GkhVosBBiqn3GefcvuDuctJLX1dKMfXuzrs7Hqj0jIVivsMA3wd+CPSOZXGDFFtrtSim3iuB25xzOwGcc1vHuMbBwn7GlwD3htmBz4FfzJW0stu49Kqdu4CpY1Ld/oq+8leVCFvv5cBjB7SikRVVr5ldY2brgf8NfHOMahtOwZrN7GRgtnPu38eysCEU+334fNDie9DMZg/x/Fgppt6jgKPM7HkzW25m54xZdUMLc3XAOcBc4OkwO/A58KVKmNmXgTbg/1S6lkKcc7c55z4C/AXwPypdz0jMLAL8LfCtStdSpEeBVufcCcAT5H67rlYx0m2dM0mPlv/RzCZVtKLiXQw86JxLhnmRz4FfzJW0stuYWQyYCGwfk+r259uVv4qq18zOBr4NnO+ci49RbUMJ+/neB1xwQCsqrFDNzcBxwLNmthE4DVhaoQO3BT9f59z2vO/AncDCMaptKMV8HzqApc65fufcO8BbpP8BqJQw3+GLCdnOAbw+aBsDNpD+tSZzgONjg7a5hoEHbR+o5nrztv0JlT9oW8znexLpg0zzPfk+zM+7/cdAe7XXPGj7Z6ncQdtiPt9Zebc/Byyv5s8WOAe4O7g9jXQ7ZWo11xxs91FgI8GJs6H2Uam/3Ch9QOeR/ld5PemLqgPcRHq0CdAI/BxYB/wOmFfl9X6c9KhjL+nfRF6r8nqfBLYAK4M/S6u83r8DXgtqfWakcK2WmgdtW7HAL/Lz/UHw+a4KPt+PVvNnCxjpltlaYA1wsQ/fB+BG4OZS3l9LK4iI1Aife/giIhKCAl9EpEYo8EVEaoQCX0SkRijwRURqhAJfRKRGKPBFRGrE/weAFXRIbaLa7AAAAABJRU5ErkJggg==\n", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYoAAAEWCAYAAAB42tAoAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3deZhcZZn38e9dXb0l6c6ekEBIAgn7asIOiq/IoOOwCAKOC6iAoAzOgKMwviNu4/LOjDqOjBqREdARlEXDDI6IsgcwDQQCQSCEJQuQfe+tqu73j3OqurpTXX2q0tXVXef3ua66qDp1qs7dRaV//TzPOc9j7o6IiEh/EtUuQEREhjcFhYiIFKWgEBGRohQUIiJSlIJCRESKUlCIiEhRCgqpaWbmZjann+c+ZGb3DHVNhZjZT83sa2W8bruZ7VOJmvoc51UzO6XSx5HhSUEhI4aZXWNmv+2z7aV+tp0/0Pu5+8/d/dS81/UbKsOVu49x9xXVrkNqm4JCRpIHgePNrA7AzKYB9cCRfbbNCfcVkUGgoJCRZDFBMBwRPj4JuA94oc+2l919Td7rTglbGZvN7DozMwAzu9DMHg7vZ4Pl6bA757xw+/vMbEn42kVmdlj2TcPumM+a2TNmtsXMbjWzprzni732SDN70sy2mdmtQO51fZnZHDN7IDzG+nD/7HO5VpCZTTSzu8xsq5ktNrOvZX++vH0v7eez2NfM/mhmG8Jj/NzMxkX6vyI1T0EhI4a7dwGPA28PN70deAh4uM+2vq2J9wFHAYcB5wJ/UeC9s68/POzOudXMjgRuAD4JTAR+BCw0s8a8l54LnAbMDt//QgiCoL/XmlkD8GvgZmAC8Cvg7CI/+leBe4DxwF7Av/ez33XADmAP4ILw1ld/n4UB3wCmAwcCM4AvFalJYkRBISPNA/SEwkkEQfFQn20P9HnNN919s7u/TtACOYJoLgF+5O6Pu3va3W8EOoFj8/b5nruvcfeNwF15713stccStIy+6+7d7n4bQWupP93ATGC6u3e4+8N9dwi73s4GrnX3ne6+DLixwHsV/Czcfbm7/97dO919HfBt4B2RPiWpeQoKGWkeBE40swnAZHd/CVhEMHYxATiEXVsUb+bd3wmMiXismcBVYTfNZjPbTPCX9vQI713stdOB1d57Rs7XitTxOYK/+P9kZs+Z2ccL7DMZSAIr87atLLBfwXrNbKqZ3WJmq81sK/AzYFKRmiRGFBQy0jwKjAUuBh4BcPetwJpw2xp3f2WQjrUS+Cd3H5d3G+Xuv9jN174B7JkdHwjt3d8bufub7n6xu08n6Mr6jwJnZ60DUgRdU1kzovyQoa8DDhzq7q3AhwnCSURBISOLu7cDbcCVBF1OWQ+H23bnbKe3gPxrEn4MXGpmx1hgtJn9pZm1RHivYq99lOCX+hVmVm9m7weO7u+NzOwDZpYNgE0Ev9Az+fu4exq4A/iSmY0yswOAj0b8uQFagO3AFjPbE/j7El4rNU5BISPRA8AUgnDIeijctjtB8SXgxrCr6Fx3byNopXyf4Bf0csLB6oEUe204KP/+8PFG4DyCX/L9OQp43My2AwuBz/Rz7cTlBK2tNwkGyn9BMC4SxZeBtwFbgP8ZoB6JGdPCRSK1ycy+Bezh7oXOfhKJTC0KkRphZgeY2WFhV9fRwCeAO6tdl4x8yWoXICKDpoWgu2k6wXjLvwK/qWpFUhPU9SQiIkWp60lERIqqua6nSZMm+axZs6pdhojIiPLEE0+sd/fJhZ6ruaCYNWsWbW1t1S5DRGREMbN+ZwdQ15OIiBSloBARkaIUFCIiUpSCQkREilJQiIhIUQoKEREpSkEhIiJFKShEYuS3S99gw/aoM4+LBBQUIjHR3pXmsp8/yZ1Pra52KTLCKChEYqIrHSyK153WRKBSGgWFSEykM0FAZDRjtJRIQSESE6lM0KLIZBQUUhoFhUhMhDmBckJKpaAQiYlci0JdT1IiBYVITGiMQsqloBCJiZSCQsqkoBCJiUwuKKpciIw4CgqRmFCLQsqloBCJidwYhZoUUiIFhUhMpNX1JGVSUIjEhLqepFwKCpGYyLYolBNSKgWFSExkL7hLq+9JSqSgEImJnik8FBRSGgWFSEz0TOFR5UJkxFFQiMREzxiFkkJKo6AQiQnN9STlUlCIxEQ2KMKF7kQiq2pQmNkNZrbWzJ7t53kzs++Z2XIze8bM3jbUNYrUipS6nqRM1W5R/BQ4rcjz7wHmhrdLgB8MQU0iNUldT1KuqgaFuz8IbCyyyxnATR54DBhnZtOGpjqR2qIpPKRc1W5RDGRPYGXe41Xhtl7M7BIzazOztnXr1g1ZcSIjSW6MQi0KKdFwD4pI3H2Bu8939/mTJ0+udjkiw5LGKKRcwz0oVgMz8h7vFW4TkRJlWxIZnfUkJRruQbEQ+Gh49tOxwBZ3f6PaRYmMROl09spstSikNMlqHtzMfgGcDEwys1XAtUA9gLv/ELgbeC+wHNgJfKw6lYqMfJpmXMpV1aBw9w8O8LwDnx6ickRqWjYgdNaTlGq4dz2JyCBRi0LKpaAQiYl0Wi0KKY+CQiQmci0KJYWUSEEhEhM9YxQKCimNgkIkJjRGIeVSUIjEhOZ6knIpKERiIq0xCimTgkIkJjTNuJRLQSESE6lMdgqPKhciI46CQiQm0po9VsqkoBCJCQ1mS7kUFCIxkT09Nq2kkBIpKERiQoPZUi4FhUhM9IxRVLkQGXEUFCIxoRaFlEtBIRITuTEKBYWUSEEhEhPqepJyKShEYkJdT1IuBYVITCgopFwKCpGYyE3hkalyITLiKChEYiIbEGpRSKkUFCIx0TMpoIJCSqOgEIkJzfUk5VJQiMRE9voJLVwkpVJQiMREKq2znqQ8CgqRmFDXk5RLQSESE7qOQsqloBCJCY1RSLkUFCIx0TNGUeVCZMRRUIjEhLqepFwKCpGYyHY9KSekVAoKkZhIaz0KKZOCQiQmUmlN4SHlUVCIxET+wkWusJASKChEYiK/y0k5IaVQUIjERDrvvFh1P0kpqhoUZnaamb1gZsvN7OoCz19oZuvMbEl4u6gadYrUglTGqUsYoAFtKU2yWgc2szrgOuDdwCpgsZktdPdlfXa91d0vH/ICRWpIJuO4Q33SSIf3RaKqZoviaGC5u69w9y7gFuCMKtYjUrOyLYj6uuCfvLqepBTVDIo9gZV5j1eF2/o628yeMbPbzGxGoTcys0vMrM3M2tatW1eJWkVGtOz4REMuKKpZjYw0w30w+y5glrsfBvweuLHQTu6+wN3nu/v8yZMnD2mBIiNBKtO7RZFWUkgJqhkUq4H8FsJe4bYcd9/g7p3hw+uBeUNUm0hNyQZDsi4YzNZ1FFKKagbFYmCumc02swbgfGBh/g5mNi3v4enA80NYn0jNUNeT7I6qnfXk7ikzuxz4HVAH3ODuz5nZV4A2d18IXGFmpwMpYCNwYbXqFRnJUplg+o6GpAazpXRVCwoAd78buLvPti/m3b8GuGao6xKpNWFO9Jz1lHG60xm60xlGNVT114CMAMN9MFtEBkG2RZEdo8g4/Md9L3PWdYuqWZaMEAoKkRhIZ3a9juLNrR2s2dxezbJkhFBQiMTAroPZTjqToTOVqWZZMkJE6pw0s3HAR4FZ+a9x9ysqU5aIDKaeFkXY9ZSBdAa60hnSeXNAiRQSdRTrbuAxYCmgP0FERpi+F9xlWxQAXakMzQ11VatNhr+oQdHk7ldWtBIRqZhciyLv9Nh0eIZsZyqtoJCioo5R3GxmF5vZNDObkL1VtDIRGTSFLrjLtig0TiEDidqi6AL+GfgCkL1Sx4F9KlGUiAyuVN8xCndSYZOis1tBIcVFDYqrgDnuvr6SxYhIZRQ6PTZ7dXZHKl21umRkiNr1tBzYWclCRKRydgmKTE8rQy0KGUjUFsUOYImZ3QdkZ3PV6bEiI0RujCJ/MDsbFGpRyACiBsWvw5uIjEDZKTzqc9OMkxcUalFIcQMGRbi29anu/qEhqEdEKqBv11PavafrSS0KGcCAYxTungZmhmtGiMgIVHAwW2MUElHUrqcVwCNmtpBgvAIAd/92RaoSkUHV9zoKz2tR6KwnGUjUoHg5vCWAlsqVIyKVsOt1FD2LF6lFIQOJFBTu/mUAMxsTPt5eyaJEZHBlQyE7hUc6k3fBnQazZQCRrqMws0PM7CngOeA5M3vCzA6ubGkiMliyoVCf0OmxUrqoF9wtAK5095nuPpPgSu0fV64sERlMPZMC5p0eq64niShqUIx29/uyD9z9fmB0RSoSkUGXyg1mB7PE5rcoNJgtA4l81pOZ/SNwc/j4wwRnQonICJBtPWQHs9OZvK4ntShkAFFbFB8HJgN3ALcDk8JtIjICpNPZK7Ozp8fqymyJrmiLwsxudvePAB/VvE4iI1ehFe5SufUo1PUkxQ3UophnZtOBj5vZ+PxFi7RwkcjIkT09tmdSwGDNbFCLQgY20BjFD4E/ECxQ9GSf57RwkcgIkW1RJHuNUYQtCo1RyACKtijc/XvufiBwg7vP7nNTSIiMEOm0pvCQ8kUazHb3y8zsRDP7GICZTTKz2ZUtTUQGS/asp2T+FB4660kiinpl9rXA54Frwk0NwM8qVZSIDK50xkkYJBN5a2brymyJKOrpsWcBpxPOHOvua9DkgCIjRirj1CWMhPUERW5SQA1mywCiBkWXuzvBADZmpquyRYaJju40z6zaXHSfdIGgSOk6CokoalD80sx+BIwzs4uBe9FcTyLDwk2Pvsrp33+El97a1u8+6YyTTCRyQZFKO2GDgo5udT1JcVEHs/8FuI3gquz9gS+6+79XsjARieahl9YDcPuTq/vdJztGEeZErjUBalHIwKIOZl8JLHP3v3f3z7r77ytcl4hE0JXK0PbqJgDufGpVblqOvlKZDMm6BHXhYHZ3uiccOtWikAFE7XpqAe4xs4fM7HIzm1rJokQkmqdXbaa9O80ZR0znra2d3P/CWn6/7C3WbevstV86Q68xiq6wFZGwoEWRSmf46n8vY9WmnUP+M8jwF7Xr6cvufjDwaWAa8ICZ3VvRykSkIHfnghv+xNfvfp5FyzdgBl9474G0NiW56KY2Lr6pjesf7j25czqTIZkwwgYF3eEFeKMbknSmMqxYv4OfPPwKCx7UpNCyq6gtiqy1wJvABmDK7h7czE4zsxfMbLmZXV3g+UYzuzV8/nEzm7W7xxQZqdq7gi6iJSs388CL61jw4ApuevRVDp7eypTWJq5411zedcBUWpuSbNze1eu1qYyTMMOsd9dTc0OwPsUbWzoAuOvpNbnWBsC9y97itQ07AHhzSwe/WdIzDpJKZ3Av3NUltSXqGMWnzOx+gnmfJgIXu/thu3NgM6sDrgPeAxwEfNDMDuqz2yeATe4+B/gO8K3dOabIcLetozt3f9Hy9SxZGZz2evsTqzj8K/fw+IoN/OqJVTTVJzh0z7Fs2NHFcftMBOCik/bh+gvmM21sM1vz3gfCs57qLDdGkQqDYnRjMN3bms3tAGza2c0DL64DYOXGnVxycxtfuWsZAP96zwt85pYlvPDmNra0d3PM1//ATx5+Jaj15fX86z0v4O64O79d+gZbdgY1bO3ozr0/wM6uVO5+dn8Z3qIuXDQD+Ft3XzKIxz4aWO7uKwDM7BbgDGBZ3j5nAF8K798GfN/MzCvwzepKZXhsxYaSX1dOIeWUX9YPXMaLvLwjUc7/kbJeU/pLyv5F1JHKsHlnF6MbkkxpbeSNzR1s7ejmkD3H4h6MD8yeNJp5M8fz4IvrWL+9kzOP2JMt7d0sfHoNx+0zkaNnT+AXf3qdtds6+dTJc3hl/Q7+/Y8vccHxszhixjiuuWMp7V1pvnPeEdy6+HW+/N/L+KczD+XwGWO58D8X01if4I7Ljueb//tnulIZrrlzKeu2dfKeQ6bxmXfN5VM/f5L3HTa9V92tzUm2tO8aFHV5XU9dYdfTqLBFkf1F3lxfx51PreLdB03lxkWvknF44MV1rN7czm+ffROAX7WtZPq4Zjbs6OIH97/MB+bP4O9/9QyrN7ez/x4t7OxM87nbn+GkuZP4wYfncc4PFrF6Uzu/vPQ4nn9jG5+//RkuPmkfLj5pNp/+rydp70rzvQ8eyeOvbOQnD73CBcfP4sQ5k7juvuXUJ42PnTCbJa9v5vFXNvB/DpjCXuNH8cc/r2VMY5KjZ09g1aadrNrUzn5TWxjdWMeqTe00JhOMH9VAVzpDVyoTtqYgYRbeyvpKlPn9K/dYpb+wpameI2aMK++ARVgp/4jMbArQlH3s7q+XfWCzc4DT3P2i8PFHgGPc/fK8fZ4N91kVPn453Gd9n/e6BLgEYO+995732muvlVzPhu2dzPuahl1k9yQT1uvU00ljGlgfdgPtNb6Ztds66UoF4wVzpozhhfDah/2ntvDCW9sY3ZCkozvN1NYmutIZdnamMDO2d6b4u1P24zv3vgjAf110DMfPmVSwhotuXMyazR3c/ZmTctsu+9kTLF+7ndsuO57Dv3wPnzhxNj95+BWOmjWexa9u4uy37cUdT63iw8fM5NbFK/nRR+ZxxS+eYtak0SxdvYVj95nAYys2stf4Ztq70rQ217OjM8XabZ0cPmMcT6/czKQxjdTXGTu70tTXJVi/vZNZE0fx2sadTBzdQMZh884uJrc08tbWTloak3SlMzQmE7R3p+lOO5PGNLJ+ezAQ31SfIJOBrrD101SfoEPzUhV1xIxx/PrTJ5T1WjN7wt3nF3ouUovCzP4K+DYwnWCcYibwPHBwWRUNMndfACwAmD9/fln53dpcz+2XHVdmBaX/eWJl/EVTzh9BVs6Byj5WOccZms+uHI3JBONGNbCto5u12zqZNraJ0Y1Jnl65GTM4fK9xPP/GNp56fRPHz5nIhNGN3Lp4JWMa6zhn3gx+s2Q1D720no8cN5OWpiT/cMdS3j53Mv/4vgO5duFzPLJ8Pd897wgy7lz1y6c5csY4fviReXzox4/z0trt/OwTx7BmSzufu+0Z/vKwaXzmlLms2dzOk69v4tiwu6mQ1qZ6/tzR++K7VJ8WRXaMYlRDT9fTuOZ6Lj15Xx5Zvp6P/XQxAF898xCuvv0ZHluxkWljm7j2rw7m4pva2LCji++cdzg3P/oaT76+mWNmT+DKd+/HeQseo6Euwe2fOZ5v/e+f+f2yt/jsqftxykFT+cAPH+XIvcdz08eP5ocPvMyvl6zmu+cdyZSWRr581zKOmT2Bj584m4VPr+bFt7ZzwXGzALjjqVUcuudYjttnIg8tX8+G7V2cvP9kdnameeL1jcycOJq9J4zixbe20dGdZsb4UXSnnU07u2hIJmioS+AEV6O7O+lMcL/cr1E5/6bK/c6W+rLs/8/BFvVdvwYcC9zr7kea2TsJ1s3eHasJurSy9gq3FdpnlZklgbEEA+mDrr4uwbyZWotJdjW5pZF9Jo/JPX7XgT1nh584t5ET5/b8ZX/1ew7I3f/kO/blk+/YN/f4vs+enPsl858XHsXW9hRjR9UDcOieY5k+rplRDUluueRYVqzfwVGzJuDujGuu55jZQTB84/2H4kCiSN9Ja3P9Ll1PmT5jFN25MYrsYHY7E0Y3sOe4Zn59+Qn8wx1LgeAv1LOO3JNv/PbPnH7EdN65/2QmtzSSyTjvPXQa40c1cPFNbXzutP2ZN3MCnz/tAKaPa2LOlDF8+9zDeWT5Bk49aCqJhPHQ597JmMYkyboEV526P1edun+uvusv6PlD9qwj9+pV+6dOnpO7/879886hGQN7TxyVezhpTGO/n4nsnqhB0e3uG8wsYWYJd7/PzL67m8deDMwNpytfDZwP/HWffRYCFwCPAucAf6zE+ITIUMj/S9TMciEBMGdKzxybE8c0MjH8pWdmnHrwHrnnigVEVmtzPds7U2Qynts/aFEk8q6jyI5RhC2KLR25vu3Wpnq+/9dvy73fOfP2YvGrG/nwMTNJ1iX47nlHANCYrOPk/afwzLV/kTt76rKTe4Kxpame0w7pqX3cqIYBa5fhKWpQbDazMcBDwM/NbC3hTLLlcveUmV0O/A6oI1gc6Tkz+wrQ5u4LgZ8AN5vZcmAjQZiISBGtTUncYVtHT4slnXHq8qbw6Ol6Cn7Bd6UyTBxd+Bf5xDGNXH/BUbnHJ/QZG8mGhNSuokFhZn8LLALOBHYCfwt8iKAL6Cu7e3B3vxu4u8+2L+bd7wA+sLvHEYmTsc1BOGzt6M4FRSqT6T0pYKb3GAXAhH6CQmSgFsVewHeBA4ClwCMEwXGXu2+scG0iUobWMCi2tHfnBgEz4RQedX26nkbntQYmqo9f+lE0KNz9swBm1gDMB44HPgYsMLPN7t73AjkRqbJciyJvQLs7k6Ehmdy166mx51dAf11PIlGn8GgGWgm6nMYCa4DHK1WUiJSvtamn6ymrvStNc0MdFl54ljvrKa9Foa4n6c9AYxQLCK6V2EYQDIuAb7v7piGoTUTKkB2XyD9Fdntnipaw9ZAw22WuJ4CJYxQUUthALYq9gUaCiQBXA6uA4msuikhVtTYFgbC1vWdOpR2dqdy8TgnrPXts1sTRGqOQwgYaozjNgpO/DyYYn7gKOMTMNgKPuvu1Q1CjiJRgdEOShPXuetreKygsN0PsqEZ1PcnABryOIrzA7Vkz2wxsCW/vI5jUT0EhMswkEtbr6uzOVDCPUktTga6n+rrcAPf4vAsARfINNEZxBUFL4nigm2CMYhFwA8HpsiIyDLU21efOetrRGaxjkR24TljPmtnJRILGZILm+jqSdaUuTyNxMVCLYhbwK+Dv3P2NypcjIoNhbF6LYntHMFZRqOspkQim4lC3kxQz0BjFlUNViIgMntbmJFvDgNjeGfw31/WUsNyV2dkWhS62k2LU1hSpQWOb87qeuvq2KMi1KOoSwXxPkxUUUkRlJi8XkapqbcrreurcteupMxcUCb5+1qFMalFQSP8UFCI1qLW5Pnd6bHaMInvBneWd9ZRMWL8r5YlkqetJpAaNba6noztDZyrNjj4tirpEzxQeUda3EFFQiNSg/KuzC3U9ZZf1TiooJAIFhUgNas1bkyIbFGPygiIrMVQLkMuIpqAQqUH5a1Ls6EzRXF+XWy87PxvUopAoFBQiNWh8uD715p1dbO9M57qdgFxgANTVKShkYAoKkRo0OTzdde3WzmCK8aaeoMjvbqpT15NEoKAQqUHZC+jWbusMpxjvmSU2Pxvq1PUkESgoRGpQQzLB+FH1rN3WEUwx3tBPi0JBIREoKERq1JSWpqDrqaN311Odup6kRAoKkRo1pbWRt7Z1sqMr1WswO5sNZrrgTqJRUIjUqCktTazb2tFrGVTo6XrSqbESleZ6EqlRU1obWbe9EzPLzfMEwRoUoPEJiU5BIVKjprQ00p12wHtfRxG2KDQ+IVGp60mkRk1pacrd7z1GEQaFWhQSkYJCpEZNae1ZY6JX11OYDwoKiUpBIVKjpuQtRlRoMLsuoX/+Eo2+KSI1Kr/raUyBKTx01pNEpaAQqVHNDXW5C+3G5E3hobOepFQKCpEalu1+Ktz1pKCQaBQUIjUs2/00RkEhu0FBIVLDsmc+jSkwhYeCQqJSUIjUsEJdT9mA0AV3EpWuzBapYefOn8HU1ibq63r+JlTXk5SqKkFhZhOAW4FZwKvAue6+qcB+aWBp+PB1dz99qGoUqQVzp7Ywd2pLr23ZfEhqGVSJqFpdT1cDf3D3ucAfwseFtLv7EeFNISEyCLJTeCTU9SQRVSsozgBuDO/fCJxZpTpEYqdOF9xJiaoVFFPd/Y3w/pvA1H72azKzNjN7zMz6DRMzuyTcr23dunWDXqxILclecKdFiySqio1RmNm9wB4FnvpC/gN3dzPzft5mpruvNrN9gD+a2VJ3f7nvTu6+AFgAMH/+/P7eS0To6XpSi0KiqlhQuPsp/T1nZm+Z2TR3f8PMpgFr+3mP1eF/V5jZ/cCRwC5BISLR6awnKVW1up4WAheE9y8AftN3BzMbb2aN4f1JwAnAsiGrUKRG1emCOylRtYLim8C7zewl4JTwMWY238yuD/c5EGgzs6eB+4BvuruCQmQ3afZYKVVVrqNw9w3AuwpsbwMuCu8vAg4d4tJEap5Oj5VSaQoPkZjRBXdSKgWFSMxkxybUopCoFBQiMaPTY6VUCgqRmEnkznrSP3+JRt8UkZjpuY6iyoXIiKGvikjMqEUhpdI3RSRmsnM8qUUhUemrIhIzPRfc6Z+/RKNvikjMJDSFh5RIQSESM5oUUEqloBCJmZ4xCgWFRKOgEImZXNeTrsyWiBQUIjGjricplYJCJGY0hYeUSkEhEjPZLietmS1RKShEYiY3zbiCQiJSUIjEjM56klIpKERixnTBnZRIQSESM3U660lKpKAQiRmdHiulUlCIxIxpMFtKpKAQiZlsi0JrZktUCgqRmMl2OSXrFBQSjYJCJGayPU5qUUhUCgqRmDEtXCQl0jdFJGZ01pOUSkEhEjPZtbIVFBKVgkIkZjR7rJRKQSESMwnNHislUlCIxIxmj5VSKShEYkYX3EmpFBQiMZPQBXdSIgWFSMwkNM24lEhBIRIzueso1PUkESkoRGJGF9xJqaoSFGb2ATN7zswyZja/yH6nmdkLZrbczK4eyhpFatVx+07kk+/Yh/33aKl2KTJCVKtF8SzwfuDB/nYwszrgOuA9wEHAB83soKEpT6R2jW2u55r3HEh9nToUJJpkNQ7q7s9DzxWi/TgaWO7uK8J9bwHOAJZVvEAREckZzn9S7AmszHu8Kty2CzO7xMzazKxt3bp1Q1KciEhcVKxFYWb3AnsUeOoL7v6bwTyWuy8AFgDMnz/fB/O9RUTirmJB4e6n7OZbrAZm5D3eK9wmIiJDaDh3PS0G5prZbDNrAM4HFla5JhGR2KnW6bFnmdkq4Djgf8zsd+H26WZ2N4C7p4DLgd8BzwO/dPfnqlGviEicVeuspzuBOwtsXwO8N+/x3cDdQ1iaiIj0MZy7nkREZBgw99o6ScjM1gGv7cZbTALWD1I5Q0H1VpbqrSzVW1ml1DvT3ScXeqLmgmJ3mVmbu/c7rchwo3orS/VWluqtrMGqV11PIiJSlIJCRESKUlDsauVvkIkAAAV4SURBVEG1CyiR6q0s1VtZqreyBqVejVGIiEhRalGIiEhRCgoRESkqlkEx0Mp5ZtZoZreGzz9uZrOGvspe9QxU79vN7EkzS5nZOdWosU89A9V7pZktM7NnzOwPZjazGnXm1TNQvZea2VIzW2JmDw+HBbSirv5oZmebmRdbSbLSIny+F5rZuvDzXWJmF1Wjzrx6Bvxszezc8Dv8nJn911DXWKCegT7j7+R9vi+a2eaSDuDusboBdcDLwD5AA/A0cFCffT4F/DC8fz5w6zCvdxZwGHATcM4I+HzfCYwK7182Aj7f1rz7pwP/O9w/43C/FoJVJB8D5g/XWoELge9X8zMtsd65wFPA+PDxlOFec5/9/wa4oZRjxLFFkVs5z927gOzKefnOAG4M798GvMsGWI6vggas191fdfdngEw1CuwjSr33ufvO8OFjBFPIV0uUerfmPRwNVPsMkCjfYYCvAt8COoayuD6i1jpcRKn3YuA6d98E4O5rh7jGvkr9jD8I/KKUA8QxKKKsnJfbx4NZbLcAE4ekul1FXulvmCi13k8Av61oRcVFqtfMPm1mLwP/D7hiiGrrz4A1m9nbgBnu/j9DWVgBUb8PZ4ddkbeZ2YwCzw+VKPXuB+xnZo+Y2WNmdtqQVVdYKauBzgRmA38s5QBxDAoZJszsw8B84J+rXctA3P06d98X+Dzwf6tdTzFmlgC+DVxV7VoiuguY5e6HAb+npzU/XCUJup9OJvjr/MdmNq6qFUV3PnCbu6dLeVEcgyLKynm5fcwsCYwFNgxJdbsaaSv9RarXzE4BvgCc7u6dQ1RbIaV+vrcAZ1a0ooENVHMLcAhwv5m9ChwLLKzSgPaAn6+7b8j7DlwPzBui2gqJ8n1YBSx09253fwV4kSA4qqWU7/D5lNjtBMRyMDsJrCBofmUHfg7us8+n6T2Y/cvhXG/evj+l+oPZUT7fIwkG3+aOkO/D3Lz7fwW0Dfea++x/P9UbzI7y+U7Lu38W8Nhw/myB04Abw/uTCLp9Jg7nmsP9DgBeJbzQuqRjVOuHq+aNYHGkF8NfVl8It32F4K9bgCbgV8By4E/APsO83qMI/srZQdDyeW6Y13sv8BawJLwtHOb1/hvwXFjrfcV+KQ+XmvvsW7WgiPj5fiP8fJ8OP98DhvNnCxhB194yYClw/kj4PgBfAr5ZzvtrCg8RESkqjmMUIiJSAgWFiIgUpaAQEZGiFBQiIlKUgkJERIpSUIgMwMwm5s28+aaZrc57vGgQj3OmmX2xyPOHmtlPB+t4IlHp9FiREpjZl4Dt7v4vFXjvRQTnva8vss+9wMfd/fXBPr5If9SiENkNZrY9/O/JZvaAmf3GzFaY2TfN7ENm9qdwLYt9w/0mm9ntZrY4vJ0Qbt8P6MyGhJl9wMyeNbOnzezBvEPeRTBbgMiQUVCIDJ7DgUuBA4GPAPu5+9EE8xf9TbjPvwHfcfejgLPD5wBOAJ7Me68vAn/h7ocTrIGR1QacVLGfQKSAZLULEKkhi939DYBwSvJ7wu1LCRZrAjgFOChveZNWMxsDTAPW5b3XI8BPzeyXwB1529cC0ytTvkhhCgqRwZM/C24m73GGnn9rCeBYd++1mJCZtRPMUgyAu19qZscAfwk8YWbz3H0DwTxk7RWqX6QgdT2JDK176OmGwsyOCO8+D8zJ276vuz/u7l8kaGlkp5HeD3h2iGoVARQUIkPtCmB+uJrbMoIxDQjWtj4yb8ndfw4HwZ8FFhHMrApBF1a1V62TmNHpsSLDhJn9G3CXu9/bz/ONwAPAiR4s0SsyJNSiEBk+vg6MKvL83sDVCgkZampRiIhIUWpRiIhIUQoKEREpSkEhIiJFKShERKQoBYWIiBT1/wFuGNzaJw+4CAAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] @@ -205,7 +244,10 @@ } ], "source": [ - "plt.plot(times_L1,whiten_signal_L1)" + "plt.plot(times_L1,whitened_signal_L1)\n", + "plt.title('Whitened signal')\n", + "plt.xlabel('Time(s)')\n", + "plt.ylabel('Waveform')" ] } ],