daviddoji/Kalman-and-Bayesian-Filters-in-Python
1
0 Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Moved chapter 4 as part of reorg

Browse Source
This commit is contained in:
Roger Labbe committed 2015-01-27 09:58:39 -08:00
1 parent b3b3b33ea6
commit bb531dde86
4 files changed
+23 -26
Download Patch File Download Diff File Expand all files Collapse all files

No files matched your search

animations/04_gaussian_animate.gif
BIN
View File
Binary file not shown.
Side by Side

After

Width:  |  Height:  |  Size: 624 KiB

animations/Gaussians_Animations.ipynb
+324
View File
@@ -0,0 +1,324 @@
{
"metadata": {
"name": "",
"signature": "sha256:aada3223e5395f37eadc453245d1072d753bcb8ccf236d890bad5a586a4ed406"
},
"nbformat": 3,
"nbformat_minor": 0,
"worksheets": [
{
"cells": [
{
"cell_type": "code",
"collapsed": false,
"input": [
"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='..')"
],
"language": "python",
"metadata": {},
"outputs": [
{
"html": [
"<style>\n",
"@import url('http://fonts.googleapis.com/css?family=Source+Code+Pro');\n",
"@import url('http://fonts.googleapis.com/css?family=Vollkorn');\n",
"@import url('http://fonts.googleapis.com/css?family=Arimo');\n",
"\n",
" div.cell{\n",
" width: 850px;\n",
" margin-left: 0% !important;\n",
" margin-right: auto;\n",
" }\n",
" div.text_cell code {\n",
" background: transparent;\n",
" color: #000000;\n",
" font-weight: 600;\n",
" font-size: 11pt;\n",
" font-style: bold;\n",
" font-family: 'Source Code Pro', Consolas, monocco, monospace;\n",
" }\n",
" h1 {\n",
" font-family: 'Open sans',verdana,arial,sans-serif;\n",
"\t}\n",
"\t\n",
" div.input_area {\n",
" background: #F6F6F9;\n",
" border: 1px solid #586e75;\n",
" }\n",
"\n",
" .text_cell_render h1 {\n",
" font-weight: 200;\n",
" font-size: 30pt;\n",
" line-height: 100%;\n",
" color:#c76c0c;\n",
" margin-bottom: 0.5em;\n",
" margin-top: 1em;\n",
" display: block;\n",
" white-space: wrap;\n",
" } \n",
" h2 {\n",
" font-family: 'Open sans',verdana,arial,sans-serif;\n",
" }\n",
" .text_cell_render h2 {\n",
" font-weight: 200;\n",
" font-size: 20pt;\n",
" font-style: italic;\n",
" line-height: 100%;\n",
" color:#c76c0c;\n",
" margin-bottom: 0.5em;\n",
" margin-top: 1.5em;\n",
" display: block;\n",
" white-space: nowrap;\n",
" } \n",
" h3 {\n",
" font-family: 'Open sans',verdana,arial,sans-serif;\n",
" }\n",
" .text_cell_render h3 {\n",
" font-weight: 300;\n",
" font-size: 18pt;\n",
" line-height: 100%;\n",
" color:#d77c0c;\n",
" margin-bottom: 0.5em;\n",
" margin-top: 2em;\n",
" display: block;\n",
" white-space: nowrap;\n",
" }\n",
" h4 {\n",
" font-family: 'Open sans',verdana,arial,sans-serif;\n",
" }\n",
" .text_cell_render h4 {\n",
" font-weight: 300;\n",
" font-size: 16pt;\n",
" color:#d77c0c;\n",
" margin-bottom: 0.5em;\n",
" margin-top: 0.5em;\n",
" display: block;\n",
" white-space: nowrap;\n",
" }\n",
" h5 {\n",
" font-family: 'Open sans',verdana,arial,sans-serif;\n",
" }\n",
" .text_cell_render h5 {\n",
" font-weight: 300;\n",
" font-style: normal;\n",
" color: #1d3b84;\n",
" font-size: 16pt;\n",
" margin-bottom: 0em;\n",
" margin-top: 1.5em;\n",
" display: block;\n",
" white-space: nowrap;\n",
" }\n",
" div.text_cell_render{\n",
" font-family: 'Arimo',verdana,arial,sans-serif;\n",
" line-height: 135%;\n",
" font-size: 125%;\n",
" width:750px;\n",
" margin-left:auto;\n",
" margin-right:auto;\n",
" text-align:justify;\n",
" text-justify:inter-word;\n",
" }\n",
" div.output_subarea.output_text.output_pyout {\n",
" overflow-x: auto;\n",
" overflow-y: scroll;\n",
" max-height: 50000px;\n",
" }\n",
" div.output_subarea.output_stream.output_stdout.output_text {\n",
" overflow-x: auto;\n",
" overflow-y: scroll;\n",
" max-height: 50000px;\n",
" }\n",
" code{\n",
" font-size: 70%;\n",
" }\n",
" .rendered_html code{\n",
" background-color: transparent;\n",
" }\n",
" ul{\n",
" margin: 2em;\n",
" }\n",
" ul li{\n",
" padding-left: 0.5em; \n",
" margin-bottom: 0.5em; \n",
" margin-top: 0.5em; \n",
" }\n",
" ul li li{\n",
" padding-left: 0.2em; \n",
" margin-bottom: 0.2em; \n",
" margin-top: 0.2em; \n",
" }\n",
" ol{\n",
" margin: 2em;\n",
" }\n",
" ol li{\n",
" padding-left: 0.5em; \n",
" margin-bottom: 0.5em; \n",
" margin-top: 0.5em; \n",
" }\n",
" ul li{\n",
" padding-left: 0.5em; \n",
" margin-bottom: 0.5em; \n",
" margin-top: 0.2em; \n",
" }\n",
" a:link{\n",
" font-weight: bold;\n",
" color:#447adb;\n",
" }\n",
" a:visited{\n",
" font-weight: bold;\n",
" color: #1d3b84;\n",
" }\n",
" a:hover{\n",
" font-weight: bold;\n",
" color: #1d3b84;\n",
" }\n",
" a:focus{\n",
" font-weight: bold;\n",
" color:#447adb;\n",
" }\n",
" a:active{\n",
" font-weight: bold;\n",
" color:#447adb;\n",
" }\n",
" .rendered_html :link {\n",
" text-decoration: underline; \n",
" }\n",
" .rendered_html :hover {\n",
" text-decoration: none; \n",
" }\n",
" .rendered_html :visited {\n",
" text-decoration: none;\n",
" }\n",
" .rendered_html :focus {\n",
" text-decoration: none;\n",
" }\n",
" .rendered_html :active {\n",
" text-decoration: none;\n",
" }\n",
" .warning{\n",
" color: rgb( 240, 20, 20 )\n",
" } \n",
" hr {\n",
" color: #f3f3f3;\n",
" background-color: #f3f3f3;\n",
" height: 1px;\n",
" }\n",
" blockquote{\n",
" display:block;\n",
" background: #fcfcfc;\n",
" border-left: 5px solid #c76c0c;\n",
" font-family: 'Open sans',verdana,arial,sans-serif;\n",
" width:680px;\n",
" padding: 10px 10px 10px 10px;\n",
" text-align:justify;\n",
" text-justify:inter-word;\n",
" }\n",
" blockquote p {\n",
" margin-bottom: 0;\n",
" line-height: 125%;\n",
" font-size: 100%;\n",
" }\n",
"</style>\n",
"<script>\n",
" MathJax.Hub.Config({\n",
" TeX: {\n",
" extensions: [\"AMSmath.js\"]\n",
" },\n",
" tex2jax: {\n",
" inlineMath: [ ['$','$'], [\"\\\\(\",\"\\\\)\"] ],\n",
" displayMath: [ ['$$','$$'], [\"\\\\[\",\"\\\\]\"] ]\n",
" },\n",
" displayAlign: 'center', // Change this to 'center' to center equations.\n",
" \"HTML-CSS\": {\n",
" styles: {'.MathJax_Display': {\"margin\": 4}}\n",
" }\n",
" });\n",
"</script>\n"
],
"metadata": {},
"output_type": "pyout",
"prompt_number": 1,
"text": [
"<IPython.core.display.HTML at 0x7fde2b0ceb38>"
]
}
],
"prompt_number": 1
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"import numpy as np\n",
"import matplotlib.pyplot as plt\n",
"from gif_animate import animate\n",
"from stats import gaussian\n",
"\n",
"def plt_g (mu, variance):\n",
" xs = np.arange(2,8,0.05)\n",
" ys = [gaussian (x, mu, variance) for x in xs]\n",
" plt.plot (xs, ys)\n",
" plt.ylim((0,1))\n",
"\n",
" \n",
"mu = 2\n",
"sigma = 0.6\n",
"def ganimate(frame):\n",
" global mu, sigma\n",
" if frame < 25:\n",
" mu += .2\n",
" elif frame == 25:\n",
" mu = 5\n",
" elif frame < 37:\n",
" sigma -= 0.05\n",
" else:\n",
" sigma += 0.05\n",
" \n",
" plt.cla()\n",
" plt_g(mu,sigma)\n",
"\n",
"animate('04_gaussian_animate.gif', ganimate, frames=80, interval=50)"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAcEAAAGrCAYAAAC41CzUAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl4VdW9//HPyXySnMwzhECYhwBCAA0oDpgrihYrINiq\nINaKI1LvVSyt7RWn0h8qiqiVgqLeoq3VqqigKBZBJIbIFMYwZ4DMc0LO2b8/jsabK5IAgZNkvV/P\nw/PIOnsn3yzi+Zy99lpr2yzLsgQAgIG8PF0AAACeQggCAIxFCAIAjEUIAgCMRQgCAIxFCAIAjEUI\nAgCM1WwIfvHFF7rmmmvUuXNneXl56ZVXXmn2i27ZskWjR49WYGCgOnfurEceeaRVigUAoDU1G4JV\nVVUaOHCgnnnmGdntdtlstpMeX15erssvv1zx8fHKyMjQM888o3nz5mn+/PmtVjQAAK3Bdio7xjgc\nDi1cuFA33XTTTx6zaNEizZ49WwUFBfL395ckPfroo1q0aJEOHz585hUDANBKWv2e4Pr163XhhRc2\nBqAkpaenKzc3VwcOHGjtbwcAwGnzae0vmJ+fry5dujRpi42NbXwtKSmpsb2srKy1vz0AAI1CQ0NP\n+nqrXwk2d88QAIC2otVDMC4uTvn5+U3aCgoKGl8DAKCtaPXh0AsuuEAPPPCA6urqGu8Lrlq1Sp06\ndWoyFPp/NXfJaoqMjAylpqZ6uow2gb5oiv5oiv5oiv74wancamvREomsrCxlZWXJ5XLpwIEDysrK\n0qFDhyRJs2fP1pgxYxqPv+GGGxQYGKipU6dq27Ztevvtt/Xkk09q1qxZp/GjAABw9jQbghs3btSQ\nIUM0ZMgQ1dbW6uGHH9aQIUP08MMPS3JPdsnJyWk8PiQkRKtWrVJubq5SU1N199136/7779d99913\n9n4KAABOQ7PDoRdffLFcLtdPvr5kyZIftQ0YMEBr1qw5s8oAADjL2DsUAGAsQhAAYCxCEABgLEIQ\nAGAsQhAAYCxCEABgLEIQAGAsQhAAYCxCEABgLEIQAGAsQhAAYCxCEABgLEIQAGAsQhAAYCxCEABg\nLEIQAGAsQhAAYCxCEABgLEIQAGAsQhAAYCxCEABgLEIQAGAsQhAAYCxCEABgLEIQAGAsQhAAYCxC\nEABgLEIQAGAsQhAAYCxCEABgLEIQAGAsQhAAYCxCEABgLEIQAGAsQhAAYCxCEABgLEIQAGAsQhAA\nYCxCEABgLEIQAGAsQhAAYCxCEABgLEIQAGAsQhAAYCxCEABgLEIQAGAsQhAAYCxCEABgLEIQAGAs\nQhAAYCxCEABgLEIQAGAsQhAAYCxCEABgLEIQAGAsQhAAYCxCEABgLEIQAGAsQhAAYCxCEABgLEIQ\nAGAsQhAAYCxCEABgLEIQAGAsQhAAYCxCEABgLEIQAGCsFoXg888/r27duslutys1NVVr16496fEr\nVqzQ+eefr5CQEEVHR2v8+PHavXt3qxQMAEBraTYEly9frpkzZ2rOnDnKyspSWlqaxo4dq0OHDp3w\n+D179mj8+PG6+OKLlZWVpU8++US1tbW68sorW714AADORLMhOH/+fE2bNk3Tp09X7969tWDBAsXH\nx2vRokUnPD4rK0sul0uPP/64kpOTNWjQID3wwAPau3eviouLW/0HAADgdJ00BOvr65WZman09PQm\n7enp6Vq3bt0Jzxk5cqSCg4P1l7/8RU6nUxUVFVq6dKmGDx+uiIiI1qscAIAzdNIQLCwslNPpVGxs\nbJP2mJgY5efnn/Cc+Ph4rVixQnPmzFFAQIDCwsK0bds2vffee61XNQAArcCntb9gTk6Oxo8fr2nT\npumGG25QeXm5fv/732vSpElavXq1bDbbCc/LyMho7VLaLfriB/RFU/RHU/RHU/SHW8+ePVt87ElD\nMCoqSt7e3iooKGjSXlBQoPj4+BOe8+KLLyoxMVFPPvlkY9trr72mxMRErV+/XmlpaSc8LzU1tcVF\nd2QZGRn0xXfoi6boj6boj6bojx+UlZW1+NiTDof6+flp6NChWrlyZZP2VatW/WSYWZYlL6+mX/b7\nv7tcrhYXBgDA2dbs7NBZs2Zp6dKlWrx4sbKzs3XvvfcqPz9ft99+uyRp9uzZGjNmTOPx11xzjTIz\nM/XII49o9+7dyszM1LRp09SlSxcNHTr07P0kAACcombvCU6aNElFRUWaO3eu8vLylJKSohUrVigx\nMVGSlJ+fr5ycnMbjR40apeXLl+uJJ57Qn/70JwUGBuqCCy7QRx99JLvdfvZ+EgAATlGLJsbMmDFD\nM2bMOOFrS5Ys+VHbhAkTNGHChDOrDACAs4y9QwEAxiIEAQDGIgQBAMYiBAEAxiIEAQDGIgQBAMYi\nBAEAxiIEAQDGIgQBAMYiBAEAxiIEAQDGIgQBAMYiBAEAxiIEAQDGIgQBAMYiBAEAxiIEAQDGIgQB\nAMYiBAEAxiIEAQDGIgQBAMYiBAEAxiIEAQDGIgQBAMYiBAEAxiIEAQDGIgQBAMYiBAEAxiIEAQDG\nIgQBAMYiBAEAxiIEAQDGIgQBAMYiBAEAxiIEAQDGIgQBAMYiBAEAxiIEAQDGIgQBAMYiBAEAxiIE\nAQDGIgQBAMYiBAEAxiIEAQDGIgQBAMYiBAEAxiIEAQDGIgQBAMYiBAEAxiIEAQDGIgQBAMYiBAEA\nxiIEAQDGIgQBAMYiBAEAxiIEAQDGIgQBAMYiBAEAxiIEAQDGIgQBAMYiBAEAxiIEAQDGIgQBAMYi\nBAEAxiIEAQDGIgQBAMYiBAEAxmpRCD7//PPq1q2b7Ha7UlNTtXbt2mbPefrpp9WnTx8FBAQoISFB\ns2fPPuNiAQBoTT7NHbB8+XLNnDlTixYt0qhRo7Rw4UKNHTtW27dvV2Ji4gnPmTVrlj744AP9+c9/\nVkpKisrKypSXl9fqxQMAcCaaDcH58+dr2rRpmj59uiRpwYIF+uijj7Ro0SI99thjPzp+586deu65\n57Rlyxb17t27sX3QoEGtWDYAAGfupMOh9fX1yszMVHp6epP29PR0rVu37oTnvPvuu0pOTtaKFSuU\nnJysbt26aerUqTp27FjrVQ0AQCs4aQgWFhbK6XQqNja2SXtMTIzy8/NPeE5OTo4OHDigN998U6++\n+qqWLVumHTt26Oqrr5ZlWa1XOQAAZ6jZ4dBT5XK5VFdXp2XLlqlHjx6SpGXLlql3797KyMjQsGHD\nTnheRkZGa5fSbtEXP6AvmqI/mqI/mqI/3Hr27NniY08aglFRUfL29lZBQUGT9oKCAsXHx5/wnPj4\nePn4+DQGoCT16NFD3t7eOnjw4E+GYGpqaouL7sgyMjLoi+/QF03RH03RH03RHz8oKytr8bEnHQ71\n8/PT0KFDtXLlyibtq1atUlpa2gnPGTVqlBoaGpSTk9PYlpOTI6fTqaSkpBYXBgDA2dbsOsFZs2Zp\n6dKlWrx4sbKzs3XvvfcqPz9ft99+uyRp9uzZGjNmTOPxY8aM0ZAhQ3TLLbcoKytLmzZt0i233KLz\nzz+fTykAgDal2XuCkyZNUlFRkebOnau8vDylpKRoxYoVjWsE8/Pzm1z12Ww2vf/++7rnnnt00UUX\nyW63Kz09XfPnzz97PwUAAKehRRNjZsyYoRkzZpzwtSVLlvyoLS4uTm+++eaZVQYAwFnG3qEAAGMR\nggAAYxGCAABjEYIAAGMRggAAYxGCAABjEYIAAGMRggAAYxGCAABjEYIAAGMRggAAYxGCAABjEYIA\nAGMRggAAYxGCAABjEYIAAGMRggAAYxGCAABjEYIAAGMRggAAYxGCAABjEYIAAGMRggAAYxGCAABj\nEYIAAGMRggAAYxGCAABjEYIAAGMRggAAYxGCAABjEYIAAGMRggAAYxGCAABjEYIAAGMRggAAYxGC\nAABjEYIAAGMRggAAYxGCAABjEYIAAGMRggAAYxGCAABjEYIAAGMRggAAYxGCAABjEYIAAGMRggAA\nYxGCAABjEYIAAGMRggAAYxGCAABjEYIAAGMRggAAYxGCAABjEYIAAGMRggAAYxGCAABjEYIAAGMR\nggAAYxGCAABjEYIAAGMRggAAYxGCAABjEYIAAGMRggAAYxGCAABjtSgEn3/+eXXr1k12u12pqala\nu3Zti7747t275XA45HA4zqhIAADOhmZDcPny5Zo5c6bmzJmjrKwspaWlaezYsTp06NBJz6uvr9fk\nyZM1evRo2Wy2VisYAIDW0mwIzp8/X9OmTdP06dPVu3dvLViwQPHx8Vq0aNFJz3vggQc0ePBgTZw4\nUZZltVrBAAC0lpOGYH19vTIzM5Went6kPT09XevWrfvJ8z744AN98MEHevbZZwlAAECb5XOyFwsL\nC+V0OhUbG9ukPSYmRvn5+Sc8Jzc3V7fddpveeecdBQYGtl6lAAC0spOG4Om48cYbNWPGDA0bNuyU\nzsvIyGjtUtot+uIH9EVT9EdT9EdT9Idbz549W3zsSUMwKipK3t7eKigoaNJeUFCg+Pj4E57z2Wef\n6YsvvtAf//hHSZJlWXK5XPL19dWiRYt06623nvC81NTUFhfdkWVkZNAX36EvmqI/mqI/mqI/flBW\nVtbiY08agn5+fho6dKhWrlyp6667rrF91apVmjhx4gnP2bp1a5O/v/POO3r00Ue1ceNGJSQktLgw\nAADOtmaHQ2fNmqUbb7xRw4cPV1paml544QXl5+fr9ttvlyTNnj1bGzdu1CeffCJJ6tevX5Pzv/76\na3l5ef2oHQAAT2s2BCdNmqSioiLNnTtXeXl5SklJ0YoVK5SYmChJys/PV05Ozkm/BusEAQBtUYsm\nxsyYMUMzZsw44WtLliw56blTp07V1KlTT7kwAADONvYOBQAYixAEABiLEAQAGIsQBAAYixLine truncated
"text": [
"<matplotlib.figure.Figure at 0x7fde0e233e80>"
]
}
],
"prompt_number": 2
},
{
"cell_type": "code",
"collapsed": false,
"input": [],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 2
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"<img src='gaussian_animate.gif'>"
]
}
],
"metadata": {}
}
]
}