Bonjour,
je cherche a résoudre numériquement l'équation de van der Pol en utilisant des dérivées décentrées, mais je ne sais pas du tout comment choisir le pas.
Forcement je peux le prendre arbitrairement très petit et être sur que ma méthode converge, mais je voudrai savoir si il y a une autre méthode que d'expérimenter a l'aveuglette.
//
l'oscillateur de VDP est défini 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
La méthode des dérivées décentrées étant
x'=(x(i+1)-x(i)) * dt, on obtient :
x(i+1)=y(i)*dt-x(i)
Y(i+1)= (-X(i)+mhu*(1-(X(i)*X(i)))*Y(i))*dt+Y(i)
mais comment fait on pour choisir dt ? En experimentant, ou je peux le calculer ?
Merci bien
-----