{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# EE16B: Homework 1" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "collapsed": false }, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib as plt\n", "from pylab import *\n", "import os\n", "import struct\n", "import wave\n", "\n", "# Plots graphs in the \n", "%matplotlib inline" ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "collapsed": false }, "outputs": [], "source": [ "def save_wav(buff, wavfile):\n", " # Write wav file\n", " print(\"Writing to \" + wavfile + \"...\")\n", " if os.path.isfile(wavfile):\n", " os.remove(wavfile)\n", " data = []\n", " for line in buff:\n", " d = max(min(int(float(line)*32767*0.99), 32767), -32767) \n", " data.append(struct.pack('h',d))\n", " w = wave.open(wavfile, 'wb')\n", " w.setparams((2, 2, 6600, len(data), 'NONE', 'not compressed'))\n", " for d in data:\n", " w.writeframes(d)\n", " w.writeframes(d)\n", " w.close()\n", " " ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "##Frequency Domain Filtering" ] }, { "cell_type": "code", "execution_count": 46, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "(0, 1.2)" ] }, "execution_count": 46, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXsAAAEACAYAAABS29YJAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAADthJREFUeJzt3W+sn2V9x/H3xxY2GTokXXBrazCzWTCpk5GVbrrxczo9\nayYsezDtlDj8x5NuS5ZoVx7ISZb9YU/mjAnrEIkzhj5Q50pGQBz+EkOgSEahSIttsElbEN3GmNMH\na9PvHpy7PT9+nJ7fOT336b/r/UpOet/3dfW6r/Ol+Zyr1/27S6oKSdKF7VVnewKSpOVn2EtSAwx7\nSWqAYS9JDTDsJakBhr0kNWBi2Cf5fJIXkuw5RfsHkjyR5MkkDyV5S//TlCQtxUJW9ncBU/O0Pwv8\nZlW9BfgL4B/7mJgkqT8Tw76qvgW8OE/7w1X1Une6C1jT09wkST3pe8/+I8C9PY8pSVqilX0NlOQd\nwIeBt/U1piSpH72EffdQ9g5gqqrm3PJJ4j/CI0mnoaqy1DGWvI2T5A3AV4EPVtWB+fpWlV9V3Hrr\nrWd9DufKl7WwFtZi/q++TFzZJ7kbuA5YleQQcCtwURfe24FPAa8Dbk8CcLSqNvQ2Q0nSkk0M+6ra\nPKH9o8BHe5uRJKl3vkF7FgwGg7M9hXOGtZhlLWZZi/6lzz2heW+U1Jm6lyRdKJJQ58IDWknSuc+w\nl6QGGPaS1ADDXpIaYNhLUgMMe0lqgGEvSQ0w7CWpAYa9JDXAsJekBhj2ktQAw16SGmDYS1IDDHtJ\naoBhL0kNMOwlqQGGvSQ1wLCXpAYY9pLUAMNekhpg2EtSAwx7SWqAYS9JDTDsJakBhr0kNcCwl6QG\nTAz7JJ9P8kKSPfP0+UyS/UmeSHJ1v1OUJC3VQlb2dwFTp2pMsgl4U1WtAz4O3N7T3CRJPZkY9lX1\nLeDFebpcD3yh67sLuCzJFf1MT5LUhz727FcDh0bODwNrehhXktSTlT2Nk7HzmqvTe9/b090kSYvS\nR9gfAdaOnK/prr3CpZdOnzxev37A+vWDHm4vSReOPXuG7Nkz7H3cVM25CH95p+RK4J6qWj9H2yZg\nS1VtSrIR+HRVbZyjXy3kXpKkWUmoqvHdk0WbuLJPcjdwHbAqySHgVuAigKraXlX3JtmU5ADwY+Cm\npU5KktSvBa3se7mRK3tJWrS+Vva+QStJDTDsJakBhr0kNcCwl6QGGPaS1ADDXpIaYNhLUgMMe0lq\ngGEvSQ0w7CWpAYa9JDXAsJekBhj2ktQAw16SGmDYS1IDDHtJaoBhL0kNMOwlqQGGvSQ1wLCXpAYY\n9pLUAMNekhpg2EtSAwx7SWqAYS9JDTDsJakBhr0kNcCwl6QGGPaS1ICJYZ9kKsm+JPuTbJ2jfVWS\n+5LsTvJUkj9alplKkk5bqurUjckK4BngXcAR4NvA5qraO9JnGvipqtqWZFXX/4qqOjY2Vs13L0nS\nKyWhqrLUcSat7DcAB6rqYFUdBXYAN4z1eR54bXf8WuA/x4NeknR2rZzQvho4NHJ+GLh2rM8dwINJ\nngNeA/xBf9OTJPVhUtgvZN/lFmB3VQ2S/CLwQJJfrqofjXecnp4+eTwYDBgMBouYqiRd+IbDIcPh\nsPdxJ+3ZbwSmq2qqO98GHK+q20b63Av8ZVU91J3/G7C1qh4bG8s9e0lapDO1Z/8YsC7JlUkuBt4H\n7Bzrs4+ZB7gkuQL4JeDZpU5MktSfebdxqupYki3A/cAK4M6q2pvk5q59O/BXwF1JnmDmh8cnq+q/\nlnnekqRFmHcbp9cbuY0jSYt2prZxJEkXAMNekhpg2EtSAwx7SWqAYS9JDTDsJakBhr0kNcCwl6QG\nGPaS1ADDXpIaYNhLUgMMe0lqgGEvSQ0w7CWpAYa9JDXAsJekBhj2ktQAw16SGmDYS1IDDHtJaoBh\nL0kNMOwlqQGGvSQ1wLCXpAYY9pLUAMNekhpg2EtSAwx7SWrAxLBPMpVkX5L9Sbaeos8gyeNJnkoy\n7H2WkqQlSVWdujFZATwDvAs4Anwb2FxVe0f6XAY8BLynqg4nWVVV/zHHWDXfvSRJr5SEqspSx5m0\nst8AHKiqg1V1FNgB3DDW5w+Br1TVYYC5gl6SdHZNCvvVwKGR88PdtVHrgMuTfDPJY0lu7HOCkqSl\nWzmhfSH7LhcBvwK8E7gEeDjJI1W1f7zj9PT0yePBYMBgMFjwRCWpBcPhkOFw2Pu4k/bsNwLTVTXV\nnW8DjlfVbSN9tgKvrqrp7vxzwH1V9eWxsdyzl6RFOlN79o8B65JcmeRi4H3AzrE+/wK8PcmKJJcA\n1wJPL3VikqT+zLuNU1XHkmwB7gdWAHdW1d4kN3ft26tqX5L7gCeB48AdVWXYS9I5ZN5tnF5v5DaO\nJC3amdrGkSRdAAx7SWqAYS9JDTDsJakBhr0kNcCwl6QGGPaS1ADDXpIaYNhLUgMMe0lqgGEvSQ0w\n7CWpAYa9JDXAsJekBhj2ktQAw16SGmDYS1IDDHtJaoBhL0kNMOwlqQGGvSQ1wLCXpAYY9pLUAMNe\nkhpg2EtSAwx7SWqAYS9JDTDsJakBE8M+yVSSfUn2J9k6T79fTXIsye/3O0VJ0lLNG/ZJVgCfBaaA\nNwObk1x1in63AfcBWYZ5SpKWYNLKfgNwoKoOVtVRYAdwwxz9/hj4MvDDnucnSerBpLBfDRwaOT/c\nXTspyWpmfgDc3l2q3mYnSerFygntCwnuTwN/XlWVJMyzjTM9PX3yeDAYMBgMFjC8JLVjOBwyHA57\nHzdVp87zJBuB6aqa6s63Acer6raRPs8yG/CrgJ8AH6uqnWNj1Xz3kiS9UhKqasnPQieF/UrgGeCd\nwHPAo8Dmqtp7iv53AfdU1VfnaDPsJWmR+gr7ebdxqupYki3A/cAK4M6q2pvk5q59+1InIElafvOu\n7Hu9kSt7SVq0vlb2vkErSQ0w7CWpAYa9JDXAsJekBhj2ktQAw16SGmDYS1IDDHtJaoBhL0kNMOwl\nqQGGvSQ1wLCXpAYY9pLUAMNekhpg2EtSAwx7SWqAYS9JDTDsJakBhr0kNcCwl6QGGPaS1ADDXpIa\nYNhLUgMMe0lqgGEvSQ0w7CWpAYa9JDXAsJekBiwo7JNMJdmXZH+SrXO0fyDJE0meTPJQkrf0P1VJ\n0ulKVc3fIVkBPAO8CzgCfBvYXFV7R/r8GvB0Vb2UZAqYrqqNY+PUpHtJkl4uCVWVpY6zkJX9BuBA\nVR2sqqPADuCG0Q5V9XBVvdSd7gLWLHVikqT+LCTsVwOHRs4Pd9dO5SPAvUuZlCSpXysX0GfBey9J\n3gF8GHjbac9IktS7hYT9EWDtyPlaZlb3L9M9lL0DmKqqF+caaHp6+uTxYDBgMBgsYqqSdOEbDocM\nh8Pex13IA9qVzDygfSfwHPAor3xA+wbgQeCDVfXIKcbxAa0kLVJfD2gnruyr6liSLcD9wArgzqra\nm+Tmrn078CngdcDtSQCOVtWGpU5OktSPiSv73m7kyl6SFu1MfvRSknSeM+wlqQGGvSQ1wLCXpAYY\n9pLUAMNekhpg2EtSAwx7SWqAYS9JDTDsJakBhr0kNcCwl6QGGPaS1ADDXpIaYNhLUgMMe0lqgGEv\nSQ0w7CWpAYa9JDXAsJekBhj2ktQAw16SGmDYS1IDDHtJaoBhL0kNMOwlqQGGvSQ1wLCXpAZMDPsk\nU0n2JdmfZOsp+nyma38iydX9T1OStBTzhn2SFcBngSngzcDmJFeN9dkEvKmq1gEfB25fprleMIbD\n4dmewjnDWsyyFrOsRf8mrew3AAeq6mBVHQV2ADeM9bke+AJAVe0CLktyRe8zvYD4B3mWtZhlLWZZ\ni/5NCvvVwKGR88PdtUl91ix9apKkvkwK+1rgODnN3ydJOgNSdepcTrIRmK6qqe58G3C8qm4b6fMP\nwLCqdnTn+4DrquqFsbH8ASBJp6GqxhfUi7ZyQvtjwLokVwLPAe8DNo/12QlsAXZ0Pxz+ezzo+5qs\nJOn0zBv2VXUsyRbgfmAFcGdV7U1yc9e+varuTbIpyQHgx8BNyz5rSdKizLuNI0m6MCz7G7QLeSnr\nfJfk80leSLJn5NrlSR5I8t0kX09y2Ujbtq4e+5K8e+T6NUn2dG1/f6a/jz4kWZvkm0m+k+SpJH/S\nXW+uHkl+OsmuJLuTPJ3kr7vrzdXihCQrkjye5J7uvMlaJDmY5MmuFo9215a3FlW1bF/MbP0cAK4E\nLgJ2A1ct5z3PxhfwG8DVwJ6Ra38LfLI73gr8TXf85q4OF3V1OcDs37AeBTZ0x/cCU2f7ezuNWrwe\neGt3fCnwDHBVw/W4pPt1JfAI8PZWa9HN/c+ALwE7u/MmawF8D7h87Nqy1mK5V/YLeSnrvFdV3wJe\nHLt88mWz7tff645vAO6uqqNVdZCZ/3DXJvl54DVV9WjX759Gfs95o6q+X1W7u+P/BfYy8y5Gq/X4\nSXd4MTOLnxdptBZJ1gCbgM8x+3HtJmvRGf/QyrLWYrnDfiEvZV2orqjZTyW9AJx4q/gXmKnDCSdq\nMn79COd5rbpPcV0N7KLReiR5VZLdzHzP36yq79BoLYC/Az4BHB+51motCvhGkseSfKy7tqy1mPTR\ny6Xy6S9QVdXaewZJLgW+AvxpVf0omV3EtFSPqjoOvDXJzwL3J3nHWHsTtUjyu8APqurxJIO5+rRS\ni87bqur5JD8HPNC9n3TSctRiuVf2R4C1I+dreflPogvZC0leD9D9desH3fXxmqxhpiZHePk/M7Gm\nu3beSXIRM0H/xar6Wne52XoAVNVLwL8C19BmLX4duD7J94C7gd9K8kXarAVV9Xz36w+Bf2Zmy3tZ\na7HcYX/ypawkFzPzUtbOZb7nuWIn8KHu+EPA10auvz/JxUneCKwDHq2q7wP/k+TazCyDbxz5PeeN\nbu53Ak9X1adHmpqrR5JVJz5RkeTVwG8Dj9NgLarqlqpaW1VvBN4PPFhVN9JgLZJckuQ13fHPAO8G\n9rDctTgDT51/h5lPZBwAtp3tp+DL9D3ezcwbxv/HzDOKm4DLgW8A3wW+Dlw20v+Wrh77gPeMXL+m\n+49+APjM2f6+TrMWb2dmT3Y3M8H2ODP/RHZz9QDWA//e1eJJ4BPd9eZqMVaX65j9NE5ztQDe2P2Z\n2A08dSIXl7sWvlQlSQ3wf0soSQ0w7CWpAYa9JDXAsJekBhj2ktQAw16SGmDYS1IDDHtJasD/A6a3\nflrsFjH3AAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "#Create a delta wave in time, by filling frequency domain with ones\n", "X = np.ones(5000, dtype=float)\n", "\n", "plot(X)\n", "ylim([0,1.2])\n" ] }, { "cell_type": "code", "execution_count": 47, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Writing to delta.wav...\n" ] }, { "name": "stderr", "output_type": "stream", "text": [ "/home/emily/anaconda3/lib/python3.4/site-packages/numpy/core/numeric.py:462: ComplexWarning: Casting complex values to real discards the imaginary part\n", " return array(a, dtype, copy=False, order=order)\n", "/home/emily/anaconda3/lib/python3.4/site-packages/IPython/kernel/__main__.py:8: ComplexWarning: Casting complex values to real discards the imaginary part\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYMAAAEACAYAAABRQBpkAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAD3BJREFUeJzt3H+s3Xddx/Hni9stbjKcy8xwbWEYqnbGOSR2Cww5iMKl\n6mr8g1EVdQrsD6tGEyjDxN0FFGdwAtmyVBxkErPFwNSi68aQnUBkbFbX/YAWeh1N2g4GKD+UH7HN\n3v5xvi2nh3u7237POffce56P5OZ+f3zu9/M57zb3db/nfD+fVBWSpOn2jOUegCRp+RkGkiTDQJJk\nGEiSMAwkSRgGkiRahkGS9yZ5MsmjJ2nz7iT7kzyc5AVt+pMkjUbbO4P3AbOLnUyyGXh+VW0A3gDc\n0rI/SdIItAqDqvo48JWTNLkSuK1p+wBwbpIL2vQpSRq+UX9msBY42Ld/CFg34j4lSadoHB8gZ2Df\n9S8kacKsGfH1DwPr+/bXNcdOkMSAkKTTUFWDf3CfllHfGewEfh0gyeXAV6vqyYUaVtV3fb31rcVL\nXnLdguem7eu666yDNbAO1uHEr2FqdWeQ5HbgpcD5SQ4C1wFnAFTVjqq6K8nmJPPAN4Cr2w5YkjR8\nrcKgqrYuoc22Nn1IkkZv4mcgP/e5neUewkTodDrLPYRlZw16rEOPdRiuDPt9p9MaRFILjeNtb4Nv\nf7v3XZJ0oiTUCvkAWZK0AhgGkiTDQJJkGEiSWAFhMAGfb0vSqjfRYZChfEYuSXo6Ex0GkqTxMAwk\nSYaBJMkwkCRhGEiSMAwkSRgGkiQMA0kShoEkiRUQBi5HIUmjN9Fh4HIUkjQeEx0GkqTxMAwkSYaB\nJMkwkCRhGEiSMAwkSRgGkiQMA0kShoEkiSGEQZLZJPuS7E+yfYHz5ye5O8meJI8l+c1Tub7LUUjS\n6LUKgyQzwE3ALHAxsDXJxoFm24CHqupSoAP8RZI1S7t+m9FJkpaq7Z3BJmC+qg5U1RHgDmDLQJvP\nA89qtp8F/FdVHW3ZryRpiJb0F/pJrAUO9u0fAi4baPMe4KNJngDOAV7dsk9J0pC1vTNYyjv6bwH2\nVNWFwKXAzUnOadmvJGmI2t4ZHAbW9+2vp3d30O9FwJ8AVNV/Jvkc8CPA7v5Gc3Nzx7c7nQ6dTqfl\n0CRpdel2u3S73ZFcO9XicZ3mg+DPAC8HngAeBLZW1d6+NjcCX6uq65NcAPw7cElV/Xdfm1poHG9/\nO3z9673vkqQTJaGqhvKoTas7g6o6mmQbcA8wA9xaVXuTXNOc3wH8KfC+JA/Te1vqTf1BIElafm3f\nJqKqdgG7Bo7t6Nv+MvCLbfuRJI2OM5AlSYaBJGkFhIHLUUjS6E10GLgchSSNx0SHgSRpPAwDSZJh\nIEkyDCRJGAaSJAwDSRKGgSQJw0CShGEgSWIFhIHLUUjS6E10GLgchSSNx0SHgSRpPAwDSZJhIEky\nDCRJGAaSJAwDSRKGgSQJw0CShGEgSWIFhIHLUUjS6E10GLgchSSNx0SHgSRpPAwDSZJhIEkaQhgk\nmU2yL8n+JNsXadNJ8lCSx5J02/YpSRquNW1+OMkMcBPws8Bh4N+S7KyqvX1tzgVuBl5ZVYeSnN+m\nT0nS8LW9M9gEzFfVgao6AtwBbBlo8yvAB6vqEEBVfblln5KkIWsbBmuBg337h5pj/TYA5yW5L8nu\nJK9t2ackachavU0ELGVK2BnATwIvB84G7k/yyara399obm7u+Han06HT6bQcmiStLt1ul263O5Jr\ntw2Dw8D6vv319O4O+h0EvlxV3wK+leRjwE8Ai4aBJOm7Df6hfP311w/t2m3fJtoNbEhyUZIzgauA\nnQNt/hG4IslMkrOBy4BPL7UDl6OQpNFrdWdQVUeTbAPuAWaAW6tqb5JrmvM7qmpfkruBR4CngPdU\n1ZLCwOUoJGk82r5NRFXtAnYNHNsxsP8O4B1t+5IkjYYzkCVJhoEkyTCQJGEYSJIwDCRJGAaSJAwD\nSRKGgSSJFRAGLkchSaM30WHgchSSNB4THQaSpPEwDCRJhoEkyTCQJGEYSJIwDCRJGAaSJAwDSRKG\ngSSJFRAGLkchSaM30WHgchSSNB4THQaSpPEwDCRJhoEkyTCQJGEYSJIwDCRJGAaSJIYQBklmk+xL\nsj/J9pO0+6kkR5P8cts+JUnD1SoMkswANwGzwMXA1iQbF2l3A3A34FQySZowbe8MNgHzVXWgqo4A\ndwBbFmj3u8AHgC+dagcuRyFJo9c2DNYCB/v2DzXHjkuyll5A3NIcWvKvd5ejkKTxWNPy55fyi/2d\nwJurqpKERd4mmpubO77d6XTodDothyZJq0u326Xb7Y7k2qkW78MkuRyYq6rZZv9a4KmquqGvzeN8\nJwDOB74JvL6qdva1qYXGceONcOhQ77sk6URJqKqhvIfS9s5gN7AhyUXAE8BVwNb+BlX1Q8e2k7wP\n+FB/EEiSll+rMKiqo0m2AfcAM8CtVbU3yTXN+R1DGKMkacTa3hlQVbuAXQPHFgyBqrq6bX+SpOFz\nBrIkyTCQJBkGkiQMA0kSKyAMXI5CkkZvosPA5SgkaTwmOgwkSeNhGEiSDANJkmEgScIwkCRhGEiS\nMAwkSRgGkiQMA0kSKyAMXI5CkkZvosPA5SgkaTwmOgwkSeNhGEiSDANJkmEgScIwkCRhGEiSMAwk\nSRgGkiRWQBg4A1mSRm+iw8AZyJI0HhMdBpKk8WgdBklmk+xLsj/J9gXO/2qSh5M8kuRfk1zStk9J\n0nC1CoMkM8BNwCxwMbA1ycaBZo8DP11VlwBvBf6qTZ+SpOFre2ewCZivqgNVdQS4A9jS36Cq7q+q\nrzW7DwDrWvYpSRqytmGwFjjYt3+oObaY3wbuatmnJGnI1rT8+SU/+JnkZcBvAS9e6Pzc3Nzx7U6n\nQ6fTaTk0SVpdut0u3W53JNdOtXiQP8nlwFxVzTb71wJPVdUNA+0uAe4EZqtqfoHr1ELjeNe74PHH\ne98lSSdKQlUN5SH8tm8T7QY2JLkoyZnAVcDO/gZJnkMvCH5toSCQJC2/Vm8TVdXRJNuAe4AZ4Naq\n2pvkmub8DuCPge8HbklvFtmRqtrUbtiSpGFq+5kBVbUL2DVwbEff9uuA153+9U9/bJKkpZnoGcgu\nRyFJ4zHRYSBJGg/DQJJkGEiSDANJEoaBJAnDQJKEYSBJwjCQJGEYSJJYAWHgchSSNHoTHQYuRyFJ\n4zHRYSBJGg/DQJJkGEiSDANJEoaBJAnDQJKEYSBJwjCQJGEYSJJYAWHgchSSNHoTHQYuRyFJ4zHR\nYSBJGg/DQJJkGEiSDANJEoaBJIkhhEGS2ST7kuxPsn2RNu9uzj+c5AVt+5QkDVerMEgyA9wEzAIX\nA1uTbBxosxl4flVtAN4A3NKmT0nS8LW9M9gEzFfVgao6AtwBbBlocyVwG0BVPQCcm+SClv1Kkoao\nbRisBQ727R9qjj1dm3Ut+5UkDdGalj+/1MUiBucSL+nnzjoLbr4Z7rzz1AYlSTo1bcPgMLC+b389\nvb/8T9ZmXXPsBHNzc8e3O50OnU6Hq6+GV72q5QglaZX4xCe63H9/9/j+jTcO79qpFivBJVkDfAZ4\nOfAE8CCwtar29rXZDGyrqs1JLgfeWVWXD1yn2oxDkqZREqpqKKu4tbozqKqjSbYB9wAzwK1VtTfJ\nNc35HVV1V5LNSeaBbwBXtx61JGmoWt0ZDG0Q3hlI0ikb5p2BM5AlSYaBJMkwkCRhGEiSMAwkSRgG\nkiQMA0kShoEkCcNAkoRhIEnCMJAkYRhIkjAMJEkYBpIkDANJEoaBJAnDQJKEYSBJwjCQJGEYSJIw\nDCRJGAaSJAwDSRKGgSQJw0CShGEgScIwkCRhGEiSaBEGSc5Lcm+Szyb5cJJzF2izPsl9ST6V5LEk\nv9duuJKkUWhzZ/Bm4N6q+mHgX5r9QUeAP6iqHwMuB34nycZT6aTb7bYY4uphHazBMdahxzoMV5sw\nuBK4rdm+DfilwQZV9YWq2tNs/y+wF7jwVDrxH7zHOliDY6xDj3UYrjZhcEFVPdlsPwlccLLGSS4C\nXgA80KJPSdIIrDnZyST3As9e4NQf9e9UVSWpk1znmcAHgN9v7hAkSRMkVYv+Dj/5Dyb7gE5VfSHJ\nDwL3VdWPLtDuDOCfgF1V9c5FrnV6g5CkKVdVGcZ1Tnpn8DR2Ar8B3NB8/4fBBkkC3Ap8erEggOG9\nGEnS6WlzZ3Ae8HfAc4ADwKur6qtJLgTeU1U/n+QK4GPAI8Cxjq6tqrtbj1ySNDSnHQaSpNVjomcg\nJ5lNsi/J/iTbl3s8w5TkvUmeTPJo37FFJ/Ilubapw74kr+g7/sIkjzbn3jXu19HGYpMSp7AO35Pk\ngSR7knw6ydub41NVB4AkM0keSvKhZn8aa3AgySNNHR5sjo2+DlU1kV/ADDAPXAScAewBNi73uIb4\n+l5C71HbR/uO/TnwpmZ7O/BnzfbFzes/o6nHPN+5q3sQ2NRs3wXMLvdrO4UaPBu4tNl+JvAZYOO0\n1aEZ89nN9zXAJ4ErprQOfwj8LbCz2Z/GGnwOOG/g2MjrMMl3BpuA+ao6UFVHgDuALcs8pqGpqo8D\nXxk4vNhEvi3A7VV1pKoO0PsHv6x5iuucqnqwafc3LDD5b1LVwpMS1zJldQCoqm82m2fS+0PoK0xZ\nHZKsAzYDfw0ce6hkqmrQZ/ChmpHXYZLDYC1wsG//UHNsNVtsIt+F9F7/McdqMXj8MCu0RgOTEqeu\nDkmekWQPvdd7X1V9iumrw18CbwSe6js2bTWA3sM2H0myO8nrm2Mjr0ObR0tHbao/2a46+US+1aSZ\nlPhBepMS/6f3RHLPtNShqp4CLk3yfcA9SV42cH5V1yHJLwBfrKqHknQWarPaa9DnxVX1+SQ/ANzb\nzOk6blR1mOQ7g8PA+r799ZyYdKvRk0meDdDc5n2xOT5Yi3X0anG42e4/fngM4xyaZlLiB4H3V9Wx\nuSpTV4djquprwD8DL2S66vAi4MoknwNuB34myfuZrhoAUFWfb75/Cfh7em+Zj7wOkxwGu4ENSS5K\nciZwFb2JbqvZsYl8cOJEvp3Aa5KcmeR5wAbgwar6AvD1JJel9+f0a1lg8t+kasa80KTEaavD+cee\nDklyFvBzwENMUR2q6i1Vtb6qnge8BvhoVb2WKaoBQJKzk5zTbH8v8ArgUcZRh+X+5PxpPlV/Fb0n\nTObpTVZb9jEN8bXdDjwB/B+9z0auBs4DPgJ8FvgwcG5f+7c0ddgHvLLv+Aub/yzzwLuX+3WdYg2u\noPf+8B56v/weAmansA4/DvxHU4dHgDc2x6eqDn2v4aV852miqaoB8Lzm/8Ee4LFjv/fGUQcnnUmS\nJvptIknSmBgGkiTDQJJkGEiSMAwkSRgGkiQMA0kShoEkCfh/Nry60DICmNwAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "#Check that time is actually a delta\n", "x = np.fft.ifft(X)\n", "plot(x)\n", "xlim([-100,5000])\n", "\n", "save_wav(x, \"delta.wav\")" ] }, { "cell_type": "code", "execution_count": 41, "metadata": { "collapsed": false }, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "/home/emily/anaconda3/lib/python3.4/site-packages/numpy/core/numeric.py:462: ComplexWarning: Casting complex values to real discards the imaginary part\n", " return array(a, dtype, copy=False, order=order)\n" ] }, { "data": { "text/plain": [ "(0, 1.2)" ] }, "execution_count": 41, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXsAAAEACAYAAABS29YJAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAEBVJREFUeJzt3W+MZXV9x/H3h11IpWiR0GALGGzdNJBApaYLqVau1eqU\niDR9UNyisfgPTbZtUiNbfCCTaP/QpKk1NHRVJNY08EBbu0b+qMWbGKLgWhYW3QU2SLILgtpSa/VB\nl/DtgznMXqbL3GHn3J3Z83u/ksk9f36c+53vbD5z5vzOuaSqkCQN23FrXYAkafYMe0lqgGEvSQ0w\n7CWpAYa9JDXAsJekBkwN+ySfSvJEkt3Psf/yJPcmuS/JnUnO679MSdJqrOTM/kZgbpn9DwOvqarz\ngA8DH++jMElSf6aGfVV9DXhymf1fr6ofdat3AWf0VJskqSd9X7N/J3BLz8eUJK3Sxr4OlOS1wDuA\nV/V1TElSP3oJ+25S9hPAXFUd9pJPEj+ER5KOQFVltcdY9WWcJC8F/hl4a1XtW25sVflVxTXXXLPm\nNayXL3thL+zF8l99mXpmn+Qm4CLg1CT7gWuA47vw3g58CHgxcH0SgINVtbm3CiVJqzY17Ktqy5T9\n7wLe1VtFkqTe+QTtGhiNRmtdwrphLw6xF4fYi/6lz2tCy75RUkfrvSRpKJJQ62GCVpK0/hn2ktQA\nw16SGmDYS1IDDHtJaoBhL0kNMOwlqQGGvSQ1wLCXpAYY9pLUAMNekhpg2EtSAwx7SWqAYS9JDTDs\nJakBhr0kNcCwl6QGGPaS1ADDXpIaYNhLUgMMe0lqgGEvSQ0w7CWpAYa9JDXAsJekBhj2ktSAqWGf\n5FNJnkiye5kxH0vyUJJ7k5zfb4mSpNVayZn9jcDcc+1McjHw8qraBLwHuL6n2iRJPZka9lX1NeDJ\nZYa8Gfh0N/Yu4OQkp/VTniSpD31csz8d2D+xfgA4o4fjSpJ6srGn42TJeh1u0CWX9PRukmZiNIL3\nv3+tq9As9BH2jwJnTqyf0W37f046aX5x+dxzR5x77qiHt5fUh3vvhS9+0bBfa+PxmPF43PtxU3XY\nk/BnD0rOAr5QVeceZt/FwNaqujjJhcBHq+rCw4yrlbyXpLVxxx3wkY8svGr9SEJVLb168rxNPbNP\nchNwEXBqkv3ANcDxAFW1vapuSXJxkn3AT4ArVluUJKlfU8O+qrasYMzWfsqRJM2CT9BKUgMMe0mL\nnFYbLsNeEgBZ9RSg1jPDXpIaYNhLUgMMe0lqgGEvaZETtMNl2EsCnKAdOsNekhpg2EtSAwx7SWqA\nYS9JDTDsJS3ybpzhMuwlAd6NM3SGvSQ1wLCXpAYY9pLUAMNe0iInaIfLsJcEOEE7dIa9JDXAsJek\nBhj2ktQAw16SGmDYSwIWJmi9G2e4DHtJaoBhL0kNMOwlqQGGvSQ1wLCXtMgJ2uGaGvZJ5pLsTfJQ\nkm2H2X9qktuS7Epyf5I/nEmlkmbKj0sYtmXDPskG4DpgDjgH2JLk7CXDtgL3VNUrgBHwN0k2zqBW\nSdIRmnZmvxnYV1WPVNVB4Gbg0iVjvge8qFt+EfAfVfVUv2VKklZj2hn46cD+ifUDwAVLxnwCuCPJ\nY8ALgd/vrzxJUh+mhf1Kpms+COyqqlGSXwa+nORXq+rHSwfOz88vLo9GI0aj0fMoVdKsOUG79sbj\nMePxuPfjppb56Sa5EJivqrlu/Wrg6aq6dmLMLcCfV9Wd3fq/AduqaueSY9Vy7yVpbd15J1x11cKr\n1o8kVNWqp8+nXbPfCWxKclaSE4DLgB1LxuwFXt8VdRrwK8DDqy1MktSfZS/jVNVTSbYCtwMbgBuq\nak+SK7v924G/AG5Mci8Lvzyuqqr/nHHdkqTnYeotklV1K3Drkm3bJ5Z/CFzSf2mSpL74BK0kNcCw\nlwT4efZDZ9hLUgMMe0lqgGEvSQ0w7CWpAYa9pEVO0A6XYS8J8PPsh86wl6QGGPaS1ADDXpIaYNhL\nUgMMe0mAH5cwdIa9JDXAsJekBhj2ktQAw16SGmDYSwKcoB06w16SGmDYS1IDDHtJaoBhL0kNMOwl\nAX7E8dAZ9pIWeTfOcBn2ktQAw16SGmDYS1IDDHtJasDUsE8yl2RvkoeSbHuOMaMk9yS5P8m49yol\nzZwflzBsG5fbmWQDcB3weuBR4JtJdlTVnokxJwN/D7yxqg4kOXWWBUuSnr9pZ/abgX1V9UhVHQRu\nBi5dMuYPgM9V1QGAqvph/2VKklZjWtifDuyfWD/QbZu0CTglyVeT7Ezytj4LlCSt3rKXcYCVXME7\nHvg14HXAicDXk3yjqh5aOnB+fn5xeTQaMRqNVlyoJLVgPB4zHo97P25qmRmZJBcC81U1161fDTxd\nVddOjNkGvKCq5rv1TwK3VdVnlxyrlnsvSWtr50543/vgm99c60o0KQlVteoPs5h2GWcnsCnJWUlO\nAC4DdiwZ86/Aq5NsSHIicAHwndUWJuno83xsuJa9jFNVTyXZCtwObABuqKo9Sa7s9m+vqr1JbgPu\nA54GPlFVhr0krSPLXsbp9Y28jCOtazt3wnvfu/Cq9eNoXcaRJA2AYS8J8PPsh86wl7TIK63DZdhL\nUgMMe0lqgGEvSQ0w7CWpAYa9JMC7cYbOsJe0yLtxhsuwl6QGGPaS1ADDXpIaYNhLApygHTrDXtIi\nJ2iHy7CXpAYY9pLUAMNekhpg2EsCnKAdOsNekhpg2Eta5N04w2XYS1IDDHtJaoBhL0kNMOwlAd6N\nM3SGvaRFTtAOl2EvSQ0w7CWpAYa9JDVgatgnmUuyN8lDSbYtM+7XkzyV5Pf6LVHS0eAE7bAtG/ZJ\nNgDXAXPAOcCWJGc/x7hrgdsA/8lIxygnaIdr2pn9ZmBfVT1SVQeBm4FLDzPuj4DPAj/ouT5JUg+m\nhf3pwP6J9QPdtkVJTmfhF8D13SbPDSRpndk4Zf9KgvujwJ9VVSUJy1zGmZ+fX1wejUaMRqMVHF6S\n2jEejxmPx70fN7XMRbokFwLzVTXXrV8NPF1V106MeZhDAX8q8FPg3VW1Y8mxarn3krS2du+Gyy+H\n++5b60o0KQlVteq50Gln9juBTUnOAh4DLgO2TA6oql+aKOpG4AtLg16StLaWDfuqeirJVuB2YANw\nQ1XtSXJlt3/7UahR0lHiH9/DNe3Mnqq6Fbh1ybbDhnxVXdFTXZKkHvkErSQ1wLCXBPgE7dAZ9pLU\nAMNe0iInaIfLsJekBhj2ktQAw16SGmDYSwK8G2foDHtJaoBhL2mRd+MMl2EvSQ0w7CWpAYa9JMAJ\n2qEz7CWpAYa9pEVO0A6XYS9JDTDsJakBhr0kwAnaoTPsJakBhr0kNcCwl7TIu3GGy7CXpAYY9pLU\nAMNeEuDdOENn2EtSAwx7SYucoB0uw16SGmDYS1IDVhT2SeaS7E3yUJJth9l/eZJ7k9yX5M4k5/Vf\nqqRZcoJ22KaGfZINwHXAHHAOsCXJ2UuGPQy8pqrOAz4MfLzvQiVJR24lZ/abgX1V9UhVHQRuBi6d\nHFBVX6+qH3WrdwFn9FumJGk1VhL2pwP7J9YPdNueyzuBW1ZTlKS14d04w7VxBWNW/ONP8lrgHcCr\njrgiSVLvVhL2jwJnTqyfycLZ/bN0k7KfAOaq6snDHWh+fn5xeTQaMRqNnkepkmbJCdr1YTweMx6P\nez9uasrfbUk2Ag8ArwMeA+4GtlTVnokxLwXuAN5aVd94juPUtPeStHYefBAuuQQeeGCtK9GkJFTV\nqn8VTz2zr6qnkmwFbgc2ADdU1Z4kV3b7twMfAl4MXJ+F04ODVbV5tcVJkvox9cy+tzfyzF5a1x58\nEN70poVXrR99ndn7BK0kNcCwl6QGGPaSAO/GGTrDXpIaYNhLUgMMe0mLvGFuuAx7SWqAYS8JcIJ2\n6Ax7SWqAYS9JDTDsJS1ygna4DHtJaoBhLwlwgnboDHtJaoBhL0kNMOwlLXKCdrgMe0lqgGEvCXCC\ndugMe0lqgGEvSQ0w7CWpAYa9pEXejTNchr0kNcCwlwR4N87QGfaS1ADDXpIaYNhLWuQE7XAZ9pLU\ngKlhn2Quyd4kDyXZ9hxjPtbtvzfJ+f2XKWnWnKAdtmXDPskG4DpgDjgH2JLk7CVjLgZeXlWbgPcA\n18+o1sEYj8drXcK6YS8OsReH2Iv+TTuz3wzsq6pHquogcDNw6ZIxbwY+DVBVdwEnJzmt90oHxH/I\nh9iLQ+zFIfaif9PC/nRg/8T6gW7btDFnrL40SVJfNk7Zv9K5+aVX+5zTl44xxx0Hjz8Ol1yy1pXA\nAw/At7611lUMS2qZe62SXAjMV9Vct3418HRVXTsx5h+AcVXd3K3vBS6qqieWHMtfAJJ0BKpq1dPn\n087sdwKbkpwFPAZcBmxZMmYHsBW4ufvl8F9Lg76vYiVJR2bZsK+qp5JsBW4HNgA3VNWeJFd2+7dX\n1S1JLk6yD/gJcMXMq5YkPS/LXsaRJA3DzJ+gXclDWce6JJ9K8kSS3RPbTkny5SQPJvlSkpMn9l3d\n9WNvkjdMbH9lkt3dvr872t9HH5KcmeSrSb6d5P4kf9xtb64fSX4myV1JdiX5TpK/7LY314tnJNmQ\n5J4kX+jWm+xFkkeS3Nf14u5u22x7UVUz+2Lh0s8+4CzgeGAXcPYs33MtvoDfBM4Hdk9s+2vgqm55\nG/BX3fI5XR+O7/qyj0N/Yd0NbO6WbwHm1vp7O4JevAR4Rbd8EvAAcHbD/Tixe90IfAN4dau96Gr/\nU+CfgB3depO9AL4LnLJk20x7Mesz+5U8lHXMq6qvAU8u2bz4sFn3+rvd8qXATVV1sKoeYeEHd0GS\nXwBeWFV3d+P+ceK/OWZU1eNVtatb/h9gDwvPYrTaj592iyewcPLzJI32IskZwMXAJzl0u3aTvegs\nvWllpr2Yddiv5KGsoTqtDt2V9ATwzFPFv8hCH57xTE+Wbn+UY7xX3V1c5wN30Wg/khyXZBcL3/NX\nq+rbNNoL4G+BDwBPT2xrtRcFfCXJziTv7rbNtBfTbr1cLWd/gaqq1p4zSHIS8DngT6rqx5n4lK2W\n+lFVTwOvSPJzwO1JXrtkfxO9SPIm4PtVdU+S0eHGtNKLzquq6ntJfh74cvd80qJZ9GLWZ/aPAmdO\nrJ/Js38TDdkTSV4C0P259f1u+9KenMFCTx7l2R8zcUa37ZiT5HgWgv4zVfX5bnOz/QCoqh8BXwRe\nSZu9+A3gzUm+C9wE/FaSz9BmL6iq73WvPwD+hYVL3jPtxazDfvGhrCQnsPBQ1o4Zv+d6sQN4e7f8\nduDzE9vfkuSEJC8DNgF3V9XjwH8nuSALp8Fvm/hvjhld7TcA36mqj07saq4fSU595o6KJC8Afhu4\nhwZ7UVUfrKozq+plwFuAO6rqbTTYiyQnJnlht/yzwBuA3cy6F0dh1vl3WLgjYx9w9VrPgs/oe7yJ\nhSeM/5eFOYorgFOArwAPAl8CTp4Y/8GuH3uBN05sf2X3Q98HfGytv68j7MWrWbgmu4uFYLuHhY/I\nbq4fwLnAv3e9uA/4QLe9uV4s6ctFHLobp7leAC/r/k3sAu5/Jhdn3QsfqpKkBvi/JZSkBhj2ktQA\nw16SGmDYS1IDDHtJaoBhL0kNMOwlqQGGvSQ14P8As3OmKHjLXhcAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "#Remove some of the zeros\n", "X[3500:5000] = np.zeros(1500, dtype=float)\n", "\n", "plot(X)\n", "ylim([0,1.2])" ] }, { "cell_type": "code", "execution_count": 42, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Writing to delta_filtered1.wav...\n" ] }, { "name": "stderr", "output_type": "stream", "text": [ "/home/emily/anaconda3/lib/python3.4/site-packages/numpy/core/numeric.py:462: ComplexWarning: Casting complex values to real discards the imaginary part\n", " return array(a, dtype, copy=False, order=order)\n", "/home/emily/anaconda3/lib/python3.4/site-packages/IPython/kernel/__main__.py:8: ComplexWarning: Casting complex values to real discards the imaginary part\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYMAAAEACAYAAABRQBpkAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAEKpJREFUeJzt3X+wXGV9x/H3JzcJNBWlSEUJYeJMaQudUizT4CBt19Jq\noK04/UOkrdPa2vJHaTt2RlHr1DAd2+qoVQfHSSl1/NGR6Uh/QDEgFq52OghmSgQ0gaSYGRJ+aCql\nvxBzud/+cXbJ5npvSjhn9+7Nvl8zO3vOs8+e59knZ+9nz9l9TlJVSJKm26rl7oAkafkZBpIkw0CS\nZBhIkjAMJEkYBpIkOgiDJJuT7EqyO8mVizx+cpKbk+xIcl+SX2/bpiSpW2kzzyDJDHA/8LPAfuDL\nwGVVtXOozhbguKp6e5KT+/VPqaq5Nh2XJHWn7ZHBJmBPVe2tqoPAdcAlC+o8Ajy/v/x84N8NAkma\nLKtbPn898NDQ+j7gvAV1rgFuS/IwcALwupZtSpI61vbI4NmcY3oHsKOqTgXOAT6S5ISW7UqSOtT2\nyGA/sGFofQPN0cGw84F3A1TVvyX5OvBDwPZBhSReIEmSnoOqShfbaXtksB04I8nGJGuBS4EbFtTZ\nRfMFM0lOoQmCBxduqKoWvb3rXe9a8rFpujkOjoHj4DgsvHWp1ZFBVc0luQK4BZgBrq2qnUku7z++\nFfgT4GNJvkITPm+tqm+17LckqUNtTxNRVduAbQvKtg4tHwB+sW07kqTRmfgZyL1eb7m7MBEcB8dg\nwHFoOA7dajXprLNOJDUJ/ZCklSQJNSFfIEuSjgGGgSTJMJAkGQaSJAwDSRKGgSQJw0CShGEgScIw\nkCRhGEiSMAwkSRgGkiQMA0kShoEkCcNAksSEh8HNN8ONNy53LyTp2DfRYbB9O9x553L3QpKOfRMd\nBpKk8TAMJEmGgSTJMJAkYRhIkjAMJEkYBpIkDANJEoaBJIkVEAZVy90DSTr2TXQYJMvdA0maDhMd\nBpKk8TAMJEmGgSTJMJAkYRhIkjAMJEl0EAZJNifZlWR3kiuXqNNLcneS+5LMtm1TktSt1W2enGQG\nuBr4WWA/8OUkN1TVzqE6JwIfAV5dVfuSnNymTUlS99oeGWwC9lTV3qo6CFwHXLKgzi8D11fVPoCq\nOtCyTUlSx9qGwXrgoaH1ff2yYWcAJyW5Pcn2JG84mga8HIUkjV6r00TAs/lTvQb4ceBCYB1wR5Iv\nVdXu4Upbtmx5ZrnX69Hr9bwchSQNmZ2dZXZ2diTbbhsG+4ENQ+sbaI4Ohj0EHKiqJ4Enk3wR+DFg\nyTCQJH23wQflgauuuqqzbbc9TbQdOCPJxiRrgUuBGxbU+QfggiQzSdYB5wFfa9muJKlDrY4Mqmou\nyRXALcAMcG1V7Uxyef/xrVW1K8nNwD3APHBNVRkGkjRB2p4moqq2AdsWlG1dsP4+4H1t25IkjYYz\nkCVJhoEkyTCQJGEYSJIwDCRJGAaSJFZAGHhtIkkavYkOA69NJEnjMdFhIEkaD8NAkmQYSJIMA0kS\nhoEkCcNAkoRhIEnCMJAkYRhIklgBYeDlKCRp9CY6DLwchSSNx0SHgSRpPAwDSZJhIEkyDCRJGAaS\nJAwDSRKGgSQJw0CShGEgSWIFhIGXo5Ck0ZvoMPByFJI0HhMdBpKk8TAMJEmGgSTJMJAkYRhIkugg\nDJJsTrIrye4kVx6h3k8kmUvyS23blCR1q1UYJJkBrgY2A2cBlyU5c4l67wFuBvzBqCRNmLZHBpuA\nPVW1t6oOAtcBlyxS73eBzwDfbNmeJGkE2obBeuChofV9/bJnJFlPExAf7Rc5p1iSJkzbMHg2f9g/\nCLytqormFNFRnSbychSSNHqrWz5/P7BhaH0DzdHBsHOB69JcW+Jk4KIkB6vqhuFKW7ZseWa51+vR\n6/W8HIUkDZmdnWV2dnYk2061+OidZDVwP3Ah8DBwF3BZVe1cov7HgBur6m8XlNdi/Xjve+HAgeZe\nknS4JFRVJx+bWx0ZVNVckiuAW4AZ4Nqq2pnk8v7jWzvooyRpxNqeJqKqtgHbFpQtGgJV9ca27UmS\nuucMZEmSYSBJMgwkSRgGkiQMA0kShoEkiRUQBl6OQpJGb6LDwMtRSNJ4THQYSJLGwzCQJBkGkiTD\nQJKEYSBJwjCQJGEYSJIwDCRJGAaSJAwDSRIrIAy8NpEkjd5Eh4HXJpKk8ZjoMJAkjYdhIEkyDCRJ\nhoEkCcNAkoRhIEnCMJAkYRhIkjAMJEmsgDDwchSSNHoTHQZejkKSxmOiw0CSNB6GgSTJMJAkGQaS\nJAwDSRIdhEGSzUl2Jdmd5MpFHv+VJF9Jck+Sf0lydts2JUndahUGSWaAq4HNwFnAZUnOXFDtQeCn\nqups4I+Bv2jTpiSpe22PDDYBe6pqb1UdBK4DLhmuUFV3VNUT/dU7gdNatilJ6ljbMFgPPDS0vq9f\ntpTfBD7bsk1JUsdWt3z+s75YRJJXAr8BvOKoGvByFJI0cm3DYD+wYWh9A83RwWH6XxpfA2yuqscX\n29CWLVueWe71evR6PS9HIUlDZmdnmZ2dHcm2Uy0+eidZDdwPXAg8DNwFXFZVO4fqnA7cBvxqVX1p\nie3UYv34wAdg377mXpJ0uCRUVScfm1sdGVTVXJIrgFuAGeDaqtqZ5PL+41uBPwK+D/homo/6B6tq\nU7tuS5K61PY0EVW1Ddi2oGzr0PKbgDe1bUeSNDrOQJYkGQaSJMNAkoRhIEnCMJAkYRhIklgBYeDl\nKCRp9CY6DLwchSSNx0SHgSRpPAwDSZJhIEkyDCRJGAaSJAwDSRKGgSQJw0CShGEgSWIFhIGXo5Ck\n0ZvoMPByFJI0HhMdBpKk8TAMJEmGgSTJMJAkYRhIkjAMJEkYBpIkDANJEoaBJAnDQJLEhIfB9dfD\nhz603L2QpMnzqU91u73UBFwJLkkt1o/BtYkmoIuSNFHOPx/uuCNUVSdXcZvoIwNJ0uLm57vdnmEg\nSStQ12dMDANJWoG6DgO/M5CkFaj5++h3BpKkDrUOgySbk+xKsjvJlUvU+XD/8a8kednRtrF7d9te\nSpKOpFUYJJkBrgY2A2cBlyU5c0Gdi4EfqKozgN8GPnq07dx0U5teStKx5dFHu99m2yODTcCeqtpb\nVQeB64BLFtR5DfBxgKq6EzgxySlH08ib3wzvfCc8/XTL3krSClYFX/gCvOQl3W97dcvnrwceGlrf\nB5z3LOqcBjx2NA29+93NbZxOOglOOw1e8AJYuxZWr4Y1a2BmprlftQpOPx1e9KKmrKqpU9U8dvzx\ncNxxzfL8fHM/N3do+0lTN2mCbu1aOHiweawK1q079FvimZnRvtZxfEk/rh8CHCuv5VhpAxbffwf7\nPzT7/eB9khx6vwzWB+8TaN5X8/NN2eD5g/ff/HyzraQp+/a3m/pzc8376zvfae7n5prb6tVw4MCh\n91kVPPkk7NrVbOfpp5vb/HxTf3A/NwePPAIPPtj8fXj88dGP4ai1DYNnuyst/Lb7u563ZcuWZ5Z7\nvR69Xu85d6or3/pWcxsY7GCDQDj+eHjhC5uddu3aZqc57rhDO966dU2YDHbk+fnmucO/khrs6KtW\nNTvqmjWH2nryycN38lFLJ79JWP42xtWObTw78/OLH9UP7/8zM4cHwML7hR+cVq1qlletasqHy1b3\n/6oNyp5+uil76qnm/fnUU4fanJtr/pAPPqQN3mf33NM8bxAAg1BYLDzHGwSz/Vv32obBfmDD0PoG\nmk/+R6pzWr/sMMNhsJT3vx9e+1rYuLH5R5akaXHgANx2W49LL+0NlV7V2fZbzTNIshq4H7gQeBi4\nC7isqnYO1bkYuKKqLk7ycuCDVfXyBds54jwDcK6BJA1s2wYXXwxdzjNodWRQVXNJrgBuAWaAa6tq\nZ5LL+49vrarPJrk4yR7gf4A3Hm07DzzQppeSdGy56KLut+kMZElagZyBLEnqnGEgSTIMJEmGgSSJ\nCQ+Ds89e7h5I0mT6xCe63d5Eh4EkaXFdT7w1DCRpBer6MiKGgSTJMJAkTXgYOPNYkhbnaSJJUucM\nA0lagabqyMDTRJI0HhMdBpKk8TAMJGkF8jSRJKlzEx0GkqTxMAwkaQWaqtNEkqTxMAwkaQWaqiMD\nv0CWpPGY6DCQJI2HYSBJK5CniSRJnZvoMJAkjYdhIEkrkKeJJEmdm+gwkCQtbqqODNatW+4eSNJ0\nWL3cHTiSm26CJ55Y7l5I0uQ555xut5eagBPzSWoS+iFJK0kSqqqTE0YTfZpIkjQehoEkyTCQJBkG\nkiRahEGSk5LcmuSBJJ9LcuIidTYkuT3JV5Pcl+T32nVXkjQKbY4M3gbcWlU/CPxTf32hg8Cbq+pH\ngJcDv5PkzKNpZHZ2tkUXjx2Og2Mw4Dg0HIdutQmD1wAf7y9/HHjtwgpV9WhV7egv/zewEzj1aBrx\nH7zhODgGA45Dw3HoVpswOKWqHusvPwaccqTKSTYCLwPubNGmJGkEjjgDOcmtwIsXeegPh1eqqpIs\nOWssyfOAzwC/3z9CkCRNkOc8AznJLqBXVY8meQlwe1X98CL11gD/CGyrqg8usS2nH0vSc9DVDOQ2\n1ya6Afg14D39+79fWCFJgGuBry0VBNDdi5EkPTdtjgxOAv4GOB3YC7yuqv4jyanANVX180kuAL4I\n3AMMGnp7Vd3cuueSpM5MxIXqJEnLa6JnICfZnGRXkt1Jrlzu/nQpyV8leSzJvUNlS07kS/L2/jjs\nSvKqofJzk9zbf+xD434dbSw1KXEKx+H4JHcm2ZHka0n+tF8+VeMAkGQmyd1JbuyvT+MY7E1yT38c\n7uqXjX4cqmoib8AMsAfYCKwBdgBnLne/Onx9P0nzU9t7h8reC7y1v3wl8Gf95bP6r39Nfzz2cOio\n7i5gU3/5s8Dm5X5tRzEGLwbO6S8/D7gfOHPaxqHf53X9+9XAl4ALpnQc/gD4a+CG/vo0jsHXgZMW\nlI18HCb5yGATsKeq9lbVQeA64JJl7lNnquqfgccXFC81ke8S4NNVdbCq9tL8g5/X/xXXCVV1V7/e\nJ1hk8t+kqsUnJa5nysYBoKr+t7+4luaD0ONM2TgkOQ24GPhLYPCjkqkagyELf1Qz8nGY5DBYDzw0\ntL6vX3YsW2oi36k0r39gMBYLy/ezQsdowaTEqRuHJKuS7KB5vbdX1VeZvnH4c+AtwPxQ2bSNATQ/\ntvl8ku1JfqtfNvJxmOT/9nKqv9muOvJEvmNJf1Li9TSTEv8rQ//T97SMQ1XNA+ckeQFwS5JXLnj8\nmB6HJL8AfKOq7k7SW6zOsT4GQ15RVY8k+X7g1v6crmeMahwm+chgP7BhaH0DhyfdseixJC8G6B/m\nfaNfvnAsTqMZi/395eHy/WPoZ2f6kxKvBz5ZVYO5KlM3DgNV9QRwE3Au0zUO5wOvSfJ14NPAzyT5\nJNM1BgBU1SP9+28Cf0dzynzk4zDJYbAdOCPJxiRrgUtpJrodywYT+eDwiXw3AK9PsjbJS4EzgLuq\n6lHgP5Ocl+bj9BtYZPLfpOr3ebFJidM2DicPfh2S5HuAnwPuZorGoareUVUbquqlwOuB26rqDUzR\nGAAkWZfkhP7y9wKvAu5lHOOw3N+c/z/fql9E8wuTPTST1Za9Tx2+tk8DDwPfoflu5I3AScDngQeA\nzwEnDtV/R38cdgGvHio/t7+z7AE+vNyv6yjH4AKa88M7aP743Q1snsJx+FHgX/vjcA/wln75VI3D\n0Gv4aQ79mmiqxgB4aX8/2AHcN/i7N45xcNKZJGmiTxNJksbEMJAkGQaSJMNAkoRhIEnCMJAkYRhI\nkjAMJEnA/wFI+2+RdDGu5QAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "#Look at time domain, there is some spread\n", "x_new1 = np.fft.ifft(X)\n", "plot(x_new1)\n", "xlim([-100,5000])\n", "\n", "save_wav(x_new1, \"delta_filtered1.wav\")" ] }, { "cell_type": "code", "execution_count": 43, "metadata": { "collapsed": false }, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "/home/emily/anaconda3/lib/python3.4/site-packages/numpy/core/numeric.py:462: ComplexWarning: Casting complex values to real discards the imaginary part\n", " return array(a, dtype, copy=False, order=order)\n" ] }, { "data": { "text/plain": [ "(0, 1.2)" ] }, "execution_count": 43, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXsAAAEACAYAAABS29YJAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAERtJREFUeJzt3W2sZVddx/Hvj5k2ioKV1IB2RovamDYRQeJAAOUgKkNj\nqfGFOApR5KFvRkwkWOGFnYnPRiOSmjpiaZBI+0J8KKGlYtoTCYHWamdaoFM6QpPpAAW0oqKJ0/D3\nxTlz5/Q4vee2d+/uPWd9P8nN3Q/r7rXumtPfWXetvU9TVUiS1ttThm6AJKl/hr0kNcCwl6QGGPaS\n1ADDXpIaYNhLUgNWhn2Sdyd5KMk9j3H+Z5McSXJ3ko8meU73zZQkbcdWRvbXAXs3Of8Z4Ieq6jnA\nrwN/2kXDJEndWRn2VfUR4OFNzn+sqr4y370d2NVR2yRJHel6zv71wE0dX1OStE07u7pQkpcBvwC8\nuKtrSpK60UnYzxdl3wXsraozTvkk8UN4JOkJqKps9xrbnsZJ8u3AXwGvqapjm5WtKr+quOqqqwZv\nw1i+7Av7wr7Y/KsrK0f2Sa4HXgqcn+Q4cBVwzjy8DwG/BnwzcE0SgJNVtaezFkqStm1l2FfVvhXn\n3wC8obMWSZI65xO0A5hMJkM3YTTsi9Psi9Psi+6lyzmhTStK6smqS5LWRRJqDAu0kqTxM+wlqQGG\nvSQ1wLCXpAYY9pLUAMNekhpg2EtSAwx7SWqAYS9JDTDsJakBhr0kNcCwl6QGGPaS1ADDXpIaYNhL\nUgMMe0lqgGEvSQ0w7CWpAYa9JDXAsJekBhj2ktQAw16SGmDYS1IDDHtJaoBhL0kNMOwlqQErwz7J\nu5M8lOSeTcq8M8n9SY4keV63TZQkbddWRvbXAXsf62SSS4HvrqqLgDcB13TUNklSR1aGfVV9BHh4\nkyKvAt4zL3s7cF6SZ3bTPElSF7qYs78AOL6w/yCwq4PrSpI6srOj62Rpvzq6rtSE970Prr9+6FbA\nZAJvecvQrVAfugj7E8Duhf1d82P/z4EDBza2J5MJk8mkg+qls99tt8F3fAe84hXDteHIEfjgBw37\noU2nU6bTaefXTdXqQXiSC4EPVNX3nuHcpcD+qro0yQuBd1TVC89QrrZSl9SiN74R9uyZfR/KrbfC\nb/zG7LvGIwlVtTx78ritHNknuR54KXB+kuPAVcA5AFV1qKpuSnJpkmPAV4HXbbdRkqRurQz7qtq3\nhTL7u2mOJKkPPkErSQ0w7KURGMty1ljaoe4Z9tJIZNtLcGd3/eqXYS9JDTDsJakBhr0kNcCwl0Zg\nLAujY2mHumfYSyMx9ALp0PWrX4a9JDXAsJekBhj2ktQAw14aARdG1TfDXtIG33TWl2EvjcTQd8MM\nXb/6ZdhLUgMMe0lqgGEvSQ0w7KURGMvC6Fjaoe4Z9tJIDL1AOnT96pdhL0kNMOwlqQGGvSQ1wLCX\nRsCFUfXNsJdGYugF0sQ3nXVm2EtSAwx7SWqAYS9JDTDspRFwrlx9M+wlbfBNZ32tDPske5McTXJ/\nkivPcP78JB9KcjjJJ5L8fC8tldbcGO7G0fraNOyT7ACuBvYClwD7kly8VGw/cFdVPReYAH+QZGcP\nbZUkPUGrRvZ7gGNV9UBVnQRuAC5fKvN54Onz7acD/1pVj3TbTEnSdqwagV8AHF/YfxB4wVKZdwG3\nJvkc8DTgp7prntQG58rVt1Vhv5WX4NuBw1U1SfJdwIeTfF9V/edywQMHDmxsTyYTJpPJ42iqpL75\npjO86XTKdDrt/Lqrwv4EsHthfzez0f2iFwG/CVBV/5Lks8D3AHcuX2wx7CU92tALpEPXr5nlgfDB\ngwc7ue6qOfs7gYuSXJjkXODVwI1LZY4CPwKQ5JnMgv4znbROktSJTUf2VfVIkv3ALcAO4NqqujfJ\nFfPzh4DfAq5LcoTZm8evVNW/9dxuSdLjsPIWyaq6Gbh56dihhe0vA5d13zRJUld8glYaARdG1TfD\nXhqJoRdI/Tz79WbYS1IDDHtJaoBhL0kNMOylEXCuXH0z7CVt8E1nfRn20kiM4W4crS/DXpIaYNhL\nUgMMe2kEnCtX3wx7SWqAYS+NxNALpH5cwnoz7CWpAYa9JDXAsJekBhj20gg4V66+GfbSSLhAqz4Z\n9pLUAMNekhpg2EtSAwx7aQScK1ffDHtpJMawQKv1ZdhL2uBfGOvLsJekBhj2ktQAw14aAadP1DfD\nXpIasDLsk+xNcjTJ/UmufIwykyR3JflEkmnnrZQaMPTdMH5cwnrbudnJJDuAq4EfAU4A/5jkxqq6\nd6HMecAfA6+oqgeTnN9ngyVJj9+qkf0e4FhVPVBVJ4EbgMuXyvwM8P6qehCgqr7cfTMlSduxKuwv\nAI4v7D84P7boIuAZSW5LcmeS13bZQEnS9m06jQNsZQbvHOD7gZcDTwU+luTjVXX/csEDBw5sbE8m\nEyaTyZYbKq0z58p1ynQ6ZTqddn7dVWF/Ati9sL+b2eh+0XHgy1X1P8D/JPkH4PuATcNe0qONYYFW\nw1seCB88eLCT666axrkTuCjJhUnOBV4N3LhU5m+BlyTZkeSpwAuAT3XSOklPKv/CWF+bjuyr6pEk\n+4FbgB3AtVV1b5Ir5ucPVdXRJB8C7ga+Bryrqgx7SRqRVdM4VNXNwM1Lxw4t7f8+8PvdNk2S1BWf\noJVGwOkT9c2wl0Zi6AXSoetXvwx7SRv8C2N9GfaS1ADDXpIaYNhLI+D0ifpm2EtSAwx7aSSGvhtm\n6PrVL8Ne0gank9aXYS9JDTDspRFwRK2+GfaS1ADDXhqJoRdIh65f/TLsJW1wOml9GfaS1ADDXpIa\nYNhLI+D0ifpm2EsjMfQC6dD1q1+GvSQ1wLCXtMHppPVl2EtSAwx7aQQcUatvhr0kNcCwl0Zi6Lth\nhq5f/TLsJW1wOml9GfaS1ADDXhoBR9Tqm2EvSQ1YGfZJ9iY5muT+JFduUu4HkjyS5Ce7baLUhqEX\nSIeuX/3aNOyT7ACuBvYClwD7klz8GOV+F/gQ4EtGOks5nbS+Vo3s9wDHquqBqjoJ3ABcfoZyvwj8\nJfCljtsnSerAqrC/ADi+sP/g/NiGJBcwewO4Zn7IsYEkjczOFee3EtzvAH61qipJ2GQa58CBAxvb\nk8mEyWSyhctL68/pE50ynU6ZTqedX3dV2J8Adi/s72Y2ul/0fOCGWc5zPvDKJCer6sbliy2GvaRH\nG3qBdOj6NbM8ED548GAn110V9ncCFyW5EPgc8Gpg32KBqvrOU9tJrgM+cKaglyQNZ9Owr6pHkuwH\nbgF2ANdW1b1JrpifP/QktFHSk8TppPW1amRPVd0M3Lx07IwhX1Wv66hdkqQO+QStNAKOqNU3w14a\niaEXSIeuX/0y7CWpAYa9pA1OJ60vw16SGmDYSyPgiFp9M+wlqQGGvTQSQ98NM3T96pdhL0kNMOwl\nbXDtYH0Z9pLUAMNeGgFH1OqbYS+NxNALpEPXr34Z9pLUAMNe0gank9aXYS9JDTDspRFwRK2+GfbS\nSAy9QDp0/eqXYS9JDTDsJakBhr2kDa4drC/DXhoBQ1Z9M+wlqQGGvTQSQ98NM3T96pdhL0kNMOwl\nbXDtYH0Z9tIIGLLqm2EvSQ3YUtgn2ZvkaJL7k1x5hvM/m+RIkruTfDTJc7pvqrTehl4gHbp+9Wtl\n2CfZAVwN7AUuAfYluXip2GeAH6qq5wC/Dvxp1w2VJD1xWxnZ7wGOVdUDVXUSuAG4fLFAVX2sqr4y\n370d2NVtMyVJ27GVsL8AOL6w/+D82GN5PXDTdholaRguFK+vnVsos+V//iQvA34BePETbpHUIENW\nfdtK2J8Adi/s72Y2un+U+aLsu4C9VfXwmS504MCBje3JZMJkMnkcTZXW29ALpEPXr5npdMp0Ou38\nuqkVQ4okO4H7gJcDnwPuAPZV1b0LZb4duBV4TVV9/DGuU6vqklr1ylfCm988+z6UT38aLrsM7rtv\nuDbo/0tCVW37rXjlyL6qHkmyH7gF2AFcW1X3Jrlifv4Q8GvANwPXZDY8OFlVe7bbOElSN7YyjUNV\n3QzcvHTs0ML2G4A3dNs0SU82//heXz5BK42AIau+GfaS1ADDXhqJoe+GGbp+9cuwl6QGGPaS1ADD\nXhqBsSzQjqUd6p5hL0kNMOylkRh6gXTo+tUvw16SGmDYS1IDDHtJG1ygXV+GvTQChqz6ZthLIzH0\nAunQ9atfhr0kNcCwl6QGGPaSNrh2sL4Me2kEDFn1zbCXRmLoBdKh61e/DHtJaoBhL0kNMOwlqQGG\nvTQCY1mgHUs71D3DXpIaYNhLIzH03TBD169+GfaS1ADDXpIaYNhL2uAC7foy7KURMGTVt5Vhn2Rv\nkqNJ7k9y5WOUeef8/JEkz+u+mdL6G3qBdOj61a9Nwz7JDuBqYC9wCbAvycVLZS4FvruqLgLeBFzT\nU1vXxnQ6HboJo2FfnHb48HToJoyGr4vurRrZ7wGOVdUDVXUSuAG4fKnMq4D3AFTV7cB5SZ7ZeUvX\niC/k0+yL044cmQ7dhNHwddG9VWF/AXB8Yf/B+bFVZXZtv2mSpK7sXHF+q8tGy7N9Z/y5yy7b4tXW\n3H33wT/909CtGAf7YubIEdg18BDpKU+BL3xhHP+d+rroXmqT2wCSvBA4UFV75/tvA75WVb+7UOZP\ngGlV3TDfPwq8tKoeWrqW9xtI0hNQVdtePl81sr8TuCjJhcDngFcD+5bK3AjsB26Yvzn8+3LQd9VY\nSdITs2nYV9UjSfYDtwA7gGur6t4kV8zPH6qqm5JcmuQY8FXgdb23WpL0uGw6jSNJWg+9P0G7lYey\nznZJ3p3koST3LBx7RpIPJ/l0kr9Lct7CubfN++Nokh9bOP78JPfMz/3Rk/17dCHJ7iS3Jflkkk8k\nefP8eHP9keTrktye5HCSTyX57fnx5vrilCQ7ktyV5APz/Sb7IskDSe6e98Ud82P99kVV9fbFbOrn\nGHAhcA5wGLi4zzqH+AJ+EHgecM/Csd8DfmW+fSXwO/PtS+b9cM68X45x+i+sO4A98+2bgL1D/25P\noC+eBTx3vv2NwH3AxQ33x1Pn33cCHwde0mpfzNv+y8BfADfO95vsC+CzwDOWjvXaF32P7LfyUNZZ\nr6o+Ajy8dHjjYbP595+Yb18OXF9VJ6vqAWb/cC9I8q3A06rqjnm5P1/4mbNGVX2hqg7Pt/8LuJfZ\nsxit9sd/zzfPZTb4eZhG+yLJLuBS4M84fbt2k30xt3zTSq990XfYb+WhrHX1zDp9V9JDwKmnir+N\nWT+ccqpPlo+f4Czvq/ldXM8DbqfR/kjylCSHmf3Ot1XVJ2m0L4A/BN4KfG3hWKt9UcDfJ7kzyRvn\nx3rti1W3Xm6Xq79AVVVrzxkk+Ubg/cAvVdV/ZuFTtlrqj6r6GvDcJN8E3JLkZUvnm+iLJD8OfLGq\n7koyOVOZVvpi7sVV9fkk3wJ8eP580oY++qLvkf0JYPfC/m4e/U60zh5K8iyA+Z9bX5wfX+6TXcz6\n5ASP/piJXfNjZ50k5zAL+vdW1d/MDzfbHwBV9RXgg8DzabMvXgS8KslngeuBH07yXtrsC6rq8/Pv\nXwL+mtmUd6990XfYbzyUleRcZg9l3dhznWNxI/Bz8+2fA/5m4fhPJzk3ybOBi4A7quoLwH8keUFm\nw+DXLvzMWWPe9muBT1XVOxZONdcfSc4/dUdFkq8HfhS4iwb7oqreXlW7q+rZwE8Dt1bVa2mwL5I8\nNcnT5tvfAPwYcA9998WTsOr8SmZ3ZBwD3jb0KnhPv+P1zJ4w/l9maxSvA54B/D3waeDvgPMWyr99\n3h9HgVcsHH/+/B/9GPDOoX+vJ9gXL2E2J3uYWbDdxewjspvrD+B7gX+e98XdwFvnx5vri6V+eSmn\n78Zpri+AZ89fE4eBT5zKxb77woeqJKkB/m8JJakBhr0kNcCwl6QGGPaS1ADDXpIaYNhLUgMMe0lq\ngGEvSQ34P9knxWP9F8bwAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "#Zero out more frequencies\n", "X[0:3000] = np.zeros(3000, dtype=float)\n", "\n", "plot(X)\n", "ylim([0,1.2])" ] }, { "cell_type": "code", "execution_count": 44, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Writing to delta_filtered2.wav...\n" ] }, { "name": "stderr", "output_type": "stream", "text": [ "/home/emily/anaconda3/lib/python3.4/site-packages/numpy/core/numeric.py:462: ComplexWarning: Casting complex values to real discards the imaginary part\n", " return array(a, dtype, copy=False, order=order)\n", "/home/emily/anaconda3/lib/python3.4/site-packages/IPython/kernel/__main__.py:8: ComplexWarning: Casting complex values to real discards the imaginary part\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYoAAAEACAYAAACtVTGuAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAFCFJREFUeJzt3X+MXWWdx/HPZ2Y6lNIqsGArpW4x1BQSE4imNOqGu+uP\nLWiA/UfkD2Ma45IYXJM1bv2RrGX3D9dN3EVDYoiLigYh6iopWV1ahInGsJQKlaIdpIYqlF9Fi1IK\nZX5894/nHO7pMPM4M+feub/er2Ryfj3nOc995nY+57n3nFNHhAAAmMtQpxsAAOhuBAUAIIugAABk\nERQAgCyCAgCQRVAAALJqB4XtLbbHbT9ie9ss2zfavsf2S7Y/MWPbQdsP2n7A9u66bQEAtN5InZ1t\nD0u6XtK7JB2SdJ/tHRGxv1Ls95I+JumKWaoISY2I+EOddgAA2qfuiGKTpAMRcTAiJiTdKunyaoGI\nOBwReyRNzFGHa7YBANBGdYNiraTHKsuPF+vmKyTdaXuP7Y/UbAsAoA1qffSk9Ie+jrdHxJO2z5S0\ny/Z4RPy0Zp0AgBaqGxSHJK2rLK9TGlXMS0Q8WUwP2/6B0kdZJwSFbR5GBQCLEBEt+Wi/7kdPeyRt\nsL3e9qikKyXtmKPsCQ22vcL2qmL+FEnvkbRvth0jYtafz33uc3NuG5Qf+oB+oB/og9l+WqnWiCIi\nJm1fI+kOScOSboyI/bavLrbfYHuNpPskvUbStO2PSzpf0uskfd922Y6bI2JnnfYAAFqv7kdPiogf\nSfrRjHU3VOaf0okfT5WOSrqg7vEBAO3V03dmNxqNTjeh4+iDhH5I6Af6oB3c6s+yWs12dHsbAaDb\n2FZ0yZfZAIA+R1AAALIICgBAFkEBAMgiKAAAWQQFACCLoAAAZBEUAIAsggIAkEVQAACyCAoAQBZB\nAQDIIigAAFkEBQAgi6AAAGQRFACALIICAJBFUAAAsggKAEAWQQEAyCIoAABZBAUAIIugAABkERQA\ngCyCAgCQRVAAALIICgBAFkEBAMgiKAAAWQQFACCLoAAAZBEUAIAsggIAkEVQAACyageF7S22x20/\nYnvbLNs32r7H9ku2P7GQfQEAneeIWPzO9rCkhyW9S9IhSfdJuioi9lfKnCnpLyVdIelIRHxxvvsW\n5aJOGwFgENlWRLgVddUdUWySdCAiDkbEhKRbJV1eLRARhyNij6SJhe4LAOi8ukGxVtJjleXHi3Xt\n3lf33CO98MJ8SwMAFmuk5v51PhOa977bt29/Zb7RaKjRaGjrVum226SNG2u0AAD6xNjYmMbGxtpS\nd92gOCRpXWV5ndLIoKX7VoMCAPBq5Ul06dprr21Z3XU/etojaYPt9bZHJV0pacccZWd+qbKQfQEA\nHVJrRBERk7avkXSHpGFJN0bEfttXF9tvsL1G6Yqm10iatv1xSedHxNHZ9q3THgBA69W6PHYpzHV5\n7MaNfEcBAHPppstjAQB9jqAAAGQRFACArJ4Oii7/egUA+kLPBoVb8hUNAODP6dmgAAAsDYICAJBF\nUAAAsggKAEAWQQEAyCIoAABZPRsU4+PS977X6VYAQPe58MLW1tezDwW0pTe+UfrNbzrQKADoYuk+\nMx4KCABYIgQFACCLoAAAZBEUAIAsggIAkNXTQcETZAGg/Xo6KLr8yl4A6AsEBQAgq6eDAgDQfgQF\nACCLoAAAZPV0UHDVEwC0X08HBQCg/Xo6KLjqCQDar6eDAgDQfj0dFIwoAKD9ejooAADt19NBwVVP\nANB+PR0UAID2IygAAFkEBQAgi6AAAGTVDgrbW2yP237E9rY5yny52P4L2xdW1h+0/aDtB2zvrtsW\nAEDrjdTZ2fawpOslvUvSIUn32d4REfsrZS6VdG5EbLB9kaSvSNpcbA5JjYj4w2KOz30UANB+dUcU\nmyQdiIiDETEh6VZJl88oc5mkmyQpIu6VdKrt1ZXtXOQKAF2sblCslfRYZfnxYt18y4SkO23vsf2R\nmm0BALRBrY+elP7Qz8dco4Z3RMQTts+UtMv2eET8tGabAAAtVDcoDklaV1lepzRiyJU5u1iniHii\nmB62/QOlj7JeFRTbt29/Zb7RaKjRaNRsNgD0l7GxMY2NjbWlbkeNb4Rtj0h6WNI7JT0habekq2b5\nMvuaiLjU9mZJ10XEZtsrJA1HxPO2T5G0U9K1EbFzxjFitjba0vr10qOPLrr5ANCX0uONrIhoyXfA\ntUYUETFp+xpJd0galnRjROy3fXWx/YaI+KHtS20fkPSCpK3F7mskfd/pFY1IunlmSAAAOq/WiGIp\nMKIAgIVp9Yiip+/M7vKMA4C+0NNBAQBoP4ICAJDV00Hx1FOdbgEA9L+eDorjxzvdAgDoLlNTra+z\np4MCANB+BAUAIIugAIA+4jY8j5ugAIA+0o77ywgKAEAWQQEAyCIoAABZBAUA9BG+owAALDmCAgCQ\nRVAAALIICgDoI3xHAQBYcgQFACCLoAAAZBEUANBH+I4CALDkCAoAQFbPB0U7hlkA0Ku+9rXW1+no\n8r+0tmO2Npb/OcfkpDQ8vMSNAoAutXKl9MILkmRFREv+G6OeH1EAANqr54Pij3/sdAsAoHuk0URr\n9XxQfPvbnW4BAPS3ng8KAEB7ERQAgKyeD4r77+90CwCgO0xPt6fenr88VuJeCgCQpDvvlN797nKp\ndZfHEhQA0Cfe+lbp5z8vlwiKE4Li0Uel9euXrk0A0I18QiwQFDM6hFEFALQrKHr+y+zS73/f6RYA\nQOcsX96+umsHhe0ttsdtP2J72xxlvlxs/4XtCxey73ydcYZ0+HCdGgCgN334w9Lx4+2rv1ZQ2B6W\ndL2kLZLOl3SV7fNmlLlU0rkRsUHS30v6ynz3XajXvU66+GLp2LE6tQBAb7j33vRxUzueGFtVd0Sx\nSdKBiDgYEROSbpV0+Ywyl0m6SZIi4l5Jp9peM899F+wnP5FOOSV1ni1t3Srdfrv05JN8jwGgNx0/\nLu3bJ33pS9KGDc2/b5s3L83xR2ruv1bSY5XlxyVdNI8yayWdNY99a/vGN9JPHaedJi1bJq1Zk6an\nnSadfLK0YkUKpWXL0qN9h4fTdGgobRsaSuVs6aST0vLoaCq3bFlajpBGRtL6k05KN8yU09NOS49R\nHx1NyyMjqfzQ0In7Tk6mfV5+OU2PH0/1Hz+ePrcspy+/3JyOjqb9RkZS3WV9Q0Np+/R0Wq5+OWY3\n2zE5mY4xMZHqKus8fjy1YWIibS+PMTWV6p6eTvWUdZd1Dg2lMtW6yzpnTstHy5f7TU+n5eXLm8er\ntqvsk3Jqp2OVdZSvMSLtPzXVXC7bVW17tS+qxy/rnK3u6usu95mtP2ZOS3OtL9tXPW61z8v1M9tf\nNfO4Q0PNOqq/53Jari+PVb53ZtZTVbZj5u911ar0O6n+rl56qfm7mppKxyyPUe2Psp3le/2kk6QX\nX0zvg2PH0r+9Y8eay8uXp7pHR9OD86anU93l+2Fiovnen55O6158MR1neDjN26kOKR1zejrtU7b1\nueeko0elI0dSG595JpX93e/S+mefXdjfnm5RNyjme47ekm/eO+XIkTQtf+kAMEjqBsUhSesqy+uU\nRga5MmcXZZbNY19J0vbt21+ZbzQaajQai21vW61Ykc4uzjgjnV289rXpDGTlyjQ6GBpKIxApne1I\nzeXR0bR9ZCT9lKMOKU3t5nSk+K2NjDSX7bSP3TwDKstMT6f6p6bSWVc5AinP8KqjlbKO8mfVqubZ\nm9Ssf3q6eYZZjiDKUU151l+OKJYta57dVs+IZ07L+fJ41bPvVauaZ8ZlHdVRSDkSq45KqmfE1TPj\n6qhr5oinPKMsj1uOTqr7Vkcp5Yit+trLs+7qaEHKj5qOH0/viRdffPW0eozy7HtionlGXPZbdWRW\nnu1Xj3HsWNpWtmvmNKLZt2V/lNOZI4ipqea+UnPf8kdq1js1lV5D+b48erTZ9qGhdBYuNUdy5bQ8\nsy9HtzOXJyeb5SOa08nJE7e//PLc06mp5ihycjKti2iOHo4dS23805/S7+3IkdSnzz6bpocPp74t\nX0NnjRU/rVfrPgrbI5IelvROSU9I2i3pqojYXylzqaRrIuJS25slXRcRm+ezb7H/vO6jyNm8Wbri\nijQ991xp9ermH1sA6EYRKah++1vpoYeku+6Sbr45Be38dNENd7YvkXSdpGFJN0bE521fLUkRcUNR\npry66QVJWyPi/rn2naX+BQfFzp3V550AQP95/vl0Wex3vztXiS4KinZbSFD87GfS2962BI0CgC5R\n/Wj4RNyZ/Sp33UVIABg8s11l1vJj9MuIostfBgC01dGj6aKPJkYUJxgf73QLAKCzVq5sX919MaLo\n8pcAAEuivNQ7YUQBAJhh7iug6un5EcWmTenBWAAw6E68AooRxStWr+50CwCgOyzkRuSF6PmgaFfH\nAACSng+Kiy/udAsAoL/1fFC8732dbgEA9LeeDwoAQHv1fFCsXdvpFgBA9/joR1tfZ89fHtvlzQeA\nJTU+Lp13nsTlsQCAWbXjSlCCAgCQRVAAALIICgBAFkEBAMgiKAAAWQQFACCLoACAPsLlsQCAJUdQ\nAACyCAoAQBZBAQDIIigAoI/wZTYAIGt0tPV18phxAOgz6e8jjxkHACwRggIAkEVQAACyCAoAQBZB\nAQDIIigAAFkEBQAgi6AAAGQtOihsn257l+1f295p+9Q5ym2xPW77EdvbKuu3237c9gPFz5bFtgUA\n0D51RhSfkrQrIt4k6cfF8glsD0u6XtIWSedLusr2ecXmkPQfEXFh8fO/NdoCAGiTOkFxmaSbivmb\nJF0xS5lNkg5ExMGImJB0q6TLK9vb8PgqAEAr1QmK1RHxdDH/tKTVs5RZK+mxyvLjxbrSx2z/wvaN\nc310BQDorJHcRtu7JK2ZZdNnqwsREbZnezxf7pF9X5H0L8X8v0r6oqQPz1Zw+/btr8w3Gg01Go1M\ntQAweMbGxjQ2NtaWuhf99Fjb45IaEfGU7ddLujsiNs4os1nS9ojYUix/WtJ0RHxhRrn1km6PiDfP\nchyeHgsAC9BNT4/dIelDxfyHJN02S5k9kjbYXm97VNKVxX4qwqX0d5L21WgLAKBN6owoTpf0HUlv\nkHRQ0vsj4jnbZ0n6akS8tyh3iaTrJA1LujEiPl+s/6akC5Q+nnpU0tWV7zyqx2FEAQAL0OoRBf9x\nEQD0mW766AkAMAAICgBAVk8Hxcknd7oFAND/ejoozjmn0y0AgP7X00EBAGi/7J3Z3WzbNun88zvd\nCgDoPhdcIO3d27r6evbyWADA3GwujwUALBGCAgCQRVAAALIICgBAFkEBAMgiKAAAWQQFACCLoAAA\nZBEUAIAsggIAkEVQAACyCAoAQBZBAQDIIigAAFkEBQAgi6AAAGQRFACALIICAJBFUAAAsggKAEAW\nQQEAyCIoAABZBAUAIIugAABkERQAgCyCAgCQRVAAALIICgBA1qKDwvbptnfZ/rXtnbZPnaPc12w/\nbXvfYvYHAHRWnRHFpyTtiog3SfpxsTybr0vaUmP/OY2NjS10l75DHyT0Q0I/0AftUCcoLpN0UzF/\nk6QrZisUET+VdGSx++fwhqAPSvRDQj/QB+1QJyhWR8TTxfzTklYv8f4AgCUwkttoe5ekNbNs+mx1\nISLCdiy2EXX3BwC0jyMW9/fZ9rikRkQ8Zfv1ku6OiI1zlF0v6faIePNC9ydAAGBxIsKtqCc7ovgz\ndkj6kKQvFNPb2rF/q14oAGBx6owoTpf0HUlvkHRQ0vsj4jnbZ0n6akS8tyh3i6SLJf2FpGck/XNE\nfH2u/eu9HABAqy06KAAAg6En78y2vcX2uO1HbG/rdHtaabYbFHM3J9r+dNEP47bfU1n/Ftv7im1f\nWurXUZftdbbvtv1L2w/Z/odi/UD1he3ltu+1vdf2r2x/vlg/UP0gSbaHbT9g+/ZieRD74KDtB4t+\n2F2sa38/RERP/UgalnRA0npJyyTtlXRep9vVwtf3V5IulLSvsu7fJf1TMb9N0r8V8+cXr39Z0R8H\n1Bwl7pa0qZj/oaQtnX5tC+yHNZIuKOZXSnpY0nkD2hcriumIpP+T9I4B7Yd/lHSzpB3F8iD2waOS\nTp+xru390Isjik2SDkTEwYiYkHSrpMs73KaWidlvUJzr5sTLJd0SERMRcVDpjXBRcRXZqojYXZT7\nphZxQ2MnRcRTEbG3mD8qab+ktRrMvjhWzI4qnSgd0YD1g+2zJV0q6b8klRe4DFQfVMy8wKft/dCL\nQbFW0mOV5ceLdf1srpsTz1J6/aWyL2auP6Qe7qPi8uoLJd2rAewL20O29yq93rsj4pcavH74T0mf\nlDRdWTdofSBJIelO23tsf6RY1/Z+qHN5bKcM9LfvEYN1c6LtlZL+W9LHI+J5u3kyNSh9ERHTki6w\n/VpJd9j+6xnb+7ofbL9P0jMR8YDtxmxl+r0PKt4eEU/aPlPSruJ+tFe0qx96cURxSNK6yvI6nZiO\n/ehp22skqRg2PlOsn9kXZyv1xaFivrr+0BK0s6VsL1MKiW9FRHmfzUD2hSRFxB8l/Y+kt2iw+uFt\nki6z/aikWyT9je1vabD6QJIUEU8W08OSfqD0UXzb+6EXg2KPpA2219selXSl0s17/ay8OVE68ebE\nHZI+YHvU9jmSNkjaHRFPSfqT7YucTsE/qIXfENlRRbtvlPSriLiusmmg+sL2GeVVLLZPlvRuSQ9o\ngPohIj4TEesi4hxJH5B0V0R8UAPUB5Jke4XtVcX8KZLeI2mflqIfOv0t/iK/+b9E6SqYA5I+3en2\ntPi13SLpCUkvK30Xs1XS6ZLulPRrSTslnVop/5miH8Yl/W1l/VuKN9EBSV/u9OtaRD+8Q+nz6L1K\nfxgfUHpc/UD1haQ3S7q/6IcHJX2yWD9Q/VB5DReredXTQPWBpHOK98FeSQ+Vf/uWoh+44Q4AkNWL\nHz0BAJYQQQEAyCIoAABZBAUAIIugAABkERQAgCyCAgCQRVAAALL+H/LqEA54604DAAAAAElFTkSu\nQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "#Even more spread\n", "x_new2 = np.fft.ifft(X)\n", "plot(x_new2)\n", "xlim([-100,5000])\n", "\n", "save_wav(x_new2, \"delta_filtered2.wav\")" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "##Bandwidth in Frequency Domain" ] }, { "cell_type": "code", "execution_count": 99, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "[]" ] }, "execution_count": 99, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX4AAAENCAYAAAAIbA6TAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztfXuMJ9lV3nemp2em57GzL+/gfZg19trYTuLYSYwDIbQT\nO1ocxQRiYhwSCFhgERkhK0qABOJdoSQgWQghYHEQIPMPjiAEFuwFL4QBC2MnFn6CV/YClnfX9r53\ndrbn1d1z80f1ddf+pqru+e49p+r2r+tIo+6eOXP7+50699zvfPdWlYQQMNtss8022/6xA1MDmG22\n2WabbVybC/9ss8022z6zufDPNttss+0zmwv/bLPNNts+s7nwzzbbbLPtM5sL/2yzzTbbPrPiwi8i\nvyQiD4vIJwd8flpEPisiHxeRV5T+ztlmm2222fLNgvH/MoDb+/5RRF4P4IUhhNsAfC+Auwx+52yz\nzTbbbJlWXPhDCB8A8OSAyxsAvHvH98MArhaRU6W/d7bZZptttjwbQ+O/CcADrZ8fBHDzCL93ttlm\nm222Dhtrc1cWfp6fEzHbbLPNNpEdHOF3PATgltbPN+/83bNMRObFYLbZZpstw0IIi+R60MZg/HcD\n+A4AEJFXA3gqhPBwl+M99wQAAb/1WwEhDP+5//7G953vTPtubTW+b3xj2rd5aF3AgQM639e+tvHX\n+P7IjzS+n/vcsN873vEO/PqvN7733pse9xOfaHzvuivtu7HR+L7lLfpYnDql8335y/WxePvbG99H\nHknH4ld+pfH94AfT437oQ43vu9+d9n3sscb3B35AH4uXvUzne+ON+lh8z/c0vmfPpmPx8z/f+H7s\nY+lx/+APGt9f+7W07+c/3/j+8A/rY7G+rvM9dEgfize9qfG9dCkdi5/8ycb3M59Jj/vbv934vve9\nad/77mt8f+zH0r6XL/O1ZW1N5/u61zX+ly+nxuTN4jjnrwL4IIAXi8gDIvLdIvJWEXkrAIQQ3gfg\nr0TkfgDvAvDv+sZ6aKcPeOKJ9O+NPl/4Qtr36aebr888k/aNduSIzi/G/cKFtO/jjz/765DFz6WJ\nRRxPE4snn3z2V40dPqzzW11tvl6+nPaNn4uJhcaXiUXEoIlx/EyHDqV9AX3+AMBTTzVfNdckzhFN\nLJg5wsQi5rooOSYTizhH45wdMiYWXnkRr51m/ke7dEnnt7XFj621YqknhPBmhc/bNGM99ljzlbmQ\njz6a9o1JpEmmWMgPKJfEOFkfewy4ObFlfeaMHgcTi5igmlgwGGKCbm+nfYHdSfvUU8C119rh8MqL\niCF+HTKWPGxuNl8vXEgXvjaOW24Z9o2fjyEE1rGIv1tLHuJcCiG9WLTz4rrrhn2ZWHjNEQZDNA0x\nWsSxtqYfX2NV3bn7zDMNa9RO8EOHgI2NtC+T1BcvAisrwPnzuoIXcWgKQkykFI719fWsWFhiAHav\nhzapvWOhneCrq/Z5wcbiiSfs8zM3L7xiocFw+TJw9ixw8KCOuS5zLM6da3xFdMWfwcFaVYV/Y6Nh\nzVoGdtNNuiITmZT2Qp48CRw92iSsZmwGh+bzxaS+6SYdZg8MQDPeDTc0zF/Tnp45A9x4I3dNUqyq\nHYvYVqfG1cbi6ad1GIBmvJtv1mHY3m6Iw6lTtjjYWDz9tP56eMViY6PpeK6+mstPbSy0OJhYaDEA\nTY1gatbVVzfxOH/eFgdrVRX+WGjOnUv7bmzYT6yI4fhx4NgxHY5z5/Q4zpwBnvc8PY5Tp/QYbrjB\nB8OJE7pYbG01f667zmcxPnVKN1nY68FguPbahq1FGWcIw9GjTeyWMRZnzzb5psVw7Fgzp7Rz9XnP\n88HB1AstBqZm5cZiLvwtY4pdXO3Pnt3V8PtsY6O5MEeP6ovdtdfqcJw927A1bdH1iAWL4cQJXSzO\nn+eKHYvj1Cldmx5joenWWAwxL1I4YuHXTnAGx9mz+lhsbHBzhIlF3MNJLYIbG/piF0JeXngQRQ8M\nTOG/eLGRg57znH0g9Xix3HPnmmCvrqZ1xvYE1xY77QQ/dw64/nrd52OTmonFddc12FOLIBOLjY28\nWGjYmlcnyGA4e5aLBcPsvGLBsNzz5zkMHovg5mbTUV11FRcLhhB4xOKaa5rvU3IoM0dibTl6VIeD\nteoKv/ZCMhP8/PlmV3xtLR3EtryRwsFO8PPnG6ak2eRiF0EmFsePN5ttqURlZC+2jY2x0E4uD0Jw\n/nwzaa0XQbbYAU3OeRR+JhbM9fCIBTNPIw6vWMS80GDwiMW5c1wsWKuq8G9sNK2N9kI+5zm6RSKu\nnpogxlWZuZCaRSLi0E6ujY2GdWjZDBMLbUK1Y6FldppYbG01bayW2XnG4sSJ5rihVrfXLoIxbtbF\nLn4+bfdz/fX2hZ+JRTuHPGKh7aBjbbGORU690MynGAvtRjBrVRX+CxcaGUJb+K+5ptHCUsYk1MWL\nzQ1LzIU8fDiNY3u7KXgnT+ou5MWLTfLtlVjECaCJBTvBL15sPp92gnvHQtMJesdCuwhee60+FidO\nNN+nFkEmFswcYYvdhQs+cyTGTXs9Dh/mumKPvGCtusKvvZAbGw1j3N5On4nNmeAMm1n2YsewXK9Y\nMBOcKXZehMBrgnsTAk3RrSEWly83C9TVV+sXQSYWJ082RC11L08thIC16gq/tthF+cY6iPFCrq3p\nNnfX1qYv/OfONRPg0qW0Vu0Vi6hJ1hKLqQs/Gwsty13mWDC69uZmc1jj+HF+jqTMKxae9YK16gr/\nyZO6RI0B9yr8mnG9MMSxT5zQxeLSpWbclZV0m+4Vi0uX6omF1jeH5R4+nC4eXrGIz245elRXwC5d\narpitvCnDh/UEAsGQ8SREwtLHJ71grUqC7/2QrJB1ByNqiGpI46rrtLFwnsRPHRouqSOLf3x4/pY\nHDvWtOipNj03FhpCcOgQHwttwdVgiP6exc4zFpYYor9HLGIN0OC4dImLhfZASo5VV/i1F6eGhPLC\nEALPcj1joV0EPTDEcY8c8SMEy8xyT5zwkwDZWKR8vUhJxMF2gtY42Lw4cmSfFP7I7DY304nKBPHi\nxSaIjHxjfSEjhiNH0r5bW80Rw7U1+2LH4GAWwRowtP1ZHJZsmyEEXhjaY2vu2aghFpcu8Rg04+YQ\nqRwcmhh7xCLHqir8hw83BW91VXek7NAhrsViJQvrC+mBoRYcNWBoj83i0B5htO5+cjAwLJfZo/HA\nUQOG7e2mrhw+3EiHKQkwF4clOWIw5FhVhT8+t9wziNpCUwOGlRV9oi57LLQstwYcDBv1WgTj3sjq\n6kwIIgaROnBYY8ixqgr/4cO7X6dMqDju1EkdE3VqHGzcPDGkJMAaYpHL+C0Xn1j0tTm0zIQgYgDq\nwDEX/gVjGT8zuVZX986FjBiAOnBYb2jmYDhwoNGq45HGLoud0cGDyx0LLRnwLnY1kBIGA7Cc5CjH\nqir88Z2tTGIv44VsJ2rKP4TdhW1qecNzggNpHGNN8KljERfAoTvW47jA9Cx36u5nDMY/JSHIsaoK\nf7vYTZ1Qnpt4LFsb8o8tfdy8WrakZmKx7BO8LQGm/MeSN2ohaEMSIEMI2kRqrxCCHKuq8EfGn/qw\n8fk8OS295S49m9RsMgF7j+Uyi/ayx2LKRXCMWNRAjlZWGtIzdACC6X7iCaCVlb1DCHKsqsLfvjgp\nNnPo0PS79DlJrTmpw7JczwkeP9/ULBfYOyx3rFgMjc3EIoTdztFLZplaAtTGoj2f9gohyLEqC79l\nUkd/T+2QwRDb9KHOo/35LJO6jaOGpF5d1bXpe43lsgvm1Iy/fQLIMxZTSoDaWJQU/qm6nxyrqvBr\nN3drnOBaDEAaB8tyx5A3vLqfVJvevtaWhCDKhWxLv1dYLhMLT5a7X2KhrRe1SD0H7YfMNy/Gv7nJ\nJ9T2tk9SA1zht9zQzN248pJ6gF0cB3sy0ZPltjuwKW/NX+x+RIYxALax8GS5YzB+zXwaIxZM3k9d\n+Kti/ExhrIHljsX4rTY0t7eboqJlubGYe0k9gG0sxmrpPQiB5j4Fr/0OT5Y7BuO3JEd7MRY5VlXh\n10o9y8JyUxiA5U5qDQ6W5XrGgilKUy+CtbDcnFho7lOoIRbt04VesdgXhV8r9TAXcmuruTAHDvgl\n9epq+hWQ8wR/Ng7t0d32wsbud2h9PRZB78I/FSGIRCp+Pg9CoDkAUUMsIgbtpvjM+HssZ9LWwHLj\nhddoyhocNSR1G0dOUg+d1PGUesaIhWYR9O5+piIEcdGORGovSIBesWDmaRsHE4uVlebn1IMaWauq\n8LdZ4F6Z4N4JpfGtYRGMOFZW0q+A9JR6PGKxtdUsZPGGwSkZ/9TkiMnNOPYYi+CUjD/6euz9aHDk\nWFWFX5vUY8kb2glujaMGlttu6ZmktsZRQyyYzq7tHyVA7V2lrOyVYrlesdDO04jDYxFkcHgStBpi\nkWNVFn5mAqRe2sIGsP0Mcy2bscYRMQBcUqcwMHJTe2+ESeo4tmUscia4VrJgMWjJQ9SqrfMipzBa\n+i5isNSqvXAsLphWsdja4mORQyqXvvBrpR6m2LEsNxY8VuqxxLG1xcleXhji2XqW8U+FYwwMUxMC\nLY4aMAC7nYdGq/bC4TVP2xg0taV9GGTo2C6LI8eqKvyMfumR1PEBTVqWW8PkGgMDw3KnxDEWhikJ\ngRZHDaRkcexUwfMkR3GOpO6VyMWQmqfth78dPJh+WOS+KvxTM/72xNKeJZ56co2BgdnQ1OBgJaec\n+zumjkUbxxSFxrPYMSy3fa1TBc8LRw2EoI0hPqjR6vh3jlVV+Kdm/O2JxT7z3JPZ1YBBc0TTo9Dk\nTq6pMHThmLrQTNmNxoe/TY3DmygyGKz3fnKsysIfn1vSZzGhoq9HMmn82zg8C83UGFZWmmRN6bMe\n14RZBNub4p4d2JDv5cu7Rz+1OLRjMzjaTHuqbrSNwRPH6qoeRw0YrHHkWFWFvx3EoQLWDriXZsfi\n8EyoqTEA6TiPgWNKDIsSoMZ3ShzthduLHDEYPHGk5sgYsWAwWOPIsaoKv7Y9bk8Ar4k1JY7FyZXS\ny70xpMYOoekG4smNGiaXF4b4Gfv02VqK3RhzhMHggaMmcsSQEg3mfVn4a5jgU+LInVxTMf6IIUfL\nTcl6i5MrtRBrJ3juBnMKR9cE9yo0VrEYA0PE4bHJzBA0r43uVIxZUrmvCn9NSZ3C0T7z74ljylho\n2dpY8sZUiyBDCFjyMEb347n3M3Q6ZTEWXhvdNRBFS9nr8mVubyvHqir8Xoy/XRgvXeo/ncJM2vaZ\n/1wcfeaZ1DkYgOHEzil2zCZljuylKXbaTUqGrTHkoQSHNSFgMERfkXQnyHY/OThqiIWG8WtjEQ8p\nxA46hSPHqiz8ltp6+6RHLNR9p1NKW/pUgc5JEiapNawjRxdN+bPFLheHpZa7eNLDMhZM97O1lYfD\nsvvJHTflz7Lc+CpMFoclOWJisUg0trc5UsmQhyEcOVZV4c+Z4JabKkxS52wEM5g95A12QdF+PnZD\nk8HhtdHNsGdGssiJRW5HU3ssmEVwe/vZ+0SMdGK996PtGhe7n5UVm+6H7aBzrKrCr2X8uUkdx9Yy\nFK8JzjIJj1jknE6xaulLFjZLxp8zbsqfiUXXPtFU3U9OsUvhKCl2bDG3JEc5m8bR36ITZElljlVZ\n+L2SGrCTb9hiV8JQaoiFFbNbxMF2SlYS4OIEZ+QNK0IQj8BqWW4uM/eSN1I4SuQNL3LELGxsLNiF\njaktS834vSd49LeQb0rkjSkZvzahGMZfKvUw+qyHlruywuuzUzC73C7FckFh55NWqy6RTi3JERsL\n5vOVdD/VFX4RuV1E7hORz4rID3b8+7qInBGRj+78+ZG+sXIKI1vs2GJuMcEX/ZmJ6NXGpnCwjF+b\n1FFaap+GYgqYR/vPnk5hOyWLxXVxbEtyVCJ7WXWCObHwIEdMd1DC+KeWeg6mXfpNRFYA/AyA1wJ4\nCMD/E5G7QwifXnD9oxDCG1Lj5U7w2osd4CtveLTppZ2SttixBcyjGLT929e0z9eq++mK2/nzOsye\n3c/WVtP9RAmqb9wUjhKpx7ITZMlRDnlIYWbqxV6Qel4F4P4QwudCCJsA3gPgmzr8OlLoSstNanZV\n9pA3hi46e1Qt9/OxrSn7+faClssQAmY/h1kwcye4tSSTQ468jzxbsdwScuQl9TCfb2rGX1r4bwLw\nQOvnB3f+rm0BwNeKyMdF5H0i8tK+wXLa2BwdjpmIFknddVRt6qROjc0WxpIJzuizNcheXoyf1aqt\nGH+pxGkRi1KpxyMWOR30FLHIsSKpB01RT9mfAbglhHBORL4RwG8CeFGX41133YHrrgMefRQ4c2Yd\nwHrngIutaWTUBzqWsa0t4MiR3Z+tkrpkBWcLmFWxYxKVWSRKkpqRZLy03NTYpRLZFLEoXQQ3N589\nb3Iwl3Y/Gxv9mHO7O2ZxjfVkqLbkkiNWcm77nj59GqdPn+7/EAorLfwPAbil9fMtaFj/ly2EcLb1\n/T0i8nMicm0I4YnFwd7+9jvwVV8FfOYzwPvf3/9L28nXfqlBfBHIou/x47s/W+nPpcXOUt7IneBW\nspenfskugl5Sj1f3w8Yip4Axvin/UonMwnfR33I+9S2CXbXFixylxl1fX8f6+vqXf77zzjv7P1CP\nlUo9HwFwm4jcKiKHALwJwN1tBxE5JdKIHCLyKgDSVfSBsknrJd94THBWvtFO2tRDszxlLy/9kpW9\ncjql1NglE9yz2Hl0BykcY5GjoesXX+mYQ3gsZS8vclS91BNC2BKRtwH4PQArAH4xhPBpEXnrzr+/\nC8AbAXyfiGwBOAfg23rBKFdwL5llrBWc0Q6ZCR67n62tZ7/tJwdHyQRgtU5Gn035ahcJL8wlE5zt\n7qy6n1Ipy4ocaTFE6SXKL5p6kUsqU7FYVBOsyBFDFHOsVOpBCOEeAPcs/N27Wt//LICfVYEpWJUt\nJlfJJpcVS1ocm23T49hdhZ9loyWdUp8+m8Nycxj/wYO7N2V1HUss1Wc9Jjhzra0lQI9YsOQoV0Ji\n6sVUcmip1GPN+Ku6c5dh/LkbJZYTvETq8Sh2OTg8OiVr2Yth/O3ux0u+8ZzgQyy3vcnI7o3k6Npd\n5rn3U0KkPLufsevFGFJPVYU/9+J4JZTXBLeWNzxwsJ1S7gS3lL28Jhe7H+CJIfdIMPNICsti59GZ\nMzFexOG592Mlh44h9VRV+BnGP0ZCeU1wdtJ6FTv281kwfgZDfJdv7qS1IgRMLErkDSsMizhyXpjC\nXGuPYmclsy7iyJlPFuSoJDeXnvGXaJJTSBaeWq5XsfNIas9it/j0SstOMFerZsmDRSdRQh40mD1w\nlEo9Fh3mIo6Y80z3Y0WOSgjaUhf++EiD+EG93mZjEXAvuWnxqJp1sWMYWImWazUuw3yYoluqVXt0\nP1YYasHB5Gap9MZ0P6kXpuTWAEtyxCzyOVZV4Y/Mjn1eiNXqyTIUD/a8eFTNs9hZyl5jMEYNy/WY\nXCWfbwoMOTi8cmgMDEPjLr7shsVhSWI8akuuVVX422a14pcwTM8JYDFul7/XYpXyHWPB1OjaU+Mo\nmeBWHUoOjpIFyCoWHnLaIpHyxmFBCNgakGPVFn5PTdmDrVlhWBw39YrEMSeXx8YVUzi8Gb8FDk95\ng4nFWDjYcT2kHmbRtsThJU+xEmeOVVv4LVlVyeTymgDacTU4SpiElxZvJfWwjN8DR2lR8pKQhlgu\nsEsaNDg89sFY2YTBwJASdj7l1gAGh1UnmGvVFn7PtslCOyzVcj0Zytj7EmNJLKlXJHrpqCWdEktg\nrIrMYg5ZLvK5hMdSQvJi/CWyF9sdWO3n5Fi1hd9TR/WatB4Yov/UOuoULHcRw9B59Hjmv81y2UIz\n9oJZS7Er6QS9ulGvhdgSB3utPRSCXKu28E8xYUpWey8M0d+ro/Fg/AcP9p/IKsEwhCPGrf1cntVV\nPQ4mblMsEn14u7qfvmJnkfdeHWaJ1JMzn7zkG4+uamb8O+ZZwEpaNw8M0b+P5ZYsQFMtmIyGqp20\nffLG2JKM1yLR1f1E6SuFwROH5YLpkW/MfGJxWCzcGgxDvrlWbeG3LB5eOmrJJpDWN2LuwtF3VM3r\n8409aftkryHG3zbP7m6KWDDdT42Mn10kLEgXk0MsDs+OZt9KPbVLMl5shmEoOczOq6PxSOo+rbov\nFiUT3CoWJcXOKhZ9i4RF0a2xa4zER9v99OHo66AtCEFpDu0bqafvw3a9X9ezHfMooiyz62NVLLMr\nkSH2IuP36sIsc6ikE+wbe0yph91/8tDLh/yZTontoL0IATOfcq3awj90IVdXn72JN4WWWyKbDN2U\nxRTz0mJnuRmV2x4zE2AIRxfj95Q3PGLMYBgae0ypxyvfUhgYQtDla9VB55IjqxjnWrWFn2Eznklt\nVQw85JtSqYdt/7UT0Vv2Yhi/x/Wrof0fGpuRC9kOupTxW7FnhhCMJZF5dgf7Ruph2IxlQo1xuiH6\nl8o3Y0o9Vt3BIoahm7JKux8vyYmNsVXX6CFv9HXQXt2dx0bw0NisXOjVKZVc630l9VgxOzZRPZid\nhXwztpZbyvi1xWDopiy2+9ESgsUXvAB1LJhTkYfSTqmGvR9Pxu8xn6zyLdeqLfwWFxLgVlrP1q2k\nmFtq/LkLm5cUMoTZi+XGO3wZluu1YHYV/tLup1Qu9Op0vaS3obFLZVbPwxIWHXSuVVv42dbNIohj\ntqaslFXK7BZf8JLC7NWaerbp2hhb7I14MP6h91DkyDdts8g3wPckC0PQPD5fTm3x2Lxm4pZr1RZ+\nVsu12LAtKWBWLJdZJJgC5nlUjd3c9WjT+6QeCy23dCNYWziGxrbIi9Kuqsvfi8UzZG4Ih7fU4xGL\nfS31MBeSZRIexdzqRIYnsyvRL71koaGxLeI2diy8CIFXDpXOJy+ZNXXk2aODLpV6LLuffSv1WFxI\ni6NqWt+huwdLE9VC48/Rcj2SupSBMYzfghF3+VvlUKkWzywSXpu7U5AjtrvzWFy7/C3jph0316ot\n/Jatm/Zmr5LNqCEcFhO81JeZAF1jeyX1EA6Lz2ch9ZS09FayV5/vmItElz8rWTCdoBUh8Bi3a8/M\nkkgtdj+RxFpZtYWfZfxjJ7UXjjF9maJklahjSllespdXSz+Ew4u5jrnRzRIpK0JQGosuDF17Zl6x\niEee+x4xnmPVFn4rZsfql6X67Jj6cykjttpgWmy9hxLVs1Mq6SS8NvJZqYdlrqXXurQT9NLAI45S\nQsAy/hIJaYpY5Fq1hd/zdIPHzrsVDotTCBaxKGlNh8Yek/GzLNBT3tCezWdZvAcj7vt8sYvzesuZ\nR8fGxtirtpR2gtYbvNUWfqsLWSL1sPqsl9RjVRjHjIWVVDf1ItHlz8Ri6IUppQvmEOP3WjAPHrzy\nLWcWsldpd2fx+bwWzK6xrWSvXKu28LNabim77PK3OpFRmnwWrGPMlp7FMXZ3x4zb9Yz2nFh4fT4P\nX6sFs2Tz0wqHhdQzRQe9b6UedlWeQsst0RnHZrkWralFm+61MccsbMy4fY93sNr7Ke3uPMb16qBT\nR549JJnSHLI4TNCFw2o+5Vq1hd8qqRmWuzi2pa7tIcmUFjvm86UWTC/Zi1lQvMYtmeAsDgsWPyZ7\nZiULKxzaGuDZQTOsvG/vR+ObGjvHqi38Vq1p6cWZIqm9NEmtb2RkzFE1D8lpTDnNogNjMZfmch85\nGhsDS448JCevPSULQtAlF6aOPDOHD3Ks2sLPMn6Pdszq4nhNLk+5qXSCe7Eqz07JYhOvFLOnlOWF\ngT0AUcqgGcZvIb2Vxq3rzP/Q48j3tdRjwcCYItP1eAf24lhNWi/G77UfwE5EDxxeC5uFvOHJMC1Y\nvJfsZdEVj4m5dO4xGFj/fSP1WCU12x20N/GGxvaSnMacLDkLVd95dC8cXizek/GPtVhZsFELDH1E\nCnj2mX8Wh1fH5iU3MRiGcOx7xj+m1DN0ccZOqLHYM4Mh9ax4LxxjdjSlGFgcFgzT4nRKKSnRYohj\nW+Dw2O+wmHtaDDk49kXhX13lbvn32PyMY5e23mO2kFYaqldrWrpgsr5jLhLxlY5dLLdUcmK6Ki+5\nqRQDO/bYslCphHTwoJ4YsTj2jdSTI0N0+ZZMcHbssSf4mOOm/Kdm8aWdIIOh71nxcdxFuXBsyclC\n6vHAkIPDI5fHlr2Gup9S8pdrVRd+j4DnFDuLSTtWQnnpl304uh5P64nDYlyLjbkuHFZdo0f3M3Yn\nkcP4S6Uei1iUzD1vUrkvGL/nZPGSN7xwWEkspbHoStSuo2opHF5s1ENuYhZBTy3XghB4sPj4HKLF\nTf8pukaPTtBqn6hU6pkZv6PEwibqWDi8xmVj0eVvkdRTfL6Swtg3tlUsvDZsPTDEh9Bpux+LRb6L\n8VuNO5aEZIUj16ou/MwE8JBj+vy7zvyzY3sxcy9ZqG9si6T2kje8ZKG+sdmWvgapxyIWXf5jb+5a\n+VrsB2i7ny4cXY++HsKRa9UWfq8JEAO6uPvOJGpMppIz/1asYyxZqA8H49vnbyELeY6r3fthfC1w\n1OAb/bXdj+e19pJ6tBii5Nm36b9oXTi2t/try74o/F4SS/RfLPxWSc3IEKVn4seWvazkjRJJZopx\nWUKg8e3zH7ujmaITXBw7PstGewy2Vqmnb2yvuJVY1YWfkQpK5RvGd6g7GFO+KWXxfRtzVvKGVyzG\n7KpqkXpKFysrDB6EYHu7+4DAFISg5LRQn79F3Kor/CJyu4jcJyKfFZEf7PH56Z1//7iIvEIzridD\nYTRJC8bvxcBKE7XvWUQW8obXaQhWKvCQm/r8p5B6vOLmKQF6xc3iWlsw/hJyxHaNuVZU+EVkBcDP\nALgdwEsBvFlEXrLg83oALwwh3AbgewHcpRnbKqlLk49p6b0mrdci0efPJB/Lcr0kmVLmGvFqu5/S\nDc2oA2tZbmlH4xXjPhw5e2Zd444Zi9Jx+/yZz8d20LlWyvhfBeD+EMLnQgibAN4D4JsWfN4A4N0A\nEEL4MICrReRUamDP1rS0Tbc4sTC2vGEhZZUyuz5/i1iUMru+ZxGxLNdD9+0b2/t0imbcPhwWc8/z\npI7FuKVhddEXAAAgAElEQVR5YbFnlmulhf8mAA+0fn5w5+9SPjenBvZq8wA+UZnW1IuNlhYDC02y\ntFPqe5bNMsXCi+Wy3Z2XnObBXFmpp29z14IcleQbwHfFJXuNJdbxKyjreEhvpy0cTur+f3fccceX\nv7/55nVsba1f4dO3Knu1Y+zpjVKWW9p6s4la2qbnxLjrWTalC6bViQyPTpBlgefPX/n3Y0o97WcR\ntaWoKeaTRXc3lrQYcZSweM24p0+fxunTp690Iqy08D8E4JbWz7egYfRDPjfv/N0V1i78n/oU8FM/\ndaVPV2Bici4mqnXAc3y7XrvWh6FvbKaADSX1kSM6zFbyhoUs5LWhWcpGS7uDmqWetv/hw8Pj9uGw\nisWY3Y8VefCQTtsxXl9fx/r6+pf/7c4777zyPySsVOr5CIDbRORWETkE4E0A7l7wuRvAdwCAiLwa\nwFMhhIdTA1swV4sJXir1DD3LxqOAMUkNjCtvjL03YiX1lBICC5ZbuvfD7oONSY4spB6vrnjoWUQW\nOVQqe+VaEeMPIWyJyNsA/B6AFQC/GEL4tIi8deff3xVCeJ+IvF5E7gewAeC7NGPnTto2Q+lKEMCv\npV9dBZ55Rj+ux+ZuO1HbkgqTfJ5SDzvBS5mdl+zFEo2STqlvbKtOyaNAs92BBSGw6IoXx20feW6P\n4yX1sLHItVKpByGEewDcs/B371r4+W3suFaJevTolb6lAWflDWaXnpGFusbuS1SrNn3MpLaQN7yK\nnVenZCX1lHRKLOZSQmAhe7HkIScWmvlU+vmGfC9cuPLvc23P3bk7NBG9pJ7S1s0rqdkC5oHZq/3v\nw8F0P9HX42x+qdTDyl6eUo9HLCyI1NDJMOYu3xLyANjIeqVEylrqqbrwM1KP10T0kkK8JjiLubRA\ney6Cpfrs0Nn8mveJLAhBia8FZisi1eerPRnGykJexZwZdwypp9rCzzIUdvX0kCxKMfT5W7A1C8nC\nSyKzYmuL9weMjXlz0wdD30PMLMiDVdc4ltQztG/nLfW0zaorLpW9cq3aws8URsBG6mGS2gNDn/9Q\nonowV7Y78NA6+/yHup8xN/JLF3lm3L6HmDFzpO89wV7dqyd5YBixhdQzdiwYcpRrVRd+dlUeS+qx\nYD5WbM0rUS3kDY9JW5NkUaNv9B9bnvLKoRIMAC8XeuS9lXS6Lwp/zgaMhz479maNhdQzpuzltcll\n8RCzIRwevl4Y2E3xKXCUdI1eGPrGXlnZfdNVLmaWKHp1xblWbeGP7/EsOZFhUcBYFl/SHfSNzU7w\nMaUeKzmtBIPl2KWxsGCB2rh5nk4pJQRec88iFvHIs/aBfGOSh30v9eScyChpm8Zu87ouJPu+zSmY\neUmn5IUh5e8xEWvtMAE/NsrMJ08MtcRizEV+3zB+gF89S5idVTEvuZB979vMSWqPAmbB1jwwDOEo\nLaRee0peRXRo7DGlHnZh88AQpRztSSvPPUEPolhiVRd+i9a09EKWTloLDCsrjeTV9RJnppCWToDS\nYr652Y2h9Drn4PCQWTwxMO1/DTjG7g4YDH042HntUczZjibXqi/8paun14kFCwzacYf8mUXFY0Pa\na8+FlXqmwOExwUuLaN87D/r891oxt+jA+nCw9cKCKM6bux1mIbNY7Acwxc5DQ+3DPPbJAouktmCi\nzNMSS+M2dktv1aEsyoWAXyEds5gz+TZEpMY8aWVBNPaV1DMmAxtbI2YwsDjGlrLG3FQV2X0CqRaz\nx8Zq6SJo0WHmyBseG5qeGrjH4uqNg1m4S+JWYlUXfrZ4TJ1QpbJJSuopSaixk9pCbtJ2P9vbzYKw\neOZ/CIcFwxyLaFjJG2PisGDPXgStD/MURLH0UEOuVV/4PVZPzxMLpYWRkXqmYB2lLLdU3ujyZxdM\nq7iNxRj74hYXunb3wxY7T8wee0psN2rRQY/Z6TLkqMR6wlKHeV4cr80aD7xd/vHMv5bleskbm5vA\n2pretxRDlz+7YHptwnqREk0s4mZuzoamR95b7WF4kJ0cHB4kpnTcj34U+PCHr/TVWPWMv0TqsUoo\nZtJ6XPSusaNv1yZe6US0YvEl3QFTzC2ZXWkOWbT/JZ0gs2Ba4RibaHRh6HoMg1Usxt4T1GL4y78E\nfv/3r/TVWPWFvz0R+15cHn292NpYG3NWE5yVp0qlHqaIemxSsrGwYrlTS0hd/kzcrHCwi7zH3Ot6\nDAMzn0Lwi4VXvg19vpRVXfgXAxPPJ3uxXI/N3SGGsngskZF6rDbxvBYJL4msa2w2FmMv3BaMWFvA\nrPaJSnNo7HzrGptZBKNsanFAoJQQWEjDKau68LPFriShLNp0BkP73bgpDF1jWzH+0oVtbEbchYMp\njENjW0g9FjFmYrGYn6y8YSEBeu6NsHJobixyFsxSwuPlq7HqCz9T7DyYuVcb2+Vv5evJzEsYSirG\nTPfTxtGHoQ8zSwhKDwh4Yegam2HE7LNsvBY2iw6sy5/tDmogUkPziVkEU1Z94c9lxNGfCWIpGy2d\ntExSsyzXIhZenVLXk1hZZucpe2k/H5MXFqSky5+JBctya5ALLYu5B5EqrS0W+aaxqgt/V6JaMH7P\nXfqSNt2qpS9drDw3dz2KEhOLvkdfd407hMNrcfWUN0riNuTvtZ/DYOjCYRm3sTZ3Wca/LzR+b5br\nJQtZMX5tS++1oend/eROxBxGrH2WjUUsmBzqejeul7yRk0MWc6Q0h7TXeopYlEqn7IK5lIyfvTge\nRYllPmNtUjLFIHUMtnQTz0qyyJ2IVrJQaTFnGbGHbs8yYouu0aI7KGXPXTgsY+El9WgZ/77a3C29\nOF7tmEUb2zW25UZwO/mGjsFa6M9jt+ldsbDYCGaKksWmuIUMUbKwsRo/Ewur7oc58jzWfBrC4Vlb\nmBxKWdWFn2HEi0Gc4oaMrkS1ZCiLvtpYWCa112ZUiexlyfg9OhoL2cuqgOUuKENjs3HT5lvXkWev\nvR+rWHiqCUwOpaz6wq9la4u+7A0ZXonqxXJrSurSuEV/rwme2x0M4fAiD11je8kbOYTAQ+phP5+F\nryc50s4RJhYrK43vYvezlIy/tP1nLqRFa1qKw2rckpYX4Fnu0ATQdj8l8oZV4SiVeiwWia6xPWPh\nJXtNEYuSuceMO/SAxNLa0oeDPfKcsqoLPyv1jJXUXonKSj0W43pJPTFRmdMpDLPLlQBZxs8sgl6y\nl5d8wxICr41uL8JjtUj0EdAxT4Z1jb1vNP7UpB0zqS0m7VibUQyGIX+G8ZfiYD+fx37AEA6LBdNK\n4y9ZMLVxszoZxiwSXWN7FnMPDEM4cvZ+tDhSVn3hHyOph8ZmV+USHFbFoCSph/wZjb8UB9tVWTF+\nD+Z64EBTONvdTw2yF4NhaM+MJVIlLNdyERyDlEQcHox/aaUeZtIuSiEsy7WUenJxWElZJcxnCEdO\n9+PV0ZTIacyktZA3RMoWQS+pxwqD53wqweEVC7ZrZCRAq1xOWdWFn520JSzXa3PXqkCXSD1WrWlf\n3MZm8Wxh1MpCXSwX6H+8g7ZrZDF7LhJeGLQ5xB559tzcHYNoDI2dU1u01y9l1Rf+3JWW1fiZdsyr\n9WZXew8MQ2OXSj1WsSiReqykAnaT0othepIjjxwS2T2a6I3DS1r0lHqsupSUVV34S05kMAGM/kw7\nNnYB85oAixhC2L3TN4Uh4phicjF5wciFHuShD4fFtR5rP8BK3rDA4bX344FhaOx4nUtu+Fxaxj9G\nYRzyL5V6LFtvr+5ncdyht5yx8kZJO20xAUrYs1UO9eFgCrSH1GOJoUSy8CQEU2O4fLkp7F2b4uyR\nZwZHyqov/Iw+W8Jy+/z77pizaMdKJpeVLFTCiOPYHvKGFTP36hpZxr/XJUAmbrXg8OzAWAxdRKrL\n3wpHyqov/F6Mf3HcPpbbdcecFYsvlSysOolSZudxIsNyESzplKwWwVIczOLK5L2HbDIkF7I4vKQ6\nL1LCzFNPHCmruvCXsGfLALJje7HcEq3TKqkZqWcsHDVgiP5ei/zYhKAEwxCRYnF4LWxeGNgN2FIc\nS1n4a2C5XThSF34Mluu1OTiEoeslIVaYS2MxRUvfJRdaFTAPeWPMhdiSSO1lQjCEwQLHvtD4vVZ7\ny1XZi/GzshejdTLdjyeO3JbeCkOJ1LO93X9naxcOVrLw0ss9ZK9UQaoBx5ikJBULD2KSsqoLv9dF\nZxO1RHJiFwltJ+GFQVP4PbTqElZlda0ZDIs3IuUwO6tYTC1DMuPWgmNMUsLOJwtCkLKqC38N7LLP\nf2yGORYGpoDFo2pDm3heMoSHL3M9Fm9EYoudFeP3XIg9MIyJYwpCwDJ+Lxwpq77wlxRzqwDWgKO0\ngFnGol3sLNtYj802VqbLLWCsvMHIll5SjzUpid1PTiw8yNEUEmCpjDwzfvDyjVcA2wmVOqpWKgtZ\nsBkvlruIw1siG7v1LukEUxPci8XXQAgWb0TKiYUXIfDqoC1IySKOoRe89OFYysLPruDaC7moz6Ym\nePtibm3ZHVVjJpcny/Uqdl442MLosZG4OLZmgpdIMlMTAubzaRh/CQ6LYj4mKdGSo4hhqLYwYw9Z\nduEXkWtF5F4R+YyIvF9Eru7x+5yIfEJEPioi/5f5HV4a8eK7cT2ZnRXD9JwsjNSzmKhWjL9kEbRk\ngbmxYIsdW8ynIAQlEqB2wWRv9rKSWbxISex8tN1PSQc9ldTzQwDuDSG8CMAf7PzcZQHAegjhFSGE\nVzG/wEsjXhzbUt7wYpilkzaV1NruJ3fBTOEo1XI9WKA3IRg7Fp6yVy4hiL5MB22xCHoRgsV3L+Qw\n/j5jF+MhKyn8bwDw7p3v3w3gnw/49lzWYStJVCaIbKIyxc6L5VrFYvGRFN7yhkUsSqQey2KXSwhS\nN3uVdgdTyV5M9+PBclmiOBYhYAjoXmD8p0IID+98/zCAUz1+AcDvi8hHROR7mF9QenG0QbS8kKUy\nhFeiWi1suWwm5c8Uu5IYW8obuXGzvNmrFglwkRB47hN5SD1TxiInbhocQza4NSAi9wL4io5/+s/t\nH0IIQURChx8AfF0I4Ysi8hwA94rIfSGED3Q53nHHHV/+fn19HQcPrrswu0V/T6nHU+O3jEWOVs1O\n2lRROn/+2b7aRYKd4IcOdfu2H0lx4IAvs2NkIUuN34vl5hICS/LgSQi0MV4c2zIWMZdPnz6N06dP\n45lngHe+s99/yAYLfwjhdX3/JiIPi8hXhBC+JCLPBfBIzxhf3Pn6qIj8bwCvApAs/ADwJ3+Sz2Ys\nNf4SqceKYbJFlI1FziLIFLsoJVmx3JKNuaNHhzFvbgKHD3NSFhNja4ksN9+sZS9mEfRguaWEgJlP\nHh2NtmtcX1/H+vo6fuIngB/9UeDHf/zO/v/UYyVSz90AvnPn++8E8JuLDiJyVERO7Hx/DMA/AfBJ\n7S/wSmogX+phGYpVQpUyFKtY5PoyGFJjl8TYKy+8xo3+FnnP5BD7SIqcApbCsDhu9Lfa+/GSAHPJ\nkWUsUlZS+H8cwOtE5DMA/tHOzxCRG0XkvTs+XwHgAyLyMQAfBvA7IYT3a38Bk9RdZ/O9WK7Wl70h\nw4vNWMcih9lZT/Bc2curE7Rk/GMSgj5f9sgzG4tcQmC1sHkRNCCfKDKxiMdgh8Yessz/BoQQngDw\n2o6//wKAf7rz/V8B+Nu5v4O5kDFRYzBqYPwRQ99RNSahmALGslxmgjPFvFQimyoWXotgTvufwlFK\nCDQ4Vlf95A3NItHe+5mCHLXfwidiLw0zsVisQ321JWVLc+du9M9hmOyFtMTgIW9EDNrux0u+KZHI\nLCc4s9/BEgIPeaOEEAzlZ1cOecXCi+VadXfMfCo58mxNjrQ1K2XVF35Gt2dZhxdby2V2TAEbwnzg\nQMMEcu8e9JB6LFkuK2+wOZRDCFipp4ZYsJ2uRyxYWSj1uBRGAswlBF51yJKUpKz6wl/C1jzYjBeG\nlD+rVefisOxoPFmuh+7bhdlqgudq4JHlesSCvdZepMtjIU6NXUoIPOYTgyE1bsqqLvwl+qxXQnlh\nSOGoJRa5k5YpBtGfYfFD3Q/77JTcz2fVNS7mm+VDAbWL6yKOGsgRgyE19liEYKpYpKzqwt91Ifse\n5gTsPcbPTERPxs+29CWbUUO+cdwQ0ougNhZdp1MsO5rcPSWGBWoxpPyZBbMLs+Ui6EVKSmJhJXt5\nxmJfFP5Dh66c4IcP9/t7BbGNY3Oz/87PxXEZZpfC0VXsPHCwLHAIQ+4mXuqRBl1x88BRQyxYeWMI\nx8GDZY8jHypghw7t+l66pI8Fq/Gn5mk7FkM4WMbPLEDtenHpUrpm5ZKSpS78ly7t/pxKKK8gtnFo\nMOQwuxTLzYlFLhudOhYMhlpwWGNgSIk2FvF1kfF0CiNPMTimIiVMLJi9kRIczCJoiSFlS1X4GYZS\n24VMsdzSWCxTUnvGIrfYMaTEUt4oXQQ9cHguxKmugyn8m5t5N3xqOo92LKxk5H3D+FdXm8DFi7Ps\nST0ly2WK3cWLegwe+wyeneAYi6AXBhbHXiNHlrEoeV0kswB5kZLUtUtZ1YX/wIFnf1ivpNbol17M\nLgeDBw7vRbAWljvHYn+Qo+3tXWlraGxvYjIVKUlZ1YUfKAviFEnNMjsthjhuZCheOCwXwTEweOLY\nC3nRxhBCWlqoIRbMPPXCEHEw86+9h+G1CKZqlhZDyvZc4dcW89Ru+uHD/hcyhWFxQ2wIQ3ylW3ts\nLY6LF+1w5MbCC4MnDq9O0CsWEW/fPtEiDsv8zF0ENRi8Cv8YOFIb3QyGw4f3icYP5CfUxYvDvm2t\n2iupUxiYxWfRn51cVjhyi53lQryy0nQ+29vNV+Z0yl4gBG0MlosPi8NzEWwvxFPM01pweGFIWcH2\nwDhWA9v2KnbshawBRw0LscguA4vfDz2lcK8RAiY32wcgWJZrSUxKGL8Wg/UiuMzkKGVLxfgX22mv\nC6mVm6wZf8QRz2GnNq48cJQsxB6x8GS5U0qA2sVnZWX3yZEesfBiuVOTkhIcXvtgDClZ6lM9wG4Q\ntRtXUye1F4Y2Di2bmWNxJY4aFkEvQmAtb+y1WMTrrO1+WJllanLEYEjZnin88QanIZbLBJzRy5lx\nczGkEq89tiap2ziYpGY2HplYWGKI/hcv8rHQYI4LkGUs2HzTYmjj0MRiEYd2MbbOi/a9IFYY2gcg\ntHlx8eLuXtEQg16cqx5zhJF6NHkxZHum8GuTup0kU01w7bgxSS9f5ie4ZiNYO7mYBSg3FtZJnZsX\nlji88mKMWITQfB3Ko7FiYU0I2FhE36F9orHqhSUhGLKlLfw1TPCU7yJD8Sp2tSyCU7LcmmIxNeO/\neLHRiFdW0h30GLGYehHUzL2x6sXM+HfMc4JPzWai/6VLy89yUxhim729bS97RRzx6OfULFcjb3gX\nO2ZcDY69TAi0+bbXyNGQLW3h99IvUzdkMBjaOLxbekvNtZ3UqSLKLILesdhLLf2lS7vXzjMWKd+p\nu+KZ8edhSNlSF36rhNrLjF/zzBImFm1fy6QGdj/flBN8keVOsd8Rn1G1uTl9LKaWembGn4chZUtZ\n+DXHuWpI6jaOqVv6qZM6+sfPx8RCu9HNYIiYtUVprxECj1ho8+LIEX0sDh7cvUt7yu6n3YWxXfFU\nEuCQLWXhj2dch1r6GtrY9thMUqfkpjYOBgMwfSy8pZ4aFkGWEGiv3zLGIt6ZPXUnGOvF5ubuTXMp\nDIAuFhcuNN9bS4BDtpSFf0qWG1+MvbWlw+HN7GrRL2uIBctyWWZnud9RQyziuNvbzR+ru1Vz9n5q\niAVLSgBbchQXHG1tGbKlLPxTstwcHN5Sz1T65aFDu/cp1BALBoN3S19DLNiFOLUp7tX91ESOvAu/\nJY4hW8rCz7BLwK/w18ByPRl/SnJqt+l7ZYK35TTrlj4uKOzn03ave5Hlsjg8ZK+9Gou58O9YDrMD\nZsYfMWhuXS9ZBPdCLFi5kJEAY5uuxZyzCGo3uhkMzIIJLD/jZ9QE61gwsteQ7dvCHy+k9gSQV+Gv\nIakXk8mqpS/FkRo3MvOpMHjHwnsR9MAApHGsru6SjL1CCFgMzFNp2b2fmfEvWE5SR5Zk1dLn4lhG\nlpuDw3MR9MIA1IGjBlICpDuPKAF6dh4enTkbixDST9FsS4AsMUlhHrKlLPzWF5J92QWT1DkJVUNS\nA3uP5XphAOrAMTWGtmSxl2LBdJhMcdYcK4/PS4o361kuxkO2lIWfuZAaX/ZlF144aojFXma5noyf\nkZz2SixyCIHm6OcYOFJdB7B7IxlD5qyvB5C3J6hZgIZsqQo/y7QZrcyTbXszO0sM7MsuaovFzPj1\nGHIWn7gADrHcMXDUEAvN4tPGYU3ShmxPFH5Wq7a+kLXgqAFDfJS0VmfM2T/wioU1u4zjxlfnDT0P\nyRMHq5czGFipR6s9e+LIyc0aYuFBWPus+sK/tgacP9/8WVsb9s29kJqLc+SIP47UhTxypLm9e0oM\nEcfGRlPwLBk/s0jkxMJ6Q/PgwWYhPHs2jcETx9paXiwsCQGDwRsHO0+nxMDiYAjBkC1l4fdg/Gtr\nTbHb3NQzQetdeq9YsMm0tgY8+WTz1bKlZ3B4xoJZBNfWgCee4CY4yxo9YmFNCFZXm+OZ7CJojaOG\nWDAkMeI4f57fG1n6wu/BcqNcceGCfoI/9VRzUbXFzloK8UrqnEUwp9hZs9y9GgsPxp+zCFpiEHk2\nIUgZK29MHQtWAtzcbIiiNhZnz6ZPALE4hmxPFH72Ql640BToITtwoAni00/7TfALF9L+UbLwSmot\nhvPn9bKXZywijhpi4ZkXqfz0joXG1zsWZ882DDe1NxI/nzXb9oqFSOOvXQSPHGlIpdaXuSZ9tpSF\n/9w54Nix9NjHjjUB95rg584BR48O+x492vh5TXANhpWVhm2cOeMbi5R/OxYestfGRjoWa2v66xH9\n2WJ34ED6tEcNecFgAJqixMTiiSfSGFgcXvXi2DHfvHjySX3Nmgv/gjETHGh8nnzSr9hpWr1jxxo/\njf7sGYtjx/xioS12zOTyKnYHD+4ugh7dj0exy2G5mmLHFn6vYrcXF8GcWGjzYmNjLvydvkzhf+yx\ndNsdfbUXcm2tKRwMs2M+H+OrSeqI4/HH9Z+PicVjj+kn+MaGTpLJjYWWVWljwbDctbVmXC2Gc+d0\nshATi7hIaDvBQ4d297ZSxuTFkSNNLLSkJM6RqWLRnqdTxuLo0Uae3t7WEZM+W6rCH28ceeopfRAf\neQQ4cUKHQ3shjx8HHn2UY3bPPNP8vxQGbSxiJ8Esgo88ksYQcXjHInVNmFgcP97EgVkEa8iLjY06\nYvHww9PH4qmnmmJnWfiZWMRF8NFHfWLxyCP6RfCRR5qvqY3gIVuqwg/wQXz4YX2x07LA48ebcbUY\n4gS3LPyrq41soW0hY0Jpk5qJBXM9YixSrLimRdAjFrHYAbb7HSdONPFlC78HIdDOkfb1SBU7tvA/\n80yTF9qO1IscxWJuiWHIqi/88eJoCmP0ZyaXNojHjwNf+hKHgZE3NMmXGwstDu0Ez4mF9no88cSu\nzp4aVxuLlZVGD2VkFq9YPPywD4aNjWYvJeV/9GhTFDWLa8TBzhENeThxQp+bLIaYFykcOYsgg+OL\nX9TFwqtmDdmeKPwbGw1zveoqnb/H6nnVVcADD3AYGG398OH0zRsnTzZ7B2fONN9rcOSwqpSdPNnE\ngsFgfT1Onmy0Tq9YaIsuGwsmLxgMZ8408UjhOHCgYaGMzMLGwmOOaDFcdVXTKWklsosXmwXTWiFg\n6kVcBK0xDFn1hf/AgeZDPvigfnJ98Yu6QhN9tUn9+c9Pi2Ftrbkx5LHH9JPrS1+yx8Eugl4Y4iI4\nJQ628HvF4vHHmyKmHdszL6xjceyYvqs6eRJ46KGmiKbuDxDhCimLg4mFNsYMhiHLLvwi8q0i8uci\nsi0irxzwu11E7hORz4rID+b8rquuatoxbRCfeAK49tq077XXNp2Eph3LwXDddXoMmgspsotDW+w8\ncJw82XRhU8bi0KHmWUGbm7pTFiyOp57SF7uNDZ/rwWA4d06ngXviiHlRQ25qMLRxMPXCIxZPPmmP\nYchKGP8nAXwzgD/ucxCRFQA/A+B2AC8F8GYReQn7i+KH1LCDWIw0CRV9NBcnjsv4ai8k8OzF5/Tp\n073+UQ5KyUIsDiYW0UdT+EsxDMWCsRwczOfzikV73L5YRGY79Ba5tsX4MvnpFQtmnmpiEWuEplYA\nu5uvmrz3yov4u3NjkWPZhT+EcF8I4TMJt1cBuD+E8LkQwiaA9wD4ptzfqUnsG29svjJBjP9nyJhi\nx2CIJzbaiTpU7EJIjxntuc/V42BiwSyCDIY43vXX7/6dVeF/znOar5rrF4tdxD5kXrHowpCKhTY3\nIsnQ3ADknRfMItjG0BeL2PFoYxG7RU2nlBML63rBYBgyb43/JgAPtH5+cOfvKLt8We8bi8bVV6d9\nmSA+//nN11tvTfvecEPzVdtuAroXNgDNWWatxQmuOVIWJ6B1LOJ411yT9o2Tr+TGlD6LRU5DHrzy\n4qadzNcUu4hXs+EXTTtPtLkG7M6jU6fSvjEGMSZDxhQ7plOK5hGLeC00uRxjwcwRZhHUXI8hG5wG\nInKviHyy488/U45P8NN+0yRStFOnmsmd2tgBgK//+uarhtnddlvz9SUKoSrKMP/wH6Z9AeCVr9SN\nG3Fo9iSAhuVqtcDXvGb3/6QsYn3Zy9K+sZPRxuK224AXvEDn+8IXAl/5lTrfG2/Us6Rv+IbmqybO\nMQYvfWnaN07a9XUdjhtuAJ73PJ3vddfp43brrcBXf7XON143TYGMsXjRi9K+cc5p8iKeRNLMU6Ah\nEJqCCzT59jVfo/ONc0TTHbz4xc1XDY6bb26+ft3XpX1jh6KNRZ9JYLSDrgFE/hDAvw8h/FnHv70a\nwB0hhNt3fv5hAJdDCD/R4WuySMw222yz7TcLIVD38Sq2CFXW90s/AuA2EbkVwBcAvAnAm7scWeCz\nzZzpHdAAAAUhSURBVDbbbLPlWclxzm8WkQcAvBrAe0Xknp2/v1FE3gsAIYQtAG8D8HsA/gLA/wwh\nfLoc9myzzTbbbLlWLPXMNttss822t2zUO3c1N3OJyE/v/PvHReQVY+Ib01KxEJFv34nBJ0TkT0Tk\nb02BcwzT3uQnIn9PRLZE5FvGxDemKefIuoh8VEQ+JSKnR4Y4minmyPUi8rsi8rGdWPzbCWC6m4j8\nkog8LCKfHPDh6mYIYZQ/AFYA3A/gVgCrAD4G4CULPq8H8L6d778GwIfGwjfmH2Us/j6Akzvf376f\nY9Hy+z8AfgfAv5ga94R5cTWAPwdw887P10+Ne8JY3AHgv8c4AHgcwMGpsTvE4usBvALAJ3v+na6b\nYzJ+zc1cbwDwbgAIIXwYwNUiUnhitUpLxiKE8KchhDM7P34YwM0jYxzLtDf5fT+AXwfw6JjgRjZN\nLP4VgP8VQngQAEIIj42McSzTxOKLAOLdMlcBeDw0+4pLZSGEDwB4csCFrptjFn7NzVxdPstY8Ngb\n294C4H2uiKazZCxE5CY0k/6unb9a1o0pTV7cBuBaEflDEfmIiPyb0dCNa5pY/AKAl4nIFwB8HMAP\njIStNqPrptVxTo1pJ+visc5lnOTqzyQirwHw3QAUt3fsSdPE4qcA/FAIIYiIoP/48F43TSxWAbwS\nwD8GcBTAn4rIh0IIn3VFNr5pYvGfAHwshLAuIi8AcK+IvDyEcNYZW41G1c0xC/9DAG5p/XwLmpVp\nyOfmnb9bNtPEAjsbur8A4PYQwlCrt5dNE4u/A+A9Tc3H9QC+UUQ2Qwh3jwNxNNPE4gEAj4UQzgM4\nLyJ/DODlAJat8Gti8bUA/isAhBD+UkT+GsCL0dw/tJ+MrptjSj1fvplLRA6huZlrceLeDeA7gC/f\n9ftUCOHhETGOZclYiMjzAPwGgH8dQrh/AoxjWTIWIYSvCiE8P4TwfDQ6//ctYdEHdHPktwD8AxFZ\nEZGjaDbz/mJknGOYJhb3AXgtAOxo2i8G8FejoqzD6Lo5GuMPIWyJSLyZawXAL4YQPi0ib93593eF\nEN4nIq8XkfsBbAD4rrHwjWmaWAD4LwCuAXDXDtPdDCG8airMXqaMxb4w5Ry5T0R+F8AnAFwG8Ash\nhKUr/Mq8+G8AfllEPo6GxP7HEMITk4F2MhH5VQDfAOD6nZtm34FG8suum/MNXLPNNtts+8yqf/Xi\nbLPNNttstjYX/tlmm222fWZz4Z9tttlm22c2F/7ZZptttn1mc+GfbbbZZnMwzcPViLFes/Ngvvjn\nvIi8IXu8+VTPbLPNNpu9icjXA3gGwK+EEP6m4bjXoHmA3U0hhAs5Y8yMf7bZZpvNwboeriYiLxCR\ne3aes/THIvLijKG/Fc3TOLOKPjAX/tlmm222Me1/APj+EMLfBfAfAPxcxhjfBuBXS0CM+aye2Wab\nbbZ9ayJyHM17Nn5t5258ADi082/fAuDOjv/2YAjhG1tjPBfA30BzR3O2zYV/ttlmm20cO4DmOTpX\nvCErhPAbaJ7NlbJ/CeA3QgjbpUBmm2222WZzthDC0wD+WkTeCADSGPtK1TejUOYB5sI/22yzzeZi\nOw9X+yCAF4vIAyLyXQC+HcBbRORjAD6F5u1Z2vFuRXOS54+Ksc3HOWebbbbZ9pfNjH+22WabbZ/Z\nXPhnm2222faZzYV/ttlmm22f2Vz4Z5ttttn2mc2Ff7bZZpttn9lc+GebbbbZ9pnNhX+22WabbZ/Z\nXPhnm2222faZ/X8/As7hgGvCPQAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "#Carrier wave 300MHz, sampling freq 10GHz\n", "n = np.linspace(0, .00001, 100000, dtype=float)\n", "c = sin((2*pi*300000000)*n)\n", "\n", "plot(n[:1000], c[:1000])" ] }, { "cell_type": "code", "execution_count": 157, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "(0, 1.2)" ] }, "execution_count": 157, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYsAAAEACAYAAABCl1qQAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAEdNJREFUeJzt3X+MZWV9x/H3113QGH8g2QTr7tJtZWOgdq2YLgQbuUZS\nRhKhaRPpVhuqojRmTUzTiPiHDFFbf8caWrpSQJtWt001FlNW1OgN1vKjWASkLO6KxN3F0lKVGIp1\nt3z7xz3sXmdn99k5M7P3OzPvVzKZ8+M593zPw3A/9znnuRCZiSRJR/O0SRcgSarPsJAkNRkWkqQm\nw0KS1GRYSJKaDAtJUlOvsIiI6yPikYi49wj7XxcRd0fEPRHxjYjYNL8yJUmT1HdkcQMwdZT9DwKv\nyMxNwHuAT/Q8jySpgF5hkZlfB350lP23ZuZj3ertwLo+55Ek1XA8nlm8CbjpOJxHkrRIVi/mi0fE\nK4E3Ai9fzPNIkhbXooVF91D7WmAqMw+7ZRUR/kepJKmHzIzjfc5FuQ0VEacCnwNen5m7j9QuM/3J\n5Morr5x4DVV+7Av7wr44+s+k9BpZRMRngHOBNRGxB7gSOAEgM7cB7waeB1wTEQD7M3PzglQsSTru\neoVFZm5p7L8UuLRXRZKkcvwGdwGDwWDSJZRhXxxiXxxiX0xeTOoeWETkJO+/SdJSFBHkcnnALUla\nXgwLSVKTYSFJajIsJElNhoUkqcmwkCQ1GRaSpCbDQpLUZFhIkpoMC0lSk2EhSWoyLCRJTYaFJKnJ\nsJAkNRkWkqQmw0KS1GRYSJKaDAtJUpNhIUlqMiwkSU2GhSSpybCQJDUZFpKkJsNCktRkWEiSmgwL\nSVLTnMMiIq6PiEci4t6jtPl4ROyKiLsj4qXzK1GSNGl9RhY3AFNH2hkRFwCnZeZG4C3ANT1rkyQV\nMeewyMyvAz86SpMLgU91bW8HToqIU/qVJ0mqYDGeWawF9oyt7wXWLcJ5JEnHyepFet2YsZ6LdB4t\nI08+CZdcAj/+8WTriID3vhc2bZpsHVIlixEW+4D1Y+vrum2HmZ6ePrg8GAwYDAaLUI6WigMH4NOf\nhs9/frJ1fOhDcN99hoVqGA6HDIfDSZdBZM79Q39EbAC+kJm/Osu+C4CtmXlBRJwNfCwzz56lXfY5\nt5avn/0MnvWs0e9J2rIFLrxw9FuqJiLIzJl3bxbdnEcWEfEZ4FxgTUTsAa4ETgDIzG2ZeVNEXBAR\nu4HHgTcsZMGSpONvzmGRmc3PW5m5tV85kqSK/Aa3JKnJsJAkNRkWKqPSfIdKtUgVGBYqJY77HI+a\nNUjVGBaSpCbDQpLUZFhIkpoMC0lSk2GhMirNQKpUi1SBYSFJajIsVEqFaasVapCqMSwkSU2GhSSp\nybCQJDUZFiqj0gykSrVIFRgWkqQmw0KlVJiJVKEGqRrDQpLUZFhIkpoMC0lSk2EhSWoyLFRGpemq\nlWqRKjAsVEqFmUgVapCqMSwkSU2GhSSpybCQJDUZFpKkJsNCZVSagVSpFqkCw0KS1NQrLCJiKiJ2\nRsSuiLh8lv1rIuKLEfGtiPh2RPzBvCvVilBh2mqFGqRq5hwWEbEKuBqYAs4AtkTE6TOabQXuysxf\nAwbARyJi9TxrlSRNSJ+RxWZgd2Y+lJn7ge3ARTPa/AB4Trf8HOC/M/NA/zIlSZPU59P+WmDP2Ppe\n4KwZba4FvhoRDwPPBl7brzxJUgV9wuJY5om8C/hWZg4i4oXAlyPiJZn5k/FG09PTB5cHgwGDwaBH\nOZK0fA2HQ4bD4aTL6BUW+4D1Y+vrGY0uxp0DvA8gM78bEd8DXgTcOd5oPCykStNVK9WilW3mB+mr\nrrpqInX0eWZxJ7AxIjZExInAxcCNM9rsBM4DiIhTGAXFg/MpVCtDhZlIFWqQqpnzyCIzD0TEVuBm\nYBVwXWbeHxGXdfu3AX8C3BARdzMKpHdk5g8XsG5J0nHUazprZu4AdszYtm1s+VHgNfMrTZJUhd/g\nliQ1GRaSpCbDQmVUmoFUqRapAsNCktRkWKiUCtNWK9QgVWNYSJKaDAtJUpNhIUlqMixURqUZSJVq\nkSowLCRJTYaFSqkwE6lCDVI1hoUkqcmwkCQ1GRaSpCbDQpLUZFiojErTVSvVIlVgWKiUCjORKtQg\nVWNYSJKaDAtJUpNhIUlqMiwkSU2GhcqoNAOpUi1SBYaFJKnJsFApFaatVqhBqsawkCQ1GRaSpCbD\nQpLUZFhIkpp6hUVETEXEzojYFRGXH6HNICLuiohvR8RwXlVqRag0XbVSLVIFq+d6QESsAq4GzgP2\nAf8aETdm5v1jbU4C/hw4PzP3RsSahSpYy1uFmUgVapCq6TOy2AzszsyHMnM/sB24aEab3wM+m5l7\nATLz0fmVKUmapD5hsRbYM7a+t9s2biNwckR8LSLujIjf71ugJGny5nwbCjiWu7knAGcCrwKeCdwa\nEbdl5q7xRtPT0weXB4MBg8GgRzmStHwNh0OGw+Gky+gVFvuA9WPr6xmNLsbtAR7NzCeAJyLiFuAl\nwBHDQpJ0uJkfpK+66qqJ1NHnNtSdwMaI2BARJwIXAzfOaPOPwG9ExKqIeCZwFvDv8ytVy12lGUiV\napEqmPPIIjMPRMRW4GZgFXBdZt4fEZd1+7dl5s6I+CJwD/AkcG1mGhaStET1uQ1FZu4AdszYtm3G\n+oeBD/cvTStRhWmrFWqQqvEb3JKkJsNCktRkWEiSmgwLlVFpBlKlWqQKDAtJUpNhoVIqzESqUINU\njWEhSWoyLCRJTYaFJKnJsJAkNRkWKqPSdNVKtUgVGBYqpcJMpAo1SNUYFpKkJsNCktRkWEiSmgwL\nSVKTYaEyKs1AqlSLVIFhIUlqMixUSoVpqxVqkKoxLCRJTYaFJKnJsJAkNRkWkqQmw0JlVJquWqkW\nqQLDQqVUmIlUoQapGsNCktRkWEiSmgwLSVJTr7CIiKmI2BkRuyLi8qO0+/WIOBARv92/REnSpM05\nLCJiFXA1MAWcAWyJiNOP0O4DwBcBHxmqqdIMpEq1SBX0GVlsBnZn5kOZuR/YDlw0S7u3Af8A/Nc8\n6pMkFdAnLNYCe8bW93bbDoqItYwC5Jpuk5/TdEwqTFutUINUzeoexxzLG//HgHdmZkZEcITbUNPT\n0weXB4MBg8GgRzmStHwNh0OGw+GkyyByjjdnI+JsYDozp7r1K4AnM/MDY20e5FBArAH+B3hzZt44\n1ibnem4tb488Aps2jX5P0lvfCi9+8ei3VE1EkJnHffzbZ2RxJ7AxIjYADwMXA1vGG2TmLz+1HBE3\nAF8YDwpJ0tIy57DIzAMRsRW4GVgFXJeZ90fEZd3+bQtco1aISgPNSrVIFfQZWZCZO4AdM7bNGhKZ\n+YY+55Ak1eE3uFVKhZlIFWqQqjEsJElNhoUkqcmwkCQ1GRaSpCbDQmVUmq5aqRapAsNCpVSYiVSh\nBqkaw0KS1GRYSJKaDAtJUpNhIUlqMixURqUZSJVqkSowLCRJTYaFSqkwbbVCDVI1hoUkqcmwkCQ1\nGRaSpCbDQpLUZFiojErTVSvVIlVgWKiUCjORKtQgVWNYSJKaDAtJUpNhIUlqMiwkSU2GhcqoNAOp\nUi1SBYaFJKnJsFApFaatVqhBqsawkCQ19QqLiJiKiJ0RsSsiLp9l/+si4u6IuCcivhERm+ZfqiRp\nUuYcFhGxCrgamALOALZExOkzmj0IvCIzNwHvAT4x30IlSZPTZ2SxGdidmQ9l5n5gO3DReIPMvDUz\nH+tWbwfWza9MrQSVZiBVqkWqoE9YrAX2jK3v7bYdyZuAm3qcR5JUxOoexxzzZ66IeCXwRuDlPc6j\nFajCTKQKNUjV9AmLfcD6sfX1jEYXP6d7qH0tMJWZP5rthaanpw8uDwYDBoNBj3IkafkaDocMh8NJ\nl0HkHG/ORsRq4AHgVcDDwB3Alsy8f6zNqcBXgddn5m1HeJ2c67m1vO3ZA+ecM/o9SW9/O2zYMPot\nVRMRZOZxH//OeWSRmQciYitwM7AKuC4z74+Iy7r924B3A88DronRmH5/Zm5euLIlScdTn9tQZOYO\nYMeMbdvGli8FLp1faZKkKvwGt8qodFeyUi1SBYaFSqkwE6lCDVI1hoUkqcmwkCQ1GRaSpCbDQpLU\nZFiojEozkCrVIlVgWEiSmgwLlVJh2mqFGqRqDAtJUpNhIUlqMiwkSU2GhSSpybBQGZWmq1aqRarA\nsFApFWYiVahBqsawkCQ1GRaSpCbDQpLUZFhIkpoMC5VRaQZSpVqkCgwLSVKTYaFSKkxbrVCDVI1h\nIUlqMiwkSU2GhSSpybBQGZVmIFWqRarAsJAkNRkWKqXCTKQKNUjV9AqLiJiKiJ0RsSsiLj9Cm493\n+++OiJfOr0xJ0iTNOSwiYhVwNTAFnAFsiYjTZ7S5ADgtMzcCbwGuWYBal63hcDjpEsp44onhpEso\nw7+LQ+yLyeszstgM7M7MhzJzP7AduGhGmwuBTwFk5u3ASRFxyrwqXcb8F+GQn/50OOkSyvDv4hD7\nYvL6hMVaYM/Y+t5uW6vNuh7nkiQVsLrHMcc6qXDmY8LDjnvNa3qcfRl64AH45jcnXcXkPf54jYfL\nT3saXH893HLLZOvw7+IQ+2LyIuc4oTwizgamM3OqW78CeDIzPzDW5i+BYWZu79Z3Audm5iNjbZzJ\nLkk9ZOZx/1jVZ2RxJ7AxIjYADwMXA1tmtLkR2Aps78Llx+NBAZO5WElSP3MOi8w8EBFbgZuBVcB1\nmXl/RFzW7d+WmTdFxAURsRt4HHjDglYtSTqu5nwbSpK08jRnQ83nC3hHOjYiTo6IL0fEdyLiSxFx\n0ti+K7r2OyPiN8e2vywi7u32/dnY9qdHxN9122+LiF/s0xHHYgn0xR9FxH3dub8SEacufC8c/Xpm\ntJlYX4zt/52IeDIizly4qz/sHOX7IiJe2/1tfDsi/nZhe+DnzlO6LyLitIj4ekTc1Z3/1QvfC0e/\nnhltjkdfvC8ivh8RP5lx7rm9d2bmEX8Y3WbaDWwATgC+BZw+o80FwE3d8lnAba1jgQ8C7+iWLwfe\n3y2f0bU7oTtuN4dGP3cAm7vlm4CpbvmtwF90yxcD2492TX1/lkhfDIBndMt/uJL7olt/NnAL8C/A\nmSu1L4CNwL8Bz+3W16zgvvgkcFm3fDrwvRXQF5uB5wM/mXH+Ob13tkYWfb+A9/zGsQeP6X7/Vrd8\nEfCZzNyfmQ91F3xWRPwC8OzMvKNr99djx4y/1meBVzWuqa/yfZGZw8z8abf9dhbvuy3l+6LzHuD9\nwP9y+FTuhbIU+uLNwNWZ+VhXw6MLc+mHWQp98QPgud3yScC++V/2rEr0Rffad2Tmf8xS45zeO1th\n0fcLeGuBFxzl2FPy0OyoR4Cnvt39gq7dbK81vn3f2GsdPH9mHgAei4iTG9fVx1Loi3FvYvSJajGU\n74vuttPazHyqDxbr4Vz5vmA0snhRRPxzRNwaEecf26XN2VLoiz8FLomIPcA/AW87lgvroUpfHFON\nx/Le2QqLvl/AO1Kbw14vR2OgpfCUfcn0RUS8HjgT+NB8X+sISvdFRATwUeCP51hLH6X7onMCcBpw\nLqNp7tdGxHOPfkgvS6EvPgr8VWauZ3Qb6G/m8VpHU6EvFvR9tRUW+4D1Y+vr+fn0mq3Nuq7NbNuf\nGvI90g236IaM/3kMr7Vulu1PHXNq91qrGd2X/WHjuvqo3BcHh9IRcR7wLuDCbgi7GCr3xV5Gzyp+\nBRhGxPeAs4EbF+khd/W+gNGnxy9k5v91tyi+wyg8FtpS6ItzgL8HyMzbgGdExJpju7w5qdAXrVts\nc3vvbDykWQ18l9EDkxNpP6Q5m0MPaY54LKOHNJd3y+/k8Ic0JwK/1B3/1EOa2xndgwsOf8B9Tbf8\nuyzeQ92l0BcvZXSv8oWL0QdLqS9m1PI1Fu8Bd/m+AM4HPtktrwG+DzxvhfbF54BLuuXTgX3L/e9i\n7HyzPeA+5vfOY7noVwMPMHoTuqLbdhndjIJu/epu/92M/Us527Hd9pOBrzD6hPMl4KSxfe/q2u8E\nzh/b/jLg3m7fx8e2P53RJ4VdwG3AhsX4h79E+uLLjB7g3dX9fH6l9sWMWhctLJZKXwAfAe4D7gFe\nu1L7AnghMGT0xnoXcN4K6IsPMhpdHuh+v7vbPqf3Tr+UJ0lq8n+rKklqMiwkSU2GhSSpybCQJDUZ\nFpKkJsNCktRkWEiSmgwLSVLT/wP8H5psQfIYVQAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "#Data Rate = 1MHz\n", "d = np.zeros(100000, dtype=float)\n", "d[40000:50000] = np.ones(10000, dtype=float)\n", "\n", "plot(n, d)\n", "ylim(0, 1.2)\n" ] }, { "cell_type": "code", "execution_count": 158, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "[]" ] }, "execution_count": 158, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYkAAAEACAYAAABGYoqtAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAADcBJREFUeJzt3G+sZHV9x/H3p7tAWm1YkWb5t4qtmMKTFm2A2jbeNpSu\n2EBNbIGkSmxSTRvSNjUVsE3YPmptYtoaU0SjhhgrGGvoEqmI1tuaNKJEJFtlC6uSLlhWhaKtPijI\ntw/mAOP1fu+u9+ydGXber2TCzJkzc875MXPfc86Z2VQVkiSt50fmvQKSpMVlJCRJLSMhSWoZCUlS\ny0hIklpGQpLUGh2JJO9NcijJvg3meXuS+5Pck+TcscuUJM3G0diTeB+wu7szycXAi6vqLOANwPVH\nYZmSpBkYHYmq+jTw3xvMcglw4zDvncCOJDvHLleStPVmcU7idODg1O0HgTNmsFxJ0kizOnGdNbf9\nt0Ak6Vlg+wyW8RCwa+r2GcO075PEcEjSJlTV2g/iR80s9iT2Aq8DSHIB8FhVHVpvxqryUsV11103\n93VYhMtkh/O64b/zvcx7LHxdOBbdZauN3pNI8kHgFcDJSQ4yeVcfB1BVN1TVbUkuTnIA+A7w+rHL\nlCTNxuhIVNUVRzDPVWOXI0maPX9xvYBWVlbmvQoLZGXeK7AwfF08w7GYnczimNaRSFKLsi5aDNmy\nU3E/PF+aWlRJqGf5iWtJ0rOUkZAktYyEJKllJCRJLSMhSWoZCUlSy0hIklpGQpLUMhKSpJaRkCS1\njIQkqWUkJEktIyFJahkJSVLLSEiSWkZCktQyEpKklpGQJLWMhCSpZSQkSS0jIUlqGQlJUstISJJa\nRkKS1DISkqSWkZAktYyEJKllJCRJLSMhSWoZCUlSy0hIklpGQpLUMhKSpJaRkCS1jIQkqWUkJEkt\nIyFJahkJSVLLSEiSWkZCktQyEpKklpGQJLWMhCSpZSQkSS0jIUlqjY5Ekt1J9ie5P8nV69y/kuRb\nSe4eLn82dpmSpNnYPubBSbYB7wAuBB4CPpdkb1Xdu2bWf6mqS8YsS5I0e2P3JM4DDlTVA1X1OHAT\ncOk682XkciRJczA2EqcDB6duPzhMm1bAy5Pck+S2JOeMXKYkaUZGHW5iEoDD+Tywq6q+m+SVwC3A\nS9abcc+ePU9fX1lZYWVlZeTqSdKxZXV1ldXV1ZktL1VH8ne+eXByAbCnqnYPt68Fnqyqt27wmK8C\nL6uqR9dMrzHromNPFuggpS9NLaokVNWWvVvGHm66CzgryZlJjgcuA/ZOz5BkZzJ5uyc5j0mYHv3B\np5IkLZpRh5uq6okkVwG3A9uA91TVvUneONx/A/Aa4PeSPAF8F7h85DpLkmZk1OGmo8nDTVrLw03S\n4S364SZJ0jHMSEiSWkZCktQyEpKklpGQJLWMhCSpZSQkSS0jIUlqGQlJUstISJJaRkKS1DISkqSW\nkZAktYyEJKllJCRJLSMhSWoZCUlSy0hIklpGQpLUMhKSpJaRkCS1jIQkqWUkJEktIyFJahkJSVLL\nSEiSWkZCktQyEpKklpGQJLWMhCSpZSQkSS0jIUlqGQlJUstISJJaRkKS1DISkqSWkZAktYyEJKll\nJCRJLSMhSWoZCUlSy0hIklpGQpLUMhKSpJaRkCS1jIQkqWUkJEmt0ZFIsjvJ/iT3J7m6meftw/33\nJDl37DIlSbMxKhJJtgHvAHYD5wBXJDl7zTwXAy+uqrOANwDXj1mmJGl2xu5JnAccqKoHqupx4Cbg\n0jXzXALcCFBVdwI7kuwcuVxJ0gxsH/n404GDU7cfBM4/gnnOAA6tfbJbbx25NtIW8bWpZTU2EnWE\n8+VIHvemN+15+vrzn7/CySevbGqlpKPtXe+a9xpIE9/85iqPPLI6s+WNjcRDwK6p27uY7ClsNM8Z\nw7QfcN99e0aujo4lWfvRYo7ck9DiWBkuE8mfb+nSxp6TuAs4K8mZSY4HLgP2rplnL/A6gCQXAI9V\n1Q8capIkLZ5RexJV9USSq4DbgW3Ae6rq3iRvHO6/oapuS3JxkgPAd4DXj15rSdJMpOpITytsrSS1\nKOuixbBIh5t8aWpRJaGqtuzd4i+uJUktIyFJahkJSVLLSEiSWkZCktQyEpKklpGQJLWMhCSpZSQk\nSS0jIUlqGQlJUstISJJaRkKS1DISkqSWkZAktYyEJKllJCRJLSMhSWoZCUlSy0hIklpGQpLUMhKS\npJaRkCS1jIQkqWUkJEktIyFJahkJSVLLSEiSWkZCktQyEpKklpGQJLWMhCSpZSQkSS0jIUlqGQlJ\nUstISJJaRkKS1DISkqSWkZAktYyEJKllJCRJLSMhSWoZCUlSy0hIklpGQpLUMhKSpJaRkCS1tm/2\ngUlOAm4GXgg8APxWVT22znwPAN8Gvgc8XlXnbXaZkqTZGrMncQ1wR1W9BPjkcHs9BaxU1bkGQpKe\nXcZE4hLgxuH6jcBvbDBvRixHkjQnYyKxs6oODdcPATub+Qr4RJK7kvzuiOVJkmZsw3MSSe4ATlnn\nrj+dvlFVlaSap/mFqvqvJD8B3JFkf1V9er0Z9+zZ8/T1lZUVVlZWNlo9SVo6q6urrK6uzmx5qer+\nth/mgcl+JucaHk5yKvCpqvrpwzzmOuB/q+pt69xXm10XHZuyQAcpfWlqUSWhqrbs3TLmcNNe4Mrh\n+pXALWtnSPJjSX58uP4c4CJg34hlSpJmaMyexEnAh4AXMPUV2CSnAe+uqlcl+UngI8NDtgMfqKq/\naJ7PPQl9H/ckpMPb6j2JTUfiaDMSWstISIe3yIebJEnHOCMhSWoZCUlSy0hIklpGQpLUMhKSpJaR\nkCS1jIQkqWUkJEktIyFJahkJSVLLSEiSWkZCktQyEpKklpGQJLWMhCSpZSQkSS0jIUlqGQlJUstI\nSJJaRkKS1DISkqSWkZAktYyEJKllJCRJLSMhSWoZCUlSy0hIklpGQpLUMhKSpJaRkCS1jIQkqWUk\nJEktIyFJahkJSVLLSEiSWkZCktQyEpKklpGQJLWMhCSpZSQkSS0jIUlqGQlJUstISJJaRkKS1DIS\nkqTWpiOR5DeTfDHJ95K8dIP5difZn+T+JFdvdnmSpNkbsyexD3g18K/dDEm2Ae8AdgPnAFckOXvE\nMpfC6urqvFdhgazOewUWhq+LZzgWs7PpSFTV/qq67zCznQccqKoHqupx4Cbg0s0uc1n4Bpi2Ou8V\nWBi+Lp7hWMzOVp+TOB04OHX7wWGaJOlZYPtGdya5AzhlnbveUlW3HsHz16bWSpK0EFI17u94kk8B\nb6qqz69z3wXAnqraPdy+Fniyqt66zrwGRZI2oaqyVc+94Z7ED6FbwbuAs5KcCXwNuAy4Yr0Zt3Ij\nJUmbM+YrsK9OchC4APhokn8app+W5KMAVfUEcBVwO/Al4Oaqunf8akuSZmH04SZJ0rFr7r+4PhZ/\nbJdkV5JPDT82/PckfzBMPynJHUnuS/LxJDumHnPtMAb7k1w0Nf1lSfYN9/3t1PQTktw8TP9MkhfO\ndit/OEm2Jbk7ya3D7aUciyQ7knw4yb1JvpTk/CUei2uH98i+JH8/rPtSjEWS9yY5lGTf1LSZbHuS\nK4dl3JfkdYdd2aqa2wXYBhwAzgSOA74AnD3PdTpK23UK8LPD9ecC/wGcDfwV8OZh+tXAXw7Xzxm2\n/bhhLA7wzF7eZ4Hzhuu3AbuH678P/N1w/TLgpnlv92HG5I+BDwB7h9tLORbAjcDvDNe3Aycu41gM\n2/MV4ITh9s3AlcsyFsAvAecC+6ambfm2AycBXwZ2DJcvAzs2XNc5D9TPAx+bun0NcM28/wduwXbe\nAlwI7Ad2DtNOAfYP168Frp6a/2NMzvWcCtw7Nf1y4J1T85w/XN8OfGPe27nB9p8BfAL4ZeDWYdrS\njQWTIHxlnenLOBYnMfnw9LxhPW8FfnWZxoLJH/zpSGz5tjP54tD1U495J3D5Rus578NNx/yP7YZv\ndp0L3MnkBXBouOsQsHO4fhqTbX/KU+OwdvpDPDM+T49dTb4g8K0kJx39LTgq/hr4E+DJqWnLOBYv\nAr6R5H1JPp/k3UmewxKORVU9CrwN+E8m33x8rKruYAnHYspWb/vzN3iu1rwjcUyfNU/yXOAfgD+s\nqv+Zvq8mGT+mtx8gya8DX6+qu2m+Kr0sY8HkE91LmRwGeCnwHSZ7z09blrFI8lPAHzH5NH0a8Nwk\nvz09z7KMxXoWadvnHYmHgF1Tt3fx/ZV71kpyHJNAvL+qbhkmH0pyynD/qcDXh+lrx+EMJuPw0HB9\n7fSnHvOC4bm2AycOn84WzcuBS5J8Ffgg8CtJ3s9yjsWDwINV9bnh9oeZROPhJRyLnwP+raoeGT7p\nfoTJ4edlHIunbPV74pF1nuuwf3PnHYmnf2yX5HgmJ1j2znmdRksS4D3Al6rqb6bu2svk5BzDf2+Z\nmn55kuOTvAg4C/hsVT0MfHv4BkyA1wL/uM5zvQb45JZt0AhV9Zaq2lVVL2JyzPSfq+q1LOdYPAwc\nTPKSYdKFwBeZHI9fqrFgcvz9giQ/OmzDhUx+S7WMY/GUWbwnPg5clMm37J7H5DzQ7Ruu1QKcvHkl\nkxNYB4Br570+R2mbfpHJ8fcvAHcPl91MTtZ9Arhv+J+1Y+oxbxnGYD/wa1PTX8bkn2U/ALx9avoJ\nwIeA+4HPAGfOe7uPYFxewTPfblrKsQB+BvgccA+TT88nLvFYvJlJJPcx+dbXccsyFkz2qr8G/B+T\ncwevn9W2D8u6f7hcebh19cd0kqTWvA83SZIWmJGQJLWMhCSpZSQkSS0jIUlqGQlJUstISJJaRkKS\n1Pp/kglqtpDY6VIAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# Data and carrier\n", "x = d * c\n", "\n", "plot(x)" ] }, { "cell_type": "code", "execution_count": 159, "metadata": { "collapsed": false }, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "/home/emily/anaconda3/lib/python3.4/site-packages/numpy/core/numeric.py:462: ComplexWarning: Casting complex values to real discards the imaginary part\n", " return array(a, dtype, copy=False, order=order)\n" ] }, { "data": { "text/plain": [ "[]" ] }, "execution_count": 159, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYgAAAENCAYAAAAVPvJNAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAGQJJREFUeJzt3X+wX3V95/HnCwIYayRlsAghCNKrJf5Am0pa19XrUNPo\nuEC3rsB2aVpTpyW2OjudtgE75TodFXdnqjgOmepSCcyQykoHscRASrnV2TVELUgwZkm6DUMiREUB\nf7WbkPf+8T03+eZ6yL253+/9kS/Px8ydfM7nfM75vr8n935f33PO95xvqgpJksY7brYLkCTNTQaE\nJKmVASFJamVASJJaGRCSpFYGhCSpVV8CIsnxSe5P8vlm+pQkm5I8nOTuJAu7xl6VZEeS7UmWd/Uv\nTbK1mXddP+qSJE1dv/Yg3gdsA8YuqlgDbKqqlwH3NNMkWQJcCiwBVgDXJ0mzzFpgVVUNAUNJVvSp\nNknSFPQcEEnOBN4G/A9g7MX+ImBd014HXNK0LwbWV9W+qtoF7ASWJTkdWFBVW5pxN3UtI0maBf3Y\ng/go8MfAga6+06pqb9PeC5zWtM8AdneN2w0saunf0/RLkmZJTwGR5O3At6vqfg7tPRymOvfy8H4e\nknSMmdfj8q8HLkryNuB5wAuT3AzsTfLiqnq8OXz07Wb8HmBx1/Jn0tlz2NO0u/v3jH+wJAaNJE1B\nVbW+iT+SnvYgqurqqlpcVecAlwH/UFVXAHcAK5thK4Hbm/YdwGVJTkxyDjAEbKmqx4GnkyxrTlpf\n0bXM+Mf0p4prrrlm1muYKz9uC7eF2+LIP1PV6x7EeGOVXAvcmmQVsAt4J0BVbUtyK51PPO0HVteh\n6lcDNwLzgQ1VtbHPtUmSjkLfAqKq/hH4x6b9PeBXn2Xch4APtfR/DXhVv+qRJPXGK6mPUcPDw7Nd\nwpzhtjjEbXGI26J36eX41ExLUsdSvZI0FyShZvoktSRpcA1cQPzwh7DR09uSBtyDD8LDD0/vYwxc\nQKxdC29962xXIUnT6/zz4fWvn97HGLiAkKTnihz1WYWjY0BIkloZEJKkVgMXENO9yyVJzxUDFxCS\npP4wICRJrQYuIDzEJEn9MXABIUnqDwNCktTKgJAktRq4gPAchCT1x8AFxCOPzHYFkjQYBu77IBYs\n6NzR9Rh6WpJ01BI49VT4zncmM9bvgwBg377ZrkCSBkNPAZHkeUnuS/JAkm1JPtz0n5JkU5KHk9yd\nZGHXMlcl2ZFke5LlXf1Lk2xt5l031Zr+7d96eUaSpDE9BURV/Svw5qp6DfBq4M1J3gCsATZV1cuA\ne5ppkiwBLgWWACuA65ODp5XXAquqaggYSrKil9okSb3p+RBTVf24aZ4IHA98H7gIWNf0rwMuadoX\nA+ural9V7QJ2AsuSnA4sqKotzbibupaRJM2CngMiyXFJHgD2AvdW1TeA06pqbzNkL3Ba0z4D2N21\n+G5gUUv/nqZfkjRL5vW6gqo6ALwmycnAXUnePG5+JfEzRZJ0jOk5IMZU1VNJ7gSWAnuTvLiqHm8O\nH327GbYHWNy12Jl09hz2NO3u/j1tjzMyMnKwPTw8zPDwcL+egiQNhNHRUUZHR3teT0/XQSQ5Fdhf\nVU8mmQ/cBXwA+DXgiar6SJI1wMKqWtOcpL4FuIDOIaS/B36+2cu4D3gvsAW4E/h4VW0c93gTXgcx\ndsrb6yAkDbKZuA6i1z2I04F1SY6jcz7j5qq6J8n9wK1JVgG7gHcCVNW2JLcC24D9wOquV/zVwI3A\nfGDD+HCQJM2sgbuS2j0ISc8FXkktSZo1BoQkqZUBIUlqZUBIkloZEJKkVgaEJKmVASFJamVASJJa\nGRCSpFYGhCSplQEhSWplQEiSWhkQkqRW3s1Vko5BR/Na591cJUl9ZUBIkloZEJKkVgaEJKmVASFJ\natVTQCRZnOTeJN9I8lCS9zb9pyTZlOThJHcnWdi1zFVJdiTZnmR5V//SJFubedf1UpckqXe97kHs\nA/5rVb0C+GXgPUnOA9YAm6rqZcA9zTRJlgCXAkuAFcD1ydiHtVgLrKqqIWAoyYoea5Mk9aCngKiq\nx6vqgab9Q+CbwCLgImBdM2wdcEnTvhhYX1X7qmoXsBNYluR0YEFVbWnG3dS1jCRpFvTtHESSs4HX\nAvcBp1XV3mbWXuC0pn0GsLtrsd10AmV8/56mX5I0S+b1YyVJXgDcBryvqn5w6KgRVFUl6dt1zSMj\nIwfbw8PDDA8P92vVkjQQRkdHGR0d7Xk9Pd9qI8kJwN8BX6iqjzV924Hhqnq8OXx0b1X9QpI1AFV1\nbTNuI3AN8Egz5rym/3LgTVX1++Mey1ttSBLHwK02mhPMNwDbxsKhcQewsmmvBG7v6r8syYlJzgGG\ngC1V9TjwdJJlzTqv6FpGkjQLetqDSPIG4IvAg8DYiq4CtgC3AmcBu4B3VtWTzTJXA+8C9tM5JHVX\n078UuBGYD2yoqve2PJ57EJLEzOxBeDdXSToGzflDTJKkwWVASJJaGRCSpFYGhCSplQEhSWplQEiS\nWhkQkqRWBoQkqZUBIUlqZUBIkloZEJKkVgaEJKmVASFJamVASJJaGRCSpFYGhCSplQEhSWplQEiS\nWhkQkqRWPQdEkr9OsjfJ1q6+U5JsSvJwkruTLOyad1WSHUm2J1ne1b80ydZm3nW91iVJ6k0/9iA+\nDawY17cG2FRVLwPuaaZJsgS4FFjSLHN9MvbV26wFVlXVEDCUZPw6JUkzqOeAqKovAd8f130RsK5p\nrwMuadoXA+ural9V7QJ2AsuSnA4sqKotzbibupaRJM2C6ToHcVpV7W3ae4HTmvYZwO6ucbuBRS39\ne5p+SdIsmTfdD1BVlaT6tb6RkZGD7eHhYYaHh/u1akkaCKOjo4yOjva8nlT1/tqd5Gzg81X1qmZ6\nOzBcVY83h4/urapfSLIGoKqubcZtBK4BHmnGnNf0Xw68qap+f9zj1ET1jp3R6MPTkqQ562he65JQ\nVZl45OGm6xDTHcDKpr0SuL2r/7IkJyY5BxgCtlTV48DTSZY1J62v6FpGkjQLej7ElGQ98Cbg1CSP\nAn8OXAvcmmQVsAt4J0BVbUtyK7AN2A+s7tolWA3cCMwHNlTVxl5rkyRNXV8OMc0UDzFJUsexfIhJ\nknSMMyAkSa0MCElSKwNCktTKgJAktTIgJEmtDAhJUisDQpLUyoCQJLUyICRJrQwISVIrA0KS1MqA\nkCS1MiAkSa0MCElSKwNCktTKgJAktRrYgPAb5SSpNwMbEAcOzHYFknRsm1MBkWRFku1JdiT5017W\n9eij/apKkp6b5kxAJDke+ASwAlgCXJ7kvKmu78/+rF+VSdLc8thjh9rf//70PU5qjhysT/IrwDVV\ntaKZXgNQVdd2jalnq/fAATj++MP7nnkGjpszEShJ/ZEcPv3jH8P8+UcaH6oqzz6i3byjXWAaLQK6\nDwztBpaNHzR+wxzJ+MCYbgsXwsgIrF8P99330/Pnz4dPfKITZt/+Njz9NLzqVTBvHvzMz3TGJPCj\nH8Hu3Z2AW7QIXvCCydeQdE7QP9t26s7X8WPG5h3NNpb67Ui/vxPN755X1XmDeOBAp+/AgUNvGH/4\nw86L6sknw09+Al/6Umf6wgvhO9+Bz362M+7cczt/h1Wd5Q8cODQ9fz68/e2wcmX/nvtUPf/507Pe\nuRQQk9yVGelqDzc/c8OTT8Kdd7aHA3R+EW+/Hfbvhy98odP36ld3/nMXLuz88j39NGzeDC99aSc4\nzj238+9kVB05ANrGdxsLl4mWn+gPWIOl1zcOk33DMl73m53xL/wTLTMWDs880+nrDoqqzhu0+fNh\n2zY46yy4//7O8l/7WucN2tg5zK985cjPbcOGI8+fPaPNT2/mUkDsARZ3TS+msxdxmKqRw6a/+lV4\n3euefaUPPQSveEVf6pOkWbVpEyxf3j7v7LPhttvgla+EE04YJhk+OC/5wJQeby6dg5gH/B/gQuBb\nwBbg8qr6ZteYZz0HAT99HmKOPDVJ6qvuPbJPfhLe/e6Jxh/j5yCqan+SPwDuAo4HbugOh8noPiH9\nwQ/2tTxJmjOeeebQm+GJwqEXc2YPYjIm2oPojOn8u3cv/NzPzUBRkjQLJnNO5tDYqe1BDOyHQBcs\nmO0KJOnYNrAB8bznzXYFknRsG9iA8GOYktSbgQ0ISVJvDAhJUisDQpLUyoCQJLUyICRJrQwISVIr\nA0KS1MqAkCS1MiAkSa0MCElSKwNCktTKgJAktTIgJEmtDAhJUisDQpLUyoCQJLWackAk+U9JvpHk\nmSS/OG7eVUl2JNmeZHlX/9IkW5t513X1n5TkM03/5iQvmWpdkqT+6GUPYivw68AXuzuTLAEuBZYA\nK4Drk4Pf77YWWFVVQ8BQkhVN/yrgiab/o8BHeqhLktQHUw6IqtpeVQ+3zLoYWF9V+6pqF7ATWJbk\ndGBBVW1pxt0EXNK0LwLWNe3bgAunWpckqT+m4xzEGcDurundwKKW/j1NP82/jwJU1X7gqSSnTENt\nkqRJmnekmUk2AS9umXV1VX1+eko6spGRkYPt4eFhhoeHZ6MMSZqzRkdHGR0d7Xk9qareVpDcC/xR\nVf1TM70GoKqubaY3AtcAjwD3VtV5Tf/lwBur6spmzEhVbU4yD3isql7U8lg1Ub1jZzt6fFqSNKcd\nzWtdEqoqE488XL8OMXU/8B3AZUlOTHIOMARsqarHgaeTLGtOWl8BfK5rmZVN+x3APX2qS5I0RUc8\nxHQkSX4d+DhwKnBnkvur6q1VtS3JrcA2YD+wuutt/2rgRmA+sKGqNjb9NwA3J9kBPAFcNtW6JEn9\n0fMhppnkISZJ6jiWDjFJkgaMASFJamVASJJaGRCSpFYGhCSplQEhSWplQEiSWhkQkqRWBoQkqZUB\nIUlqZUBIkloZEJKkVgaEJKmVASFJamVASJJaGRCSpFYGhCSplQEhSWplQEiSWk05IJL89yTfTPL1\nJH+b5OSueVcl2ZFke5LlXf1Lk2xt5l3X1X9Sks80/ZuTvGTqT0mS1A+97EHcDbyiqs4HHgauAkiy\nBLgUWAKsAK5Pxr5em7XAqqoaAoaSrGj6VwFPNP0fBT7SQ12SpD6YckBU1aaqOtBM3gec2bQvBtZX\n1b6q2gXsBJYlOR1YUFVbmnE3AZc07YuAdU37NuDCqdYlSeqPfp2DeBewoWmfAezumrcbWNTSv6fp\np/n3UYCq2g88leSUPtUmSZqCeUeamWQT8OKWWVdX1eebMe8H/l9V3TIN9f2UkZGRg+3h4WGGh4dn\n4mEl6ZgxOjrK6Ohoz+tJVU194eS3gXcDF1bVvzZ9awCq6tpmeiNwDfAIcG9Vndf0Xw68saqubMaM\nVNXmJPOAx6rqRS2PVxPVO3a2o4enJUlz3tG81iWhqjLxyMP18immFcAfAxePhUPjDuCyJCcmOQcY\nArZU1ePA00mWNSetrwA+17XMyqb9DuCeqdYlSeqPKe9BJNkBnAh8r+n6clWtbuZdTee8xH7gfVV1\nV9O/FLgRmA9sqKr3Nv0nATcDrwWeAC5rTnCPf0z3ICSJmdmD6OkQ00wzICSpY04fYpIkDTYDQpLU\nyoCQJLUyICRJrQwISVIrA0KS1MqAkCS1MiAkSa0GLiBOOGG2K5CkwTBwASFJ6g8DQpLUyoCQJLUy\nICRJrQwISVIrA0KS1MqAkCS1MiAkSa0MCElSqykHRJK/SPL1JA8kuSfJ4q55VyXZkWR7kuVd/UuT\nbG3mXdfVf1KSzzT9m5O8ZOpPSZLUD73sQfy3qjq/ql4D3A5cA5BkCXApsARYAVyfjH17KmuBVVU1\nBAwlWdH0rwKeaPo/Cnykh7okSX0w5YCoqh90Tb4A+G7TvhhYX1X7qmoXsBNYluR0YEFVbWnG3QRc\n0rQvAtY17duAC6dalySpP+b1snCSDwJXAD8BLmi6zwA2dw3bDSwC9jXtMXuafpp/HwWoqv1Jnkpy\nSlV9r5f6JElTd8Q9iCSbmnMG43/+A0BVvb+qzgI+DXxsJgqWJM2MI+5BVNVbJrmeW4ANTXsPsLhr\n3pl09hz2NO3x/WPLnAV8K8k84ORn23sYGRk52B4eHmZ4eHiSJUrSc8Po6Cijo6M9rydVNbUFk6Gq\n2tG0/xC4oKquaE5S30LnkNMi4O+Bn6+qSnIf8F5gC3An8PGq2phkNfCqqroyyWXAJVV1Wctj1kT1\nnngi7NsHU3xaknRMGPvoz2Re65JQVZl45OF6OQfx4SQvB54B/hm4EqCqtiW5FdgG7AdWd72qrwZu\nBOYDG6pqY9N/A3Bzkh3AE8BPhYMkaWZNeQ9iNrgHIUkdM7EH4ZXUkqRWBoQkqZUBIUnHqJNPnt71\nGxCSpFYGhCQdo3LUp52PjgEhSWplQEiSWhkQkqRWBoQkHaM8B3GUfvd3Z7sCSRoMAxcQL/HLSiWp\nLwYuICRJ/WFASNIxynMQkqRZMXABcdzAPSNJanf88dO7/l6+MGhOes974Jd+abarkKTptW0bnHTS\n9D7GwH1hkCTpcH5hkCSprwwISVKrngMiyR8lOZDklK6+q5LsSLI9yfKu/qVJtjbzruvqPynJZ5r+\nzUm83E2SZllPAZFkMfAW4JGuviXApcASYAVwfXLw07prgVVVNQQMJVnR9K8Cnmj6Pwp8pJe6ngtG\nR0dnu4Q5w21xiNviELdF73rdg/hL4E/G9V0MrK+qfVW1C9gJLEtyOrCgqrY0424CLmnaFwHrmvZt\nwIU91jXw/OU/xG1xiNviELdF76YcEEkuBnZX1YPjZp0B7O6a3g0saunf0/TT/PsoQFXtB57qPmQl\nSZp5R7wOIskm4MUts94PXAUs7x7ex7okSbNsStdBJHklcA/w46brTDp7BMuA3wGoqmubsRuBa+ic\np7i3qs5r+i8H3lhVVzZjRqpqc5J5wGNV9aKWx/UiCEmagqlcBzGlK6mr6iHgtLHpJP8CLK2q7yW5\nA7glyV/SOXQ0BGypqkrydJJlwBbgCuDjzSruAFYCm4F30Amftsd1L0WSZki/brVx8J19VW1Lciuw\nDdgPrO66/Hk1cCMwH9hQVRub/huAm5PsAJ4ALutTXZKkKTqmbrUhSZo5c/JK6iQrmovsdiT502cZ\n8/Fm/teTvHama5wpE22LJL/ZbIMHk/yvJK+ejTpnwmR+L5pxr0uyP8l/nMn6ZtIk/0aGk9yf5KEk\nozNc4oyZxN/IqUk2Jnmg2Ra/PQtlTrskf51kb5KtRxhzdK+bVTWnfoDj6Vw7cTZwAvAAcN64MW+j\nc4gKOifGN8923bO4LX4FOLlpr3gub4uucf8A/B3wG7Nd9yz+XiwEvgGc2UyfOtt1z+K2GAE+PLYd\n6BzGnjfbtU/Dtvj3wGuBrc8y/6hfN+fiHsQFwM6q2lVV+4C/oXPxXbeDF9ZV1X3AwiSnMXgm3BZV\n9eWqeqqZvI/OJ8oG0WR+LwD+EPgs8J2ZLG6GTWZb/GfgtqraDVBV353hGmfKZLbFY8ALm/YL6dy1\nYf8M1jgjqupLwPePMOSoXzfnYkAcvGiuMXah3URjBvGFcTLbotsqYMO0VjR7JtwWSRbReXFY23QN\n6gm2yfxeDAGnJLk3yVeTXDFj1c2syWyLTwGvSPIt4OvA+2aotrnmqF835+IXBk32j3r8R14H8cVg\n0s8pyZuBdwH/bvrKmVWT2RYfA9ZUVTX3/xrUj0VPZlucAPwindvWPB/4cpLNVbVjWiubeZPZFlcD\nD1TVcJJzgU1Jzq+qH0xzbXPRUb1uzsWA2AMs7ppezOG36GgbM3ah3qCZzLagOTH9KWBFVR1pF/NY\nNpltsRT4m+bekKcCb02yr6rumJkSZ8xktsWjwHer6ifAT5J8ETgfGLSAmMy2eD3wQYCq+ufmuq2X\nA1+dkQrnjqN+3ZyLh5i+SudOr2cnOZHOnWHH/4HfAfwWQJJfBp6sqr0zW+aMmHBbJDkL+Fvgv1TV\nzlmocaZMuC2q6qVVdU5VnUPnPMSVAxgOMLm/kc8Bb0hyfJLn0zkpuW2G65wJk9kW24FfBWiOub8c\n+L8zWuXccNSvm3NuD6Kq9if5A+AuOp9QuKGqvpnk95r5f1VVG5K8LclO4Ec0t/cYNJPZFsCfAz8L\nrG3eOe+rqgtmq+bpMslt8Zwwyb+R7c0tbB4EDgCfqqqBC4hJ/l58CPh0kq/TeVP8J1X1vVkrepok\nWQ+8CTg1yaN0bnF0Akz9ddML5SRJrebiISZJ0hxgQEiSWhkQkqRWBoQkqZUBIUlzwGRuttc19o1J\n/inJviS/MW7eyiQPNz+/1UtNBoQkzQ2fpnPDzcl4hM6XrN3S3ZnkFDoffb+g+bkmycKpFmRASNIc\n0HazvSTnJvlCcz+tLyZ5eTP2karaSucal26/BtxdVU9W1ZPAJiYfOj9lzl0oJ0k66JPA71XVzubr\nmq+nc3+tZ3MGh99qZKIbfB6RASFJc1CSF9D5vpf/2dwlAeDEmazBgJCkuek4OvdLmuib37pvh7EH\nGO6aXkznC7SmXIAkaY6pqqeBf0nyDoB0jP9K4fG3tb8LWJ5kYZKfBd7S9E2JASFJc0Bzs73/Dbw8\nyaNJfgf4TWBVkgeAh+h8K9zY964/CrwD+Kuxj8Y2t/v/C+ArwBbgA83J6qnV5M36JElt3IOQJLUy\nICRJrQwISVIrA0KS1MqAkCS1MiAkSa0MCElSKwNCktTq/wM7uyJgG6oBSgAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "X = np.fft.fft(x)\n", "k = np.linspace(0, 10000000000, size(X), dtype=float)\n", "\n", "plot(k, X)\n", "#xlim(0,600000000)" ] }, { "cell_type": "code", "execution_count": 160, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "[]" ] }, "execution_count": 160, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY4AAAEACAYAAACkvpHUAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAADg5JREFUeJzt3W2sZdVdx/HvTwbaaEiRYobCDEULVfAVJRlIfdGbSMoU\nE2hjhZKYksa0qME31kipRsYY0wfTJhJi09ZqMTalTWzIEDAtNVzFKFTSFrEywphOCmM7psXW2kbD\nw98XZzOcuZy5c/9335lz7tzvJzmZ/bD2XmuvnHN+Z+11zp1UFZIkrdWPzLsBkqTNxeCQJLUYHJKk\nFoNDktRicEiSWgwOSVLL6OBI8mdJDiV5dJUytyV5IskjSS4ZW6ckaX42YsTx58Duo+1MchVwQVVd\nCLwL+MgG1ClJmpPRwVFVDwD/tUqRq4E7hrIPAWck2T62XknSfJyIOY5zgSen1p8CdpyAeiVJx8GJ\nmhzPinX/zokkbVLbTkAdB4GdU+s7hm1HSGKYSNI6VNXKD+fH1YkYcewF3g6Q5HLgu1V1aFbBqvJR\nxa233jr3NizCYzIwvXX4d/KYd5t8XizGw7548TEPo0ccST4NvAE4K8mTTF7ppwJU1Uer6t4kVyXZ\nD/wAeMfYOiVJ8zM6OKrq+jWUuWlsPZKkxeAvxxfQ0tLSvJuwQJbm3YCF4fPiRfbFfGVe98hWSlKL\n0hYthsyY7vMpIh0pCXUSTo5Lkk4iBockqcXgkCS1GBySpBaDQ5LUYnBIkloMDklSi8EhSWoxOCRJ\nLQaHJKnF4JAktRgckqQWg0OS1GJwSJJaDA5JUovBIUlqMTgkSS0GhySpxeCQJLUYHJKkFoNDktRi\ncEiSWgwOSVKLwSFJajE4JEktBockqcXgkCS1GBySpBaDQ5LUYnBIkloMDklSi8EhSWoxOCRJLQaH\nJKnF4JAktRgckqQWg0OS1GJwSJJaDA5JUovBIUlqMTgkSS0GhySpxeCQJLUYHJKkFoNDktQyOjiS\n7E6yL8kTSW6esX8pyfeSfGV4/O7YOiVJ87NtzMFJTgFuB64ADgL/lGRvVT22oujfVtXVY+qSJC2G\nsSOOXcD+qjpQVc8AdwLXzCiXkfVIkhbE2OA4F3hyav2pYdu0Al6f5JEk9ya5eGSdkqQ5GnWrikko\nHMuXgZ1V9cMkbwLuAl47q+CePXsOLy8tLbG0tDSyeZJ0clleXmZ5eXmubUjVWt77j3Jwcjmwp6p2\nD+u3AM9X1QdWOebrwKVV9fSK7TWmLTr5ZMYNTp8i0pGSUFUndDpg7K2qh4ELk5yf5DTgOmDvdIEk\n25PJW0CSXUzC6umXnkqStBmMulVVVc8muQn4PHAK8ImqeizJjcP+jwJvBX4tybPAD4G3jWyzJGmO\nRt2q2kjeqtJK3qqSjm0z3qqSJG0xBockqcXgkCS1GBySpBaDQ5LUYnBIkloMDklSi8EhSWoxOCRJ\nLQaHJKnF4JAktRgckqQWg0OS1GJwSJJaDA5JUovBIUlqMTgkSS0GhySpxeCQJLUYHJKkFoNDktRi\ncEiSWgwOSVKLwSFJajE4JEktBockqcXgkCS1GBySpBaDQ5LUYnBIkloMDklSi8EhSWoxOCRJLQaH\nJKnF4JAktRgckqQWg0OS1GJwSJJaDA5JUovBIUlqMTgkSS0GhySpxeCQJLUYHJKkFoNDktRicEiS\nWgwOSVLL6OBIsjvJviRPJLn5KGVuG/Y/kuSSsXVKkuZnVHAkOQW4HdgNXAxcn+SiFWWuAi6oqguB\ndwEfGVOnJGm+xo44dgH7q+pAVT0D3Alcs6LM1cAdAFX1EHBGku0j65Ukzcm2kcefCzw5tf4UcNka\nyuwADq082d13j2yNTno+R6T5GxsctcZyWctx7373nsPLr3zlEmedtbSuRunk9bGPzbsF0nx9+9vL\nfOc7y3Ntw9jgOAjsnFrfyWREsVqZHcO2l3j88T0jm6OTSVZ+3MARhwRLw2Mi+f0T3oKxcxwPAxcm\nOT/JacB1wN4VZfYCbwdIcjnw3ap6yW0qSdLmMGrEUVXPJrkJ+DxwCvCJqnosyY3D/o9W1b1Jrkqy\nH/gB8I7RrZYkzU2q1jpNcXwlqUVpixbDrFtVPkWkIyWhqma8Wo4ffzkuSWoxOCRJLQaHJKnF4JAk\ntRgckqQWg0OS1GJwSJJaDA5JUovBIUlqMTgkSS0GhySpxeCQJLUYHJKkFoNDktRicEiSWgwOSVKL\nwSFJajE4JEktBockqcXgkCS1GBySpBaDQ5LUYnBIkloMDklSi8EhSWoxOCRJLQaHJKnF4JAktRgc\nkqQWg0OS1GJwSJJaDA5JUovBIUlqMTgkSS0GhySpxeCQJLUYHJKkFoNDktRicEiSWgwOSVKLwSFJ\najE4JEktBockqcXgkCS1GBySpBaDQ5LUYnBIklq2rffAJGcCnwFeDRwArq2q784odwD4b+A54Jmq\n2rXeOiVJ8zdmxPEe4L6qei3wN8P6LAUsVdUlhoYkbX5jguNq4I5h+Q7gzauUzYh6JEkLZExwbK+q\nQ8PyIWD7UcoV8MUkDyd554j6JEkLYNU5jiT3AWfP2PU70ytVVUnqKKf5uar6ZpKfAO5Lsq+qHphV\ncM+ePYeXl5aWWFpaWq15krTlLC8vs7y8PNc2pOpo7/fHODDZx2Tu4ltJXgXcX1U/c4xjbgX+p6o+\nNGNfrbctOjllxg1OnyLSkZJQVSd0OmDMraq9wA3D8g3AXSsLJPnRJKcPyz8GvBF4dESdkqQ5GzPi\nOBP4LHAeU1/HTXIO8PGq+oUkPwV8bjhkG/CpqnrfUc7niENHcMQhHds8RhzrDo6NZnBoJYNDOrbN\ndqtKkrQFGRySpBaDQ5LUYnBIkloMDklSi8EhSWoxOCRJLQaHJKnF4JAktRgckqQWg0OS1GJwSJJa\nDA5JUovBIUlqMTgkSS0GhySpxeCQJLUYHJKkFoNDktRicEiSWgwOSVKLwSFJajE4JEktBockqcXg\nkCS1GBySpBaDQ5LUYnBIkloMDklSi8EhSWoxOCRJLQaHJKnF4JAktRgckqQWg0OS1GJwSJJaDA5J\nUovBIUlqMTgkSS0GhySpxeCQJLUYHJKkFoNDktRicEiSWgwOSVKLwSFJall3cCT5pSRfS/Jcktet\nUm53kn1Jnkhy83rrkyQthjEjjkeBtwB/d7QCSU4Bbgd2AxcD1ye5aESdW8Ly8vK8m7BAlufdgIXh\n8+JF9sV8rTs4qmpfVT1+jGK7gP1VdaCqngHuBK5Zb51bhS+KacvzbsDC8HnxIvtivo73HMe5wJNT\n608N2yRJm9S21XYmuQ84e8au91bV3Ws4f62rVZKkhZWqce/tSe4H3l1VX56x73JgT1XtHtZvAZ6v\nqg/MKGvISNI6VFVOZH2rjjgajtboh4ELk5wP/AdwHXD9rIIn+sIlSesz5uu4b0nyJHA5cE+Svx62\nn5PkHoCqeha4Cfg88K/AZ6rqsfHNliTNy+hbVZKkraU14ljLj/mS3DbsfyTJJcc6NsmZSe5L8niS\nLyQ5Y2rfLUP5fUneOLX90iSPDvv+eGr7y5J8Ztj+YJJXd66vaxP0x28OP9J8JMkXk5y38b2w+P0w\ntf8Xkzy/2g9Wx9oMfZHk2uF58S9JPrWxPXBEPQvdF0kuSPJAkq8M9b9p43th9etZUeZE9MUfJvlG\nku+vqLv33llVa3oApwD7gfOBU4GvAhetKHMVcO+wfBnw4LGOBT4I/PawfDPw/mH54qHcqcNx+3lx\nhPQlYNewfC+we1j+deBPhuXrgDvXen3dxybpjyXg5cPyrx6P/tgM/TCsn87kx6r/ALxuCz8nLgS+\nDLxiWD9rC/fFJ4Ebh+WLgK9vgb7YxeSbst9fUX/rvbMz4ljLj/muBu4AqKqHgDOSnH2MYw8fM/z7\n5mH5GuDTVfVMVR0YLv6yJK8CTq+qLw3l/mLqmOlz/RXw843r61r4/qiq5ar632H7Q8COjbn0Iyx8\nPwz+AHg/8H8c/cscY22GvngncHtVfW9ow7c35tJfYjP0xTeBVwzLZwAHx1/2TAvRF8O5v1RV35rR\nxtZ7Zyc41vJjvqOVOWeVY7dX1aFh+RCwfVg+Zyg361zT2w9Onetw/TWZmP9ekjPXcG3rsRn6Y9qv\nMPm0tdEWvh+GW1PnVtUL13+8JvYWvi+YjDh+OsnfJ/nHJFeu7dLaNkNfvA+4IZMv+dwD/MZaLmwd\nFqUv1tTGtbx3doJjrS+2tXyay6zz1WSctFlm6zdNfyT5ZeB1wB+NPdcMC90PSQJ8GPitZlvWY6H7\nYnAqcAHwBiZfjf94klesfsi6bIa++DDwp1W1k8mtor8cca7VLEJfbOj7aic4DgI7p9Z3cmSqzSqz\nYygza/sLw8JDw5CMYVj5n2s4144Z21845rzhXNuY3Md9em2X17bI/XF4yJ3kCuC9wNXDUHejLXI/\nPMVkbuNngeUkX2fy9fG9x2mCfNH7AiafKu+uqueG2xiPMwmSjbYZ+uL1wGcBqupB4OVJzlrb5bUs\nQl8c6zZc772zMcGzDfh3JpMtp3HsCZ7LeXGC56jHMpnguXlYfg8vneA5DfjJ4fgXJngeYnLPLrx0\ncvwjw/LbOL6T45uhPy5hcn/zNVu5H1a05X6O3+T4wvcFcCXwyWH5LOAbwI9v0b74HHDDsHwRcPBk\nf15M1TdrcnzN753dDngT8G9M3oxuGbbdyPDNhGH99mH/I0y9QGcdO2w/E/gik08+XwDOmNr33qH8\nPuDKqe2XMvmz7vuB26a2v4zJJ4gngAeB84/HE2ET9cd9TCYAvzI87tqK/bCircctODZLXwAfAr4G\n/DNw7VbtC+A1TP788leH18cVW6AvPshk1Pns8O/vDdtb753+AFCS1OJ/HStJajE4JEktBockqcXg\nkCS1GBySpBaDQ5LUYnBIkloMDklSy/8Dj+9aLSoGIjYAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "#Data Rate = 10MHz\n", "d2 = np.zeros(100000, dtype=float)\n", "d2[40000:41000] = np.ones(1000, dtype=float)\n", "x2 = d2 * c\n", "\n", "plot(n, x2)\n" ] }, { "cell_type": "code", "execution_count": 161, "metadata": { "collapsed": false }, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "/home/emily/anaconda3/lib/python3.4/site-packages/numpy/core/numeric.py:462: ComplexWarning: Casting complex values to real discards the imaginary part\n", " return array(a, dtype, copy=False, order=order)\n" ] }, { "data": { "text/plain": [ "[]" ] }, "execution_count": 161, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYEAAAENCAYAAADpK9mHAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAFw5JREFUeJzt3X+0XWV95/H3l0QkCBiQNglJXERMlFCgwvCjnbbcUcGM\nyyEsSkVsLRV0tc0UXS6nmsjqkK7pUNBppzItLMsUxFlDOhFbQUUkdXJ1loMErUAkZEiEpCSVQBAq\n2lCS3O/8sfflnlzuj3PPvufX3e/XWmfds5/9nP08d6979+fsX8+OzESSVE+HdbsDkqTuMQQkqcYM\nAUmqMUNAkmrMEJCkGjMEJKnGKodARMyNiDsi4tGI2BIR50TEcRGxISIei4h7I2JuQ/01EbEtIrZG\nxAVV25cktW469gQ+DdydmScDpwFbgdXAhsxcBny9nCYilgOXAsuBFcCNEeHeiCR1SaUNcES8Fvjl\nzLwFIDMPZOY/ARcCt5XVbgMuKt+vBNZl5v7M3AFsB86u0gdJUuuqfgtfAjwTEbdGxN9HxM0R8Rpg\nXmbuKevsAeaV708AdjV8fhewsGIfJEktqhoCs4EzgBsz8wzgp5SHfoZlMS7FRGNTOG6FJHXJ7Iqf\n3wXsyswHyuk7gDXAUxExPzOfiogFwNPl/N3A4obPLyrLDhERBoMktSAzYyr1K+0JZOZTwJMRsaws\nejvwCPAl4PKy7HLgi+X7u4D3RMThEbEEWApsGmfZvjK55pprut6HXnm5LlwXrouJX62ouicAcBXw\nPyPicOAHwPuBWcD6iLgS2AG8u9ywb4mI9cAW4ACwKlvtuSSpssohkJkPAWeNMevt49S/Fri2aruS\npOq8Rr/HDQwMdLsLPcN1McJ1McJ1UU304tGYiPAokSRNUUSQnTwxLEnqb4aAJNWYISBJNWYISFKN\nGQKSVGOGgCTVmCEgSTVmCEhSjRkCklRjhoAk1ZghIEk1ZghIUo31bQg89BA880y3eyFJ7bVxIwwN\ntW/5fRsCP//zcOWV3e6FJLXXW98K3/pW+5bftyEA8OKL3e6BJLXfwYPtW3Zfh0A7d5EkqQ76OgTa\nmY6S1Cva+Yytvg4BSVI1hoAk9Tj3BCRJbWEISFKPc09AktQWhoAk9Tj3BCRJbWEISFKP++d/bt+y\nDQFJ6nE9fzgoImZFxPci4kvl9HERsSEiHouIeyNibkPdNRGxLSK2RsQF09G+JM1kPR8CwIeBLcBw\nV1cDGzJzGfD1cpqIWA5cCiwHVgA3RoR7I5LUJZU3wBGxCHgn8N+BKIsvBG4r398GXFS+Xwmsy8z9\nmbkD2A6cXbUPkqTWTMe38P8K/D7QOKbnvMzcU77fA8wr358A7GqotwtY2GrDg4OtflKSBBVDICLe\nBTydmd9jZC/gEJmZjBwmGrNKlT5I0kzXznMCsyt+/heBCyPincARwDER8T+APRExPzOfiogFwNNl\n/d3A4obPLyrLXmHt2rUvvx8YGGBgYKBiVyWpP40XAoODgwxWPCQSOU0RExHnAf8hM/9dRHwSeDYz\nr4+I1cDczFxdnhi+neI8wELg74A35qhORMToojHaK362MyElqdsi4AtfgIsvbqZukJljHpUZT9U9\ngdGGN8nXAesj4kpgB/BugMzcEhHrKa4kOgCsmnRrL0k1186t5LTtCUwn9wQkqRABn/88XHJJM3Wn\nvifgNfqS1ONuv719yzYEJKnHPf305HVaZQhIUo8bGpq8TqsMAUnqcZ4YHrNO8bMHuy9J0yITDjsM\n5sxpbjhpTwxL0gy0b1/7lm0ISFKNGQKSVGOGgCT1qP3729+GISBJPWrv3va3YQhIUo0ZApJUY4aA\nJNWYISBJNWYISFKNGQKSVGOGgCTVmCEgSTVmCEhSj9q1q/1tGAKS1KNeeqn9bRgCklRjhoAk9aiY\n0uNhWmMISFKNGQKS1KPcE5AktZUhIEk96v7729+GISBJPWr37va3YQhIUo/at6/9bRgCktSjNm9u\nfxuVQiAiFkfExoh4JCK+HxEfKsuPi4gNEfFYRNwbEXMbPrMmIrZFxNaIuKDqLyBJal3VPYH9wEcy\n8xTgXODfR8TJwGpgQ2YuA75eThMRy4FLgeXACuDGiHBvRJK6pNIGODOfyswHy/c/AR4FFgIXAreV\n1W4DLirfrwTWZeb+zNwBbAfOrtIHSZqphoba38a0fQuPiBOBtwD3A/Myc085aw8wr3x/AtA4Lt4u\nitCQJI1y333tb2P2dCwkIo4CvgB8ODNfiIbb3DIzIyIn+PiY89auXfvy+4GBAQYGBqajq5LUNybb\nExgcHGRwcLBSG5E50fa5iQVEvAr4MvDVzPyzsmwrMJCZT0XEAmBjZr45IlYDZOZ1Zb17gGsy8/5R\ny8zJ+jWcMxW7L0k9q3HYiGa2dRFBZk5psImqVwcF8FfAluEAKN0FXF6+vxz4YkP5eyLi8IhYAiwF\nNlXpgySpdZX2BCLil4BvAg8zclhnDcWGfT3wemAH8O7MfL78zCeAK4ADFIePvjbGct0TkFR7ndgT\nqHw4qB0MAUnqg8NBkqT+ZghIUo0ZApJUY4aAJNWYISBJNWYISFKNGQKSVGOGgCTVmCEgSTVmCEhS\njRkCklRjfRkCe/ZMXkeSNLm+DIF/+Idu90CSZoa+DAFJ0vQwBCSpxvoyBB58sNs9kKSZoS9DYOfO\nbvdAkmaGvgyBoaFu90CSZoa+DIEvf7nbPZCkmaEvnzE81eduSlI/8hnDTTAEJKl1fR8Cf/RH3e6B\nJPWvvg+BBx7odg8kqf0OHGjPcvs+BH70o273QJLa71Ofas9y++7E8NAQzJp1aFkP/gqSVMmLL8Kc\nOSPTs2fD/v0Tf6YWJ4YHB7vdA0lqv3vvPXS6XYeD+m5PIMbIuB78FSSpkrG2dT/5CbzmNRN9Zobv\nCWzePHb5unWd7YckdcNRR03/MrsSAhGxIiK2RsS2iPh4M5+55RY47bSx5733vdPZO0nqro9+dPx5\nEfDCC9PX1uzpW1RzImIW8OfA24HdwAMRcVdmPtpY7+Mfh40bm78ENAL27YMjjpjuHktS5/zar8Ed\nd0xc55hjRt5fdlnxBfmss1prr+PnBCLiF4BrMnNFOb0aIDOva6iTUK1ff/AHcPXV8OpXV1qMJLXV\njh1w3nnT9cTEqZ8T6EYIXAK8IzM/WE7/BnBOZl7VUKdyCNTFvHnFt4J582D+/OLy2ZUr4TOfgXPP\nhRNPhOefh4UL4frri8vOfvADuOEGOHgQtm2D++6Dj3ykqPczPwMvvQSPP14sb8GCQ9vbvr04MXXk\nkSMnqQ4eLPbCIoo9saOPHtkje+YZeOwxOOOM4nK3Z5+FY48t+nHkkUWdiOLk/tAQ/Mu/FJ8dPik2\n/PPAgVdeGjxaZlF/rBNqk31u9DLGm+6Uydodb14rn+nk8jrV1s6dcMUV439m5pp6CHT8cBBNb93X\nNrwfKF8abc+e4rVt20jZU0/BN75RvE4+GR59tNhVfOSRkTp/8idw2GFF3SOPLKYff7wIk507i+uR\nTz0VliwZ+UxmcXPevn3w05/C614Hzz0Hc+cWZUNDRYgccUQRKHPnwhNPFBv9O++EN76xCJgDB4q6\nw1c5HDxYbOAjivLhgBjeOA8NFddIT/R9Zbq+ywwH0ngDdw2XD5eNt3FqZhmN9Ub3f6KN3ni/a6vr\np1PL62Rbzz03/mdmlsHy1bpuhMBuYHHD9GJg1yurrZ3WRtesgfPPh5NOguOOKzZUwxseSWqngwdh\n797iC9bDD8OGDbB+/XQseYBDvyD/4ZSX0I3DQbOB/we8DfhHYBNwWeOJ4SqHg/bsgZ/92enoqSR1\nztAQ/MVfwIc+VGUpfXCfQGYeAH4P+BqwBfhfo68MKuqNvF58cfLlfuUrRV0DQFI/OuwwuOqq5p6c\neOuth24jh1+t6Ks7hv/2b+Hii8f+TA/+GpLUkr17i/NrY5loW9fKHcN9FQLFvFeW/fCHxZUskjRT\njLWte+GFie8anvHDRgBcfvkrywwASTPNWCMhtGPYiL7bEyjmHzrdg7+CJFWybRssW3Zo2WTbulrs\nCUhSHbz+9YdO33xze9rp+xAY6/CQJPW70UPefOAD7Wmn70PglFO63QNJ6l99HwKnntrtHkhS/+r7\nE8M92H1JmhZT3dZ5YliSNCWGgCTVmCEgSTXWlyFw+und7oEkzQx9GQJnn93tHkjSzNCXVwdt3Vo8\nMQu8OkjSzOXVQeN43eu63QNJmhn6ck+gqFP87MHuS9K0cE9AktRWhoAk1ZghIEk1ZghIUo0ZApJU\nY4aAJNWYISBJNWYISFKNGQKS1KOOOab9bRgCktSjRj9svh0MAUnqUZ0YFscQkKQaazkEIuJTEfFo\nRDwUEX8TEa9tmLcmIrZFxNaIuKCh/MyI2FzO+3TVzkuSqqmyJ3AvcEpmng48BqwBiIjlwKXAcmAF\ncGPEy2Ph3QRcmZlLgaURsaJC+5KkiloOgczckJlD5eT9wKLy/UpgXWbuz8wdwHbgnIhYABydmZvK\nep8DLmq1fUma6frpnMAVwN3l+xOAXQ3zdgELxyjfXZZLkrpk9kQzI2IDMH+MWZ/IzC+Vda4GXsrM\n26ezY2vXrn35/cDAAAMDA9O5eEnqecuWwX33jT9/cHCQwcHBSm1UerJYRPwW8EHgbZn5Ylm2GiAz\nryun7wGuAXYCGzPz5LL8MuC8zPydMZbrk8Uk1d5FF8Gddxbve+7JYuVJ3d8HVg4HQOku4D0RcXhE\nLAGWApsy8yngxxFxTnmi+H3AF1ttX5JU3YSHgybx34DDgQ3lxT/3ZeaqzNwSEeuBLcABYFXD1/pV\nwGeBOcDdmXlPhfYlaUYb3nJ+8IPta6PlECgv8xxv3rXAtWOUfxc4tdU2JamOjjqqfcv2jmFJ6nHt\nPPdpCEhSjRkCklRjhoAk9ajFi9vfhiEgST1q/li36k4zQ0CSaswQkKQaMwQkqcfFlAaCmBpDQJJ6\nVD8NJS1J6kOGgCT1qDlz2t+GISBJPeqkk9rfhiEgST3OsYMkqca8OkiS1BaGgCT1OPcEJKmGXv3q\n9rdhCEhSj1q0qP1tGAKSVGOGgCTVmCEgST2qnSeEhxkCktSjFiwofnp1kCTV0KxZxU9DQJJq6Nhj\n29+GISBJPc6xgyRJbWEISFKNGQKSVGOVQyAiPhoRQxFxXEPZmojYFhFbI+KChvIzI2JzOe/TVduW\nJFVTKQQiYjFwPrCzoWw5cCmwHFgB3Bjx8gVONwFXZuZSYGlErKjSviSpmqp7An8KfGxU2UpgXWbu\nz8wdwHbgnIhYABydmZvKep8DLqrYviTNeMP3C7RDyyEQESuBXZn58KhZJwC7GqZ3AQvHKN9dlkuS\numT2RDMjYgMwf4xZVwNrgAsaq09jvyRJHTBhCGTm+WOVR8TPAUuAh8rD/YuA70bEORTf8Bc3VF9E\nsQewu3zfWL57vLbXrl378vuBgQEGBgYm6qok1c7g4CCDg4OVlhE5DbeiRcQTwJmZ+aPyxPDtwNkU\nh3v+DnhjZmZE3A98CNgEfAW4ITPvGWN5OVm/hk81t/NOOknqtgj42Mfg+uubqRtk5pSOyky4JzAF\nL2+KM3NLRKwHtgAHgFUNW/RVwGeBOcDdYwWAJKlzpiUEMvMNo6avBa4do953gVOno01JUnXeMSxJ\nNWYISFKNGQKS1OPe9a72LdsQkKQeN3+su7WmiSEgSTVmCEhSjfV1CCxePHkdSdL4+joEjj662z2Q\npP7W1yEQDlknSZUYApJUY4aAJNVYX4eAJKma6RpFtOM++Ul485u73QtJaq93vAMWtvEZjNPyPIHp\n1szzBCRJh2rleQIeDpKkGjMEJKnGDAFJqjFDQJJqzBCQpBozBCSpxgwBSaoxQ0CSaswQkKQaMwQk\nqcYMAUmqMUNAkmrMEJCkGjMEJKnGKoVARFwVEY9GxPcj4vqG8jURsS0itkbEBQ3lZ0bE5nLep6u0\nLUmqruUQiIh/A1wInJaZPwf8l7J8OXApsBxYAdwY8fKDIG8CrszMpcDSiFhRpfN1MDg42O0u9AzX\nxQjXxQjXRTVV9gR+F/jjzNwPkJnPlOUrgXWZuT8zdwDbgXMiYgFwdGZuKut9DrioQvu14B/4CNfF\nCNfFCNdFNVVCYCnwKxHx7YgYjIh/VZafAOxqqLcLWDhG+e6yXJLUJRM+YzgiNgDzx5h1dfnZYzPz\n3Ig4C1gPvGH6uyhJapeWnzEcEV8FrsvMb5TT24FzgQ8AZOZ1Zfk9wDXATmBjZp5cll8GnJeZvzPG\nsn3AsCS1YKrPGJ5wT2ASXwTeCnwjIpYBh2fm3oi4C7g9Iv6U4nDPUmBTZmZE/DgizgE2Ae8DbpiO\nX0KS1JoqIXALcEtEbAZeAn4TIDO3RMR6YAtwAFiVI7sbq4DPAnOAuzPzngrtS5IqavlwkCSp/3X1\njuGIWFHeULYtIj4+Tp0byvkPRcRbOt3HTplsXUTEr5fr4OGI+FZEnNaNfrZbM38TZb2zIuJARFzc\nyf51UpP/HwMR8b3yhs3BDnexY5r4/zg+Iu6JiAfLdfFbXehmR0TELRGxpzwKM16d5rebmdmVFzCL\n4h6CE4FXAQ8CJ4+q806Kw0YA5wDf7lZ/e2Bd/ALw2vL9ipm4LppZDw31/jfwZeBXu93vLv5NzAUe\nARaV08d3u99dXBdrKe5bAjgeeBaY3e2+t2l9/DLwFmDzOPOntN3s5p7A2cD2zNyRxQ1nf01xo1mj\nC4HbADLzfmBuRMzrbDc7YtJ1kZn3ZeY/lZP3A4s63MdOaOZvAuAq4A7gmTHmzRTNrIv3Al/IzF0A\nmbm3w33slGbWxQ+BY8r3xwDPZuaBDvaxYzLz/wDPTVBlStvNbobAQuDJhunhm8omqzMTN37NrItG\nVwJ3t7VH3THpeoiIhRQbgJvKopl6UquZv4mlwHERsTEivhMR7+tY7zqrmXVxM3BKRPwj8BDw4Q71\nrRdNabtZ5eqgqpr95x19uehM/Kdv+ncqx2y6AvjX7etO1zSzHv4MWJ2ZWY5JNVMvJ25mXbwKOAN4\nG3AkcF9EfDszt7W1Z53XzLr4BPBgZg5ExEnAhog4PTNfaHPfelXT281uhsBuYHHD9GIOHVZirDqL\nyrKZppl1QXky+GZgRWZOtDvYr5pZD2cCf12OSXg88G8jYn9m3tWZLnZMM+viSWBvZu4D9kXEN4HT\ngZkWAs2si18E/jNAZv4gIp4A3gR8pyM97C1T2m5283DQdyhGEj0xIg6nGHl09D/yXZT3H0TEucDz\nmbmns93siEnXRUS8Hvgb4Dcyc3sX+tgJk66HzHxDZi7JzCUU5wV+dwYGADT3/3En8EsRMSsijqQ4\nCbilw/3shGbWxVbg7QDl8e83AY93tJe9Y0rbza7tCWTmgYj4PeBrFGf//yozH42I3y7nfyYz746I\nd5ZDUvwUeH+3+ttOzawL4D8CxwI3ld+C92fm2d3qczs0uR5qocn/j63lsCwPA0PAzZk540Kgyb+L\na4FbI+Ihii+3H8vMH3Wt020UEeuA84DjI+JJimF5XgWtbTe9WUySaszHS0pSjRkCklRjhoAk1Zgh\nIEk1ZghIUoc0M/hbQ91fiYi/j4j9EfGro+ZdHhGPla/frNInQ0CSOudWigEgm7ETuBy4vbEwIo6j\nuGT87PJ1TUTMbbVDhoAkdchYg79FxEkR8dVy/KdvRsSbyro7M3MzxT0gjd4B3JuZz2fm88AGmg+W\nV+jmsBGSJPhL4Lczc3v5+N0bKcaDGs8JHDpsxmQDTk7IEJCkLomIoyieFfL5ciQAgMM72QdDQJK6\n5zCKsX0me2pi49AOu4GBhunFFA9ZarkDkqQuyMwfA09ExCUAURj96NjRQ6Z/DbggIuZGxLHA+WVZ\nSwwBSeqQcvC3/wu8KSKejIj3A78OXBkRDwLfp3gy2PBztJ8ELgE+M3xZaTmM/H8CHgA2AX9YniBu\nrU8OICdJ9eWegCTVmCEgSTVmCEhSjRkCklRjhoAk1ZghIEk1ZghIUo0ZApJUY/8fDP4T//0C2FMA\nAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "X2 = np.fft.fft(x2)\n", "\n", "plot(k, X2)\n", "#xlim(0,600000000)" ] }, { "cell_type": "code", "execution_count": 162, "metadata": { "collapsed": false }, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "/home/emily/anaconda3/lib/python3.4/site-packages/numpy/core/numeric.py:462: ComplexWarning: Casting complex values to real discards the imaginary part\n", " return array(a, dtype, copy=False, order=order)\n" ] }, { "data": { "text/plain": [ "[]" ] }, "execution_count": 162, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXsAAAENCAYAAADuRcXXAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAF4NJREFUeJzt3XuQHOV57/HvA0LBcgpkFTEyNxsbIQIuMBALcHC8NjdZ\nCZACfDs2IQ4Qlx1T/oMTG0JVpBTxOTjGiU+c2GVjoDjEQDAkXBIulsFLTmJxswGbSBwhjCgJIxmb\nWxACdHnyR8+yq2W1OzszPT0z/f1UTW13T0/3s107v32n5+23IzORJA22HaouQJJUPsNekmrAsJek\nGjDsJakGDHtJqgHDXpJqoO2wj4jZEXFdRKyIiOURcUREzImIpRGxMiK+FxGzO1GsJKk1nWjZ/x/g\nlsz8TeBg4BHgPGBpZu4P3NGYlyRVJNq5qCoidgUeyMy3j1v+CPC+zFwfEXOB4cw8oL1SJUmtardl\nvy/wdERcHhE/johLIuKNwO6Zub6xznpg9zb3I0lqQ7thPwM4DPh6Zh4GbGDcKZssPjo4JoMkVWhG\nm69fC6zNzPsa89cB5wPrImJuZq6LiLcAvxj/wojwH4AktSAzY7qvaatln5nrgDURsX9j0bHAfwI3\nA2c0lp0B3LCd1/vIZPHixZXX0CsPj4XHwmMx+aNV7bbsAc4BvhMRM4HHgE8COwLXRsSZwGrgwx3Y\njySpRW2HfWY+BLx7gqeObXfbkqTO8AraHjA0NFR1CT3DYzHKYzHKY9G+tvrZt7XjiKxq35LUryKC\n7PYXtJKk/mDYS1INGPaSVAOGvSTVgGEvSTVg2EtSDRj2klQDhr0k1YBhL0k1YNhLUg0Y9pJUAz0d\n9rfdBjHtESAkqb9s2gQzOjHg/CR6Ouw/+MGqK5Ck8i1dClu2wPPPl7ePng57SaqDO+8sfj73XHn7\n6IuwdyRkSYPsoYeKn3fdVd4++iLsX3ml6gokqb/1Rdj7Ja2kQfb97xc/yzyN09N3qhoJ+Zdfhl/7\ntS4UJUkVGMm6uXPhqaemWneA71Q11S8vSYOgzLZ3X4T98uVVVyBJ/a0vwl6S6qD2LXu7Xkqqg1/8\norxt90XYS5La0xdhf889VVcgSf2tL8L+ttuqrkCS+ltfhL3n7CWpPR0J+4jYMSIeiIibG/NzImJp\nRKyMiO9FxOx2tn///Z2oUpJ6z8svd2c/nWrZfw5YDoy0wc8Dlmbm/sAdjXlJ0jhveEN39tN22EfE\nXsAi4NvAyCW8JwFXNKavAH6/3f1IklrXiZb93wB/Cmwds2z3zFzfmF4P7N6B/UiSWtRW2EfE7wG/\nyMwHGG3Vb6Mx2lnbX7Fu3Tr1OpKkibV718P3ACdFxCJgZ2CXiLgSWB8RczNzXUS8BZjwurAlS5a8\nNj00NMTQ0NBr8+OHNX70UZg/v81qJanPDA8PMzw83PZ2OjbEcUS8D/ifmXliRPwV8KvM/FJEnAfM\nzszzxq0/6RDH48N+xQo44ICOlCpJPWN81q1dC3vuOdn6vTHE8Uh6XwQcFxErgQ805iVJUyhrfJyO\nhX1m3pWZJzWmn8nMYzNz/8w8PjPbvv+KwxxLqoPVq8vZbl9cQQvwwx9WXYEkle+mm8rZbt+EvUMm\nSKqDsrLOsJekHvLKK+Vst2/CfsuWqiuQpPJdc0052+1Y18tp73iaXS/B1r2kwTPdrOuVrpeSpB5k\n2EtSDRj2klQDhr0k1YBhL0k1YNhLUg0Y9pJUA4a9JNWAYS9JNWDYS1INGPaSVAOGvSTVgGEvSTVg\n2EtSDRj2klQDhr0k1UBPhr03KZGkzurJsN+8ueoKJGmw9GTYv/BC1RVI0mDpybB/6qmqK5CkwdKT\nYf+zn1VdgSQNlp4M+4cfrroCSSpfNzuj9GTY33df1RVIUvlefrl7++rJsP/BD6quQJLKt2xZ9/YV\n2cbniIjYG/i/wJuBBL6VmX8bEXOAfwTeCqwGPpyZz417bW5v3xET78/+95IGSStZFxFk5nZeOcnr\n2gz7ucDczHwwIn4d+BHw+8AngV9m5l9FxBeAN2XmeeNea9hLqrVuhn1bp3Eyc11mPtiYfhFYAewJ\nnARc0VjtCop/AG3bsKETW5Gk+unYOfuIeBtwKHAPsHtmrm88tR7YvRP7sGUvSa2Z0YmNNE7hXA98\nLjP/K8Z8NsnMjIgJY3rJkiWvTQ8NDTE0NDTpflauhMMO60DBktQnhoeHGR4ebns7bZ2zB4iInYB/\nAW7NzK82lj0CDGXmuoh4C/CDzDxg3Oumfc7+2mvhQx9qq1xJ6hl9c84+iib8pcDykaBvuAk4ozF9\nBnBDO/sZ4WkcSWpNu+fsfxv4BPD+iHig8VgIXAQcFxErgQ805tt2992d2Iok1U/bp3Fa3nELp3He\n8Q5YtarEoiSpi7aXdc8/D7vssr3XVHAap9u2bKm6Akkq389/3vlt9lXYb9pUdQWSVL5nn+38Nvsq\n7J98suoKJKl8X/lK57fZV2EvSXXw0kud36ZhL0k9Zt26zm+zZ8I+s/hmevnyqiuRpGo98ADMm7f9\n3jqt6Jmul2vWwD77FN0rH3ts+6/zwipJg6KZMB+feX3f9fKGxjW2kwW9JKk1PdOyb/bjii17SYOi\nmdzbtAlmjBmysu9b9pKk1+vUxaSGvST1sE4NEVNp2H/pS539tlmSBs0TTxQ/I+CFF1rfTqVhf955\nU68jSXW2Zg08/HAxfdpprW+nI3eqkiSV47LLRk/lLF3a+nZ64pz9iy9WXYEk9a6LL25/Gz0R9qtX\nN7/uZZfBRz9aWimSVLpvfxs2b25u3Xvv7cw+K+1nD8W+r78eTj11eq+3v72kfrRxI8yaBXPmwDPP\ntLKFPu5nf911VVcgSd0xMnxxa0Hfup5o2bfClr2kftR+d/M+btm3wv75ktS8vg17SVLzDHtJ6pIq\n76Nt2EtSl3z2s9Xtu6+voN2ypRj6c/wQoJLUS269FT70Idiwoboa+rY3zlizZlV7ECVpMp3tUNJa\nb5yBCHuwK6ak3tT5noM163o53o9+VBxUu2RK6gU77ABXXVV1FaMGpmW/887w8svF9AMPwLve1bFN\nS9K0lNvo7LGWfUQsjIhHIuLRiPhCWfsZMRL0AIceOtrKdygGSd0wkjk33lh1JRMrpWUfETsC/x84\nFngSuA/4WGauGLNOR1v2zdq6tfh4BZ7nlzR9EcVIlAsWVFZBT7XsFwCrMnN1Zm4CrgFOLmlf07LD\nmN945D/xyOPMM0enN2wYnV6/HlauLKY3bixee/XV8Oyz1fwOklqzatXo+/ZrX4OjjiqmR97rt9yy\nbSYcdNDo9KJFo6dnqgv61pXVsj8NOCEzz27MfwI4IjPPGbNOJS17dcenPgXf/ObU6914I6xbV6wP\n8Cd/An//91O/7vTT4cori+lZs4rTdYsWFfPXXw///M/wD/8A3/gG7LHH6Ke4zNHptWvh3HOL6zSu\nvBJ23XX0jb14cfGlP8Bf/iV861vFtsY3ECLg+OOL9T7+cfjOd4rpk0/e9uP8UUfBsmUwfz7stx/8\n679O/Tsedhj8+MfF9MyZ8OqrxfTf/V21F+eoaj3U9TIiTgUWTh32i8e8aqjxUJ2ceCI8/vjoPTZb\nddxxo7dsO/FEuPnmYvrgg+Gtby2mR1plIz9XrCg+sY343d8d/Wdw662jy/fbr2gRHnNMcfHeyDqZ\nxWnBO+9sr3ZpcsONx4i/aCnsycyOP4AjgdvGzJ8PfGHcOrnt26Y/H88+Ozp9++2j0+eem/nqq8X0\nKadkbt2aeeGFmeeck/nSS5mbN2fefXfm+vWZmcXzGzbka7ZuLR5Snb30UuaWLcX0xRePvl/OPjvz\n1FMzX3wx89Zbi/fZDTdkvvLK6Htw48bR6S9/edv37ec/Pzp94YXV58j0HmS2kMtltexnUHxBewzw\nc+BeeuQL2ueeg9mzi+l//3c48siitbZsWTE9cjh3GJgrECRN16ZNsNNOxfTIJ8HM0el3vrP9T6Ot\na+00TikjymTm5oj4LHA7sCNw6dig77bM4iP6rrsWj/H/38bOe2GWpJGgh23zYXx2RMANNxTf0Yzk\nxtlnwyWXlF/jdA3MRVVjnXVWcbAXLIDh4eILPEkq07HHwh13bPsJoBw99AVtUzvucNg/80xxA18o\nRsP0NIykquywQxH6t98OJ5zQ6a3XPOwr+jUkaVIOhNYBl11WdQWSNLnHHqu6gkJft+zHnhuzZS+p\nV43k1B13FNdrtLm13umN002GvKReN5JTVeZVX5/GkaR+UmW3bsNekmrAsJekGjDsJakG+jbsDz+8\n6gokqXV//Mfd3V/fhf3IAWpmrHRJ6lXdvidB3/Wz37oVFi4sLkOWpH7z3e/Chz+87cia09PHwyWc\ncELz4W2/ekmDorWumH0c9g8+CO96V3OvM+wlDYpuhn2l5+wz4f774ZBDqqxCknrbM88Up3/aaexW\n2rIfu+9m/8PZspc0KFrJvYg+bNlLkiZ30EGd2Y5hL0k97JxzOrMdT+NIUkWayb0XX4Q3vnHsa/r8\nNI4hLkmvNzbo29EzYQ/w6U/D5ZdXXYUk9YavfrVz2+qZ0zjbPjfxa+bNg5UrSyxKkrpostM424vm\nvj+N04xOfZyRpLrpq7CXpDo45ZTOb9Owl6Qec8ABnd+mYS9JPaaMe9X2VdifdVbVFUhSf+qrsN97\n76orkKT+1HLYR8SXI2JFRDwUEf8UEbuOee78iHg0Ih6JiOM7U6ok1cN++3V+m+207L8HHJSZhwAr\ngfMBIuJA4CPAgcBC4OsR0VefICSpSkcf3fltthzCmbk0M7c2Zu8B9mpMnwxcnZmbMnM1sApY0FaV\nkqS2dKrF/UfALY3pPYC1Y55bC+zZiZ3ss08ntiJJ9TNjsicjYikwd4Kn/iwzb26scwHwamZeNcmm\nJrzwd8mSJa9NDw0NMTQ0NGmx8+dP+rQkDZzh4WGGh4fb3k5bY+NExB8CZwPHZObLjWXnAWTmRY35\n24DFmXnPuNdOe2ycl16CN7yh5XIlqadsL+vWrIG99pr4ua6PjRMRC4E/BU4eCfqGm4CPRsTMiNgX\nmAfc2+p+JKluthf07Zj0NM4UvgbMBJZG8e9pWWZ+JjOXR8S1wHJgM/CZ7TbhJUld0VdDHG/cCDvv\nXGJRktRF28u6yWK5FkMcG/SS1Jq+CntJUmsMe0mqAcNekmrAsJekGujJsJ/iQlpJ0jT1ZNiffnrV\nFUhS+WbN6t6+ejLs3//+qiuQpPJ97GPd21dPhr2jW0qqg0WLurevngz7HXesugJJKt+hh3ZvXz05\nXELx/OuXOcKOpEGyaRPMnPn65bUfLkGSBsmMdoainCbDXpIqsr2B0Mpg2EtSDzn11HK2a9hLUg+Z\nM6ec7Rr2klQDhr0k1YBhL0k1YNhLUg30TdhfemnVFUhS+cq6/WrfXEHr1bOSBtH4rLv/fjj88MnW\n9wpaSep7ZV1Va9hLUg0Y9pJUA4a9JNWAYS9JNWDYS1IP2XXXcrZr10tJqtB0s66yrpcRcW5EbI2I\nOWOWnR8Rj0bEIxFxfLv7kCS1p60enRGxN3Ac8MSYZQcCHwEOBPYEvh8R+2fm1nb2JUlqXbst+78G\nPj9u2cnA1Zm5KTNXA6uABW3uR5LUhpbDPiJOBtZm5k/GPbUHsHbM/FqKFr4kaZx587qzn0lP40TE\nUmDuBE9dAJwPjD0fP9kXBn69KkkTOOkk+MpXyt/PpGGfmcdNtDwi3gnsCzwUxVfJewE/iogjgCeB\nvcesvldj2essWbLktemhoSGGhoZem//iF+GCC5r4DSSpjx1zzORhPzw8zPDwcNv76UjXy4h4HDg8\nM59pfEF7FcV5+j2B7wP7je9nOVXXy2XL4D3vGZ2366WkQfT00/DmN4/Ol9X1slPjq71WXmYuj4hr\ngeXAZuAzk6a6JNXYb/xGd/bTsxdV2bKXVBdjL6zq2YuqJEm9r2fDft99q65AkgZHz57GKdYZnfY0\njqRB5WkcSVJHGPaSVAOGvSTVgGEvSTVg2EtSDRj2klQDPR32l1xSdQWS1D3r1pW37Z7uZ1+sV/y0\nn72kQTWdnBvYfvZHHw333FN1FZJUrve+t9zt93zLXpI0amBb9pKk9hn2klQDhr0k1YBhL0k1YNhL\nUg0Y9pJUA4a9JNWAYS9JNWDYS1INGPaSVAOGvSTVgGEvSTVg2EtSDRj2klQDhr0k1UBbYR8R50TE\nioh4OCK+NGb5+RHxaEQ8EhHHt1+mJKkdLYd9RLwfOAk4ODPfCVzcWH4g8BHgQGAh8PWI8BPEJIaH\nh6suoWd4LEZ5LEZ5LNrXTgh/GvjfmbkJIDOfbiw/Gbg6Mzdl5mpgFbCgrSoHnH/IozwWozwWozwW\n7Wsn7OcBvxMRd0fEcET8VmP5HsDaMeutBfZsYz+SpDbNmOzJiFgKzJ3gqQsar31TZh4ZEe8GrgXe\nvp1NebNZSapQyzccj4hbgYsy867G/CrgSOAsgMy8qLH8NmBxZt4z7vX+A5CkFrRyw/FJW/ZTuAH4\nAHBXROwPzMzMX0bETcBVEfHXFKdv5gH3dqJYSVJr2gn7y4DLIuKnwKvAHwBk5vKIuBZYDmwGPpOt\nfnyQJHVEy6dxJEn9o/T+7xGxsHFx1aMR8YXtrPO3jecfiohDy66pKlMdi4j4eOMY/CQi/iMiDq6i\nzm5o5u+isd67I2JzRJzSzfq6qcn3yFBEPNC4gHG4yyV2TRPvkd0i4raIeLBxLP6wgjJLFxGXRcT6\nxpmT7a0zvdzMzNIewI4U/ezfBuwEPAj85rh1FgG3NKaPAO4us6aqHk0ei6OAXRvTC+t8LMasdyfw\nL8CpVddd4d/FbOA/gb0a87tVXXeFx2IJxfU9ALsBvwJmVF17CcfivcChwE+38/y0c7Pslv0CYFVm\nrs7i4qtrKC66Gusk4AqALHrszI6I3UuuqwpTHovMXJaZzzdm7wH26nKN3dLM3wXAOcB1wNMTPDco\nmjkW/wO4PjPXAmTmL7tcY7c0cyyeAnZpTO8C/CozN3exxq7IzP8HPDvJKtPOzbLDfk9gzZj5iS6w\nmmidQQy5Zo7FWGcCt5RaUXWmPBYRsSfFG/0bjUWD+uVSM38X84A5EfGDiLg/Ik7vWnXd1cyxuAQ4\nKCJ+DjwEfK5LtfWaaedmO71xmtHsG3R8N8xBfGM3/Ts1xh36I+C3yyunUs0ci68C52VmRkTw+r+R\nQdHMsdgJOAw4BpgFLIuIuzPz0VIr675mjsWfAQ9m5lBEvANYGhGHZOZ/lVxbL5pWbpYd9k8Ce4+Z\n35tth1KYaJ29GssGTTPHgsaXspcACzNzso9x/ayZY3E4cE2R8+wGfDAiNmXmTd0psWuaORZrgF9m\n5kZgY0T8G3AIMGhh38yxeA/wRYDMfCwiHgfmA/d3pcLeMe3cLPs0zv3AvIh4W0TMpBgNc/yb9SYa\nffQj4kjgucxcX3JdVZjyWETEPsA/AZ/IzFUV1NgtUx6LzHx7Zu6bmftSnLf/9AAGPTT3HrkRODoi\ndoyIWRRfyC3vcp3d0MyxeAQ4FqBxjno+8LOuVtkbpp2bpbbsM3NzRHwWuJ3im/ZLM3NFRHyq8fw3\nM/OWiFjUGG5hA/DJMmuqSjPHAvhz4E3ANxot2k2ZOXAjhjZ5LGqhyffII41hR34CbAUuycyBC/sm\n/y7+F3B5RDxE0Vj9fGY+U1nRJYmIq4H3AbtFxBpgMcXpvJZz04uqJKkGvKmIJNWAYS9JNWDYS1IN\nGPaSVAOGvSR1WDMDmY1Z93ci4scRsSkiTh333BkRsbLx+IN2ajLsJanzLqcYzLAZTwBnAFeNXRgR\ncyi6Yy9oPBZHxOxWCzLsJanDJhrILCLeERG3NsY3+reImN9Y94nM/CnFNRRjnQB8LzOfy8zngKU0\n/w/kdcoeLkGSVPgW8KnMXBURRwBfpxjvaHv2YNvhIqYaPHFShr0klSwifp3ifhXfbVwdDzCzmzUY\n9pJUvh0oxq+Z6o5SY4c0eBIYGjO/N8XNfFouQJJUosx8AXg8Ik4DiML4246OH8r7duD4iJgdEW8C\njmssa4lhL0kd1hjI7IfA/IhYExGfBD4OnBkRDwIPU9xtauQ+y2uA04BvjnTXbAxxfiFwH3Av8BeN\nL2pbq8mB0CRp8Nmyl6QaMOwlqQYMe0mqAcNekmrAsJekGjDsJakGDHtJqgHDXpJq4L8B7Z/2vjyG\nuOcAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "#Data Rate = 100MHz\n", "d3 = np.zeros(100000, dtype=float)\n", "d3[40000:40100] = np.ones(100, dtype=float)\n", "x3 = d3 * c\n", "\n", "X3 = np.fft.fft(x3)\n", "\n", "plot(k, X3)\n", "#xlim(0,600000000)\n" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": true }, "outputs": [], "source": [] } ], "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.4.3" } }, "nbformat": 4, "nbformat_minor": 0 }