Discussion :

[ramrouma]

Salut tout le monde

Je voudrais connaître comment résoudre une équation différentielle non linéaire de 2ème ordre sachant que j'ai utilisé ode45 et que ça ne marche pas

Mon équation est la suivante :

Code :

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

md2z/dt+a dz/dt=fm sin(2x/l)-fo signe(dz/dt)

Aidez-moi
Merci

[salseropom]

C'est normal que ça ne marche pas. Les solveurs de matlab ne savent résoudre que des équa diff d'ordre 1. Pour cela : tu augmentes la dimension de ton équation différentielle mais tu diminues son ordre

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
 
EDO à résoudre : 
 
m z'' + a z' = fm sin(2x/l)-fo signe(z')   (1)
 
c'est la même chose que de dire
 
 
tu poses y=z et tu obtiens le système suivant
 
y'=z'
y''=(fm sin(2x/l)-fo signe(y') - a y')/m
 
tu poses X=(y, y')
tu arrives alors au système suivant
 
X'=(z',(fm sin(2x/l)-fo signe(y') - a y')/m)

qui est un système de d'ordre 1 mais de dimension 2 et ça, ode45 c'est faire