16_09_29 deconv window.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "%matplotlib inline"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[<matplotlib.lines.Line2D at 0x7359358>]"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX0AAAEACAYAAABfxaZOAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XeYVEW+xvHvj6yogAkUJQgi6spiQoJeRlgkrXlXQQRF\nEVBQ1DWh64ruelfXsKASFTCDInLFNeKw4woKooAIDgxIEMHFCEZy3T+q0XEIEzjd1X36/TzPPE43\nZ3recuDXNVV1qsw5h4iIZIdyoQOIiEjqqOiLiGQRFX0RkSyioi8ikkVU9EVEsoiKvohIFim26JvZ\naDNbY2bzdnHNA2a22MzmmlnTaCOKiEhUStLTHwu039kfmllHoIFz7nCgDzAiomwiIhKxYou+c24a\n8M0uLjkTeDxx7UygmpnVjCaeiIhEKYox/drAykKPVyWeExGRNKOJXBGRLFIhgtdYBRxa6PEhiee2\nY2ba6EdEpAyccxbF65S0p2+Jjx2ZDPQAMLPmwFrn3JqdvZBzLrYft912W/AMmdS+NWsc117rqF7d\n0bWrY+pUx+bNJf/6n35yTJjg6NjRUauW4667HN99lz7tS6ePOLctG9oXpZIs2XwaeBtoZGafmFlP\nM+tjZr0TRfxlYJmZLQFGAldEmlBiZ9MmuOsuOPJI2LABFiyAp5+GU0+F8uVL/jpVqsAf/gAvvwyv\nvQazZ/vXHD8eIv53IhIbxQ7vOOcuKME1/aOJI3H3/vtwySVw0EEwaxYcdlg0r9ukCTzzDEybBlde\nCWPGwKOPwsEHR/P6InGhidwI5eTkhI6QVLvTPudgyBDo2BFuuAFeeSW6gl/YySf7N5NWreDYY2HS\npJJ/bZx/fnFuG8S/fVGyqMeLdvnNzFwqv5+khx9/hB49YMUK3xtPRrHfkRkz4Pzz4aKLYNAgKKcu\njmQoM8OleCJXpEzWrIGcHNhzTz/0kqqCD9C8Obz7Lkyd6sf+f/opdd9bJF2p6EvSLF4MLVr4IZ3H\nHoPKlVOfoWZNyM31bzodO8K336Y+g0g6UdGXpCgo8KtxBg6E228Hi+QX07KpXBkef9yv7GnbFr76\nKlwWkdBU9CVyBQXQpg3ccQdcdlnoNF65cjBsmH8jatcO1q4NnUgkDE3kSqQ++cSvoLn9dujZM3Sa\n7TkHAwb4Nf2vv+6HfUTSXZQTuSr6EplvvvEFv1cvuOaa0Gl2butW/4a0Zg28+CJUrBg6kciuqehL\n2lm/Hk47DU44Ae6/P3Sa4m3eDGedBbVrw4gRYeccRIqjoi9pxTno1g22bIFx4zJnPfy33/qbuC69\nFK6+OnQakZ2LsuhHscumZLn774dFi/w6/Ewp+AD77OOHd1q0gEaNoFOn0IlEkk89fdktU6f6Xv6M\nGVC3bug0ZfP2236o5513oEGD0GlEtqc7ciUtrFgBF1wATz2VuQUfoGVLuPVWOO88PzchEmfq6UuZ\nbN4MrVvDmWf6DdQynXO+6B9wgF/PL5JO1NOX4P72N6haFa67LnSSaJjBI4/AlCl+P36RuFJPX0pt\n2jT44x/9DU4HHRQ6TbTmzIH27f2+/4ceWvz1Iqmgnr4Es3YtXHghPPxw/Ao++D34r77a37y1dWvo\nNCLRU09fSqVHD9h7bxg6NHSS5Nm8GU45xU9SX3ll6DQiWqcvgbz0EkyfDvPmhU6SXBUq+F05W7Tw\nm7M1bhw6kUh01NOXElm3Dn7zG78vfps2odOkxrBhvvhPn166A9tFoqYxfUm566+Hzp2zp+AD9O3r\nN2MbMSJ0EpHoqKcvxcrN9ROb8+f7rQuySX6+H9//4AO/OZtICOrpS8ps2ACXX+6HOrKt4IM/batf\nP03oSnyo6Msu3XMPHHUU/P73oZOEM3AgfPQRTJoUOonI7tPwjuzU8uV+f/z33oN69UKnCevNN6F7\ndz/cU7Vq6DSSbTS8IykxYIA/ASvbCz74fYZatYK77w6dRGT3qKcvO/Svf8G118KHH0LlyqHTpIeV\nK6FpU/+bT/36odNINtHJWZJUGzb4cfxhw/w+NPKLv/3N788zcWLoJJJNNLwjSfXQQ77oq+Bv77rr\nfNHPzQ2dRKRs1NOXX/nyS79M8a23tP3AzkyaBH/+s1+7X0EbmUgKqKcvSXPHHf4wERX8nTvrLKhV\nC0aPDp1EpPTU05efFRT4owPz8/0JUrJz778Pp58OixdrCackn3r6khQ33OD32FHBL97xx/tlnPff\nHzqJSOmopy+Av/noootg4UKoUiV0msywdCk0a+bv1j3wwNBpJM60ZFMi5ZzfO75/f38qlpTcgAGw\nZYtf8SSSLCr6EqkXX4Sbb4a5c7VvfGl98YVf7fTOO3D44aHTSFylfEzfzDqY2UIzKzCzG3fw5/uY\n2WQzm2tmH5rZxVGEk+TbutUvP/zrX1Xwy+KAA/xWFX/+c+gkIiVTbE/fzMoBBUBbYDUwC+jinFtY\n6JqBwD7OuYFmtj+wCKjpnNtc5LXU008zzzwD990HM2eCRdKPyD4//AANGsCUKXDMMaHTSByluqff\nDFjsnFvhnNsEjAfOLHKNA/ZOfL438FXRgi/pZ/Nm+Mtf/NYCKvhlV7WqX/k0aFDoJCLFK0nRrw2s\nLPT408RzhT0EHGVmq4EPgAHRxJNkevJJf5NRu3ahk2S+vn39uP6cOaGTiOxaVDeRtwfmOOfamFkD\nYIqZNXHOfV/0wkGFukM5OTnk5OREFEFKY+NGuP12eOIJ9fKjsOeecNNNvrf/wguh00imy8vLIy8v\nLymvXZIx/ebAIOdch8TjmwDnnLu70DX/Av7unJueeJwL3Oice6/Ia2lMP00MH+5X7bz8cugk8bF+\nPTRsCP/3f/7wGZGopHpMfxbQ0MzqmlkloAswucg1K4DfJcLVBBoBS6MIKNHbuBHuugtuuy10knip\nUsUfraj/r5LOii36zrktQH/gdWABMN45l29mfcysd+KyvwEtzWweMAW4wTn3dbJCy+554gk44gg4\n6aTQSeKnVy9/8MyMGaGTiOyYbs7KMps3+x00x46FU04JnSaeRo70QzyvvBI6icSFNlyTMhs3Dg45\nRAU/mS6+2Pf2Z88OnURke+rpZ5EtW+Doo2HoUGjbNnSaePvnP+Htt2HChNBJJA7U05cymTgRatSA\nNm1CJ4m/3r39zqX5+aGTiPyain6W2LrV33l7661al58KVavCVVfB3XcXf61IKumEzywxeTJUqgQd\nO4ZOkj369fPr9pcvh3r1QqcR8dTTzwLOwd//Drfcol5+KtWo4Yd57rkndBKRX2giNwu8+aYvPvn5\nUE5v8ym1Zo3fb/+jj/w+RyJloYlcKZV77oHrrlPBD6FmTX8amc7SlXShnn7MzZ/vd9Fctkxn34by\nySfQtKk/U7d69dBpJBOppy8ldu+9cOWVKvgh1akDnTrBww+HTiKinn6sffopNGkCH3/sJxUlnDlz\n4PTTfW+/UqXQaSTTqKcvJTJkiN8SQAU/vGOP9XsePfNM6CSS7dTTj6m1a/25rXPm+OEFCe/VV/1B\nK3PmaOmslI56+lKskSP9OLIKfvpo397vf5SbGzqJZDP19GNowwaoX9/3LJs0CZ1GCnv0URg/3v9s\nREpKPX3Zpaef9sVeBT/9dO0K8+b5D5EQVPRjxjkYPBiuvTZ0EtmRypX9ElrdrCWhaMO1mMnLg02b\n/A1Zkp769PEbsa1aBbVrh04j2UY9/ZgZPBiuvlqrQ9LZvvtC9+7w4IOhk0g20kRujCxZAi1awIoV\nsOeeodPIrnz8MTRvrp+VlIwmcmWHHnwQevVSEckEDRr4N+gnnwydRLKNevoxsW6dX6Y5b54/+FzS\nX24uDBjgD1HXcJzsinr6sp2xY/3NPyr4mWPbWcVTp4bNIdlFRT8GtmyBBx7wE7iSOcz8ObpDhoRO\nItlERT8GXnzRH9Zx0kmhk0hpXXghvPOOn9gVSQUV/RgYPNiPDUvm2XNPuPRSeOih0EkkW2giN8Nt\n26d92TKoWDF0GimLTz7xWy8vXw577x06jaQjTeTKz4YMgf79VfAzWZ06flL30UdDJ5FsoJ5+Bvvi\nC2jUyN+Utd9+odPI7pg2DS65BBYu1AH2sj319AWAMWPgrLNU8OOgVSvYay9tuSzJp6KfobZs8Qel\nXHFF6CQSBTM/Ga/lm5JsKvoZ6rXX/MZdJ54YOolEpUsXmDsXCgpCJ5E4U9HPUMOHq5cfN5Ur+3H9\n4cNDJ5E400RuBlq+HI4/Hlau1OZqcbPtZ/vJJ1C1aug0ki40kZvlRo2CHj1U8OOoXj0/qTtuXOgk\nElclKvpm1sHMFppZgZnduJNrcsxsjpnNN7N/RxtTttmwAUaPhr59QyeRZOnXD4YO9UdfikSt2KJv\nZuWAh4D2wNFAVzNrXOSaasBQ4PfOud8Af0xCVgGefx6OOQaOOCJ0EkmWdu3g++9hxozQSSSOStLT\nbwYsds6tcM5tAsYDZxa55gJgonNuFYBz7stoY8o2w4bB5ZeHTiHJVK6c/xkPHRo6icRRSYp+bWBl\nocefJp4rrBGwr5n928xmmVn3qALKLz78EJYuhTPOCJ1Ekq1nT3jpJfj889BJJG6imsitABwHdAQ6\nALeaWcOIXlsShg+Hyy7TPjvZoEYNOOccP38jEqUKJbhmFVCn0ONDEs8V9inwpXNuPbDezP4D/BZY\nUvTFBg0a9PPnOTk55OTklC5xlvruOxg/3vf2JTv06wdnnw033ADly4dOI6mUl5dHXl5eUl672HX6\nZlYeWAS0BT4D3gW6OufyC13TGHgQ38uvDMwEznfOfVTktbROv4yGD4c33oCJE0MnkVRq0QIGDtSQ\nXrZL6Tp959wWoD/wOrAAGO+cyzezPmbWO3HNQuA1YB4wAxhVtOBL2Tnni74mcLPPFVdoQleipTty\nM8C0af50pfx8bbubbdavh7p14a23/Dbakp10R26WGT7c34ylgp99qlTx+/GMGBE6icSFevpp7vPP\nfQ9v2TK/okOyz/LlcMIJfj8ebb2RndTTzyJjxvileyr42atePWjZEp5+OnQSiQMV/TS2ZYv/tV5b\nKMsVV/i7sfWLsuwuFf009uqrcMAB/ld7yW6nnQbr1sG774ZOIplORT+NDRumXr542/bjGTYsdBLJ\ndJrITVPLlvmjEDV5J9t89RU0bAhLlsB++4VOI6mkidwsMHKkDkqRX9tvP39n7tixoZNIJlNPPw1t\n2AB16uiGHNnezJnQrZs/PF33bWQP9fRj7rnnoEkTFXzZXrNmUK0aTJkSOolkKhX9NKR9dmRnzDSh\nK7tHwztpZt486NTJ34VZoSQbX0vW+eEHP/w3Z47/r8SfhndibNtBKSr4sjNVq0L37jBqVOgkkonU\n008j337rb7mfPx8OPjh0GklnCxdCTo5f0lupUug0kmzq6cfUk09CmzYq+FK8xo3h6KPh+edDJ5FM\no6KfJrYdlKI7cKWktu3HI1IaKvppYto02LQJTj01dBLJFGecAR9/7IcDRUpKRT9NbFumaZGM2kk2\nqFjRT/oPHx46iWQSTeSmgTVr/BjtsmVQvXroNJJJVq2CY46BFStg771Dp5Fk0URuzIwZA+eeq4Iv\npVe7tp/8f+qp0EkkU6inH9iWLXDYYX4VxvHHh04jmSg3F665Bj74QMODcaWefoy88grUrKmCL2XX\npg1s3AjTp4dOIplART8wHZQiu2vbfjya0JWS0PBOQEuX+l0TV66EPfYInUYy2dq1UL8+LFoEBx4Y\nOo1ETcM7MTFyJFx0kQq+7L7q1f1igNGjQyeRdKeefiDr1/sdEqdPh8MPD51G4mD2bDj7bP8bZPny\nodNIlNTTj4HnnoOmTVXwJTrHHQe1avnFASI7o6IfiA5KkWTQfjxSHA3vBPDBB9C5sw5Kkej99JMf\nNpw509//IfGg4Z0MN3w49O6tgi/R22MPuPhiv0hAZEfU00+xb7+FunVhwQLtmy/JsWQJtGzpD1ip\nUiV0GomCevoZ7PHH4Xe/U8GX5GnY0E/qTpgQOomkIxX9FHLOT7L16xc6icSdJnRlZ1T0U+jNN/1/\nW7cOm0Pir3Nnv+3ynDmhk0i6UdFPoW377GgnREm28uWhTx/txyPb00Ruiqxe7Q+yXrEC9tkndBrJ\nBv/9Lxx5pF8aXK1a6DSyO1I+kWtmHcxsoZkVmNmNu7juRDPbZGbnRBEuTh5+GLp0UcGX1KlVC9q3\n94sHRLYptqdvZuWAAqAtsBqYBXRxzi3cwXVTgJ+AMc6553fwWlnZ09+0CerV87fHN2kSOo1kk//8\nxw/zfPSRhhUzWap7+s2Axc65Fc65TcB44MwdXHcl8BzweRTB4mTyZH93pAq+pNopp/jx/by80Ekk\nXZSk6NcGVhZ6/GniuZ+Z2cHAWc654YD6E0UMHaplmhKGDliRoqJavTMYKDzWr8KfkJ/vf7U+R7Mc\nEkj37jBlil9MIFKS3V9WAXUKPT4k8VxhJwDjzcyA/YGOZrbJOTe56IsNGjTo589zcnLIyckpZeTM\nMnw49OoFlSqFTiLZap99/CKCRx6Bv/wldBopiby8PPKSNCZXkonc8sAi/ETuZ8C7QFfnXP5Orh8L\nvKiJXPj+e7/j4dy5/r8iocybB506aWfXTJXSiVzn3BagP/A6sAAY75zLN7M+ZtZ7R18SRbA4eOop\nf/etCr6E1qSJX0H24ouhk0houjkrSZzzJ2Pdey+0axc6jQg8/TSMHevH9yWzaJfNDPD22/5Ai7Zt\nQycR8c491w/zFBSETiIhqegnybBhfqlcOf0fljRRuTJccgmMGBE6iYSk4Z0kWLMGGjeGpUuhRo3Q\naUR+sXw5nHCCP2Blzz1Dp5GS0vBOmhs1Cv7wBxV8ST/16kGLFjB+fOgkEop6+hHbuNH/w3rtNTjm\nmNBpRLb38stw220wa1boJFJS6umnseee80M7KviSrtq3h6++UtHPVir6EXvgAbjqqtApRHaufHno\n21fHKWYrDe9EaOZMf7v7kiX+H5ZIuvrySzj8cPj4Y9h339BppDga3klTDz4I/fur4Ev6239/+P3v\n4dFHQyeRVFNPPyKffQZHHaVlmpI53nkHevSARYt0P0m6U08/DY0Y4Yd2VPAlUzRv7nfgfOWV0Ekk\nldTTj8CGDVC3Lkyd6nv7IpniiSfgscfgjTdCJ5FdUU8/zTz7rN/FUAVfMs355/tDfubNC51EUkVF\nfzc5B0OGaJmmZKZKlfxRnoMHh04iqaLhnd309tt+MqygQJNhkpm2Ld9cuBBq1gydRnZEwztpZMgQ\n31NSwZdMtf/+cN55Ojw9W6invxuWL4fjj4dly/wqCJFMlZ8Pp57q/05XqRI6jRSlnn6aGDIEevZU\nwZfMd+SRcOyx/nQtiTf19Mto7Vo47DD44AM49NDQaUR23+uvw3XX+b/TFkmfUqKinn4aePhh6NhR\nBV/io1072LoVcnNDJ5FkUk+/DDZu9L38yZPhuONCpxGJziOPwKRJ8NJLoZNIYerpB/bss9CokQq+\nxE+3bvDee375psSTin4pOQf33gt/+lPoJCLR22MPuPxyuP/+0EkkWTS8U0q5uXDllTB/vtbmSzx9\n+aX/TXbBAjjooNBpBDS8E9R998G116rgS3ztv78f5hkyJHQSSQb19EthwQJo21Y3sEj8bbvxcOlS\nqFYtdBpRTz+Qe+7xWy6o4Evc1asHHTrAyJGhk0jU1NMvoRUr/GqdJUt0UIpkhw8+8PeiLFsGlSuH\nTpPd1NMP4N57oVcvFXzJHr/9rf944onQSSRK6umXwJo10Lix35SqVq3QaURSJy8P+vTxB62ULx86\nTfZSTz/FhgyBrl1V8CX7tG7tJ3JfeCF0EomKevrFWLfOb7nw3ntQv37oNCKpN3Ei3H03zJypjdhC\nUU8/hYYNg86dVfAle519Nvzwg9+FUzKfevq78OOPvpefmwtHHx06jUg448bBgw/C9Onq7Yegnn6K\njB4NLVqo4Iucdx58/TVMnRo6iewu9fR3Yv16f1j088/DiSeGTiMS3pNP+nMk3nwzdJLsk/Kevpl1\nMLOFZlZgZjfu4M8vMLMPEh/TzOyYKMKF9Mgj0LSpCr7INl26wOrVKvqZrtievpmVAwqAtsBqYBbQ\nxTm3sNA1zYF859w6M+sADHLONd/Ba2VET3/9emjY0C9TO/740GlE0sejj/qbtXS6VmqluqffDFjs\nnFvhnNsEjAfOLHyBc26Gc25d4uEMoHYU4UIZNcoXexV8kV/r1s1vyzBtWugkUlYlKfq1gZWFHn/K\nrot6L+CV3QkV0k8/wV13waBBoZOIpJ+KFeHmm+GOO0InkbKqEOWLmdmpQE/g5J1dM6hQNc3JySEn\nJyfKCLttxAg46SQ49tjQSUTSU48ecOed8NZbcMopodPEU15eHnl5eUl57ZKM6TfHj9F3SDy+CXDO\nubuLXNcEmAh0cM59vJPXSusx/R9/hAYN4NVX/UZTIrJjjz/uh0Hfekvr9lMh1WP6s4CGZlbXzCoB\nXYDJRQLVwRf87jsr+Jlg+HBo2VIFX6Q43brBN9/AKxk7kJu9SrROP7EiZwj+TWK0c+4uM+uD7/GP\nMrOHgXOAFYABm5xzzXbwOmnb01+3zq/LnzoVfvOb0GlE0t+kSXD77TB7to4PTbYoe/q6OSvhllvg\ns89gzJjQSUQyg3PQvDlcc41fwy/Jo6IfsdWr4ZhjYO5cOPTQ0GlEMkduLvTt6/fbr1gxdJr40t47\nEbv9drj0UhV8kdJq2xbq1vU3bUlmyPqe/sKFftnZokWw776h04hknnffhXPO8f+GqlYNnSae1NOP\n0C23wPXXq+CLlFWzZnDyyf4caUl/Wd3Tf/ttOP98KCiAPfYInUYkc61YAccdB/PmQe2M3oQlPamn\nH4GtW2HAAL/lggq+yO6pW9cfoH7zzaGTSHGytug/9phfbXDBBaGTiMTDwIEwZYo/T1rSV1YO73z7\nLTRu7LdO1n75ItF55BG/kkfbM0RLwzu76c47oX17FXyRqPXsCd9/DxMmhE4iO5N1Pf3Fi/25t/Pn\nQ61aQaOIxNKbb0L37v6Grb32Cp0mHnRHbhk5B507Q04O3HBDsBgisXfRRXDAAVrGGRUV/TJ69ln4\n61/9BlG6ZVwkeT7/3G9c+MYb0KRJ6DSZT0W/DNauhaOP9mONLVsGiSCSVUaN8qvk3npLu3DuLk3k\nlsHNN8Ppp6vgi6RKr17+fhjtXJtesqKn/847cO65sGAB1KiR8m8vkrXmzvUr5ebOhYMOCp0mc6mn\nXwrr1/sexz//qYIvkmpNm0Lv3n775TS4RUfIgqJ/661w1FFw3nmhk4hkp1tvhWXL4KmnQicRiPnw\nzrRpvtjPmwf775+ybysiRcyeDR06+GGegw8OnSbzaHinBL7/Hi6+GEaMUMEXCe244/wQT+/eGuYJ\nLbY9/b59/Xi+TvQRSQ8bN/rVcxdfDP37h06TWaLs6VeI4kXSzYQJ/qaQ998PnUREtqlUCcaP99ug\nnHyyn+SV1Ivd8M7HH0O/fvDMM1CtWug0IlJYw4YwZIg/vOj770OnyU6xGt7ZsAFatYIePeCqq5L2\nbURkN116KWze7IdftQVz8bQNw05cfjmsWQMTJ+ovkkg6++EHOOkkuOIK/yG7pjH9HRg2DP7zH3/3\nrQq+SHqrWtUfYtSqlb+PJicndKLsEYue/tSp/tjD6dOhQYPIX15EkiQ3F7p18521+vVDp0lfWqdf\nSEEBdO0K48ap4ItkmrZt/dm6p58O33wTOk12yOiiv2qV38zpf/8XTj01dBoRKYurroJ27Xzh//HH\n0GniL2OHd77+Gk45xZ/Qo1OwRDLb1q3+3/I338CkSTrkqKisX73zzTe+h9+6NdxzTwTBRCS4TZvg\nrLNgn33g8cdV+AvL6jH9L76ANm38HX3/+EfoNCISlYoV4bnnYN06f/PWxo2hE8VTRhX91av90q5O\nneC++7Q0UyRu9tjDD+8AnH22X88v0cqYoj97tr+Zo0cPuPNOFXyRuKpc2W+jcuCB8D//4xdsSHQy\nouhPmODH8AcPhhtvDJ1GRJKtYkV/tu755/vO3rvvhk4UH2k9kfvjj3DNNf4Gjmef9Xtyi0h2mTwZ\nLrsMrr0Wrr8eymVEVzVaKZ/INbMOZrbQzArMbId9bTN7wMwWm9lcM9vtTVNzc+HYY/2Y3uzZKvgi\n2eqMM2DWLHjpJb+IIz8/dKLMVmzRN7NywENAe+BooKuZNS5yTUeggXPucKAPMKKsgZYs8XfYXnqp\nX4755JN+CVcmyMvLCx0hqdS+zJXpbatTB/79bzj3XD/Of9118OWXv/x5prcvlUrS028GLHbOrXDO\nbQLGA2cWueZM4HEA59xMoJqZ1SxpCOf8O3mPHtC8ORxxBCxY4N/hM0nc/+KpfZkrDm0rXx6uvBLm\nz4fvvoNGjeBPf4KlS+PRvlQpyS6btYGVhR5/in8j2NU1qxLPrdnZi/70kz/Z6o03/ETt+vW+d//A\nA1C9egnTi0jWqVkTRo6EW2+F++/3HcXKlf2IQE4ONGkCFWKzf3D0Uj4l0rIlHHYY1KgBAwb44j9i\nhB/WuflmFXwRKZlDDvFFf9UqvyXL4sVw4YWw115+tKB1a+jcGT7/PHTS9FLs6h0zaw4Mcs51SDy+\nCXDOubsLXTMC+Ldz7pnE44VAa+fcmiKvlbqlQiIiMZLKQ1RmAQ3NrC7wGdAF6FrkmslAP+CZxJvE\n2qIFH6ILLSIiZVNs0XfObTGz/sDr+OGg0c65fDPr4//YjXLOvWxmncxsCfAD0DO5sUVEpCxSenOW\niIiElbKJ3JLc4JXOzOwQM5tqZgvM7EMzuyrxfA0ze93MFpnZa2ZWrdDXDEzcsJZvZqeFS19yZlbO\nzGab2eTE49i0z8yqmdmERN4FZnZSXNpnZteY2Xwzm2dmT5lZpUxvm5mNNrM1Zjav0HOlbpOZHZf4\n/1JgZoNT3Y4d2Unb/pHIPtfMJprZPoX+LLq2OeeS/oF/c1kC1AUqAnOBxqn43hG2oRbQNPH5XsAi\noDFwN3BD4vkbgbsSnx8FzMEPodVLtN9Ct6ME7bwGeBKYnHgcm/YBjwI9E59XAKrFoX3AwcBSoFLi\n8TPARZkCwyrFAAADDUlEQVTeNuBkoCkwr9BzpW4TMBM4MfH5y0D7NG3b74Byic/vAv6ejLalqqdf\nkhu80ppz7r/OubmJz78H8oFD8O14LHHZY8BZic/PAMY75zY755YDi9n+/oa0YmaHAJ2ARwo9HYv2\nJXpNpzjnxgIkcq8jJu0DygNVzawCsAf+XpmMbptzbhpQ9OTcUrXJzGoBezvnZiWue7zQ1wSzo7Y5\n595wzm1NPJyBry8QcdtSVfR3dINX7RR978iZWT38u/QMoKZLrFRyzv0XODBx2c5uWEtn/wSuBwpP\n9MSlffWBL81sbGL4apSZ7UkM2uecWw3cB3yCz7nOOfcGMWjbDhxYyjbVxtebbTKl9lyC77lDxG3L\nwv3qdo+Z7QU8BwxI9PiLzoRn5My4mXUG1iR+m9nV0tqMbB/+V+PjgKHOuePwq8xuIgY/PzOrju8B\n18UP9VQ1s27EoG0lELs2mdktwCbn3LhkvH6qiv4qoE6hx4cknssoiV+dnwOecM69kHh6zbZ9hhK/\nbm27/28VcGihL0/3NrcCzjCzpcA4oI2ZPQH8Nybt+xRY6Zx7L/F4Iv5NIA4/v98BS51zXzvntgCT\ngJbEo21FlbZNGdVWM7sYP8R6QaGnI21bqor+zzd4mVkl/A1ek1P0vaM0BvjIOTek0HOTgYsTn18E\nvFDo+S6JVRT1gYZA2h4F4Zy72TlXxzl3GP7nM9U51x14kXi0bw2w0swaJZ5qCywgHj+/T4DmZlbF\nzAzfto+IR9uMX//mWao2JYaA1plZs8T/mx6Fvia0X7XNzDrgh1fPcM5tKHRdtG1L4Wx1B/yKl8XA\nTaFnz8uQvxWwBb/yaA4wO9GmfYE3Em17Hahe6GsG4mfa84HTQrehFG1tzS+rd2LTPuC3+A7IXOB5\n/OqdWLQPuC2Rcx5+grNiprcNeBpYDWzAv7H1BGqUtk3A8cCHidozJHS7dtG2xcCKRG2ZDQxLRtt0\nc5aISBbRRK6ISBZR0RcRySIq+iIiWURFX0Qki6joi4hkERV9EZEsoqIvIpJFVPRFRLLI/wNGseKp\nV4N5IgAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x70527f0>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "w = np.hanning(1024)\n",
    "plt.plot(w)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "ename": "TypeError",
     "evalue": "'builtin_function_or_method' object has no attribute '__getitem__'",
     "output_type": "error",
     "traceback": [
      "\u001b[1;31m\u001b[0m",
      "\u001b[1;31mTypeError\u001b[0mTraceback (most recent call last)",
      "\u001b[1;32m<ipython-input-18-0dc6fb0ead23>\u001b[0m in \u001b[0;36m<module>\u001b[1;34m()\u001b[0m\n\u001b[0;32m      3\u001b[0m     \u001b[0mnp\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mreal\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mnp\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mfft\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mifft\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mnp\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mlog\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mnp\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mabs\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mnp\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mfft\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mfft\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mx\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m      4\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m----> 5\u001b[1;33m \u001b[1;32mprint\u001b[0m \u001b[0mnp\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mlog\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mnp\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0marray\u001b[0m\u001b[1;33m[\u001b[0m\u001b[1;36m0.0\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;36m0.0\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m      6\u001b[0m \u001b[0mWmp\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mrceps\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mw\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m      7\u001b[0m \u001b[1;32mprint\u001b[0m \u001b[0mWmp\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;31mTypeError\u001b[0m: 'builtin_function_or_method' object has no attribute '__getitem__'"
     ]
    }
   ],
   "source": [
    "#https://www.mathworks.com/help/signal/ref/rceps.html\n",
    "def rceps(x):\n",
    "    np.real(np.fft.ifft(np.log(np.abs(np.fft.fft(x)))))\n",
    "\n",
    "print np.log(np.array[0.0, 0.0])\n",
    "Wmp = rceps(w)\n",
    "print Wmp"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python [Root]",
   "language": "python",
   "name": "Python [Root]"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 2
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython2",
   "version": "2.7.12"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}