{ "cells": [ { "cell_type": "markdown", "id": "3a8dfc0e", "metadata": {}, "source": [ "# Analisi della molla statica con calibro 2" ] }, { "cell_type": "markdown", "id": "2ff20b69", "metadata": {}, "source": [ "## Import delle librerie e set di variabili gloabali" ] }, { "cell_type": "code", "execution_count": 82, "id": "f34c5b88", "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import pandas as pd\n", "import scipy as sc\n", "from scipy.stats import chi2\n", "import seaborn as sns\n", "import matplotlib.pyplot as plt\n", "import matplotlib as mpl\n", "import statsmodels.api as sm\n", "\n", "\n", "g = 9.807\n", "ug = 0.001" ] }, { "cell_type": "markdown", "id": "39d2ffea", "metadata": {}, "source": [ "## Lettura dei dati e calcolo delle deviazioni standard campionarie\n", "- Lettura del CSV\n", "- Creazione del data frame\n", "- Deviazioni standard\n", "- Permutazioni e calcolo dei Delta" ] }, { "cell_type": "code", "execution_count": 83, "id": "08efb2be", "metadata": {}, "outputs": [], "source": [ "df = pd.read_csv(r'statica2.csv')\n", "\n", "def calcola_stats(df, prefix, err_arbitrario):\n", " cols = [col for col in df.columns if col.startswith(prefix)]\n", "\n", " def riga_stats(row):\n", " valori = row[cols].dropna()\n", " n = len(valori)\n", "\n", " if n == 0:\n", " return pd.Series({prefix: np.nan, f\"u{prefix}\": np.nan})\n", " elif n == 1:\n", " return pd.Series({prefix: valori.iloc[0], f\"u{prefix}\": err_arbitrario})\n", " else:\n", " media = valori.mean()\n", " sigma = valori.std(ddof=1)\n", " return pd.Series({prefix: media, f\"u{prefix}\": sigma})\n", "\n", " stats = df.apply(riga_stats, axis=1)\n", " df[prefix] = stats[prefix]\n", " df[f\"u{prefix}\"] = stats[f\"u{prefix}\"]\n", "\n", " return df\n", "\n", "df = calcola_stats(df, \"Dx\", err_arbitrario=0.01)\n", "df = calcola_stats(df, \"m\", err_arbitrario=0.0028867513)\n" ] }, { "cell_type": "code", "execution_count": 84, "id": "5494409f", "metadata": {}, "outputs": [ { "data": { "text/html": [ "
| \n", " | m1 | \n", "Dx1 | \n", "Dx2 | \n", "Dx3 | \n", "Dx4 | \n", "Dx5 | \n", "Dx6 | \n", "Dx | \n", "uDx | \n", "m | \n", "um | \n", "
|---|---|---|---|---|---|---|---|---|---|---|---|
| 0 | \n", "49.246667 | \n", "212.10 | \n", "211.64 | \n", "212.00 | \n", "212.18 | \n", "212.52 | \n", "212.04 | \n", "212.080000 | \n", "0.284816 | \n", "49.246667 | \n", "0.002887 | \n", "
| 1 | \n", "69.276667 | \n", "150.92 | \n", "150.26 | \n", "150.02 | \n", "150.16 | \n", "150.40 | \n", "150.18 | \n", "150.323333 | \n", "0.317847 | \n", "69.276667 | \n", "0.002887 | \n", "
| 2 | \n", "88.966667 | \n", "90.34 | \n", "90.34 | \n", "90.38 | \n", "90.52 | \n", "90.26 | \n", "90.28 | \n", "90.353333 | \n", "0.092664 | \n", "88.966667 | \n", "0.002887 | \n", "
| 3 | \n", "108.610000 | \n", "29.82 | \n", "30.18 | \n", "30.10 | \n", "30.20 | \n", "30.10 | \n", "30.40 | \n", "30.133333 | \n", "0.188750 | \n", "108.610000 | \n", "0.002887 | \n", "