function basic_plot(X,nam0) cud=cd; numObjects=size(X,1); for i=1:numObjects % Par rapport au bleu X1(i,1)=X(i,9)-X(i,13); X1(i,2)=X(i,10)-X(i,14); X1(i,3)=X(1,11)-X(i,15); X1(i,4)=(X1(i,1)^2+X1(i,2)^2+X1(i,3)^2)^(0.5); X1(i,5)=X(i,5)-X(i,13); X1(i,6)=X(i,6)-X(i,14); X1(i,7)=X(i,7)-X(i,15); X1(i,8)=(X1(i,5)^2+X1(i,6)^2+X1(i,7)^2)^(0.5); X1(i,17)=real(acos([X1(i,1) X1(i,2) X1(i,3)]*transpose([X1(i,5) X1(i,6) X1(i,7)])/(X1(i,8)*X1(i,4)))*180/pi); % par rapport au vert X1(i,1)=X(i,5)-X(i,9); X1(i,2)=X(i,6)-X(i,10); X1(i,3)=X(1,7)-X(i,11); X1(i,4)=(X1(i,1)^2+X1(i,2)^2+X1(i,3)^2)^(0.5); X1(i,5)=X(i,13)-X(i,9); X1(i,6)=X(i,14)-X(i,10); X1(i,7)=X(i,15)-X(i,11); X1(i,8)=(X1(i,5)^2+X1(i,6)^2+X1(i,7)^2)^(0.5); X1(i,18)=real(acos([X1(i,1) X1(i,2) X1(i,3)]*transpose([X1(i,5) X1(i,6) X1(i,7)])/(X1(i,8)*X1(i,4)))*180/pi); % par rapport au rouge X1(i,1)=X(i,9)-X(i,5); X1(i,2)=X(i,10)-X(i,6); X1(i,3)=X(1,11)-X(i,7); X1(i,4)=(X1(i,1)^2+X1(i,2)^2+X1(i,3)^2)^(0.5); X1(i,5)=X(i,13)-X(i,5); X1(i,6)=X(i,14)-X(i,6); X1(i,7)=X(i,15)-X(i,7); X1(i,8)=(X1(i,5)^2+X1(i,6)^2+X1(i,7)^2)^(0.5); X1(i,19)=real(acos([X1(i,1) X1(i,2) X1(i,3)]*transpose([X1(i,5) X1(i,6) X1(i,7)])/(X1(i,8)*X1(i,4)))*180/pi); end AA(:,1)=X1(:,19);AA(:,2)=X1(:,18);AA(:,3)=X1(:,17); y1=0:180/10:180; h=figure; hist([X1(:,17) X1(:,18) X1(:,19)],y1);legend('AngB','AngG', 'AngR',1); grid on cd(nam0) saveas(h,'histogramme_angle','tiff'); close(h); cd(cud) h=figure; y=0:2000/10:2000; hist([X(:,3) X(:,2) X(:,4)],y);legend('disBG','disRG', 'disBR',1); grid on; cd(nam0) saveas(h,'histogramme_dist','tiff');close(h); cd(cud) X01=X(:,3);X2=X(:,2);X3=X(:,4); t1=size(X01,1);t2=size(X2,1);t3=size(X3,1); for i=1:t2 FF(i)=X01(i); GG(i,1:6)=char('distBG'); end for i=1:t1 FF(t2+i)=X2(i); GG(t2+i,1:6)=char('distRG'); end for i=1:t3 FF(t1+t2+i)=X3(i); GG(t1+t2+i,1:6)=char('distBR'); end h=figure('Name','Boxplot without color'); boxplot(FF,GG); cd(nam0) saveas(h,'Boxplot','tiff');close(h); cd(cud) h = msgbox('boxplots and histograms saved','Loci analysed','warn');