Vérification de programation juste et interpretation graphique
bonjour a tous
je dois resoudre sous matlab 3 equations differentielles . Prenons la premeire l'equation est :x’(t)=-tx(t)+cos(at) , une equation du premier ordre non lineaire .
Sous matlab j'ai programé:(en utilisant la methode d 'euler)
Code:
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25
| clear all
n=100
a=0;%borne inferieure%
b=5 ;%borne superieure%
a1=0;
b1=-3;
t(1)=0.1; %condition initiale%
x(1)=0.1;%condition initiale%
u(1)=0.1; %condition initiale%
v(1)=0.1;
h=(b-a)/n%pas%
h1=(b1-a1)/n
for i=2:n
t(i)=t(i-1)+h%
x(i)=x(i-1)+h*(cos(t(i-1))-t(i-1)*x(i-1))%dérivée de x par rapport par à t%
u(i)=u(i-1)+h1%incrementation de t par le pas%
v(i)=v(i-1)+h1*(cos(u(i-1))- u(i-1)*v(i-1))%dérivée de x par rapport par à t%
end
plot (t,x,u,v) |
Pouvez vous me dire deja si cette resolution est juste et m'expliquer l'interpretation graphique de resultat(point singuliers,points max, chalmps des tangeantes)
Merci , j'attens vos suggestions.