Discussion :

[langv26]

Bonjour à tous

Dans le cadre d'un exercice d'optimisation statistique, nous sommes amenés à exécuter la méthode de newton Raphson sur la fonction de Rosenbrock.
J'ai donc tapé ce code là :

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
 
import numpy as np
from numpy import exp,log
from scipy.special import factorial
from scipy import optimize
 
## Fonction automatique de Newton-Raphson
 
def f(x):
  x = x[1:]
  xp = x[:-1]
  return 200 * (x - xp ** 2)
 
x = np.array([0.4, 0.6])
optimize.newton(f, x)

Pourtant à l'exécution, j'obtiens ce message d'erreur :

Code :

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

1
2
 
ValueError: operands could not be broadcast together with shapes (0,) (2,)

J'ai beau me creuser la tête, mais je ne vois absolument pas comment résoudre ça.

Si vous pouviez me donner un coup de pouce, je vous en serait vraiment reconnaissant

Merci !

[wiztricks]

J'ai beau me creuser la tête, mais je ne vois absolument pas comment résoudre ça.

avant de résoudre un problème, il faut essayer de le comprendre.
Ajoutez des print dans f:

Code :

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

1
2
3
4
5
6
7
8
def f(x):
   print(x)
   x = x[1:]
   xp = x[:-1]
   print(x, xp)
   a = 200 * (x - xp ** 2)
   print (a)
   return a

Et posez vous la question de savoir si ce que çà affiche est ce que vous attendiez.

- W