Molto lavoro per negare il modello, ora si inizia con LaTeX perché altrimenti non finiamo
This commit is contained in:
@@ -18,7 +18,7 @@
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"execution_count": 58,
|
||||
"execution_count": 20,
|
||||
"id": "f34c5b88",
|
||||
"metadata": {},
|
||||
"outputs": [],
|
||||
@@ -51,7 +51,7 @@
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"execution_count": 59,
|
||||
"execution_count": 21,
|
||||
"id": "08efb2be",
|
||||
"metadata": {},
|
||||
"outputs": [],
|
||||
@@ -86,7 +86,7 @@
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"execution_count": 60,
|
||||
"execution_count": 22,
|
||||
"id": "5494409f",
|
||||
"metadata": {},
|
||||
"outputs": [
|
||||
@@ -215,7 +215,7 @@
|
||||
"4 168.53 0.002887 "
|
||||
]
|
||||
},
|
||||
"execution_count": 60,
|
||||
"execution_count": 22,
|
||||
"metadata": {},
|
||||
"output_type": "execute_result"
|
||||
}
|
||||
@@ -226,7 +226,7 @@
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"execution_count": 61,
|
||||
"execution_count": 23,
|
||||
"id": "976d5531",
|
||||
"metadata": {},
|
||||
"outputs": [
|
||||
@@ -282,7 +282,7 @@
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"execution_count": 62,
|
||||
"execution_count": 24,
|
||||
"id": "2ad19283",
|
||||
"metadata": {},
|
||||
"outputs": [
|
||||
@@ -290,10 +290,10 @@
|
||||
"name": "stdout",
|
||||
"output_type": "stream",
|
||||
"text": [
|
||||
"[23.12238655 23.79957321 23.89802584 23.89480155 24.50322786 24.30058505\n",
|
||||
" 24.15933014 24.09836266 23.99026349 23.88530016]\n",
|
||||
"[0.22419984 0.16285763 0.10115386 0.15188559 0.31275403 0.14045171\n",
|
||||
" 0.2012065 0.36234381 0.31901124 0.61345112]\n"
|
||||
"[0.04459466 0.05636189 0.07164148 0.08906597 0.04476779 0.05643232\n",
|
||||
" 0.07206497 0.04463879 0.05651695 0.04482161]\n",
|
||||
"[0.08074652 0.11183321 0.10315038 0.2075251 0.10230673 0.09273618\n",
|
||||
" 0.20255041 0.12077527 0.21682558 0.21247745]\n"
|
||||
]
|
||||
}
|
||||
],
|
||||
@@ -305,13 +305,13 @@
|
||||
"uK = np.sqrt((1/este)**2 * uF**2 + (F / este**2)**2 * ueste**2 )\n",
|
||||
"\n",
|
||||
"\n",
|
||||
"print(K)\n",
|
||||
"print(uK)"
|
||||
"print(uF)\n",
|
||||
"print(ueste)"
|
||||
]
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"execution_count": 63,
|
||||
"execution_count": 25,
|
||||
"id": "5f59d6c9",
|
||||
"metadata": {},
|
||||
"outputs": [
|
||||
@@ -360,7 +360,7 @@
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"execution_count": 64,
|
||||
"execution_count": 26,
|
||||
"id": "7e75ec05",
|
||||
"metadata": {},
|
||||
"outputs": [],
|
||||
@@ -387,7 +387,7 @@
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"execution_count": 65,
|
||||
"execution_count": 27,
|
||||
"id": "aefe7756",
|
||||
"metadata": {},
|
||||
"outputs": [
|
||||
@@ -400,8 +400,8 @@
|
||||
"Dep. Variable: y R-squared: 1.000\n",
|
||||
"Model: OLS Adj. R-squared: 1.000\n",
|
||||
"Method: Least Squares F-statistic: 2.238e+04\n",
|
||||
"Date: Sat, 04 Apr 2026 Prob (F-statistic): 4.46e-15\n",
|
||||
"Time: 19:37:55 Log-Likelihood: -27.238\n",
|
||||
"Date: Tue, 07 Apr 2026 Prob (F-statistic): 4.46e-15\n",
|
||||
"Time: 15:51:49 Log-Likelihood: -27.238\n",
|
||||
"No. Observations: 10 AIC: 58.48\n",
|
||||
"Df Residuals: 8 BIC: 59.08\n",
|
||||
"Df Model: 1 \n",
|
||||
@@ -449,7 +449,7 @@
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"execution_count": 66,
|
||||
"execution_count": 28,
|
||||
"id": "1d42b009",
|
||||
"metadata": {},
|
||||
"outputs": [
|
||||
@@ -515,7 +515,7 @@
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"execution_count": 67,
|
||||
"execution_count": 29,
|
||||
"id": "986ff4a6",
|
||||
"metadata": {},
|
||||
"outputs": [
|
||||
@@ -577,7 +577,7 @@
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"execution_count": 68,
|
||||
"execution_count": 30,
|
||||
"id": "2d4b7144",
|
||||
"metadata": {},
|
||||
"outputs": [
|
||||
@@ -653,7 +653,7 @@
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"execution_count": 69,
|
||||
"execution_count": 31,
|
||||
"id": "e2407a04",
|
||||
"metadata": {},
|
||||
"outputs": [
|
||||
@@ -722,7 +722,7 @@
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"execution_count": 70,
|
||||
"execution_count": 32,
|
||||
"id": "32e9948f",
|
||||
"metadata": {},
|
||||
"outputs": [
|
||||
@@ -785,7 +785,7 @@
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"execution_count": 71,
|
||||
"execution_count": 33,
|
||||
"id": "bfb895c6",
|
||||
"metadata": {},
|
||||
"outputs": [
|
||||
@@ -868,7 +868,7 @@
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"execution_count": 72,
|
||||
"execution_count": 34,
|
||||
"id": "202de438",
|
||||
"metadata": {},
|
||||
"outputs": [
|
||||
@@ -942,7 +942,7 @@
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"execution_count": 73,
|
||||
"execution_count": 35,
|
||||
"id": "caf23dbe",
|
||||
"metadata": {},
|
||||
"outputs": [
|
||||
@@ -1000,7 +1000,7 @@
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"execution_count": 74,
|
||||
"execution_count": 36,
|
||||
"id": "8f5c8bb7",
|
||||
"metadata": {},
|
||||
"outputs": [
|
||||
@@ -1010,8 +1010,8 @@
|
||||
"text": [
|
||||
"u_strum_m = 0.004082\n",
|
||||
"u_strum_Dx = 0.020412\n",
|
||||
"uF_strum = 0.089066\n",
|
||||
"uK_strum = 0.662828\n"
|
||||
"uF_strum = 0.058091\n",
|
||||
"uK_strum = 0.273764\n"
|
||||
]
|
||||
}
|
||||
],
|
||||
@@ -1028,10 +1028,10 @@
|
||||
"umasse_strum = np.maximum(umasse, u_strum_m)\n",
|
||||
"\n",
|
||||
"# Worst-case scalare: prendi il massimo anche di ueste_strum\n",
|
||||
"uF_strum = np.max(np.sqrt( (g * umasse_strum)**2 + (masse * ug)**2 ))\n",
|
||||
"uDx_strum = np.max(ueste_strum)\n",
|
||||
"uF_strum = np.average(np.sqrt( (g * umasse_strum)**2 + (masse * ug)**2 ))\n",
|
||||
"uDx_strum = np.average(ueste_strum)\n",
|
||||
"\n",
|
||||
"uK_strum = np.max(np.sqrt( (1/este)**2 * uF_strum**2 + (F/este**2)**2 * uDx_strum**2 ))\n",
|
||||
"uK_strum = np.average(np.sqrt( (1/este)**2 * uF_strum**2 + (F/este**2)**2 * uDx_strum**2 ))\n",
|
||||
"\n",
|
||||
"print(f\"u_strum_m = {u_strum_m:.6f}\")\n",
|
||||
"print(f\"u_strum_Dx = {u_strum_Dx:.6f}\")\n",
|
||||
@@ -1049,7 +1049,7 @@
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"execution_count": 75,
|
||||
"execution_count": 37,
|
||||
"id": "a1dc24c9",
|
||||
"metadata": {},
|
||||
"outputs": [],
|
||||
@@ -1096,12 +1096,12 @@
|
||||
"id": "2e57c7d8",
|
||||
"metadata": {},
|
||||
"source": [
|
||||
"## Risutltati della propagazione dell'errore strumentale massimo"
|
||||
"## Risultati della propagazione dell'errore strumentale massimo"
|
||||
]
|
||||
},
|
||||
{
|
||||
"cell_type": "code",
|
||||
"execution_count": 76,
|
||||
"execution_count": 38,
|
||||
"id": "c8e264fa",
|
||||
"metadata": {},
|
||||
"outputs": [
|
||||
@@ -1110,22 +1110,22 @@
|
||||
"output_type": "stream",
|
||||
"text": [
|
||||
"RISULTATI CON ERRORE STRUMENTALE INCLUSO:\n",
|
||||
"Media pesata K = 23.95113 ± 0.66530\n",
|
||||
"Media pesata K = 23.95113 ± 0.27970\n",
|
||||
"\n",
|
||||
"RISULTATI REGRESSIONE OLS:\n",
|
||||
"Aols = 23.96871 ± 0.68192\n",
|
||||
"Bols = 0.09993 ± 2.98913\n",
|
||||
"Chi² OLS = 3.95368 | rid = 0.49421 | P = 0.13872\n",
|
||||
"Aols = 23.96871 ± 0.31720\n",
|
||||
"Bols = 0.09993 ± 2.92755\n",
|
||||
"Chi² OLS = 7.21389 | rid = 0.90174 | P = 0.48626\n",
|
||||
"\n",
|
||||
"RISULTATI REGRESSIONE Carpi:\n",
|
||||
"AC = 24.04738 ± 0.67577\n",
|
||||
"BC = -1.56983 ± 2.18409\n",
|
||||
"Chi² Carpi = 4.03420 | rid = 0.50428 | P = 0.14598\n",
|
||||
"AC = 24.04738 ± 0.30375\n",
|
||||
"BC = -1.56983 ± 2.09901\n",
|
||||
"Chi² Carpi = 7.28743 | rid = 0.91093 | P = 0.49404\n",
|
||||
"\n",
|
||||
"RISULTATI REGRESSIONE York:\n",
|
||||
"AY = 24.06480 ± 0.67582\n",
|
||||
"BY = -1.81758 ± 2.18788\n",
|
||||
"Chi² York = 4.06701 | rid = 0.50838 | P = 0.14897\n"
|
||||
"AY = 24.06480 ± 0.30386\n",
|
||||
"BY = -1.81758 ± 2.10295\n",
|
||||
"Chi² York = 7.33375 | rid = 0.91672 | P = 0.49891\n"
|
||||
]
|
||||
}
|
||||
],
|
||||
@@ -1163,8 +1163,8 @@
|
||||
"\n",
|
||||
"In generale con il Chi² vale:\n",
|
||||
"- \\~50%: Errori ottimamente stimati\n",
|
||||
"- \\<5% : Fit troppo grande\n",
|
||||
"- \\>95%: Fit troppo piccolo\n",
|
||||
"- \\<5% : Fit troppo grande (errori sovrastimati)\n",
|
||||
"- \\>95%: Fit troppo piccolo (errori sottostimati)\n",
|
||||
"\n",
|
||||
"In generale mi sembra che nei paper se è presente un Chi² si riporta la probabilità complementare (probabilità di ottenere gli stessi risultati o peggio)"
|
||||
]
|
||||
|
||||
Reference in New Issue
Block a user