Bonjour
je débute le cours des différences finies pour résoudre des équations différentielles mais je suis un peu bloquée et j ai besoin d un coup de pousse ..
prenant cet exemple qui va m aider pour bien comprendre la stratègie de la méthode et son utilité:
on cherche à déterminer la solution exacte et la solution approchée de cette équations par la méthode des différences finies :et par la suite on compare les résultats:
-U" +4x *U=2 e^(-x²)
U(0)=1
U(1)=1/e
PS : ils ont trouvé comme solution exacte U(x)=e^(-x²) mais je ne sais pas comment ils ont trouvé cette solution !
mes réponses:
concernant la méthode des différences finies
on remplace U" par -Uj+1-2Uj+Uj-1 / h avec h=b-a/N+1 = 1/e(N+1)
et on obtient :
-Uj+1+2Uj-Uj-1 +4xj *h* Uj=2 *h * e^(-xj²)
alors :
-Uj-1 +Uj[2+4xj*h]-Uj+1 =2 *h * e^(-xj²)
en utilisant la notation matricielle on obtient ;
AU=F
ou :
F=[2 *h * e^(-(x1)²) ;2 *h * e^(-(x2)²).....;2 *h * e^(-(xN)²)]
U=[U1 ;U2........;UN]
A=[-1 2+4x1*h -1 0..0;0 -1 2+4*x2*h -1..0 ; ....2+4*xN*h -1]
(F , U et A sont des matrices , les lignes sont séparées par des point virgules et les colonnes par un vide )
si A est inversible alors U=A^(-1)*F
je suis bloquée ici ,comment résoudre ce système S: pour trouver U
merci de m avoir aidée
-----