Discussion :

[bilou_12]

Bonjour,

J'ai une fonction main4 qui appelle la fonction descente1 dans le but de coder la méthode des descentes. Problème: le résultat est NaN...

J'ai l'impression que le code ne rentre pas dans la fonction descente1...

Pourriez-vous m'indiquez pourquoi svp?

Ci-dessous le code de main4:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
clear all % efface tout en mémoire
close all %ferme fenêtres graphiques
clc %efface console
format long
 
tstart = tic; %start du temps
 
x0=[2.2;2.2]; 
eps=1e-6;
[X,nb,XX]=descente1(x0,eps)
plot(XX(1,:),XX(2,:), 'o-')
 
telapsed = toc(tstart) %stop du temps

et celui de descente1:

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
function [X,nb,XX]=descente1(X0,eps)
%méthode de Newton
X=X0;
XT=X0; 
XX=[X];
nb=0; 
err=1;
maxiter = 200
%choix du pas
alpha=0.1
 
while (err>eps && nb<maxiter)
        %direction de descente -gradient
        d=-df(X);
        %calcul de x(k+1)
        X=X+alpha*d;
        %calcul de l'erreur
        err=norm(X-XT);
        %iteration
        XT=X;
        %compte nb iteration
        nb=nb+1;
        %complète la matrice
        XX=[XX,X];
end

Merci d'avance!

[Jean Dumoncel]

Bonjour,

si tu obtiens un résultat, même nan, c'est que le code rentre dans la fonction.

Que fait la fonction df?

[tachmou]

Bonjour,

Je plussoie le commentaire de Magelan : peux-tu nous montrer la fonction df stp ?

Un NaN peut venir d'une division par zéro... qui peut se produire lors du calcul d'une dérivée mal codée par exemple...

Cordialement,

[Jean Dumoncel]

Un NaN peut venir d'une division par zéro...

Il faut aussi que le numérateur soit nul sinon, matlab renvoie la limite (Inf).

These operations produce NaN:

• Any arithmetic operation on a NaN, such as sqrt(NaN)
• Addition or subtraction, such as magnitude subtraction of infinities as (+Inf)+(-Inf)
• Multiplication, such as 0*Inf
• Division, such as 0/0 and Inf/Inf
• Remainder, such as rem(x,y) where y is zero or x is infinity

[bilou_12]

Salut tout le monde!

Merci pour votre aide!

Voici ci-dessous df, dérivée de ma fonction:

Code :

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

1
2
3
4
5
function [a] = df(X)
a = zeros(2,1);
x=X(1,1);
y=X(2,1);
a = [-400*x*(y-x.^2)-2*(1-x); 200*(y-x.^2)];

et d2f, dérivée seconde:

Code :

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

1
2
3
4
5
function [a]=d2f(X)
a=zeros(2,2);
x=X(1,1);
y=X(2,1);
a=[1200*x.^2-400*y+2, -400*x; -400*x, 200];

[Jean Dumoncel]

En fait les nan sont obtenus car ton algo diverge : réduit ton pas (valeur de alpha).

[bilou_12]

Merci tout le monde!

Ca marche avec un pas de 0.001!

Par contre, si vous avez testé mon programme, que signifie les points sur la courbe? Que cette méthode n'est pas trop adaptée?

[Jean Dumoncel]

Par contre, si vous avez testé mon programme, que signifie les points sur la courbe? Que cette méthode n'est pas trop adaptée?

Ce sont les valeurs trouvées à chaque itération. Je n'ai pas vu de problèmes particulier, pourquoi dis-tu qu'elle est mal adaptée?