{ "cells": [ { "cell_type": "markdown", "id": "a80529e4", "metadata": {}, "source": [ "# Analisi dei dati con il sonar\n", "Minimo Indice:\n", "- Analisi dei dati statici\n", "- Analisi dei dati dinamici\n", " - Sonar\n", " - Cronometro" ] }, { "cell_type": "markdown", "id": "32702b8f", "metadata": {}, "source": [ "## Import delle librerie e set di variabili gloabali" ] }, { "cell_type": "code", "execution_count": 140, "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\n", "\n", "m_mol = 15.43\n", "um_mol = 0.01" ] }, { "cell_type": "markdown", "id": "fd0b8b1d", "metadata": {}, "source": [ "## Lettura dei dati e calcolo delle deviazioni standard campionarie\n", "- Lettura del CSV\n", "- Creazione del data frame\n", "- Deviazioni standard\n", "\n", "ATTENZIONE: Linea cursed ~17" ] }, { "cell_type": "code", "execution_count": 141, "id": "08efb2be", "metadata": {}, "outputs": [], "source": [ "df = pd.read_csv(r'dinamica2.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", "\n", "df = calcola_stats(df, \"w\", err_arbitrario=0.0002)\n", "df = calcola_stats(df, \"m\", err_arbitrario=0.0028867513)\n", "df = calcola_stats(df, \"c\", err_arbitrario=0.01)\n", "df = calcola_stats(df, \"a\", err_arbitrario=0.01)\n", "df = calcola_stats(df, \"t\", err_arbitrario=0.01)" ] }, { "cell_type": "code", "execution_count": 142, "id": "5494409f", "metadata": {}, "outputs": [ { "data": { "text/html": [ "
| \n", " | m1 | \n", "a1 | \n", "ua1 | \n", "w1 | \n", "uw1 | \n", "c1 | \n", "uc1 | \n", "t1 | \n", "a2 | \n", "ua2 | \n", "w2 | \n", "uw2 | \n", "c2 | \n", "uc2 | \n", "t2 | \n", "a3 | \n", "ua3 | \n", "w3 | \n", "uw3 | \n", "c3 | \n", "uc3 | \n", "t3 | \n", "a4 | \n", "ua4 | \n", "w4 | \n", "uw4 | \n", "c4 | \n", "uc4 | \n", "t4 | \n", "w | \n", "uw | \n", "m | \n", "um | \n", "c | \n", "uc | \n", "a | \n", "ua | \n", "t | \n", "ut | \n", "
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 0 | \n", "49.25 | \n", "9.7171 | \n", "0.016 | \n", "7.65650 | \n", "0.00040 | \n", "484.4550 | \n", "0.011 | \n", "15.62 | \n", "8.9110 | \n", "0.015 | \n", "7.65690 | \n", "0.00040 | \n", "484.516 | \n", "0.011 | \n", "15.58 | \n", "10.446 | \n", "0.027 | \n", "7.66030 | \n", "0.00050 | \n", "485.082 | \n", "0.019 | \n", "15.76 | \n", "8.377 | \n", "0.016 | \n", "7.65820 | \n", "0.00040 | \n", "484.752 | \n", "0.011 | \n", "15.87 | \n", "7.657975 | \n", "0.001711 | \n", "49.25 | \n", "0.002887 | \n", "484.70125 | \n", "0.284314 | \n", "9.362775 | \n", "0.908254 | \n", "15.7075 | \n", "0.133010 | \n", "
| 1 | \n", "69.28 | \n", "9.8600 | \n", "0.016 | \n", "6.55968 | \n", "0.00029 | \n", "423.3520 | \n", "0.011 | \n", "18.31 | \n", "10.3900 | \n", "0.012 | \n", "6.55891 | \n", "0.00022 | \n", "423.154 | \n", "0.009 | \n", "18.27 | \n", "10.491 | \n", "0.013 | \n", "6.56002 | \n", "0.00024 | \n", "423.697 | \n", "0.010 | \n", "18.34 | \n", "10.968 | \n", "0.019 | \n", "6.56000 | \n", "0.00030 | \n", "423.465 | \n", "0.014 | \n", "18.16 | \n", "6.559652 | \n", "0.000519 | \n", "69.28 | \n", "0.002887 | \n", "423.41700 | \n", "0.226641 | \n", "10.427250 | \n", "0.454472 | \n", "18.2700 | \n", "0.078740 | \n", "
| 2 | \n", "88.97 | \n", "11.5840 | \n", "0.014 | \n", "5.84417 | \n", "0.00020 | \n", "363.2290 | \n", "0.010 | \n", "20.27 | \n", "10.1763 | \n", "0.017 | \n", "5.84585 | \n", "0.00028 | \n", "363.354 | \n", "0.012 | \n", "20.44 | \n", "12.044 | \n", "0.018 | \n", "5.84500 | \n", "0.00026 | \n", "363.183 | \n", "0.013 | \n", "20.54 | \n", "11.224 | \n", "0.016 | \n", "5.84513 | \n", "0.00025 | \n", "363.233 | \n", "0.011 | \n", "20.49 | \n", "5.845038 | \n", "0.000689 | \n", "88.97 | \n", "0.002887 | \n", "363.24975 | \n", "0.073109 | \n", "11.257075 | \n", "0.794837 | \n", "20.4350 | \n", "0.117331 | \n", "
| 3 | \n", "108.61 | \n", "11.5420 | \n", "0.026 | \n", "5.32780 | \n", "0.00030 | \n", "303.5502 | \n", "0.019 | \n", "22.49 | \n", "8.4240 | \n", "0.017 | \n", "5.32820 | \n", "0.00030 | \n", "303.581 | \n", "0.012 | \n", "22.27 | \n", "10.501 | \n", "0.022 | \n", "5.32960 | \n", "0.00030 | \n", "303.842 | \n", "0.016 | \n", "22.55 | \n", "9.959 | \n", "0.014 | \n", "5.32822 | \n", "0.00020 | \n", "303.445 | \n", "0.010 | \n", "22.15 | \n", "5.328455 | \n", "0.000787 | \n", "108.61 | \n", "0.002887 | \n", "303.60455 | \n", "0.168669 | \n", "10.106500 | \n", "1.299853 | \n", "22.3650 | \n", "0.187172 | \n", "
| 4 | \n", "128.64 | \n", "11.5740 | \n", "0.020 | \n", "4.92663 | \n", "0.00023 | \n", "242.9620 | \n", "0.014 | \n", "24.33 | \n", "11.5920 | \n", "0.023 | \n", "4.92537 | \n", "0.00029 | \n", "242.876 | \n", "0.017 | \n", "24.38 | \n", "10.264 | \n", "0.023 | \n", "4.92430 | \n", "0.00030 | \n", "242.789 | \n", "0.017 | \n", "25.09 | \n", "9.118 | \n", "0.021 | \n", "4.92610 | \n", "0.00030 | \n", "243.115 | \n", "0.015 | \n", "24.26 | \n", "4.925600 | \n", "0.001009 | \n", "128.64 | \n", "0.002887 | \n", "242.93550 | \n", "0.138954 | \n", "10.637000 | \n", "1.188344 | \n", "24.5150 | \n", "0.386480 | \n", "