Bonjour,
je cherche a résoudre l'oscillateur de van der pol, defini par
x''+x-x'*mhu*(1-x^2) avec x'=dx/dt
que l'ont peut ramener a deux équations,
x'=y
y'=-x+mhu(1-x^2)*y
J'ai beaucoup de mal a comprendre l'algorithme de runge-kutta,
est-il correct d'écrire :
x(i+1)=x(i)+0.5(k1+k2)*h
avec k1=y(i) et k2=y(i)+h
y(i+1)=y(i)+0.5(k1+k2)*h
avec k1=-x(i)+mhu(1-x(i)^2)*Y(i) et k2=-(x(i)+h)+mhu(1-(x(i)+h)^2)*(Y(i)+k1*h)
J'utilise comme définition k1=f(xi,yi) et k2=f(xi+h,yi+k1*h) mais je ne suis pas sur d'appliquer l'algorithme correctement.
merci !
Partager