Discussion :

[dahman59155]

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

[phryte]

Bonjour,

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

[dahman59155]

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 :

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

[phryte]

Bonjour,

Il y a plusieurs solutions.
Un exemple :

Code :

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

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)

[dahman59155]

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

[phryte]

Bonjour,
Autre méthode :

Code :

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

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)