la géométrie dans matlab

je débute en matlab, je me suis bloqué :cry:
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:

---

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 :cry: ya t'il quelqu'un qui peut m'aider :cry::cry::cry:

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

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:

---

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')

---

Salut!

Citation:

---

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

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:

---

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+

Citation:

---

Envoyé par Kangourou

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..

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

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..

Merci souviron mais tout ce que vous m'avez expliquer c'est déja fait

Code:

---

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

---

mais mon problème c'est comment teminer le code et tracer le segment OE
je suis franchement très nul en géomètrie et j'éspère que vous pouvez m'aider :cry::cry::cry:

le point est sur un cercle de centre O e de rayon D(O-E), à l'angle calculé..

Franchement, si on sait pas faire ce genre de calculs, c'est même idiot de se servir de MatLab..

Salut!

Citation:

---

comment teminer le code et tracer le segment OE

---

Tu n'as pas l'air de te rendre compte que ton problème comporte deux aspects totalement distincts:

1. La détermination des coordonnées des deux extrémités de ton segment, ce qui est un problème de calcul, dont Souviron34 t'a donné la solution; un tel problème peut se résoudre à l'aide d'un logiciel de calcul ou en faisant les calculs à la main ou avec une calculette.
2. Le dessin d'un segment de droite, connaissant les coordonnées de ses extrémités, ce qui est un problème graphique; un tel problème peut se résoudre à l'aide d'un logiciel graphique ou sur une feuille de papier millimétré avec un crayon et une règle graduée.

Alors, quel est ton problème?
Jean-Marc Blanc

Désolé de te décevoir mais un copier/coller de mon code fonctionne...

La variable "A" est l'angle permettant d'orienter le cercle afin d'obtenir la représentation de toutes les solutions.

Aide-toi de la solution de souviron34 et tu pourras faire une représentation graphique.

Christophe