Comment afficher la valeurs d'un point d'une courbe précise ?

Je bug un peu sur un problème simple je pense ,

J'ai tracé la trajectoire d'une balle à une vitesse V(0) et je voudrais avoir le point d’impact au sol c'est à dire sortir la valeurs de x à y=0 ( et non la première valeur initiale y=0 à x=0 ).

Quelle est la commande s'il vous plait ?

Bien à vous
Da

Bonjour,

Tu as tracé la trajectoire à partir d’une équation ?

Citation:

---

Envoyé par phryte

Bonjour,

Tu as tracé la trajectoire à partir d’une équation ?

---

Exactement ;

Oui j'ai utilisé 2 méthodes Ode45 et par l’équation de la trajectoire théorique.
Je vous voudrais juste afficher le point sur X de l'impact à Y=0

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

clear close all   % Methode ODE45   t0=0; tf=10.5; samp=200;   t=linspace(t0,tf,samp); v0 = input(' Quelle est la Vitesse initiale V0 ? ') y0=[0,v0*cos(pi/6),0,v0*sin(pi/6)]; % Conditions initiales sur y'   [t,y]=ode45('f5test',t,y0);     subplot(2,1,1) plot(y(:,1),y(:,3)); grid title('Trajectoire de la balle'); xlabel('Y1'); ylabel('Y3');     % Methode théorique   g=9.81; % gravité V0=100; % Vitesse initiale en m/s   x=0:1000;   y=-(g)*x.^2/(2*(V0^2)*cos(pi/6)^2)+ tan(pi/6).*x; % Trajectoire subplot(2,1,2) plot(x,y,'r') % représentation graphique grid title('Trajectoire de la balle Théorique'); xlabel('Y1'); ylabel('Y3');

---

Bonjour,

Il y a plusieurs solutions.
Un exemple :

Code:

---

1 2 3 4 5 6

xlabel('Y1'); ylabel('Y3'); n=x(2:end); Ind=find(y(n).*y(n+1)<0); hold on plot(x(Ind),y(Ind),'.b','MarkerSize',15)

---

Merci de votre réponse vous n'auriez pas une autre méthode ( par exemple d’intersection ) pour afficher ce point ?

Bonjour,
Autre méthode :

Code:

---

1 2 3 4

hold('on') m=abs(min(diff(y))) idx=find(abs(y(2:end))<=m/2) plot(x(idx),y(idx),'*b','MarkerSize',13)

---