Lugovtsov/4.1.1-centring-optical-system/processing.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 14,
   "id": "a21511dc-d48d-4171-bf38-7ed030488d14",
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "# %load ../processing_tools.py\n",
    "import numpy as np\n",
    "from scipy.optimize import curve_fit\n",
    "import pandas as pd\n",
    "\n",
    "def line(x, k, b):\n",
    "    return k * x + b\n",
    "\n",
    "import matplotlib.pyplot as plt\n",
    "import matplotlib\n",
    "import scienceplots\n",
    "\n",
    "plt.style.use(['science', 'russian-font'])\n",
    "\n",
    "matplotlib.rcParams.update({\n",
    "    'figure.figsize': [6, 4],\n",
    "    'savefig.facecolor': 'white',\n",
    "    'figure.dpi': 150.0,\n",
    "    'font.size': 12.0,\n",
    "})\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "id": "ebe7c967-eed8-4701-80a5-3f276e07347b",
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "h = np.array([10, 20, 40]) # mm\n",
    "ds = np.array([0, 1, 3]) # mm"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "id": "549aff15-dfb2-48b0-9802-f33cbfd7c3f5",
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "(k, b), _ = curve_fit(line, h**2, ds)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "id": "c6234bfd-5756-4422-a7c9-557be49f0e67",
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "x = np.linspace(0, 1700, 100)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "id": "2160928d-976d-4dd2-9074-ce026fb1f801",
   "metadata": {
    "tags": []
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAxgAAAIbCAYAAABonjh+AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy88F64QAAAACXBIWXMAABcSAAAXEgFnn9JSAACA/0lEQVR4nO3deVwTd/4/8BeJiceoBPDoaWXoYW8NePRuBXru1TbRdtvdXgLaY69WkLq2ttkWw3a32+1aBW277fZS0G53e5pgz21RIGqrrbYyau3hBQxoRAnJ/P7gl3xBQEMyYSbwej4ePB6STD7zzmfGYd7zuRIURVFARERERESkAoPWARARERERUd/BBIOIiIiIiFTDBIOIiIiIiFTDBIOIiIiIiFTDBIOIiIiIiFQzQOsA4kkgEIDX64XBYEBCQgIMBgOMRqPWYRERERERRcTv9yMQCEBRFAQCAQiCAIMhujYIJhg9cOjQIQwfPlzrMIiIiIiIYsLr9WLIkCFRlcEEowcGDGirLlmWYTKZwmrB8Pl8EAQBXq8XJpMpqv3rsSw9xqRmWXqMSc2y9BiTmmXpMSY1y9JjTGqWxZh6vyw9xqRmWXqMSc2y9BiTmmXpMSY1y9IqpmALhs/ng8ViCd3vRoMJRg8kJCQAAIYMGdLjA28ymaI+WfRclh5jUrMsPcakZll6jEnNsvQYk5pl6TEmNctiTL1flh5jUrMsPcakZll6jEnNsvQYk5pl9XZMwfd9Ph+A/7vfjQYHeceYwWDAww8/HHVfNrWpFZea30+vZalFj99PjzGpSa/fT4/HT016/H56jElNevx+ei1LLXr8fnqMSU16/X56PH5aS1AURdE6iHjh8/lgNpvR0tKiWmYZT/vvj1jn2mC99z7Wee9jnfc+1rk2WO+9L5I6V/M4xX+KREREREREusEEI470paazeME61wbrvfexznsf67z3sc61wXrvfVrXObtI9QCb+IiIiIioL2IXKSIiIiIi0iUmGEREREREpBomGEREREREpBomGEREREREpBqu5B2B4EqHBoMBRqNR42iIiIiIiCLj9/sRCARC97dqYAtGBARBgNlshsPh0DoUIiIiIqKIORwOmM1mCIKgWpmcprYHgtN3eb1emEwmtmAQERERUVxr34IhCIIq09Syi1QETCYT18EgIiIiorhnNBpVf2DOLlJERERERKQaJhhERERERKQaJhhERERERKQaJhhERERERDrQ0urXOgRVMMEgIiIiItKQrzWAP/9nIyYVvoXGgy1ahxM1ziJFRERERKSRjd82YNbSSqzf3gAAeO1/25CXfYbGUUWHCQYRERERUS/ztQbwlzc3ofiNTfD5A0gSzCi+NR3TLxyrdWhRY4JBRERERNSLPt/RgFlLKvH5t22tFtdZT8Lfbp+I4yyDNY5MHUwwiIiIiIh6QUurH3/+zyY88d9NaPUrSB46EE/8Kh22KacgISFB6/BUwwQjAj6fDwBgMBhUX/mQiIiIiPqe9dvrMWtJJTbulAEAP8s4GU/eloFRidq2Wvj9fgQCgdD9rRo4i1QEBEGA2WyGw+HQOhQiIiIi0rHDPj8eLd+Ay+e/h407ZaQMG4h/3n0RXrrvYs2TCwBwOBwwm80QBEG1MhMURVFUK62P8/l8MJvN8Hq9MJlMbMEgIiIiom7VSHWYtaQSX33fCAC4YdIYPPHrDIwcPkjjyP5P+xYMQRDQ0tICk8kUVZnsIhUBk8kUdcUTERERUd90qMWPon9/gb+99RUCioKRwwfhydsy8POJY7QOrROj0aj6A3MmGEREREREKlm7dR9mLanE1z82AQCmXXAKnLemY8Qw/bRaxBoTDCIiIiKiKDW3tOJPK77AP97djICiYFTiIPzt9on4afrJWofW65hgEBERERFF4bOv9+LupZXYums/AODmi8ZiwS3pSB46UOPItMEEg4iIiIgoAgcPt+LR8g14ZtUWKApwfNJgPHX7JFwz4UStQ9MUEwwiIiIioh76ZPMe3PtsJWp3HwAA3HKJiKJfWpEkmDWOTHtMMIiIiIiIwuQ93Ir5y9djsetrAMAJSYPx9J2TceX5J2gcmX4wwSAiIiIiCsPHX+3GPc+uwbY9ba0Wv74sDY/fPAGJQ9hq0R4TDCIiIiKiozhwyIeHlq3HkopvAAAnJQ/B03dOQtZ5bLXoChMMIiIiIqJufLBpF+59dg127PMCAO644lT86aYJGD6Yiy53hwlGBHw+HwDAYDCovvIhEREREWmvqdmHea+tw3PvbwUAjBkhYOFdk3H52cdpHJm6/H4/AoFA6P5WDQbVSupHBEGA2WyGw+HQOhQiIiIiUlnFFz9icuFboeQiJ/M0VD52bZ9LLgDA4XDAbDZDEATVykxQFEVRrbQ+zufzwWw2w+v1wmQysQWDiIiIqA9pPNiCB19dhxc/rAUAjB0pYOFdU3DpWaM1jix22rdgCIKAlpYWmEzRdf9iF6kImEymqCueiIiIiPRj1YYf8Jvn1+L7+oMAgJnZp2P+tPEQBvbt22Wj0aj6A/O+XWNEREREREfR4G1B4SsevPyxBAAQRw3FwhlTcPG4URpHFr+YYBARERFRv/T2uu/w2+ersEtuRkICcPeVZ+Ah2/kY0sdbLWJN97VXWlqKmpoaAEB9fT0kSUJeXh5yc3N7XJYkSXA6nbBYLKHfCwsLYbVa1QyZiIiIiHSs/sBhzHm5Bq/+bzsA4NTjhuGZGVNwwekjtQ2sj9D1IO+CggJMnz69QwLg8XiQnp6OrKwsuFyusMvyeDyw2+1wuVwQRREAIMsy0tPT4XQ6YbPZjllGcJC3GoNfiIiIiKh37WlsxoOvevD62p1oaQ0AAC44fSSW5l2AMSOHahydttS8z9VtgiFJEtLS0pCbm4uSkpIO72VnZ8PtdsPlciErKyus8tLS0pCXl4f8/PwOr5eXlyMnJwfbtm0LtWx0hwkGERERUfxpbmnFb55bi2WfbUdXd74mowG3Xiqi+JZ0DDL3zxlC1bzP1e06GLIsAwCqq6s7vde+i1M43G43JEnqMhmx2WyQZRmlpaURx0pERERE+tTc0opLH3oXr33adXIBAD5/AM+/vxXXP/E+mltaezfAPki3CYbVakVtbS0qKio6vRdMLDIyMsIqK9iVKtg16kgWiwXLli2LMFIiIiIi0qO9TYcwqfAtbP6hKaztP9m8BwUve2IcVd+n2wQDaEsIjuy2JEkSPB4PsrKywh6c7Xa7AaDbLlCiKMLj4clERERE1BcoioKVa3YgveBNbN/r7dFnX/pIwp7G5hhF1j/oOsHoSl5eHqxWK8rKysL+jCRJxxxfQURERETxb09jM371j09w28L/ocHb0uPP+/wBvPBheN3wqWtxkWAUFxfDbrcjLS0NVqsVNTU1PUoYguM51NrO5/N1+vH7/WHHQ0RERETqUhQFZZ9tx8TCt/FG1U4MMCZg7EghorI+2bxb5ej0xe/3d3k/q5a4SDDy8/NRVlYGl8sFSZKQnp4e9gDvnqivrw9rO0EQYDabO/w4HA7V4yEiIiKiY9slN+Pmpz7GnYs+Rf2Bwzh3jAUfzr8aI4cPiqi8/c3q3WzrkcPh6HQvKwiRJWNd0f1Ce+2JooiysjKkp6cjLS0NtbW13Q7cbk/t7lFer7fT9F0GQ1zkakRERER9hqIoWPbpduS/VIMGbwtMRgPyf3427v/J2TANMGDY4MimW430c/Fi3rx5mDt3bofXfD6faklGXN4V5+XlAQDsdntY2ycnJ6u6nclk6vRjNPbPOZOJiIiItPBjw0FM/9tHyCn5DA3eFowfm4SPHrkKc35xLkwD2m5xLx43KqKyLx43Ws1QdcdoNHZ5P6sW3SYYkiShvLy8y/eC09OGO/OTxWI56viK4HscCE5ERESkb4qi4OWPJUwqfAvvrPse5gEGPGQ7H6sfugrnjEnqsO2vL02Dydiz212T0YDbLjt2Dxnqnm4TjOzsbNjt9m6TjKBwBmYHE5LutpUkKewpb4mIiIhIG9/XH4TtLx9g5pJKyAd9sKYm4+NHr8bsn50darVob7RlMG65JLVH+7j1UhGjEge
      "text/plain": [
       "<Figure size 900x600 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.scatter(h**2, ds)\n",
    "plt.plot(x, line(x, k, b))\n",
    "\n",
    "plt.xlabel(r'Квадрат диаметра диафрагмы $h^2, \\text{мм}^2$')\n",
    "plt.ylabel(r'Отклонение $\\delta s$, мм')\n",
    "\n",
    "plt.savefig('images/dSonH2.svg')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 32,
   "id": "2f749ff7-6a94-4a4a-8339-ea843247e055",
   "metadata": {
    "tags": []
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "1.1904761904772945"
      ]
     },
     "execution_count": 32,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "line(625, k, b) - line(0, k, b)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "5889dd05-8f4c-42a1-8c27-fe2b8fbcbea7",
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "f_chr = 0.5e-3\n",
    "f_D = 5e-2"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "e350de8c-412c-4732-a65c-33220c837c19",
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "nu = f_D / f_chr"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "4e359270-52c7-47e4-ac2d-8685fb69d1d9",
   "metadata": {
    "tags": []
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "100.0"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "nu"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "4cf6e3cb-a91f-4dd7-b1ac-efd0b15b32ab",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "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.10.6"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}