bonjour,
je cherche à résoudre l'équation différentielle suivante: (il s'agit de l'équation régissant le modèle de Dahl)
R'(t ) = b ×u'(t )×[h(u) - R×sgn(u'(t ))]
avec
h (u) = 1/2 [(hu + hl) *sgn(u'(t)) + (hu - hl)]
hu=a1*u+b1
hl=a2*u+b2
R:force
u: déplacement de forme sinusoïdale (frequence f connue)
t:temps
les constantes (b, a1, a2,b1, b2) sont déterminées.
le morceau de pgrm :
function [y]=sol(T,Te,beta,a1,a2,b1,b2,R0, f);
T0=[0:Te:T];
options=odeset('RelTol',1e-5);
[t,y]=ode45(@sec,T0,R0,options,beta ,a1,a2,b1,b2,f);
le tracé de R en fonction de u est nettement divergent pour la plupart des fréquences f
que faire pour y remédier?
merci d'avance
-----