[matlab/octave] TPE équations différentielles

[invite44ec1809]

Bonjour nous fesons un TPE sur la trajectoire d'une balle de tennis. Nous avons grâce à un document trouvé un système d'équation différencielle de premier ordre permettant de créer une courbe de la trajectoire de cette balle.
Nous sommes en 1ère S et nous avons aucune idée de comment résoudre ce système. Nous avons essayé de télécharger plusieurs logiciels permettant de résoudre les équations jointes.
Nous aimerions savoir avec quel logiciel les résoudres, où le télécharger? Puis comment faire pour que ce logiciel résolve les équations
Merci

Ps: Kl est la force de magnus et Kd est la force de frottement.
Ils nous restent quelques jours svp de vite répondre merci .
-----

[invite44ec1809]

Nous avons trouver ça, mais ce code est normalement pour Mtlab marche t-il pour octave?

%programme principal représentant la trajectoire de la balle et le terrain
clear;
theta=0.7;
global vo;
vo=50;
e=0.7;
Dx=23.78;Dy=10.97;Dn=1.067;Ds=5.49;%% dimension du terrain de tennis
figure ('NumberTitle','off','Name','projet tennis')
axis([-1 25 -2 12 -2 16]);
hold on;
xlabel('position de la balle');
zlabel('hauteur de la balle');
axis equal;
options = odeset('event',@rebond);
[t,x] = ode45(@tennis_fct,[0:0.001:10], [0 vo*cos(theta) 5 0 1.5 vo*sin(theta)],
options);
plot3(x(:,1),x(:,3),x(:,5))
% lignes
hold on
plot3(Dx*[0 1 1 0 0],Dy*[0 0 1 1 0],Dn*[0 0 0 0 0],'r'); % court
plot3(0.5*Dx*[1 1 1 1 1],Dy*[0 0 1 1 0],Dn*[0 1 1 0 0],'r'); % filet
plot3(Ds*[1 1],Dy*[0 1],Dn*[0 0],'r'); % service ligne 1
plot3((Dx-Ds)*[1 1],Dy*[0 1],Dn*[0 0],'r'); % service ligne 2
plot3([Ds Dx-Ds],0.5*Dy*[1 1],Dn*[0 0],'r'); % lignes intermédiaires
x;
y=rot90(rot90(x));
y=abs(y);
y(1,:)
z=[y(1,6) y(1,5) y(1,4) y(1,3) y(1,2) y(1,1)]
options = odeset('event',@rebond);
[t,x] = ode45(@tennis_fct,[0:0.001:10], [z(1,1) z(1,2) z(1,3) z(1,4) z(1,5) e*z
(1,6)],options);
plot3(x(:,1),x(:,3),x(:,5))
hold on;
hold off
%Fonction 1 : Résolution de l'équation différentielle du mouvement
function xdot= tennis_fct (t,x)
global vo;
xdot=zeros(6,1);
km=0.01352;kd=0.0011;g=9.81; % kd et km sont les coefficients de trainée et de portance
w=350;% Rotation propre de la balle de tennis -> Effet Magnus
xdot(1)=x(2);
xdot(3)=x(4);
xdot(5)=x(6);
xdot(2)= -kd*vo*x(2)+km*w*x(6);
xdot(4)= 0;
xdot(6)=-g-kd*vo*x(6)-km*w*x(2);
%% Fonction 2 : Permet de modéliser le rebond
function [value,isterminal,isdirection] = rebond(t,x)
value = x(5);
isterminal = 1;
isdirection = 0;

On a essayer ça marche jusque là:
options = odeset('event',@rebond);

[Flyingsquirrel]

Salut,

J'ai déplacé la discussion dans le forum « mathématiques du supérieur » car il est plus adapté à ta question.

Quand tu veux joindre du code à l'un de tes messages, il vaut mieux utiliser la balise code (icône ), ça rend ton message plus lisible.

Pour la modération, Flyingsquirrel.

[invite408a730a]

On a découvert que il y a des différences entre Matlan et Octave mais on ne trouve pas de guide les expliquant toute j'ai vu par exemple que odeset et en fait odeget sur Octave..

[A voir en vidéo sur Futura]

[invite408a730a]

Up!
Personne n'a une idée ou un bon lien? On bloque sur

@rebond

qui ne marche pas sur Octave