Moved chapter 5 as part of reorg.

This commit is contained in:
Roger Labbe committed 2015-01-27 10:18:25 -08:00
1 parent bb531dde86
commit 2a651f6c5b
6 files changed
+2426 -2438

No files matched your search

Binary file not shown.

After

Width:  |  Height:  |  Size: 179 KiB

Binary file not shown.

After

Width:  |  Height:  |  Size: 383 KiB

+270
View File
@@ -0,0 +1,270 @@
{
"metadata": {
"name": "",
"signature": "sha256:42fcc233db58d607109d988acbd2c05f37460792df171ef7155393554328c969"
},
"nbformat": 3,
"nbformat_minor": 0,
"worksheets": [
{
"cells": [
{
"cell_type": "code",
"collapsed": false,
"input": [
"\n",
"from __future__ import division, print_function\n",
"%matplotlib inline\n",
"import sys\n",
"sys.path.insert(0,'../code') # allow us to format the book\n",
"\n",
"# use same formatting as rest of book so that the plots are\n",
"# consistant with that look and feel.\n",
"import book_format\n",
"book_format.load_style(directory='..')\n",
"\n",
"\n",
"\n",
"def update(mu1, var1, mu2, var2):\n",
" mean = (var1*mu2 + var2*mu1) / (var1+var2)\n",
" variance = 1 / (1/var1 + 1/var2)\n",
" return (mean, variance)\n",
"\n",
"def predict(pos, variance, movement, movement_variance):\n",
" return (pos + movement, variance + movement_variance)"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 3
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"This notebook creates the animations for the Kalman Filter chapter. It is not really intended to be a readable part of the book, but of course you are free to look at the source code, and even modify it. However, if you are interested in running your own animations, I'll point you to the examples subdirectory of the book, which contains a number of python scripts that you can run and modify from an IDE or the command line. This module saves the animations to GIF files, which is quite slow and not very interactive. "
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"def plot_3d_covariance(ax, mean, cov):\n",
" \"\"\" plots a 2x2 covariance matrix positioned at mean. mean will be plotted\n",
" in x and y, and the probability in the z axis.\n",
"\n",
" Parameters\n",
" ----------\n",
" mean : 2x1 tuple-like object\n",
" mean for x and y coordinates. For example (2.3, 7.5)\n",
"\n",
" cov : 2x2 nd.array\n",
" the covariance matrix\n",
"\n",
" \"\"\"\n",
"\n",
" # compute width and height of covariance ellipse so we can choose\n",
" # appropriate ranges for x and y\n",
" o,w,h = stats.covariance_ellipse(cov,3)\n",
" # rotate width and height to x,y axis\n",
" wx = abs(w*np.cos(o) + h*np.sin(o))*1.2\n",
" wy = abs(h*np.cos(o) - w*np.sin(o))*1.2\n",
"\n",
"\n",
" # ensure axis are of the same size so everything is plotted with the same\n",
" # scale\n",
" if wx > wy:\n",
" w = wx\n",
" else:\n",
" w = wy\n",
"\n",
" minx = mean[0] - w\n",
" maxx = mean[0] + w\n",
" miny = mean[1] - w\n",
" maxy = mean[1] + w\n",
"\n",
" xs = np.arange(minx, maxx, (maxx-minx)/40.)\n",
" ys = np.arange(miny, maxy, (maxy-miny)/40.)\n",
" xv, yv = np.meshgrid (xs, ys)\n",
"\n",
" zs = np.array([100.* stats.multivariate_gaussian(np.array([x,y]),mean,cov) \\\n",
" for x,y in zip(np.ravel(xv), np.ravel(yv))])\n",
" zv = zs.reshape(xv.shape)\n",
"\n",
" ax = plt.figure().add_subplot(111, projection='3d')\n",
" ax.plot_surface(xv, yv, zv, rstride=1, cstride=1, cmap=cm.autumn)\n",
"\n",
" ax.set_xlabel('X')\n",
" ax.set_ylabel('Y')\n",
"\n",
" ax.contour(xv, yv, zv, zdir='x', offset=minx-1, cmap=cm.autumn)\n",
" ax.contour(xv, yv, zv, zdir='y', offset=maxy, cmap=cm.BuGn)\n"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 4
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"from __future__ import print_function, division\n",
"import matplotlib.pyplot as plt\n",
"import numpy.random as random\n",
"import math\n",
"import stats\n",
"from gif_animate import animate\n",
"\n",
"\n",
"\n",
"sensor_error = 1.2**2\n",
"movement_error = .2\n",
"movement = 0\n",
"voltage = (25,20) #who knows what the first value is?\n",
"\n",
"volts = [14.44, 18.74, 17.21, 15.98, 16.76, 14.8, 16.01, 16.90, 15.41, 17.31,\n",
" 16.39, 17.96]\n",
"ps = [voltage[0]]\n",
"\n",
"\n",
"i = 0\n",
"pred = 0\n",
"def volt_animate(frame):\n",
" global i, ps, voltage, Z, x, pred\n",
" \n",
" step = frame % 4\n",
" plt.subplot(211)\n",
" plt.ylim([14,26])\n",
"\n",
" if step == 0:\n",
" prev = voltage[0]\n",
" voltage = predict(voltage[0], voltage[1], movement, movement_error)\n",
" pred = voltage[0]\n",
" plt.plot ([i, i+1], [prev, pred], c='g') \n",
" \n",
" elif step == 1:\n",
" Z = volts[i]\n",
" i += 1\n",
" \n",
" plt.scatter(i, Z, marker='+', s=64, color='r')\n",
"\n",
" elif step == 2:\n",
" plt.plot ([i,i], [pred, Z], c='r', alpha=0.3)\n",
"\n",
" else:\n",
" voltage = update(voltage[0], voltage[1], Z, sensor_error)\n",
" ps.append(voltage[0])\n",
" plt.plot(ps, c='b')\n",
" \n",
" plt.subplot(212)\n",
" plt.cla()\n",
" stats.plot_gaussian(voltage[0], voltage[1], xlim=[0,32])\n",
" plt.ylim([0,1])\n",
" plt.tight_layout()\n",
" \n",
"N=12\n",
"animate('05_volt_animate.gif', volt_animate, N*3, 350)"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAcYAAAHGCAYAAADuYispAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3Xl8FPX9P/DX7iYh93LmDhAgJIBASQKFqIgCKUFA/FIo\nqFxSUQQbiP5EFBuqCHJWkEP0azECVqB+C4IpBiWAaWhNhCgGuQyGQLKRMxc5d+f3x7AXOTd7zB6v\n5+OxD2ZnZ3bes5vkxWfmM5+RCYIggIiIiAAAcqkLICIisicMRiIiIgMMRiIiIgMMRiIiIgMMRiIi\nIgMMRiIiIgMMRiIiIgMmB+PKlSsxePBgKJVKBAQEYMKECcjLy2uw3Pnz5/E///M/6NChA3x8fBAb\nG4uzZ89apGgiIiJrMTkYjx07hgULFuDEiRM4cuQI3NzcMGrUKNy6dUu3zKVLl3D//fejZ8+eyMjI\nQF5eHt566y34+vpatHgiIiJLk5k78k1lZSWUSiX279+PRx99FADwxBNPQKFQYMeOHRYpkoiIyFbM\nPsdYVlYGjUaDDh06AAA0Gg0OHjyIPn36YMyYMQgICMCQIUOwZ88es4slIiKyNrNbjFOmTMHPP/+M\nnJwcyGQyqFQqhISEwNvbG8uXL8cjjzyCr7/+Gi+//DL279+PsWPHWqp2IiIii3MzZ+Xk5GRkZWUh\nMzMTMpkMgNhiBICJEydi4cKFAIABAwYgJycHmzZtMgrG0tJSczZPRETUIqVSadLybT6UumjRIuze\nvRtHjhxB9+7ddfM7d+4MNzc39O3b12j56OhoXL58ua2bIyIisok2tRiTkpKwd+9eZGRkoHfv3kav\neXh4YPDgwQ0uzTh//rxRgBIREdkjk4Nx/vz52LlzJ/bt2welUgmVSgUA8PPzg4+PDwDg5ZdfxpQp\nU/Dggw/i4YcfRkZGBnbv3o39+/c3+b6mNnWdVU5ODuLi4qQuw27w8zDGz8MYPw9j/Dz0zDlVZ/Kh\n1K1bt6KiogIjR45ESEiI7rFu3TrdMo899hjef/99rF27FgMGDMDmzZuxY8cOJCYmtrlQIiIiWzC5\nxajtXNOSmTNnYubMmSYXREREJCWOlUpERGSAwUhERGSAwUhERGSAwUhERGSAwUhERGSAwUhERGSA\nwUhERGSAwUhERGSAwUhERGSAwUhERGTArPsxkuVkZgIJCUBV1SBgmUzqciCkmHX/aiIih8UWo524\ndAmoqgIAOXBmjKS1jDsH4MABSWsgIpIKW4x2Yvp0YMYMAJChS8a/8OtuCYthKBKRC2OL0Q5duyZ1\nBURErovBaEcUCqkrICIiBqMdMbzxtni+kYiIbI3BaEfS0gBA7A365JOSlkJE5LIYjHakY0f99MGD\n0tVBROTKGIx2R2wx1tVJXAYRkYtiMNoZb2+11CUQEbk0BqOd+eMfr+qmMzMlLISIyEWZHIwrV67E\n4MGDoVQqERAQgAkTJiAvL6/J5Z999lnI5XKsW7fOrEJdxfTp13XTU6dKWAgRkYsyORiPHTuGBQsW\n4MSJEzhy5Ajc3NwwatQo3Lp1q8Gy//jHP5CdnY2QkBDIZNKP/+lorl5teRkiIrIsk4eEO3TokNHz\nHTt2QKlUIisrC48++qhufkFBARYuXIivv/4aY8ZIO/ano5HJAIFjeBMRScLsc4xlZWXQaDTo0KGD\nbl59fT2mTZuG119/HVFRUeZuwuX06yd1BURErsvsYExKSsKgQYMwbNgw3byUlBQEBATg2WefNfft\nXZJho3zOHOnqICJyRWbdXSM5ORlZWVnIzMzUnUM8evQoUlNTkZuba7Ss0MKxwZycHHNKcSrFxTkA\nYgHIsGOHGvPmnbLp9pUXLgAASu3kO+HPhjF+Hsb4eRjj5yGKjIxs87oyoaXEasKiRYuwZ88eZGRk\noHfv3rr5f/nLX/DGG29ALtc3RtVqNeRyOUJCQnD58mXd/NLSUt20UqlsSxlOJycnB3FxcTDsq2Tz\n843a206NH2/jDTek/TxIxM/DGD8PY/w89MzJlza1GJOSkrB3794GoQgAzz//PCZPnqx7LggCfve7\n3+GJJ57AM88805bNuSQvLw4kTkQkBZODcf78+di5cyf27dsHpVIJlUoFAPDz84OPjw+6dOmCLl26\nGK3j7u6OoKAgs5q2ruaFF4DVq8Xp06eB/v2lrYeIyFWY3Plm69atqKiowMiRIxESEqJ78AJ+y1q1\nSj/92GPS1UFE5GpMbjFqNBqTN3Lp0iWT1yG9X36RugIiItfBsVLtmLYDDi/2JyKyHQajHeveXeoK\niIhcD4PRjqWn66cXL5auDiIiV8JgtGO9eumn331XujqIiFwJg9FB8JpGIiLbYDDaOQ8PqSsgInIt\nDEY799RT+mnen5GIyPoYjHbuww/10wkJ0tVBROQqGIwO5KefpK6AiMj5MRgdCC/0JyKyPgajAwgN\nlboCIiLXwWB0AJ9+qp9eu1a6OoiIXAGD0QE88IB++o03pKuDiMgVMBgdTHm51BUQETk3BqODcHeX\nugIiItfAYHQQ48bpp2/elK4OIiJnx2B0ELt26afHjpWuDiIiZ8dgdBBeXvrp7Gzp6iAicnYMRgek\n0UhdARGR82IwOpAuXaSugIjI+ZkcjCtXrsTgwYOhVCoREBCACRMmIC8vT/d6fX09Fi9ejIEDB8LX\n1xchISF48sknUVhYaNHCXdG2bfrpHTukq4OIyJmZHIzHjh3DggULcOLECRw5cgRubm4YNWoUbt26\nBQCorKzEqVOnsHTpUpw6dQr79+9HYWEhxowZA7VabfEdcCWPP66fTkqSrg4iImfmZuoKhw4dMnq+\nY8cOKJVKZGVl4dFHH4VSqUR6errRMtu2bUO/fv1w9uxZ9OvXz7yKCQBw9/8hRERkYWafYywrK4NG\no0GHDh2aXKa0tBQAml2GWkehkLoCIiLnZnYwJiUlYdCgQRg2bFijr9fW1uLFF1/EhAkTEBISYu7m\nXN6DD+qnq6qkq4OIyFmZFYzJycnIysrCZ599BplM1uD1+vp6PPXUUygrK8P27dvN2RTdlZamnx4/\nXro6iIiclUwQ2nb720WLFmHPnj3IyMhA7969G7xeX1+PadOmIS8vD0ePHkVAQECDZbSHWAHgwoUL\nbSnDJQ0eHAtABoVCg//856TF3195/DgAoHT4cIu/NxGRLURGRuqmlUqlSeu2KRiTkpKwd+9eZGRk\nICoqqsHrdXV1mDp1Ks6cOYOjR48iMDCw0fcxDEZTC3dWOTk5iIuLa3YZw8Z52/5b04IDB8R/7aBJ\n2prPw5Xw8zDGz8MYPw89c/LF5F6p8+fPx86dO7Fv3z4olUqoVCoAgJ+fH3x8fKBWqzF58mTk5OTg\nwIEDEARBt0z79u3h6elp6ibpHu3bA7dvS10FEZFzMvkc49atW1FRUYGRI0ciJCRE91i3bh0AoLCw\nEJ9//jmKi4sRGxtrtMyePXssvgOuaO1a/bThOUciIjKfyS1GTQsDdXbv3r3FZcg8c+YAf/yjOD17\nNlBSIm09RETOhGOlOrhff5W6AiIi58JgdFByfnNERFbBP68OKjZW6gqIiJwTg9FBGQ5ZO3mydHUQ\nETkbBqOD6thRP/3559LVQUTkbBiMTqC2VuoKiIicB4PRgfn6Sl0BEZHzYTA6sNde009nZkpXBxGR\nM2EwOrBXXtFPT5smXR1ERM6EwegkrlyRugIiIufAYHRwjdwGk4iIzMBgdHB9+khdARGRc2EwOrj0\ndP30s89KVwcRkbNgMDq40FD9dGqqdHUQETkLBqMTqamRugIiIsfHYHQCXl5SV0BE5DwYjE5g3jz9\n9MWL0tVBROQMGIxOYN06/XRionR1EBE5Awajk/n5Z6krICJybAxGJ6G90F8QpK2DiMjRmRyMK1eu\nxODBg6FUKhEQEIAJEyYgLy+vwXLLli1DaGgovL298fDDD+PMmTMWKZga162b1BUQETkHk4Px2LFj\nWLBgAU6cOIEjR47Azc0No0aNwq1bt3TLrFq1CuvXr8emTZuQnZ2NgIAAjB49GhUVFRYtnvQMb1b8\n5z9LVwcRkaMzORgPHTqEmTNnom/fvrjvvvuwY8cOXLt2DVlZWQAAQRDwzjvvYMmSJXj88cfRr18/\npKamory8HJ988onFd4BE/fvrpw074xARkWnMPsdYVlYGjUaDDh06AAAuXbqEkpISJCQk6Jbx9PTE\n8OHDdeFJ1nXnjtQVEBE5LrODMSkpCYMGDcKwYcMAACqVCgAQGBhotFxAQIDuNbIODw+pKyAicnxm\nBWNycjKysrLw2WefQdaK+x+1Zhlqu6lT9dNXr0pXBxGRI5MJQts6+C9atAh79uxBRkYGevfurZuf\nn5+PXr16ITs7G7Gxsbr5jz76KAICArB9+3bdvNLSUt30hQsX2lIG3WPw4FgAMvTsWYFPPz3bpvdQ\nHj8OACgdPtyClRER2U5kZKRuWqlUmrSuW1s2mJSUhL179zYIRQCIiIhAUFAQ0tPTdcFYXV2NzMxM\nrF27tsn3jIuLa0spTicnJ8cin8WlS75tf5/iYvFfO/hOLPV5OAt+Hsb4eRjj56Fn2PAylcnBOH/+\nfOzcuRP79u2DUqnUnTf08/ODj48PZDIZFi5ciBUrViA6OhqRkZFYvnw5/Pz88MQTT7S5UDKNRiN1\nBUREjsnkYNy6dStkMhlGjhxpNH/ZsmX4890L6F5++WVUVVVh/vz5uHXrFoYOHYr09HT4+PhYpmpq\nUnCwvsFHRESmMzkYNa1siqSkpCAlJcXkgsg8e/YADz4oTm/YACQlSVsPEZGj4VipTuaBB/TTr78u\nXR1ERI6KwejEysulroCIyPEwGJ2QW5v6GhMREcBgdEqGNyuuqpKuDiIiR8RgdEK7d+un7+Line truncated
"text": [
"<matplotlib.figure.Figure at 0x7fa8ae47c7f0>"
]
}
],
"prompt_number": 5
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"<img src='volt_animate.gif'>"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"from __future__ import print_function, division\n",
"import matplotlib.pyplot as plt\n",
"import numpy.random as random\n",
"import math\n",
"import stats\n",
"from gif_animate import animate\n",
"\n",
"# assume dog is always moving 1m to the right\n",
"movement = 1\n",
"movement_error = .05\n",
"sensor_error = 4.5\n",
"pos = (0, 100) # gaussian N(0,100)\n",
"\n",
"# this is the recorded output of a run of the dog sensor with an initial\n",
"# seed of 200\n",
"ZS = [-2.07, 6.05, 4.51, 3.47, 5.76, 5.93, 6.53, 9.01, 7.53, 11.68, \n",
" 11.15, 14.76, 13.45, 16.15, 19.05, 14.87, 20.90, 15.75, \n",
" 17.16, 20.50]\n",
"zs = []\n",
"ps = []\n",
"\n",
"N=20\n",
"\n",
"def dog_animate(frame):\n",
" global pos, zs, ps, N, ZS\n",
" pos = predict(pos[0], pos[1], movement, movement_error) \n",
" Z = ZS[frame]\n",
" zs.append(Z)\n",
" \n",
" pos = update(pos[0], pos[1], Z, sensor_error)\n",
" ps.append(pos[0])\n",
"\n",
"\n",
" plt.subplot(211)\n",
" plt.plot(zs,c='r', linestyle='dashed')\n",
"\n",
" plt.xlim([0,N*1.2])\n",
" plt.ylim([0,N*1.2])\n",
"\n",
" if len(ps) > 1:\n",
" plt.plot(ps, c='#8EBA42')\n",
" \n",
" plt.subplot(212)\n",
" plt.cla()\n",
" stats.plot_gaussian(pos[0], pos[1], xlim=[0,N*1.2])\n",
" plt.ylim(0, 1)\n",
" plt.tight_layout()\n",
"\n",
"\n",
"animate('05_dog_track.gif', dog_animate, N, 200)"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAcYAAAHFCAYAAABo9lmHAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3XlYVXX+B/D3uRcuuygim6BiIuaWCpqSuUPaYjku2a5Z\npqOTyzQZM85gaZrVOGqmWZmSbdjYr00zMU0jLLk6lCGuuMNFNkH25Z7fHwfuApf9Hu7C+/U8PJz7\nPduH0318d5bv9wiiKIogIiIiAIDC0gUQERFZEwYjERGRAQYjERGRAQYjERGRAQYjERGRAQdL7jw/\nP9+SuycionbO09OzThvPGImIiAwwGImIiAxY9FKqIVOns9RyarUa4eHhli7DLvHYyoPHVT48tsYa\nu43HM0YiIiIDDEYiIiIDDEYiIiIDDEYiIiIDDEYiIiIDDEYiIiIDDEYiIiIDDEYiIiIDDEYiIiID\nDEYiIiIDDEYiIiIDDEYiIiIDDEYiIiIDDEYiIiIDDEYiIiIDDEYiIiIDDEYiIiIDDEYiIiIDDEYi\nIiIDDEYiIiIDDEYisg+pqUBgIBAXZ+lKyMYxGInIPoSFAdevAzNnWroSsnEMRiKyP9u2WboCsmHN\nDsY1a9Zg6NCh8PT0hI+PDyZPnoyUlJQ6y61YsQJdu3aFq6srxo4di1OnTpmlYCIik4qLATc3afql\nlyxbSyuVlOfjl7SduHHrvKVLaZeaHYyHDx/GwoULcfToURw8eBAODg6YMGEC8vLydMusXbsW69at\nw6ZNm5CUlAQfHx9ERkaisLDQrMUTERn5xz+k39nZwIULlq2lBfJLMnDo9Ft4/6dHcfTCDqj3rwS0\nWkuX1e44NHeFffv2GX3euXMnPD09kZiYiPvuuw+iKGL9+vWIjo7GlClTAACxsbHw8fHBJ598grlz\n55qnciKi2qKjgVWrgOBgQGE7d4oyC85CfSkO5zKPQIQUhD2ScjDgq1PA0GtAt24WrrB9aXYw1lZQ\nUACtVotOnToBAC5evIjMzExERUXplnF2dsaoUaOQmJjIYCQieRUVWbqCJhFFEZdz1FBfjsPV3P8B\nABSCEn38IhF+PRDem14EXlzGULSAVgfjokWLMHjwYIwYMQIAoNFoAAC+vr5Gy/n4+CA9Pb21uyMi\nsmlV2kqc0RzC8cu7kF2YBgBQKV0xIPA+DO72J3g4+wD9AaQ+ALi7W7bYdqpVwbh06VIkJiYiISEB\ngiA0unxTliEiarKDB4EJE4CQEODMGUtX06DyymKcvL4HJy7vRmFZFgDATdUZg7tNwYDAB+DsWCsE\nPT0tUCUBrQjGJUuWYNeuXTh06BB69Oiha/fz8wMAZGZmIjAwUNeemZmpm2eKWq1uaSlUDx5T+fDY\nyqO5x7X/44/DWRQhnj2L41b636RcW4DrpT8jo/wXVImlAAAXhQ8CnUfBRzUYQo4D/sg5LXsd/M7q\nhYSENDi/RcG4aNEifP755zh06BB69+5tNC84OBh+fn7Yv38/wsLCAAClpaVISEjAm2++We82w8PD\nW1IK1UOtVvOYyoTHVh4tOq4ZGQAAoXdv0+uWlAAffQQ8+6wZKmyenMLLOH55F05n/IAqsQIA0LXj\nAIT3eBjB3ndCENru4SB+Z43l5+c3OL/ZwbhgwQJ89NFH+PLLL+Hp6am7p+jh4QE3NzcIgoDFixdj\n9erV6NOnD0JCQrBq1Sp4eHjg0UcfbdlfQURU29tv66c/+qju/O+/ByZOlKanTgW8vGQvKafwMn69\n+BEu56hRWlFQ3Sqgl8/dCO/xMPw9bze9YlISsG8fsGwZoFLJXic1rNnBuGXLFgiCgPHjxxu1r1ix\nAv/6178AAC+++CJKSkqwYMEC5OXlYfjw4di/fz/cajrfEhG1VvW/N1AqgaFD686fMAEQBEAUgUce\nkYJSBjVheCk7CWWVt3TtCkGJ/l3vxZBu09DJLbD+DZSXA3PmACdPAk5OwIsvylInNV2zg1HbxM6m\nMTExiImJaXZBRERNUjOoyF13mZ5fE5jHjgEHDph11/WFIQC4OHqiZ5cRGBb8GDq6BjS+sddek0Kx\nZ09g4UKz1kkt0+ruGkREFpGbCzz8MPDxx/Uv89lnUuBotcCaNdIAAC2UU3gJv178GJeyj6Gs0ngU\nLxfHjujZZQTuDH4Mnq7+Td9oSoo0IAEAvP8+4Ora4vrIfBiMRGSbOnZs/PJocDDQpQuQlSWdmTUz\nGGUJwxpVVdIl1IoKYO5cYOzY5m+DZMFgJCL79vrrUgCNHNmkxW8UnIf6clz9YegTgTt7PNqyMDRU\nUCCFe0CAVCNZDQYjEbXMt98C999v6SoaN2uW9FOPqqoqqC9/hj/Sv8OtEg1EiEbzzRqGhjp1Ar77\nTnqHJDvzWxUGIxE137/+BaxcCTg4SA/BuLsDX34JfPEF8OGHlq6uUTeLM/DTua24knMC5VV1x1Z1\ncfRET5+7cGePx+DpWv/AJK0mCEBgA0+skkUwGImoeTQaKRQBoLJS+j1lihSMAPDWW/KdAaWkAFFR\nwObNwIMPNmvV1IwfoL70GXKLLkMrVtWZr3JwR3evIRgZMhcdzXlmSDaHwUhEzdO1q35640bpbHHD\nBn0w9u0rXR6Uw8yZQHo68NBDUv/EBpSUF+KXCx/g7I0jKC7PqzNfgAAPFz8MCLgXYd0fhlKplKdm\nsjkMRiJqul699C/OHTgQ+MtfpOlu3QAfH+DGDSm4MjIAfxnOulJSpN9BQSZnp988hZ/Pb4MmPxWV\n2rI68x0UTvBXdEOEajwC7ppu/voaIopSl5FHHwUMxpcm68NgJKKmKSwELlyQph0cgN9+M56fkiJ1\njQCAAQOA7Gzz7v+LL/Rnie++azQr4+ZpfH58Caq05XVWc1V1Qi+fkYi47Rm4THkY2LsV8Pk/ILON\ng3H3buAf/wA2bQIuXpRGuSGrxGAkoqZxdwdiYqT7i3l1L03C21s6c7xyBcjJkYLR29t8+685OxUE\n3RiolZWV+DRpvu69hoA0FJuXW3eEdZ+BvgGRxtsYMQLYu1c6s714Uern2BZyc4EFC6TpmBiGopVr\nu+Hdicj2rVghdUyv7wW6qanSUGxTp5o3FAHpEi0ADBoEADh6YSfeOnSPUSjeETQFiybsxxMj3qsb\nigCwfDmgqP5nb+ZM89bXkCVLpDAePdoib/qg5uEZIxGZj6ur/klVc3v0UeCLL5Dz4QZ8evB+VFSV\n6GZ1cPHDU8O3wcHBufHtjBsnjZ2alCSFvNwP3ezbJ3VhcXYG3ntPH8xktfhfiIhsw8cfI+7IXHyY\nsUIXigpBiYcGrcGckR83LRQBIC5O+i2KwF//KlOxBi5dkl4l9corQCMvyCXrwDNGIjItN1e6B3fx\nYpu8y7AhyVe+xqEzGwGDUWlCfcfh3oH/aP7GvLyA7t2lJ2fNfbnXlHnzpLPUnj3l3xeZBYORiEzr\n0kXqmtG5c6N9BuVSWJqNnb88a/DSX8BV5YWnhm+Hs1M99zmb4sIF+S+hGurdu+32Ra3GS6lEVFef\nPvr+irfX89b5pnjtNSmAGno1VD2+/N8/8N5PD+tCUYAC9/SLxnOjP29dKAJtG4pkc3jGSETG3nwT\nOHNGmlYqgVOnWr6tmtc8zZoFPPZYk1Y5nXEI+/5YDRFSMCsqtOjmNxxThqxpeR1tqaBAuqfo3MR7\nnmR1eMZIRHq5ucDf/qb/fONG67Y3dar0u7ISePvtBhetrCzEu4dn4Ls/VulCMehkARY8dBhThv+7\ndXW0lfffl0YA+tOfLF0JtQKDkYj04uP102vWtP6hm//+Vz/9/PMmF8kv1iA2cTaO3lqJovKc6lYB\nY0IXYNqL/4NDhSj/a5mSk1u3/o0bQFiY1EexrAw4dAgoLTVPbdTmGIxEpPfww9LTmk89Bbz0knm2\n+cwz0m+tVhogAEBFRTE++/V5/Cd+Aj74+THkFl3RLe7v2Q9LIg9g8HV/qZ8hIN2rlEvnzsDgwcDa\ntS1b/403pIHVT5yQPg8fLg2izkupNksQRQs9bgYgPz9fN+3JF3WalVqtRnh4uKXLsEs8ti2gUEAr\nikh5MBg/LOgFEXVf+6SAEx4b8Ta83auHaQsJAc6fl4aAq3kQSA5dukjD13l6AjdvNm9drVbq8pGX\np+/A//jj8tTZCvzOGmsse/jwDRHJ6mDqJlz49G4UeyigdVQABqGoEBzQx2887un/ItRqtT4UAf2A\n5bfdJm+Br70mndXm5wPHj0uXRJtKoZAuF7/2mvTaLVdX+eqkNsNgJCKzO37pcyRe2IFKbfV9Ni/9\nPzUCFOjmNQQPDFwFR0dH0xuorJRCpqhIGk5NTnPmAPPnAxUVwBNPNP8p3HHjpB+yGwxGovZs4EDg\n5EnphcM1b69ooTOawziYuh6llQUm5grw8eiFKYPWwtW5CbdNHByk11xVVkrTcpsyBdi1SxoEvaQE\ncHExnq/VAosWAU8+CQwdKn89ZFEMRqL26rXXpFAEgKVLWxSM5eXl+OjYHOSXpJuc7+kSgMmDXoW3\ne7eW1dgWoQgAO3YAn38OeHhIl3D799fPS0oC7r9fevJ01y7p4SQOBG7XGIxE7dGFC/rO9wBw9Wqz\nN3Ho9GYkX91dp91V1Qn39H0JPbo08rDHhQvy3z9sKhcX6bKt4ZliZaX0dO4nn+jb7r237WujNsdg\nJGqPevXST0dHA35+TV71ZqEGsb88Ba2of72UUqHC2N7PY0DQpMY3sG6d/q0Wlnsovq7aoejrKw14\nAEjHZ88eYMgQy9RGbYrBSNTe+Pvrp/v1A1avbvKqcccWIT3/D6O2kb3mYmjww03ff2iowcojgYQE\n/ee2uqfYGAcH6dgkJACLF0vD5PHyabthBd9AImpTv/8uDVvm5AT88UfjywO4cOMovv5tuVLine truncated
"text": [
"<matplotlib.figure.Figure at 0x7fa8a6b732e8>"
]
}
],
"prompt_number": 6
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"<img src='05_dog_track.gif'>"
]
}
],
"metadata": {}
}
]
}