codice parziale creato da giulia

osc_attr.ipynb

@@ -2,7 +2,7 @@
  "cells": [
   {
    "cell_type": "code",
-   "execution_count": null,
+   "execution_count": 1,
    "id": "9c25e271",
    "metadata": {},
    "outputs": [],
@@ -24,7 +24,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": null,
+   "execution_count": 2,
    "id": "e1efad4a",
    "metadata": {},
    "outputs": [],
@@ -45,7 +45,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": null,
+   "execution_count": 3,
    "id": "81d935ac",
    "metadata": {},
    "outputs": [],
@@ -68,10 +68,25 @@
   },
   {
    "cell_type": "code",
-   "execution_count": null,
+   "execution_count": 6,
    "id": "7b2c2af7",
    "metadata": {},
-   "outputs": [],
+   "outputs": [
+    {
+     "data": {
+      "application/vnd.jupyter.widget-view+json": {
+       "model_id": "a5b577c48b9e41549f4804af6032b481",
+       "version_major": 2,
+       "version_minor": 0
+      },
+      "text/plain": [
+       "Sheet(cells=(Cell(column_end=0, column_start=0, row_end=6, row_start=0, squeeze_row=False, type='numeric', val…"
+      ]
+     },
+     "metadata": {},
+     "output_type": "display_data"
+    }
+   ],
    "source": [
     "# Calcolo w0 varie cnfigurazioni\n",
     "\n",
@@ -85,7 +100,10 @@
     "param_sistema = param_sistema.calc_var(w0_1_CD, \"w0_1_CD\")\n",
     "param_sistema = param_sistema.calc_var(w0_1_rete, \"w0_1_rete\")\n",
     "param_sistema = param_sistema.calc_var(w0_2_CD, \"w0_2_CD\")\n",
-    "param_sistema = param_sistema.calc_var(w0_2_rete, \"w0_2_rete\")"
+    "param_sistema = param_sistema.calc_var(w0_2_rete, \"w0_2_rete\")\n",
+    "\n",
+    "sheet = ipysheet.from_dataframe(param_sistema.analisi_stat)\n",
+    "display(sheet)"
    ]
   },
   {
@@ -100,19 +118,146 @@
   },
   {
    "cell_type": "code",
-   "execution_count": null,
+   "execution_count": 7,
    "id": "c040f35e",
    "metadata": {},
    "outputs": [],
    "source": [
-    "config = [molla1_leggera_rete_inizio, molla1_leggera_rete_fine, molla1_leggera_rete_mid,\n",
-    "          molla1_pesante_rete_inizio, molla1_pesante_rete_fine, molla1_pesante_rete_mid,\n",
-    "          molla1_leggera_CD_inizio,   molla1_leggera_CD_fine,   molla1_leggera_CD_mid,\n",
-    "          molla1_pesante_CD_inizio,   molla1_pesante_CD_fine,   molla1_pesante_CD_mid,\n",
-    "          molla2_leggera_rete_inizio, molla2_leggera_rete_fine, molla2_leggera_rete_mid,\n",
-    "          molla2_pesante_rete_inizio, molla2_pesante_rete_fine, molla2_pesante_rete_mid,\n",
-    "          molla2_leggera_CD_inizio,   molla2_leggera_CD_fine,   molla2_leggera_CD_mid,\n",
-    "          molla2_pesante_CD_inizio,   molla2_pesante_CD_fine,   molla2_pesante_CD_mid]"
+    "# config = [molla1_leggera_rete_inizio, molla1_leggera_rete_fine, molla1_leggera_rete_mid,\n",
+    "#           molla1_pesante_rete_inizio, molla1_pesante_rete_fine, molla1_pesante_rete_mid,\n",
+    "#           molla1_leggera_CD_inizio,   molla1_leggera_CD_fine,   molla1_leggera_CD_mid,\n",
+    "#           molla1_pesante_CD_inizio,   molla1_pesante_CD_fine,   molla1_pesante_CD_mid,\n",
+    "#           molla2_leggera_rete_inizio, molla2_leggera_rete_fine, molla2_leggera_rete_mid,\n",
+    "#           molla2_pesante_rete_inizio, molla2_pesante_rete_fine, molla2_pesante_rete_mid,\n",
+    "#           molla2_leggera_CD_inizio,   molla2_leggera_CD_fine,   molla2_leggera_CD_mid,\n",
+    "#           molla2_pesante_CD_inizio,   molla2_pesante_CD_fine,   molla2_pesante_CD_mid]"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 8,
+   "id": "c9611e53",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "h = [0.51415, 0.4983, 0.485, 0.4735]\n",
+    "df = pd.DataFrame()\n",
+    "df[\"h\"] = h\n",
+    "molla2_leggera_rete_inizio = Data(df)\n",
+    "molla2_leggera_rete_inizio.analisi_stat = df\n",
+    "\n",
+    "h = [0.5121, 0.498, 0.4856, 0.47389]\n",
+    "df = pd.DataFrame()\n",
+    "df[\"h\"] = h\n",
+    "molla2_leggera_CD_inizio = Data(df)\n",
+    "molla2_leggera_CD_inizio.analisi_stat = df"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 9,
+   "id": "3014c983",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "config = [molla2_leggera_rete_inizio,\n",
+    "          molla2_leggera_CD_inizio]"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "id": "28734bca",
+   "metadata": {},
+   "source": [
+    "# Calcolo A"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 13,
+   "id": "00b5d09c",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "molla2_leggera_CD_inizio.analisi_stat[\"h1\"] = molla2_leggera_CD_inizio.analisi_stat[\"h\"][1]\n",
+    "molla2_leggera_CD_inizio.analisi_stat[\"h2\"] = molla2_leggera_CD_inizio.analisi_stat[\"h\"][2]\n",
+    "\n",
+    "molla2_leggera_rete_inizio.analisi_stat[\"h1\"] = molla2_leggera_rete_inizio.analisi_stat[\"h\"][1]\n",
+    "molla2_leggera_rete_inizio.analisi_stat[\"h2\"] = molla2_leggera_rete_inizio.analisi_stat[\"h\"][2]"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "2e3e28c0",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "# Calcolo numerico x e y\n",
+    "\n",
+    "for data in config:\n",
+    "    data.analisi_stat[\"h1\"] = data.analisi_stat[\"h\"][1]\n",
+    "    data.analisi_stat[\"h2\"] = data.analisi_stat[\"h\"][2]\n",
+    "    data.analisi_stat[\"uh1\"] = 0.001\n",
+    "    data.analisi_stat[\"uh2\"] = 0.001\n",
+    "    data.analisi_stat[\"uh\"] = 0.001\n",
+    "    data.analisi_stat[\"n\"] = [1, 2, 3, 4]\n",
+    "    data.analisi_stat[\"un\"] = 0.000000001"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "857214ff",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "h, h1, h2 = sp.symbols('h h1 h2')\n",
+    "h0 = (h*h2-h1**2) / (2*h1 + h + h2)"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 22,
+   "id": "4c030372",
+   "metadata": {},
+   "outputs": [
+    {
+     "data": {
+      "application/vnd.jupyter.widget-view+json": {
+       "model_id": "b83544b403c04d0493ef734f9dc039ab",
+       "version_major": 2,
+       "version_minor": 0
+      },
+      "text/plain": [
+       "Sheet(cells=(Cell(column_end=0, column_start=0, row_start=0, squeeze_row=False, type='numeric', value=[0.51415…"
+      ]
+     },
+     "metadata": {},
+     "output_type": "display_data"
+    },
+    {
+     "data": {
+      "application/vnd.jupyter.widget-view+json": {
+       "model_id": "ab432f011593444aa2d42ae9fe4f3601",
+       "version_major": 2,
+       "version_minor": 0
+      },
+      "text/plain": [
+       "Sheet(cells=(Cell(column_end=0, column_start=0, row_start=0, squeeze_row=False, type='numeric', value=[0.5121,…"
+      ]
+     },
+     "metadata": {},
+     "output_type": "display_data"
+    }
+   ],
+   "source": [
+    "# Calcolo numerico x e y\n",
+    "\n",
+    "for data in config:\n",
+    "    data = data.calc_var(h0, \"h0\")\n",
+    "    sheet = ipysheet.from_dataframe(data.analisi_stat)\n",
+    "    display(sheet)"
    ]
   },
   {