physics update

Author: NF-coder

physics/sem1/lab1/avg-brightness-calculator/experement.ipynb

@@ -36,7 +36,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 2,
+   "execution_count": 1,
    "id": "6acd30a5",
    "metadata": {},
    "outputs": [],
@@ -50,7 +50,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 3,
+   "execution_count": 10,
    "id": "564b8910",
    "metadata": {},
    "outputs": [],
@@ -709,30 +709,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 11,
-   "id": "1b091193",
-   "metadata": {},
-   "outputs": [
-    {
-     "data": {
-      "text/plain": [
-       "50"
-      ]
-     },
-     "execution_count": 11,
-     "metadata": {},
-     "output_type": "execute_result"
-    }
-   ],
-   "source": [
-    "n = 2\n",
-    "\n",
-    "len(list(filter(lambda x: 2.89 - 0.38*n <= x <= 2.89 + 0.38*n, data)))"
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": 207,
+   "execution_count": 4,
    "id": "647e0661",
    "metadata": {},
    "outputs": [],
@@ -752,7 +729,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 4,
+   "execution_count": 12,
    "id": "dbea9be6",
    "metadata": {},
    "outputs": [
@@ -771,7 +748,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 5,
+   "execution_count": 14,
    "id": "553175bd",
    "metadata": {},
    "outputs": [
@@ -802,18 +779,25 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 6,
+   "execution_count": null,
    "id": "10681818",
    "metadata": {},
    "outputs": [
+    {
+     "name": "stdout",
+     "output_type": "stream",
+     "text": [
+      "0.08132897477598496\n"
+     ]
+    },
     {
      "ename": "ValueError",
      "evalue": "could not convert string to float: ''",
      "output_type": "error",
      "traceback": [
       "\u001b[31m---------------------------------------------------------------------------\u001b[39m",
       "\u001b[31mValueError\u001b[39m                                Traceback (most recent call last)",
-      "\u001b[36mCell\u001b[39m\u001b[36m \u001b[39m\u001b[32mIn[6]\u001b[39m\u001b[32m, line 4\u001b[39m\n\u001b[32m      1\u001b[39m f = \u001b[38;5;28;01mlambda\u001b[39;00m x: max_ro*np.exp(-(x-mean)**\u001b[32m2\u001b[39m/(\u001b[32m2\u001b[39m*dispersion))\n\u001b[32m      3\u001b[39m \u001b[38;5;28;01mwhile\u001b[39;00m \u001b[38;5;28;01mTrue\u001b[39;00m:\n\u001b[32m----> \u001b[39m\u001b[32m4\u001b[39m     \u001b[38;5;28mprint\u001b[39m(f(\u001b[38;5;28;43mfloat\u001b[39;49m\u001b[43m(\u001b[49m\u001b[38;5;28;43minput\u001b[39;49m\u001b[43m(\u001b[49m\u001b[43m)\u001b[49m\u001b[43m)\u001b[49m))\n",
+      "\u001b[36mCell\u001b[39m\u001b[36m \u001b[39m\u001b[32mIn[15]\u001b[39m\u001b[32m, line 4\u001b[39m\n\u001b[32m      1\u001b[39m f = \u001b[38;5;28;01mlambda\u001b[39;00m x: max_ro*np.exp(-(x-mean)**\u001b[32m2\u001b[39m/(\u001b[32m2\u001b[39m*dispersion))\n\u001b[32m      3\u001b[39m \u001b[38;5;28;01mwhile\u001b[39;00m \u001b[38;5;28;01mTrue\u001b[39;00m:\n\u001b[32m----> \u001b[39m\u001b[32m4\u001b[39m     \u001b[38;5;28mprint\u001b[39m(f(\u001b[38;5;28;43mfloat\u001b[39;49m\u001b[43m(\u001b[49m\u001b[38;5;28;43minput\u001b[39;49m\u001b[43m(\u001b[49m\u001b[43m)\u001b[49m\u001b[43m)\u001b[49m))\n",
       "\u001b[31mValueError\u001b[39m: could not convert string to float: ''"
      ]
     }
@@ -825,6 +809,91 @@
     "    print(f(float(input())))"
    ]
   },
+  {
+   "cell_type": "code",
+   "execution_count": 24,
+   "id": "22e1e1a0",
+   "metadata": {},
+   "outputs": [
+    {
+     "name": "stdout",
+     "output_type": "stream",
+     "text": [
+      "1 - from 2.74 to 3.05: 40 0.80\n",
+      "2 - from 2.58 to 3.21: 46 0.92\n",
+      "3 - from 2.42 to 3.37: 48 0.96\n"
+     ]
+    }
+   ],
+   "source": [
+    "for n in [1,2,3]:\n",
+    "    _ = len(list(filter(lambda x: mean - std*n <= x <= mean + std*n, data)))\n",
+    "\n",
+    "    print(f\"{n} - from {mean - std*n:.2f} to {mean + std*n:.2f}: {_} {_/50:.2f}\")"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 28,
+   "id": "1f848cc7",
+   "metadata": {},
+   "outputs": [
+    {
+     "name": "stdout",
+     "output_type": "stream",
+     "text": [
+      "0.022271515623398092\n",
+      "0.04476574640303016\n"
+     ]
+    }
+   ],
+   "source": [
+    "mean_std = np.sqrt(1/(50*49)*sum(map(lambda x: (x-mean)**2, data)))\n",
+    "\n",
+    "print(mean_std)\n",
+    "\n",
+    "trust_range = 2.01 * mean_std\n",
+    "\n",
+    "print(trust_range)"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 30,
+   "id": "b6fe708c",
+   "metadata": {},
+   "outputs": [
+    {
+     "name": "stdout",
+     "output_type": "stream",
+     "text": [
+      "0.04525943543024871\n"
+     ]
+    }
+   ],
+   "source": [
+    "ttl_pogreshnost = np.sqrt(trust_range**2 + (2/3*0.01)**2)\n",
+    "print(ttl_pogreshnost)"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 35,
+   "id": "c06117ce",
+   "metadata": {},
+   "outputs": [
+    {
+     "name": "stdout",
+     "output_type": "stream",
+     "text": [
+      "1.56%\n"
+     ]
+    }
+   ],
+   "source": [
+    "print(f\"{ttl_pogreshnost/mean*100:.2f}%\")"
+   ]
+  },
   {
    "cell_type": "code",
    "execution_count": 8,
@@ -1040,7 +1109,7 @@
    "name": "python",
    "nbconvert_exporter": "python",
    "pygments_lexer": "ipython3",
-   "version": "3.12.4"
+   "version": "3.13.7"
   }
  },
  "nbformat": 4,