Discussion :

[slaima15]

je débute en matlab, je me suis bloqué
je vous explique tout d'abords mon besoin puis mon problématique: j'ai deux points A et B et je veux tracé la perpendiculaire OE a ces deux points qui passe par leur milieu le point O, le segment OE doit avoir un longeur de 1/3 de AB.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
 
xA=100
yA=300
xB=200
yB=500
DAB=sqrt(((xA-xB)^2)+((yA-yB)^2))
DAB=DAB/3
xO1=(xA+xB)/2
yO1=(yA+yB)/2

je me suis bloquer maintement comment tracer le segment OE qui doit être perpendiculaire a AB au point O et de longueur AB/3 ya t'il quelqu'un qui peut m'aider

[phryte]

Bonjour.
Tu peux t'inspirer de la discussion là :
http://www.developpez.net/forums/d80...trice-segment/

[christophe_halgand]

Ce n'est pas un problème matlab mais de géométrie...

Pour ton segment OE, il y a une infinité de solution puisque ta seule contrainte est la "perpendicularité" au segment AB. (ton point E se balade sur le cercle centré en O, exécutes le code ci-dessous pour visualiser).

Si tu rajoutes la contrainte d'être dans le plan où se situe les points A et B, tu n'as plus que deux solutions (si tu visualises le segment A vers B, le point E peut se trouver à droite comme à gauche du segment AB)...

Par rapport à ta manière de présenter ton problème, la solution est :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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
xA=100;
yA=300;
xB=200;
yB=500;
DAB=sqrt(((xA-xB)^2)+((yA-yB)^2));
DAB=DAB/3;
xO1=(xA+xB)/2;
yO1=(yA+yB)/2;
 
A = - acos((yB-yO1)/(DAB*3/2));
 
teta = -pi:pi/180:pi;
xE = xO1 + DAB*cos(teta);
yE = yO1*ones(length(teta),1);
zE =  0  + DAB*sin(teta);
 
xE0 = DAB*cos(teta);
yE0 = zeros(size(xE0,1),size(xE0,2));
zE0 = DAB*sin(teta);
 
Rot_Z = [cos(A)   -sin(A)   0   0; 
         sin(A)    cos(A)   0   0;
             0         0    1   0;
             0         0    0   1];
 
for q=1:length(xE0)
    a(q,:)=Rot_Z*[xE0(q);yE0(q);zE0(q);1];
end
 
for w=1:length(a)
    xE1(w)=a(w,1)+xO1;
    yE1(w)=a(w,2)+yO1;
    zE1(w)=a(w,3);
end
 
figure
hold on
line([xA xB],[yA yB]);
plot3(xE1,yE1,zE1,'Color','r')

[Kangourou]

salut,

j'ai développé geom2d, une bibliothèque de calculs géométrique, qui te permettra de simplifier certains calculs.
http://www.mathworks.fr/matlabcentral/fileexchange/7844

Par exemple, une fois la lib téléchargée et installée, pour tracer un segment perpendiculaire à une autre et commençant a son milieu tu peux faire :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
 
a = [2 1];
b = [8 5];
edge = createEdge(a, b);
cla
drawPoint([a;b], 'ro');
hold on;
axis([0 10 0 10]);
axis equal;
drawEdge(edge);
line = edgeToLine(edge);
o = centroid([a;b]);
perp = orthogonalLine(line, o);
dist = edgeLength(edge);
e = pointOnLine(perp, dist);
oe = createEdge(o, e);
drawEdge(oe, 'color', 'k')

Au cas ou ça aide...
A+

[souviron34]

salut,

j'ai développé geom2d, une bibliothèque de calculs géométrique, qui te permettra de simplifier certains calculs.

merci pour ton message et lien, mais pas merci par rapport à la question du PO.

Dans un cas de géométrie simple, mieux vaut apprendre à s'y retrouver en réfléchissant et avec les équations (le problème est simple) que prendre un outil qui "fait tout"..

Ce n'est pas lui rendre service..

[slaima15]

Merci christophe ,mais ton code n'as pas marcher et elle me donne cette erreur "Subscripted assignment dimension mismatch" sur la ligne "a(q,=Rot_Z*[xE0(q);yE0(q);zE0(q);1];"
concernant les contraintes, j'ai deux contraintes qui sont :
le segment OE qui doit être perpendiculaire à AB en O et la longeur de OE doit être le 1/3 de AB

[souviron34]

bon, prenons logiquement, et mathématiquement :


A (Xa,Ya)
B (Xb,Yb)

O = milieu de (A-B)

Xo = (Xa+Xb)/2
Yo = (Ya+Yb)/2


Angle formé par le segment AB et l'horizontale : sa tangente est la pente . Donc l'angle est l'arc tangente.

Donc l'angle perpendiculaire est cet angle + PI/2


d(A_B) = racine ( (Xa-Xb)^2 + (Ya-Yb)^2)

d(O-E) = 1/3 d(A-B)


Et là tu finis tout seul..

[FR119492]

Salut!

j'ai deux points A et B et je veux tracé la perpendiculaire OE a ces deux points qui passe par leur milieu le point O, le segment OE doit avoir un longeur de 1/3 de AB.

Un marteau est très utile pour planter des clous; une brouette est tout aussi utile pour transporter du gravier. Mais si tu essaies de planter des clous avec une brouette ou de transporter du gravier avec un marteau, tu verras que ça ne va pas très bien. De manière analogue, Matlab est un excellent outil pour effectuer des calculs, mais ce n'est pas un outil de dessin.

D'autre part, comme te le fait remarquer Christophe Halgand, c'est un problème de géométrie. En conséquence, je déplace cette discussion dans le forum algo/maths.
Jean-Marc Blanc