Scribd
Importer
4 vues2 pages

Simulation d'une roue libre dynamique

Le document présente un modèle de simulation d'une roue libre en utilisant des équations différentielles. Il définit les paramètres du modèle, les conditions initiales et résout les équations avec la fonction ode45. Les résultats sont ensuite tracés pour montrer l'évolution de la position et de la vitesse angulaire au cours du temps.

Transféré par

deogracias kawanda
Copyright
© All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats DOCX, PDF, TXT ou lisez en ligne sur Scribd
4 vues2 pages

Simulation d'une roue libre dynamique

Le document présente un modèle de simulation d'une roue libre en utilisant des équations différentielles. Il définit les paramètres du modèle, les conditions initiales et résout les équations avec la fonction ode45. Les résultats sont ensuite tracés pour montrer l'évolution de la position et de la vitesse angulaire au cours du temps.

Transféré par

deogracias kawanda
Copyright
© All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats DOCX, PDF, TXT ou lisez en ligne sur Scribd

% Paramètres du modèle

R = 1; % Rayon de la roue libre

K = 1; % Constante de raideur

b = 0.1; % Coefficient de frottement

% Conditions initiales

theta0 = 0; % Position initiale

omega0 = 0; % Vitesse angulaire initiale

% Intervalle de temps

tspan = [0 10];

% Fonction des équations différentielles

function dXdt = roue_libre(t, X)

% X(1) = theta (position)

% X(2) = omega (vitesse angulaire)

% Les équations du système

dXdt = [X(2); -b * X(2) - K * X(1)];

end

% Résolution des équations différentielles

[t, X] = ode45(@roue_libre, tspan, [theta0 omega0]);

% Tracé des résultats


figure;

plot(t, X(:,1), '-o', 'DisplayName', 'Position');

hold on;

plot(t, X(:,2), '-x', 'DisplayName', 'Vitesse angulaire');

xlabel('Temps (s)');

ylabel('Valeur');

title('Comportement d\'une roue libre');

legend;

grid on;

Vous aimerez peut-être aussi