Exercice sur fsolve de scipy

Bonjour,

Je suis en train de faire un exercice sur scipy et je ne comprends pas l'erreur que j'ai faite:

Voici l'exercice:
Pièce jointe 579897

Voici mon code:

Code:

---

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 26 27 28 29 30

from scipy import linalg   ##Question 1: def f(x): M = np.array([[1, 2*x - 3, 3] , [-2, -1, x] , [3, 2, 1]])   return linalg.det(M) #La fonction retourne le déterminant de M     ##Question 2: y = np.linspace (-3, 5, 1000)   #On crée une liste nommée ff qu'on remplit avec les images de y par la fonction f   ff=[] for i in y: ff.append(f(i))   plt.plot(y,ff) plt.title("Tracé de f dans l'intervalle [-3,5]") plt.xlabel("x") plt.ylabel("f(x)")     ## Question 3: #Les zéros: #Imporation du module fsolve: from scipy.optimize import fsolve   zeros = fsolve(f,[-3,5])

---

et là il m'affiche :

Citation:

---

ValueError: object arrays are not supported

---

Merci d'avance et bonne journée :)

Salut,

fsolve() envoie à f() un array de la même taille que celle de l'argument x0. Or la fonction f n'accepte pas ce type car numpy n'est pas en mesure de calculer le déterminant d'un truc du genre :

Code:

---

1 2 3 4

mat = np.array([ [1, 2*[0, 1], 3], [0, 0, 1], ])

---

Il y a plusieurs façon de résoudre le problème :

1. Modifier la fonction f()
2. Créer une fonction intermédiaire
3. Utiliser une fonction jetable via l'emploi de lambda

1-
Pour moi pas la bonne approche, mais ce serait quelque chose du genre :

Code:

---

1 2 3 4

def f(x): _x = x[0] mat = np.array([[1, 2*_x - 3, 3] , [-2, -1, _x] , [3, 2, 1]]) return linalg.det(mat)

---

2-

Code:

---

1 2 3 4 5

def obj_func(x): """Objective function : fonction a pour laquelle on cherche les racines""" return f(x[0])   zeros = fsolve(obj_func, 0.0)

---

3-

Code:

---

1 2	zeros = fsolve(lambda x: f(x[0]), 0.0)

---

Dans tous les cas, il définir un nombre comme valeur de départ (argument x0 de fsolve), pas les bornes. D'ailleurs cette fonction n'accepte pas de bornes. Il faut savoir à l'avance à peu près où se trouvera la racine.

J

Aaahh ouiiii je comprends très bien!!

Je vous remercie énormément pour votre réponse détaillée et votre aide :)