🔄 synced local 'matlabfiles/' with remote 'matlabMfiles/'

Author: UniprJRC

matlabfiles/capAlgebra/normaEucl.m

@@ -1,6 +1,3 @@
-% Numeri dalla Chi2 con 5 gradi di libertà.
-% Per replicabilità dei risultati porre $n=100000$, $p=3$ e fissare il seed
-% dei numeri casuali a 25.
 rng(25)
 n=100000;
 p=3;
@@ -19,7 +16,7 @@
 
 Xtilde=X-medie;
 Z=zscore(X);
-%% Calcolo delle norme richieste tramite ciclo for
+%% Calcolo varianze e sigma tramite le norme
 for j=1:p
     Xtildej=Xtilde(:,j);
     normXtildej=sqrt(sum(Xtildej.^2));
@@ -28,19 +25,20 @@
     diff=abs(normXtildej-normXtildejCHK);
     assert(diff<1e-12,"Errore di programmazione " + ...
         "nell'implementazione della norma euclidea")
-    diff1=abs(normXtildej-sigma(j)*sqrt(n-1));
+    % Verifica uguaglianza varianze
+    diff1=abs(normXtildej^2/(n-1)-varianze(j));
     assert(diff1<1e-12,"Errore di programmazione " + ...
-        "nell'implementazione della norma euclidea")
+        "nell'implementazione della varianza tramite la norma")
     Zj=Z(:,j);
     normZj=sqrt(sum(Zj.^2));
     % Di seguito chiamiamo direttamente la funzione norm
     normZjCHK=norm(Zj);
     diff2=abs(normZj-normZjCHK);
     assert(diff2<1e-12,"Errore di programmazione " + ...
         "nell'implementazione della norma euclidea")
-    diff3=abs(normZj-sqrt(n-1));
+    diff3=abs(normZj/sqrt(n-1)-1);
     assert(diff3<1e-12,"Errore di programmazione " + ...
-        "nell'implementazione della norma euclidea")
+        "nell'implementazione di sigma tramite la norma")
 end
 
 

matlabfiles/capAlgebraBase/sistemalin.m

@@ -16,4 +16,4 @@
 x=A\b;
 disp("Soluzione del sistema di equazioni lineari " + ...
     "tramite istruzione A\b")
-disp(x)
+disp(x)

matlabfiles/capGraficiMult/ballonPLOT.m

@@ -16,7 +16,7 @@
     'Repairs' 'Holidays'};
 colslab={'Wife'	'Alternating'	'Husband'	'Jointly'};
 tableN=array2table(N,'VariableNames',colslab,'RowNames',rowslab);
-% In this case a table is supplied
+%    In this case a table is supplied
 balloonplot(tableN);
 
 % print -depsc balloonPLOT.eps;

matlabfiles/capGraficiMult/graf2Dpersonalizzati.m

@@ -4,12 +4,8 @@
 X=stars{:,:};
 plot(X(:,1),X(:,2),'o','Marker','p','MarkerEdgeColor','r', ...
     'MarkerFaceColor','k','MarkerSize',12)
-% Per invertire la direzione dell'asse x
-set(gca,'XDir','reverse')
-xlabel('Temperatura della superficie')
-ylabel('Intensità della luce')
-title('Diagramma Hertzsprung-Russell')
-% print -depsc diagrammaHR.eps;
+set(gca,"XDir","reverse")
+
 
 %% Aggiunta delle regioni verticali 
 % xregion([3.45; 4.6],[3.6; 4.7])