From 46cca564f48fab8902b562d6ea4d7346b802bb8d Mon Sep 17 00:00:00 2001
From: mirimatcode <mirimatcode@gmail.com>
Date: Fri, 21 Feb 2025 13:59:47 +0100
Subject: [PATCH] lavoro ancora alle nuove formule

---
 classificatore_singolo   | Bin 21304 -> 21392 bytes
 classificatore_singolo.c |  23 ++++++++++++++++++++++-
 percettroni.h            |  37 ++++++++++++++++++++++++++++++++++++-
 visualizzatore           | Bin 22448 -> 22448 bytes
 4 files changed, 58 insertions(+), 2 deletions(-)

diff --git a/classificatore_singolo b/classificatore_singolo
index 0719ddb268ebc821e169b759a505e1683a13f97f..3d44c6e0a095a33597167d287db3eb12d6b621ed 100755
GIT binary patch
delta 4084
zcmZu!4^R}>8Q<LlJdApGfZUxR94a3ENAMra)IiPY;-82ZT1@R23WzA8qC}lWY3DqX
z&3a0Fp_8Q3Or5FiG;ONR#a62b*u;ZLP0(tjRwF2yB^-)~=HghlzxQ?z0lPD|-|qM4
z_kQnt-}m0$iQn@zf8c9U1oP5&Cf|r0%du<OSQ5ogkGvGFy<YO=@H6-%N1FCJRG!w<
z+-cwYf{Bm$cEg5nUrV0h!L9a)yPuago4WE|@trmt87X`g?b~QE^OovFj_Wn8YV0)?
zH1wLxBe(nZM@8`bW?yaed_l;FCznS%M;ha)Q$Z6P!k>`H*kAMO$l=(OQFG&->8Hf`
z{u6sl;CJ}E_KO04*!Nm|Enl4_O2HhdN0iRhQ)WvgzbI9;h|;cBQOa);rDYAWtJSRt
zl1u4K<~dH3T`h9HA{Jc~8*he*e2biGlwFE+PV|PSN5L-nPf#G`>-(zecLnHk=jXV=
zNqkTA#(YLKKxsm(C|7fGen0~=Kmq;~?f*DF=VzvK&5Ne<TYXkbkj2%5VDAoIltCM~
zGG+q$IZ~5!L3%7o-(i|c1;$@1nszsck}Fs}bhHLW>$%ZfK@8#Q=`~q6w<~}FhGrql
znBfYT_SsZbb#GL7_b!DqNf(}~YE4b^Kk`NCw&)I|vRvPPmg_IYYN@IX$yTKNYtpjT
zrIM=!d9+zr9+fjas=C&656f<G1u<lsoPTXd;;*s9KQs?W{5nMWU41moU2n*)YlD*6
z5}IDC24oLq{W*A&I5t7ly>us}n%1Ft8In<e2e5Lx`e-tSFr&%nGwmxyGT`nCHk0nL
zwz=7i<yWW><@<Z8iu5ZY%~O>fb9QP|P}YEbJ?e~84pSfGne0x86`Vcfm4qze_C4~e
zgdE|6d!#+V;n<;h1V!%-BRmyRYWzD5J!7Qtuihg`iD}glaC)``Sx9yhe49k?=*OdI
z))=36?lL}Y;M1n@sb@yN@*OO!pgSiB99IvPoU&5$9>D&j?%mw){YdK-g<MRGwPw=*
z5R|iwSdGTEFq(rW-9$}{8~$TGwDc~S;8-l2>L!~VV}<f=a=;Nc+&@yQge%?T6UR#H
zxNe3Ml(WnX$3Ste2S`Fv-0*q2tJqBzBsqop04Yz37Y+r;+evA{Zvx~@Qmi>Sz@{0j
zA3>*K4Umqc%ZrQ|JXb8rVaX{_(nXosz_tSAn=Vxa{$>cRrEj2|Mag6KKUGC{9Tl8u
z&-pM!$sY5;@o^KG1tE{Vr98vc$r}3!PgtByi?fmJ@v*jiIu*&2OLbCh0Yq0OU~FX5
z__%8337j`6jyIfpLBmPGtL({wFROaHrq;LxI#KnWzH{J&tMeNJ$N4=${Z(%wj}f7M
z&y$eT&*7x4=u<?4jRSTlu;VJ-Omh~Ls8W>7qU;s%BqN-pOqlG9L7|Frt!5Sc9u7uY
ze4Rd;Mw!h_uI0$FWLq`0(SGxexJ=oF!o2Yox`2z>PzISn=<DjuMZ;@D(`^!?Yk*!2
zs#$Do$L8B4_rH5*o!6kfjFC5flqjPqP;k<h#GD*9aSr7x8QdCCcOBzv+Mk6U9XPej
z<S{=zfJ%SM^8jfIad0;P3%u#ZS<#@J1eSEx^jGGSan`iw7-p_zvP@x_Ou>x(o-o`g
ztFD(r{A<HlLa1SILTG@fCnSntfGTz4-E~at{-ZvgwwRL}26Lk)2TlUID|u2WmYzM0
zIkgVZRAVMku_I)9<e($OGfxX<eENeKBz8Y$rH24g#ys=_kdL!MmOZqP^>gM_AI+?1
zFc2NqYs`d~zJ;=y3>Ab9>ebug3tAm(*(Ud__V!LptY%9rAE(gSC$Ly(Tfg#O9JpP>
z+^p!4G+O&%RS7cd3Dz&~4r7j5zbw21Z@Zf{1Ftnz2758>+e+J2X2jTXqXW0-p|rr2
z7qSo=3ZcVni+8h8_+Qr11~fT2&a^)gC|zMSU3?AQ`HICU+6J>NTv>gG6g%gdf)NAg
z?7dAsaxNKtRM)X_rrxGJTnZ_Svy;nCGf7XGEqoY3)~7gy>m7uoWDDzWldmAIbP)Tb
zRN-(3Su*L<v?Fv!3Vnt|-2&JA0AN@S{x;C|=q5WFIL`Ot<VN9OYSki4Z&{gRTji$L
zH#sUw${po}J4-4Z$;DG0>o!3+HkVWum288(qHH*~qN>bMROTqGtSozd6Sq`%&2+rH
z##fyB6Q1~Ij;yAid$IU&U#n^lMCuJyZH1iFrm8l~ay#Ta$mn)eT@ASblzsu30@}ZW
zd;<9?<oOO&?SMRnoEWj<wOE3w^gx1K3i&?vkaEZ`uvzbeoQLi8&jL4xTY;nJpHS4H
zi~^98dR27<_Pv+~0DxQ^RMlL_Rgjw@<?oRbvNO%e=aNHdR^CWX;yIGErzP@p37<Y*
zIC6tbOdluQY9ou&$MJ8Fjp>PeD)FW#3hg(^U!XcfK2DD#fpjBZNFKttfsCG&$S0H8
zvl6Q_!E-(4!DWv!?-dq|j4+h&BZe6{e7|x%(C@((Zy#ij>9?Cfe;4LgKFD6$Z%>82
z87Evh<TLzIA%veRaErLn$U_(G{SLw%YUdnKZAM@l>~rX+e2AdvHW(&_GapLwUHJ6&
zk*l-DM0|qB_CE3$kC*WH;Gu8K>_@}M(I2ehr}>7@pKcHwX(TgyK3_q0W{-^?#TDRS
zp%<GVWHGUXS?!gWP82ta%f(k8y(Vgc+s>?e@n+#{+;4L=HT|@q*W2WNww1Jo+xU-|
z-zn&#1QVAvjhXvzjnv9tW_N|&XFi6FA|4%r7^d@(W3=Ep9V3=l7FhW-(*8qhcsBJR
z$AoA;k8167Da^pC(Xm3q?j+f;z6mQGsbR=S2gw##Kh~`$hFDLLL$Lmfj)bhlv+m)z
z3p(^Ia$$irCrD!j8`m^s+%T*Ijr!9R&HW(!=w*SP`oA~CKY`4A$;zjYB$u@?Q;*%^
zU?ChUGdwQQ^~Y_R9v?OA`jM_*uPaQtz7%>hXMc{{R?>*E+;lAdLvsHmtGz)-xS=Co
z#XG&OP#=4zZq1pXX<60c@iAEe1eM`Lk#?6g$BDkBvrfmqA$k?)v(C`<J9UO>x}Nqs
z^}n8^`Ll+><IgA>U&ZdN6o2Mh@TBWIG|vZgguD8GVAj$=IK*Plu={zkWucXirO)<m
z-G5fHHX`dVJob{43yZ8ph2=$?$_v+&Ru&eQl~imiS;O{SU)-Xp21{|-)}oTF1B||f
zi;wcY;F8x2LU$Z_b@_})z2vi|k%P+@3J!r>T|V9N+#hdZ&z)gu)Bl!-WS?(U?u&w8
Ziy$B5&56YRz^)+qcivn+NGvOz{{xmd>>~gG

delta 3188
zcmZWr4RBP|6~6cFhGa7&Z#RUzY<|e*2a;}3@`oTpW0w5*5(8w=I7nqQhEdu^g$|hl
zg+^SKar?S7c-X;G#({~}C|XF@b`Xe3z+vMs*u{cvM=c2jvLTTTi!y1oyq<gS+ZgDb
z*?aCi-#Ons_uO;u+s?P>vES3qBI(#8dz_an8HBgx8IVHnw^YQ5X+@x(E}*%gQZc<(
zzAMz^-N0`pnWk?0=EB4CkG!}2`QG0wyxMa5WcT@}+5#s{N0Ovvw!rf?E43Bo5E7N2
z7>LTN`=hcoDJQTyWd^0`0spK@NgB<Bi?a)pj%G5W)q`jd_Xb$fen->c&9tKAOBiEJ
z|JRJb=V?6>bp(9&PbIoM@N|}+!e)ncQCkWjs`oS18w_uY>bmN!4Xf_$VbX3_eQ~Qe
z7z|%Vq>p0s;UZ?E)rNrPrb+Ka;ct#AX?+yzx%E<Z6xQWB9N*sL7L%%Pa{^k1Rqw!Y
z=K5g*^Z)WDypdbh;YFj)k54zN-m9v6yh`<MA4J<h)i-;tZ7QBr#^kLX!gz*6JlCxZ
z42HKNa*QQ&v_K-H549xx*~GokUoos3r>t*dAS@<m&P!YR5p%$aR+XT>rVb358FL3`
zIVA$Oa-k8`)6diKvO#!#0^Z7NsH+qN6O5oWq<TXt?N@0CsZ5M25@7)`stMK0s3z3!
zsiLyj$(jjRl>eymCbQLQC$%aokA)LF9f3pnnTf52>DdUJ$uE?aMIe%&B^5*<sh~`9
zM8H*$X6?SgH$T}olci$+4S25LqBOY}D(5O4KS2z`x%HKiv+$)m%-G`y-qBTeFB3%G
zxysB^Ba~Vbo*X9HQ^icYew7IWL#8drq-t?5-zT<E@ExZs5_C+E_JZShF62Yjb@Z`5
z{CA!WPm^F=jbdyGa&oV!^`@aQ4A)C!>-pk}lf-(5$8}xfp?S}sx<;sZIfP$CEa)+_
z7_g1DnFv`O^duPbXXOZ9Bbz53=RwcRm=!@QQA0b>9k*SAh)`_c3Oqxu<HzVz5!KKe
z4$pHI?&J5t&pcGG6!chez^Z*TfE~KA9cg|md^az>gJEjxr>bw}MHZc^2{&Jsjg`x}
zF}Aiz;}k{__?jkTOJiMWSVV%*C#XJF2Zqfm9n_}gzKU~#+E>CB^JOKg8lEK<sP`{U
z)q5uV0l(uWw+PhT%ZRAFyGqB!wePpqTbQvn70>MzkBRXZ+jsZ#y3ok8s{BHKc<C)H
ztW)xzxy^&8<ed-T&PCqFuk%DMM(3a>6?+!d*os*E{#2e22Hc$x=73T2g!nlJhRV#S
zyO)b!zHO{0&eGX`cWn3uVu&72BchGUPvY}wV`p)f)=HusvelQ4-1tZ`6J8dvL(zO9
z6r*52VJEbXpY>FdSnShoJ;m2crN$1YX#T`#T*ufIU=-het?Og)ZeEeJ!`;g&)7(Bs
z<0Z?>c}XI-6y;nqisp7VFI6gP{y{#Ge|>}9I*rH7PXlj)Vz{#U3tP}N=D}9kl;!4&
z&>g+(ddP)@%#nBRM3ymg=NRo|Pu)7T*gUrgy*I^$*Nnn5&PDR;mTBSa3`4i`k=c(M
zI=+rG%;LZ&wI$OItBPkqMRA!_V}bR>h0^CEu&3B1tq;Qm6c<N8nO`DxjX?eUk4l>3
zrx9DPVX)2PNW*imb~^U;RsN|6JapfH)O9e?^lC|a;9N-qr39>HsU7T_nu@=sf9m=W
z%9k$b`XI_JL%N=T5e88fp==%2^+uE@F6;URl$-yh>w8f)UeWdUQD);Z7f@0RWfJ98
zlv#Le@`$c4MLCN3d=ljxd=a*x{1rY*2P86&+{TCPBr1mSYw!ii1K;R6dr=OfbfEn5
zwysyfKgtT}QK01tT@Uv1Li%HHmFLiX&{Cc+$wRQKJWr|~hIh*I=ofGfkxuvudl!^e
zWF|UrJ1lg11(GW=ps^x>4nT874m}I)*n7ZVk<+ms<7-1B;!3velB$ztn4Y6E;}b|0
zswPl>4R2E$hBpnDVz?FcWw@U)gr{)n8iR=gRV9f37x5GLjhMn!x5@MeGq=Gt;>lxZ
zv^NS{7iP@ZghO`}_Tz9b4*j<R#~0p?hl7a+;ph*apw2*I<^3ipR0>O6mGp7g;>w0J
za|&Gp11^O=1vaxnUx1q~CGKaqo&$QMS#$?%G%v7s;XFl>NtF`^)7snIr<eC3-d~2E
zWlHuq+N6+VqOl6W)++?LgOc;wGcYwag5UKL0kIvj_5-k0D|87AFH`7a;Hp;aEojNO
z)Y-z~f)M0I*jTN^ZACPXN9u<|Xx(L4@56sktQ`=7oPhI)o;A=9@UI(d3o%Qq4FRoD
z>c<WAJ$4nby^3?z{*BTmrIXj1#5~r|KyMaVy|LB`1F%9%jbhI;&>Pu8&5Hzw<**0Q
zN&{`f_kcD0uff|&7(nz#h}uXx={8oxdf@G+oOzewdtMRRD+byPxwQ)IV>YGd4VyE$
z0^}uz!^^0*5<3aqvzHO(MgPNl8G34!(p<blY<IsfA}lkqvDnbB!T-pt@#03=0Jb`X
zuH}6{=h*-*)cl-l_VT_94x#1_L+z&xo}N5TLwJ$na1ci8S_A7IE;a=YH5{RV6^+lB
zq)0maapeMwQGQiUxW2NEhG64L1=OZ=+yDLHhmB1OY%i4w=9qjx(9`sgB!%N3y6OQ-
Ws7|1tK<T54EU!0+SsL72p7TG9vZhV|

diff --git a/classificatore_singolo.c b/classificatore_singolo.c
index 15edcbf..e64b157 100644
--- a/classificatore_singolo.c
+++ b/classificatore_singolo.c
@@ -71,6 +71,7 @@ void main() {
 
                 //gradienti è un array bidimensionale, la prima dimensione identifica il layer, la seconda il percettrone nel layer
                 //gradienti[indice_layer][indice_percettrone]
+                //Questo vettore identifica i gradienti dei percettroni
                 double **gradienti = (double**)malloc(sizeof(double*) * NUM_LAYERS);
 
                 //Alloco la dimensione per ogni layer
@@ -88,7 +89,27 @@ void main() {
                     - gradiente dell'errore retropropagato = peso del ne
                     
                 */
-                gradienti[NUM_LAYERS-1][0] = (output_corretto - sigmoidi[NUM_LAYERS-1][0]);
+                
+                double gradiente_errore = (output_corretto - sigmoidi[NUM_LAYERS-1][0]);
+                double derivata_sigmoide_out = sigmoidi[NUM_LAYERS-1][0] * (1 - sigmoidi[NUM_LAYERS-1][0]);
+                
+                gradienti[NUM_LAYERS-1][0] = gradiente_errore * derivata_sigmoide_out;
+
+                //Ricorda di partire dal penultimo layer in quanto l'ultimo è già fatto
+                discesa_gradiente(rete_neurale, sigmoidi, gradienti);
+
+                /* A questo punto ho tutti i gradienti dei percettroni, non mi resta che trovare i gradienti dei pesi e correggerli
+                */
+               
+               //Applico la correzione dal penultimo layer andando indietro fino al secondo (il primo si fa diverso)
+               for(int indice_layer = NUM_LAYERS - 2; indice_layer > 0; indice_layer--) {
+                    //Applico la correzione a tutti i percettroni del layer dal primo a seguire
+                    for(int indice_percettrone = 0; indice_percettrone <= rete_neurale.layers[indice_layer].size; indice_percettrone++) {
+                        correggi_pesi_percettrone();
+                    }
+                }
+                
+                //gradienti[NUM_LAYERS-1][0] = (output_corretto - sigmoidi[NUM_LAYERS-1][0]);
                 errore_totale += gradienti[NUM_LAYERS-1][0];
 
                 correggi_layer_interni(&rete_neurale, gradienti, sigmoidi);
diff --git a/percettroni.h b/percettroni.h
index 4012f68..ff653d4 100644
--- a/percettroni.h
+++ b/percettroni.h
@@ -44,6 +44,9 @@ double *funzioni_attivazione_layer_double(Layer, double*);
 void correggi_layer_interni(ReteNeurale*, double**, double**);
 void correggi_layer_input(Layer*, double**, double**, byte*, int);
 
+double calcola_gradiente_layer(ReteNeurale, int, int, double**);
+void discesa_gradiente(ReteNeurale, double**, double**);
+
 int previsione(double);
 
 void salvaReteNeurale(const char*, ReteNeurale*);
@@ -211,7 +214,6 @@ void correggi_layer_interni(ReteNeurale *rete, double **gradienti, double **sigm
         }
     }
 }
-
 //Questa funzione prende tutti i parametri della precedente + gli input passati dal dataset per correggere il layer di ingresso
 void correggi_layer_input(Layer *layer, double **gradienti, double **sigmoidi, byte *inputs, int n_layers) {
     //L'indice del layer d'ingresso che prende byte per input
@@ -227,6 +229,39 @@ void correggi_layer_input(Layer *layer, double **gradienti, double **sigmoidi, b
     }
 }
 
+
+void discesa_gradiente(ReteNeurale rete, double **sigmoidi, double **gradienti) {
+    //For che scorre i layer dal penultimo al primo QUINI SIZE -2
+    for(int indice_layer = rete.size -2; indice_layer >= 0; indice_layer--) {
+        //printf("Mi trovo nel layer %d, ho %d percettroni\n", indice_layer, rete.layers[indice_layer].size);
+
+        //For che scorre i percettroni del layer partendo dal primo
+        //Per ogni percettrone mi devo prendere il gradiente disceso dal livello sopra e moltiplicarlo per la derivata di attivazione
+        for(int indice_percettrone = 0; indice_percettrone < rete.layers[indice_layer].size; indice_percettrone++) {
+            
+            double derivata_attivazione = sigmoidi[indice_layer][indice_percettrone] * (1 - sigmoidi[indice_layer][indice_percettrone]);
+            //Passo anche l'indice del percettrone perchè corrisponde all'indice del peso del livello sopra
+            double gradiente_disceso = calcola_gradiente_layer(rete, indice_layer + 1, indice_percettrone, gradienti);
+            
+
+            gradienti[indice_layer][indice_percettrone] = gradiente_disceso * derivata_attivazione;
+        }
+    }
+}
+
+double calcola_gradiente_disceso(ReteNeurale rete, int livello, int indice_peso, double **gradienti) {
+    //printf("Qui ci arrivo\n");
+    double sommatoria = 0.0;
+    //printf("Layer %d: N_percettroni: %d\n", livello, rete.layers[livello].size);
+    //Calcolo la sommatoria dei gradienti dei percettroni per i pesi
+    for(int indice_percettrone = 0; indice_percettrone < rete.layers[livello].size; indice_percettrone++) {
+        sommatoria += (gradienti[livello][indice_peso] * rete.layers[livello].percettroni[indice_percettrone].pesi[indice_peso]);
+    }
+
+    return sommatoria;
+}
+
+
 //Una volta finito il ciclo delle epoche viene salvato lo stato della rete neurale
 void salvaReteNeurale(const char *filename, ReteNeurale *rete) {
     FILE *file = fopen(filename, "wb");
diff --git a/visualizzatore b/visualizzatore
index cf86a31fbe13b88a2c5a5897ae12250913b06288..9395e1eb87302ef57172f4727ef35397054da250 100755
GIT binary patch
delta 54
zcmdn6o^b;ZtzZ@rp7Z6ytv>PIqYV$1I_dt@OLW@Dzj+UHr!A|gfuWJb<c&;nlNSWA
KY|e3X5d#3Lm==Km

delta 54
zcmdn6o^b;ZtzZ_JQT6TMM7dRM;lYKPsj6@19W0&rXY(HBPFq$J14~1L$s3vECNBtJ
K*_`9(A_f4eoEFpo