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 46 47 48 49 50 51 52 53 54 55 56 57 58 59
| x1 = A(:,1);
y1 = A(:,2);
x2 = A(:,3);
y2 = A(:,4);
% Réduction des points
x11 = A(1:47:end,1);
y11 = A(1:47:end,2);
x12 = A(1:47:end,3);
y12 = A(1:47:end,4);
% Maillage
dt = DelaunayTri(x11,y11);
numtri = size(dt,1);
MVT=zeros(numtri,7);
for i = 1:numtri
% Indexation linéaire des points des triangles
nT =dt(i,:); % ième triangle
np1=dt(i,1);
np2=dt(i,2);
np3=dt(i,3);
% Coordonnées des 3 points de chacun des triangles
x_p1=x11(np1);
y_p1=y11(np1);
x_p2=x11(np2);
y_p2=y11(np2);
x_p3=x11(np3);
y_p3=y11(np3);
p1=[x_p1,y_p1];
p2=[x_p2,y_p2];
p3=[x_p3,y_p3];
% Centre de gravité de chacun des triangles du maillage
x_CG=(x_p1+x_p2+x_p3)/3;
y_CG=(y_p1+y_p2+y_p3)/3;
% Points appartenant à chacun des triangles
tol=1e-3;
pp=[x1,y1];
node = [x11(nT) y11(nT)];
[in, on] = inpoly(pp, node, [], tol);
%Obtenir les indices des points appartenant à chacun des triangles
indices1 = find(in);
end |
Partager