Discussion :

[Fleur du Pays]

Tout d'abord bonjour à tous, j'aimerais avoir votre aide concernant la programmation de la méthode d'Euler Implicite, je vous montre ce que j'ai effectué.

EDIT: j'ai supprimé mes programmes suite à la copie par certains élèves de ma promo (ils se reconnaitrons :p), demandez donc moi plutôt que de copier des programmes à moitié faux que vous ne comprenez pas ^^

[phryte]

Bonjour.

function y=dF(x);
y=[0 1;-1 0];

Tu n'utilises pas l'argument d'entrée x ?

[Fleur du Pays]

dF correspond à la Jacobienne de F, (sauf erreur de ma part sa donne celà ^^) donc je n'ai pas utilisé x pour cela, néanmoins je m'étais dit qu'en y mettant sous cette forme se serait plus simple si je devais le refaire avec une autre fonction.
Le fait que je ne l'utilise pas pourrait donc créer un soucis? =/

EDIT:

J'ai modifié une ligne de ma fonction de Newton par celle ci

Code :

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

Nw=Nw-inv(feval(dF,Nw))*feval(F,Nw);

Cela m'enlève déjà la seconde erreur, mais la première persiste =/...

[Jerome Briot]

Je ne vois pas très bien l'intérêt de FEVAL ici

A mon avis, tu devrais reprendre ton algorithme depuis le début

[Fleur du Pays]

Je me suis inspiré justement de ce que j'ai pu trouver afin d'essayer de construire cette petite chose bancale qui ne marche pas ^^'... Comme je l'ai souligné, je suis plutôt débutant en la matière. Il me semble avoir réussi pour Euler Explicite, mais même en reprenant tout du début pour Euler Implicite, le problème est que je ne vois pas par où commencer différemment...

Concernant feval, je dois bien avouer ne pas trop connaitre les moments précis où sont utilisation est strictement nécessaire, mais les fonctions f et df sont des arguments d'entrée de la fonction newton et "il est donc
nécessaire d'utiliser la commande feval pour les appeler" à ce que j'ai pu recevoir comme informations jusqu'à présent...

Mis à part cela, il me semble par contre m'être rendu compte que je n'applique pas newton aux bonnes fonction F et dF pour construire un schéma d'Euler Implicite. Il me faudra relire quelque peu mon cours là dessus pour en être certain.

Mais j'ai quelque peu l'impression d'être dans une impasse,ne sachant plus trop que faire pour m'en sortir...

[salseropom]

Salut, dans ta méthode de Newton, tu peux pas plutôt faire plus simplement :

Code :

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

1
2
 
Nw=Nw-F(Nw)\dF(Nw);

Ca m'a l'air plus simple.

Autre probème : il faut tester que ta matrice jacobienne soit inversible (donc 0 ne doit pas être une valeur propre...)
Si c'est le cas, tu peux t'inspirer d'une méthode de Newton modifiée... Si J(t) non inversible alors tu poses J(t) = J(t - 1) ou J(t) est la jacobienne à l'instant t