1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45
| clc;
disp('La matrice de données X')
X=[6 6 5 5.5 ; 8 8 8 8 ; 6 7 11 9.5 ; 14.5 14.5 15.5 15 ; 14 14 12 12.5 ; 11 10 5.5 7 ; 5.5 7 14 11.5 ; 13 12.5 8.5 9.5 ; 9 9.5 12.5 12]
disp('Centrage des donneés de la matrice X')
ind=size(X)
i=mean(X)
for j=1:ind(2)
Y(:,j)=X(:,j)-i(j)
end
disp('La matrice X centrée réduite')
ect=std(X,1)
for j=1:ind(2)
X(:,j)=Y(:,j)/ect(j)
end
disp('La matrice de variance covariance')
C=X'*X
disp('Les valeurs et les vecteurs propres')
vlp=eig(C)
[vct,u]=eig(C)
disp('Classement par ordre decroissant des valeurs propres')
[vpor,posi]=sort(-vlp);
VlpOrd=-vpor
vect=vct(:,posi)
disp('Coordonées des individus sur les axes')
F=[vect(:,1) vect(:,2)];
M=X*F
plot(M(:,1),M(:,2),'.g');
axis([-5 5 -5 5]);
title('Projection des individus sur les axes principaux','edgecolor','red','BackgroundColor','white','fontsize',15);
xlabel('Axe 1 = vect1','edgecolor','red','BackgroundColor','white','fontsize',14);
ylabel('Axe 2 = vect2','edgecolor','red','BackgroundColor','white','fontsize',14); |
Partager