Discussion :

[juliettesch]

Bonjour, j'essaye de tracer le graphe de ma fonction g définie ci dessous mais ça me met une erreur que je ne comprends pas... ma fonction marche mais c'est le traçage du graphe qui ne marche pas..

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
def fact(n):
    if n<2:
        return 1
    else:
        return n*fact(n-1)
 
def caga(a,x):
    if x>0 and x<1 :
        return (1/fact(a-1))*x**(a-1)
    elif x>=1 :
        return (1/fact(a-1))*np.exp(-x)
    else :
        return 0
 
caga(2,3)
 
x_grid = np.linspace (0,5,100)
plt.plot(x_grid,caga(a,x_grid))

ça me met :

ValueError: The truth value of an array with more than one element is ambiguous. Use a.any() or a.all()

lorsque je remplace par

Code :

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

1
2
3
4
5
6
7
def caga(a,x):
    if x.any()>0 and x.any()<1 :
        return (1/fact(a-1))*x**(a-1)
    elif x.any()>=1 :
        return (1/fact(a-1))*np.exp(-x)
    else :
        return 0

ça me met l'erreur suivante :

x and y must have same first dimension, but have shapes (100,) and (1,)

Je ne comprends vraiment pas.. Merci d'avance pour votre aide

[Sve@r]

Bonjour

Dans la fonction caga() (qui est appelée inutilement à la ligne 15), tu utilises les arguments "a" et "x" comme des nombres individuels. Et lors de ton appel en ligne 18, tu lui passes x_grid qui n'est pas un nombre mais un tableau de nombres. Fatalement il ne comprend rien mais tente quand-même de t'aider à sa façon.
Et ce n'est pas en remplaçant x par x.any() (fonction renvoyant vrai si un de ses éléments de "a" est vrai) que ça va améliorer les choses (surtout qu'à la ligne 3 de ta fonction, tu as laissé "x" tel quel).

[juliettesch]

Bonjour,
merci pour votre réponse.. Je pensais que le plot allait tracer ma fonction en chaque point de x_grid. Du coup il faut que je fasse une boucle ?
Dans le cours nous avions pourtant fait quelque chose de semblable et cela avait marché

Code :

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

1
2
3
4
5
6
7
def f(x):
    return((2/np.pi)*np.sqrt(1-(x-1)**2))
 
x_grid = np.linspace(0,2,100)
#print(x_grid)
plt.figure(figsize=(8,5))
plt.plot(x_grid,f(x_grid))

[wiztricks]

Salut,

Du coup il faut que je fasse une boucle ?

Vous faites comme vous voulez à partir du moment où vous comprenez ce que vous cherchez à faire et que vous obtenez le résultat attendu.

Dans le cours nous avions pourtant fait quelque chose de semblable et cela avait marché

Si 2 choses qui se ressemblent ne se comportent pas de la même façon, c'est qu'il faut essayer de comprendre ce qui est différent.
En général, le message d'erreur est significatif... et utile pour faire une recherche sur Internet si on n'a pas d'idée.

- W

[Julien N]

Salut,

Lorsque vous faites ceci :

Code :

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

1
2
3
4
5
def f(x):
    return((2 / np.pi) * np.sqrt(1 - (x - 1)**2))
 
x = np.array([1, 2, 3])
f(x)

Vous faites une opération vectorielle. A la place de faire une opération sur chaque élément de votre array x, vous faites l'opération sur tous les éléments en une seule fois. Numpy est ok avec ça, et c'est même l'un de ses grands atouts, mais dès que vous demandez ceci :

Code :

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

1
2
x = np.array([1, 2, 3])
x > 2

Il ne comprend pas ce que vous lui demandez. Et c'est normal. Vous avez deux possibilités :

1. Appliquer votre fonction à chaque terme de votre array (une boucle)
2. transformer votre opération de sorte à ce qu'elle fonctionne sur l'ensemble à la fois

Si vous débutez, partez sur la première.

J

[juliettesch]

Bonsoir,
Merci pour toutes vos réponses. J'essaye de faire une boucle, en effet, je m'approche peu à peu du résultat.