{ "cells": [ { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAEHCAYAAAC+1b08AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nO3deXhU5fn/8fedhLAISZAoVqMsElCKFGoqWlRUpGJtoVapS7VgUVrr0qpVca3a1qL+6lZxoQoosggoSC1+UYsWbVXEuhQRK6tEZIlC2IRsz++PJ4EQsswkMzlzZj6v6zrXZGbOOXOfLPfceeZZzDmHiIiEX1rQAYiISGwooYuIJAkldBGRJKGELiKSJJTQRUSShBK6iEiSyGhoBzMbD/wA2OCc61XL8z8Frq+8uw241Dn3QUPnzc3NdZ07d44u2kplZWVkZDQYelLRNacGXXNqaMo1v/vuu0XOuQNqey6SM04EHgKequP5lcAA59wmMzsdGAf0a+iknTt3ZtGiRRG8/L6KiorIzc1t1LFhpWtODbrm1NCUazaz1XU912BCd84tMLPO9Tz/72p33wLyoglORERiI9Zt6COBF2N8ThERiUDMGq7M7GR8Qj++nn1GAaMA8vLyKCoqatRrFRcXN+q4MNM1pwZdc2qI1zXHJKGbWW/gceB059yXde3nnBuHb2OnoKDANaXdLNXa3EDXnCp0zakhHtfc5CYXMzsMeA640Dn3v6aHJCIijRFJt8WpwElArpkVAr8DWgA45x4FbgU6AA+bGUCZc64gXgGLiEjtIunlcl4Dz18MXByziEREpFFSqze/NJ8tW2DpUli7Fr74AoqLobwcKiqgdWvYf3+/dekC3brBfvsFHbFI6CmhS2x8/jm89BK8/DIsXAjLl0d3/CGHwLe/DcceC/36+U1EoqKELo23dStMnw5PPgmvv+4fO+ggOP54uOgi6NUL8vLg4IMhJwfS0yEtDXbsgE2boKjIJ/5ly3w1/8478Le/+fNkZpJ17LEwZAiccQYccURw1ykSEkroEr0NG+CBB2DsWN+U0r07/OEP8MMfwlFHgf9wvG5ZWX7r1AmOPnrv5zZt8hX+K6+Q9sIL8Nvf+u2oo+C88+Ccc6Br1/hdm0iIabZFidyOHXD77dC5M/zpT3DqqfCvf/nq+qaboHfvhpN5Q9q3h9NOg3vuYfPrr8Nnn8Ff/gLt2sGNN8Lhh8N3vwsTJvh4RGQ3JXSJzOzZvtnjttvgBz+AJUtg5kyfXJuaxOtz6KFw+eX+jWPVKhgzxlfxP/+5b8q54gpYvDh+ry8SIkroUr/iYhg+HM480/dKWbDAt5sH0abdqRNcf71/M1mwwL+x/PWvvjlm4ECYO9f3ohFJUUroUrd33/XNKJMnw623+g8tTzgh6Kj8fwQnnABPP+1719x1F3zyif/wtFcvePxx2Lkz6ChFmp0SutRu0iTfWwV8c8ftt0OLFsHGVJsOHeC662DFCh9zy5ZwySVw2GFw553+PwyRFKGELntzDm6+GX72M98nfNGicPQJz8yECy6A//wH5s+HggL/QW2nTv6/iy/rnDNOJGkoocse5eXwy1/CH//oq9yXXoIDal3pKnGZwckn+/b0RYt82/rvf+8T+3XXwbp1QUcoEjdK6OKVlcH558O4cb6yfeyxxGxiicbRR8Ozz8J//wtDh8Kf/+ynGrjySlizJujoRGJOCV18ZT58uO+9cs89fpBQPLsiNrdevfwHu0uX+jetRx7x/dl/8QtYuTLo6ERiRgk91Tnnm1mmTPFNLb/9bdARxU9+PjzxBHz6KYwcCRMn+scuugj+p6n8JfyU0FPdLbf4bn433ui3VNC5s6/SV6zwg5amTYMjj/TV+0cfBR2dSKMpoaeyp57yVfnIkb6ZJdUccgjcf78fgXrNNTBnjm+eOftseP/9oKMTiZoSeqp6/XW4+GLfI+Thh5OrzTxaHTvC3Xf7xH7TTX4K4L59/UyPCxcGHZ1IxJTQU9Fnn/mh/F27+l4gmZlBR5QYcnP9fyqrV8Mdd8Abb/g++Ked5r8WSXBK6KmmpMRPQVtS4psY2rcPOqLEk5PjP1tYvdpPBvbee36qgZNP9oOWnAs6QpFaKaGnmhtugLfe8r09uncPOprE1q6dnwxs5Uq4914/X8zAgXDMMTB1KpSWBh2hyF6U0FPJ88/7xHTZZTBsWNDRhMd++8FVV/leMY8+6ldqOv9832R1zz2weXPQEYoASuip44sv/BziRx/tR0xK9Fq18oORliyBF17wfdivu84vs/frX0e/jqpIjCmhpwLnYNQov8LP5Ml+RkJpvLQ0P1Xv/Pl+MrAf/9j3FOrWDb73Pf9Bs5pjJABK6Klg4kRfUf7pT9CjR9DRJJe+fX1//tWr/RTDH3/s+7EfdpiftXL16qAjlBSihJ7sVq/2zQEDBvhJqSQ+Dj7YT9O7ahX87W++aevOO/1kYIMGwZNP+rZ3kThqMKGb2Xgz22BmtS7caN6DZrbMzD40s2/HPkxpFOfg0kv9smwTJvimAomv9HS/NN4LL/jkfsstvm19xAg/gOm88/xzapKROIjkL3wiMLie508H8iu3UcAjTQ9LYuK55+DFF/184F26BB1N6jnsMN8Ms3w5/PvffhKwl1+GH/4QDjrIJ/nZs/1nGyIx0GBCd84tAL6qZ5ehwFPOewvIMbNvxCpAaaStW30Ty7e+BVdcEXQ0qc0MjjsOxo6FtWv9gK7vf993Iz3zTD9C9Uc/ggkTSPvii6CjlVgrL4ctW/ziKitXwscfY+vXx+WlMmJwjkOA6qsFFFY+FpffTOccl1xyCWtSbIGCsrIyMjKi+HGtWQMbNkBWll9KLoSivuaw6doVtm3z/djnzvUJHnz3yKwsP7CpXTvfjJPEEvbn7Jxf+KW83N9WbdXvV1TsvZWX7/tYLSOLzxswgGteey3mIcfiu1jbrE61jo02s1H4Zhny8vIoKiqK+sV27tzJ7Nmz6datG11SqBkhml9627qVFhs2UJ6XR/mhh8Y5svhJ2D/0WOrQwS+Ph/+5uY0bydi8mbSiIv+GDLh27ajIzsbl5OCys3Ft2gQZcczF/efsHJSWYqWle93W9thez5WX131OM1xGBmRk+M+m0tNxLVr4N+P0dH+/8pb0dEhL23M/LQ3r27dR+a8hsfguFgLVs0YesLa2HZ1z44BxAAUFBS43NzfqF9u5cycAI0eOZPTo0VEfH1ZFRUVE9P1yDk49FfbfHz78MNRztUR8zUlk9zXv2uWnaPjHP3z7+8KFUFjod8rN9d0ljzrKb717+/ncW7cONvhGivjnXFHh/5v56iu/6Hekt8XFdZ8zLc3/jXTo4P9mIr1t165JM5TG63c7Fgl9DnC5mU0D+gHFzrm4NQQ6TYxUv7//3Q94efDBUCfzlNeype9qOmCAv19e7keovvWW3z74wA9mqixwSEvzC3d07eqX1+va1W+dOvkPYDt2TJxZNUtKYNOm3VuLVat8st60qf7kvGlT3ROjmflJ1aqSbm6uH3NRWzKu/nV2dlL1/mowoZvZVOAkINfMCoHfAS0AnHOPAnOB7wPLgB3ARfEKtkZczfEy4VJa6peQ697dLysnySM9fU9Ffskl/rHycli2zC+C/d//+mX0VqzwI1Vr+3e+Q4c9yT0nx7fTZ2X5pJaVBW3b+oXBa9tqtieXl+/5+uuvYfv2PduOHXvf37JlrwTO11/vFVZ2zTizsvZOwl261J6Mq9/m5CT9Zw2RaDChO+fOa+B5B1wWs4gaoAq9Ho895mcEnDPH/xFKcktP91Vojx5+dGp1xcW+R8Xq1bB+vd/WrduzffKJT7RVW6z+rlq29JOZVd+ysnyM7dv7xNu+/Z4tJ4fNZuR06bLnMf3uNlqSf+KUQjZvhttug1NO8QNbJLVlZ0OfPn5rSEWFr6q3bt39geA+W+UHf2Rk7Pmgr+rr1q194m7TplFVcllRkW8ikSYLXUKvqtDV5FLDPff4dsY//zm1l5OT6KWl+eaWtm2DjkSaKHk+DUhlGzfCAw/4lYgiqchEJCmFLqGrQq/FXXf5D5puuy3oSEQkQKFL6FLD2rV+SPmFF2pqXJEUF7qErgq9hj/9yXcdu/XWoCMRkYCFLqFLNZ99BuPG+aXlunYNOhoRCVhoE7oqdODuu33/4ZtuCjoSEUkAoUvoGlhUaf16eOIJGD7cz7stIikvdAm9SspX6Pff7+fEuO66oCMRkQQRuoSuCh0/KnTsWD/cOz8/6GhEJEGELqFXSekK/eGH/TDtG24IOhIRSSChS+gpX6Hv2AH33Qenn65RoSKyl9Al9JQ3YYKfGlXVuYjUELqEntIDiyoq/Ieh/frBCScEHY2IJJjQJfSU9ve/+wUNrroq6EhEJAGFLqGndIV+331w6KFw1llBRyIiCSh0CT1VpS9eDK++Cpdf7hcWEBGpIbQJPdUq9NaPPeZXhKlaT1JEpIbQJfSU7La4fj0tn3sOLrrIr7koIlKL0CX0KilVoT/yCFZSAldeGXQkIpLAQpfQU65CLymBRx+lZNAg6N496GhEJIGFLqFXSZkK/fnnYf16vr7ooqAjEZEEF7qEnnIV+iOPQOfOlJ5yStCRiEiCC11CTylLl/quir/4BaSnBx2NiCS4iBK6mQ02s0/MbJmZja7l+cPM7FUze8/MPjSz78c+VC+lBhY9+ii0aOGXmBMRaUCDCd3M0oGxwOlAT+A8M+tZY7ebgenOub7AucDDsQ405ezYAU8+6UeFHnhg0NGISAhEUqEfAyxzzq1wzpUA04ChNfZxQFbl19nA2tiFWOOFUqVCf+YZv5DFpZcGHYmIhEQkY8gPAdZUu18I9Kuxz23AS2Z2BbAfcGpMoktljzwCPXtqVkURiVgkCb22UrhmV5PzgInOuT+b2XHAJDPr5Zyr2OtEZqOAUQB5eXkUFRVFHXDVMdu3b2/U8WGQ/sEHtH/nHbbdeSc7v/wSgOLi4oCjan665tSga46dSBJ6IXBotft57NukMhIYDOCce9PMWgG5wIbqOznnxgHjAAoKClxubm7UAZeVlQHQtm1bGnN8KMycCa1b0/bSS2mbk7P74aS93nromlODrjk2ImlDfwfIN7MuZpaJ/9BzTo19PgMGApjZkUArYGMsA60padvQd+yAKVP8AtDVkrmISEMaTOjOuTLgcmAe8DG+N8tHZnaHmQ2p3O0a4BIz+wCYCoxwcRoBlPQDi2bNgi1b1FVRRKIW0cTazrm5wNwaj91a7eslQP/Yhla/pK3Qx4+Hrl3hxBODjkREQiZ0I0WTukJfuRLmz/fT5KaF7kcjIgELbdZIygp94kQwg+HDg45EREIodAk9aSv08nKYMAG+9z2/bqiISJRCl9CT1vz5sGaNb24REWmE0CX0pB36P368X15uaM1ZFUREIhO6hJ6UNm3y3RV/+lNo1SroaEQkpEKX0JOyQp82DXbtUnOLiDRJ6BJ6Upo0CXr1gr59g45EREIstAk9aSr05cvhzTfhggt8l0URkUYKXUJPum6Lkyf7RH7++UFHIiIhF7qEXiUpKnTnfHPLSSep77mINFnoEnpSVegLF8KyZb65RUSkiUKX0KskRYX+9NO+m+JZZwUdiYgkgdAl9KSp0EtLfXfFIUMgOzvoaEQkCYQuoSeNefOgqEjNLSISM6FL6EkzsOjpp6FDBzjttKAjEZEkEbqEnhS2bIHnn4dzz4XMzKCjEZEkEbqEnhQV+rPPws6dam4RkZgKXUJPClOnwuGHQ79+QUciIkkktAk9tBV6UZGf+/wnP9FQfxGJqdAl9NB3W5w1y69O9JOfBB2JiCSZ0CX0KqGt0KdPh27d4FvfCjoSEUkyoUvooa7QN26EV19Vc4uIxEXoEnqVUFboVc0tw4YFHYmIJKHQJfRQV+gzZqi5RUTiJqKEbmaDzewTM1tmZqPr2OcnZrbEzD4ysymxDbPW14v3S8SWmltEJM4yGtrBzNKBscAgoBB4x8zmOOeWVNsnH7gB6O+c22RmB8Yr4NBW6GpuEZE4i6RCPwZY5pxb4ZwrAaYBQ2vscwkw1jm3CcA5tyG2YSaBGTMgP1/NLSISN5Ek9EOANdXuF1Y+Vl13oLuZ/cvM3jKzwbEKsKZQDv3fuNEPJho2TM0tIhI3DTa5ALVloJrtHhlAPnASkAe8bma9nHOb9zqR2ShgFEBeXh5FRUVRB7xp0yYAtm7d2qjjg9By0iTaVVSwadAgyhsZc3FxcYyjSny65tSga46dSBJ6IVB9wcs8YG0t+7zlnCsFVprZJ/gE/071nZxz44BxAAUFBS43NzfqgKuSeFZWFo05PhBz50J+Pu0HDGhShR6a640hXXNq0DXHRiRNLu8A+WbWxcwygXOBOTX2mQ2cDGBmufgmmBWxDDS01LtFRJpJgwndOVcGXA7MAz4GpjvnPjKzO8xsSOVu84AvzWwJ8CpwrXPuy3gFDSFqQ581Cyoq1LtFROIukiYXnHNzgbk1Hru12tcOuLpyi6vQdVucPt33bundO+hIRCTJhW6kaJVQVOhqbhGRZhS6hB6qCv2559TcIiLNJnQJvUooKvQZM6B7dzW3iEizCF1CD02FvmGDb27RYCIRaSahS+hVEr5Cr+rdopWJRKSZhC6hh6ZCr2puOeqooCMRkRQRuoQeCmpuEZEAhC6hh2JyLjW3iEgAQpfQQ2H6dDW3iEizC21CT9gKfcMGeO01DSYSkWYXuoSe8B+KajCRiAQkorlcElHCVugzZkCPHoE1t1RUwLZtUFwMW7fCzp17tq+/3vt+1WMlJVBW5rfy8j1fN/RYebl/Pef23aJ5vK59S0qyyYjwNzSa9/loa4J4nbu2fcvKcmq95kSvY5qivDyH9PSgo2hew4a14uabY3/e0CX0hK7Qq5pbbrwxZs0tW7bAqlXwySct2LUL1q/3L1O1bdrk9yku9rdbtzbtjz8jY+8tPb3u+2lpfjPbe6vtsboeT0+ve9/SUkfLlpHHHs23PNofT7zOXXPfkpJyMjNr/7NM1Bqmqeq75mS1//7xyWOh/S4mZIXeyOaW7dthyRJYvNjfrljhk/jKlT5he9m798/MhI4d4YADYP/94RvfgOxsyMrac5uVBe3aQevWfmvVau+t+mOZmdCihU+iiaSoaEvKLXxQVLSV3Nwo3sWSQGpe8y6gXczPG7qEntAVegTNLaWl8MEH8O9/+23RIp/Aqy6rZUvo0gU6d4Z+/fxt587Qps1munfPoWNHn6wT8f1MRIIVuoReJeEq9HqaWz79FF580W///KdvtwY49FA45hj42c/8e0CvXtC1K7W2JxYVlZFixaqIRCl0CT1hK/Sq5pbKwUTLl8OUKX5butTv0qMHXHIJ9O8Pxx3nE7qISKyELqEnrOnTcd17MGNJL/7yK3jjDf/wgAHwq1/B978Phx8ebIgiktxCl9ATcej/9pUbaP3aP3mw7U1cda7RtSuMGQPnn68qXESaT+gSeiLZtQv++lf47MbnuNtV8FHPYcy5yVfjqdavVkSCF9qEHnSF/o9/+KaU//0P3s2Zzo5OR/DXN3tB4vzjICIpJsF6Hjcs6A9FN2+GCy+EU0/1IyVfmbyevlv+SZufaapcEQlW6BJ6lSAq9H//G/r0gWnT4JZb4L//hYHFz2Gau0VEEkDoEnoQFbpzcO+9cOKJfjTlG2/AHXf40ZbMmAFHHOE7kYuIBCh0Cb1Kc1XopaVw6aVwzTUwdCi8954fwQn4iVX++U+tTCQiCSGihG5mg83sEzNbZmaj69nvbDNzZlYQuxD31pwV+o4d8MMfwmOPwejRvhjPzq62Q43BRCIiQWqwl4uZpQNjgUFAIfCOmc1xzi2psV874Erg7XgEWktccT3/9u0+mb/2mu+aePHFtexU1dzyzW/GNRYRkUhEUqEfAyxzzq1wzpUA04Chtez3e+BuYGcM49tHc1ToX38NP/iBb0156qk6knlVc4tWJhKRBBFJQj8EWFPtfmHlY7uZWV/gUOfcCzGMrV7xqtDLy323xKpkfsEFdeyolYlEJMFEMrCotsy5u0w2szTgPmBEgycyGwWMAsjLy6OoqCiyKKvZVDlBeHFxcaOOb8gtt+zHs8+25o47tnHaaTup6yWyJk8mLT+fzR07UudOMVRcXBz310g0uubUoGuOnUgSeiFQfUaSPGBttfvtgF7Aa5VV80HAHDMb4pxbVP1EzrlxwDiAgoIC15jFC3JycgDIzs6O+eIHEybAo4/ClVfCzTe3xaxt7TuuWwdvvgk330zuAQfENIb6pNpiD6BrThW65tiIpMnlHSDfzLqYWSZwLjCn6knnXLFzLtc519k51xl4C9gnmcdarJtcPvjAD+UfOND3Oa/39GpuEZEE1GBCd86VAZcD84CPgenOuY/M7A4zGxLvAGuJJ+bn3LLF5+b27WHy5Agm1poxA448Ur1bRCShRDQ5l3NuLjC3xmO31rHvSU0Pq2GxrNCvucYvSPHaa36tznqtWwcLFsDNN6t3i4gklNCNFI11hf7SS/D443DttXDCCREc8OyzGkwkIgkpdAm9Siwq9C1b/JJwRxwBt90W4UHTp/umFjW3iEiCCV1Cj2WFfuutUFjoe7e0ahXBAWvXwuuvqzoXkYQUuoRepakV+scfw0MP+Qr92GMjPGjmTD/1ohK6iCSg0CX0WFTozsFVV0HbtvD730dx4DPPQO/evo1GRCTBhC6hV2lKhT53LsybB7/7HUQ8LmjNGr/CxTnnNPp1RUTiKXQJvakVekUF3HAD5OfDZZdFceCMGf5Wg4lEJEGFdpHoxpo1yy8d9/TTkJkZxYHTp0Pfvv6dQEQkAYW2Qm9Mk0tFBdx+O3TvDueeG8WBq1bB22+ruUVEElpKVejVq/MGh/dXp+YWEQmB0FXoVaKt0J2DP/4RevSIsjoH37vlO9+Brl2jPFBEpPmELqE39kPRBQv8As/XXBNldb58Obz7rvqei0jCC11CrxJthX7ffdChQz0rENVl+nR/q+YWEUlwoUvojanQly2DOXPgl7+E1q2jPHj6dD+UtFOnqF9XRKQ5hS6hV4mmQh87FjIy/AIWUfn4Y3j/ffVuEZFQCF1Cj7ZC37XLL/b8ox/BwQdH+WKTJ0NaWiM+RRURaX6hS+hVIq3QZ8+Gr77yk3BFxTmYMsWvSXfQQdEHKCLSzEKX0KOt0B9/3Dd/DxwY5Qu9+SasXAk//WmUB4qIBCN0Cb1KJBX6ypXwyiswcqRvOYnK5Ml+kvQzz2xcgCIizSx0CT2aCn3SJL/s54gRUb5Iaanv3TJkCGRlRXmwiEgwknbov3MwdSqceCIcemiUB7/0EhQVqblFJAZKSkpYvnw5O3bsqHOf1atXN2NEiaGha27Tpg2HH344mVHMIhjahN5Qk8uHH8LSpfDrXzfi5JMnw/77w+DBjQtORHZbvnw5OTk59OjRg7So2z5TU0VFBevWrePDDz+kS5cudOjQIaLjQvfdjbTJZepU3/f87LOjfIFt2+D55/3I0Kjm1xWR2uzYsYOOHTsqmUchLS2Ngyp7182YMYPNmzdHdlw8g4qn+ip052DaNBg0CHJzozzx7NmwY4eaW0RiSMk8emlpaZgZO3fuZO3atZEdE+eYYi6SCn3hQli9upEDPCdNgsMOg/79G3GwiCSazZs38/DDD+++v3btWs6O+l/32rVt27ZRx61cuZJ+/fqRn5/POeecQ0lJSZ37OucoLy+P6LwRJXQzG2xmn5jZMjMbXcvzV5vZEjP70Mz+YWZxn/ikvgr9+ef9jIpDhkR50jVr4OWXYfjwRvRzFJFEVDOhH3zwwcycOTPAiOD666/nqquu4tNPP6V9+/Y88cQTMTlvg1nLzNKBscDpQE/gPDPrWWO394AC51xvYCZwd0yiq0UkFfqcOb53S/v2UZ78qad8e03U/RxFJFGNHj2a5cuX06dPH6699lpWrVpFr169AJg4cSJDhw5l8ODB9OjRg9tvv333cffeey+9evWiV69e3H///fW+RlFREccddxx///vfG4zHOcf8+fN3/5cwfPhwZs+e3YQr3COSXi7HAMuccysAzGwaMBRYUi3AV6vt/xYQ7SS1UaurQl++HD76qJFD/SdOhAEDtJCFSJz85jd+vrtY6tMH6su3Y8aMYfHixbxf+cKrVq3a6/mFCxeyePFi2rRpw3e+8x3OOOMMzIwJEybw9ttv45yjX79+DBgwgL59++5z/vXr1zNkyBD+8Ic/MGjQILZu3coJJ5xQayxTpkzhwAMPJCcnh4wMn37z8vL4/PPPG3fxNUSS0A8B1lS7Xwj0q2f/kcCLTQmqPg1V6HPm+Nuom1v+9S8/z+7NNzcuMBEJpUGDBu3uFvjjH/+YN954AzPjzDPPZL/99tv9+Ouvv75PQi8tLWXgwIGMHTuWAQMGANCuXbvdbx612bhx4z6PNWaN5NpEktBre6Vas6qZXQAUAAPqeH4UMAr8u1JRUVGEYe5RXFy8+7a242fOzKZnT6Ndu81Ec/q2jzxC5n778dVJJxHVgc2k6rpTia45+TTQchGImsnUzCLuHp2RkcHRRx/NvHnzdif0hir0I488ks2bN1NWVkZGRgaFhYUcXM9UsLt27aoz3+0TTwQxFwLVx1rmAfv0oTGzU4GbgAHOuV21ncg5Nw4YB1BQUOByo+5TCFmVQ/FzcnKoefzmzfD223D99ezzXL22b/el/TnnkJvAC1k05vsVdrrm8At6FGi7du3YunVrnc+//PLLfPXVV7Ru3ZrZs2czfvx40tLSGDFiBKNHj8Y5x6xZs5g0adI+x5oZ48ePZ9iwYYwZM4bRo0c3WKEDnHzyycycOZNzzz2XJ598kqFDh9a5b8uWLcnOzo7o9yKSrhzvAPlm1sXMMoFzgTk1Lqov8BgwxDm3IYJzNllt/6LMnw/l5XD66VGebOZMP6BIH4aKJJ0OHTrQv39/evXqxbXXXrvP88cffzwXXnghffr04ayzzqKgoIBvf/vbjBgxgmOOOYZ+/fpx8cUX19p+DpCens60adN49dVX9+pNU5+77rqLe++9l27duvHll18ycuTIJl1jlQYrdOdcmZldDswD0oHxzrmPzOwOYJFzbg5wD9AWmFGZaDw2Eb4AAAp4SURBVD9zzkXbit1kL70E7dpBv/pa+Gvz+OPQrRscf3xc4hKRYE2ZMmWv+4sXL9799YEHHshDDz20zzFXX301V199db3n3bZtGwCZmZnMmzcv4ni6du3KwoULI94/UhHN5eKcmwvMrfHYrdW+PjXGcdUXSx2Pw7x5cMop0KJFFCdcvBjeeAPuucdPzSgiElKhHT1Ts8ll+XJYtcoP94/KY49By5ZqbhFJQSNGjKi1Og+r0CX0uir0l1/2t9/7XhQn27bNDyYaNqwRk76IiCSW0CX0KjUr9Jdegs6dfVN4xKZNgy1b4Je/jGlsIiJBCF1Cr61CLy+H116DU0+NohncOXjkETjqKPjud2Mao4hIEEKX0KtUr9AXL/Z90AfUOpypDosWwX/+46tzfRgqIkkgdAm9tgp9wQJ/e+KJUZzogQd8H8cL4j7tjIgEKBGnz33ooYfo1q0bZtaoEfN1CV1Cr1K9Ql+wADp18tOYR6SwEJ55Bi6+WItAiyS5RJw+t3///rzyyit0ivHI9NAl9JoVunM+oUdVnf/lL1BRAVdeGdvgRCThJNr0uQB9+/alc+fOjb6muoR+kej//Q82bIgioW/b5vuen3WW7xYjIs0ngPlzE2363J49ay4nETuhS+g1K/TXX/e3ESf0CROguBgaGNIrIqmhuafPjafQJfSaFiyAAw+E/PwIdi4thfvug+OOg2OPjXtsIlJDAs6f29zT58azQg9dG3qVqh/CW2/5buQR9TycPBlWroQbb4xvcCKSMCKdPvfrr79m9uzZ9O/fnxNPPJHZs2ezY8cOtm/fzqxZs2pN0lXT5y5dupQxY8bsfr3333+/1i2eyRxCmNCrv3N+9RV8+mmEsyuWlcEf/wh9+8IZZ8QvQBFJKIk4fe6DDz5IXl4ehYWF9O7dm4svvrhJ11gltE0uZkbV7JMRJfSpU/0Sc7NmaSCRSIpJtOlzr7zySq6MQy+7UFfob7/tc3NBQQMHlZXBH/4AvXs3YrFREZFwCHWF/vbb0LOnH/BZr8cf9/0bZ82CtNC9h4lInIwYMYIRSTR1duiyW1WF7hwsXBhBc8vWrfC738EJJ0A96/aJiIRdaCv0wkLjyy8jSOh33+1HHv3tb2o7F5GkFtoKveozjXoT+qpV8Oc/w3nnwTHHxD02EZEghS6hV1myxGjdGr75zTp2cA4uvRTS06Gyf6iISDILXUKvqtCXLvWdVjLqajR65hn4v//zfc8jnoZRRJKNps8NgU8+Mero5w9r18IVV/hmlssua9a4RCSxaPrcENi2zU+yto/ycrjwQtixAyZO9E0uIpKyNH1uAqs+sKjWCv3WW2H+fHjiCTjyyOYLTEQa9Jvf/CbmMxH26dOn3oSr6XNDwMw46qgaD44fD3fe6VciuuiiQOISkXBJuelzzWww8ACQDjzunBtT4/mWwFPA0cCXwDnOuVWxDdWrqtC7dIHWras9MWECXHIJDBoEDz+sPuciCaihposgpNT0uWaWDowFTgd6AueZWc2IRgKbnHPdgPuAu2IdaE1HHln5Qygr8yNBf/5zOPVUP7y/RYt4v7yIhISmz93bMcAy59wK51wJMA2oOYZ+KPBk5dczgYFW820vRoqL/TvnN/NL4PnnfU+WO+6A4cNhzhyo/BdJRAQ0fW5NhwBrqt0vBGqOz9y9j3OuzMyKgQ5A7DpYVip65T0AfvrQcXB/iV8X9Jln4Cc/ifVLiUiSSJXpcyNJ6LVV2jUbmCLZBzMbBYwCyMvLa1SH+uwDvsGAlodgQ0+h+EenUzpgAGRmQgw75yei4uLioENodrpmEdi1axfFxcUR5ctIEnohcGi1+3nA2jr2KTSzDCAb+KrmiZxz44BxAAUFBS43NzeCl9/bLx+8grNvPY/GHBt2uubUkGzXvHr16qBDqFMYps9t2bIl2dnZEf1eRNKG/g6Qb2ZdzCwTOBeYU2OfOcDwyq/PBua7SD8mFhGRmGiwQq9sE78cmIfvtjjeOfeRmd0BLHLOzQGeACaZ2TJ8ZX5uPIMWkXCpqKggTYvLRKWioiLqYyLqh+6cmwvMrfHYrdW+3gkMi/rVRSTptWnThvXr19OxY0cl9QhVVFSwbt06SktLozoutCNFRSQcDj/8cJYtW8bnn3++zyAeqVtpaSmrV6/GOUfrvUZR1k0JXUTiKjMzk549e/Lpp5/y4osv7vP8zp07adWqVQCRBSfSa66oqOCII46IeCIvJXQRaRb5+fl07NiRrVu37jW0ftOmTbRv3z7AyJpfpNecmZlJhw4dSI9w1lgldBFpNllZWWRlZe31WKtWrZKuq2ZD4nXN+oRCRCRJWFDdxc1sI9DYEQe5xGFagQSna04NuubU0JRr7uScO6C2JwJL6E1hZouccwVBx9GcdM2pQdecGuJ1zWpyERFJEkroIiJJIqwJfVzQAQRA15wadM2pIS7XHMo2dBER2VdYK3QREakhdAndzAab2SdmtszMRgcdT7yZ2aFm9qqZfWxmH5nZr4OOqTmYWbqZvWdmLwQdS3Mxsxwzm2lmSyt/3scFHVM8mdlVlb/Ti81sqpkl5fh/MxtvZhvMbHG1x/Y3s5fN7NPK25gMlQ1VQo9wwepkUwZc45w7EjgWuCwFrhng18DHQQfRzB4A/s85dwTwLZL4+s3sEOBKoMA51ws/NXeyTrs9ERhc47HRwD+cc/nAPyrvN1moEjqRLVidVJxzXzjn/lP59Vb8H/khwUYVX2aWB5wBPB50LM3FzLKAE/FrC+CcK3HObQ42qrjLAFpXrnLWhn1XQksKzrkF7LuC21DgycqvnwR+FIvXCltCr23B6qRObtWZWWegL/B2sJHE3f3AdUD0M/yHV1dgIzChsqnpcTPbL+ig4sU59znw/4DPgC+AYufcS8FG1aw6Oue+AF+0AQfG4qRhS+gRLUadjMysLfAs8Bvn3Jag44kXM/sBsME5927QsTSzDODbwCPOub7AdmL0b3giqmwzHgp0AQ4G9jOzC4KNKvzCltAjWbA66ZhZC3wyn+ycey7oeOKsPzDEzFbhm9ROMbOngw2pWRQChc65qv++ZuITfLI6FVjpnNvonCsFngO+G3BMzWm9mX0DoPJ2QyxOGraEHsmC1UnF/BIvTwAfO+fuDTqeeHPO3eCcy3POdcb/fOc755K+cnPOrQPWmFmPyocGAksCDCnePgOONbM2lb/jA0niD4FrMQcYXvn1cOD5WJw0VPOh17VgdcBhxVt/4ELgv2b2fuVjN1au8yrJ5QpgcmWxsgK4KOB44sY597aZzQT+g+/J9R5JOmLUzKYCJwG5ZlYI/A4YA0w3s5H4N7eYrMmskaIiIkkibE0uIiJSByV0EZEkoYQuIpIklNBFRJKEErqISJJQQheppnLGw18FHYdIYyihi+wtB1BCl1BSQhfZ2xjgcDN738zuCToYkWhoYJFINZUzWr5QOUe3SKioQhcRSRJK6CIiSUIJXWRvW4F2QQch0hhK6CLVOOe+BP5VuXCxPhSVUNGHoiIiSUIVuohIklBCFxFJEkroIiJJQgldRCRJKKGLiCQJJXQRkSShhC4ikiSU0EVEksT/B5RRjNWCw4PVAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "import matplotlib.pyplot as plt\n", "import numpy as np\n", "from scipy.signal import lsim\n", "import control as ctl\n", "from control import (TransferFunction)\n", "\n", "\n", "y0 = [0]\n", "\n", "# malha aberta\n", "L1 = TransferFunction(1, [1,1]) # tipo k=0\n", "L2 = TransferFunction(1, [1,1,0]) # tipo k=0\n", "\n", "plant_tf1 = ctl.feedback(L1, 1)\n", "plant_tf2 = ctl.feedback(L2, 1)\n", "\n", "t_ = np.linspace(0, 10, 1001)\n", "u_ = step_input = t_>0\n", "\n", "plt.figure(1)\n", "tout1, yout1 = ctl.step_response(plant_tf1, T=t_)\n", "tout2, yout2 = ctl.step_response(plant_tf2, T=t_)\n", "plt.plot(tout1, yout1, 'b', linewidth=1.5, label='tipo k=0')\n", "plt.plot(tout2, yout2, 'r', linewidth=1.5, label='tipo k=1')\n", "plt.plot(t_, u_, 'k', linewidth=1.5, label='tipo k=1')\n", "plt.legend(loc='best', shadow=True, framealpha=1)\n", "plt.grid(alpha=0.3)\n", "plt.xlabel('t')\n", "plt.show()\n", "\n" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.7.6" } }, "nbformat": 4, "nbformat_minor": 4 }