2023-09-15 23:46:38 +03:00
|
|
|
|
{
|
|
|
|
|
|
"cells": [
|
|
|
|
|
|
{
|
|
|
|
|
|
"cell_type": "code",
|
2023-10-07 19:38:12 +03:00
|
|
|
|
"execution_count": 9,
|
2023-09-15 23:46:38 +03:00
|
|
|
|
"id": "3e9c91d2-a8a3-4cf4-a3bd-682f4979776a",
|
|
|
|
|
|
"metadata": {},
|
|
|
|
|
|
"outputs": [],
|
|
|
|
|
|
"source": [
|
|
|
|
|
|
"# %load /home/glebi/git/experiment-automation/processing_tools.py\n",
|
|
|
|
|
|
"import numpy as np\n",
|
|
|
|
|
|
"from scipy.optimize import curve_fit\n",
|
|
|
|
|
|
"import pandas as pd\n",
|
|
|
|
|
|
"\n",
|
|
|
|
|
|
"import matplotlib.pyplot as plt\n",
|
|
|
|
|
|
"import matplotlib\n",
|
|
|
|
|
|
"import scienceplots\n",
|
|
|
|
|
|
"\n",
|
|
|
|
|
|
"plt.style.use(['science', 'russian-font'])\n",
|
2023-09-29 12:46:12 +03:00
|
|
|
|
"plt.style.use([\"bright\"])\n",
|
2023-09-15 23:46:38 +03:00
|
|
|
|
"\n",
|
|
|
|
|
|
"matplotlib.rcParams.update({\n",
|
|
|
|
|
|
" 'figure.figsize': [6, 4],\n",
|
|
|
|
|
|
" 'savefig.facecolor': 'white',\n",
|
2023-09-29 13:25:15 +03:00
|
|
|
|
" 'figure.dpi': 200.0,\n",
|
2023-09-15 23:46:38 +03:00
|
|
|
|
" 'font.size': 12.0,\n",
|
2023-09-29 13:25:15 +03:00
|
|
|
|
"})\n",
|
|
|
|
|
|
"\n",
|
|
|
|
|
|
"def line_func(x, arg1, arg2):\n",
|
|
|
|
|
|
" return x * arg1 + arg2"
|
2023-09-15 23:46:38 +03:00
|
|
|
|
]
|
|
|
|
|
|
},
|
|
|
|
|
|
{
|
|
|
|
|
|
"cell_type": "code",
|
2023-10-07 19:38:12 +03:00
|
|
|
|
"execution_count": 10,
|
2023-09-15 23:46:38 +03:00
|
|
|
|
"id": "5e765906-9269-4ba2-81ba-fb2fd3d798f1",
|
|
|
|
|
|
"metadata": {},
|
|
|
|
|
|
"outputs": [
|
|
|
|
|
|
{
|
|
|
|
|
|
"data": {
|
2023-09-29 13:41:43 +03:00
|
|
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAABBgAAALUCAYAAAChG4GcAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAB7CAAAewgFu0HU+AACiPElEQVR4nOzdeXxU5aH/8e8kBBJUMhBEx7qQCUsHVybBLra2moyi2LiQgNXYRU2mdrm3v7Ykprc1SXtvMeBd21uZoN7aRm9hBlGRimSCtrV2STLgUlKQDKitoyhhwIXEwMzvD26miWSZJCc5mcnn/Xr5ek3nnOecb6Ao853neY4lGo1GBQAAAAAAMAIpZgcAAAAAAACJj4IBAAAAAACMGAUDAAAAAAAYMQoGAAAAAAAwYhQMAAAAAABgxCgYAAAAAADAiFEwAAAAAACAEaNgAAAAAAAAI0bBAAAAAAAARoyCAQAAAAAAjNgkswOMpuLiYi1atEjl5eVDGhcOh1VXV6d169YpHA7LarVKkiorK1VUVDTgWJ/PJ4/H0+s9l8s15AwAAAAAACSSpCoYwuGwgsGg1q1bp7q6OoXDYdnt9iFdo66uThUVFSorK1NjY2OsXAiHwyotLVV7e7vKysr6HOtyudTe3i6v1xu7b/e4nJwcNTQ0DDkPAAAAAACJIGmWSFgsFmVnZ6uiokI5OTnD+iC/atUqud1ueb1e1dbWxsoFSfL7/fL7/aqoqOhzbHFxsZqbm9XY2Njr3larVV6vV1arVS6Xa8iZAAAAAABIBElTMESjUR08eFANDQ0qKyvTjBkzhjTe5/OpoqJC5eXlKigoOOF4Q0ODwuGwwuHwCcf8fr98Pp/Kysp6lRI9VVZWKhgM9ltQAAAAAACQyJKmYBiJ7mUM0vEioC/dyya8Xu8Jx7r3XBhohkL33g11dXUjjQsAAAAAwLiTVHswDNfKlSsVDodVVFTU7wwEu91+wuaN3Xw+X+ycgVitVoXDYfn9/j5nSQAAAAAAkKiYwaC/zyoYzh4Jfr8/9nqwgqH7eCAQGPJ9AAAAAAAYzyb8DIZgMBjbVyEvLy/2ns/nU1tbm6xWq5YvXy6n09nn+O6yoL+ZDz117wvR1NQ07Lxvv/22nnrqKc2ePVsZGRnDvg4AAAAAAPE4cuSI9u3bpyuvvFIzZ87s97wJXzD0nIFgtVpVV1entrY2ud1u2e12BQIBlZaWym63a+3atScUCQcOHBjyPfvaKDJeTz31lEpKSoY9HgAAAACA4aivr9fNN9/c7/EJXzC0tbXFXgcCAbW1tam2tjb2ntPpVGNjo7Kzs5Wdna29e/f2KhmGUhZ0j2tvbx923tmzZ0uSfvjDHyo7O7vPc2bOnKlTTz110Gu1traqpKRE9fX1cjgcw840GsZrtvGaSyLbcI3XbOM1l0S24RivuSSyDdd4zTZec0lkG47xmksi23CN12zjNZdEtuEaTra33npLb7/9dp/H9u7dq+9///uxz6P9mfAFQ8+CwOPxqKGh4YRzrFaramtr5Xa7VVFR0Wuzx+GUBSOZwdC9LOL73/9+v+dUVVWpuro67ms6HI5+l4CYbbxmG6+5JLIN13jNNl5zSWQbjvGaSyLbcI3XbOM1l0S24RivuSSyDdd4zTZec0lkG66hZKuurlZNTc2A5wy2TH/CFww9DfQL3/3Uh7q6Ornd7ti53fsqjLWBmiibzTbGaQAAAAAAicztdquwsLDPY90zIgYz4QuGnssdBnqKRM8nRHg8ntgshng2dxwNRrRkNptNVVVV47KQGM/Zxqvx/GtGtqEbr7nGu/H66zZec0lkSzbj+deMbEM3XnONd+P51228ZhuvuaTxnW08G86vm81mG/mvczRJFRQURCVFy8vLBzyvtrY2KikqKdrW1jbguVarNSop6nQ6TxhvtVrjztRz/FC1tLREJUVbWlqGfQ0MH7/+AH8OgGiUPwcAfwaAifXnIN6fNWVk9UTi6zkzYTDdyyF67qHQPYMhnn0VuvdrGMo9AQAAAABIBBO+YOi5zGCwkqC7IOi5LKJnWTDY+O7jFAwAAAAAgGQz4QsGu90eKwyCwWDcY7rl5eXFXg/2RInu6w+01wPGN9aAAfw5ACT+HAD8GQD4c9CXCV8wSNKyZcskDVwwhMPh2AyE5cuXx963Wq2xWRCDje8+v/uJFEg8NptN1dXV/EsEExp/DgD+HAD8GQD4c9AXCgZJFRUVkqR169b1e47f75d0fElFUVFRr2OVlZWSJK/X2+/49evXS5LKyspGlBUAAAAAgPEo6QuGeDZftNvtKi8vVyAQkM/n6/Oc7hKirxKhqKhIRUVFqqur63cWQ21traxWq2pra+MPDwAAAABAgkjKgiEQCKi5uVnS8ZkH8eytUFtbq7KyMhUXF6uuri72fjAYjO2Z0NLS0u8GjWvXrlVBQYFcLlev+4XD4V7jAQAAAABIRpPMDmAUl8ul5ubmE2YsBINB5eTkSDq+/8GyZcvk8Xj6vIbH41FxcbFqa2tjMxZmzJihoqIiNTQ0DHh/q9WqhoYG+Xw+ud3uXo+vdLlcKi8vH9kPCAAAAADAOJY0BcNgBUC8CgoKRrQJY/dyCQAAAAAAJpKkXCIBAAAAAADGFgUDAAAAAAAYMQoGAAAAAAAwYhQMAAAAAABgxJJmk8eJprW1NfbaZrPJZrOZmAYAAAAAkExCoZBCoZCk3p8/B0LBkKBKSkpir6uqqlRdXW1eGAAAAABAUvF4PKqpqRnSGAqGBFVfXy+HwyFJzF4AAAAAABjK7XarsLBQ0vEZDD2/5O4PBUOCcjgccjqdZscAAAAAAAxTJBpRR1en0tOmKMUyvrZIHM5SfAoGAAAAAADGUOv+3bq/6SFt2dWoI10dykhL1+L5+bpt0c1yzJpndrxhG18VCQAAAAAASezxnVt03YO3aONLm3Wkq0OSdKSrQxtf2qzrHrxFj+/cYnLC4aNgAAAAAABgDLTu360Vm6t0NHKsz+NHI8e0YnOVWvfvHuNkxqBgAAAAAABgDNzf9FC/5UK3o5FjeqDp4TFKZCwKBgAAAAAARlkkGtGWXY1xnfvkLr8i0cgoJzIeBQMAAAAAAKOso6sztufCYI50daijq3OUExmPggEAAAAAgFGWnjZFGWnpcZ2bkZau9LQpo5zIeBQMAAAAAACMshRLihbPz4/r3KvmFyjFkngf1xMvMQAAAAAACei2RTdrUkrqgOdMSknVrYtuGqNExqJgAAAAAABgDDhmzVOVq7zf45NSUrV6SY0cs+aNYSrjUDAAAAAAADBGTjv5VGVMSteSj14R25MhIy1dN5x3jR794i9UuGCxyQmHb5LZAQAAAAAAmCjy51yqP32jQVMnZygSjaijq1PpaVMScs+FD0v8nwAAAAAAgATw9nvt6jzaqamTMyQd3/hx6uSMpCgXJAoGAAAAAADGxMqn/0Off7jM7BijhoIBAAAAAIBRdqjjsJ7c1agr5l1mdpRRwx4MCaq1tTX22mazyWazmZgGAAAAADCQx/78pI5FjmrpedeYHSUuoVBIoVBIUu/PnwOhYEhQJSUlsddVVVWqrq42LwwAAAAAoF/RaFTrnn9Ul8+5VKeePNPsOHHxeDyqqakZ0hgKhgRVX18vh8MhScxeAAAAAIBx7HDnO0qfNEXLLrjW7Chxc7vdKiwslHR8BkPPL7n7Q8GQoBwOh5xOp9kxAAAAAACDyEyfpg1f+Jmi0ajZUeI2nKX4bPIIAAAAAMAoebfzPQX+9oKi0agsFovZcUYVBQMAAAAAAKPkidatWlZ/m9589y2zo4w6CgYAAAAAAEbJ+hce1aX2T+j0U2aZHWXUUTAAAAAAADAK/rL/ZT0f+rNuvPAGs6OMCQoGAAAAAABGwS+f36hTT8rSZTmfMjvKmKBgAAAAAABgFJwx7XSVXnyL0lInxgMcJ8ZPCQAAAADAGCv72BfMjjCmmMEAAAAAAIDB1j//qN54Z7/ZMcYUBQMAAAAAAAZqO7BPlVv+WU2vbTc7ypiiYAAAAAAAwEDrnn9U0zMydcW8y8yOMqYoGAAAAAAAMEjn0Q+08aUndP25SzR
|
2023-09-15 23:46:38 +03:00
|
|
|
|
"text/plain": [
|
2023-09-29 13:25:15 +03:00
|
|
|
|
"<Figure size 1200x800 with 1 Axes>"
|
2023-09-15 23:46:38 +03:00
|
|
|
|
]
|
|
|
|
|
|
},
|
|
|
|
|
|
"metadata": {},
|
|
|
|
|
|
"output_type": "display_data"
|
|
|
|
|
|
}
|
|
|
|
|
|
],
|
|
|
|
|
|
"source": [
|
|
|
|
|
|
"df = pd.read_csv(\"VAH.csv\")\n",
|
|
|
|
|
|
"cols = df.columns\n",
|
|
|
|
|
|
"\n",
|
2023-09-29 13:41:43 +03:00
|
|
|
|
"plt.plot(df[cols[0]] * 1e-3, df[cols[1]]*1e-1, marker=\".\", lw=.5, color=\"C2\", ls=\"--\")\n",
|
2023-09-15 23:46:38 +03:00
|
|
|
|
"plt.xlabel(\"Напряжение, кВ\")\n",
|
|
|
|
|
|
"plt.ylabel(\"Сила тока, мкА\")\n",
|
|
|
|
|
|
"\n",
|
2023-09-29 12:46:12 +03:00
|
|
|
|
"plt.savefig(\"hand-VAH.png\")\n",
|
2023-09-15 23:46:38 +03:00
|
|
|
|
"plt.show()"
|
|
|
|
|
|
]
|
|
|
|
|
|
},
|
|
|
|
|
|
{
|
|
|
|
|
|
"cell_type": "markdown",
|
|
|
|
|
|
"id": "d084ce30-37aa-440c-ab41-591ab094d21f",
|
|
|
|
|
|
"metadata": {},
|
|
|
|
|
|
"source": [
|
|
|
|
|
|
"сопротивление нашей цепи"
|
|
|
|
|
|
]
|
|
|
|
|
|
},
|
|
|
|
|
|
{
|
|
|
|
|
|
"cell_type": "code",
|
2023-10-07 19:38:12 +03:00
|
|
|
|
"execution_count": 11,
|
2023-09-15 23:46:38 +03:00
|
|
|
|
"id": "5b0cd085-4398-4b10-b6cc-3147d4e616f8",
|
|
|
|
|
|
"metadata": {},
|
|
|
|
|
|
"outputs": [
|
|
|
|
|
|
{
|
|
|
|
|
|
"data": {
|
|
|
|
|
|
"text/plain": [
|
|
|
|
|
|
"13991140.642303433"
|
|
|
|
|
|
]
|
|
|
|
|
|
},
|
2023-10-07 19:38:12 +03:00
|
|
|
|
"execution_count": 11,
|
2023-09-15 23:46:38 +03:00
|
|
|
|
"metadata": {},
|
|
|
|
|
|
"output_type": "execute_result"
|
|
|
|
|
|
}
|
|
|
|
|
|
],
|
|
|
|
|
|
"source": [
|
|
|
|
|
|
"1263.4 / (903e-7)"
|
|
|
|
|
|
]
|
|
|
|
|
|
},
|
|
|
|
|
|
{
|
|
|
|
|
|
"cell_type": "markdown",
|
|
|
|
|
|
"id": "cd1ec375-3bdd-4bcb-88a0-fa6dcf173ec4",
|
|
|
|
|
|
"metadata": {},
|
|
|
|
|
|
"source": [
|
|
|
|
|
|
"максимальное напряжение на резисторе, который снимает ток"
|
|
|
|
|
|
]
|
|
|
|
|
|
},
|
|
|
|
|
|
{
|
|
|
|
|
|
"cell_type": "code",
|
2023-10-07 19:38:12 +03:00
|
|
|
|
"execution_count": 12,
|
2023-09-15 23:46:38 +03:00
|
|
|
|
"id": "cfe69bd0-e5c9-45af-826b-22e759bc6eca",
|
|
|
|
|
|
"metadata": {},
|
|
|
|
|
|
"outputs": [
|
|
|
|
|
|
{
|
|
|
|
|
|
"data": {
|
|
|
|
|
|
"text/plain": [
|
|
|
|
|
|
"8.32"
|
|
|
|
|
|
]
|
|
|
|
|
|
},
|
2023-10-07 19:38:12 +03:00
|
|
|
|
"execution_count": 12,
|
2023-09-15 23:46:38 +03:00
|
|
|
|
"metadata": {},
|
|
|
|
|
|
"output_type": "execute_result"
|
|
|
|
|
|
}
|
|
|
|
|
|
],
|
|
|
|
|
|
"source": [
|
|
|
|
|
|
"1600e-7 * 52e3"
|
|
|
|
|
|
]
|
|
|
|
|
|
},
|
|
|
|
|
|
{
|
|
|
|
|
|
"cell_type": "markdown",
|
|
|
|
|
|
"id": "6a53d766-e438-4378-b7b0-db3c46578f2c",
|
|
|
|
|
|
"metadata": {},
|
|
|
|
|
|
"source": [
|
|
|
|
|
|
"аналитическая функция вах"
|
|
|
|
|
|
]
|
|
|
|
|
|
},
|
|
|
|
|
|
{
|
|
|
|
|
|
"cell_type": "code",
|
2023-10-07 19:38:12 +03:00
|
|
|
|
"execution_count": 13,
|
2023-09-15 23:46:38 +03:00
|
|
|
|
"id": "8a4273c4-db67-40b9-b887-926d8e6866e1",
|
|
|
|
|
|
"metadata": {},
|
|
|
|
|
|
"outputs": [],
|
|
|
|
|
|
"source": [
|
|
|
|
|
|
"# def vah_func(u):\n",
|
|
|
|
|
|
"# return S_e * (1.537e10 / )"
|
|
|
|
|
|
]
|
|
|
|
|
|
},
|
2023-09-29 13:25:15 +03:00
|
|
|
|
{
|
|
|
|
|
|
"cell_type": "markdown",
|
|
|
|
|
|
"id": "59f5605a-6a22-49de-98a2-bc83e9d8d49b",
|
|
|
|
|
|
"metadata": {},
|
|
|
|
|
|
"source": [
|
|
|
|
|
|
"обработка данных"
|
|
|
|
|
|
]
|
|
|
|
|
|
},
|
2023-09-15 23:46:38 +03:00
|
|
|
|
{
|
|
|
|
|
|
"cell_type": "code",
|
2023-10-07 19:38:12 +03:00
|
|
|
|
"execution_count": 14,
|
2023-09-29 13:25:15 +03:00
|
|
|
|
"id": "87f15973-cb08-477e-a81d-2a81e1536105",
|
|
|
|
|
|
"metadata": {},
|
|
|
|
|
|
"outputs": [],
|
|
|
|
|
|
"source": [
|
|
|
|
|
|
"U_min = 900\n",
|
|
|
|
|
|
"R = 26.99e3 # Ом\n",
|
|
|
|
|
|
"\n",
|
|
|
|
|
|
"column_names = [\"Voltage\", \"Current\"]\n",
|
|
|
|
|
|
"\n",
|
|
|
|
|
|
"experiment_names = [\"Первая серия измерений\", \"Вторая серия измерений\"]\n",
|
|
|
|
|
|
"file_idxs = [[1, 2, 3], [4, 5]]\n",
|
|
|
|
|
|
"\n",
|
|
|
|
|
|
"dots_result = [None, None]\n",
|
|
|
|
|
|
"ex_result = [None, None]\n",
|
|
|
|
|
|
"\n",
|
|
|
|
|
|
"for num, (ex_name, f_idxs) in enumerate(zip(experiment_names, file_idxs)):\n",
|
|
|
|
|
|
" result = pd.DataFrame(columns=column_names)\n",
|
|
|
|
|
|
"\n",
|
|
|
|
|
|
" for i in f_idxs:\n",
|
|
|
|
|
|
" df = pd.read_csv(f\"RISING-{i}.csv\", header = None)\n",
|
|
|
|
|
|
" df = df.loc[lambda x : x[4] > U_min]\n",
|
|
|
|
|
|
"\n",
|
|
|
|
|
|
" U = df[4] # В\n",
|
|
|
|
|
|
" U_cur = df[10]\n",
|
|
|
|
|
|
" I = U_cur / R # А\n",
|
|
|
|
|
|
"\n",
|
|
|
|
|
|
" result = pd.concat([result, pd.concat(dict(zip(column_names, [U, I])), axis=1)], ignore_index=True)\n",
|
|
|
|
|
|
" \n",
|
|
|
|
|
|
" result[\"Current\"] -= min(result[\"Current\"])\n",
|
|
|
|
|
|
" ex_result[num] = result.groupby(\"Voltage\").mean().reset_index()\n",
|
|
|
|
|
|
" dots_result[num] = result"
|
|
|
|
|
|
]
|
|
|
|
|
|
},
|
|
|
|
|
|
{
|
|
|
|
|
|
"cell_type": "code",
|
2023-10-07 19:38:12 +03:00
|
|
|
|
"execution_count": 15,
|
2023-09-15 23:46:38 +03:00
|
|
|
|
"id": "24817123-a6cd-47a5-a153-7d4e6ab68135",
|
|
|
|
|
|
"metadata": {},
|
|
|
|
|
|
"outputs": [
|
|
|
|
|
|
{
|
|
|
|
|
|
"data": {
|
2023-09-29 12:46:12 +03:00
|
|
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAABBgAAALLCAYAAABTmzHSAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAB7CAAAewgFu0HU+AAEAAElEQVR4nOz9eXwb550n+H8AAiRAUSJA+ghlJ5FAOQkdJ7ZAOs7pHARtJ51WDoNmulu9OzuxiHZ6emd+rx2TVm/32t7uaRlK+jfHznYEyD1X1JNYgO1OOp3EImgnzm0RsJw4URKLkBw7oqwDKFIkAeKq/QOoIgrEjSqAx+f9eulV4IOq+hYKoM3ni+f7PDpRFEUQERERERERETVA3+oLICIiIiIiIqKNjwkGIiIiIiIiImoYEwxERERERERE1DAmGIiIiIiIiIioYUwwEBEREREREVHDmGAgIiIiIiIiooYxwUBEREREREREDWOCgYiIiIiIiIgaZmj1BVBtLl++jGeeeQa7du2C2Wxu9eUQERERERHRJheLxXDu3DncfffduOaaa0ruxwTDBvO1r30Nf/Znf9bqyyAiIiIiIqIt5tixY/ijP/qjks8zwbDBdHZ2AgD+6q/+Cp/4xCc0jXX69Gns378fx44dw8DAAOMwDuMwjuZxmhmLcRiHcdZ/nGbGYhzGYRzGaUWsjRJHOn7Xrl1l92OCYYPp6OgAAOzevRt2u70pMQcGBpoSi3EYh3EYpxWxGIdxGGf9x2lmLMZhHMZhnFbE2ihxKpXpc5JHIiIiIiIiImoYEwwbjDShRrmJNdTS19eHhx9+GH19fZsiTrNstvu22eI0y2a7b818f/hZYJxmxmmWzXbf+N+E+m2292izxWmWzXbfNlucZtps904niqKoaQRSVSgUwuDgIILBYNOGBW0WvHcE8HNAq/hZIICfA1rFzwIB/BzQKn4WlKq9HxzBQEREREREREQNY4KBiIiIiIiIiBrGBANtGZuxZotqx88BSfhZIICfA1rFzwIB/BzQKn4W6sM5GDYY1gIRERERERFRM3EOBiIiIiIiIiJqGiYYiIiIiIiIiKhhTDAQERERERERUcOYYCAiIiIiIiKihjHBQEREREREREQNY4KBiIiIiIiIiBrGBAMRERERERERNYwJBiIiIiIiIiJqGBMMRERERERERNQwJhiIiIiIiIiIqGFMMBARERERERFRw5hgICIiIiIiIqKGMcFARERERERERA0ztPoCqD6nT5+WH/f19aGvr6+FV0NERERERESbydzcHObm5gAo+5/lMMGwQe3fv19+/PDDD+ORRx5p3cUQERERERHRpuLxePDoo4/WdAwTDBvUsWPHMDAwAAAcvUBERERERFuCKIrQ6XStvowtweVyYd++fQCyIxjyv+QuhQmGDWpgYAB2u73Vl0FERERERKQ5URQBUVx9rNMx0aCxekrxmWAgIiIiIiKi9U0UV5MMOh10AMAEw7rDVSSIiIiIiIho3RJzIxdQsJXbad1ggoGIiIiIiIjWLbkUomDLEon1hwkGonXM5XJBl6svU/Nff39/q18aEREREVH1pL9l9fpsYoHJhXWJCQaidczj8UAURczOziom9ZTaK/0LBoOK8wWDQfl8REREREQbhZRcQH6SgdYdJhiINgCbzYaxsTH556GhoaqOy09KOBwOrjxCRERERBsaEwvrGxMMRERERHlCoRD6+/sxMjKy5rmRkRH09/cjFAq14MqIiIjWNyYYiIiIiPJ4PB6Ew2EEAgEEAgG5PRQKIRAIIBwOw+PxtPAKiYiI1icmGIiIiIjyDA4OAsiWp+WXpNlstjX7EBER0SpDqy+AiIiIaD0ZHx9HMBiE1+uF1WqFxWIBAAiCID8/Pj7eugskIiJap5hgICIiIirg8XgwOTkJv98vr7zT398Pp9OpGMlAREREq5hgoLIyGRELy4lWX8a6taOzHXo9Z7IlItrIBEHA5ORkyeel+RZcLhfcbnfRfdxutzzSgYiIaKtigoHKWlhOYPSvv9Hqy1i3fH+xD5aujlZfhmqkP7JnZmZgs9nk4cAulwtOp1Ox7+joKMLhMARBQCQSwfj4ONxuN/x+v/zHeCQSAQCMjY1hYmKiaEzpPPkzsjudTvT09CASiUAQBITDYTgcDrhcrrJLbQYCAcUf/5FIBA6HAwcPHqzqD/9QKCRP7gZAcUwgEJDvB5Bd9nN0dFSVYdKHDx/GE088AQDy6x4bG8P4+HjJ6272eyUIAg4cOLDmvbLb7bDZbPD5fIpYxfY5evTomtfTrHuu1vUDwPDwMHw+HwKBAHw+n2L4fCQSkZeVLXwfXC4XZmZm5HNbLBacPXu24mczEAhgdHRUvhfl7kMrfofria/m+6HW59Lj8SAUCmF4eBiCIMjH5/83R/qsjoyMIBwOw2KxYHp6mksAExER5TDBQEQAAL/fjwMHDsi1x5JAIICRkRGMj48rZk0fGxtDOByWv/WTOhb9/f2YmpqS95P+YPd4PGv+WAcgdwBGRkYQCATgdDrlNokgCBgcHITX64Xb7S7a0XG5XACgiC0IAkZHR2G1WjE1NQWHw1H0tUv7BQIBOBwORacxfx+r1SrHKHWuWoTDYbnjODU1pRh2HQgEMDg4WPSeteK9slgs8vsyODiIUCgEh8OhOB5Y+37a7XbFNUqafc/Vun5BEBAKhTA4OIijR4+uOR4AJicnMTo6CpvNpnhf878FlxINhw4dKvmNeOE1AUAwGCzZmW3V73A98dX8PKn5ubTb7XA4HPD7/XA4HEVfq81mg91ulxOfTC4QERGt4ioSRAS/34/R0VE4HI41nR2HwwGPxwOv1wuv1yu3O51OTExMyN+ier1e+Vv3fHa7HdPT0wiHwxgcHJS/qS4k/ZHe09Oz5rn8zsjk5OSac0hLyUlJhvzjpA6e9I1jMZOTk/JSdFNTU0W/Uc5vy59Vvl7hcBj9/f0Ih8NrkgtSDLvdjgMHDija18N7Jb3+cnXo0rGl7lUr7nnhuRq5/unp6TUjAiRut1vugI6Ojq55vr+/X/6s5r9Ppa5jcHBQ/r0odc2t/lzUE19SzfshPVfuc6DG+woU/29QvkAgAL/fX9W+REREWw0TDERbnDRUGUDJb1KlDkexGmWpE1jumzy73S6POijW4apGfqchf5gzkO2gSp2f/CH1EqkjWGrdemkYeDPrp6UOptvtLtohOnToEPx+v+K1bpT3anR0tOK9bMU9r1al65+YmKj4rbU02qJUkqanpwfj4+MQBEHurBbj8XgqloS0+nPRaPxqqPE5qeZzWYkgCHC73VxBgoiIqASWSFBZOzrb4fuLfa2+jHVrR2d7qy+hYZOTkxAEQa5HLsXhcCAQCCAUCpUcNlyOy+XC4cOHEQqF4Pf7S377W8qhQ4fkOIXHjo2Nwev1YmhoqGgHor+/HwDkb8wLjY2Nwe/3QxAECIKgeadXGnEBAPfdd1/Ja5LKByQb4b06fPgwRkZGEAwG1ySC8jX7nler0vVbLJaKJQ0A5KTBwYMHS+4zOTkJr9eLQ4cOFb3HxZJlpc7Tys+FWvG1VO3nspLR0VF4PJ6qPgNERFuNKIrQ6bSf/Dx6+TKs11yjeRwAWF5eRmdnZ1NibRYcwUBl6fU6WLo6+K/Ev82wgoTU0a00BF3qAM7MzNQVx2azyecoVrdejCAI8gR30vwL0nJx+ex2O6LR6JrzShNEFjsmnzRUHIBiMj2tSCMp8u9JIalGPL8js57fKyA7smR2draqb3ebfc+rUcv1lzvH4OCgXJ5TbmJEm80Gh8OBUChUtNPr9XrXlP0U0+rPRTPiN/L5UON9BbKJFJfLxSUqiYgKiKIIMZMBcltRFDWJE718GdFf/QrIbaOXL2sSB8gmFpYvXwZy2+XlZc1ibTYcwUC0xUlDuGdmZsoOX7bZbJiYmGhooj2bzYZQKFS2g1E4l0I4HEYkEqlqWHI4HIbf78fU1BQikQiGhoZKlk0UcrvdcLlc8Hg88hBwrToSxVZMqOW49fJe5ZMmCKwlIdHMe15JPdcPZN8Tt9uNSCSCUCgkr8ZQ7col0lwUhw4dWjO
|
2023-09-15 23:46:38 +03:00
|
|
|
|
"text/plain": [
|
2023-09-29 13:25:15 +03:00
|
|
|
|
"<Figure size 1200x800 with 1 Axes>"
|
2023-09-15 23:46:38 +03:00
|
|
|
|
]
|
|
|
|
|
|
},
|
|
|
|
|
|
"metadata": {},
|
|
|
|
|
|
"output_type": "display_data"
|
|
|
|
|
|
}
|
|
|
|
|
|
],
|
|
|
|
|
|
"source": [
|
2023-09-29 13:25:15 +03:00
|
|
|
|
"for num, (ex_name, dots, line) in enumerate(zip(experiment_names, dots_result, ex_result)):\n",
|
2023-09-29 12:46:12 +03:00
|
|
|
|
" plt.plot(dots[\"Voltage\"]*1e-3, dots[\"Current\"]*1e6, alpha=0.01, ls=\"\", marker=\".\", markersize=2, color=f\"C{num}\")\n",
|
|
|
|
|
|
" plt.plot(line[\"Voltage\"]*1e-3, line[\"Current\"]*1e6, lw=1.5, alpha=1, color=f\"C{num}\", label=ex_name)\n",
|
|
|
|
|
|
"else:\n",
|
|
|
|
|
|
" df = pd.read_csv(\"VAH.csv\")\n",
|
|
|
|
|
|
" cols = df.columns\n",
|
|
|
|
|
|
" plt.plot(df[cols[0]]*1e-3, df[cols[1]]*1e-1, marker=\".\", lw=1, markersize=3, ls=\"--\", color=\"C2\", label=\"Ручные измерения\")\n",
|
|
|
|
|
|
" \n",
|
|
|
|
|
|
"plt.xlabel(fr\"Напряжение, кВ\")\n",
|
|
|
|
|
|
"plt.ylabel(fr\"Сила тока, мкА\")\n",
|
2023-09-29 13:25:15 +03:00
|
|
|
|
"# plt.yscale(\"log\")\n",
|
|
|
|
|
|
"\n",
|
|
|
|
|
|
"plt.legend()\n",
|
2023-09-29 12:46:12 +03:00
|
|
|
|
"plt.savefig(\"all-VAH.png\")\n",
|
2023-09-29 13:25:15 +03:00
|
|
|
|
"plt.show()"
|
|
|
|
|
|
]
|
|
|
|
|
|
},
|
2023-09-29 12:46:12 +03:00
|
|
|
|
{
|
|
|
|
|
|
"cell_type": "markdown",
|
|
|
|
|
|
"id": "8c61ec35-5368-4a1c-8708-9f5b08d969f7",
|
|
|
|
|
|
"metadata": {},
|
|
|
|
|
|
"source": [
|
|
|
|
|
|
"### Прямые с осциллографа"
|
|
|
|
|
|
]
|
|
|
|
|
|
},
|
2023-09-29 13:25:15 +03:00
|
|
|
|
{
|
|
|
|
|
|
"cell_type": "code",
|
2023-10-07 19:38:12 +03:00
|
|
|
|
"execution_count": 16,
|
2023-09-29 13:25:15 +03:00
|
|
|
|
"id": "f964e34e-f019-49d3-9815-86764aad6a64",
|
|
|
|
|
|
"metadata": {},
|
|
|
|
|
|
"outputs": [
|
|
|
|
|
|
{
|
|
|
|
|
|
"name": "stdout",
|
|
|
|
|
|
"output_type": "stream",
|
|
|
|
|
|
"text": [
|
2023-09-29 12:46:12 +03:00
|
|
|
|
"0: k= -10.50\tb= -15.26\n",
|
|
|
|
|
|
"1: k= -7.50\tb= -17.57\n"
|
2023-09-29 13:25:15 +03:00
|
|
|
|
]
|
|
|
|
|
|
},
|
|
|
|
|
|
{
|
|
|
|
|
|
"data": {
|
2023-09-29 13:50:41 +03:00
|
|
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAABE0AAAL1CAYAAAAy1hoDAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAB7CAAAewgFu0HU+AADrrUlEQVR4nOz9a1BbZ74u+j5j6MLdiLu52CDJmGDj2AjfYrsvicXM7Nmkq9ZsEc9Tdq19PsxAen48u9u01zrVSbpqTzfu3vtrItJ7Ve1V9qndhtXrnBlm90yDk07HdnxBAts4xMaS8AVjGwPC2AaENMb5oEjhIgkJhCTg+VVRIYNXY/w1NDAaj96LIMuyDCIiIiIiIiIimkOMdwFERERERERERImIoQkRERERERERUQAMTYiIiIiIiIiIAmBoQkREREREREQUAEMTIiIiIiIiIqIAGJoQEREREREREQXA0ISIiIiIiIiIKABlvAug9eXp06f49NNPUVZWhpSUlHiXQ0RERERERGvc5OQkBgYG8OabbyI3NzeixzI0oZj69NNPcezYsXiXQUREREREROvM6dOncfTo0Ygew9CEYqqsrAyA92KtrKyc87O+vj4cO3Ys4M8SQSLXx9qWLpHrS+TagMSuj7UtXSLXl8i1AYldH2tbukSuL5FrAxK7Pta2dIlcXyLXBiR2fYlcG7D8+nyP992PRoKhCcWUb0hOZWUlDAZDwDahfpYIErk+1rZ0iVxfItcGJHZ9rG3pErm+RK4NSOz6WNvSJXJ9iVwbkNj1sbalS+T6Erk2ILHrS+TagOXXt5QpIjgRLBERERERERFRAAxNKGEUFhbivffeQ2FhYbxLCSiR60vk2hJdIp+7RK4NSOz6Erm2RJfI5y6RawMSu75Eri3RJfK5S+TagMSuL5FrS3SJfO4SuTYgsetL5NriTZBlWY53EbR+WK1W1NTUwGKxJHS3L4oevuZEc/F3gmgu/k4QzcXfCaKFlvt7sZzHs6cJEREREREREVEADE2IiIiIiIiIiAJgaEJEK4rjI4nm4u8E0Vz8nSCai78TRAvF8/eCc5pQTHGMJhEREREREcUS5zQhIiIiIiIiIooyhiZERERERERERAEwNCEiIiIiIiIiCoChCRERERERERFRAAxNiIiIiIiIiIgCYGhCRERERERERBQAQxMiIiIiIiIiogAYmhARERERERERBcDQhIiIiIiIiIgoAGW8C6D1qa+vz/99YWEhCgsL41gNERERERERrSVDQ0MYGhoCMPf+M1IMTSgujh075v/+vffew/vvvx+/YoiIiIiIiGhNMZvN+OCDD5a9H4YmFBenT59GZWUlALCXCREREREREUVVY2MjfvKTnwDw9jSZ/cF9JBiaUFxUVlbCYDDEuwwiIiIiIiJag6I1DQRDE6IEJksS3P0OyOMTEDIzoCzXQhA5fzMREREREVEsMDQhSlCu7l5Mnm2HNOr0bxNzNEipr4O6uip+hREREREREa0T/MiaKAG5unvxwnxmTmACANKIEy/MZ+Dq7o1PYUREREREROsIQxOiFSJLEmZu2eC60oOZWzbIkhT24ybPtgOyHKSBjMmz7WHvj4iIiIiIiJaGw3OIVsByhta4+x0LepjMJ4064e53QFWhj0K1REREREREFAh7mhBF2XKH1sjjE2EdJ9x2REREREREtDQMTYiiKBpDa4TMjLCOFW47IiIiIiIiWhqGJkRRFMnQmmCU5VqI2ZqQ+xBzNFCWayOqbalzrBAREREREa1XnNOEKIqiMbRGEEWkvF2HF+YzgXusCAJS6usgiOFnnly+mIiI1qLGxka0tLREfb86nQ42my3q+yUiotWHPU2I5llOj4xoDa1RV1chrfHogh4nYo4GaY1HIwo6uHwxERGtVWazGbIsw2azwWAwLNi+2JfFYpmzP4vF4t8fERERwJ4mRHMst0eGb2hNqCE64Q6tUVdXQbVzG9z9DsjjExAyM6As10bUwyTcOVZUO7dFtF8iIqJEotPpcOTIEVitVgDA7t27w3rc7KDFaDTO+X8iIiKAPU2I/EL2yPjodFg9MnxDayAIQRpENrRGEEWoKvRQ790FVYU+4mAjGnOsEBERERER+VitVuj1etTW1i74WW1tLfR6vT/EXgsYmhAhjB4ZAF7+f/5/kNzuRfcVzaE1y8Xli4mIiIiIKJrMZjPsdjs6OzvR2dnp3261WtHZ2Qm73Q6z2RzHCqOLw3OIEF6PDPnZBCbb/oTk7+2FWJgfstdHNIbWRAOXLyYiIiIiomiqqakB4B0aOXs4pE6nW9BmLWBoQoTwe1oIsgSX9QaE9HSoKnTe8CTIUBzf0Jp4iuYcK0RElLg8koxexzBGJqaQk5GMKm0eFGKQoaJERETL0NDQAIvFgpaWFmRlZUGj0QAAnE6n/+cNDQ3xKzDKGJoQIfyeFirDDijyczFzyw6X5TqEjHSotoYOT+JpJZYvJiKixHK+dxAftvfgifOlf1uBJhXv1u3CoariOFZGRERrldlsRlNTE9ra2vwrjun1ephMpjk9TtYChiZEiKxHhiCKSNpfDWnUiZnb3vBE3JABZYUOYkFewoUn6uoqoPHoslYFIiKixHS+dxC/PnNxQS7+2PkSvz5zEb86eoDBCRERLYvT6URTU1PQn/vmL2lsbERzc3PANs3Nzf4eKasNQxMiLK1HhpitQdJ+gzc8uWWD6+o1iJkboNyqg1iQm1DhSaLMsUJERNHjkWR82N4TalV5fNjeg9e2FXGoTpT4bhy6urqg0+n8XdEbGxthMpnmtK2vr4fdbofT6cTo6CgaGhrQ3NyMtrY2/w3G6OgoAODIkSM4fvx4wGP69jN7JQqTyYTs7GyMjo7C6XTCbrfDaDSisbEx5LLJnZ2dc25oRkdHYTQaceLEibBuZqxWq38CSABzHtPZ2ek/H4B3Cef6+vqodNE/deoU/vCHPwCA/3kfOXIEDQ0NQeuO9WvldDrxzjvvLHitDAYDdDodWltb5xwrUJuPP/54wfOJ1TmPVv0AcPjwYbS2tqKzsxOtra1zhm6Mjo76lwif/zo0Njaiq6vLv2+NRgOHw7HotdnZ2Yn6+nr/uQh1HuLxO7yU40fz9YjWdWk2m2G1WnH48GE4nU7/42f/m+O7Vmtra2G326HRaHDu3LnVv5y7TBRDFotFBiBbLJZ4lxLQtPWG7PzlSXm0ocn/5TxxUp623lj0se6no/LUhavyy3/7izz1xSXZ/WhYliQpBlUTEdF61HPnsWxsOrvoV8+dx/EuNSaam5tlABG/z/A9xmg0hmzX2toqazQa+fjx43O2d3R0yADkhoaGBe1n19TQ0CAfP35cNpvNc9pZLBZZo9HIOp0uZN1Go1EGIJtMpgU/Gxsbk3U6nQxAbm5uDvj4hoaGBTWOjY3599vR0RH02LPbGY1GeWxsLGAb33MNta9I2Gw22WAwyDqdTrbZbHN+1tHREfScxfu1MhgMi15TvvNpMBgC/jxe51yWl1e/ryaNRiO3trYGfOzx48dlAAFfV1n2Xqu+Gua/hoE0NDTIGo1m0d/9eF8XkR7fJxrXUzT3YzKZQtY7u02gf6/iZTn3oQxNKKYSPTSRZVmWPB7Z9c0defpyt+z65o4seTwRPd49PCJPnb/iDU++vCy7HzM8ISKi6DvXfTes0ORc9914lxoTKxmatLa2hrwBMJvNMoAFN1Oy7L2hW6wu3/sjAAFvImX5uxvNYDcqofYRKmCQZdkfuAQ79uznEIqvTaAb/EjZbDb/zXegusbGxmSTybTg5i4RXivffkLdVPrOebA28Tjn84+9lPp9ocliv4O+G/hAN+fNzc3+10mj0YTcj81mk81ms7+eYOch3tdFNI4f6vWIVpvFrstw9uMLgRbbT6wt5z6UffOJ5vGteqPeuwuqCn3EQ1gUudlQH9gN9X4DAAGuy91wXeiCZ3gEcrA+1ERERBHKyUiOajsKzNdNHkDQsfq+YQCBxvz7hhYYjcagXdQNBoO/a399ff2S6pw98eLsLvYA0NHRAbvdjpqamjnDOXx8wwJ8Qw7m8w1BiOV8BI2NjQC85zzQpJI
|
2023-09-29 13:25:15 +03:00
|
|
|
|
"text/plain": [
|
|
|
|
|
|
"<Figure size 1200x800 with 1 Axes>"
|
|
|
|
|
|
]
|
|
|
|
|
|
},
|
|
|
|
|
|
"metadata": {},
|
|
|
|
|
|
"output_type": "display_data"
|
|
|
|
|
|
}
|
|
|
|
|
|
],
|
|
|
|
|
|
"source": [
|
|
|
|
|
|
"opt_arr = []\n",
|
|
|
|
|
|
"for num, (ex_name, line) in enumerate(zip(experiment_names, ex_result)):\n",
|
2023-09-29 12:46:12 +03:00
|
|
|
|
" x = 1e3 / line[\"Voltage\"]\n",
|
2023-09-29 13:25:15 +03:00
|
|
|
|
" y = np.log(line[\"Current\"]) - 2 * np.log(line[\"Voltage\"])\n",
|
2023-09-15 23:46:38 +03:00
|
|
|
|
" \n",
|
2023-09-29 13:25:15 +03:00
|
|
|
|
" popt, pcov = curve_fit(line_func, x, y)\n",
|
|
|
|
|
|
" print(f\"{num}: k= {popt[0]:.2f}\\tb= {popt[1]:.2f}\")\n",
|
|
|
|
|
|
" opt_arr.append(popt)\n",
|
2023-09-15 23:46:38 +03:00
|
|
|
|
" \n",
|
2023-09-29 13:25:15 +03:00
|
|
|
|
" plt.plot(x, y, ls=\"\", marker=\".\", markersize=5, color=f\"C{num}\", label=ex_name)\n",
|
|
|
|
|
|
" plt.plot(x, line_func(x, *popt), lw=0.5, alpha=.5, color=f\"C{num}\")\n",
|
|
|
|
|
|
"\n",
|
2023-09-29 13:41:43 +03:00
|
|
|
|
"plt.xlabel(r\"$\\dfrac {10^3} U$, у.е.\")\n",
|
|
|
|
|
|
"plt.ylabel(r\"$\\log \\dfrac I {U^2}$, у.е.\")\n",
|
2023-09-15 23:46:38 +03:00
|
|
|
|
"\n",
|
|
|
|
|
|
"plt.legend()\n",
|
2023-09-29 12:46:12 +03:00
|
|
|
|
"plt.savefig(\"osc-NORDGAME-FAULER.png\")\n",
|
|
|
|
|
|
"plt.show()"
|
|
|
|
|
|
]
|
|
|
|
|
|
},
|
|
|
|
|
|
{
|
|
|
|
|
|
"cell_type": "markdown",
|
|
|
|
|
|
"id": "bb204b4f-dc40-4da5-bb0f-465162d6460f",
|
|
|
|
|
|
"metadata": {},
|
|
|
|
|
|
"source": [
|
|
|
|
|
|
"### Прямая ручная"
|
|
|
|
|
|
]
|
|
|
|
|
|
},
|
|
|
|
|
|
{
|
|
|
|
|
|
"cell_type": "code",
|
2023-10-07 19:38:12 +03:00
|
|
|
|
"execution_count": 17,
|
2023-09-29 12:46:12 +03:00
|
|
|
|
"id": "6b0c2a09-1033-43c5-baa4-02616c6bdb91",
|
|
|
|
|
|
"metadata": {},
|
|
|
|
|
|
"outputs": [
|
|
|
|
|
|
{
|
|
|
|
|
|
"name": "stdout",
|
|
|
|
|
|
"output_type": "stream",
|
|
|
|
|
|
"text": [
|
|
|
|
|
|
"1: k= -11.41\tb= -14.52\n"
|
|
|
|
|
|
]
|
|
|
|
|
|
},
|
|
|
|
|
|
{
|
|
|
|
|
|
"data": {
|
2023-09-29 13:50:41 +03:00
|
|
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAABEkAAAL2CAYAAAC9qcjXAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAB7CAAAewgFu0HU+AACqWElEQVR4nOz9e1DcZb73/X5+QM4EGkIOHciBhhxQc6CbaMZDNKERxxnjxEDiKFO1q+4JrHW79x+Pa4TJOPsx1mMtJaNP1fNUzVo2WrX3vkVnQqNxMicjTaIxHoEmiQdMAg05EnOADsZDTKD3Hy56hYRj0vCjm/erKjVt9/W7+NDdo+lvf6/rMgKBQEAAAAAAAABjXJTZAQAAAAAAAEYDiiQAAAAAAACiSAIAAAAAACCJIgkAAAAAAIAkiiQAAAAAAACSKJIAAAAAAABIokgCAAAAAAAgiSIJAAAAAACAJCnG7AAYW86ePaudO3dq/vz5mjRpktlxAAAAAAAR7rvvvlNLS4tyc3OVlJTU71iKJBhRO3fuVEFBgdkxAAAAAABjTHl5uR599NF+x1AkiVAlJSWqrKyUz+eTzWaT3W7XihUrVFxcPOC1fr9fzz77rLxerySpra1NkuR0OrV582ZZLJbrzjV//nxJP745MzIy+h3b0NCggoKCQY0dCuYd/rmZl3mZd2TnZl7mZd6Rm3c452Ze5mXekZt3OOdm3tE3b/c13Z9H+0ORJML4fD7l5+erqKhIdXV1slgs8nq9evbZZ1VSUiKXyyW32y273d7r9R6PRyUlJdq8ebNKS0uD93u9XmVnZ2vr1q1yuVwqLCy8rnzdS2wyMjL6zHC1oYwdCuYd/rmZl3mZd2TnZl7mZd6Rm3c452Ze5mXekZt3OOdm3tE372C2fGDj1gji9/uVk5Mjt9utwsLCYMeH3W6X2+1WcXGxfD6fHA6HfD7fNdf7fD6VlJSorq5OeXl5PR7rnkOSioqK5PF4hv33AQAAAABgJFEkiSAlJSUqLS2VzWbr9fHS0tJg4SQ/P/+ax10ul7xer3Jycnq93ul09pgrXFmtVj311FOyWq1jet7hnns4hNtzHG7zDpdwex74/9x/C7fnONzmHS7h9jyE27zDKdyei3Cbd7iE2/MQbvMOp3B7LsJt3uEy7HkDiBgWi2XAMaWlpQFJAUmBqqqqHo85nc7gY01NTb1eb7PZApIG9bN6U1dXF5AUqKurC+lYjC68duGJ1y188dqFJ1638MVrF5543cITr1v4Gk2v3VCy0EkSIbxer/x+vwzD6HcpzJVrtqqqqno8VlRUJIvFIqfT2Wc3Svcmrn09DgAAAABAuKJIEiGu3GPE5XL1OS4xMTF42+/393gsLy9P7e3t1xRPrhzffc2VS28AAAAAAIgEFEkihNPplN1ul8ViUVFRUZ/jujtBJCktLW1IP6OiokKSZLFYtHnz5usLOgThtjYO/43XLjzxuoUvXrvwxOsWvnjtwhOvW3jidQtf4fraGYFAIGB2CIycrVu3qqSkRNKPy20G2xHi9/uVmpoqv98/pOuu5vV65XA4VFdXN2xHeAEAAAAA0G0on0NjRigTRonupTh2u31QhQ6/36/a2loVFRXJZrPppZdeCklxo6Ghoc/HrFZr2FUbAQAAAADmaW1tVWtra6+P9ff582oUScYQj8cT3LvkpZde6ndsZWWlXC6X2tra5PV6ZbfbVVpaGrLuj4KCgj4fe+qpp7Rly5aQ/BwAAAAAQORzuVx6+umnb3geltuMIWlpafL5fCotLVVxcfGQrvV6vcrPz1dbW5vcbvcNL7cpLy9XRkZGr2PoJAEAAAAADMVAnSQFBQUst8F/Kyoquu4CifTj8py6ujolJCQoJyfnuufplpGRwZ4kAAAAAICQCNWX7ZxuMwL8fr8MwwjZH4fDMaSfX1ZWprKyshsubFgsFhUWFkqSSkpK5PV6r3suAAAAAABGGzpJRoDFYlFVVZX8fn9I5rPZbIMe6/F4VFRUJJfLFSxw3IicnByVlZVJ+nHNV/dGsAAAAAAAhDuKJCPkevfwuBHd+4i43W7l5eUNaqzFYpHb7e6zEGOxWIK3a2trQxkXAAAAAABTsdwmQvl8PmVnZ/dZIPH7/fJ4PMF/drlc8vl88nq9dIcAAAAAAMYkiiQRyO/3KycnR9XV1X12sNTW1vbYU+TKDpEVK1b0OXf3EcKSOd0xAAAAAAAMF4okEcbv9wc7SPo7PaaqqqrH4xs3bpTFYlFTU1O/S3PcbnfwdlFRUWhCAwAAAAAwCrAnSYTJzs6W0+mUz+fr0fVxJZ/Pp7Kysh5FDrvdrsLCQpWWlva53Mbj8QSX6PS3bwkAAAAAAOGIIkkEycnJkdfrHfTRvFcXOUpLS1VSUiKHw6HNmzfL6XQGl+FUVlZq06ZNslgsKi0tHXAjWAAAAAAAwg1FkghRWVnZYyPWgfS1FKe0tFQ+n0+lpaV69tln5ff71dbWpsTERBUWFmrz5s099i8BAAAAACBSUCSJEHl5eQoEAiGZy2azjekTbjq7OlVzvF6nL5zVjNgkrUjJVHRUtNmxAAAAAADDjCIJcIWdh3brmeoXdLLjVPC+5Dirnsx+XLkLV5uYDAAAAAAw3DjdBvgvOw/t1mPbi3sUSCTpREerHtterJ2HdpuUDAAAAAAwEiiSAPpxic0z1S8ooN6XLAUU0DPVL6izq3OEkwEAAAAARgpFEkBSzfH6azpIrnay45RqjtePUCIAAAAAwEijSAJIOn3hbEjHAQAAAADCD0USQNKM2KSQjgMAAAAAhB+KJICkFSmZmh03q98xyXFWrUjJHKFEAAAAAICRRpEEkBQdFa3fZ/+bDBm9Pm7I0JPZjys6Kvqaxzq7OvXR0Vrt+OItfXS0ls1dAQBARCgqKpJhGIP+k5aWppycHJWVlZkdHQCuG0US4L/kLlytP67bek1HSXKcVX9ct1W5C1dfc83OQ7t1j+tBPfqnf9H/9tff69E//YtWu37BccEAACDsuVwuBQIBNTU1yW63X3P/1X/cbrfsdrtKSkqUkJCgrVu3mpgeAK6PEQgEej/zFBgGXq9XDodDdXV1Pf5jO5p0dnWq5ni9Tl84qxmxSVqRktlrB8nOQ7v12PbiXo8NNmT0WVgBAAAIN1u3blVJSYkkqampSTabrc+xfr9fqamp8vv9KiwslMvlGqmYANCroXwOpZMEuEp0VLRWzs3S2pvu08q5WX0usXmm+oVeCySSFFBAz1S/wNIbAAAQcRITE/t93GKx6KWXXpIklZWVBYsrABAOKJLAFA0NDfJ6vfJ6vWptbTU7zpDVHK/XyY5T/Y452XFKNcfrRygRAADA6JGXlxe8zbIbAGZobW0NfuZsaGgY9HUUSWCKgoICORwOORyOsGzBPH3hbEjHAQAARBqLxRK87fV6zQsCYExyuVzBz5wFBQWDvi5mGDMBfSovL1dGRoYkyWq1mpxm6GbEJoV0HAAAGL0Gu18ZekpMTJTf7zc7BoAxqqioSGvXrpX040qGwRZKKJLAFBkZGaN249bBWJGSqdlxs/pdcpMcZ9WKlMwRTAUAAEJt56Hdeqb6hR7/zU+Os+rJ7MfZoH0APp8veLu/jV4BYDhYrdbr+kKe5TZAL37ovKRPTzWoK9DV6+PRUdH6ffa/yZDR6+OGDD2Z/TjfMgEAEMa6T7K7+kuREx2temx7sXYe2m1SstHP4/EEbxcWFvZYegMAoxlFEqAXR9uPa1fjHlUe2KFz37b3OiZ34Wr9cd1WzY6b1eP+5Dgrx/8CABDmOMnuxnSfaON0OnvsP+fxeJSTkyPDMHr8ycnJ6dF50j3HlWPS0tJUWVkpv9+v/Px8ORyOHo87HA7l5+cHr+9vjN/vH/F5erN169Zrno+EhATl5+ersrLymvF+v19FRUXBuXNycpSTk9Pr2KEI5XPhcDjk8/lUVlamnJwc5efnB7N2j+0tb/fvdeXzMJjlWh6PRwkJCT3eS2VlZX3+noN9/rp/17S0NCUkJATf05WVlcHruve7GOzmxIP9+aPtPX41r9eroqKiYP7u1zg/P7/HazHQ6zFqBYARVFd
|
2023-09-29 12:46:12 +03:00
|
|
|
|
"text/plain": [
|
|
|
|
|
|
"<Figure size 1200x800 with 1 Axes>"
|
|
|
|
|
|
]
|
|
|
|
|
|
},
|
|
|
|
|
|
"metadata": {},
|
|
|
|
|
|
"output_type": "display_data"
|
|
|
|
|
|
}
|
|
|
|
|
|
],
|
|
|
|
|
|
"source": [
|
|
|
|
|
|
"df = pd.read_csv(\"VAH.csv\")\n",
|
|
|
|
|
|
"cols = df.columns\n",
|
|
|
|
|
|
"\n",
|
|
|
|
|
|
"volt_h = df[cols[0]]\n",
|
|
|
|
|
|
"curr_h = df[cols[1]]*1e-7\n",
|
|
|
|
|
|
"\n",
|
|
|
|
|
|
"x = 1e3 / volt_h\n",
|
|
|
|
|
|
"y = np.log(curr_h) - 2 * np.log(volt_h)\n",
|
|
|
|
|
|
"\n",
|
|
|
|
|
|
"popt, pcov = curve_fit(line_func, x, y)\n",
|
|
|
|
|
|
"print(f\"{num}: k= {popt[0]:.2f}\\tb= {popt[1]:.2f}\")\n",
|
|
|
|
|
|
"opt_arr.append(popt)\n",
|
|
|
|
|
|
"\n",
|
|
|
|
|
|
"plt.plot(x, y, ls=\"\", marker=\".\", markersize=5, color=f\"C2\", label=\"Ручные измерения\")\n",
|
|
|
|
|
|
"plt.plot(x, line_func(x, *popt), lw=0.5, alpha=.5, color=f\"C2\")\n",
|
|
|
|
|
|
"\n",
|
2023-09-29 13:41:43 +03:00
|
|
|
|
"plt.xlabel(r\"$\\dfrac {10^3} U$, у.е.\")\n",
|
|
|
|
|
|
"plt.ylabel(r\"$\\log \\dfrac I {U^2}$, у.е.\")\n",
|
2023-09-29 12:46:12 +03:00
|
|
|
|
"\n",
|
|
|
|
|
|
"plt.legend()\n",
|
|
|
|
|
|
"plt.savefig(\"hand-NORDGAME-FAULER.png\")\n",
|
|
|
|
|
|
"plt.show()"
|
|
|
|
|
|
]
|
|
|
|
|
|
},
|
|
|
|
|
|
{
|
|
|
|
|
|
"cell_type": "markdown",
|
|
|
|
|
|
"id": "a212083d-2d53-4522-860d-8eba58766461",
|
|
|
|
|
|
"metadata": {},
|
|
|
|
|
|
"source": [
|
|
|
|
|
|
"### Все вместе"
|
|
|
|
|
|
]
|
|
|
|
|
|
},
|
|
|
|
|
|
{
|
|
|
|
|
|
"cell_type": "code",
|
2023-10-07 19:38:12 +03:00
|
|
|
|
"execution_count": 18,
|
2023-09-29 12:46:12 +03:00
|
|
|
|
"id": "d668aa7a-0a40-475b-9876-698f4bc62c0a",
|
|
|
|
|
|
"metadata": {},
|
|
|
|
|
|
"outputs": [
|
|
|
|
|
|
{
|
|
|
|
|
|
"data": {
|
2023-09-29 13:50:41 +03:00
|
|
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAABEkAAAL1CAYAAAA7Pbp5AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAB7CAAAewgFu0HU+AAEAAElEQVR4nOzdeXzT9f0H8Nc36ZGepKUClUNI8SgKalI8QAEltfXCwVLxqJtOaXVz89hsBTeBbYqpm7rNTVt0P7fh2YoOpxwJcih4kEQBFTwakKscpU3vpG3y/f1Rv1+aNknTNm16vJ6PRx7i9/M93t+jbT7v7+cQRFEUQUREREREREQ0zCnCHQARERERERER0UDAJAkREREREREREZgkISIiIiIiIiICwCQJEREREREREREAJkmIiIiIiIiIiAAwSUJEREREREREBIBJEiIiIiIiIiIiAEBEuAOg4aWyshLr16/HxIkTERMTE+5wiIiIiIiIaIhramrC/v37kZWVhZSUlIDrMklC/Wr9+vXIzc0NdxhEREREREQ0zKxatQq33nprwHWYJKF+NXHiRABtD2d6enrAdffs2YPc3Nyg1u1rjGXgxjGQYhkocTCWgR3HQIploMQxkGIZKHEwloEdx0CKZaDEwVgGdhwDKZaBEgdjGdhxhDoWaV9SfTQQJkmoX0ldbNLT06HVaoPapjvr9jXGMnDjAAZOLAMlDoCxDOQ4gIETy0CJAxg4sQyUOADGMpDjAAZOLAMlDoCxDOQ4gIETy0CJA2AsAzkOILSxBDPkAwduJSIiIiIiIiICkyQ0gKWmpmLp0qVITU0NdyiMZYAbKNdkoMTBWAa+gXJNBkocAymWgRLHQDNQrstAiWMgxTJQ4hhoBsp1GShxDKRYBkocA81AuS4DJY5wxiKIoij26xFpWLPZbNDpdLBarQOm+Rb1DO8lDTZ8Zmmw4TNLgw2fWRps+MwOH92512xJQkREREREREQEJkmIiIiIiIiIiAAwSUJEPTSQ+isSBYPPLA02fGZpsOEzS4MNn1nyhWOSUL9ivz8iIiIiIiLqTxyThIiIiIiIiIiom5gkISIiIiIiIiICkyRERERERERERACYJCEiIiIiIiIiAsAkCRERERERERERACZJiIiIiIiIiIgAMElCRERERERERASASRIiIiIiIiIiIgBMkhARERERERERAWCShIiIiIiIiIgIABAR7gBoeNqzZ4/879TUVKSmpoYxGiIiIiIiIhpKKioqUFFRAcC7/tkVJkkoLHJzc+V/L126FMuWLQtfMERERERERDSkFBcXY/ny5d3ejkkSCotVq1YhPT0dANiKhIiIiIiIiEIqPz8f8+bNA9DWkqT9i/pAmCShsEhPT4dWqw13GERERERERDQE9XRYByZJiAYQ97FKNG+3wF1ZBWVKMqJmZEA5OiXcYREREREREQ0LTJIQDRCubRY0rloNeDwAgBYAzg1bEZu7ANEzM8IbHBERERER0TDAJAlRiPWkNYj7WKVXgkTm8aBx1WpETJ7IFiVERERERER9jEkSohDqaWuQ5u2WzgkSiceD5u0WxMzP7oOIiYiIiIiISKIIdwBEQ0VXrUHcxyr9b1tZFXjfXZQTERERERFR7zFJQhQiwbQG8UeZkhxw312VExERERERUe+xuw1RiPSmNUjUjAw4N2z1nWRRKBA1I7iBWzk7DhERERERUc8xSULUgftYJZzvboTbUYvISeODTjQoU5LR0kW537LRKYjNXdC5u45CgdjcBUEdn7PjEBHRcJOfn4+SkpKQ71ej0aC8vDzk+yUiooGP3W2I2nFts6B26Z/R/MlncH9dDue6zahd+mc4P9zR5bZRMzIAhZ8fqSBag0TPzEDisgehyp6DyIxpUGXPQeKyB4NKcPRmPBQiIqLBqri4GKIoory8HFqtttPyrj5Wq9Vrf1arVd4fERENT0ySEP1ATjSIoneBKKLpP2/Caf4Qor8xR3CqNUinREk3WoMoR6cgZn424hfdgpj52UF3lenNeChERESDnUajwcKFC+X/z8gIrgVl+8SKXq/3+n8iIhqe2N2G6AcBEw0Amkr/B9fmj6C6ahaiLtVBiOz84xM9MwMRkyf2+7ggnB2HiIiIiIj6g81mQ05ODjQaDUwmk1dZZmYm7HY7SktLB23imUkSoh8Ek0jwnDiJxpffQtP/zFBlXo7oyy+GoIr2WkdqDdKfejMeChERERERUbCKi4tht9tht9thNpuh1+sBtCVPzGazvE5xcXE4w+wxdrch+kFXiQTFSLX8b7GmDk1l76Fm8RNoWmOCp76hj6MLrLfjoRAREREREQVDp9MBaOvq2L57o0aj6bTOYMQkCdEPuko0xN93FxIW/wKRF5wrLxYbm+B8dyNqFj+Bxjfegae6pp+i9RaK8VCIiIiIiIi6kpeXh7y8PNjtdiQlJXl92pcPVuxuQ/SDYKfhjb/nNrgrjsO5fjOaP/m8bd3mFrg2boNr88eIuuRCqLJmQzn6tH6NP1zjoRARUfgdqqzDest+HK1qwJjkOGRlTMS4lIRwh0VERENUcXExCgsLUVZWJs8IlpaWBoPB4NWiZDBikoSonWATDcrUUYi7/Uaors+Ey/QBXB/uAFpaALcbzdssaN5uReSF50F19RxETBjbb/GHYzwUIiIKr3WWfXh6tRUez6nZ2d7Y+jUeWKBDdsakMEZGRERDgcPhQGFhod9yaeyR/Px8GI1Gn+sYjUao1eq+CC/kmCQh6qA7iQblyCTE3jQPqmuuhOv9bXBt/ghikxMQRbTYdqPFthsRU85qS5acOQmCIPRx9ERENJwcqqzrlCABAI9HxNOrrThvYgpblISQVFGwWCzQaDRwOBwA2ioGBoPBa92cnBzY7XY4HA5UVVUhLy8PRqMRZWVlcoWiqqpt0PiFCxeioKDA5zGl/dhsNnmZwWBAcnIyqqqq4HA4YLfbodfrkZ+fH3A2CbPZ7FWBqaqqgl6vx+LFi4OqvNhsNnnARgBe25jNZvl6AG1TKufk5ISkyX1RURFef/11AJDPe+HChcjLy/Mbd3/fK4fDgUWLFnW6V1qtFhqNBqWlpV7H8rXOypUrO51Pf13zUMUPAHPnzkVpaSnMZjNKS0vlmKXrK03Z3fE+5Ofnw2KxyPtWq9XYt29fl8+m2WxGTk6OfC0CXYdw/Az35PihvB+hei6Li4ths9kwd+5cOBwOefv2v3OkZ1Wa4UatVmPjxo2Db5YbkagfWa1WEYBotVrDHUqf8DQ2iU1rN4nVv/mDWJVX6PWpeeLvouvzL0WP2x3uMImIaIh4Ye0uUV/4ht/PC2t3hTvEfmM0GkUA3f6eIW2j1+sDrldaWiqq1WqxoKDAa7nJZBIBiHl5eZ3Wbx9TXl6eWFBQIBYXF3utZ7VaRbVaLWo0moBx6/V6EYBoMBg6lVVXV4sajUYEIBqNRp/b5+XldYqxurpa3q/JZPJ77Pbr6fV6sbq62uc60rkG2ld3lJeXi1qtVtRoNGJ5eblXmclk8nvNwn2vtFptl8+UdD21Wq3P8nBdc1HsXfxSTGq1WiwtLfW5bUFBgQjA530VxbZnVYqh4z30JS8vT1Sr1V3+7If7ueju8SWheJ5CuR+DwRAw3vbr+Pp9FS7dqYcySUL9Sno4V61aJVqtVtFqtYpHjhwJd1gh53E1i87NH4mOxU90SpY4lj0lOj+2iZ7W1nCHSUREg9wfX/4oYJLkjy9/FO4Q+01fJklKS0sDfuEvLi4WAXSqPIliWwWuq7ik70cAfFYaRfFUxdJfxSTQPgIlFERRlBMs/o7d/hwCkdbxVaHvrvLycrmy7Suu6upq0WAwdKrMDYR7Je0nUCVSuub+1gnHNe947J7ELyVJuvoZlCrsvirjRqNRvk9qtTrgfsrLy8Xi4mI5Hn/XIdzPRSiOH+h+hGqdrp7LYPYjJX262k9/OHLkiFznXLVqVdB/Hzi7DYVFbm4udDoddDrdoJ0/OxAhKhLRsy9B4h9+g9ifLYTi9NFymefIMTT+83XUPvpnuLZ8DLGlJYyRduY+Vommt9ahfuUraHprHdzHKsMdEhER+TE
|
2023-09-29 12:46:12 +03:00
|
|
|
|
"text/plain": [
|
|
|
|
|
|
"<Figure size 1200x800 with 1 Axes>"
|
|
|
|
|
|
]
|
|
|
|
|
|
},
|
|
|
|
|
|
"metadata": {},
|
|
|
|
|
|
"output_type": "display_data"
|
|
|
|
|
|
}
|
|
|
|
|
|
],
|
|
|
|
|
|
"source": [
|
|
|
|
|
|
"for num, (ex_name, line) in enumerate(zip(experiment_names, ex_result)):\n",
|
|
|
|
|
|
" x = 1e3 / line[\"Voltage\"]\n",
|
|
|
|
|
|
" y = np.log(line[\"Current\"]) - 2 * np.log(line[\"Voltage\"])\n",
|
|
|
|
|
|
" \n",
|
2023-09-29 13:41:43 +03:00
|
|
|
|
" plt.plot(x, y, ls=\"\", marker=\".\", markersize=4, color=f\"C{num}\", label=ex_name)\n",
|
|
|
|
|
|
" plt.plot(x, line_func(x, *(opt_arr[num])), lw=1, color=f\"C{num}\")\n",
|
2023-09-29 12:46:12 +03:00
|
|
|
|
"else:\n",
|
|
|
|
|
|
" x = 1e3 / volt_h\n",
|
|
|
|
|
|
" y = np.log(curr_h) - 2 * np.log(volt_h)\n",
|
|
|
|
|
|
"\n",
|
2023-09-29 13:41:43 +03:00
|
|
|
|
" plt.plot(x, y, ls=\"\", marker=\".\", markersize=4, color=f\"C2\", label=\"Ручные измерения\")\n",
|
|
|
|
|
|
" plt.plot(x, line_func(x, *(opt_arr[-1])), lw=1, color=f\"C2\")\n",
|
2023-09-29 12:46:12 +03:00
|
|
|
|
" \n",
|
2023-09-29 13:41:43 +03:00
|
|
|
|
"plt.xlabel(r\"$\\dfrac {10^3} U$, у.е.\")\n",
|
|
|
|
|
|
"plt.ylabel(r\"$\\log \\dfrac I {U^2}$, у.е.\")\n",
|
2023-09-29 12:46:12 +03:00
|
|
|
|
"\n",
|
|
|
|
|
|
"plt.legend()\n",
|
|
|
|
|
|
"plt.savefig(\"osc-NORDGAME-FAULER.png\")\n",
|
2023-09-15 23:46:38 +03:00
|
|
|
|
"plt.show()"
|
|
|
|
|
|
]
|
|
|
|
|
|
},
|
2023-09-29 13:25:15 +03:00
|
|
|
|
{
|
|
|
|
|
|
"cell_type": "markdown",
|
|
|
|
|
|
"id": "31fee377-ff25-4382-8f41-228790745a73",
|
|
|
|
|
|
"metadata": {},
|
|
|
|
|
|
"source": [
|
|
|
|
|
|
"Угол наклона полученной прямой получается по формуле\n",
|
|
|
|
|
|
"\n",
|
|
|
|
|
|
"$\\tan \\alpha = -0.683 \\cdot \\dfrac {\\phi^\\frac 3 2} \\beta$\n",
|
|
|
|
|
|
"\n",
|
|
|
|
|
|
"зная $\\phi = 4.5$ эВ -- работу выхода, находим $\\beta$ -- форм-фактор острия:"
|
|
|
|
|
|
]
|
|
|
|
|
|
},
|
2023-09-15 23:46:38 +03:00
|
|
|
|
{
|
|
|
|
|
|
"cell_type": "code",
|
2023-10-07 19:38:12 +03:00
|
|
|
|
"execution_count": 19,
|
2023-09-29 13:25:15 +03:00
|
|
|
|
"id": "21ac3719-9a57-4f80-9858-fc3c018fb1ad",
|
|
|
|
|
|
"metadata": {},
|
|
|
|
|
|
"outputs": [
|
|
|
|
|
|
{
|
|
|
|
|
|
"name": "stdout",
|
|
|
|
|
|
"output_type": "stream",
|
|
|
|
|
|
"text": [
|
2023-09-29 13:41:43 +03:00
|
|
|
|
"beta:\t6.21e-01\n",
|
|
|
|
|
|
"beta:\t8.69e-01\n",
|
|
|
|
|
|
"beta:\t5.71e-01\n"
|
2023-09-29 13:25:15 +03:00
|
|
|
|
]
|
|
|
|
|
|
}
|
|
|
|
|
|
],
|
|
|
|
|
|
"source": [
|
|
|
|
|
|
"for opt in opt_arr:\n",
|
|
|
|
|
|
" tan_a = opt[0]\n",
|
|
|
|
|
|
" beta = -0.683 * (4.5**(3/2)) / tan_a\n",
|
2023-09-29 12:46:12 +03:00
|
|
|
|
" print(f\"beta:\\t{beta:.2e}\")"
|
2023-09-29 13:25:15 +03:00
|
|
|
|
]
|
|
|
|
|
|
},
|
2023-09-29 13:41:43 +03:00
|
|
|
|
{
|
|
|
|
|
|
"cell_type": "markdown",
|
|
|
|
|
|
"id": "d030c573-6440-4bcc-bab6-2887bd67624d",
|
|
|
|
|
|
"metadata": {},
|
|
|
|
|
|
"source": [
|
|
|
|
|
|
"# тоже самое, но для каждого файла в отдельности"
|
|
|
|
|
|
]
|
|
|
|
|
|
},
|
|
|
|
|
|
{
|
|
|
|
|
|
"cell_type": "code",
|
2023-10-07 19:38:12 +03:00
|
|
|
|
"execution_count": 20,
|
2023-09-29 13:41:43 +03:00
|
|
|
|
"id": "9e61c01d-77b7-4984-b906-68f55031aed4",
|
|
|
|
|
|
"metadata": {},
|
|
|
|
|
|
"outputs": [
|
|
|
|
|
|
{
|
|
|
|
|
|
"data": {
|
|
|
|
|
|
"text/plain": [
|
|
|
|
|
|
"['RISING-1.csv',\n",
|
|
|
|
|
|
" 'RISING-2.csv',\n",
|
|
|
|
|
|
" 'RISING-3.csv',\n",
|
|
|
|
|
|
" 'RISING-4.csv',\n",
|
|
|
|
|
|
" 'RISING-5.csv',\n",
|
|
|
|
|
|
" 'BOTH-1.csv',\n",
|
|
|
|
|
|
" 'BOTH-2.csv',\n",
|
|
|
|
|
|
" 'BOTH-3.csv']"
|
|
|
|
|
|
]
|
|
|
|
|
|
},
|
2023-10-07 19:38:12 +03:00
|
|
|
|
"execution_count": 20,
|
2023-09-29 13:41:43 +03:00
|
|
|
|
"metadata": {},
|
|
|
|
|
|
"output_type": "execute_result"
|
|
|
|
|
|
}
|
|
|
|
|
|
],
|
|
|
|
|
|
"source": [
|
|
|
|
|
|
"file_names = [f\"RISING-{i}.csv\" for i in range(1, 6)] + [f\"BOTH-{i}.csv\" for i in range(1, 4)]\n",
|
|
|
|
|
|
"file_names"
|
|
|
|
|
|
]
|
|
|
|
|
|
},
|
|
|
|
|
|
{
|
|
|
|
|
|
"cell_type": "code",
|
2023-10-07 19:38:12 +03:00
|
|
|
|
"execution_count": 21,
|
2023-09-29 13:41:43 +03:00
|
|
|
|
"id": "cef56c78-13e9-435f-b84e-9549f195e647",
|
|
|
|
|
|
"metadata": {},
|
|
|
|
|
|
"outputs": [],
|
|
|
|
|
|
"source": [
|
|
|
|
|
|
"U_min = 900\n",
|
|
|
|
|
|
"R = 26.99e3 # Ом\n",
|
|
|
|
|
|
"\n",
|
|
|
|
|
|
"column_names = [\"Voltage\", \"Current\"]\n",
|
|
|
|
|
|
"dots_result = []\n",
|
|
|
|
|
|
"\n",
|
|
|
|
|
|
"for num, f_name in enumerate(file_names):\n",
|
|
|
|
|
|
" result = pd.DataFrame(columns=column_names)\n",
|
|
|
|
|
|
"\n",
|
|
|
|
|
|
" df = pd.read_csv(f_name, header = None)\n",
|
|
|
|
|
|
" df = df.loc[lambda x : x[4] > U_min]\n",
|
|
|
|
|
|
"\n",
|
|
|
|
|
|
" U = df[4] # В\n",
|
|
|
|
|
|
" U_cur = df[10]\n",
|
|
|
|
|
|
" I = U_cur / R # А\n",
|
|
|
|
|
|
"\n",
|
|
|
|
|
|
" result = pd.concat([result, pd.concat(dict(zip(column_names, [U, I])), axis=1)], ignore_index=True)\n",
|
|
|
|
|
|
" \n",
|
|
|
|
|
|
" result[\"Current\"] -= min(result[\"Current\"])\n",
|
|
|
|
|
|
" dots_result.append(result.groupby(\"Voltage\").mean().reset_index())"
|
|
|
|
|
|
]
|
|
|
|
|
|
},
|
|
|
|
|
|
{
|
|
|
|
|
|
"cell_type": "code",
|
2023-10-07 19:38:12 +03:00
|
|
|
|
"execution_count": 22,
|
2023-09-29 13:41:43 +03:00
|
|
|
|
"id": "ec61c33d-578b-45e4-9146-049724620a32",
|
|
|
|
|
|
"metadata": {},
|
|
|
|
|
|
"outputs": [
|
|
|
|
|
|
{
|
|
|
|
|
|
"name": "stdout",
|
|
|
|
|
|
"output_type": "stream",
|
|
|
|
|
|
"text": [
|
|
|
|
|
|
"0: k= -15.31\tb= -11.32\n",
|
|
|
|
|
|
"1: k= -10.81\tb= -14.96\n",
|
|
|
|
|
|
"2: k= -15.25\tb= -11.58\n",
|
|
|
|
|
|
"3: k= -10.82\tb= -15.03\n",
|
|
|
|
|
|
"4: k= -7.69\tb= -17.53\n",
|
|
|
|
|
|
"5: k= -7.45\tb= -17.89\n",
|
|
|
|
|
|
"6: k= -7.95\tb= -17.49\n",
|
|
|
|
|
|
"7: k= -7.96\tb= -17.50\n"
|
|
|
|
|
|
]
|
|
|
|
|
|
},
|
|
|
|
|
|
{
|
|
|
|
|
|
"data": {
|
2023-09-29 13:50:41 +03:00
|
|
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAABE0AAAL1CAYAAAAy1hoDAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAB7CAAAewgFu0HU+AAEAAElEQVR4nOzdeXxU9b3/8deZfZLJZCbJTFYCSUggEAIkoGyKCLhU3Fqx9kqtbdXo7e399ba9cr3tVbS9t8XW9ra2V8EutsVqgdaquBUQFwQUCAiRQAIJSyB7ZpJMMvuc3x+RwJBt2CN+no+HjybnfM8530lJYN75fj8fRVVVFSGEEEIIIYQQQggRRXOxJyCEEEIIIYQQQggxHEloIoQQQgghhBBCCNEPCU2EEEIIIYQQQggh+iGhiRBCCCGEEEIIIUQ/JDQRQgghhBBCCCGE6IeEJkIIIYQQQgghhBD9kNBECCGEEEIIIYQQoh+6iz0B8dnS0tLCm2++yahRozCbzRd7OkIIIYQQQgghLnFer5eDBw9y7bXXkpKSclrXSmgiLqg333yTRYsWXexpCCGEEEIIIYT4jFmxYgV33nnnaV0joYm4oEaNGgX0/GEtLCyMOldZWcmiRYv6PTccDOf5ydzO3HCe33CeGwzv+cncztxwnt9wnhsM7/nJ3M7ccJ7fcJ4bDO/5ydzO3HCe33CeGwzv+Q3nucHZz+/49cffj54OCU3EBXV8S05hYSElJSX9jhns3HAwnOcncztzw3l+w3luMLznJ3M7c8N5fsN5bjC85ydzO3PDeX7DeW4wvOcncztzw3l+w3luMLznN5znBmc/vzMpESGFYIUQQgghhBBCCCH6IaGJGDbS09N55JFHSE9Pv9hT6ddwnt9wnttwN5y/dsN5bjC85zec5zbcDeev3XCeGwzv+Q3nuQ13w/lrN5znBsN7fsN5bsPdcP7aDee5wfCe33Ce28WmqKqqXuxJiM+O8vJySktLo/aipaenyzfnJez4/+fbt28f1kv9hLhQ5HtCiGjyPSFENPmeEKKvM/m+qK+vp76+HjhR0+RMvq+kpom4KE7uoPPII4+wZMmSizcZIYQQQgghhBCXlGXLlvHoo4+e9X0kNBEXxakrTYQQQgghhBBCiHOlrKyMm266CTix0uRMSGgiLorhXpVZnDuyP1KIaPI9IUQ0+Z4QIpp8TwjR15l8X5yrMhBS00RcULJHUwghhBBCCCHEhXQ270Ole44QQgghhBBCCCFEPyQ0EUIIIYQQQgghhOiHhCZCCCGEEEIIIYQQ/ZDQRAghhBBCCCGEEKIfEpoIIYQQQgghhBBC9ENCEyGEEEIIIYQQQoh+SGgihBBCCCGEEEII0Q8JTYQQQgghhBBCCCH6IaGJEEIIIYQQQgghRD8kNBFCCCGEEEIIIYToh4QmQgghhBBCCCGEEP2Q0EQIIYQQQgghhBCiHxKaCCGEEEIIIYQQQvRDd7EnID6bKisrez9OT08nPT39Is5m+FLVMB5XBcFAK3pDMhZ7EYqivdjTEkIIIYQQQohhrb6+nvr6eiD6/efpktBEXBSLFi3q/fiRRx5hyZIlF28yw5SraSN1VU8R8DX1HjOYUskquB+7c9ZFnJkQQgghhBBCDG/Lli3j0UcfPev7SGgiLooVK1ZQWFgIIKtM+uFq2kjNrscANep4wNdIza7HyC1+WIITIYQQQgghhBhAWVkZN910E9Cz0uTkX9yfDglNxEVRWFhISUnJxZ7GeXWmW2tUNUxd1VOcGpicNIK6qqewOabLVh0hhBBCCCGE6Me5KgMhoYkQ58HZbK3xuCp6r4ug0KArpFuxE6e6SAtVokEl4Gui5dg/SEmfj6KJ7dtY6qMIIYQQQgghxOmR0ESIc+z41poI0KAbdyLw8FVSs+sxsgrux2BMwdddh6+7Dn/XETJGfw1r0mQAgoFWAGr1U9lsvhuP1tF7b0u4ieneP5AT3Mrhyp9xZO+TFM18FoPJ2TtGVVUURekzJ6mPIoQQQgghhBCnR0ITIc6h41travVTBgw8qHqqz3VeT21vaKI3JFOrn8ra+O+AEt0V3KN1sjb+O8zveoKc4FYUjR69MSVqzLEDv8PV+C5xCaMxJ+QRCftpOPg8Uh9FCCGEEEIIIU6PhCZCnCKiqlR7QrQHVRL1CvkWHZpTVm4MxOOqYF94ZEyBB4BfiaNZm4ery0xNqx9fWMUbGs078f/S5/peiobNcV9losGCTqtFOWVcd+d+/N5j+L3HcDW9O8SMpT6KEEIIIYQQQgxEQhMhTrLDFWBlXTdtgROrMpIMCrdnxTHZbugzPhz20tlajrt5M4FgN8GkebwXVzZ44GG+m4m2OFLSruCYMoo/HLFCN3Cw+6Rx5kHn6dGkEM75LvnWvnNSFC2KxoAaCcT0mgO+JjyuChKSJsY0XgghhBBCCCE+KyQ0EeITO1wBltV4UIleVdIWUFlW46Es18JkuwGft4XO1i24mzfT6dpBJBJitfUnuDWZqM1aGCAvOc6jdeBJuY381LF0e8NAxxnN972WAK/W+5mSZGCyTY9V3/Pg0ZN+iBoJ4+s+QtORl2k5+sqQ9zpeR0UIIYQQQgghxAkSmghBz5acFw65UOm7cgNAReGZmnZMqofU4B7mdf2i95wCqGhQT2N7S9A0CoBEvcL1aSZMWgWTBkxahdZAhJeP+Ya8x8GuMC2BCFWeEM8fhoIEHVPsPQFKgl6L2TKKpNTZvaHJQJ14oKeOSmfbRzQc+gspGdeT6JiGRqOP+fUIIYQQQgghxKVIQhMhgKrOAO5w/4HJcWF0dCk2XNqs3mN6YwqJKdPIUaw0RzTEaxX2ecJDPi/R0BOwWHQabsmM3ooTUVU2tgRoC0QGvN6uV3oDD+gp8bqvM8S+zp4AZUyCjilJBiYmjsNgcrIvPHLAwrRjtIex2Is4+PHjdLRupaN1Kzq9jeT0eSRnXoc5fmTUs6V1sRBCCCGEEOKzQkITIYBG92EgZchx2oiPeCVE2qgvY3dOx5wwGkVRuP+T8xFV5XsVHYMGHskGDfmWgb/1NIrC7VlmltV0ndLvpocCfHFEHJNseuq8Yba7gmxzBWj29zwzAlR2hqjsDPEcMCf1Yda3pwxYmDY1pQ7oKSB7XCjopvHwahoPryY+cRwpGddjT51NR9t2aV0shBBCCCGE+MyQ0EQIIF51EUtocpd9L9NGzwFK+z0fS+CxMMs8ZDeeyXYDZbmwss4bFcAkGzQszDL3FqUdEadjRJyOmzNMHDkeoLQFaPnkGgUo96aC0t9sAEXDa52jmImGcdOeobNtJy3HXsfd9D6qGgSgq30PXe17OLzvyX6Ly0rrYiGEEEIIIcSlSkITIYB8ixFLY3PU9pVTWcJNFNqThrxXrIFHLPeZaNPH1P5YURSy43Rkx+m4JcPEoe4w210BGn0RPmoPDvqctkCEak+IMQl6rMklWJNLCAU6aGtYT8ux1/F6agGG6MZz5q2LZbuPEEIIIYQQYriS0EQIwJpUxBX7fsrrmnv6bxesRrgi9BLWpO/EdL/TCTwGo1EUxiScXkFWRVEYFa9jVLyOD9sCQ4YmAO3B6JUoOoMVZ/atOEbcQndHFcdqV9DRsmXQewR8Tezd+v+ItxZgNKdjMKdhNKdjNKej1cX3e42raaNs9xFCCCGEEEIMWxKaCAEoipbZo2cSqvwZm81f6adg6h+ZXXj9aa2AOJPA41ThiEpFbTOtnT6SE0wU5TjQamIPXhL1sY21DvCTQFEU4hPHkJw2Z8jQBKC7Yx/dHfv6HNfqraSNXEjaqDt6j7maNlKz6zE4ZSOTbPcRQgghhBBCDBcSmoiLorKysvfj9PR00tPTL+Jsetids5gL5Fc9xuFQcm9r3mxdG9mFZRf8DfzGiqM8tWYnTe7u3mOptjjuXzCJWUWZMd0j36IjyaAZtDAtwMaWAKPi9Zi0/YcsekNy7BPvRzjYAScFTqoa5si+X3NqYHLCmW/3EUIIIYQQQoj6+nrq6+uB6Pefp0tCE3FRLFq0qPfjRx55hCVLlly8yZzE7pyFzTGd3HNQY+NsVolsrDjKY89tQj0lU2h
|
2023-09-29 13:41:43 +03:00
|
|
|
|
"text/plain": [
|
|
|
|
|
|
"<Figure size 1200x800 with 1 Axes>"
|
|
|
|
|
|
]
|
|
|
|
|
|
},
|
|
|
|
|
|
"metadata": {},
|
|
|
|
|
|
"output_type": "display_data"
|
|
|
|
|
|
}
|
|
|
|
|
|
],
|
|
|
|
|
|
"source": [
|
|
|
|
|
|
"opt_arr = []\n",
|
|
|
|
|
|
"for num, (f_name, line) in enumerate(zip(file_names, dots_result)):\n",
|
|
|
|
|
|
" x = 1e3 / line[\"Voltage\"]\n",
|
|
|
|
|
|
" y = np.log(line[\"Current\"]) - 2 * np.log(line[\"Voltage\"])\n",
|
|
|
|
|
|
" \n",
|
|
|
|
|
|
" popt, pcov = curve_fit(line_func, x, y)\n",
|
|
|
|
|
|
" print(f\"{num}: k= {popt[0]:.2f}\\tb= {popt[1]:.2f}\")\n",
|
|
|
|
|
|
" opt_arr.append(popt)\n",
|
|
|
|
|
|
" \n",
|
2023-09-29 13:50:41 +03:00
|
|
|
|
" plt.plot(x, y, ls=\"--\", marker=\".\", markersize=5, color=f\"C{num}\", label=f_name)\n",
|
|
|
|
|
|
" # plt.plot(x, line_func(x, *popt), lw=0.5, alpha=.5, color=f\"C{num}\")\n",
|
2023-09-29 13:41:43 +03:00
|
|
|
|
"\n",
|
|
|
|
|
|
"plt.xlabel(r\"$\\dfrac {10^3} U$, у.е.\")\n",
|
|
|
|
|
|
"plt.ylabel(r\"$\\log \\dfrac I {U^2}$, у.е.\")\n",
|
|
|
|
|
|
"\n",
|
|
|
|
|
|
"plt.legend()\n",
|
|
|
|
|
|
"plt.savefig(\"separ-NORDGAME-FAULER.png\")\n",
|
|
|
|
|
|
"plt.show()"
|
|
|
|
|
|
]
|
|
|
|
|
|
},
|
2023-09-29 13:25:15 +03:00
|
|
|
|
{
|
|
|
|
|
|
"cell_type": "code",
|
2023-10-07 19:38:12 +03:00
|
|
|
|
"execution_count": 23,
|
2023-09-29 13:41:43 +03:00
|
|
|
|
"id": "285e6509-6d8b-4235-8fc8-1134142e685d",
|
|
|
|
|
|
"metadata": {},
|
|
|
|
|
|
"outputs": [],
|
|
|
|
|
|
"source": [
|
|
|
|
|
|
"opt_arr = np.array(opt_arr)\n",
|
|
|
|
|
|
"k_arr = opt_arr.T[0]\n",
|
|
|
|
|
|
"b_arr = opt_arr.T[1]"
|
|
|
|
|
|
]
|
|
|
|
|
|
},
|
2023-10-07 19:38:12 +03:00
|
|
|
|
{
|
|
|
|
|
|
"cell_type": "code",
|
|
|
|
|
|
"execution_count": 53,
|
|
|
|
|
|
"id": "90c41f99-6ee6-43e1-abd5-9957a1df13ff",
|
|
|
|
|
|
"metadata": {},
|
|
|
|
|
|
"outputs": [
|
|
|
|
|
|
{
|
|
|
|
|
|
"name": "stdout",
|
|
|
|
|
|
"output_type": "stream",
|
|
|
|
|
|
"text": [
|
|
|
|
|
|
"\\begin{tabular}{rrrrrrrr}\n",
|
|
|
|
|
|
"\\hline\n",
|
|
|
|
|
|
" -15.31 & -10.81 & -15.25 & -10.82 & -7.69 & -7.45 & -7.95 & -7.96 \\\\\n",
|
|
|
|
|
|
" -11.32 & -14.96 & -11.58 & -15.03 & -17.53 & -17.89 & -17.49 & -17.50 \\\\\n",
|
|
|
|
|
|
"\\hline\n",
|
|
|
|
|
|
"\\end{tabular}\n"
|
|
|
|
|
|
]
|
|
|
|
|
|
}
|
|
|
|
|
|
],
|
|
|
|
|
|
"source": [
|
|
|
|
|
|
"from tabulate import tabulate\n",
|
|
|
|
|
|
"print(tabulate(opt_arr.T, tablefmt=\"latex\", floatfmt=\".2f\"))"
|
|
|
|
|
|
]
|
|
|
|
|
|
},
|
2023-09-29 13:41:43 +03:00
|
|
|
|
{
|
|
|
|
|
|
"cell_type": "markdown",
|
|
|
|
|
|
"id": "4b8623ef-91d3-455d-9093-8c2a60201356",
|
|
|
|
|
|
"metadata": {},
|
|
|
|
|
|
"source": [
|
2023-09-29 13:50:41 +03:00
|
|
|
|
"### Первый и второй случаи (b от ln -k)"
|
2023-09-29 13:41:43 +03:00
|
|
|
|
]
|
|
|
|
|
|
},
|
|
|
|
|
|
{
|
|
|
|
|
|
"cell_type": "code",
|
2023-10-07 19:38:12 +03:00
|
|
|
|
"execution_count": 54,
|
2023-09-29 13:41:43 +03:00
|
|
|
|
"id": "5af98fb6-4675-493e-aad6-19a8810a096a",
|
|
|
|
|
|
"metadata": {},
|
|
|
|
|
|
"outputs": [
|
|
|
|
|
|
{
|
|
|
|
|
|
"name": "stdout",
|
|
|
|
|
|
"output_type": "stream",
|
|
|
|
|
|
"text": [
|
|
|
|
|
|
"k= 8.97\tb= -36.05\n"
|
|
|
|
|
|
]
|
|
|
|
|
|
},
|
|
|
|
|
|
{
|
|
|
|
|
|
"data": {
|
2023-09-29 13:50:41 +03:00
|
|
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAABCIAAALVCAYAAAAVlkDIAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAB7CAAAewgFu0HU+AAB/IklEQVR4nO39a3Rb933n/35g2Y4tx9aGZNmBbxEBOw5sNZYBKld30lhA3JnxsBcDUmfC5slERNJHZ62JiSjrtCQz51QGm67zrA6gPMowM5XA5D+d+r/+0wBK2sZjpxUBy61j2rEBylfEskVu36/iPg8YwOAFIEhu7I3L+7UWl3DZ2PtLbALC/uC3f1+PZVmWAAAAAAAAHHCB2wUAAAAAAID+QRABAAAAAAAcQxABAAAAAAAcQxABAAAAAAAcQxABAAAAAAAcQxABAAAAAAAcQxABAAAAAAAcQxABAAAAAAAcQxABAAAAAAAcQxABAAAAAAAcQxBho3g8rsnJSdceDwAAAABAp7vQ7QK6mWmaKpfLOn78uDKZjEzTlN/vd+zxAAAAAAB0G0ZEbJLH49HAwICSyaQCgcCGA4StPh4AAAAAgG7EiIhNsixr2fVsNuvo4wEAAAAA6EaMiAAAAAAAAI4hiAAAAAAAAI4hiAAAAAAAAI5hjoge9corr+jv/u7vtGfPHl166aVulwMAAAAA6HFvv/22zpw5o7vuuktXXnllw+UIInrU3/3d32l4eNjtMgAAAAAAfWZqakpf+cpXGt5PENGj9uzZI0n6r//1v2pgYGDV/XNzc/rTP/1TTU1NKRgM2r792dlZDQ8Pt2X97Vx3u9dP7e6sn9rdWT+1O7/udq+f2t1ZP7W7s35qd2f91O7O+qndnfV3a+0vv/yyTp06pT/90z9ddbxZPc6sHo82QhDRo6qnY/zpn/5p0+WCwaBCoVDb6mjn+qndnfVTuzvrp3Z31t+t6273+qndnfVTuzvrp3Z31k/t7qyf2t1Zf7fVPj4+romJCUmNjzfXmx6AIKLHNUq/qukYAAAAAACtSiQSuummm9YcbdHqcSZBRI9rlH75fD6NjY3J5/O5UNXWtLv2dq6/m2tvN55359fdbjzvzq+73din7q2/nXjenV93u/G8u7f+duJ5d37d7datz7vP59Odd96psbEx3XnnnZtbv9VDFhYWLEm2/YRCoZa3HYlELEnW6Ojopmrf6uNXKhQKliSrUCjYsr5u2z7sxz7tTezX3sM+7T3s097Efu097NPexH7dmFafr54aEWEYhnK5nEzTtGV9fr/flvUAAAAAAIAlPRVESFIkEnG7BAAAAAAA0MAFbheA3tTN52phbezT3sR+7T3s097DPu1N7Nfewz7tTezX9vBYlmW5XUQviEajyufzGh0dVSqVcvzxKxWLRYXDYRUKhba2ggEAAAAAQGr9OJQREQAAAAAAwDEEETbb6kSZdk20CQAAAABAJyKIsEGxWNTMzIwkKZ/Pq1wuO/p4AAAAAAC6BUHEJkWjUXm9Xnk8HoXD4dpIhnK5rEAgII/HI6/Xq0Qi0ZbHAwAAAADQjXqufadTcrmcq48HAAAAAKAbMSICAAAAAAA4hiACAAAAAAA4hlMzetzs7Gztss/nk8/nc7EaAAAAAEAvqVQqqlQqkpYffzZDENHjhoeHa5fHxsY0Pj7uXjEAAAAAgJ6STqc1MTGxoccQRPS4qakpBYNBSWI0BAAAAADAVolEQkNDQ5KWRkTUfxneCEFEjwsGgwqFQm6XAQAAAADoQZuZAoDJKgEAAAAAgGMIIgAAAAAAgGM4NQMAAAAAgA50ftHSY3Mv69zr72jX5Zdo78BubbvA43ZZW0YQAQAAAABAh3nwsRd0/wOnddZ8q3bb1cZ2ff3ufbpj77UuVrZ1nJoBAAAAAEAHefCxF/SdHz60LISQpJfMt/SdHz6kBx97waXK7EEQAQAAAABAhzi/aOn+B07Lsta+37Kk+x84rfOLDRboAgQRAAAAAAB0iMfmXl41EmKls+ZbemzuZYcqsh9BBAAAAAAAHeLc6+/YulwnIogAAAAAAKBD7Lr8EluX60QEEQAAAAAAdIi9A7t1lbG96TJXG9u1d2C3QxXZjyAC6AGJREIej6fpz/T0dMPHN3pMJpNx8LfYnEQioWg0qkAgIK/Xq3g83pXbAAAAACRp2wUefePuffJ41r7f45G+fvc+bbugwQJdgCAC6AHpdFqWZalUKikUCtVuHxkZ0cLCgizLUiwWa/j4UqkkwzAkSX6/X+l0WgsLCxoZGWl36VsWDocVCoVULpdlmmbXbgMAAACoumPvtfqzr3x+1ciIq43t+rOvfF537L3WpcrscaHbBQCwj9/v16FDh1QsFiUtfZNfDRjWe5wkRSIR5XK5dpZou2pYYppm20ZwOLENAAAAoN4de6/V5265Ro/Nvaxzr7+jXZdfor0Du7t6JEQVQUSPm52drV32+Xzy+XwuVoNOFY/HNTIyolQq5XYpmxYIBHpiGwDcUSwWFY/H5ff7VwWy0WhU5XJZ2Wx22agzAADabdsFHt0WuMrtMpqqVCqqVCqSlh9/NsOpGT1ueHhY4XBY4XBY6XTa7XLQgeLxuPbv39/VIYSklkZ+dMM2ALgjnU6rXC4rn88rn8/Xbi8Wi8rn8yqXy/w/CgDAGtLpdO2Yc3h4uKXHMCKix01NTSkYDEoSoyGwSjweVzQa7Yq5IACgncLhsKSlU9UGBwdrt1dPXatfBgAAfCiRSGhoaEjS0oiIVsIIgogeFwwG+2oY6flFqyfPoWqHXg0hnJhMkgkrgd4zMjKiQqGgTCYjr9dbGwFVfb2PjIz03PslAAB22MwUAAQR6BkPPvaC7n/gtM6ab9Vuu9rYrq/fva/rZ5W1Wy+GEAQQALYqnU4rmUxqenpapVJJ0tLcMLFYbNnICAAAsDUEEegJDz72gr7zw4dkWctvf8l8S9/54UM90eLGLtFotDY5JQD0K9M0lUwmG95fnQ8ikUg0nEMnlUoxdwwAAJtAEIGud37R0v0PnF4VQlRZlnT/A6f1uVuu6fvTNKLRqPL5vOLx+JbXlclklMvlVC6XJS19qI/FYjpy5MiaH8wTiYTS6bSKxaKOHj1aG10wPz+vwcFBJZPJlr5xzGQySqfTMk1Tfr9fhmHo0KFDy5ap30a5XJZhGIpEIi1PyGn3NkzT1OHDh1Uul2utVSUpFArJ7/crm81KWhqp0miZY8eOSVJHrWflfq4e2M3MzMjv99f2cSKRUCwWa8tzYhiGvF6vTNOUYRgaHByUYRgqFou1v83qvqnuq+rtCwsLLT8Xa5mcnFQul1s2sWF1W4cOHVr1O6+n0XMSiURWvTbm5+c1PT3d8DmpPm+maWp+fr7WFWd6erp2gD0/Py9JOnTokEZHR9etr1gs1iZ0rP6uVfl8ftmIoUgksizwTCQSmpmZqf1ehmFobm5u3YP46vtVdd0r17vSRt6XDMOovScdOHCg9npf2Rmj+jtXO2cYhqGTJ0/21WmPAADYzkJPKhQKliSrUCi4XUrbnX76JSuSPLHuz+mnX3K7VEekUilL0qr9H4lELMMwavctLCxsav2FQsHy+/3WyMjIqvtGR0ctwzCsXC637PaFhQVLkuX3+63R0dFV2x4ZGbEkWdlstuF2FxYWrFAoZEmy0un0qu36/f7aNmKx2LJtFAoFyzAMy+/3N31NOLGN6vojkUjDZSKRiCXJCoVCXbOebDZrGYZhjY6OLrs9l8tZktb8e7GrBkmrtlu9XdKqv8d0Om1Jskql0obrWEs2m61ta+U6N6taSywWa7hMLBZr+Jxks9ll7wUjIyPW6Ojoqr/rVv5uFxYWas9/JBJZ872j+hpf6/muNzIyUvvd1tpnay1ffd9q9rrazPtSvepz2ezvtLpMs30CAEC/a/U4lPad6HrnXn/H1uV6UTQaVSKRWPbt7mZGReTzeYXDYcVisTXb2KVSKaVSKUWjUWUymVX3J5PJNYcyp9NpjYyMKB6Pa3JyctXjTNPUwMCAisWicrncqm9DU6lU7dvJ+fl5ZbPZZdsIhUIqFAqan59XOBxe9g22k9u
|
2023-09-29 13:41:43 +03:00
|
|
|
|
"text/plain": [
|
|
|
|
|
|
"<Figure size 1200x800 with 1 Axes>"
|
|
|
|
|
|
]
|
|
|
|
|
|
},
|
|
|
|
|
|
"metadata": {},
|
|
|
|
|
|
"output_type": "display_data"
|
|
|
|
|
|
}
|
|
|
|
|
|
],
|
|
|
|
|
|
"source": [
|
|
|
|
|
|
"x = np.log(-k_arr)\n",
|
|
|
|
|
|
"y = b_arr\n",
|
|
|
|
|
|
"\n",
|
|
|
|
|
|
"popt, pcov = curve_fit(line_func, x, y)\n",
|
|
|
|
|
|
"print(f\"k= {popt[0]:.2f}\\tb= {popt[1]:.2f}\")\n",
|
|
|
|
|
|
"\n",
|
|
|
|
|
|
"plt.plot(x, y, ls=\"\", marker=\".\", markersize=5, color=f\"C0\", label=\"Коэффициенты прямой\")\n",
|
|
|
|
|
|
"plt.plot(x, line_func(x, *popt), lw=0.5, alpha=.5, color=f\"C0\")\n",
|
|
|
|
|
|
"\n",
|
|
|
|
|
|
"plt.xlabel(r\"$b$, у.е.\")\n",
|
|
|
|
|
|
"plt.ylabel(r\"$\\log -k$, у.е.\")\n",
|
|
|
|
|
|
"\n",
|
|
|
|
|
|
"plt.legend()\n",
|
|
|
|
|
|
"plt.savefig(\"separ-1-KOEF.png\")\n",
|
|
|
|
|
|
"plt.show()"
|
|
|
|
|
|
]
|
|
|
|
|
|
},
|
|
|
|
|
|
{
|
|
|
|
|
|
"cell_type": "markdown",
|
|
|
|
|
|
"id": "9591a07f-d002-4b3c-b329-cfc1823e5080",
|
|
|
|
|
|
"metadata": {},
|
|
|
|
|
|
"source": [
|
|
|
|
|
|
"### Третий случай (b + 2 * log(-k) от exp(-b / 2))"
|
|
|
|
|
|
]
|
|
|
|
|
|
},
|
|
|
|
|
|
{
|
|
|
|
|
|
"cell_type": "code",
|
2023-10-07 19:38:12 +03:00
|
|
|
|
"execution_count": 55,
|
2023-09-29 13:41:43 +03:00
|
|
|
|
"id": "0202a50a-02d6-40dd-8d79-854d179273e4",
|
|
|
|
|
|
"metadata": {},
|
|
|
|
|
|
"outputs": [
|
|
|
|
|
|
{
|
|
|
|
|
|
"name": "stdout",
|
|
|
|
|
|
"output_type": "stream",
|
|
|
|
|
|
"text": [
|
|
|
|
|
|
"k= -0.00\tb= -6.90\n"
|
|
|
|
|
|
]
|
|
|
|
|
|
},
|
|
|
|
|
|
{
|
|
|
|
|
|
"data": {
|
|
|
|
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAABD4AAALPCAYAAACKSQBTAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAB7CAAAewgFu0HU+AACfL0lEQVR4nOz9e1hU973/f7/WDKAih+HsiArMgAYxUQdMGpP0oDNNm2an3Q3E7l+t9x+7Qvf+7752A7F/bM2+7mtH3O11/7cTSPYft7XfbxNM27TpIRk0aXOOzKjxgAoz4BEFhfF8hLn/yGaKCgg4sGB4Pq7LK7BmzVovhoGw3uvzeX+McDgcFgAAAAAAQAyymB0AAAAAAABgvFD4AAAAAAAAMYvCBwAAAAAAiFkUPgAAAAAAQMyi8AEAAAAAAGIWhQ8AAAAAABCzKHwAAAAAAICYReEDAAAAAADErDizA2B8nD17Vu+8847y8/M1a9Yss+MAAAAAAGLc1atX1d7erieffFKZmZlmx4mg8BGj3nnnHa1du9bsGAAAAACAaWbbtm364Q9/aHaMCAofMSo/P1/Sl2+44uLiQfdpbm7W2rVrh91nrDh2bBx7Kmbm2Bw7lo89FTNzbI5t9rGnYmaOzbFj+dhTMTPHHvmx+7f1X49OFhQ+YlT/9Jbi4mK5XK5h9x3JPmPFsWPj2FMxM8fm2LF87KmYmWNzbLOPPRUzc2yOHcvHnoqZOfbIjz3Z2i3Q3BQAAAAAAMQsRnzEuObm5sjHdrtddrt9Qs5rt9u1cePGcTkfx544U/X1mKrHHk9T9TWZqsceL1P19Ziqxx5PU/U1marHHi9T9fWYqsceT1P1NZmqxx4vU/X1mKrH7tfV1SW/3y/p9uvPycQIh8Nhs0Mg+vx+v0pLS2/btnHjRm3atOmufXw+37gNe8LE4HsZG/g+xg6+l7GB72Ns4PsYO/hexga+j7Gj/3tZWVmp+vr62x6bbN9fRnzEuIGNZqZStRYAAAAAMPl9//vfV1VVlaS/NzedbCh8xLjxbGIDAAAAAJjesrKyJv01J81Np7GpOG8Pg+N7GRv4PsYOvpexge9jbOD7GDv4XsYGvo+xYyp9L+nxEaOYOwcAAAAAmEiT9TqUER8AAAAAACBmUfgAAAAAAAAxi8IHAAAAAACIWRQ+AAAAAABAzKLwAQAAAAAAYhaFDwAAAAAAELMofAAAAAAAgJhF4QMAAAAAAMQsCh8AAAAAACBmUfgAAAAAAAAxi8IHAAAAAACIWRQ+AAAAAABAzKLwAQAAAAAAYlac2QEwffX2hbW/rUvnLl5TRvJMLSnIktVimB0LAAAAABBDKHzAFB/uP6mX396jztCVyLYcW6J+8vQyPb4k18RkAAAAAIBYwlQXTLgP95/Uf/zq49uKHpJ0JnRF//Grj/Xh/pMmJQMAAAAAxBoKH5NAKBTSli1bVFpaKqfTqdLSUpWWlmr79u1mR4u63r6wXn57j8LhwR8Ph6WX396j3r4hdgAAAAAAYBQofJisvr5eBQUFOnfunHbs2KFAICCfz6cdO3bo9ddfV319vdkRo2p/W9ddIz3u1Bm6ov1tXROUCAAAAAAQy+jxYaItW7aopqZGXq9Xbrf7tscaGxsj/yorK01KGH3nLl6L6n4AAAAAAAyHwodJtm/frpqaGlVXV99V9JAkr9erUCg08cHGWUbyzKjuBwAAAADAcJjqYoJQKKT169dLkjZs2DDoPjU1NaqsrFRDQ8NERht3SwqylG1LHHafHFuilhRkTVAiAAAwXqqqqmQYxrD/hutpNtRzpsJU4KqqKnk8HjmdTqWlpamiomJKngMAYgGFDxO89NJLCoVCKi8vl81mG3Qfh8Ohuro6lZeXT2y4cWa1GPqXp5fJMAZ/3DCknzy9TFbLEDsAAIApo66uTuFwWIFAQC6XK7K9srJSPT09CofDw/6tEwgEIn8r9f9t1NPTMyWmAZeWlsrlcikYDI7bKN6JOAcAxAKmupig/y6Fx+MxOYk5Hl+Sq3//4Uq9/Pae2xqd5tgS9ZOnl+nxJbkmpgMAANHmcDi0Zs0a+f1+SV+OVBjq5s+dz5Mkt9str9c7nhGjrr84EwqFxm2EykScAwBiAYWPCTawIl9WVhbZtn379shdjTVr1tx2VyQWPb4kV48unqv9bV06d/GaMpJnaklBFiM9AABAREVFhSorK1VbW2t2lDFzOp0xcQ4A5vD7/aqoqJDD4birAOzxeBQMBtXQ0BDz14/3i8LHBGtsbIx8bLPZVF9fr0AgoKqqKjkcDvn9fq1fv14Oh0OvvvrqiO6GTFVWi6GlzmyzYwAAgEmooqJCK1asUHV1tdlR7stE/C0Xy38vAtNdXV2dgsGggsGgGhsbIwtj+P3+yLVlXV2d6urqzIw56VH4mGCBQCDysd/vVyAQuO0uhsvl0o4dO1RQUKCCggK1tbXd1//Mmpubh3zMbrfLbreP+dgAAADjoaKiQh6PZ0r08gCA8VRaWirpy6l//TMG+j+/c5+J0tHRoY6OjkEfG+7600wUPibYwMZTdXV1g85Xtdlsqq2tVVVVlWpqau6rerd27dohH9u4caM2bdo05mMDAICR6e0LM71zhGK16DERzUdpcArEnsrKSvl8PtXX1ystLS1yU7z/572ysnLCf1/W1dXpxRdfnNBz3i8KHyYabh5W/xCm+vp6VVVVjXnO1rZt21RcXDzoY4z2AABg/H24/yQNvUcoFoseFDwA3K+6ujrV1NRE+kJKX/b2KS8vv23kx0SpqqrSM888M+hjzc3Nw958NwuFjwk2cNrKcKu6DHwD38+creLiYhrdAABgkg/3n9R//OpjhcO3bz8TuqL/+NXH+vcfrqT48b88Hk+kmSkATFehUEg1NTVDPt5/XVhVVTVk4+fa2tpx7f0zFVsmUPgYIBQKKS0tLWrHc7lc8vl8t23LyMiIfHyv6pzNZlMoFFJTU1PUMgEAgInR2xfWy2/vuavo0S8cll5+e48eXTx32k978Xg8amxsVEVFxX0fq76+Xl6vV8FgUNKXf9+Vl5drw4YNg14IVFVVqa6uTn6/Xy+99FJk9ER3d7fKyspUU1Mzojuq9fX1qqurUygUksPhiKzUN9DAcwSDQdlsNrnd7hGvWhPtc4RCIa1fv17BYDCy1LD05d+wDodDDQ0Nkr4ciTPUPq+++qokTarj3Pl97r+QbGpqksPhiHyPq6qqVF5ePi6vic1mU1pamkKhkGw2m8rKymSz2eT3+yPvzf7vTf/3qn97T0/PiF+LwWzZskVer/euRRXcbrfWrFlz19d8L0O9Jm63+66fje7ubm3fvn3I16T/dQuFQuru7o6s2rR9+/bIBX13d7ckac2aNSNqbuz3+yMNQPu/1n6NjY23jYhyu923FVirqqrU1NQU+bpsNtuIeiz2/77qP/adx73TaH4v2Wy2yO+k1atXR37e71y5pf9r7l/ZxWazaceOHdz0HkoYt/F6veGGhoao/PP5fHcdv6GhISwpLCkcCASGzeJwOMKSwg6HY9Rfh8/nC0saNMNkdau3L7yn9Ux4x+6j4T2tZ8K3evvMjgQAwJjtaT0Tdte8cc9/e1rPmB11QtTW1kb+Bhr494nb7Q7bbLbIYz09PWM6vs/nCzscjnBlZeVdj1VXV4dtNlvY6/Xetr2npyfyt1Z1dfVd566srAxLCjc0NAx53p6enrDL5QpLCtfV1d113oF/z5WXl992Dp/PF7bZbGGHwzHs32wTcY7+47vd7iH3cbvdYUlhl8s1ZY7T0NAQttls4erq6tu2e73esKRB3y/RyiDprvP2b5d01/uxrq5u0GuEkeQYzGiuO0aqP0t5efmQ+5SXlw/5mjQ0NNz2u6CysjJcXV191/t6JO/bnp6eyOvvdrsH/d3R/zM+2Os9UGVlZeRrG+x7Ntj+/b+3hvu5GsvvpYH6X8vh3qf9+wz3PZlIk/U6lMLHBAsEAoP+T38w/T9Mw/0yH8pkfcMN5YN9J8L/z0tv3/aH4A9fejv8wb4TZkcDAGBMduw+OqLCx47dR82OOiEGK3y43e5wQ0ND5CJ0LBd34fD
|
|
|
|
|
|
"text/plain": [
|
|
|
|
|
|
"<Figure size 1200x800 with 1 Axes>"
|
|
|
|
|
|
]
|
|
|
|
|
|
},
|
|
|
|
|
|
"metadata": {},
|
|
|
|
|
|
"output_type": "display_data"
|
|
|
|
|
|
}
|
|
|
|
|
|
],
|
|
|
|
|
|
"source": [
|
|
|
|
|
|
"x = np.exp(-b_arr / 2)\n",
|
|
|
|
|
|
"y = b_arr + 2 * np.log(-k_arr)\n",
|
|
|
|
|
|
"\n",
|
|
|
|
|
|
"popt, pcov = curve_fit(line_func, x, y)\n",
|
|
|
|
|
|
"print(f\"k= {popt[0]:.2f}\\tb= {popt[1]:.2f}\")\n",
|
|
|
|
|
|
"\n",
|
|
|
|
|
|
"plt.plot(x, y, ls=\"\", marker=\".\", markersize=5, color=f\"C0\", label=\"Коэффициенты прямой\")\n",
|
|
|
|
|
|
"plt.plot(x, line_func(x, *popt), lw=0.5, alpha=.5, color=f\"C0\")\n",
|
|
|
|
|
|
"\n",
|
|
|
|
|
|
"plt.xlabel(r\"$\\exp{-b / 2}$, у.е.\")\n",
|
|
|
|
|
|
"plt.ylabel(r\"$b + 2\\log{-k}$, у.е.\")\n",
|
|
|
|
|
|
"\n",
|
|
|
|
|
|
"plt.legend()\n",
|
|
|
|
|
|
"plt.savefig(\"separ-3-KOEF.png\")\n",
|
|
|
|
|
|
"plt.show()"
|
|
|
|
|
|
]
|
|
|
|
|
|
},
|
|
|
|
|
|
{
|
|
|
|
|
|
"cell_type": "code",
|
2023-10-07 19:38:12 +03:00
|
|
|
|
"execution_count": 56,
|
2023-09-15 23:46:38 +03:00
|
|
|
|
"id": "8dfd4820-1065-4364-abe3-498825506a74",
|
|
|
|
|
|
"metadata": {},
|
|
|
|
|
|
"outputs": [
|
|
|
|
|
|
{
|
|
|
|
|
|
"data": {
|
2023-09-29 12:46:12 +03:00
|
|
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAABBgAAALWCAYAAADs0yCXAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAB7CAAAewgFu0HU+AAEAAElEQVR4nOz9e3jc533feb9/pzkDGAxASiApiQfJEhTJsUDJthzHji3CTuKE2zZElLb8Y5snJuJ2/ex2W5NmepCUtpGppNvd9aYRKK+zz3XxShoBsbdynMQibceOHckiAcuWIko2CepEgiKBwRCY8+/0/DEABBAHzoAgQYKf13XBmpnf/Z3vbwD6j/s79/29jTAMQ0RERERERERELoO52jcgIiIiIiIiItc/FRhERERERERE5LKpwCAiIiIiIiIil00FBhERERERERG5bCowiIiIiIiIiMhlU4FBRERERERERC6bCgwiIiIiIiIictlUYBARERERERGRy6YCg4iIiIiIiIhcNhUYREREREREROSy2at9A1dCLpfj8ccfZ2hoCIBsNgvAjh072L9/P+l0etHY3t5estksDz/8MF1dXWzdunXmPYeHh/mzP/szBgYGOHz48My1hQwMDNDX1zfnte7ubvbu3XuZn05ERERERETk2mOEYRiu9k2spCNHjrBv3z7279/Prl27Zl4fGhrioYceIpfL0dfXx549exaM7+npYWBgYMkc/f39c977Yt3d3WSzWfr7++cUKD796U8zNDR0yeKEiIiIiIiIyPVmTRUYhoeH6enpYXBwcMHrR44cobu7G4DDhw+zY8eOeWOWKjDs2rWLAwcOLFkc6Onp4ciRI5w6dWrBlRLbt28nl8tx8uTJOj6RiIiIiIiIyPVhTW2R6OvrY2hoiO7ubg4fPjzv+uyCwoEDBxYsMEBthUI6nZ7ZYtHV1cX999+/5NYKqBUwBgYG2Lt376Jj9+/fT09PD/v27ePAgQP1fTARERERERGRa9yaWsHQ3d3NkSNHADh58uSCKw22bdvG8PAw6XSa8fHxedd7enrYv38/XV1dDeefXv2w2OqIaYZhLJpfRERERERE5Hq0pk6R6O3tJZ1Os2PHjkW3MUw3fLwSPRCmt1Zc6r3T6TS5XG6mGCIiIiIiIiJyvVtTBYZdu3YxPj6+4PYIqDVazOVyAEuuMFiO2cWCSxUYpq9Pb8EQERERERERud6tqQLDpTz99NNAbQXB/v37V/S9p4sFl+rTAJDJZAA4evToit6DiIiIiIiIyGpZU00el5LL5di3bx/wbhPHpQwPD8+Mz+VyZLNZduzYwf79+xeMHRsbW9Y9iYiIiIiIiKwFa77AkMvlOHbsGL29vWzdupWnnnrqkg0cH3/8cXK53LxCRE9PD9u3b+fw4cPztkE0UiyYfs/pfhCNGB0d5b/9t/9GIpEgGo0uOKa9vZ1169Y1/N4iIiIiIiJyYzp//jyjo6MLXqtUKhSLRX7jN36D9vb2Rd9jzRYYBgYG6OvrI5vNMjQ0RFdXFwcOHKj7dIiF+jj09/ezbds2tm3bxvj4+Jziw3KKBctZwfCNb3yDz372sw3HiYiIiIiIiFyO1tZW/uk//aeLXl+zBYZdu3axa9eumedDQ0P09PSQzWbp7+9ftMnjgQMHlmzSuGvXLp544gn27dtHX1/fzOvTfRWutM2bNwPwH/7Df2DLli0LjrnUCobjx4+ze/duDh06RGdn55W4TeVQjuvu/ZVDOfTvVjmU4+rlWAufQTluvBxr4TMoh3IslWOpFQynTp3i3/27fzczH13Mmi0wXKyrq4vBwUFaW1vp7u7mwIED7N27d964S50A0d3dzRNPPMHBgwc5cODAzCqGepo7roR4PA7AL//yL9e9GmMxnZ2dl/0eyqEcVzvHWvgMynHj5VgLn0E5lON6fH/lUI7r8f2VQzmuxRxDQ0P8u3/372bmo4u5oU6RSKfT7NmzB4B9+/Yt65jI2QWIY8eOzTxua2ur+z2mt0ZcraKEiIiIiIiIyJV2QxUYoLYCYdrsLQ71mr0VYnaBYrpYUE9fhel+DZdaLXE96+jo4JFHHqGjo0M5roEcV8OV/hxr5W+hv/eNlUN/b+W4Xq2F/2+slRxXw1r5Xa2Ff7dXw1r4W6ylHFfDVf0c4RoxODgYbt26Nezq6gpPnjy56LjDhw+HQAiEXV1dM6+Pj4+HXV1d4datW5eMHx8fn4nfs2fPgu87Pj6+5L1u3bo1BMK9e/fW/wGnDA4OhkA4ODjYcOxKvodcX/Q3v7Ho731j0d/7xqO/+Y1Ff+8bi/7eN57r5W9e732umRUMfX19DA8PMzQ0tKyVCUeOHGFoaIjh4WEGBgYWHTf7tIht27bNPL7//vsXHLOQ4eFhYO5qChEREREREZHr2ZopMMzuZ/DAAw8sOm56cg/MOUliutlFV1fXnNMn6o1Pp9Mz7zF7zMVm919Y7CQLERERERERkevNmikwPPzww6TTaU6ePLlkgaC/v3/mcW9v78zjrVu3zpw0sVRvhMOHDwO14yov7sC5f//+eTku9vTTTwPMNJtcDWtlL5HUT3/zG4v+3jcW/b1vPPqb31j0976x6O9941lrf3MjDMNwtW9ipezbt49cLrfoFokjR47MbEvo7++fV4g4cuQI/f39i8YPDQ2xfft20uk0p06dWvAUiJ6eHgYGBjh58uSChYpt27aRzWYZHx9v8NPNvYfBwcErfoyJiIiIiIiISL3z0DWzggHgwIEDpNNptm/fzsDAwJwTHQYGBujp6SGdTtPX17fgKocdO3awbds2uru7OXLkyMzruVyOgwcPsn37dnbs2MHg4OCiR0w+9dRT7Nixg+7u7jlbJXK53ExxY3BwcGU+sIiIiIiIiMg1Yk2tYJg2PDzMgQMHOHbsGLlcjmw2SyaTYdeuXezfv3/R4sBC8cPDw2QyGbq6unj44YeX3H4x28DAAH19fXOOr+zu7mbv3r2X9dm0gkFERERERESupnrnoWuywLCWqcAgIiIiIiIiV9MNuUVCRERERERERFaHCgwiIiIiIiIictlUYBARERERERGRy2av9g3I8hw/fnzmcUdHx5o5N1VERERERERW38jICCMjI8Dc+edSVGC4Tu3evXvm8SOPPMKjjz66ejcjIiIiIiIia0pfXx+PPfZYQzEqMFynDh06RGdnJ4BWL4iIiIiIiMiK6u3tZefOnUBtBcPsL7kXowLDdaqzs1PHVIqIiIiIiMgVsZyt+GryKCIiIiIiIiKXTSsYRERERERERGa5kK9yNjtJxfeJWhY3Z5poSUWuuxxXmwoMIiIiIiIiIsBINs9Lp0Y5N14EwA98LNMC3mF9a4J7t7TTkUld8zlWiwoMIiIiIiIicsMbHsnx3CtnyJerjOZK5EoVCAAT0vEohUqV87kiD969ga0d6Ws2x2pSgUFERERERERuaCPZPM+9cobzuSKns3ks06Q9FScasahUfcYLFXJnKmzMpHjulTPEo3bDqwyuRo7VpgKDiIiIiIiI3NBeOjVKvlzldDZPcyLCreuaprYt1NycSfDm+UlOZ/PEYzYvnRptePJ/NXKsNp0iISIiIiIiIjesC/kq58aLjOZKWKY5b+IPYJnW1Osmo7kS58aLXMhXr6kc1wIVGEREREREROSGdTY7CUCuVKE1GZ038Z9mmRatyWitb8KsuGslx7VABQYRERERERG5YVV8Hz/wIYBoZOGJ/7RoxIKgdvJDxfevqRzXAhUYRERERERE5IYVtazaigITKtWlJ/SVqg9mbaVB1Fq6UHC1c1wLVGAQERERERGRG9bNmSagdkzkeKFSW2mwAD+onfSQjkfnxF0rOa4FKjCIiIiIiIjIDaslFWF9a4L2dBw/CHjz/OS8AoAf+FOvB7Sn46xvTdCSilxTOa4FOqbyOnX8+PGZxx0dHXR0dKzi3YiIiIiIiFy/7t3SzvlckY2ZFKezeY6/NU5rMko0YlGp+lOrDgI2ZlKkYhHu3dJ+TeZYSSMjI4yMjABz559LUYHhOrV79+6Zx4888giPPvro6t2MiIiIiIjIVVItVylOFAj8ANM
|
2023-09-15 23:46:38 +03:00
|
|
|
|
"text/plain": [
|
2023-09-29 13:25:15 +03:00
|
|
|
|
"<Figure size 1200x800 with 1 Axes>"
|
2023-09-15 23:46:38 +03:00
|
|
|
|
]
|
|
|
|
|
|
},
|
|
|
|
|
|
"metadata": {},
|
|
|
|
|
|
"output_type": "display_data"
|
|
|
|
|
|
}
|
|
|
|
|
|
],
|
|
|
|
|
|
"source": [
|
2023-09-29 13:25:15 +03:00
|
|
|
|
"idx = 900\n",
|
2023-09-15 23:46:38 +03:00
|
|
|
|
"U_min = 800\n",
|
|
|
|
|
|
"\n",
|
|
|
|
|
|
"for i in range(1, 4):\n",
|
|
|
|
|
|
" df = pd.read_csv(f\"BOTH-{i}.csv\", header = None)\n",
|
2023-09-29 13:25:15 +03:00
|
|
|
|
" data = [None, None]\n",
|
2023-09-15 23:46:38 +03:00
|
|
|
|
" df = df.loc[lambda x : x[4] > U_min]\n",
|
|
|
|
|
|
" \n",
|
2023-09-29 13:25:15 +03:00
|
|
|
|
" data[0] = df.loc[lambda x : x.index < idx]\n",
|
|
|
|
|
|
" data[1] = df.loc[lambda x : x.index >= idx]\n",
|
2023-09-15 23:46:38 +03:00
|
|
|
|
" \n",
|
2023-09-29 13:25:15 +03:00
|
|
|
|
" for (i, d) in enumerate(data):\n",
|
|
|
|
|
|
" U = d[4] * 1e-3 # кВ\n",
|
|
|
|
|
|
" U_res = d[10]\n",
|
|
|
|
|
|
" R = 26.99e3 # Ом\n",
|
|
|
|
|
|
" I = U_res / R * 1e6 # мкА\n",
|
|
|
|
|
|
"\n",
|
|
|
|
|
|
" plt.plot(U, I, lw=5, alpha=0.2, ls=\"\", marker=\"o\", markersize=4, color=f\"C{i}\")\n",
|
2023-09-15 23:46:38 +03:00
|
|
|
|
" \n",
|
|
|
|
|
|
" \n",
|
|
|
|
|
|
"plt.xlabel(fr\"Напряжение, кВ\")\n",
|
|
|
|
|
|
"plt.ylabel(fr\"Сила тока, мкА\")\n",
|
|
|
|
|
|
"# plt.yscale(\"log\")\n",
|
|
|
|
|
|
"\n",
|
2023-09-29 13:41:43 +03:00
|
|
|
|
"plt.savefig(\"all-VAH-BOTH.png\")\n",
|
2023-09-15 23:46:38 +03:00
|
|
|
|
"plt.show()"
|
|
|
|
|
|
]
|
2023-10-07 19:38:12 +03:00
|
|
|
|
},
|
|
|
|
|
|
{
|
|
|
|
|
|
"cell_type": "code",
|
|
|
|
|
|
"execution_count": null,
|
|
|
|
|
|
"id": "d2dc3e7e-c3bd-4dfe-9b14-26e190916c39",
|
|
|
|
|
|
"metadata": {},
|
|
|
|
|
|
"outputs": [],
|
|
|
|
|
|
"source": []
|
2023-09-15 23:46:38 +03:00
|
|
|
|
}
|
|
|
|
|
|
],
|
|
|
|
|
|
"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.12"
|
|
|
|
|
|
}
|
|
|
|
|
|
},
|
|
|
|
|
|
"nbformat": 4,
|
|
|
|
|
|
"nbformat_minor": 5
|
|
|
|
|
|
}
|
