Discussion :

[Youyayouyou]

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:


Voici mon code:

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
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 :

ValueError: object arrays are not supported

Merci d'avance et bonne journée

[Julien N]

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 :

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

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 :

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

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 :

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

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 :

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

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

[Youyayouyou]

Aaahh ouiiii je comprends très bien!!

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