Discussion :

[copc1]

Bonjour,

je souhaite résoudre un système d'équations différentielles du type :
(1) f(t)*f''(t) = A - B g(t) - C f(t)
(2) g'(t) = D f'(t) - E g(t)

je sais résoudre des systèmes d'équations différentiels linéaires sous Matlab avec ode,
mais ici le fait d'avoir le terme f*f'' sur l'équation (1) me bloque et ne permet pas la résolution...
des idées pour linéariser le système ou résoudre autrement ?

Merci

[peaknoise]

(1) f(t)*f''(t) = A - B g(t) - C f(t)
(2) g'(t) = D f'(t) - E g(t)

Bonsoir,
Je suis un peu rouillé en equa diff mais si je ne me trompe pas tu peut sortir les dérivé de cette façon:
f'(t) = (g'(t) + E*g(t))/D
f''(t) = (A - B*g(t) - C*f(t))/f(t)
g'(t) = = D*f'(t) - E*g(t)

Il y a juste le problème avec la division par f(t)/f(t). Il faut éviter que f(t) = 0 pour la division 0/0 qui est indéterminé ou tu simplifie la division f(t)/f(t) = 1.

Il te faut ensuite 2 composantes pour f et 1 pour g pour le vecteur d'état. La suite tu doit connaître.

y = f
y = f'
y = g

y' = f'
y' = f''
y' = g'

[Beltharion]

Sinon tu peux aussi utiliser la fonction fsolve.