Permalink
Cannot retrieve contributors at this time
Name already in use
A tag already exists with the provided branch name. Many Git commands accept both tag and branch names, so creating this branch may cause unexpected behavior. Are you sure you want to create this branch?
compmech-project02/02_Analyze-data_project.ipynb
Go to fileThis commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
executable file
375 lines (375 sloc)
47.3 KB
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| { | |
| "cells": [ | |
| { | |
| "cell_type": "markdown", | |
| "metadata": {}, | |
| "source": [ | |
| "# Computational Mechanics Project #02 - Create specifications for a projectile robot\n", | |
| "\n", | |
| "On the first day of class, we threw $2\"\\times~2\"$ dampened paper (spitballs) at a target on the whiteboard. Now, we are going to analyze the accuracy of the class with some cool Python tools and design a robot that has the same accuracy and precision as the class, but we will have the robot move farther away from the target and use a simpler projectile i.e. a tennis ball so we don't need to worry about knuckle-ball physics. \n", | |
| "\n", | |
| "The goal of this project is to determine the precision of necessary components for a robot that can reproduce the class throwing distibution. We have generated pseudo random numbers using `numpy.random`, but the class target practice is an example of truly random distributions. If we repeated the exercise, there is a vanishingly small probability that we would hit the same points on the target, and there are no deterministic models that could take into account all of the factors that affected each hit on the board. \n", | |
| "\n", | |
| "<img src=\"../images/robot_design.png\" style=\"height: 250px;\"/>\n" | |
| ] | |
| }, | |
| { | |
| "cell_type": "markdown", | |
| "metadata": {}, | |
| "source": [ | |
| "Now, we ask ourselves some questions:\n", | |
| "\n", | |
| "1. How do we quantify the class accuracy and precision?\n", | |
| "\n", | |
| "2. If we design a robot, what design components can we control?\n", | |
| "\n", | |
| "3. How can we relate the controlled components to the class accuracy, and specify the component precision?\n", | |
| "\n", | |
| "The first question, we have some experience from our work in [02_Seeing_Stats](../notebooks/02_Seeing_Stats.ipynb). We can define the mean, standard deviation, measure the first, second, and third quartiles, etc. \n", | |
| "\n", | |
| "The second question is a physical question. We cannot control the placement of the robot or the target those are chosen for us. We cannot control temperature, mechanical vibrations, etc. We *can* control the desired initial velocity. The initial velocity will have some speed and direction, and both will be subject to random noise. Once the speed and direction are set, the location on the target is determined by kinematic equations for an object in freefall, as such\n", | |
| "\n", | |
| "$x_{impact} = \\frac{v_x}{v_y}d + x(0)~~~~~~~~~~~~~~~~~~~~(1.a)$\n", | |
| "\n", | |
| "$z_{impact} = d\\left(\\frac{v_z(0)}{v_y}-\\frac{g}{2v_y^2}d\\right)+ z(0)~~~~~(1.b)$.\n", | |
| "\n", | |
| "Where the location of impact is at a $y$-distance of $d$ at a point on the target with coordinates $(x_{impact},~z_{impact})$, and the initial velocity is $\\bar{v}=v_x\\hat{i}+v_y\\hat{j}+v_z(0)\\hat{k}$, the object is released at an initial location $\\bar{r}(0)=x(0)\\hat{i}+0\\hat{j}+z(0)\\hat{k}$, and the only acceleration is due to gravity, $\\bar{a}=-g\\hat{k}$. Equation (1) becomes much easier to evaluate if we assume that $v_x=0$, resulting in an evalution of the accuracy of the height of the impact, $z_{impact}$, as such\n", | |
| "\n", | |
| "$x_{impact} = x(0)~~~~~~~~~~~~~~~~~~~~(2.a)$\n", | |
| "\n", | |
| "$z_{impact} = \\frac{d}{\\cos{\\theta}}\\left(\\sin{\\theta}-\\frac{g}{2v_0^2\\cos{\\theta}}d\\right)+ z(0)~~~~~(2.b)$.\n", | |
| "\n", | |
| "Where $\\theta$ is the angle of the initial velocity and $v_0$ is the initial speed. Equation (2) restricts the analysis to height accuracy. You can incorporate the 2D impact analysis if you finish the 1D analysis. " | |
| ] | |
| }, | |
| { | |
| "cell_type": "markdown", | |
| "metadata": {}, | |
| "source": [ | |
| "The third question, is how we can relate equation (2) to the measured points of impact? For this, we can use Monte Carlo methods *(There are other methods, but Monte Carlo is one of the most straight-forward)*. Our Monte Carlo approach is as such, if we have a desired initial speed, $v_0$, and desired angle, $\\theta$, we can propagate the uncertainty of our actual speeds and angles into the $z_{impact}$ locations. Then, we can choose distributions in speed and angles that match the distributions in $z_{impact}$ locations. Here are the steps:\n", | |
| "\n", | |
| "1. Generate random $\\theta_i$ and $v_{0~i}$ variables\n", | |
| "\n", | |
| "2. Plug into eqn 2 for random $z_{impact~i}$ locations\n", | |
| "\n", | |
| "3. Compare to our measured $z_{impact}$ location statistics\n", | |
| "\n", | |
| "4. Repeat 1-3 until the predicted uncertainty matches the desired uncertainty, we can use a number of comparison metrics:\n", | |
| " \n", | |
| " - standard deviation\n", | |
| " \n", | |
| " - first, second, and third quartiles\n", | |
| " \n", | |
| " - visually, with box plots and histograms" | |
| ] | |
| }, | |
| { | |
| "cell_type": "markdown", | |
| "metadata": {}, | |
| "source": [ | |
| "# Project Deliverables\n", | |
| "\n", | |
| "1. Statistical analysis of class accuracy and precision (x- and z-locations) data is in the csv file [../data/target_data.csv](../data/target_data.csv) _Note: if you want to see how I turned the images into data check out the jupyter notebook [process_target_practice](./process_target_practice.ipynb)\n", | |
| "\n", | |
| "2. A Monte Carlo model to generate impact heights based upon uncertainty in $\\theta_0$ and $v_0$. \n", | |
| "\n", | |
| "3. The precision required to recreate the class accuracy and precision with a robot. \n", | |
| "**You must show some validation of your work**\n", | |
| "\n", | |
| "4. [BONUS] Repeat 2-3 taking into account the variation in $x_{impact}$ due to misalignment. \n", | |
| "\n", | |
| "Given constants and constraints:\n", | |
| "\n", | |
| "- $d=$3 m, distance to target\n", | |
| "\n", | |
| "- $g=$9.81 m/s$^2$, acceleration due to gravity\n", | |
| "\n", | |
| "- $z(0)=$0.3 m, the initial height is 0.3 m above the bull's eye\n", | |
| "\n", | |
| "- 4 m/s$<v_0<$12 m/s, the initial velocity is always higher than 9 mph and less than 27 mph" | |
| ] | |
| }, | |
| { | |
| "cell_type": "code", | |
| "execution_count": 282, | |
| "metadata": {}, | |
| "outputs": [], | |
| "source": [ | |
| "import matplotlib.pyplot as plt\n", | |
| "import numpy as np\n", | |
| "import pandas as pd\n", | |
| "\n", | |
| "#Import rcParams to set font styles\n", | |
| "from matplotlib import rcParams\n", | |
| "\n", | |
| "#Set font style and size \n", | |
| "rcParams['font.family'] = 'sans'\n", | |
| "rcParams['font.size'] = 16\n", | |
| "rcParams['lines.linewidth'] = 3" | |
| ] | |
| }, | |
| { | |
| "cell_type": "markdown", | |
| "metadata": {}, | |
| "source": [ | |
| "# Section 1" | |
| ] | |
| }, | |
| { | |
| "cell_type": "code", | |
| "execution_count": 283, | |
| "metadata": {}, | |
| "outputs": [ | |
| { | |
| "name": "stdout", | |
| "output_type": "stream", | |
| "text": [ | |
| "The first quartile for x is -0.16854422326722185\n", | |
| "The second quartile for x is 0.028452517775170515\n", | |
| "The third quartile for x is 0.1670894539891312\n" | |
| ] | |
| }, | |
| { | |
| "data": { | |
| "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZUAAAD9CAYAAAB0i+q4AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAbfklEQVR4nO3df5QV5Z3n8fdH/J2NkyZKdLLyQ3E8IJpM9upqwmbQg4E9GjXjqhCZ6MYJGSeiZo67kWESjegJTrLqJJuZlawbJpAlCR6NP+KgBpETVEwalFlQERSMOiZiGn8A0ij57h9VbYrq2933dj/VzY/P65x7uu9Tz1P11O3u++mnnrpVigjMzMxS2GegO2BmZnsOh4qZmSXjUDEzs2QcKmZmloxDxczMktl3oDswkA499NAYPnz4QHfDzGy3snz58tci4rB6y/bqUBk+fDitra0D3Q0zs92KpBe6WubDX2ZmloxDxczMknGomJlZMg4VMzNLxqFiZmbJOFTMzCwZh4qZmSXjUDEzs2T26g8/mvUXSf2yHd8fyQaaQ8WsHzT7Zi/JAWG7JR/+MjOzZBwqZmaWjEPFzMyScaiYmVkyDhUzM0vGoWJmZsk4VMzMLBmHipmZJeNQMTOzZCoNFUlHSrpd0huS3pR0h6ShDbaNLh4fLdXbR9J0SRskbZO0UtK51eyRmZl1p7LLtEg6GHgIaAcuAgK4Hlgs6YSI2NLAauYAt5bKni09nwlcBcwAlgOTgAWSzoyI+3q/B2Zm1qwqr/31BeAo4NiIWAcg6V+BtcAXgZsaWMfLEbGsq4WShpAFyqyI+FZevFjSSGAW4FAxM+tHVR7+OgtY1hEoABGxHngEODvRNiYA+wPzSuXzgOMljUi0HTMza0CVoXIcsKpO+WpgdIPruFRSu6Stkh6S9J/qbKMdWFcqX51/bXQ7ZmaWQJWhMhjYVKe8DWhpoP084K+B8cBU4IPAQ5LGlbbxenS+RnhbYflOJE2V1CqpdePGjQ10w8zMGlX1/VTq3RCiobsVRcRfFJ7+QtJdZCOf64GxhXU1tY2ImA3MBqjVar5hhZlZQlWOVDZRZ6RANkqpN4LpVkS8BfwMOLFQ3Aa0qPNt9VoKy83MrJ9UGSqryeY8ykYDT/VyneWRyWrgAODoOtugD9sxM7NeqDJU7gZOlnRUR4Gk4cAn8mVNkXQIcAbweKF4IbAduLBUfQqwKj/bzMzM+kmVcyrfAy4D7pL0d2QjjJnAixQ+0ChpGPAccF1EXJeXXQUcCywG/g0YRvZ5lMMpBEhEvCrpZmC6pLeAFcAFwGmkO23ZzMwaVFmoRMQWSacBNwNzyQ5dLQKujIjNhaoCBrHzqGkN8Jn88UfAm2Sfb7kkIn5Z2tQMYDNwBVnorAHOj4h7ku+UmZl1S53Pxt171Gq1aG1tHehumHUiib35b9N2bZKWR0St3jJfpdjMzJJxqJiZWTIOFTMzS8ahYmZmyThUzMwsGYeKmZkl41AxM7NkHCpmZpaMQ8XMzJJxqJiZWTIOFTMzS8ahYmZmyThUzMwsGYeKmZkl41AxM7NkHCpmZpaMQ8XMzJJxqJiZWTIOFTMzS8ahYmZmyThUzMwsGYeKmZkl41AxM7NkHCpmZpZMpaEi6UhJt0t6Q9Kbku6QNLSBdjVJsyU9I2mrpF9L+qGkEXXqbpAUdR7nVLNXZmbWlX2rWrGkg4GHgHbgIiCA64HFkk6IiC3dNJ8EHAd8G1gNfBj4KtAq6aMR8WKp/v3AtaWyNX3eCTMza0ploQJ8ATgKODYi1gFI+ldgLfBF4KZu2t4YERuLBZIeAdbn6/1aqf5rEbEsVcfNzKx3qjz8dRawrCNQACJiPfAIcHZ3DcuBkpe9AGwkG7WYmdkuqMpQOQ5YVad8NTC62ZVJGgUMAZ6us/jT+dxLu6Rlnk8xMxsYVYbKYGBTnfI2oKWZFUnaF/hfZCOV20qL7wGmAROAC4FtwJ2SpnSxrqmSWiW1btzYaUBkZmZ9UOWcCmST82XqxXr+J/Bx4IyI2CmoImLaTiuX7gSWAd8A5nXqUMRsYDZArVar1z8zM+ulKkcqm8hGK2Ut1B/B1CXpG8BU4PMR8UBP9SNiB7AA+PeSjmh0O2Zm1ndVjlRWk82rlI0GnmpkBZJmAFcDl0fE3Ca23TEa8kjEzKwfVTlSuRs4WdJRHQWShgOfyJd1S9LlZJ9rmRER32l0o/n8y3nAryPiN0322czM+qDKUPkesAG4S9LZks4C7gJeBG7tqCRpmKR3JX2tUDYJuAVYCDwk6eTCY3Sh3mRJP5L0OUmn5u0WA/8B+EqF+2ZmZnVUdvgrIrZIOg24GZhLdkhqEXBlRGwuVBUwiJ0DbmJePjF/FC0BxuXfryc7zfibZPM3W4FfARMj4v6U+2NmZj1TxN477VCr1aK1tXWgu2HWiST25r9N27VJWh4RtXrLfJViMzNLxqFiZmbJOFTMzCwZh4qZmSXjUDEzs2SqvvaX2R5n8ODBbNrU8JWGek3qzWXymtPS0kJbW1vl27G9h0PFrEmbNm3aY0737Y/gsr2LD3+ZmVkyDhUzM0vGoWJmZsk4VMzMLBmHipmZJeNQMTOzZBwqZmaWjEPFzMyScaiYmVkyDhUzM0vGoWJmZsk4VMzMLBmHipmZJeNQMTOzZBwqZmaWjEPFzMyScaiYmVkyDhUzM0um0lCRdKSk2yW9IelNSXdIGtpg2wMlfVPSK5LelvSYpE/WqbePpOmSNkjaJmmlpHPT742ZmfWksnvUSzoYeAhoBy4CArgeWCzphIjY0sMqbgPOAP4b8DzwJeB+SadExJOFejOBq4AZwHJgErBA0pkRcV/KfTIDiGsOgWv/aKC7kURcc8hAd8H2MIqIalYsXQHcBBwbEevyshHAWuC/R8RN3bT9CPAk8PmI+H5eti+wGlgTEWflZUOAF4FZEXFNof0i4LCIOKG7PtZqtWhtbe3DXtreSBJV/d30tz1pX6z/SFoeEbV6y6o8/HUWsKwjUAAiYj3wCHB2A23fAX5caPsu8CNggqQD8uIJwP7AvFL7ecDxeYiZmVk/qTJUjgNW1SlfDYxuoO36iNhap+3+wMhCvXZgXZ16NLAdMzNLqMpQGQxsqlPeBrT0oW3H8o6vr0fn8Xu53nskTZXUKql148aNPXTDzMyaUfUpxfUO1qqBdmqwbaP1/tChiNkRUYuI2mGHHdZAV8zMrFFVhsom6owUyEYp9UYhRW3dtO1Y3vG1RVI5RMr1zMysH1QZKqvJ5jzKRgNPNdB2RH5acrntdv4wh7IaOAA4uk49GtiOmZklVGWo3A2cLOmojgJJw4FP5Mt6arsfcF6h7b7ABcADEdGeFy8kC5kLS+2nAKvys83MzKyfVPbhR+B7wGXAXZL+jmzuYybZ50pu7agkaRjwHHBdRFwHEBFPSvoxcIuk/YD1wKXACAoBEhGvSroZmC7pLWAFWfCcRs+nLZuZWWKVhUpEbJF0GnAzMJds8nwRcGVEbC5UFTCIzqOm/wrcQPYp/A8AK4GJEbGiVG8GsBm4AjgcWAOcHxH3pN0jMzPrSWWfqN8d+BP11ht70qfQ96R9sf4zUJ+oNzOzvYxDxczMknGomJlZMg4VMzNLxqFiZmbJOFTMzCwZh4qZmSXjUDEzs2QcKmZmloxDxczMknGomJlZMg4VMzNLxqFiZmbJOFTMzCwZh4qZmSXjUDEzs2QcKmZmlkyV96g322NJGuguJNHS0jLQXbA9jEPFrEn9cftd3+bXdlc+/GVmZsk4VMzMLBmHipmZJeNQMTOzZBwqZmaWTGWhImkfSdMlbZC0TdJKSec20O4QSV+T9Kik30l6Pf/+nDp1r5UUdR4/rWavzMysO1WeUjwTuAqYASwHJgELJJ0ZEfd1024o8NfA9/N1/B6YDNwp6bKI+G6dNmOBHYXnbQn6b2ZmTaokVCQNIQuUWRHxrbx4saSRwCygu1BZDxwVEVsLZfdLOhL4ClAvVB6PiHcTdN3MzPqgqsNfE4D9gXml8nnA8ZJGdNUwIraUAqVDK/DH6bpoZmapVRUqxwHtwLpS+er86+herPOTwDNdLHtR0g5JL0i6UdJBvVi/mZn1UVVzKoOB16PzdSbaCssbJmkqcDIwpbRoHXA18AQQwKeALwMfA05vss9mZtZHDYWKpPHAgw1UXRIR4wCRvcl3WlXjXXtv2+OAbwNzI+KHxWURUT689qCkl4BbJI2PiJ/XWd9UYCrA0KFDm+2OmZl1o9GRyqPAqAbqdcyFtAEtklQarbQUlvdI0onA3cBDwCUN9nU+cAtwItApVCJiNjAboFar+Yp9ZmYJNRQq+cR5V/MZ9awGDgCOZud5lY65lKd6WoGk44H7gSeBcyPinSa2D/VHSmZmVqGqJuoXAtuBC0vlU4BVEbG+u8aSjiE73PY8cGZEvN3Etju2+XgTbczMLIFKJuoj4lVJNwPTJb0FrAAuAE4Dzi7WlbQIGBYRI/PnQ8gCZX/gGmB06YZIT0REe173CeAHwBqykcnpwDRgYUQsrmLfzMysa1V+on4GsBm4Ajic7I3//Ii4p1RvUKkfo4Fh+ff31lnvCGBD/v0a4DLgiHw9zwHXAX/f9+6bmVmztDffXa5Wq0Vra+tAd8OsE9/50XZlkpZHRK3eMl+l2MzMknGomJlZMg4VMzNLxqFiZmbJOFTMzCwZh4qZmSXjUDEzs2QcKmZmloxDxczMknGomJlZMg4VMzNLxqFiZmbJOFTMzCwZh4qZmSXjUDEzs2QcKmZmloxDxczMknGomJlZMg4VMzNLxqFiZmbJOFTMzCwZh4qZmSXjUDEzs2QcKmZmloxDxczMkqksVCTtI2m6pA2StklaKencBtvOkRR1HrfUqTtW0qOS3pb0G0k3SToo/R6ZmVlP9q1w3TOBq4AZwHJgErBA0pkRcV8D7TcCZ5XKXik+kXQC8CBwP3AmMAL4JvBh4II+9d7MzJpWSahIGkIWKLMi4lt58WJJI4FZQCOhsj0ilvVQ5+vAS8B5EfFOvu3twD9LujEiVvRuD8zMrDeqOvw1AdgfmFcqnwccL2lEXzcgaT9gIvCTjkDJ/QTYDpzd122YmVlzqgqV44B2YF2pfHX+dXQD6xgi6TVJ70p6VtJXJA0qLD8aOBBYVWwUEduA5xrchlm/kNTUozdtOtqZDaSq5lQGA69HRJTK2wrLu/Mk2TzMarLg+AzwDeAY4C9L69hUp31bV9uQNBWYCjB06NAeumGWRuc/BbM9U0OhImk82YR4T5ZExDhAQL2/oob+lYqI8lle90naDFyZz5WsLayrqe1ExGxgNkCtVvNfuplZQo2OVB4FRjVQb2v+tQ1okaTSaKWlsLxZ84ErgRqwlu5HPS384VCbmZn1k4ZCJSK2As80sd7VwAFk8x7FeZWOeY6nmlhXh/LI5DmyeZvjdqokHQgcBSzoxTbMzKwPqpqoX0h2BtaFpfIpwKqIWN+LdX6WLFB+BRAR2/PtnC+pGI7/hSzQ7u7FNszMrA8qmaiPiFcl3QxMl/QWsILsw4inUTrVV9IiYFhEjMyfDwPmAj8iG+UcQDZRfzFwa0Q8V2h+LfAY8BNJ3wWGk3348faIWF7FvpmZWdeq/ET9DGAzcAVwOLAGOD8i7inVG1Tqx1tk8yVfAT5ENjp5Grgc+Mdiw4h4UtIE4EbgZ8AbwA+Av029M2Zm1jPtzac61mq1aG1tHehumJntViQtj4havWW+SrGZmSXjUDEzs2QcKmZmloxDxczMknGomJlZMg4VMzNLxqFiZmbJOFTMzCwZh4qZmSXjUDEzs2QcKmZmloxDxczMknGomJlZMg4VMzNLxqFiZmbJOFTMzCwZh4qZmSXjUDEzs2QcKmZmloxDxczMknGomJlZMg4Vs13I/PnzGTNmDIMGDWLMmDHMnz9/oLtk1pR9B7oDZpaZP38+M2bM4LbbbmPs2LEsXbqUSy65BIDJkycPcO/MGqOIGOg+DJharRatra0D3Q0zAMaMGcN3vvMdTj311PfKFi9ezLRp01i1atUA9sxsZ5KWR0St7jKHikPFdg2DBg1i27Zt7Lfffu+VvfPOOxx44IHs2LFjAHtmtrPuQqWyORVJ+0iaLmmDpG2SVko6t4F2wyVFN49JhbrXdlHnp1Xtl1lVRo0axdKlS3cqW7p0KaNGjRqgHpk1r8o5lZnAVcAMYDkwCVgg6cyIuK+bdq8Ap9Qpvx4YCzxQZ9lYoPivXFuvemw2gGbMmMEll1zSaU7lhhtuGOiumTWsklCRNIQsUGZFxLfy4sWSRgKzgC5DJSLagWWl9R0MnATcExH1AuPxiHg3SefNBkjHZPy0adN4+umnGTVqFDfccIMn6W23UtVIZQKwPzCvVD4P+D+SRkTE+ibW9+fA+4F/TtQ/s13S5MmTHSK2W6tqTuU4oB1YVypfnX8d3eT6LgJeBRZ2sfxFSTskvSDpRkkHNbl+MzNLoKqRymDg9eh8allbYXlDJH0YOA34hzqHuNYBVwNPAAF8Cvgy8DHg9C7WNxWYCjB06NBGu2FmZg1oKFQkjQcebKDqkogYB4jsTb7Tqhrv2nv+gmxE1enQV0SUD689KOkl4BZJ4yPi53XazAZmQ3ZKcS/6Y2ZmXWh0pPIo0Mh5jVvzr21AiySVRistheWN+hzwZESsbLD+fOAW4ESgU6iYmVl1GgqViNgKPNPEelcDBwBHs/O8SsdcylONrETSiWRh9uUmtt3BoxAzs35W1ZzKQmA7cCHw9UL5FGBVE2d+XQS8C/zfJrZ9Yf718Z4qLl++/DVJLzSxbrP+cijw2kB3wqwLw7paUEmoRMSrkm4Gpkt6C1gBXEA24X52sa6kRcCwiBhZKt+P7AOT/xIRr9bbjqQngB8Aa8hGJqcD04CFEbG4gX4e1uy+mfUHSa1dXQbDbFdW5SfqZwCbgSuAw8ne+M+PiHtK9QZ10Y8zgQ/S/WdT1gCXAUfk63kOuA74+z713MzMemWvvqCk2a7KIxXbXfkmXWa7ptkD3QGz3vBIxczMkvFIxczMknGomJlZMg4V22tImiNpQ+H58PxGb0fVqbtB0pz+7F9XOm5GVyoLSdc2sY6P5utp5rp74/LtjCuUPSxpaTfNmtJdv5rdR9s1OFRsbzIT+Ezh+XDgGqBTqOT1ZvZDn3rrFOB/N1H/o2T72nCokH2+7JT8a1W661ez+2i7gCo/p2K2S4mI55qo+0SVfemriFjWc63ekTSI7CSeNyndMK8/VbmPVh2PVKzXJL1P0jOSfplfAaGj/FOSfi/pSz203yBpnqQvSFonaZukFZJOrVN3iqSVeZ3XJM2VdESpzmclPSFps6Q3JP0/SV8sLH/v8Fd+SKfjqgsP5oda3jvUU+/wl6STJP08X/8WSYsknVSqM0fSS5L+VNIvJG2VtFbSX/X8ikKh3TZJL0v6KnWu7l0+NCTpTyTdKenVvO2vJS2QtK+ki4Hv51XXFvZ1eGFdN0i6WtJ6skssHV/v8Fdhe2dLWiWpPf8dOL/O67ChTruHJT2cf99Iv64ttZ8o6TFJb+c/459KOrbONpZKGp//Pm3N+3pOvdfc0nKoWK9FxBZgMvAR8kNFym4l/QPg3oj4bgOr+TPgb8iuwDCJ7OZu/1J8o1B2D5y5wNNkdwG9muzuoksk/bu8zliyO4suAc4BzgO+B3ygi+2uADpC73KyQy1dHuqRdEK+7hbgYrKrZx+S9+EjpeqHkF2vbh7ZZYl+BfxTvbAsbeNQ4CGy635dlPdvIvD57trl7gU+DFxK9tpcTfZa7gP8DLg+r3deYV9fKbS/GDiD7DbgZwD/1s22RgLfBv4H2c9jHfCjnvavjkb69R5JE/M2m8ku+3QpMAZYquy+S0VHA/8A3JT38RXgdmW3NLcqRYQffvTpQXYV6d8D48kuJvoycGgD7TaQ/Vc8tFD2frJbI8zNnw8CfgssLrUdS3a9t8vz51cBbT1sbw6wofB8XL6O8V30bU7h+e3A68AHCmWH5H29o7SNAE4tlB1AdnHI2T3074Y6r8f78rZRqhvAtfn3h+bPz+pm3RfndUbWWRZkIXJQqbzj9RlXKHs4Lzu5UDaI7Crmv+jqtS61f7iJfl1beN4KrAX2LZSNAN4Bbipt4x3gmELZEGAH8LcD/feypz88UrEUbiELk3vJ7r75uYho9Aq7yyLi1x1PIuItsv9GT8mLjiV7Q/hhsVFELAVeIBvpQDYaaMkPp50pqasRSm99kmz09XqhD28Cdxf60GFrFC5oGhHtZG+GPd1q9BQ6vx5bgPL18sp+BzwPzMoPJR7T087UsTAi3m6w7otRmO+IiB3AAuAkSZW8p0h6H9kdXX8chTvARnbF80fo/DNYGxFrC/VeJbsluW/3WjGHivVZZP8KziX7j3xlRCxqovlvuyjrOJzRcVZQvUMiv+lYHhFLyA6hHAncCWzM5z9OaKIv3RncTR9aSmWb6tRrBw7sYRtH0PXr0aX89T+d7D/5bwDPSnpe0qU9bK+o7iGnLnTVx/2Bqq783UI2t9Tt70FBvRsBNvIzsD5yqFifSTqcbLSyAviIpCuaaP6hLspezr/veHM4vE69w8n+SwcgIm6PiD8jewP6DNmb9MJE/z23ddOHZu5k2p1X6Pr16FZEPB8RnyN7U/9TsrmZf5T0nxvcdjPXa+qqj9uBjfnzbWQhU/bBJrZTtImsjz3+HtjAcqhYn0gS2e0JtpP9t3wLcGMTI4STJR1ZWN/7ySaKH8uL1pD9FzyptN2Pk90oaEl5hRGxOSLuBW4lC5au3sja868HNdDPJcAZef+Kff10vT700mN0fj3el2+jIZF5kuzkB8gmsqG5fe3JkZJOLvRxENko8ZcR8fu8+AXgQ/nJBx31jiY7nFnUUL/yw4DLgfPy7XWscxjwcdL9DKyPHCrWV39DNkE/JSLayM46egqYL6mRN7DfAg9IuiA/5fMBssnpmfDe8fqvAePz+ZKJki4B7iCbp/g+gKTrJN0qaZKkT0r6LNlZXU9GxMY62wV4luzOop+X9AlJtWJolMwke+NbJOlcSX8O/Bw4mOwePincDGyh8+vR7VyHpBMkLZb0V/lptBPIAvVdshEL/OEW3l+SdEq+r/VGEo34LfBjSRdLOgO4C/gTsp9ThwVkI4sfSpog6cK8XnmurZl+fRU4BrhX0qclTQYeBN4gOxPNdgUDfaaAH7vvg+wwSztwfan8WLI3x3/qof0GstNu/5LsBmvtwBPAaXXqTgFW5nV+RzaHc0Rh+RnA/WSHkNqBF4HbgD8u1JlD6Ywk4Itkk9zvUjjTidLZX3nZfyQLks35/i0CTirVmQO8VKf/D1M466mb1+RjwC/IDh+9TPZG+nW6P/trCNlo8VlgK9nhuCXAhFKba/J17sjbDy+s6/o6fRlH/bO/lgJnAavy13oNcEGd9ufkdd7Of3afqvc69NCva0t1J5KN6N4mC5O7gGPrvNZLu/h9m1Mu9yPtw5e+twGTfzhuaURMGei+mFkaPvxlZmbJOFTMzCwZH/4yM7NkPFIxM7NkHCpmZpaMQ8XMzJJxqJiZWTIOFTMzS+b/AyNAqG+EZjOYAAAAAElFTkSuQmCC\n", | |
| "text/plain": [ | |
| "<Figure size 432x288 with 1 Axes>" | |
| ] | |
| }, | |
| "metadata": { | |
| "needs_background": "light" | |
| }, | |
| "output_type": "display_data" | |
| } | |
| ], | |
| "source": [ | |
| "data = pd.read_csv(\"../data/target_data.csv\")\n", | |
| "x_position = data[' x position (m)']\n", | |
| "x_mean = np.mean(x_position)\n", | |
| "x_std = np.std(x_position)\n", | |
| "Q1_x = np.percentile(x_position, q=25)\n", | |
| "Q2_x = np.percentile(x_position, q=50)\n", | |
| "Q3_x = np.percentile(x_position, q=75)\n", | |
| "plt.boxplot(x_position, labels=['x position distribution']);\n", | |
| "print('The first quartile for x is {}'.format(Q1_x))\n", | |
| "print('The second quartile for x is {}'.format(Q2_x))\n", | |
| "print('The third quartile for x is {}'.format(Q3_x))" | |
| ] | |
| }, | |
| { | |
| "cell_type": "code", | |
| "execution_count": 284, | |
| "metadata": { | |
| "scrolled": true | |
| }, | |
| "outputs": [ | |
| { | |
| "name": "stdout", | |
| "output_type": "stream", | |
| "text": [ | |
| "The first quartile for z is -0.21785495681409592\n", | |
| "The second quartile for z is -0.070851135369002\n", | |
| "The third quartile for z is 0.13222390283592378\n" | |
| ] | |
| }, | |
| { | |
| "data": { | |
| "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZUAAAD9CAYAAAB0i+q4AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAYiElEQVR4nO3df5SkVX3n8fdHEPwVkyFKMAYYEEKAVTmmjajEIEsWNrKgh1WIYGSj4tH158ZzDOKPCWjE6IKbqGfBZCUynqhwVEQJ6Aq4UcHYg2IYlQDO+BNldIbfMIh+94/naSlqqruru29NI/N+nVOnp+5z73PvU9VVn3nufao6VYUkSS08aLkHIEl64DBUJEnNGCqSpGYMFUlSM4aKJKmZ7Zd7AMvpUY96VK1cuXK5hyFJv1LWrFnzk6p69Kht23SorFy5kunp6eUehiT9Sknyndm2Of0lSWrGUJEkNWOoSJKaMVQkSc0YKpKkZgwVSVIzhookqRlDRZLUzDb94Udpa0myVfrx7yNpuRkq0law0Df7JAaEfiU5/SVJasZQkSQ1Y6hIkpoxVCRJzUw0VJLsmuS8JDcnuSXJx5LsNmbbmuV2wFC9ByU5Kcn6JHcluSrJ0ZM5IknSXCZ29VeShwGXAJuBFwIFvBW4NMkTqur2MXZzNnDmUNm/D90/FXgdcDKwBjgWODfJEVV14eKPQJK0UJO8pPglwJ7APlV1HUCSrwPXAi8FTh9jHz+oqitm25hkZ7pAOa2q3tUXX5pkL+A0wFCRpK1oktNfRwJXzAQKQFWtA74IHNWoj8OAHYDVQ+Wrgccn2aNRP5KkMUwyVPYHrh5RvhbYb8x9vCzJ5iR3JLkkyR+O6GMzcN1Q+dr+57j9SJIamGSo7ARsGlG+EVgxRvvVwMuBQ4ETgd8ELkly8FAfN9WWHz3eOLD9PpKcmGQ6yfSGDRvGGIYkaVyT/pqWUd8zMdaXIFXVCwbu/kuS8+nOfN4KHDSwrwX1UVVnAWcBTE1N+T0YktTQJM9UNjHiTIHuLGXUGcycqupW4NPAkweKNwIrsuW39a0Y2C5J2komGSpr6dY8hu0HfGOR+xw+M1kL7Ag8bkQfLKEfSdIiTDJUPgkcmGTPmYIkK4Gn99sWJMkjgWcBXx4ovgi4GzhuqPrxwNX91WaSpK1kkmsq7wdeAZyf5I10ZxinAt9j4AONSXYHrgdOqapT+rLXAfsAlwI/BHan+zzKLgwESFXdmOQM4KQktwJXAscAh9DusmVJ0pgmFipVdXuSQ4AzgHPopq4+B7ymqm4bqBpgO+571nQN8Jz+9uvALXSfb3lRVf3rUFcnA7cBr6YLnWuA51XVBc0PSpI0p2zLfwhoamqqpqenl3sY0hb8I126P0uypqqmRm3zW4olSc0YKpKkZgwVSVIzhookqRlDRZLUjKEiSWrGUJEkNWOoSJKaMVQkSc0YKpKkZgwVSVIzhookqRlDRZLUjKEiSWrGUJEkNWOoSJKaMVQkSc0YKpKkZgwVSVIzhookqRlDRZLUjKEiSWrGUJEkNWOoSJKaMVQkSc0YKpKkZgwVSVIzhookqRlDRZLUjKEiSWpmoqGSZNck5yW5OcktST6WZLcx2k0lOSvJt5LckeS7ST6UZI8RddcnqRG3Z0/mqCRJs9l+UjtO8jDgEmAz8EKggLcClyZ5QlXdPkfzY4H9gb8F1gKPBd4ETCc5oKq+N1T/YmDVUNk1Sz4ISdKCTCxUgJcAewL7VNV1AEm+DlwLvBQ4fY6276iqDYMFSb4IrOv3++ah+j+pqitaDVyStDiTnP46ErhiJlAAqmod8EXgqLkaDgdKX/YdYAPdWYsk6X5okqGyP3D1iPK1wH4L3VmSfYGdgW+O2Pxf+rWXzUmucD1FkpbHJENlJ2DTiPKNwIqF7CjJ9sD/pjtT+YehzRcArwQOA44D7gI+nuT4WfZ1YpLpJNMbNmxxQiRJWoJJrqlAtzg/LIvYz3uApwHPqqr7BFVVvfI+O08+DlwBvB1YvcWAqs4CzgKYmpoaNT5J0iJN8kxlE93ZyrAVjD6DGSnJ24ETgT+vqs/MV7+qfg6cC/xOkseM248kaekmeaaylm5dZdh+wDfG2UGSk4G/BF5VVecsoO+ZsyHPRCRpK5rkmcongQOT7DlTkGQl8PR+25ySvIrucy0nV9Xfjdtpv/7yXOC7VfWjBY5ZkrQEkwyV9wPrgfOTHJXkSOB84HvAmTOVkuye5J4kbx4oOxZ4N3ARcEmSAwdu+w3U+9MkH07yZ0me2be7FPh94PUTPDZJ0ggTm/6qqtuTHAKcAZxDNyX1OeA1VXXbQNUA23HfgDu8Lz+8vw36PHBw/+91dJcZv5Nu/eYO4CvA4VV1ccvjkSTNL1Xb7rLD1NRUTU9PL/cwpC0kYVt+ber+Lcmaqpoatc1vKZYkNWOoSJKaMVQkSc0YKpKkZgwVSVIzhookqRlDRZLUjKEiSWrGUJEkNWOoSJKaMVQkSc0YKpKkZib954SlB5yddtqJTZvG/uOli5Ys5i9vL8yKFSvYuHHjxPvRtsNQkRZo06ZND5hvEN4awaVti9NfkqRmDBVJUjOGiiSpGUNFktSMoSJJasZQkSQ1Y6hIkpoxVCRJzRgqkqRmDBVJUjOGiiSpGUNFktSMoSJJasZQkSQ1Y6hIkpqZaKgk2TXJeUluTnJLko8l2W3Mtg9J8s4kNyS5M8nlSZ4xot6DkpyUZH2Su5JcleTo9kcjSZrPxEIlycOAS4DfA14IvADYG7g0ycPH2MU/AC8B3gwcAdwAXJzkgKF6pwKrgPcA/xm4Ajg3yZ80OAxJ0gJM8i8/vgTYE9inqq4DSPJ14FrgpcDpszVM8kTg+cCfV9UH+rLPA2uBU4Aj+7KdgdcBp1XVu/rmlybZCzgNuHACxyVJmsUkp7+OBK6YCRSAqloHfBE4aoy2PwM+MtD2HuDDwGFJduyLDwN2AFYPtV8NPD7JHks6AknSgkwyVPYHrh5RvhbYb4y266rqjhFtdwD2Gqi3GbhuRD3G6EeS1NAkQ2UnYNOI8o3AiiW0ndk+8/Omqqp56v1SkhOTTCeZ3rBhwzzDkCQtxKQvKR5+swfIGO0yZttx6907oKqzqmqqqqYe/ehHjzEUSdK4JhkqmxhxpkB3ljLqLGTQxjnazmyf+bkiyXCIDNeTJG0FkwyVtXRrHsP2A74xRts9+suSh9vezb1rKGuBHYHHjajHGP1IkhqaZKh8EjgwyZ4zBUlWAk/vt83X9sHAcwfabg8cA3ymqjb3xRfRhcxxQ+2PB67urzaTJG0lk/ycyvuBVwDnJ3kj3drHqcD3gDNnKiXZHbgeOKWqTgGoqq8l+Qjw7iQPBtYBLwP2YCBAqurGJGcAJyW5FbiSLngOYf7LliVJjU0sVKrq9iSHAGcA59Atnn8OeE1V3TZQNcB2bHnW9N+AtwFvBX4DuAo4vKquHKp3MnAb8GpgF+Aa4HlVdUHbI5IkzSdbXo277Ziamqrp6enlHoZ+xSThgfK6eSAdi7aeJGuqamrUNr+lWJLUjKEiSWrGUJEkNWOoSJKameQlxdIDUr3lkbDq15d7GE3UWx653EPQA4yhIi1Q/uqWB8wVU0moVcs9Cj2QOP0lSWrGUJEkNWOoSJKaMVQkSc0YKpKkZgwVSVIzhookqRlDRZLUjKEiSWrGUJEkNWOoSJKaMVQkSc0YKpKkZgwVSVIzhookqRlDRZLUjKEiSWrGUJEkNWOoSJKaMVQkSc0YKpKkZgwVSVIzhookqZntl3sA0q+iJMs9hCZWrFix3EPQA8zEzlSSPCjJSUnWJ7kryVVJjh6j3SOTvDnJl5L8NMlN/b+fPaLuqiQ14vaJyRyVBFU18dvW6mfjxo3L/GjqgWaSZyqnAq8DTgbWAMcC5yY5oqounKPdbsDLgQ/0+/gF8KfAx5O8oqreO6LNQcDPB+77SpGkZTCRUEmyM12gnFZV7+qLL02yF3AaMFeorAP2rKo7BsouTrIr8HpgVKh8uaruaTB0SdISTGr66zBgB2D1UPlq4PFJ9pitYVXdPhQoM6aB3243RElSa5MKlf2BzcB1Q+Vr+5/7LWKfzwC+Ncu27yX5eZLvJHlHkocuYv+SpCWa1JrKTsBNNbPieK+NA9vHluRE4EDg+KFN1wF/CXwVKOA/Aa8FngT88Rz7OhFgt912W8gwJEnzGCtUkhwKfHaMqp+vqoOB0L3Jb7Gr8Yf2y74PBv4WOKeqPjS4raqGp9c+m+T7wLuTHFpV/3d4f1V1FnAWwNTU1KgxSpIWadwzlS8B+45Rb2YtZCOwIkmGzlZWDGyfV5InA58ELgFeNOZY/wl4N/BkYItQkSRNzlih0i+cz7aeMcpaYEfgcdx3XWVmLeUb8+0gyeOBi4GvAUdX1c8W0D+MPlOSJE3QpBbqLwLuBo4bKj8euLqq1s3VOMnedNNt3waOqKo7F9D3TJ9fXkAbSVIDE1mor6obk5wBnJTkVuBK4BjgEOCowbpJPgfsXlV79fd3pguUHYC3APsNfSXGV6tqc1/3q8AHgWvozkz+GHglcFFVXTqJY5MkzW6Sn6g/GbgNeDWwC90b//Oq6oKhetsNjWM/YPf+358asd89gPX9v68BXgE8pt/P9cApwN8sffiSpIXKllf9bjumpqZqenp6uYchbSEJ2/JrU/dvSdZU1dSobX71vSSpGUNFktSMoSJJasZQkSQ1Y6hIkpoxVCRJzRgqkqRmDBVJUjOGiiSpGUNFktSMoSJJasZQkSQ1Y6hIkpoxVCRJzRgqkqRmDBVJUjOGiiSpGUNFktSMoSJJasZQkSQ1Y6hIkpoxVCRJzRgqkqRmDBVJUjOGiiSpGUNFktSMoSJJasZQkSQ1Y6hIkpqZWKgkeVCSk5KsT3JXkquSHD1m27OT1Ijbu0fUPSjJl5LcmeRHSU5P8tD2RyRJms/2E9z3qcDrgJOBNcCxwLlJjqiqC8dovwE4cqjshsE7SZ4AfBa4GDgC2AN4J/BY4JgljV6StGATCZUkO9MFymlV9a6++NIkewGnAeOEyt1VdcU8df4K+D7w3Kr6Wd/33cA/JnlHVV25uCOQJC3GpKa/DgN2AFYPla8GHp9kj6V2kOTBwOHAR2cCpfdR4G7gqKX2IUlamEmFyv7AZuC6ofK1/c/9xtjHzkl+kuSeJP+e5PVJthvY/jjgIcDVg42q6i7g+jH7kCQ1NKk1lZ2Am6qqhso3Dmyfy9fo1mHW0gXHc4C3A3sDLx7ax6YR7TfO1keSE4ETAXbbbbd5hiFJWoixQiXJoXQL4vP5fFUdDAQYDhT68nlV1fBVXhcmuQ14Tb9Wcu3AvhbUT1WdBZwFMDU1NaqtJGmRxj1T+RKw7xj17uh/bgRWJMnQ2cqKge0L9U/Aa4Ap4FrmPutZwb1TbZKkrWSsUKmqO4BvLWC/a4Ed6dY9BtdVZtY5vrGAfc0YPjO5nm7dZv/7VEoeAuwJnLuIPiRJSzCphfqL6K7AOm6o/Hjg6qpat4h9Pp8uUL4CUFV39/08L8lgOP5XukD75CL6kCQtwUQW6qvqxiRnACcluRW4ku7DiIcwdKlvks8Bu1fVXv393YFzgA/TneXsSLdQfwJwZlVdP9B8FXA58NEk7wVW0n348byqWjOJY5MkzW6Sn6g/GbgNeDWwC3AN8LyqumCo3nZD47iVbr3k9cBv0Z2dfBN4FfC+wYZV9bUkhwHvAD4N3Ax8EHhD64ORJM0vW171u+2Ympqq6enp5R6GtIUkbMuvTd2/JVlTVVOjtvktxZKkZgwVSVIzhookqRlDRZLUjKEiSWrGUJEkNWOoSJKaMVQkSc0YKpKkZib5NS2SeslYf0poyW38FL6Wm6EibQW+2Wtb4fSXJKkZQ0WS1IyhIklqxlCRJDVjqEiSmjFUJEnNGCqSpGYMFUlSM9v036hPsgH4znKPQxrhUcBPlnsQ0ix2r6pHj9qwTYeKdH+VZLqqppZ7HNJCOf0lSWrGUJEkNWOoSPdPZy33AKTFcE1FktSMZyqSpGYMFUlSM4aKthlJLkty2cD9A5KsSrLTiLqVZNXWHN9skpydZP3A/ZX9+E5YwD4O7o917Nd8khP6flYOlK1PsnrcfSx2XIs5Rt0/GCralry8v804AHgLsEWoAE8F/n5rDGoRbqAb36cX0OZgumNdyGv+030/NyygzUIdzOhxLeYYdT/gnxPWNqOqvrGAuldMcixLUVWbgYmNL8mDgXuqagOwYVL9zGXSx6jJ8UxFS9JPX9Qst7PnaVtJ3pbk5CTfT3Jnkv+X5IChekny2iTXJLk7yQ1J3pPkkUP1Xp3km/1+NiWZTvKcge2/nP7qp1U+0G+6dmDMKwfGtmpo/4cnubzf/81JPpFkn6E6lyX5QpJDk1yZ5I4kVyd59piP53/s292V5PokLx1RZ4upoSRPTvLZJD/t+/x2kvf121bRnQ0A/GzmWIf29fIkf5Pkh8Bm4DdGTX8N9PeSJNf147wyyTNHPA6XjWi3fub3YsxxnTDU/vgkV/X9/iTJOUkeM6KP1UmO7X8fbu9/Fw4a/airJc9UtFRX0k1TDHom8NfAN8do/2fAd4FXADsCpwCfS7J3VW3s67wNOAl4L3ABsB9wKvDEJH9UVb9IchzwP/v2/wI8FHgCo6e2oJtWeSvwRuC5wPf78pFTPUkO79tcAhwDPKLv6wtJDqiqHwxUfxzwv4C3031/118A5yX5vaq6brYHIsm+wIXANHBs/3is6vv6+RztHgFcDPwrcAJwK7ASeFpf5e+B3wFeBBw0y75OBr4CnAhsB9w1W3/AHwG/37fZDLwe+OckT6yqa+ZoN2yccf1SkhOBM4GP0P0+/Dbd79lTkjypqm4bqP6HwD7Am/pjORX4VJKVVXXTAsaohaoqb96a3YC9gY3AefSfg5qjbtG96T58oGwl8DPg1P7+TnRvCmcPtT2+b39kf/89wJXz9HcZcNnA/RP6few1y9hWDdyfBq4Fth8o26Mf6+lDffwM2HugbGe6N8w3zDO+D414PHYF7gbWDz1GBZzQ35/q7z9hjn2v6utsP1Q+s68rh5+vgcdn5UDZ+n48uw2U/Vr/nJ8z22M91P7sBYxr5hi3A34MXDpU76C+3quG+tgErBgom3mMnr/cr5EH+s3pLzWTZAXwKeA64AXVv5rncWFV3T5zp6rW082lz5z9HEj3P/bhK44+DNxD979m6P6XfUCSv+unnh626AMZkuThwJOAj1TVPQNjXQd8cWAMM66tqmsH6t0I3AjsNk9XT2XLx+N7fR9zuRa4CTiznx7adZ76o3xizOcL4Iqq+u7AGG/l3kX9SdmHLpw/NFhYVV+g+6bx4efg8qraNHD/3/qf8z0HWiJDRU0k2Z7u7OQhdGcPd47Z9MezlD22//fM9NV9pqX6N/efDmz/IPAy4Cl0U0Ebk3xs1HrAIqwAMjyG3o/Ycopt44h6m+kem7k8htkfj1lV1c10U44/BN4HfLdfxzl6nv4GLeQKr/mes0kY+XvQm/c5qG7hH+Z/DrREhopaeS/wZOCIqvrRAtr91ixlM2sUM28OuwxW6EPsN+mCheqcWVV/QPe3SF4I/AHd/PtSbaKbOtllxLZdZsbQwA3M/njMqaq+VlVH0725PhW4Hvhokv8wZt8L+b6m+Z4z6KYsdxhRb7Y1rvmM/D0YKGv1HGiJDBUtWZLX0i22HltV/zZf/SF/0k8vzexrJd2U1+V90RV0/8s/dqjdMXQXmnx+eIdVtamqPgJ8FJjrTXXmf68PnWuA/XTUGuC5SbYbGOvudIvhW4xhkS5ny8djV+Dp4+6gqu6p7nLoN9G9vvftN411rGM6cHCKLcmvAc/i3ucMuimp302yw0C9Z9Ctvwwad1zX0J0N3ef3IMnTgN1p9xxoibz6S0vSv6jfRTf9tDHJgQObN1TV9fPs4k7gM0neSbd28lfALcAZAFW1McnpwElJbqe7Ompfuiu3vkD/4bgkZ9Fd9XQ53frF7wIvAD4zR98zn1v570n+kW6B/etVdfeIum/q+/pUf6nuI/qx3kx31VkLb6W7Em3m8dih72PO6a8kR9BdtfUJYB3wcOBV3Pt4wL3H+hdJ/hn4eVVNL3KcP+7HuIp7r/56ON0VVjM+3I/p//SXEO8B/A+6x2vQWOOqqp8neTPdutFqujW2x9JdGXgt914eruW23FcKePvVvnHvFUKjbmfP07bo3hTeQHdJ7110lwMfMFQvwGvp/rd6N9000XuBRw7UeSHdFUc30r3RraMLpsE6lzF0RRLd5yR+QHd11i+vdGLo6q++7HC6N+k76d4czwf2GapzGfCFEce6fr7Ho693KPDV/hi+DbwUOJu5r/7ah26ab13/GG6gC9+nDLTZrn/MbgR+0b3077OvF8/x3K4cOo7VwIvpptg29+M9ZET7l9K94d8JfInuMuT7PA5jjOuEoX0eD1zV9/tT4BzgMSMe69Wz/L6tGi731vbmV99r2fQfdHtbVb1xucciqQ3XVCRJzRgqkqRmnP6SJDXjmYokqRlDRZLUjKEiSWrGUJEkNWOoSJKa+f96rDBUwNAnSQAAAABJRU5ErkJggg==\n", | |
| "text/plain": [ | |
| "<Figure size 432x288 with 1 Axes>" | |
| ] | |
| }, | |
| "metadata": { | |
| "needs_background": "light" | |
| }, | |
| "output_type": "display_data" | |
| } | |
| ], | |
| "source": [ | |
| "y_position = data[' y position (m)']\n", | |
| "y_mean = np.mean(y_position)\n", | |
| "y_std = np.std(y_position)\n", | |
| "Q1_y = np.percentile(y_position, q=25)\n", | |
| "Q2_y = np.percentile(y_position, q=50)\n", | |
| "Q3_y = np.percentile(y_position, q=75)\n", | |
| "plt.boxplot(y_position, labels=['z position distribution']);\n", | |
| "print('The first quartile for z is {}'.format(Q1_y))\n", | |
| "print('The second quartile for z is {}'.format(Q2_y))\n", | |
| "print('The third quartile for z is {}'.format(Q3_y))" | |
| ] | |
| }, | |
| { | |
| "cell_type": "code", | |
| "execution_count": 285, | |
| "metadata": { | |
| "scrolled": true | |
| }, | |
| "outputs": [ | |
| { | |
| "name": "stdout", | |
| "output_type": "stream", | |
| "text": [ | |
| "The first quartile for tot is 0.19909935113674182\n", | |
| "The second quartile for tot is 0.2967735594999671\n", | |
| "The third quartile for tot is 0.40816954320848103\n" | |
| ] | |
| }, | |
| { | |
| "data": { | |
| "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX0AAAD9CAYAAABQvqc9AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAUyUlEQVR4nO3df7RlZX3f8ffHAWMkpZkp2KgVxqkUB4ixrsmSyCQOLF3SJAtitfFHNdA1BY31Z01a7SiiMFqjC2x+rMrgNFgwo41LM9QfLSqDdoikGYwahwkBy4Ck/hhhBEEZUb/9Y+9bD4dz79n33hnuMM/7tdZZZ85znuc837PPPZ+7z7P3PZOqQpLUhkcsdQGSpIeOoS9JDTH0Jakhhr4kNcTQl6SGHLbUBUxz1FFH1cqVK5e6DEl62Lj++uu/XVVHT7pvUOgneQJwMfBsIMCngddW1W0Dxh4DXACcChwF3A78N+AdVXXvtPErV65kx44dQ8qUJAFJbp3tvqmhn+TRwNXAPuAsoIALgW1JnjJXcCc5gu4XxOHAm4HbgF8E3gocB7xg+NOQJC3WkD39c4BVwPFVdTNAki8DNwEvAy6aY+wpdOH+nKq6qm/blmQF8DtJHl1V31tw9ZKkeRlyIPcM4LqZwAeoqluAa4Ezp4x9ZH9991j7d/q5M7BOSdJ+MCT0TwS+MqF9J3DClLGfpvtE8M4kJyT5mSSnAa8B3jtkTV+StP8MCf0VwN4J7XcCy+caWFX3AWv7eXYC3wU+A3wMeOW8KpUkLdrQ8/QnfSvb1KWZJI8CPgQ8Bngp8Ezgd+kO4P7RHOPOTbIjyY49e/YMLFF6aGzZsoWTTjqJZcuWcdJJJ7Fly5alLkkabMiB3L10e/vjljP5E8Co9cA64ElV9dW+7XNJ7gI2JXlvVX1pfFBVbQI2AaxZs8avAdVBY8uWLWzYsIHNmzezdu1atm/fzvr16wF40YtetMTVSdMN2dPfSbeuP+4E4IYpY38e2DsS+DP+d3+9esD80kFj48aNbN68mVNPPZXDDz+cU089lc2bN7Nx48alLk0aZEjoXwmcnGTVTEOSlXSnY145Zew3gOVJnjTW/vT++u+GlSkdHHbt2sXatWsf0LZ27Vp27dq1RBVJ8zMk9C8FdgNbk5yZ5AxgK/A14JKZTkmOTfLDJOeNjL2M7uDtJ5KcleTUJL8LvBu4nu60T+lhY/Xq1Wzfvv0Bbdu3b2f1aj+06uFhauj3p1WeBvwtcDnwAeAW4LSqumeka4Blo49ZVbuBk4Ev0v0V7yfo/thrE/DsqvrxfnkW0kNkw4YNrF+/nm3btnH//fezbds21q9fz4YNG5a6NGmQQd+903/HzvOm9NnNhDN6quoG4DcXUpx0sJk5WPuqV72KXbt2sXr1ajZu3OhBXD1s5GD/P3LXrFlTfuGaJA2X5PqqWjPpPr9PX5IaYuhLUkMMfUlqiKEvSQ0x9CWpIYa+JDXE0Jekhhj6ktQQQ1+SGmLoS1JDDH1JaoihL0kNMfQlqSGGviQ1xNCXpIYY+pLUEENfkhpi6EtSQwb9H7nSoS550H/vfMAc7P9FqQ5thr7EwoI4iQGuhx2XdySpIYa+JDXE0Jekhhj6ktQQQ1+SGmLoS1JDDH1JaoihL0kNMfQlqSGGviQ1xNCXpIYY+pLUEENfkhpi6EtSQwx9SWqIoS9JDTH0Jakhhr4kNcTQl6SGDAr9JE9I8uEkdyW5O8lHkhwzdJIkq5P8aZJvJ/l+khuTvGbhZUuSFmLqf4ye5NHA1cA+4CyggAuBbUmeUlX3Thm/ph9/DfCvgbuA44CfWVTlkqR5mxr6wDnAKuD4qroZIMmXgZuAlwEXzTYwySOA9wOfqarnjty1bcEVS5IWbMjyzhnAdTOBD1BVtwDXAmdOGbsOOIE5fjFIkh46Q0L/ROArE9p30gX6XNb2149Kcl2S+5N8K8nvJ/np+RQqSVq8IaG/Atg7of1OYPmUsY/rrz8EXAU8G/g9urX9P5ltUJJzk+xIsmPPnj0DSpQkDTFkTR+6g7fjMmDczC+VK6rqvP7f1yRZBvzHJCdU1Q0PmqxqE7AJYM2aNZPmliQtwJA9/b10e/vjljP5E8CoO/rrT421X9VfP3XA/JKk/WRI6O+kW9cfdwLwoL30CWPhwZ8UZj4l/HjA/JKk/WRI6F8JnJxk1UxDkpXAKf19c/kk3fn9p4+1P6e/3jGoSknSfjEk9C8FdgNbk5yZ5AxgK/A14JKZTkmOTfLDJDNr91TVHcA7gJcneXuSZyV5A3Ae8P7R00AlSQfe1AO5VXVvktOAi4HL6ZZmPgO8tqruGekaYBkP/kXyNuC7wCuA3wG+DrwLuGDR1UuS5mXQ2TtVdRvwvCl9djPhjJ6qKro/zvIPtCRpifktm5LUEENfkhpi6EtSQwx9SWqIoS9JDTH0Jakhhr4kNcTQl6SGGPqS1BBDX5IaYuhLUkMMfUlqiKEvSQ0x9CWpIYa+JDXE0Jekhhj6ktQQQ1+SGmLoS1JDDH1JaoihL0kNMfQlqSGGviQ1xNCXpIYY+pLUEENfkhpi6EtSQwx9SWqIoS9JDTH0Jakhhr4kNcTQl6SGGPqS1BBDX5IaYuhLUkMMfUlqiKEvSQ0x9CWpIYa+JDXE0JekhgwK/SRPSPLhJHcluTvJR5IcM9/JkrwxSSXZPv9SJUmLNTX0kzwauBp4MnAW8FLgOGBbkiOGTpRkFbAB+NbCSpUkLdZhA/qcA6wCjq+qmwGSfBm4CXgZcNHAuf4z8AHg+IHzSpL2syHLO2cA180EPkBV3QJcC5w5ZJIkLwaeBrxxIUVK87VixQqSHNALcMDnWLFixRJvSR1qhuxxnwhsndC+E/gX0wYnWQ5cDPy7qrpz5s0iHUh79+6lqpa6jEXz/aL9bcie/gpg74T2O4HlA8a/C/hb4LKhRSU5N8mOJDv27NkzdJgkaYqhp2xO2mWauguS5JeB3wJ+u+ax21VVm6pqTVWtOfroo4cOkyRNMWR5Zy/d3v645Uz+BDDqEmAzcHuSnx2Zc1l/+/tVtW9osZKkxRkS+jvp1vXHnQDcMGXs6v7y8gn37QVeB7xnQA2SpP1gSOhfCbw7yaqq+j8ASVYCpwBvmDL21Alt7wGWAa8Cbp5wvyTpABkS+pcCrwS2JnkT3fr+BcDX6JZvAEhyLPBV4G1V9TaAqrpm/MGSfAc4bNJ9kqQDa+qB3Kq6FziN7gycy+n+wOoW4LSqumeka+j24P0+H0k6SA36y9iqug143pQ+uxlwRk9VrRsypyRp/3OvXJIaYuhLUkMMfUlqiKEvSQ0x9CWpIYa+JDXE0Jekhhj6ktQQQ1+SGmLoS1JDDH1JaoihL0kNMfQlqSGGviQ1xNCXpIYY+pLUEENfkhpi6EtSQwx9SWqIoS9JDRn0H6NLDzf1liPh/L+/1GUsWr3lyKUuQYcYQ1+HpLz1bqpqqctYtCTU+UtdhQ4lLu9IUkMMfUlqiKEvSQ0x9CWpIYa+JDXE0Jekhhj6ktQQQ1+SGmLoS1JDDH1JaoihL0kNMfQlqSGGviQ1xNCXpIYY+pLUEENfkhpi6EtSQwx9SWrIoNBP8oQkH05yV5K7k3wkyTEDxq1JsinJ3yT5XpLbknwgyRMXX7okab6mhn6SRwNXA08GzgJeChwHbEtyxJThLwROBH4f+GfAG4CnATuSPGERdUuSFmDIf4x+DrAKOL6qbgZI8mXgJuBlwEVzjH1nVe0ZbUhyLXBL/7jnLaRoSdLCDFneOQO4bibwAarqFuBa4My5Bo4Hft92K7AHePz8SpUkLdaQ0D8R+MqE9p3ACfOdMMlq4DHArvmOlSQtzpDlnRXA3gntdwLL5zNZksOA99Lt6W+eo9+5wLkAxxwz9XixNFGSpS5h0ZYvn9dbTJpqSOgD1IS2hbyj/hB4BvBrVTXpF0k3WdUmYBPAmjVrJs0tzanqwP/YJHlI5pH2pyGhv5dub3/cciZ/ApgoyTvo9t7Pqqqrho6TJO0/Q0J/J926/rgTgBuGTJJkA93pmq+uqsuHlydJ2p+GHMi9Ejg5yaqZhiQrgVP6++aU5NXAhcCGqvqDhZUpSdofhoT+pcBuYGuSM5OcAWwFvgZcMtMpybFJfpjkvJG2FwLvAf4HcHWSk0cu8z7zR5K0OFOXd6rq3iSnARcDl9MdwP0M8Nqqumeka4BlPPAXyel9++n9ZdRngXULrlySNG+Dzt6pqtuA503ps5uxM3qq6mzg7IWVJkna3/yWTUlqiKEvSQ0x9CWpIYa+JDXE0Jekhhj6ktQQQ1+SGmLoS1JDDH1JaoihL0kNMfQlqSGGviQ1xNCXpIYY+pLUEENfkhpi6EtSQwx9SWqIoS9JDTH0Jakhhr4kNcTQl6SGGPqS1BBDX5IaYuhLUkMMfUlqiKEvSQ0x9CWpIYa+JDXE0Jekhhj6ktQQQ1+SGmLoS1JDDH1JaoihL0kNMfQlqSGHLXUB0sEgyUM2rqoWNJe0Pxj6Egax2uHyjiQ1xNCXpIYY+pLUEENfkhpi6EtSQwx9SWqIoS9JDTH0JakhOdj/KCXJHuDWpa5DmuAo4NtLXYQ0wbFVdfSkOw760JcOVkl2VNWapa5Dmg+XdySpIYa+JDXE0JcWbtNSFyDNl2v6ktQQ9/QlqSGGviQ1xNA/RCT5jST/dqnGH2hJdie5bOT2uiTnJ3nEWL+VSSrJ2Q91jZMkuSbJNSO31/X1rZvHY8z7tem3TY21VZIL5/M4C6lrIc9RDx1D/9DxG8BiQnux4w+05wIXjNxeB7yFB/8Mfx34JeDjD01Z8/YFuvq+MI8xC3lt3tfPcyDNVtdCnqMeIv53iXpYqKq/GthvH3DdAS5nwarqbg5gfUl+qqr2VdXtwO0Hap65HOjnqMVxT/8Q0C97nAU8vv9YXUl2j9x/fJKPJvlOku8nuS7J6UPHT5hvZgnlFUkuSvKtJN9L8rEkK8f6Hp7kwn555gf99YVJDh/pc1iSC5J8Ncl9Sb6dZHuStSN9/v/yTpLz6fbyAe6fqXmstrPH6nhJki+NPP7lSR471md3kiuSvDDJriT3JtkxWsdc+nF/k2Rfkp1Jnjuhz4OWPpI8J8mfJ7kryT1JbkxyXn/fZczy2ow81j9Pcmn/lSXfnNlG48s7P5kuG5Lc3v8sfC7JUydsh8smDKx+2w+ta/Q5Jsnr+uf2gyRfT/KHSY6cMMeFSV6d5JYk303y2SQnzr31NZR7+oeGC4CjgV8Ezujb9gEkeRywHfgu8ErgLuDfAB9P8utV9cm5xk/xRuCLwL8CHgO8HbgqyYlVdX/f5/3Ab/b3baf72P8mYBXw4r7PvwdeB2zoH+9IYA2wYpZ53wf8I2A9sBb40VxFJjkXuAT4UF/z4/p6np7kaVV1z0j3XwaOB94M3Ee3bT6WZGVVfWeOOZ4F/AndstLr6bbnfwIOB26cY9wq4Ergw8DbgB8Ax9FtHxj22vwB8EngpcCjZpur91vAbXQ/Cz/Vz/mZJMdV1Z1Txo6a78/MRrpt/0fAfwdO6B/jF5I8s6p+PNL3JXTb7DXAI4F3AVuTPLmqfjiPGjVJVXk5BC7AZcDtE9rfDfwQeNJI2zK6N9UXpo2fZa6VQAE3AI8YaT+lb1/f3z6pv33+2Pg39e1P6W9/DPjIlDl3A5eN3D6/f4zDZqnt7JHn+k1g21i/tX2/V4/NsRdYPtK2pu/34in1XTthezy9H3vNSNu6vm1df/v5/e0jF/DazjzWRyfcd3739n5AW9F9QdwRY9vrfuCC2bb12Pjz51HXzHNcQfcL9LKxfi/p+50xNsdNwOEjbTPb6BlL+R47VC4u7xz6fgW4rqpunmmoqh8BW4Cnjn+8nqcP18geWlVdS7eOPHMA8Vf66yvGxs3cfmZ//ZfArybZmGRtkkcuoqZxx9N9CvnAaGNVbaf79tZnjvX/fFXtHbn91/31MbNNkGQZ3R7v+Pb4C7oAncsX6UL3g0men+QxU/pP8tF59P1EVd07c6OqdtOtvx/Ig74n032qGP85+CDdDsn4a/Cp+sknRRjwGmg4Q//Qt4LujJZx3wACLF/EY39zlrbHj8zNhPm/MXb/2+nW6M8A/hdwR5I/TnLUImqbMVsNM3WMLyE9YImjugPDMPeyyVF0yzizbY9Z9b+Mn0P3Xrwc+EaSv0gyHoRzmfTcZjPtNTsQJr4G1S3V3MGU14CfLBtNW7rSAIb+oe9O4OcmtP8c3Ufm+azjjvuHs7T93cjcM3ONzw3dG56qur+q3llVPw88lm59/3l067+LNVsNM2137Ic5vk23tz7b9phTVW2rqtOBnwWe1T/Wx+fxS28+36Uy7TWDbinmAZ+2ksx2fGWIia9BksOAf8D+eQ00kKF/6NgH/PSE9s8CJ4+eVdMvR7wA+Kuq+u6U8XN5fkb+OCrJKXQHWD8/MjfAC8fG/cv++nPjD1hV36iq9wGfpjsmMJuZvb9pNd9Ityf7gBqSPAM4dqTGBeuXy/6SB2+Pp9OtmQ99nH1VdTXwe8ARwBP7uxby2szmV5McMVLjSrrll8+P9LmVB2/7X5/wWEPruq7vO/5z8AK6k0kW/RpoOM/eOXTcAKxI8tvADuC+qvpr4GLgbOBTSd4C3A28AvgnwK8NGD+Xvwf8WZJL6M7keAfdQbj/ClBVO5NsAc7v9+r+nG7t+M3Alqr6MkCSrcCX6P6YZy/wT4HT6c64mev5Arw+ySeBH1XVjvFOVfWj/vTHS5JcQbeu/Hi6s0luAv54ynMc6i3AVTxwe7yVnyxlTZTk5XTHPj4BfI1uqeiNwP8FvtJ3W8hrM5vv051h9S66dfa30v1MXDzS54PAf0lyMd1B9l+g+xkaN6iuqrozyUXAG5Pc2z/X1cCFdGd0Hax/SHdoWuojyV72z4Vuz3ALXWgWsHvkvuOBP6M7XfM+uj2v04eOnzDXyr7PK4CLgD3A9+jevE8c63s43Zv7Vrpli1v726NnZ7y+r+kOulC6ke7sk9E+u3ng2TvL6JZ/vgX8mP5MFcbO3hnp/xK6Xyz7+nkuBx471mc3cMWE5/ugM5Bm2S4v6mvfB+yk+yvia5j77J1fArbSBf4+unXvPwWOn/bajDzWsybUcv7MNhl7HhuB/0B3wP0+umMoTx3r9wjgvP61+h7wP4F/PL4dBtS1bqRv6JbtbqQ7LfXr/et35IQaL5zl5+3s8efpZf4Xv1pZ89YvCdwCnFPdUoykhwnX9CWpIYa+JDXE5R1Jaoh7+pLUEENfkhpi6EtSQwx9SWqIoS9JDfl/ytqKO5KjEn4AAAAASUVORK5CYII=\n", | |
| "text/plain": [ | |
| "<Figure size 432x288 with 1 Axes>" | |
| ] | |
| }, | |
| "metadata": { | |
| "needs_background": "light" | |
| }, | |
| "output_type": "display_data" | |
| } | |
| ], | |
| "source": [ | |
| "dist = np.sqrt(data[' x position (m)']**2 + data[' y position (m)']**2)\n", | |
| "dist_mean = np.mean(dist)\n", | |
| "dist_std = np.std(dist)\n", | |
| "Q1_tot = np.percentile(dist, q=25)\n", | |
| "Q2_tot = np.percentile(dist, q=50)\n", | |
| "Q3_tot = np.percentile(dist, q=75)\n", | |
| "plt.boxplot(dist, labels=['tot position distribution']);\n", | |
| "print('The first quartile for tot is {}'.format(Q1_tot))\n", | |
| "print('The second quartile for tot is {}'.format(Q2_tot))\n", | |
| "print('The third quartile for tot is {}'.format(Q3_tot))" | |
| ] | |
| }, | |
| { | |
| "cell_type": "markdown", | |
| "metadata": {}, | |
| "source": [ | |
| "# Section 2" | |
| ] | |
| }, | |
| { | |
| "cell_type": "code", | |
| "execution_count": 303, | |
| "metadata": {}, | |
| "outputs": [ | |
| { | |
| "name": "stdout", | |
| "output_type": "stream", | |
| "text": [ | |
| "The first quartile for z_height is -0.1789018487947081\n", | |
| "The second quartile for z_height is -0.005914843444322249\n", | |
| "The third quartile for z_height is 0.2602416017739394\n" | |
| ] | |
| }, | |
| { | |
| "data": { | |
| "text/plain": [ | |
| "array([-0.15202024, -0.26453083, -0.00591484, -0.15769379, -0.40009519,\n", | |
| " 0.16552245, 0.64461535, 0.18039022, 0.4106879 , 0.2602416 ,\n", | |
| " -0.10482995, -0.19368973, 0.16953191, 0.16857507, -0.25959889,\n", | |
| " -0.19530618, 0.11614742, -0.24414948, -0.03817304, 0.14973499,\n", | |
| " -0.0493952 , 0.4239159 , -0.05896327, -0.10622657, 0.33590505,\n", | |
| " 0.15311425, -0.17890185, -0.19615454, -0.11545261, -0.42585763,\n", | |
| " -0.30865695, 0.32180157, 0.27525837, -0.023196 , 0.26131531,\n", | |
| " -0.36320222, 0.47176472, -0.36319588, 0.01933821, 0.20214081,\n", | |
| " 0.27905946, -0.24434303, -0.16188914, 0.02955067, 0.54633433,\n", | |
| " 0.38801365, 0.44864213, 0.10254686, -0.00918668, -0.59019435,\n", | |
| " -0.03052037, 0.01370724, 0.33712062])" | |
| ] | |
| }, | |
| "execution_count": 303, | |
| "metadata": {}, | |
| "output_type": "execute_result" | |
| }, | |
| { | |
| "data": { | |
| "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZUAAAD9CAYAAAB0i+q4AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAYwElEQVR4nO3de7SkVX3m8e/D1csaku4g3mLTIA4RRmUyxyVRogyLDIwyoIMiKkYTRxwd42V0LUPw0gJGjC5wEp01oBmJtAOKS8QLikbQxEs7HlAMjRJgaNQEBO1GbtIg/uaP9z2hurrOOXXO2dUHmu9nrVrVtd+9373fqtP1nHfvt+qkqpAkqYUdlnsAkqTth6EiSWrGUJEkNWOoSJKaMVQkSc3stNwDWE677757rV69ermHIUkPKJdeeunPquoRo7Y9qENl9erVTE9PL/cwJOkBJcn1s21z+kuS1IyhIklqxlCRJDVjqEiSmjFUJEnNGCqSpGYMFUlSM4aKJKmZB/WHH6VtJck26ce/j6TlZqhI28BC3+yTGBB6QHL6S5LUjKEiSWrGUJEkNWOoSJKaMVQkSc0YKpKkZgwVSVIzhookqRlDRZLUjKEiSWrGUJEkNWOoSJKaMVQkSc0YKpKkZgwVSVIzhookqZmJhkqSxyX5ZJJfJLk1yaeSrBqzbc1yO2Co3g5JTkiyIcldSS5PcvRkjkiSNJeJ/eXHJA8DLgY2Ay8DCjgFuCTJk6vqjjF2cxZwxlDZPw49Phl4M3AicClwLHBekiOq6sLFH4EkaaEm+eeEXwnsDexbVdcAJPk+cDXwKuC0MfbxT1W1braNSfagC5RTq+p9ffElSfYBTgUMFUnahiY5/XUksG4mUACq6jrgG8BRjfo4DNgFWDtUvhZ4UpK9GvUjSRrDJENlf+CKEeXrgf3G3Merk2xOcmeSi5P8/og+NgPXDJWv7+/H7UeS1MAkQ2UlsGlE+UZgxRjt1wKvAQ4Fjgd+C7g4ycFDfdxSVTWij5ntW0hyfJLpJNM333zzGMOQJI1rkmsq0C3OD8tYDateOvDw75NcQHfmcwpw0MC+FtRHVZ0JnAkwNTU1qq0kaZEmeaayiRFnCnRnKaPOYOZUVbcBnweeOlC8EViRZDhEVgxslyRtI5MMlfV0ax7D9gOuXOQ+h89M1gO7Ao8f0QdL6EeStAiTDJXPAAcm2XumIMlq4Bn9tgVJshvwHODbA8VfBO4GXjJU/Tjgiv5qM0nSNjLJNZUPAa8FLkjyVrozjJOBHzPwgcYkewLXAidV1Ul92ZuBfYFLgH8G9qT7PMqjGAiQqropyenACUluAy4DXggcQrvLliVJY5pYqFTVHUkOAU4HzqabuvoK8Iaqun2gaoAd2fKs6Srgef3tN4Bb6T7f8oqq+r9DXZ0I3A68ni50rgKOqarPNj8oSdKcsvXVuA8eU1NTNT09vdzDkLaShAfz/03dvyW5tKqmRm3zW4olSc0YKpKkZgwVSVIzhookqRlDRZLUjKEiSWrGUJEkNWOoSJKaMVQkSc0YKpKkZgwVSVIzhookqRlDRZLUjKEiSWrGUJEkNWOoSJKaMVQkSc0YKpKkZgwVSVIzhookqRlDRZLUjKEiSWrGUJEkNWOoSJKaMVQkSc0YKpKkZgwVSVIzEw2VJI9L8skkv0hya5JPJVk1RrupJGcm+WGSO5P8KMnHkuw1ou6GJDXi9tzJHJUkaTY7TWrHSR4GXAxsBl4GFHAKcEmSJ1fVHXM0PxbYH/hLYD3wWOBtwHSSA6rqx0P1LwLWDJVdteSDkCQtyMRCBXglsDewb1VdA5Dk+8DVwKuA0+Zo+56qunmwIMk3gOv6/b59qP7Pqmpdq4FLkhZnktNfRwLrZgIFoKquA74BHDVXw+FA6cuuB26mO2uRJN0PTTJU9geuGFG+HthvoTtL8kRgD+AHIzb/p37tZXOSda6naJJWrlxJkonegIn3kYSVK1cu87Op7c0kp79WAptGlG8EVixkR0l2Av4X3ZnKXw9t/izwHbqpsUcCrwXOT/LSqlo7Yl/HA8cDrFo17zUD0lY2bdpEVS33MJqYCTCplUmGCnSL88MW81P8AeDpwHOqaougqqo/2WLnyfnAOuDdwFahUlVnAmcCTE1NbR/vDJJ0PzHJ6a9NdGcrw1Yw+gxmpCTvpjuz+OOq+tJ89avqXuA84LeTPHrcfiRJSzfJM5X1dOsqw/YDrhxnB0lOBP4UeF1Vnb2AvmfOhjwTkaRtaJJnKp8BDkyy90xBktXAM/ptc0ryOrrPtZxYVX81bqf9+ssLgB9V1Y0LHLMkaQkmGSofAjYAFyQ5KsmRwAXAj4EzZiol2TPJr5K8faDsWOD9wBeBi5McOHDbb6Dei5Kcm+QPk/z7vt0lwL8D3jLBY5MkjTCx6a+quiPJIcDpwNl0U1JfAd5QVbcPVA2wI1sG3OF9+eH9bdDXgIP7f19Hd5nxe+nWb+6kuxLs8Kq6qOXxSJLml+3l0sjFmJqaqunp6eUehh5gkmxXlxRvL8eibSfJpVU1NWqb31IsSWrGUJEkNWOoSJKaMVQkSc0YKpKkZgwVSVIzhookqRlDRZLUjKEiSWrGUJEkNWOoSJKaMVQkSc0YKpKkZgwVSVIzhookqRlDRZLUjKEiSWrGUJEkNWOoSJKaMVQkSc0YKpKkZgwVSVIzhookqRlDRZLUjKEiSWrGUJEkNWOoSJKamWioJHlckk8m+UWSW5N8KsmqMds+JMl7k9yQ5JdJvpXkmSPq7ZDkhCQbktyV5PIkR7c/GknSfCYWKkkeBlwM/A7wMuClwBOAS5I8fIxd/DXwSuDtwBHADcBFSQ4YqncysAb4APAfgXXAeUme3eAwJEkLsNME9/1KYG9g36q6BiDJ94GrgVcBp83WMMlTgBcDf1xVH+nLvgasB04CjuzL9gDeDJxaVe/rm1+SZB/gVODCCRyXJGkWk5z+OhJYNxMoAFV1HfAN4Kgx2t4DfHyg7a+Ac4HDkuzaFx8G7AKsHWq/FnhSkr2WdASSpAWZZKjsD1wxonw9sN8Yba+rqjtHtN0F2Geg3mbgmhH1GKMfSVJDkwyVlcCmEeUbgRVLaDuzfeb+lqqqeer9iyTHJ5lOMn3zzTfPMwxJ0kJMck0FYPjNHiBjtMuYbcetd9+Aqs4EzgSYmpoa1VaaU71jN1jzG8s9jCbqHbst9xC0nZlkqGxixJkC3VnKqLOQQRuBUZcerxjYPnO/IkmGzlaG60nN5J23svXJ8QNTEmrNco9C25NJTn+tp1vzGLYfcOUYbffqL0sebns3962hrAd2BR4/oh5j9CNJamiSofIZ4MAke88UJFkNPKPfNl/bnYEXDLTdCXgh8KWq2twXf5EuZF4y1P444Ir+ajNJ0jYyyemvDwGvBS5I8la6tY+TgR8DZ8xUSrIncC1wUlWdBFBV30vyceD9SXYGrgNeDezFQIBU1U1JTgdOSHIbcBld8BzC/JctS5Iam1ioVNUdSQ4BTgfOpls8/wrwhqq6faBqgB3Z+qzpj4B3AacAvwlcDhxeVZcN1TsRuB14PfAo4CrgmKr6bNsjkiTNJ9vLguNiTE1N1fT09HIPQw8wSbavhfrt5Fi07SS5tKqmRm3zW4olSc0YKpKkZgwVSVIzhookqRlDRZLUjKEiSWrGUJEkNWOoSJKaMVQkSc0YKpKkZgwVSVIzhookqRlDRZLUjKEiSWrGUJEkNWOoSJKaMVQkSc0YKpKkZgwVSVIzhookqRlDRZLUzE7LPQDpgSjJcg+hiRUrViz3ELSdMVSkBaqqifeRZJv0I7Xm9JckqRlDRZLUjKEiSWrGUJEkNTOxUEmyQ5ITkmxIcleSy5McPUa73ZK8Pck3k/w8yS39v587ou6aJDXi9unJHJUkaS6TvPrrZODNwInApcCxwHlJjqiqC+dotwp4DfCRfh+/Bl4EnJ/ktVX1wRFtDgLuHXi8scH4JUkLNJFQSbIHXaCcWlXv64svSbIPcCowV6hcB+xdVXcOlF2U5HHAW4BRofLtqvpVg6FLkpZgUtNfhwG7AGuHytcCT0qy12wNq+qOoUCZMQ08pt0QJUmtTSpU9gc2A9cMla/v7/dbxD6fCfxwlm0/TnJvkuuTvCfJQxexf0nSEk1qTWUlcEtt/ZHgjQPbx5bkeOBA4LihTdcAfwp8FyjgPwBvBH4X+IM59nU8wKpVqxYyDEnSPMYKlSSHAl8eo+rXqupgIHRv8lvtavyh/UvfBwN/CZxdVR8b3FZVw9NrX07yE+D9SQ6tqr8d3l9VnQmcCTA1NeX3YEhSQ+OeqXwTeOIY9WbWQjYCK5Jk6GxlxcD2eSV5KvAZ4GLgFWOO9Rzg/cBTga1CRZI0OWOFSr9wPtt6xijrgV2Bx7PlusrMWsqV8+0gyZOAi4DvAUdX1T0L6B9GnylJkiZoUgv1XwTuBl4yVH4ccEVVXTdX4yRPoJtu+3/AEVX1ywX0PdPntxfQRpLUwEQW6qvqpiSnAyckuQ24DHghcAhw1GDdJF8B9qyqffrHe9AFyi7AO4D9hv52xXeranNf97vAR4Gr6M5M/gD4E+CLVXXJJI5NkjS7SX6i/kTgduD1wKPo3viPqarPDtXbcWgc+wF79v/+3Ij97gVs6P99FfBa4NH9fq4FTgL+YunDlyQtVB7Mfwhoamqqpqenl3sY0lb8I126P0tyaVVNjdrmtxRLkpoxVCRJzRgqkqRmDBVJUjOGiiSpGUNFktSMoSJJasZQkSQ1Y6hIkpoxVCRJzRgqkqRmDBVJUjOGiiSpGUNFktSMoSJJasZQkSQ1Y6hIkpoxVCRJzRgqkqRmDBVJUjOGiiSpGUNFktSMoSJJasZQkSQ1Y6hIkpoxVCRJzRgqkqRmJhYqSXZIckKSDUnuSnJ5kqPHbHtWkhpxe/+Iugcl+WaSXya5MclpSR7a/ogkSfPZaYL7Phl4M3AicClwLHBekiOq6sIx2t8MHDlUdsPggyRPBr4MXAQcAewFvBd4LPDCJY1ekrRgEwmVJHvQBcqpVfW+vviSJPsApwLjhMrdVbVunjrvBH4CvKCq7un7vhv4myTvqarLFncEkqTFmNT012HALsDaofK1wJOS7LXUDpLsDBwOfGImUHqfAO4GjlpqH5KkhZlUqOwPbAauGSpf39/vN8Y+9kjysyS/SvKPSd6SZMeB7Y8HHgJcMdioqu4Crh2zD0lSQ5NaU1kJ3FJVNVS+cWD7XL5Htw6zni44nge8G3gC8F+G9rFpRPuNs/WR5HjgeIBVq1bNMwxJ0kKMFSpJDqVbEJ/P16rqYCDAcKDQl8+rqoav8rowye3AG/q1kqsH9rWgfqrqTOBMgKmpqVFtJUmLNO6ZyjeBJ45R787+fiOwIkmGzlZWDGxfqHOANwBTwNXMfdazgvum2iRJ28hYoVJVdwI/XMB+1wO70q17DK6rzKxzXLmAfc0YPjO5lm7dZv8tKiUPAfYGzltEH5KkJZjUQv0X6a7AeslQ+XHAFVV13SL2+WK6QPkOQFXd3fdzTJLBcHw+XaB9ZhF9SJKWYCIL9VV1U5LTgROS3AZcRvdhxEMYutQ3yVeAPatqn/7xnsDZwLl0Zzm70i3Uvxw4o6quHWi+BvgW8IkkHwRW03348ZNVdekkjk2SNLtJfqL+ROB24PXAo4CrgGOq6rND9XYcGsdtdOslbwEeSXd28gPgdcD/HGxYVd9LchjwHuDzwC+AjwJ/1vpgJEnzy9ZX/T54TE1N1fT09HIPQ9pKEh7M/zd1/5bk0qqaGrXNbymWJDVjqEiSmjFUJEnNGCqSpGYMFUlSM4aKJKkZQ0WS1MwkP/woqZeM9QXdS27jZ1u03AwVaRvwzV4PFk5/SZKaMVQkSc0YKpKkZgwVSVIzhookqRlDRZLUjKEiSWrGUJEkNfOg/suPSW4Grl/ucUgj7A78bLkHIc1iz6p6xKgND+pQke6vkkzP9udapfszp78kSc0YKpKkZgwV6f7pzOUegLQYrqlIkprxTEWS1IyhIklqxlDRWJI8N8l/X672k5ZkQ5KzBh4fnGRNkh2G6q1OUklevq3HOEqSryb56sDjg/vxHbyAfSz4temfmxoqqySnLGQ/ixnXYo5R246honE9F1hKKCy1/aQ9Dzh54PHBwDvY+v/IDcDvAZ/fNsNasMvoxnfZAtos5rX5cN/PJM02rsUco7YR/5ywBFTVd8estxlYN+HhLFpV3coEx5dk16raXFU/AX4yqX7mMulj1NJ4pqJ59dNCLwMe2087VJINA9v3TXJ+kluS/DLJuiSHj9t+RH8zU0yvSXJakpuS3Jnkc0lWD9XdOckp/fTV3f39KUl2HqizU5KTk1yb5K4kP0vy9SQHDdT5l+mvJGvozlIA7pkZ89DYXj40juOSXD6w/7OTPHqozoYka5Mcm+QHSe5IMj04jrn07X6YZHOS9UmeN6LOVlNDSQ5L8s0kv0hye5Krkry933YWs7w2A/v6z0k+1H+t0U9nnqPh6a/7usuJSX7S/yz8XZIDRjwPZ41oWP1zP+64Bo8xSd7YH9vdSW5I8oEku43o45Qkr0tyXZLbknwtyf5zP/sal2cqGsfJwCOApwJH9mWbAZI8Bvg6cBvwWuAXwH8DPp/kiKr6wlzt53EC8D3gj4A9gD8HvpRk/6q6p6/zN8Ax/bav002LvBXYG3hxX+ctwBuBE/v97QZMAStn6ffDwG8DrwAOAu6da5BJjgfOAD7ej/kx/XieluR3q+r2geq/D+wLvA24i+65+VyS1VV1yxx9HAr8H7pptzfRPZ//A9gZuGqOdnsDnwE+CZwE3A08ge75gfFem78CvgC8FHjIbH31/hD4Ed3Pwq59n19J8oSq2jhP20EL/Zl5F91z/0Hgs8B+/T6ekuRZVfXrgbrH0T1nrwd2Ad4LXJDkd6rqVwsYo0apKm/e5r0BZwE/GVH+PuBXwD4DZTvS/ae9bL72s/S1GijgSmCHgfJn9OWv6B//m/7xmqH2b+3Ln9w//hzwqXn63ACcNfB4Tb+PnWYZ28sHjvWnwCVD9Q7q671uqI9NwIqBsqm+3ovnGd83RjwfT+vbfnWg7OC+7OD+8fP7x7st4rWd2df5I7at6d4+tigrui/BfPjQ83UPcPJsz/VQ+zULGNfMMa6kC+izhuod19c7cqiPq4GdB8pmnqOnL+f/se3l5vSXluqZwLqqumamoKruBc4BDhiefligT9bAb5hV9Q26efyZBeJn9vdrh9rNPH5Wf/8d4NlJ3pXkoCS7LGFMw/alO4v62GBhVX2d7huwnzVU/1tVtWng8T/096tm6yDJjnS/sQ8/H9+me4Oey/fo3tTPTfL8JHvMU3+U8xdQ98KqumPmQVVtoFv/mOSi/oF0Z0XDPwfn0v3CM/wafLnuO9OFMV4Djc9Q0VKtpLsiatiNQIAVS9j3T2cpe+xA34zo/8ah7X9Ot0ZyJPD3wM+TfCTJ7ksY24zZxjAzjuEpti2mgKpb+Ie5p5V2p5vmmu35mFUf9ofR/V8/G7gxybeTDL/RzmXUsc1mvtdsEka+BtVNZf2ceV4D7ptWm29qT2MwVLRUG4FHjSh/FN2UwkLm0Yc9cpayfxroe6av4b6he0Ohqu6pqvdU1ZOAR9OtrxxNN/++VLONYabs5w36+Bnd2cZsz8ecquqSqjoc+E3g0H5fn19AqC7ku5zme82gm6ra4mwxyWzrW+MY+Rok2Qn4Ldq8BhqToaJxbQYeOqL8a8CBg1dl9dM1LwS+W1W3zdN+Ls/PwIcPkzyDbgH9WwN9Axw71O4l/f3fDe+wqm6sqg8Df0u3JjObmd9e5xvzVXS/iW8xhiRPB/YcGOOi9dOJ32Hr5+NpdGsW4+5nc1VdDPwF8HBgr37TYl6b2Tw7ycMHxriabnrqWwN1rmfr5/6IEfsad1zr+rrDPwcvpLsYacmvgcbn1V8a15XAyiSvBqaBu6rqH4DTgZcDX07yDuBW4DXAvwaeM0b7ufwr4NNJzqC7EujddIusHwWoqvVJzgHW9L+VfpNu7v5twDlV9X2AJBcAl9N9WG4T8G+Bw+mu2JrreAHelOQLwL1VNT1cqaru7S/PPSPJWrp5/cfSXY10NfCReY5xXO8AvsSWz8c7uW+qb6Qk/5Vu7elC4Md0U2knAP8MXNFXW8xrM5tf0l2h9166dY530v1MnD5Q51zgfyc5ne4iiqfQ/QwNG2tcVbUxyWnACUnu6I/1icApdFcE3l8/qLp9Wu4rBbw9MG50v9meQ/emXMCGgW37Ap+mu5z4LrrfHA8ft/2Ivlb3dV4DnAbcDNxJ9+aw11DdnenePK6nm9a5vn88eHXPm/ox/ZzuTe8ququXButsYMurv3akmx67Cfg1/ZVODF39NVD/OLrg2tz3czbw6KE6G4C1I453qyvYZnleXtSPfTOwnu5bAL7K3Fd//R5wAV2gbKZbdzgP2He+12ZgX4eOGMuamedk6DjeBfwZ3QUVd9GtYR0wVG8H4O39a3UncBHw+OHnYYxxHTxQN3TTmlfRXTZ9Q//67TZijKfM8vP28uHj9Lbwm199r/udfsrkOuCV1U1VSXqAcE1FktSMoSJJasbpL0lSM56pSJKaMVQkSc0YKpKkZgwVSVIzhookqZn/D+6Xgzgo2ZtDAAAAAElFTkSuQmCC\n", | |
| "text/plain": [ | |
| "<Figure size 432x288 with 1 Axes>" | |
| ] | |
| }, | |
| "metadata": { | |
| "needs_background": "light" | |
| }, | |
| "output_type": "display_data" | |
| } | |
| ], | |
| "source": [ | |
| "def montecarlo_heights(D,V_low,V_high,N):\n", | |
| " '''Creates random points on the target based off of variable inputs.\n", | |
| " ---------\n", | |
| " D: Distance away from the target\n", | |
| " V_low: lower bound of the initial velocity\n", | |
| " V_high: upper bound of the initial velocity\n", | |
| " N: number of projectiles thrown\n", | |
| " -------\n", | |
| " '''\n", | |
| " import random\n", | |
| " pi = 3.14159\n", | |
| " g = 9.81\n", | |
| " z0 = .3\n", | |
| " V = np.zeros(N)\n", | |
| " for i in range (0,N):\n", | |
| " V[i] = V_low + np.clip(np.random.normal(.5,.1),0,1)*(V_high-V_low)\n", | |
| " theta = np.zeros(N)\n", | |
| " for i in range (0,N):\n", | |
| " theta[i] = np.random.rand(1)*pi/11\n", | |
| " cos = np.cos(theta)\n", | |
| " sin = np.sin(theta)\n", | |
| " z = D/cos*(sin-g/(2*V**2*cos)*D)+z0\n", | |
| " return z\n", | |
| "z_height = montecarlo_heights(3,4,12,53)\n", | |
| "Q1_z_height = np.percentile(z_height, q=25)\n", | |
| "Q2_z_height = np.percentile(z_height, q=50)\n", | |
| "Q3_z_height = np.percentile(z_height, q=75)\n", | |
| "plt.boxplot(z_height, labels=['tot position distribution']);\n", | |
| "print('The first quartile for z_height is {}'.format(Q1_z_height))\n", | |
| "print('The second quartile for z_height is {}'.format(Q2_z_height))\n", | |
| "print('The third quartile for z_height is {}'.format(Q3_z_height))\n", | |
| "z_height" | |
| ] | |
| }, | |
| { | |
| "cell_type": "markdown", | |
| "metadata": {}, | |
| "source": [ | |
| "# Section 3" | |
| ] | |
| }, | |
| { | |
| "cell_type": "markdown", | |
| "metadata": {}, | |
| "source": [ | |
| "I ajusted the angle variation until the positive regon of the box plot whisker reached around .6 meters. This occured when the angle was between 0 and pi/11. I then set the mean of the velocity to be 8 m/s and adjusted the standard deviation until the whisker plot was around -.6. This occured when the standard deviation was set to .1." | |
| ] | |
| }, | |
| { | |
| "cell_type": "code", | |
| "execution_count": null, | |
| "metadata": {}, | |
| "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.7.3" | |
| } | |
| }, | |
| "nbformat": 4, | |
| "nbformat_minor": 4 | |
| } |