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
| close all; clear all; clc
R = 100/255; G = 100/205; B = 100/205; RGB = [R G B];
phi = -2*pi:0.01:2*pi;
M = xlsread('datas.xlsx');
a = M(1:8:end,1);
b = M(1:8:end,2);
c = M(1:8:end,3);
Norm = sqrt(a.^2 + b.^2 + c.^2); %Pour placer les points à la surface
a = a./Norm;
b = b./Norm;
c = c./Norm;
hF = figure('Color','w','Position',[150 250 350 300]);
set(gca,'FontSize',18,'Position',[0.17 0.17 0.8 0.9])
[X,Y,Z] = sphere(100);
hL = surf(X,Y,Z,'FaceAlpha',0.25,'FaceColor',RGB,'LineStyle','none',...
'FaceLighting','phong','AmbientStrength',0.5,...
'DiffuseStrength',0.7,'BackFaceLighting','reverselit');
axis([-1.1 1.1 -1.1 1.1 -1.1 1.1]); axis vis3d; grid off; hold on
light('Position',[1 1 1]); axis on
plot3(a,b,c,'or','LineStyle','none','MarkerEdgeColor','r','MarkerFaceColor','r','MarkerSize',4)
Xlabel('S_{1}','FontSize',18,'Position',[1.4 -0.1 -2.5])
Ylabel('S_{2}','FontSize',18,'Position',[-0.40 1.8 -2.5])
Zlabel('S_{3}','FontSize',18,'Position',[-1.6 1.5 -0.25],'rot',0)
plot3(sin(phi),cos(phi),0*phi,'Color','w')
plot3(0*phi,cos(phi),sin(phi),'Color','w')
plot3(sin(phi),0*phi,cos(phi),'Color','w')
plot3([-1 1],[0 0],[0 0],'Color','w')
plot3([0 0],[-1 1],[0 0],'Color','w')
plot3([0 0],[0 0],[-1 1],'Color','w') |
Partager