extraer_vel_del_photogate_sin_graficos_con_funcion.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# IMPORTACIÓN DE LAS BIBLIOTECAS QUE NECESITAMOS USAR\n",
    "import numpy as np\n",
    "import pandas as pd\n",
    "import matplotlib.pyplot as plt\n",
    "import matplotlib.mlab as mlab"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# Me genero unos datos de juguete\n",
    "from scipy import signal\n",
    "\n",
    "fs = 1000 #Hz\n",
    "tmax = 3 # Tiempo (segundos) de registro\n",
    "t = np.linspace(1/fs,tmax,tmax*fs)\n",
    "\n",
    "f = 13. # Hz. Acuerdense del punto para que no sean enteros!!\n",
    "t0 = 1\n",
    "ciclos = 13\n",
    "dur = ciclos/f\n",
    "A = 5 # Amplitud en Volts\n",
    "\n",
    "xo = A*(signal.square(2 * np.pi * f * t)+1)/2\n",
    "xo[t<t0] = 0\n",
    "xo[t>t0+dur] = 0\n",
    "\n",
    "ruido = np.random.normal(0, 0.01*A, len(xo)) # Le sumo ruido con el 1% de la amplitud\n",
    "x = xo + ruido"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# OPCION 2: GENERO LA FUNCION (como hice para el R^2), y despues lo calculo\n",
    "def calculo_velocidad(t,x,distancia,UMBRAL):\n",
    "    y = (x > UMBRAL)*1.\n",
    "    dy = y[1:len(y)] - y[0:(len(y)-1)]\n",
    "    ti = []\n",
    "    for item in range(0,len(dy)-1):\n",
    "        if (dy[item]==1.):\n",
    "            ti.append(t[item])\n",
    "    ti = np.asarray(ti)\n",
    "    #items = [item for item in range(0,len(dy)-1) if (dy[item]==1.)]\n",
    "    #ti = t[items]\n",
    "\n",
    "    dt = ti[1:len(ti)] - ti[0:(len(ti)-1)]\n",
    "    vel = distancia / dt # m/s\n",
    "    tprom = (ti[1:len(ti)] + ti[0:(len(ti)-1)])/2\n",
    "\n",
    "    output = {'tprom':tprom,'vel':vel}\n",
    "    return output "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXEAAAEACAYAAABF+UbAAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAADNpJREFUeJzt3X+M5PVdx/HnCw7+UBEDRoxQT6WQBhCIP+gltnTamvAj\npmcajW1NaVFjoyJN/MM2RnPXpInWxKQiKjlzYmlC0LSJHJWmmJRNg72Ss+XKz0uP2vDjWq9pAe+K\nNjnI2z9mgGXZ3Zk7ZnfmzTwfyeRmdj633/d+2Xne8Nmd3VQVkqSeTpr1AJKkE2fEJakxIy5JjRlx\nSWrMiEtSY0ZckhobG/Ek5yT5fJKHkjyQ5Po11t2Q5GCS/Ukunf6okqSVtkyw5jngj6pqf5IfAr6c\n5K6qOvDCgiRXAedW1XlJ3gjcBGzbmJElSS8Y+0y8qv67qvaPrn8PeAQ4e8Wy7cAtozX3AqcnOWvK\ns0qSVjiuPfEkPwVcCty74q6zgSeW3T7EK0MvSZqyiSM+2kr5FPDB0TNySdKMTbInTpItDAP+yaq6\nfZUlh4DXLbt9zuhtK9+PP6hFkk5AVWW1t0/6TPwfgYer6q/XuH8PcA1Akm3AM1V1eI1BvExw2bFj\nx8xn6HDxPHmeFuFcrWfsM/EkvwT8JvBAkvuAAv4E2Dpscu2qqjuTXJ3kUeBZ4Npx71eS9OqNjXhV\n/Qdw8gTrrpvKRJKkifmKzTk1GAxmPUILnqfJeJ4m1+1cZdx+y1QPltRmHk+SXguSUK/yC5uSpDlk\nxCWpMSMuSY0ZcUlqzIhLUmNGXJIaM+KS1JgRl6TGjLgkNWbEJakxIy5JjRlxSWrMiEtSY0Zckhoz\n4pLUmBGXpMaMuCQ1ZsQlqTEjLkmNGXFJasyIS1JjRlySGjPiktSYEZekxoy4JDVmxCWpMSMuSY0Z\ncUlqzIhLUmNGXJIaM+KS1JgRl6TGjLgkNWbEJakxIy5JjRlxSWrMiEtSY0Zckhoz4pLUmBGXpMaM\nuCQ1ZsQlqTEjLkmNGXFJasyIS1JjYyOeZHeSw0nuX+P+tyR5JslXRpc/nf6YkqTVbJlgzc3A3wC3\nrLPmC1X1jumMJEma1Nhn4lV1D/D0mGWZzjiSpOMxrT3xbUnuS/JvSS6Y0vuUJI0xyXbKOF8GtlbV\n/ya5CvhX4Py1Fu/cufPF64PBgMFgMIURJOm1Y2lpiaWlpYnWpqrGL0q2AndU1cUTrP0G8PNV9dQq\n99Ukx5MkvSQJVbXqtvWk2ylhjX3vJGctu34Zw38YXhFwSdL0jd1OSXIrMADOTPI4sAM4Faiq2gX8\nWpLfA44B/wf8xsaNK0labqLtlKkdzO0USTpu09hOkSTNISMuSY0ZcUlqzIhLUmNGXJIaM+KS1JgR\nl6TGjLgkNWbEJakxIy5JjRlxSWrMiEtSY0Zckhoz4pLUmBGXpunoUdi7d/intAmMuDQtR4/Cm98M\nl18+/NOQaxPMf8Rn8cxmEY65CB/jZh/zwQfhoYfguefg4YeH1zfaa/2czuqYnT7Gqtq0y/Bwx+HI\nkapLLqnasmX455Ejx/f3T8QiHHMRPsZZHPOF451yiv8dOx9zDj/GUTtX7+pad2zE5bgj/sUvDj8o\nGD4w9u49vr9/IhbhmIvwMc7qmEeODI+zGQ/8RTmnPj7Wjfh8b6dcdBFceCGccgpccMHwusfsd7xF\nOuZpp8G2bcM/N9qinFMfH+ua/1+UfPTocG/xwgs354GxKMdchI9xVsfcTItyThf88bHeL0qe/4hL\n0oLzt91L0muUEZekxoy4JDVmxCWpMSMuSY0ZcUlqzIhLUmNGXJIaM+KS1JgRl6TGjLgkNWbEJakx\nIy5JjRlxSWrMiEtSY0Zckhoz4pLUmBGXpMaMuCQ1ZsQlqTEjLkmNGXFJasyIS1JjRlySGjPiktTY\n2Ign2Z3kcJL711lzQ5KDSfYnuXS6I0qS1jLJM/GbgSvWujPJVcC5VXUe8AHgpinNJkkaY2zEq+oe\n4Ol1lmwHbhmtvRc4PclZ0xlPkrSeaeyJnw08sez2odHbJEkbbMtmH3Dnzp0vXh8MBgwGg80eQZLm\n2tLSEktLSxOtTVWNX5RsBe6oqotXue8m4O6q+ufR7QPAW6rq8Cpra5LjSZJekoSqymr3TbqdktFl\nNXuAa0YH2gY8s1rAJUnTN3Y7JcmtwAA4M8njwA7gVKCqaldV3Znk6iSPAs8C127kwJKkl0y0nTK1\ng7mdIknHbRrbKZKkOWTEJakxIy5JjRlxSWrMiEtSY0Zckhoz4pLUmBGXpMaMuCQ1ZsQlqTEjLkmN\nGXFJasyIS1JjRlySGjPiktSYEZekxoy4JDVmxCWpMSMuSY0ZcUlqzIhLUmNGXJIaM+KS1JgRl6TG\njLgkNWbEJakxIy5JjRlxSWrMiEtSY0Zckhoz4pLUmBGXpMaMuCQ1ZsQlqTEjLkmNGXFJasyIS1Jj\nRlySGjPiktSYEZekxoy4JDVmxCWpMSMuSY0ZcUlqzIhLUmNGXJIaM+KS1NhEEU9yZZIDSb6W5EOr\n3P++JN9O8pXR5bemP6okaaUt4xYkOQm4EXg78E1gX5Lbq+rAiqW3VdX1GzCjJGkNkzwTvww4WFWP\nVdUx4DZg+yrrMtXJJEljTRLxs4Enlt1+cvS2ld6ZZH+Sf0lyzlSmkySta+x2yoT2ALdW1bEkvwt8\nguH2yyvs3LnzxeuDwYDBYDClESTptWFpaYmlpaWJ1qaq1l+QbAN2VtWVo9sfBqqqPrbG+pOAp6rq\nR1a5r8YdT5L0ckmoqlW3rCfZTtkHvD7J1iSnAu9i+Mx7+QF+fNnN7cDDJzqsJGlyY7dTqur5JNcB\ndzGM/u6qeiTJR4B9VfUZ4Pok7wCOAU8B79/AmSVJI2O3U6Z6MLdTJOm4vdrtFEnSnDLiktSYEZek\nxoy4JDVmxCWpMSMuSY0ZcUlqzIhLUmNGXJIaM+KS1JgRl6TGjLgkNWbEJakxIy5JjRlxSWrMiEtS\nY0Zckhoz4pLUmBGXpMaMuCQ1ZsQlqTEjLkmNGXFJasyIS1JjRlySGjPiktSYEZekxoy4JDVmxCWp\nMSMuSY0ZcUlqzIhLUmNGXJIaM+KS1JgRl6TGjLgkNWbEJakxIy5JjRlxSWrMiEtSY0Zckhoz4pLU\nmBGXpMaMuCQ1ZsQlqTEjLkmNTRTxJFcmOZDka0k+tMr9pya5LcnBJHuT/OT0R5UkrTQ24klOAm4E\nrgAuBN6d5A0rlv028FRVnQd8HPjLaQ+6aJaWlmY9Qguep8l4nibX7VxN8kz8MuBgVT1WVceA24Dt\nK9ZsBz4xuv4p4O3TG3ExdftEmhXP02Q8T5Prdq4mifjZwBPLbj85etuqa6rqeeCZJGdMZUJJ0po2\n6gub2aD3K0laJlW1/oJkG7Czqq4c3f4wUFX1sWVrPjtac2+Sk4FvVdWPrfK+1j+YJGlVVbXqk+Mt\nE/zdfcDrk2wFvgW8C3j3ijV3AO8D7gV+Hfj88QwhSToxYyNeVc8nuQ64i+H2y+6qeiTJR4B9VfUZ\nYDfwySQHge8yDL0kaYON3U6RJM2vqX9hM8nuJIeT3L/OmhtGLwzan+TSac/QxbhzleQ9Sb46utyT\n5Gc3e8Z5MMnn1GjdLyY5luSdmzXbvJnw8TdIcl+SB5PcvZnzzYsJHns/nGTPqFEPJHn/Jo84sY34\n7pSbGb4waFVJrgLOHb0w6APATRswQxfrnivgv4DLq+oS4KPAP2zKVPNn3Hl64UVpfwF8blMmml/j\nHn+nA38L/EpVXcTwa1iLaNzn1B8AD1XVpcBbgb9KMsnXEDfd1CNeVfcAT6+zZDtwy2jtvcDpSc6a\n9hwdjDtXVfWlqvqf0c0v8crvz18IE3xOAfwhwxeafXvjJ5pfE5yr9wCfrqpDo/Xf2ZTB5swE56mA\n00bXTwO+W1XPbfhgJ2AWPwBr5YuHDrGgcTpOvwN8dtZDzKMkPwH8alX9Pb5GYZzzgTOS3J1kX5L3\nznqgOXUjcEGSbwJfBT4443nWNJf/e6CXS/JW4FrgTbOeZU59HFj+g9kM+dq2AD8HvA34QWBvkr1V\n9ehsx5o7VwD3VdXbkpwL/HuSi6vqe7MebKVZRPwQ8Lplt88ZvU2rSHIxsAu4sqrGbSksql8AbksS\n4EeBq5Icq6o9M55rHj0JfKeqvg98P8kXgEsAI/5y1wJ/DlBVX0/yDeANwH/OdKpVbOTL7td6NrQH\nuAZefDXoM1V1eIPm6GDNczX6kb6fBt5bVV/f1Knmz5rnqap+ZnT5aYb74r+/4AFf7/F3O/CmJCcn\n+QHgjcAjmzbZfFnvPD0G/DLA6Gt25zP8RoO5M/Vn4kluBQbAmUkeB3YApzJ8qf6uqrozydVJHgWe\nZfgv3kIad66APwPOAP5u9CzzWFVdNqt5Z2WC87TcQr/wYYLH34EknwPuB54HdlXVwzMbeEYm+Jz6\nKPBPy74F8Y+r6qmZDDuGL/aRpMb89WyS1JgRl6TGjLgkNWbEJakxIy5JjRlxSWrMiEtSY0Zckhr7\nf2Vl+PZyuhihAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7ffb626dbe90>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "output = calculo_velocidad(t,x,distancia,UMBRAL)\n",
    "\n",
    "plt.plot(output['tprom'],output['vel'],'r.')\n",
    "plt.ylim(0,2)\n",
    "plt.show()"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 2",
   "language": "python",
   "name": "python2"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 2
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython2",
   "version": "2.7.12"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}