Discussion :

[sososo3131]

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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.

[mr_samurai]

Salut,

tu es obliger de coder un solveur ?

Car sinon, la fonction DSOLVE permet de résoudre les systèmes différentiels.

++

[sososo3131]

bonjour
Mais qu'est ce qu un solveur?
Il me sembler que pour resoudre les equation diff il fallait utilser les methodes du style euler ou kunge kutta?

[mr_samurai]

Un solveur est une fonction qui résout une équation ou une système d'équations.

[sososo3131]

Penssez vous que la programation de mon equation differentielle est juste?

[mr_samurai]

Salut,

tu confond deux pas qui n'ont rien à avoir. Le Pas de coordonnées en x, et le pas de deplacement de la solution.

En principe, si tu ne vectorise pas ton calcul, tu devrait avoir deux boucles (une sur les points x(i), et la deuxième met a jour t(x(i)) ) .

++ bon courage

[sososo3131]

je n'ai pas bien compris votre suggestion vous penssez que a la place de mettre x(i)=x(i-1)+h*(cos(t(i-1))-t(i-1)*x(i-1))
je devrais mettre
x(i)=x(i-1)+h*(cos(x(i-1))-t(i-1)*x(i-1)) et pareille pour la deuxieme ou ce n'est pas cela dont vous me parlez?