381 lines
252 KiB
Plaintext
381 lines
252 KiB
Plaintext
|
{
|
||
|
|
"cells": [
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": 109,
|
||
|
|
"id": "f8d06107",
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [],
|
||
|
|
"source": [
|
||
|
|
"# -*- coding: utf-8 -*-\n",
|
||
|
|
"\"\"\"\n",
|
||
|
|
"Created on Tue Aug 24 16:24:52 2021\n",
|
||
|
|
"\n",
|
||
|
|
"@author: Joschka\n",
|
||
|
|
"\"\"\"\n",
|
||
|
|
"\n",
|
||
|
|
"import matplotlib.pyplot as plt\n",
|
||
|
|
"import numpy as np\n",
|
||
|
|
"import sys\n",
|
||
|
|
"sys.path.insert(0,'..')\n",
|
||
|
|
"\n",
|
||
|
|
"from src import coil_class as BC\n",
|
||
|
|
"from src import physical_constants as cs\n",
|
||
|
|
"\n",
|
||
|
|
"#from IPython import get_ipython\n",
|
||
|
|
"#get_ipython().run_line_magic('matplotlib', 'qt')\n",
|
||
|
|
"#get_ipython().run_line_magic('matplotlib', 'inline')\n",
|
||
|
|
"\n",
|
||
|
|
"\n",
|
||
|
|
"from scipy.optimize import curve_fit\n",
|
||
|
|
"plt.rcParams['axes.grid'] = True\n",
|
||
|
|
"plt.rcParams['grid.alpha'] = 0.5"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "markdown",
|
||
|
|
"id": "99eb6ef1-8920-4f7d-9b58-fd00c731d2bc",
|
||
|
|
"metadata": {},
|
||
|
|
"source": [
|
||
|
|
"## Set up coils"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": 274,
|
||
|
|
"id": "199e2602-bec6-429a-a639-6ea35ababb86",
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [],
|
||
|
|
"source": [
|
||
|
|
"HH = 1\n",
|
||
|
|
"I = 5\n",
|
||
|
|
"\n",
|
||
|
|
"d_coils = 50\n",
|
||
|
|
"Radius = 30\n",
|
||
|
|
"\n",
|
||
|
|
"layers = 4\n",
|
||
|
|
"windings = 4\n",
|
||
|
|
"wire_width = 1\n",
|
||
|
|
"wire_height = 1\n",
|
||
|
|
"\n",
|
||
|
|
"Coil = BC.BCoil(HH,d_coils,Radius, layers, windings, wire_width, wire_height)"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": 287,
|
||
|
|
"id": "40cee79d-d6ef-4df4-ae03-5e82785433e8",
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [
|
||
|
|
{
|
||
|
|
"name": "stdout",
|
||
|
|
"output_type": "stream",
|
||
|
|
"text": [
|
||
|
|
"20\n"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"data": {
|
||
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAkcAAAGwCAYAAACjPMHLAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAA9hAAAPYQGoP6dpAAEAAElEQVR4nOyddVxU2fvHPzMMjYB0qIgidq/djY3da7drrx3Y2Lq2ssbatYrdIipgYICCiqAgJSkdE8/vD35zvyAhA/eKuuf9evkS7tz5nDNnhrnPfc4TIiIiMBgMBoPBYDAAAOKSngCDwWAwGAzGjwQzjhgMBoPBYDCywYwjBoPBYDAYjGww44jBYDAYDAYjG8w4YjAYDAaDwcgGM44YDAaDwWAwssGMIwaDwWAwGIxsSEp6Aj8jCoUC4eHhKFWqFEQiUUlPh8FgMBgMRiEgIiQlJcHKygpicf7+IWYcFYHw8HCULVu2pKfBYDAYDAajCHz69AllypTJ93FmHBWBUqVKAchaXH19fd505XI5AgMDUbFiRaipqfGmy/g2bO1LDrb2JQdb+5KDrX3JkJiYiLJly3LX8fxgxlERUG6l6evr824c6enpQV9fn/2xfGfY2pccbO1LDrb2JQdb+5LlWyExLCCbwWAwGAwGIxvMOGIwGAwGg8HIBjOOGAwGg8FgMLLBjCMGg8FgMBiMbDDjiMFgMBgMBiMbzDhiMBgMBoPByAYzjhgMBoPBYDCywYwjBoPBYDAYjGww44jBYDAYDAYjG8w4+sGIj48XTJuI8OnTJ8H0ASAuLk5Q/czMTEH1GQwGg8FgxtEPxtWrVyGXywXRfvjwIc6ePSuINgD4+vriyJEjgul/+vQJx44dE0w/PT0dd+7cEUxfJpMhKChIMH0AyMjIEFSfwWAw/gsw4+gHwsvLCwEBAYiMjBREf8+ePXj69Kkg2gDg5OSE4OBgwfT/+OMPxMbGCqJNRLh48aKg+hMnTkR4eLhg+k5OTggLCxNEHwBOnz4tqPH18eNHwbQBCHbTwWAwfj2YcfQDsXr1agAQxMCIi4vD6dOn8eTJE961AeD58+f4999/ERISIoj++fPn4erqKti22u7du+Hn5wcDAwNB9J2cnODi4gJ7e3vetRUKBWbMmAEXFxfY2tryri+VSjF16lQcPnwYmpqavOunp6fj5s2bOHr0KO/aAJCSkoIVK1bg9u3bguinpaXBxcUFX758EUQ/MzMTjx49EkQbyDKshZp79jEYjJ8JZhz9IERHR8PMzAwmJiaC3P2/ePECDRo0gJGRERITE3nXDwwMRMWKFaGlpcW7NgCEh4ejXLlygnzJSqVSBAcHQ1tbG4aGhrzrf/nyBT4+PrC2toapqSnv+qGhoXB3d0erVq2+2Wm6KPj5+eHQoUPo27cv79oAcPHiRTx69AijRo0SRN/JyQmurq7o0KED79oKhQI9evRAYGCgIJ+d1NRUVKlSBUlJSbxrA4C/vz9cXFyQlpYmiP6///6LmTNnQiaTCaK/ceNGwYxehUKBFStWCBZHOWvWLMG8pYmJiejVqxeio6MF0ffw8MCECRME0QaAHTt2YM+ePYLpFwYRMZNeZRITE2FgYICEhATo6+vzpiuXyxEQEABbW1tB7tCVEJEgF1Egy9BQV1cXRJuIkJaWBh0dHd615XI53rx5gzJlygjiPSIifP78GRYWFrxrA1nxTFFRUbCyshJEPy4uDmpqaoKsjVwux+PHj9GwYUOoqanxri+TyTjjWgg+f/6MUqVKCfK5BICQkBDB5i6VSvHixQvUq1dPkLVPSEiAhoYGtLW1edcGsm4qhbjhUBIXFwcjIyNBtLdu3Yo6deqgefPmvK99amoqVq5ciQULFkBPT49XbSDLqL5y5QpmzZrFuzaQZVRra2ujc+fOvGsX9vrNjKMiILRxVKlSJUG+qBj5w9a+5GBrX3KwtS852NqXDIW9frNtNQaDwWAwGIxsMOOIwWAwGAwGIxvMOGIwGAwGg8HIBjOOGAwGg8FgMLLBjCMGg8FgMBiMbDDjiMFgMBgMBiMbzDhiMBgMBoPByAYzjhgMBoPBYDCywYyjHwRWi5PBYDAYjB8DZhz9IAjVdPNrrly5Iqj+69evBdWXyWSCdoZnMBgMBoMZRz8Ix44dw9WrVwUd482bN9iwYYNg+jKZDBMnThRMHwDWrl2LhIQEwfTDwsLw9u1bwfSBrC70DAaDwfhxYcbRD0JqaipmzJgBqVQq2Bi7d+/Gs2fPBNvC++eff3D//n1kZmYKov/q1SssW7YMYrEwH1upVIp+/fpBIpEIog9kdZv29/cXTD8yMhKRkZGC6TMYDMZ/AWYc/SCkpqbi48ePePjwoSD6KSkpOHjwIBISEvDhwwfe9dPT07F06VIAQFRUFO/6UqkUI0aMgFQqhUgk4l2fiHDhwgU8e/ZMsM72O3bswOzZs1G9enVB9J89e4aOHTvC2NhYEP2UlBSsXr1aUAM+PDxcMG0Gg8EoLMw4+kFITU0FADx48ACBgYG86584cYLbjnr27Bnv+jt37kRoaCgACOK5WLduHby9vQFAEM/Rhg0b8OrVKxgZGUFbW5t3/R07dmDKlCmoVasWNDQ0eNc/ffo0mjdvjho1akBdXZ13fQ8PD9SpUwepqamC6H/58gUjR44ULCYuLS0NLi4uePfunSD6ABASEgKFQiGYPoPB+H4w4+gHIS0tDQCgrq6OxYsX8771tWvXLqipqQHg3zhKSEjAli1buN/5No5evXqFnTt3cr/z7Tl68OABp29pacmrNgDcuHEDS5YsAQD89ttvvOufP38ew4cPR1paGnr27Mm7/tGjR9GyZUsEBQVh3LhxvOvv3bsX27dvx+PHjzF8+HBetYkIGzduRLly5XDr1i3Y29vzqi+TybB161Y0bNgQp0+f5t1wT05Oxj///IOxY8fi1atXvGoDWX+7gYGBOHbsGORyOe/6KSkpePv2LYKDg3nXVpKZmSnI3H92hApvUCJkYoxCoYBMJhNMvzAw4+gHoWnTppgxYwYcHBzwzz//8Kqdnp6OI0eOwMnJCbt370blypV51ReLxXj16hXGjh2LTZs28f6hrlq1Kh48eIDBgwdj+PDhvBtHzZs3x5EjR1CnTh3UqVOHV20A6NixI9avX48OHTqgQYMGvOs7OjpizJgxqFGjBjp37sy7ft++fdGtWzd0794d5cqV412/ffv2MDAwwNKlS3n3SolEItjY2MDAwAB//fUXr9oAIJFIEB4ejjJlymDmzJm86+vq6mLLli1o2bIlatSowbt+bGwsrl69ivr163M3T3xy9uxZzJkzB0ZGRrxrA8C4ceOwdetWQbbaFQoF+vbti5cvX/KuDQCBgYFwc3MTxMj4/Pkz9PT0EBISwrs2ABw/fhx169YVRBsAxo4di1mzZgmmXxhExArsqExiYiIMDAyQkJAAfX19XjSVdz4BAQGoVKmSIF9UvwLKbQu+79DlcjkCAgJgbW2NUqVK8aqtRKFQICEhAaVLlxZEPzIyEhYWFoJoA1mZfNbW1rzryuVyPHv2DHXq1BFkyw4AIiIiBPEKAll30Onp6TAwMBBEPz4+XrDPjFwux+vXr1G9enVBvnMyMjKgrq4uaBKFUJ8ZIGt9hPoulslkeP/+vWDf948ePULDhg0FMRxTUlIQFBSEmjVr8q4NAJ8+fYJYLBbk+6aw12/h0nIYKqGmpsZcw4VAqC9ZJTo6OoJpi8ViwS5yAAQ1jAAI8kWlpFSpUoK+t0IZRgCgqakJTU1NwfSF/MwAECQGTomQ6wJAUMMIgKA3qUIYLdlp1KiRYNq6urqCGUYAULZsWcG0CwvbVmMwGAwGg8HIBjOOGAwGg8FgMLLBjCMGg8FgMBiMbPxSxlH58uUhEoly/Zs8eTIAYMSIEbkea9y4cQnPmsFgMBgMxo/ELxWQ/eTJkxxBza9evUKHDh3Qr18/7piDgwMOHDjA/S5kMCKDwWAwGIyfj1/KODI1Nc3xu7OzMypWrIhWrVpxxzQ1NVXO6snIyMhRiyIxMRFAVponXxlm3t7eqFOnDhQKhaBZa0FBQbC1tRU0U0L
|
||
|
|
"text/plain": [
|
||
|
|
"<Figure size 640x480 with 1 Axes>"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
"metadata": {},
|
||
|
|
"output_type": "display_data"
|
||
|
|
}
|
||
|
|
],
|
||
|
|
"source": [
|
||
|
|
"z = np.arange(-100,100,10)\n",
|
||
|
|
"print(np.size(z))\n",
|
||
|
|
"r = np.arange(1e-15,100,5)\n",
|
||
|
|
"\n",
|
||
|
|
"B, B_tot = Coil.B_multiple_3d(I, r,z,raster = 4)\n",
|
||
|
|
"\n",
|
||
|
|
"z_m, r_m = np.meshgrid(z,r)\n",
|
||
|
|
"\n",
|
||
|
|
"#plt.figure(figsize=(16,10))\n",
|
||
|
|
"plt.quiver(r_m,z_m,B[:,:,1],B[:,:,0])\n",
|
||
|
|
"plt.xlabel(\"radius r [m]\")\n",
|
||
|
|
"plt.ylabel(\"z-axis [m]\")\n",
|
||
|
|
"plt.show()"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": 276,
|
||
|
|
"id": "28dcdd21-f4b6-49e6-a935-4012a941f186",
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [],
|
||
|
|
"source": [
|
||
|
|
"z = np.arange(-100,100,1)\n",
|
||
|
|
"r = np.arange(1e-3,100,1)\n",
|
||
|
|
"\n",
|
||
|
|
"B, B_tot = Coil.B_multiple_3d(I, r,z,raster = 2)\n",
|
||
|
|
"\n",
|
||
|
|
"z_m, r_m = np.meshgrid(z,r)"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": 277,
|
||
|
|
"id": "ad02e8cd-e92a-4d27-845d-cc7e22d6e7e5",
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [],
|
||
|
|
"source": [
|
||
|
|
"x = np.concatenate((-np.flip(r),r))\n",
|
||
|
|
"B_tot_x = np.concatenate((np.flip(B_tot[:,len(z)//2]),B_tot[:,len(z)//2]))\n",
|
||
|
|
"B_tot_z = B_tot[0,:]\n"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": 278,
|
||
|
|
"id": "52475c2c-e2f0-4615-a5f4-9f0a86155a3e",
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [
|
||
|
|
{
|
||
|
|
"data": {
|
||
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAjMAAAGwCAYAAABcnuQpAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAA9hAAAPYQGoP6dpAABgt0lEQVR4nO3deXxTVfo/8M9N0qRbkm50XygtWymbbIrKKiAqDjoq6gzujoqoiMvI1w0YGdRxFBVRUQZxHBVnUH7uCMqiIvu+lW60dKN0Tde0Tc7vjzTB0oWkTXuzfN6vV1+Qm5ub5/SmyZNznnuOJIQQICIiInJTCrkDICIiIuoKJjNERETk1pjMEBERkVtjMkNERERujckMERERuTUmM0REROTWmMwQERGRW1PJHUB3M5vNKCgogFarhSRJcodDREREdhBCoKqqCtHR0VAoOu578fhkpqCgAHFxcXKHQURERJ1w+vRpxMbGdriPxyczWq0WgOWXodPpunw8k8mEzMxMJCUlQalUdvl4roht9Axso2dgGz0D2+g4g8GAuLg42+d4Rzw+mbEOLel0OqclM4GBgdDpdB79gmQb3R/b6BnYRs/ANnaePSUiLAAmIiIit8ZkhoiIiNwakxkiIiJya0xmiIiIyK0xmSEiIiK3xmSGiIiI3BqTGSIiInJrTGaIiIjIrTGZISIiIrfGZIaIiIjcGpMZIiIicmtMZoiIiMitMZkhIrdkNguU1zagtLYJxiaz3OEQkYw8ftVsIvIMZ6uM+O5IIXZll2F/bgUKK+tgFtZ7c6DzVWFwrB4jEkIwZWAEUmN0dq22S0Tuj8kMEbm0nVmleP+XbPx0ohimc9mLjQRAADDUN+HXjFL8mlGKN35Mx4BILf58cQJuGhkHtYqd0ESejMkMEbmk9DNVePG7E/jxRLFt27C4IEweEI4RvYORHB4IrVqJ7KwMRMT2xpmqRuzLLcdvmaXYePwMThRV4Zn1R/Dez1l4Ylp/XD04ij01RB6KyQwRuRSzWWDVL9n4x4Y0NJjMUCokzBoVhzvH9kbfCG2LfU0mExSShCB/NUK1fkiJ1uHPFyegsrYR6/blYcWWTOSU1mLux/vxZUoBXvzjEIQEqGVqGRF1FyYzROQyymoa8NAn+/BrRikAYGL/XnjmmhQk9Qp06Dh6fx/cdVkiZo2Kw8ptWVixJQM/HDuD/ae34c1bhuPiPqHdET4RyYQDyUTkEk6V1OD6Fb/i14xS+PoosOS6VPzrjlEOJzK/F6BR4dEp/fDFnEuRHB6Is1VGzF61E//vQL4TIyciuTGZISLZHcqrwPVvb8ep0lrEBPnhy7mX4U9jEpxW45Iao8dXcy/DVYMj0WgSeOTTA3hna6ZTjk1E8mMyQ0SyOl5owOxVu1BW04DBMXp88eBY9DuvNsYZ/NRKLL/lItx7eSIA4MXvTmDVL9lOfx4i6nlMZohINllnqzF71U5U1jVieHwQPvnLxQjX+nbb8ykUEp6+OgWPXtEPAPC3r4/hk1253fZ8RNQzmMwQkSzKahpw++pdKKluQEqUDh/cMRqBmp65JuHhycm4b1wfAMD/fXEYm393+TcRuR8mM0TU4xpNZjz4n304XVaH+BB//Pvu0dD7+/TY80uShKemD8CskXEQAnj40/3IPFvdY89PRM7FZIaIetySb47jt6xS+KuVeO+2kQgN1PR4DJIk4W8zUzEyIRhV9U2498M9MNQ39ngcRNR1TGaIqEdtOFqED7afAgC8etMw9I90frGvvdQqBd7+8whE6X2RdbYGz60/IlssRNR5TGaIqMcUG+rx1LpDAID7xvXBlamRMkcE9NJq8NafLoJSIWH9gQLOQUPkhpjMEFGPEELgif8dQnltI1KidJg/tZ/cIdlcFB+MhyYlAwCeWX8E+RV1MkdERI5gMkNEPWLt7tPYevIsNCoFXr95GDQqpdwhtTB3YjKGxQWhqr4JT607BCFar9BNRK6JyQwRdbvSaiOWfncCAPDEtP6tFox0BSqlAq/NGga1SoGf00vw9aFCuUMiIjsxmSGibvf3b0+gss4yvHTH2N5yh9OuxLAAzJ1oGW5a/PUxXt1E5CaYzBBRt9qZVYp1+/IgScCS61KhUrr228594/ugT1gAzlYZ8c8NaXKHQ0R2cO13FSJya2azwKKvjgEAbhkdj+HxwTJHdGEalRJ/m5kKAPj3jhykn6mSOSIiuhAmM0TUbf7fwXwcKzRA66vCE1P7yx2O3S5NDsO0QREwC+Cl79k7Q+TqmMwQUbeobzThlQ0nAQAPTEhCcIBa5ogc8+SVA6BUSNh0/Ax2ZpXKHQ4RdYDJDBF1i3//loP8ijpE6X1x16WJcofjsKRegbh5VBwA4O/fneCl2kQujMkMETldtbEJb23JAAA8OqUffH1ca04Zez1yRV/4q5U4eLoCPxw7I3c4RNQOWZOZbdu2YcaMGYiOjoYkSVi/fn27+953332QJAnLli3rsfiIqHM+2pGDitpG9AkLwB8vipU7nE4L1/raLiVf/lMGe2eIXJSsyUxNTQ2GDh2K5cuXd7jf+vXrsXPnTkRHR/dQZETUWXUNJrz/cxYAYM7EZCgVkswRdc3dlyXCz0eJw/mV2HLyrNzhEFEbVHI++fTp0zF9+vQO98nPz8fcuXOxYcMGXH311Rc8ptFohNFotN02GAwAAJPJBJPJ1LWAm49jNpudcixXxTZ6Brna+J+dp1BS3YC4YD9cMziiW5+/J9oY5KfCraPjsOrXU3jjx3RcnhQCSeq5BI2vVc/ANnbuePaSNZm5ELPZjNmzZ+OJJ57AoEGD7HrM0qVLsWjRolbbMzMzERgY6JSYysrKkJGRAYXCM0uO2EbPIEcbG0wCb2/OAQBcPzAQp7Iyu/X5eqqNk2IEPlRI2J9bgf/9fAjDovy77bnOx9eqZ2AbHVddXW33vi6dzLz00ktQqVR4+OGH7X7MggULMH/+fNttg8GAuLg4JCUlQafTdTkmk8mEjIwMJCcnQ6l0z6LGC2EbPYMcbfzf3jyU1poQqdPgvmnDoVF175t2T7WxL4CbcwT+vTMX32Y14sZxfbvtuc7H16pnYBsdZx1ZsYfLJjN79+7F66+/jn379jnUpavRaKDRaFptVyqVTnsBKRQKpx7PFbGNnqEn2yiEwOrtll6ZOy9NhL/Gp9ufE+i5Nt4zrg8+2pWLbeklyC6tRXJ4zy2WydeqZ2AbHePIMVy2r+vnn39GcXEx4uPjoVKpoFKpkJOTg8ceewy9e/eWOzwiOs9vWaU4UVQFPx8lbh4VL3c4TpcQGoArBkYAAFb/ekreYIioBZdNZmbPno1Dhw7hwIEDtp/o6Gg88cQT2LBhg9zhEdF5/vXLKQDADSNioffvmV6Znmad/G/dvjxU1DbIHA0RWck6zFRdXY2MjAzb7ezsbBw4cAAhISGIj49HaGhoi/19fHwQGRmJ/v3dZ40XIm9wqqQGP56wTCp3x6W95Q2mG13cJwQDo3Q4XmjAJ7tO44EJSXKHRESQuWdmz549GD58OIYPHw4AmD9/PoYPH47nnntOzrCIyEEf7ciBEMCE/r2Q1KvrVw26KkmScFdzsvbRjhyYzJxEj8gVyNozM2HCBIdm1Dx16lT3BUNEnWJsMmHdvjwAwG2XJMgcTfebMTQaL3xzHPkVddiWfhYT+4fLHRKR13PZmhkicg8bjp5BeW0jovS+GN/P8z/YfX2UuG54DADg0125MkdDRACTGSLqok92Wj7QbxoZ5/ZLF9jrltGWq7V+PF6MYkO9zNEQEZMZIuq07JIa/JZVCkkCbhoVJ3c4PaZ/pBYXxQehySzw3715codD5PWYzBBRp32629IrM75fL8QE+ckcTc+6ubl3Zu3u0zCzEJhIVkxmiKhTTGaBz/flA4BHTpJ3IdcMiYJWo0JuWS12nyqTOxwir8Zkhog6ZXtmCc5WGRHk74NJAzy/8Pd8/moVpg+OBACsP5AvczRE3o3JDBF1yhf7LR/g1wyJgrqbF5R0VTObr2r6+lAh6htNMkdD5L288x2IiLqktqEJG44UAYDtMmVvdHFiKKL0vqiqb8LmE8Vyh0P
|
||
|
|
"text/plain": [
|
||
|
|
"<Figure size 640x480 with 1 Axes>"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
"metadata": {},
|
||
|
|
"output_type": "display_data"
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"data": {
|
||
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAkAAAAG1CAYAAAARLUsBAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAA9hAAAPYQGoP6dpAAByvElEQVR4nO3dd3gU1foH8O/sJtn0TipppAAhEEInSBMBQVFEBbEAigUFFZGLl3tFwQK2q9jLFQHlp6AXVOyAEooUCRA6Ib2H9J5skt35/ZFkNSbAbtjN7M5+P8+T52FnZ2bfl9lk3z3nzDmCKIoiiIiIiKyIQuoAiIiIiLobCyAiIiKyOiyAiIiIyOqwACIiIiKrwwKIiIiIrA4LICIiIrI6LICIiIjI6rAAIiIiIqvDAoiIiIisDgsgIiIisjqSFkBr1qzB0KFD4eLiAh8fH0yfPh3Jycnt9hFFEStXrkRAQAAcHBwwbtw4nDlz5orn3rp1K6Kjo6FSqRAdHY2vv/7aVGkQERGRhZG0ANqzZw8WLlyIQ4cOYefOnWhubsakSZNQW1ur2+eVV17B66+/jnfeeQdHjhyBn58fJk6ciOrq6kue9+DBg5g1axbuuecenDhxAvfccw9mzpyJw4cPd0daREREZOYEc1oMtbi4GD4+PtizZw/GjBkDURQREBCAxYsX46mnngIAqNVq+Pr64uWXX8ZDDz3U6XlmzZqFqqoq/PTTT7pt119/PTw8PPDFF19cMQ6tVov8/Hy4uLhAEATjJEdEREQmJYoiqqurERAQAIXi8m08Nt0Uk14qKysBAJ6engCAjIwMFBYWYtKkSbp9VCoVxo4diwMHDlyyADp48CCeeOKJdtsmT56MtWvXdrq/Wq2GWq3WPc7Ly0N0dPTVpEJEREQSycnJQc+ePS+7j9kUQKIoYsmSJbjmmmsQExMDACgsLAQA+Pr6ttvX19cXWVlZlzxXYWFhp8e0ne/v1qxZg1WrVnXYvnv3bjg7OxuUx5VotVqUl5fDw8PjitWppZJ7jnLPD2COciD3/AD55yj3/ADj51hTU4Px48fDxcXlivuaTQG0aNEinDx5Evv37+/w3N+7oURRvGLXlCHHLF++HEuWLNE9rqqqQlBQEGJjY+Hq6qpvCnrRaDRITU1FREQElEqlUc9tLuSeo9zzA5ijHMg9P0D+Oco9P8D4OVZVVQHoWAN0xiwKoEcffRTbt2/H3r172zVZ+fn5AWhp0fH399dtLyoq6tDC81d+fn4dWnsud4xKpYJKpeqwXalUmuRNp1AoTHZucyH3HOWeH8Ac5UDu+QHyz1Hu+QHGzdGQc0japiaKIhYtWoRt27bht99+Q1hYWLvnw8LC4Ofnh507d+q2NTY2Ys+ePYiPj7/keUeOHNnuGADYsWPHZY8hIiIi6yFpC9DChQvx+eef49tvv4WLi4uu1cbNzQ0ODg4QBAGLFy/G6tWrERkZicjISKxevRqOjo648847deeZM2cOAgMDsWbNGgDA448/jjFjxuDll1/GzTffjG+//Ra7du3qtHuNiIiIrI+kBdD7778PABg3bly77evXr8e8efMAAMuWLUN9fT0eeeQRlJeXY/jw4dixY0e7AU7Z2dntBk/Fx8dj8+bNePrpp7FixQqEh4djy5YtGD58uMlzIiIiIvMnaQGkzxREgiBg5cqVWLly5SX3SUhI6LDttttuw2233XYV0REREZFcyfO+OiIiIqLLYAFEREREVocFEBEREVkdFkBERERkdVgAERERkdVhAURERERWhwUQERERWR0WQER0RSU1ahzNKkNxtVrqUKgT5bWNOJpVjoLKeqlDIbIYZrEYKhGZn/omLV79JRnfnypEbvmfH6xeTnaYFhuAR6+NgJdzx0WEqXtUNzThwz3p+N/RXBRWNei2+7qqMLmfHx6/NlzC6IjMHwsgIupg9/ki/OubHBTVNgMABAHwcVGhqFqN0tpGbDiQia1Hc/H4dZGYf00YBEGQOGLr8r+juVjz4zmU1jbqtvm6qlBcrcbFKjU+PZiF708W4P5B7oiMlDBQIjPGAoiI2vnySA6WbT0JAAh0t8e/pkZjdJQ3XO1tUd+owaGMUrz2SzLO5FfhhR/OIeViDV68JQY2Svaom5ooinj1l2S8l5AGAOjl7YSlk3tjTFQPOKtsUNfYjMPpZVj94zmkFNXglX1F0Nin4dEJURJHTmR+WAARkc43x/Pw1LaW4mdKlCtemT0MLg5/dnM52CkxvrcPxkb2wKcHM/Hc92exJTEHxTVqvHfXINjbKqUKXfY0WhFPbT2J/x3NBQA8dm0EHp0QCdu/FJ6OdjYY38cHoyK8sXZXMt5LSMd/dqbAwc4G94/uJVXoRGaJX9mICACwL6UYT351AqII3DUsCI+N9IajXeffkRQKAfNGheGDuwdDZaPAb+eLsHzbKb0WOKaueeXn8/jf0VwoFQJevrU/lkzq3a74+Ss7GwWenBiFewZ6AABe+OEcvjme153hEpk9FkBEhIq6Rjz55QlotCJmxAVi5bRovcb1TOrnh3Vzh0KpEPD18Ty8vyetG6K1Pl8l5uDDvekAgNdnxmLW0GC9jrsz1gMPjA4DADz9zWnklteZLEYiS8MCiIjw9DenUVStRq8eTnjxlv5QKPQf1HxNpDdW3tQPAPDKz8n49dxFU4VplY5nl+NfX58CADx6bQRuHhio97GCIOAfk6IwOMQDNepmLP3qBLRattIRASyAiKzedyfy8f3JAigVAt6YORAOdoaP47lnRAjuGRECAHhq60mU1nC+IGOob9RgyZcn0KQRMbmfL564zvDBzEqFgNdnxsLRTolD6WX45PcME0RKZHlYABFZsfpGDV784RwAYNH4CMQGuXf5XP++oS+ifJ1RUtOIFd+e5nggI3j55/PIKKmFr6sKr9waa1DL3F+FeDnh6RuiAQBrd6WwQCUCCyAiq7bhQCYKqxoQ6O6AR8Zf3cR59rZK/Of2gbBRCPjxVCG2n8g3UpTW6WBaKTYcyAQAvHzrALg52l7V+e4YGoSYQFfUqJvxzu5UI0RIZNlYABFZqfLaRryX0PJB+OSkKKhsrv4W9v493bBwfAQA4LnvzqKyvumqz2mNGpu1unE/s4cFY1xvn6s+p0Ih4J/X9wUAbDqUhexSDogm68YCiMhKvZeQiuqGZvTxczFoYO2VLBwfgfAeTiitbcSbu1KMdl5rsv73DGSU1MLbWYV/Te1jtPNeE+mN0ZHeaNKI+M/OZKOdl8gSsQAiskKlNS3LJQDAU1P6QNnFsSWdsbNR4NlpLXeFbTyYiZSL1UY7tzW4WNWAt35tKRz/OaUPXOyvruvr7566vqWg+jYpH2nFNUY9N5ElYQFEZIU2HcqGulmLAT3dMC6qh9HPPyaqByZG+0KjFbHyuzMcEG2Al386j9pGDeKC3TEjzngtc21iAt1wXV9fAMAn+3lHGFkvFkBEVqahSYPPDmUCgEkXMl1xQzTsbBT4PbUUey4Um+Q15OZ0XiW2tc7YvHJavy7f9XUl97dOjvi/o7ko+8uCqkTWhAUQkZX5NikPJTWNCHCzx9T+/iZ7nWAvR8xpnRvo1V+SOQGfHl7b0TIu56bYgKuakuBKhod5on+gG9TNWvzfoSyTvQ6ROWMBRGRFRFHEx/tauj3uHRV2ybWkjOWR8RFwVtngTH4VfjhVYNLXsnSH00uRkFwMG4WAJRNNu3q7IAi6VqCNB7PQ0KQx6esRmSMWQERW5EBaKVKKauBkp8SsYUEmfz1PJzs80LoK+es7L6BJozX5a1oiURTxyi8trT+zhgYh1NvJ5K85tb8/AtzsUVKjxk+nWZyS9WEBRGRFvkzMAQBMjwuEq5HvLrqU+aPD4OVkh4ySWnzNFck7tedCMY5mlcPeVoHHJkR2y2vaKhW4Y1jLoqpfHsntltckMicsgIisRGVdE346XQgAmDnE9K0/bZxVNnhobEsr0PsJadBwLFA7oijizdbb3u8eHgJfV/tue+1bB/eEIAAH00uRVVrbba9LZA5YABFZie0n8tDYrEUfPxcM6OnWra991/AQuDvaIqOkFj9yLFA7+1N
|
||
|
|
"text/plain": [
|
||
|
|
"<Figure size 640x480 with 1 Axes>"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
"metadata": {},
|
||
|
|
"output_type": "display_data"
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"data": {
|
||
|
|
"text/plain": [
|
||
|
|
"'\\nplt.plot(x,np.gradient(B_tot_x,x))\\nplt.xlabel(\"radius r [mm]\")\\nplt.ylabel(\"total field B_tot [G]\")\\n#plt.xlim(0,0.01)\\nplt.show()\\nplt.plot(z,np.gradient(B_tot[0,:],z))\\nplt.xlabel(\"z_axis [mm]\")\\nplt.ylabel(\"total field B_tot [G]\")\\nplt.show()\\n'"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
"execution_count": 278,
|
||
|
|
"metadata": {},
|
||
|
|
"output_type": "execute_result"
|
||
|
|
}
|
||
|
|
],
|
||
|
|
"source": [
|
||
|
|
"plt.plot(x,B_tot_x)\n",
|
||
|
|
"plt.xlabel(\"radius r [mm]\")\n",
|
||
|
|
"plt.ylabel(\"total field B_tot [G]\")\n",
|
||
|
|
"#plt.xlim(-0.1,0.1)\n",
|
||
|
|
"plt.show()\n",
|
||
|
|
"plt.plot(z,B_tot_z)\n",
|
||
|
|
"plt.xlabel(\"z_axis [mm]\")\n",
|
||
|
|
"plt.ylabel(\"total field B_tot [G]\")\n",
|
||
|
|
"plt.show()\n",
|
||
|
|
"\n",
|
||
|
|
"\"\"\"\n",
|
||
|
|
"plt.plot(x,np.gradient(B_tot_x,x))\n",
|
||
|
|
"plt.xlabel(\"radius r [mm]\")\n",
|
||
|
|
"plt.ylabel(\"total field B_tot [G]\")\n",
|
||
|
|
"#plt.xlim(0,0.01)\n",
|
||
|
|
"plt.show()\n",
|
||
|
|
"plt.plot(z,np.gradient(B_tot[0,:],z))\n",
|
||
|
|
"plt.xlabel(\"z_axis [mm]\")\n",
|
||
|
|
"plt.ylabel(\"total field B_tot [G]\")\n",
|
||
|
|
"plt.show()\n",
|
||
|
|
"\"\"\""
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "markdown",
|
||
|
|
"id": "37c71478-cae2-40ea-8712-b2f751c6c54d",
|
||
|
|
"metadata": {},
|
||
|
|
"source": [
|
||
|
|
"## Fit harmonic function"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": 279,
|
||
|
|
"id": "443f51a2-a502-4d4a-8692-15610c193cc7",
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [],
|
||
|
|
"source": [
|
||
|
|
"def func(x,a,b):\n",
|
||
|
|
" return a*x**2+b\n",
|
||
|
|
"\n",
|
||
|
|
"x_SI = 1e-3 * x\n",
|
||
|
|
"z_SI = 1e-3 * z\n",
|
||
|
|
"\n",
|
||
|
|
"B_tot_x_SI = 1e-4*B_tot_x\n",
|
||
|
|
"B_tot_z_SI = 1e-4*B_tot_z\n"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": 280,
|
||
|
|
"id": "59145c9b-aac3-4bea-ba79-b2a3992d05f6",
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [],
|
||
|
|
"source": [
|
||
|
|
"nr_points = 20\n",
|
||
|
|
"a = 100-nr_points//2\n",
|
||
|
|
"b = 100+ nr_points//2\n",
|
||
|
|
"popt_x, pcov = curve_fit(func,x_SI[a:b],B_tot_x_SI[a:b])"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": 286,
|
||
|
|
"id": "49c6fb78-ba74-41e3-9a2a-c4dca7109643",
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [
|
||
|
|
{
|
||
|
|
"data": {
|
||
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAkAAAAG2CAYAAACXuTmvAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAA9hAAAPYQGoP6dpAACHfUlEQVR4nOzdd3hT1RvA8W+S7t1SuqC0Bcoue8heskUQBy6GuBBQEf2hqCiIggtERUUcgCKCylRBhrL33qOF0jJaChS6V5L7++NCsXaQQNOkzft5nj7e3HvOzRtOa9+ee4ZGURQFIYQQQgg7orV2AEIIIYQQZU0SICGEEELYHUmAhBBCCGF3JAESQgghhN2RBEgIIYQQdkcSICGEEELYHUmAhBBCCGF3JAESQgghhN2RBEgIIYQQdkcSICGEEELYHasmQFOmTKFFixZ4enoSEBBA//79OXHiRIEyiqIwYcIEQkJCcHV1pVOnThw5cuSW9160aBH16tXD2dmZevXqsWTJEkt9DCGEEEKUM1ZNgDZs2MDIkSPZvn07a9asQa/X0717dzIyMvLLfPjhh0ybNo0ZM2awa9cugoKC6NatG2lpacXed9u2bQwcOJBBgwZx4MABBg0axEMPPcSOHTvK4mMJIYQQwsZpbGkz1EuXLhEQEMCGDRvo0KEDiqIQEhLC6NGjefXVVwHIyckhMDCQDz74gGeffbbI+wwcOJDU1FRWrlyZf65nz574+vry888/l8lnEUIIIYTtcrB2AP+WkpICgJ+fHwCxsbEkJibSvXv3/DLOzs507NiRrVu3FpsAbdu2jZdeeqnAuR49ejB9+vQiy+fk5JCTk5P/2mg0kpycTKVKldBoNHfykYQQQghRRhRFIS0tjZCQELTakh9y2UwCpCgKY8aMoV27djRo0ACAxMREAAIDAwuUDQwMJC4urth7JSYmFlnnxv3+a8qUKUycOPFOwhdCCCGEjTh79ixVq1YtsYzNJECjRo3i4MGDbN68udC1//bCKIpyy54Zc+qMGzeOMWPG5L9OSUmhWrVqnDlzBi8vL1M/gs0xGAycPn2a6tWro9PprB2OXZO2sB3SFrZF2sN2VIS2SE1NJTw8HE9Pz1uWtYkE6Pnnn2f58uVs3LixQMYWFBQEqD06wcHB+eeTkpIK9fD8W1BQUKHenpLqODs74+zsXOi8r69vuU+AvLy88PX1LbffzBWFtIXtkLawLdIetqMitMWNuE0ZvmLVWWCKojBq1CgWL17MP//8Q0RERIHrERERBAUFsWbNmvxzubm5bNiwgTZt2hR739atWxeoA7B69eoS6wghhBDCfli1B2jkyJHMnz+fZcuW4enpmd9r4+3tjaurKxqNhtGjRzN58mQiIyOJjIxk8uTJuLm58eijj+bfZ/DgwVSpUoUpU6YA8OKLL9KhQwc++OAD+vXrx7Jly1i7dm2Rj9eEEEIIYX+smgB99dVXAHTq1KnA+dmzZzN06FAAxo4dS1ZWFiNGjODq1au0atWK1atXF3i+Fx8fX2C0d5s2bViwYAFvvvkm48ePp0aNGixcuJBWrVpZ/DMJIYQQwvbZ1DpAtiI1NRVvb29SUlLK/Rig6OhoIiMjy+3z3IpC2sJ2SFtYh9FoJDc3t9B5g8FAXFwcYWFh0h5WVl7awsnJqdgp7ub8/raJQdBCCCEqrtzcXGJjYzEajYWuKYqCXq8nLi5O1l2zsvLSFlqtloiICJycnO7oPpIACSGEsBhFUUhISECn0xEaGlroL3dFUcjJycHZ2dmmf+nag/LQFkajkQsXLpCQkEC1atXuKE5JgIQQQliMXq8nMzOTkJAQ3NzcCl2/MQrDxcXFZn/p2ovy0haVK1fmwoUL6PV6HB0db/s+Vp0GL4QQomIzGAwAd/y4Qogbbnwv3fjeul2SAAkhhLA4W+5REOVLaX0vSQIkhBBCCLsjCZAQQgjxH506dWL06NHWDqNUnTlzBo1Gw/79+2/7Hlu2bCEqKgpHR0f69+/P+vXr0Wg0XLt2rdTiLCsyCFoIIYSwA6GhoSQkJODv73/b9xgzZgyNGzdm5cqVeHh44ObmRkJCAt7e3gDMmTOH0aNHl4uESHqAhBBCiDJgMBiKXAuprOh0OoKCgnBwuP2+j1OnTtGlSxeqVq2Kj48PTk5OBAUFlcsxXpIACSGEEEUwGo2MHTsWPz8/goKCmDBhQoHr06ZNIyoqCnd3d0JDQxkxYgTp6en51+fMmYOPjw9//PEH9erVw9nZmbi4OMLDw3n33XcZPHgwHh4ehIWFsWzZMi5dukS/fv3w8PAgKiqK3bt3F3i/RYsWUb9+fZydnQkPD2fq1KkFroeHhzN58mSGDRuGp6cn1apVY9asWfnXi3oEduTIEfr06YOXlxdeXl7cfffdnDp1qtC/xY26V65cYdiwYWg0GubMmVPgEdj69et54oknSElJQaPRoNFoCv2b2RJJgIQQQpQZRVHIzNX/58tQxLnS/zJ356e5c+fi7u7Ojh07+PDDD3nnnXdYs2ZN/nWtVstnn33G4cOHmTt3Lv/88w9jx44tcI/MzEymTJnCt99+y5EjRwgICADgk08+oW3btuzbt48+ffowaNAgBg8ezOOPP87evXupWbMmgwcPzo95z549PPTQQzz88MMcOnSICRMmMH78eObMmVPg/aZOnUrz5s3Zt28fI0aM4LnnnuP48eNFfr7z58/ToUMHXFxc+Oeff9i9ezeDBw9Gr9cXKnvj8ZmXlxfTp08nISGBgQMHFijTpk0bpk+fjpeXFwkJCSQkJPDKK6+Y9W9elmQMkBBCiDKTlWeg3lurrPLeR9/pgZuT6b/2GjZsyNtvvw1AZGQkM2bM4O+//6Zbt24ABQZJR0REMGnSJJ577jm+/PLL/PN5eXl8+eWXNGrUqMC9e/fuzbPPPgvAW2+9xVdffUWLFi148MEHAXj11Vdp3bo1Fy9eJCgoiGnTptG1a1fGjx8PQK1atTh69CgfffRR/ubhN+47YsSI/Ht88sknrF+/njp16hT6fF988QXe3t4sWLAAR0dHFEWhWrVquLi4FCp74/GZRqPB29uboKCgQmWcnJzw9vZGo9EUed3WSA+QEEIIUYSGDRsWeB0cHExSUlL+63Xr1tGtWzeqVKmCp6cngwcP5sqVK2RkZOSXcXJyKnSf/947MDAQgKioqELnbrzfsWPHaNu2bYF7tG3blujo6AILAv77vjcSkX/H/G/79++nffv2d7SacnkmPUBCCCHKjKujjqPv9Mh/rSgK2dk5uLhYfv8pV0fzdjj/b2Kg0WjyBzHHxcXRu3dvhg8fzqRJk/Dz82Pz5s08+eST5OXl3XxPV9ciP9e/733jelHnbryfoiiF7lPUI72SYv4vV1fXIs/bC0mAhBBClBmNRlPgMZSiKGiNelycHMrVTKLdu3ej1+uZOnVq/gavv/zyi8Xer169emzevLnAua1bt1KrVi10OvMSuxsaNmzI3LlzycvLK7VeICcnpzveoqKsyCMwIYQQwkw1atRAr9fz+eefc/r0aX788Udmzpxpsfd7+eWX+fvvv5k0aRInT55k7ty5zJgx444GGY8aNYrU1FQefvhhdu/eTXR0NPPnz+fEiRO3fc/w8HDS09P5+++/uXz5MpmZmbd9L0uTBEgIIYQwU+PGjZk2bRoffPABDRo04KeffmLKlCkWe7+mTZvyyy+/sGDBAho0aMBbb73FO++8U2AAtLkqVarEP//8Q3p6Oh07dqR58+bMnj37jnqD2rRpw/Dhwxk4cCCVK1fmww8/vO17WZpGMXdeoB1ITU3F29ublJQUvLy8rB3ObTMYDERHRxMZGXnbXaSidEhb2A5pi7KVnZ1NbGwsERERRc4uUscAZePi4lKuHoFVROWlLUr6njLn97f0AAkhhBDC7kgCJIQQQgi7IwmQEEIIIeyOJEBCCCGEsDuSAAkhhBDC7kgCJIQQQgi7IwmQEEIIIeyOJEBCCCGEsDuSAAkhhBDC7kgCJIQQQgi7I7vBCyGEEBVMeno6np6e1g6jVFhqxy7pARJCCCEqmNWrV3PgwAEURTH5y2g0kpmZidF
|
||
|
|
"text/plain": [
|
||
|
|
"<Figure size 640x480 with 1 Axes>"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
"metadata": {},
|
||
|
|
"output_type": "display_data"
|
||
|
|
}
|
||
|
|
],
|
||
|
|
"source": [
|
||
|
|
"plt.plot(1e3*x_SI,1e4*func(x_SI,popt_x[0],popt_x[1]),label = \"harmonic fit\")\n",
|
||
|
|
"plt.plot(1e3*x_SI,1e4*B_tot_x_SI, label = r\"$ B_{tot} = \\sqrt{r^2 + z^2} $\" )\n",
|
||
|
|
"#plt.ylim(0,8)\n",
|
||
|
|
"#plt.xlim(-0.01,0.01)\n",
|
||
|
|
"plt.xlabel(\"x-axis [mm]\")\n",
|
||
|
|
"plt.xlim(-50,50)\n",
|
||
|
|
"plt.ylim(0,20)\n",
|
||
|
|
"plt.ylabel(\"B_tot [G]\")\n",
|
||
|
|
"plt.legend()\n",
|
||
|
|
"plt.show()"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": 282,
|
||
|
|
"id": "77c981bb-02be-4de9-8990-849b6a7fbee1",
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [],
|
||
|
|
"source": [
|
||
|
|
"nr_points = 20\n",
|
||
|
|
"a = 100-nr_points//2\n",
|
||
|
|
"b = 100+ nr_points//2\n",
|
||
|
|
"popt_z, pcov = curve_fit(func,z_SI[a:b],B_tot_z_SI[a:b])"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": 285,
|
||
|
|
"id": "f1f7cdd7-f13c-41e2-be5e-a375dc5e4d06",
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [
|
||
|
|
{
|
||
|
|
"data": {
|
||
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAjMAAAG2CAYAAACKxwc0AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAA9hAAAPYQGoP6dpAABxyElEQVR4nO3dd3gUVd/G8e9ueg8BkhBKAoQemiAICIQqgr03ECuIqFgeFHlUbKCoyKvYHwVEEXtHBEWa9N4hlFATAgRSSN+d94+BSKQlkGR2k/tzXXvpzs7O3svJJr89c+Ycm2EYBiIiIiJuym51ABEREZELoWJGRERE3JqKGREREXFrKmZERETEramYEREREbemYkZERETcmooZERERcWsqZkRERMStqZgRERERt6ZiRkRERNyapcXMqFGjsNlsRW6RkZGFjxuGwahRo4iKisLPz4/4+Hg2bNhgYWIRERFxNZb3zDRr1oykpKTC27p16wofGzt2LOPGjWPChAksW7aMyMhIevXqRUZGhoWJRURExJVYXsx4enoSGRlZeKtevTpg9sqMHz+ekSNHct111xEXF8fkyZPJyspi6tSpFqcWERERV+FpdYCEhASioqLw8fGhffv2jB49mnr16rFz506Sk5Pp3bt34b4+Pj507dqVhQsXMmjQoNMeLzc3l9zc3ML7TqeT1NRUqlatis1mK/P3IyIiIhfOMAwyMjKIiorCbj9734ulxUz79u359NNPadiwIQcOHOCll16iY8eObNiwgeTkZAAiIiKKPCciIoJdu3ad8Zhjxozh+eefL9PcIiIiUj727NlDrVq1zrqPzTAMo5zynNOxY8eoX78+w4cP55JLLqFTp07s37+fGjVqFO5z3333sWfPHmbMmHHaY/y7ZyYtLY06deqQmJiIl28Al//ffI5kFzDm2mb0bV7jtMdwRQ6Hgx07dlCvXj08PDysjlOpqS1ch9rCdagtXIu7tcfhzFz6vPU3eQVOPuh/EZfUDSM9PZ2YmBiOHj1KSEjIWZ9v+WmmkwUEBNC8eXMSEhK45pprAEhOTi5SzKSkpJzSW3MyHx8ffHx8TtlepUoVgoODuadHHONmbWXKykPc1rmJ25x6cjgcBAcHU6VKFbf4wazI1BauQ23hOtQWrsXd2uPDJZsp8PDlophQ+rSuh81mK8xdnL/Tlg8APllubi6bNm2iRo0a1K1bl8jISGbNmlX4eF5eHnPnzqVjx47n/RoDOkQT4O3B5uQMZm9OKY3YIiIicp7SsvKZssgcPvJgt9jz6mSwtJh54oknmDt3Ljt37mTJkiXccMMNpKenc+edd2Kz2Rg2bBijR4/m+++/Z/369QwcOBB/f39uu+22837NUH9v7rgkGoC3Zm/Dhc6yiYiIVDqfLkokM7eAxpFB9Ggcfl7HsPQ00969e7n11ls5dOgQ1atX55JLLmHx4sVER5vFxvDhw8nOzmbIkCEcOXKE9u3bM3PmTIKCgi7ode/tXI/JixJZs+co8xMO0aVh9dJ4OyIiIlICx3IL+OTvnQA8EF8fu/38hn5YWsxMmzbtrI/bbDZGjRrFqFGjSvV1qwf5cGu7Okz8O5G3ZyfQuUE1txk7IyIiUlF8sXQ3R7LyianqzxUtos77OC41ZqY8DepSH28PO8sSj7B4R6rVcURERCqV3AIHH87bAZi9Mh7n2SsDlbiYiQzx5eaLawPw9uwEi9OIiIhULt+s2EtKRi41Qny5tvXZ55E5l0pbzAAMjq+Pl4eNhdsPszxRvTMiIiLlocDh5P252wG4v0s9vD0vrByp1MVMzVA/rr/IrAbfmr3N4jQiIiKVw89r97MnNZuqAd7ccnGdCz5epS5mAIbEx+JhtzFv60FW7zlqdRwREZEKzek0ePcvs1fm7kvr4ud94ZP6Vfpipk5Vf65uZY6gnqCxMyIiImVq5sYDJKRkEuTrSf8O0aVyzEpfzMCJGQfhj00pbNifZnUcERGRCskwDN75yxzWcWeHGIJ9vUrluCpmgPrVAwuvb5+gsTMiIiJlYvbmFNbtS8Pf24O7OsWU2nFVzBz3UPdYAH5bn8zWAxkWpxEREalYDMNg/B/mcI4BHWKoGnjqotDnS8XMcQ0jgrg8LhJQ74yIiEhpO7lX5r7OdUv12CpmTjL0eO/Mz2v3s/1gpsVpREREKoay7JUBFTNFNIsKoWeTcAyDwgFKIiIicmH+2mL2yvh5lX6vDKiYOcVD3RsA8OPq/ew+nGVxGhEREfdWpFemY3Sp98qAiplTtKwdSpeG1XE4Dd6do94ZERGRC/HXlhTW7jV7Ze7vXK9MXkPFzGk8fHzszLcr97LvaLbFaURERNxTefTKgIqZ02obE0bH+lXJdxi8P2e71XFERETcUnn0yoCKmTM6MXbmy2V7SEpT74yIiEhJGIbB/xVewVR2vTKgYuaMLqkXRru6YeQ5nLqySUREpITmbDnImuO9Mvd1KbteGVAxc0Y2m43HejUEzN6ZvUd0ZZOIiEhxmGNltgJmr0y1MuyVARUzZ3VJvaqFY2c0K7CIiEjxlGevDKiYOafHe5u9M1+v2Muuw8csTiMiIuLaDMNg/J/mWJn+5dArAypmzqlNdBhdj88789af6p0RERE5mzlbD7Jmz1F8vezcXw69MqBiplgePT525vtVe9mhNZtERERO699rMJVHrwyomCmWVrVD6dkkHKcB/3e860xERESKOrlX5r4ynFfm31TMFNOwnmbvzE9r9pNwIMPiNCIiIq7l5Hll+l8STfWg8umVARUzxRZXM4Q+zSIxDAq70ERERMQ0d+tBVheOlalfrq+tYqYEhvVqgM0Gv65LYlNSutVxREREXIJhGLw5y5xXprx7ZUDFTIk0jgymX/MaAIWNJiIiUtnN3HiANXvT8Pf2KPdeGVAxU2LDejbAbjMbbt3eNKvjiIiIWMrhNBg30/yCf1enmHLvlQEVMyUWGx7E1a1qAvDmH+qdERGRyu3nNfvZciCDIF9P7u9c/r0yoGLmvDzcowEedhuzN6ewcvcRq+OIiIhYIt/hLPxiP7hrfUL8vSzJoWLmPNStFsB1rY/3zmjsjIiIVFJfL9/LrsNZVAv0ZmDHGMtyqJg5Tw/3aICn3cb8hEMsS0y1Oo6IiEi5ysl38NbxiWSHxMcS4ONpWRYVM+epdpg/N7atDVA48ElERKSy+GzxLpLTc6gR4stt7etYmkXFzAUY2j0Wbw87i3YcZuH2Q1bHERERKReZuQW8N2c7AI/0aICvl4eleVTMXICaoX7c0u6f3hnDMCxOJCIiUvYmLtjJ4WN5xFT15/o2tayOo2LmQj3YLRYfTzvLdx1hzpaDVscREREpU0ez8vhw/g4AHu3VEC8P60sJ6xO4uYhg38IR3K/O2IzTqd4ZERGpuD6Yt4OMnAIaRwZxZYsoq+MAKmZKxQPx9Qny9WRzcgY/rdlvdRwREZEykZKRw6S/EwF4vHcj7HabtYGOUzFTCkL9vRnc1Zz18I1ZW8grcFqcSEREpPS9+9d2svMdtKwdSs8m4VbHKaRippScWI9iT2o2XyzdbXUcERGRUrX3SBZTl5h/34Zf1gibzTV6ZUDFTKnx9/bkkR4NAHh7dgLHcgssTiQiIlJ63vozgTyHkw71qtIptprVcYpQMVOKbr64NjFV/TmUmcf/5u+0Oo6IiEip2H4wk29X7gPgicsaWZzmVCpmSpGXh53He5uN/NH8HRzOzLU4kYiIyIV7c9ZWHE6DHo3DaRNdxeo4p1AxU8r6Na9Bs6hgMnMLeOev7VbHERERuSAb96fzy9okgMIv7K5GxUwps9ttPNmnMWCuW7H3SJbFiURERM7fa79vBuCKFjVoGhVscZrTUzFTBjo3qEaHelXJczh5c1aC1XFERETOy8Lth/hry0E87TYe69XQ6jhnpGKmDNhsNp683Oyd+W7VXrYkZ1icSEREpGScToNXfjN7ZW5rX4d61QMtTnRmKmbKSKvaofRpFolhwGu/b7E6joiISIn8si6JtXvTCPD24OHjU4+4KhUzZeiJyxpht8Efmw6wYleq1XFERESKJbfAUThWZlDX+lQL9LE
|
||
|
|
"text/plain": [
|
||
|
|
"<Figure size 640x480 with 1 Axes>"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
"metadata": {},
|
||
|
|
"output_type": "display_data"
|
||
|
|
}
|
||
|
|
],
|
||
|
|
"source": [
|
||
|
|
"plt.plot(1e3*z_SI,1e4*func(z_SI,popt_z[0],popt_z[1]),label = \"harmonic fit\")\n",
|
||
|
|
"plt.plot(1e3*z_SI,1e4*B_tot_z_SI, label = r\"$ B_{tot} = \\sqrt{r^2 + z^2} $\")\n",
|
||
|
|
"plt.xlabel(\"x-axis [mm]\")\n",
|
||
|
|
"plt.xlim(-50,50)\n",
|
||
|
|
"plt.ylim(0,50)\n",
|
||
|
|
"plt.ylabel(\"B_tot [G]\")\n",
|
||
|
|
"plt.legend()\n",
|
||
|
|
"plt.show()"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": 289,
|
||
|
|
"id": "624fbda5-4658-4f9e-ba0f-33168c27e5c8",
|
||
|
|
"metadata": {},
|
||
|
|
"outputs": [
|
||
|
|
{
|
||
|
|
"name": "stdout",
|
||
|
|
"output_type": "stream",
|
||
|
|
"text": [
|
||
|
|
"T_expansion = T/4 = 66.2961644033217 ms\n",
|
||
|
|
"omega_x = 23.69362301623883\n",
|
||
|
|
"Magnification for t_tof = 20.0 ms: M = 0.47387246032477665\n",
|
||
|
|
"for z_start = 1 μm after T_expansion z_end = 1.4110642480981759 μm\n"
|
||
|
|
]
|
||
|
|
}
|
||
|
|
],
|
||
|
|
"source": [
|
||
|
|
"m = 2.69e-25\n",
|
||
|
|
"\n",
|
||
|
|
"k_x = -2*popt_x[0]*9.9*cs.mu_B\n",
|
||
|
|
"\n",
|
||
|
|
"omega = np.sqrt(k_x/m)\n",
|
||
|
|
"f = omega/(2*np.pi)\n",
|
||
|
|
"\n",
|
||
|
|
"T = 1/f\n",
|
||
|
|
"T_exp = T/4\n",
|
||
|
|
"print(f\"T_expansion = T/4 = {T_exp*1e3} ms\" )\n",
|
||
|
|
"print(f\"omega_x = {omega}\") \n",
|
||
|
|
"\n",
|
||
|
|
"t_tof = 20e-3\n",
|
||
|
|
"M = omega * t_tof\n",
|
||
|
|
"\n",
|
||
|
|
"print(f\"Magnification for t_tof = {t_tof*1e3} ms: M = {M}\")\n",
|
||
|
|
"\n",
|
||
|
|
"\n",
|
||
|
|
"start_z = 1e-6\n",
|
||
|
|
"d_t = 1e-3\n",
|
||
|
|
"def force(z):\n",
|
||
|
|
" return 2*0.248*z*9.9*cs.mu_B\n",
|
||
|
|
"z = start_z\n",
|
||
|
|
"v = 0\n",
|
||
|
|
"for t in np.arange(0,T_exp,d_t):\n",
|
||
|
|
" v = v + force(z)/m * d_t\n",
|
||
|
|
" #print(v)\n",
|
||
|
|
" z = z + v * d_t\n",
|
||
|
|
"print(f\"for z_start = 1 μm after T_expansion z_end = {z*1e6} μm\")\n"
|
||
|
|
]
|
||
|
|
},
|
||
|
|
{
|
||
|
|
"cell_type": "code",
|
||
|
|
"execution_count": null,
|
||
|
|
"id": "1561200c-5630-4c8a-b0b5-b6a535805eef",
|
||
|
|
"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.8.8"
|
||
|
|
}
|
||
|
|
},
|
||
|
|
"nbformat": 4,
|
||
|
|
"nbformat_minor": 5
|
||
|
|
}
|