Maximums à partir d'un tracé de plusieurs courbes

Version imprimable

1 pièce(s) jointe(s)

Bonjour à tous,
j'ai tracé 6 courbes correspondants aux 6 débits de piston d'une pompe hydraulique
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 import numpy as np import matplotlib.pyplot as plt #créer un tableau de 100 valeurs de t compris entre 0 et 10 secondes t=np.linspace(0,10,1000) omega=2 R=20 r=4 alpha=22 n=6 #nb de pistons #Définir une fonction de tracé def f_trace_courbe (x,y,color,xlim,titre,ylabel,xlabel): plt.plot (x,y,color) plt.title(titre) plt.xlim(xlim) plt.ylabel(ylabel) plt.xlabel(xlabel) plt.show() ##Programme #calcul du debit for i in range (1,n+1) : debiti=omega*R*np.tan(alpha)*np.sin(omega*t+(i*2*np.pi/n))*np.pi*r**2 f_trace_courbe (t,np.abs(debiti),'b',(0,10),'Débit de la pompe','debiti','temps')
ce qui me donne le tracé suivant :
Pièce jointe 503622

Maintenant, j'aimerais tracer la courbe qui passe par le maximum de chacune des 6 courbes ci-dessus.
Lorsque je récupére np.max, python me retourne bien entendu une seule valeur ; j'ai pensé récupérer les points d'intersection des courbes mais cela ne colle pas avec mon codage en i.

Auriez vous une idée de code ? Merci d'avance pour votre aide.

17/09/2019, 09h32
wiztricks

Salut,

Citation:

Envoyé par segolsed

Auriez vous une idée de code ? Merci d'avance pour votre aide.

Une courbe c'est une liste de coordonnées sur l'axe des X associée à une liste de coordonnées sur l'axe de Y. Et cette association passe par les indices i.e. X[i] est associé à Y[i] ou dit plus mathématiquement Y[i] = F(X[i]).

A partir de là, pour trouver le x associé au max des Y, il faut récupérer (aussi) l'index de ce max (le i).

- W
17/09/2019, 09h37
segolsed

et ici le i est le temps (l'axe des abscisses) mais la difficulté vient du fait que j'utilise une seule fonction de tracé (voir le code) et non pas 6.
17/09/2019, 10h34
wiztricks

Citation:

Envoyé par segolsed

et ici le i est le temps (l'axe des abscisses) mais la difficulté vient du fait que j'utilise une seule fonction de tracé (voir le code) et non pas 6.

Je veux bien que vous ayez des difficultés à stocker les x[i], y[i] correspondant aux maxima des différentes courbes dans une liste... mais je ne vois pas trop lesquelles.

- W

1 pièce(s) jointe(s)

salut,

pour commencer le code que tu donnes ne correspond pas à l'image qui l'accompagne, il faut au minimum sortir le plt.show() de la fonction et l’appeler tout à la fin après la boucle pour obtenir le même rendu.
Code:

1 2 3 4 5 6 7 8 9 10 11 12 def f_trace_courbe (x, y , color, xlim, titre, ylabel, xlabel): plt.plot(x, y, color) plt.title(titre) plt.xlim(xlim) plt.ylabel(ylabel) plt.xlabel(xlabel) for i in range (1,n+1): debiti = omega * R * np.tan(alpha) * np.sin(omega * t + (i * 2 * np.pi / n)) * np.pi * r**2 f_trace_courbe (t, np.abs(debiti), 'b', (0,10), 'Débit de la pompe', 'debiti', 'temps') plt.show()
Citation:

Envoyé par segolsed

Maintenant, j'aimerais tracer la courbe qui passe par le maximum de chacune des 6 courbes ci-dessus.

une courbe qui passe par 1 seul point ? tout ça remonte à loin mais je crois bien qu'il y en a une infinité...

Citation:

Envoyé par segolsed

Lorsque je récupére np.max, python me retourne bien entendu une seule valeur

si tu veux récupérer tous les index où ça maximise tu peux utiliser np.argmax(), mais dans le cas présent chacune de tes courbes ne présente qu'un seul maximum apparemment (mais peut-être est-ce un problème d'approximation coté Python ?)

par ailleurs, même en admettant qu'il y ait plusieurs maximums, ils ont tous la même valeur (la valeur maximum donc...) du coup la "courbe" que tu cherches ce sera toujours une droite de type f(x)=debiti.max()
Code:

1 2 3 4 5 6 7 8 9 10 11 12 13 14 def f_trace_courbe (x, y , color, xlim, titre, ylabel, xlabel): plt.plot(x, y, color) plt.title(titre) plt.xlim(xlim) plt.ylabel(ylabel) plt.xlabel(xlabel) color=' bgrcym' for i in range (1,n+1): debiti = omega * R * np.tan(alpha) * np.sin(omega * t + (i * 2 * np.pi / n)) * np.pi * r**2 f_trace_courbe (t, np.abs(debiti), color[i], (0,10), 'Débit de la pompe', 'debiti', 'temps') f_trace_courbe (t, np.full(debiti.shape, np.abs(debiti).max()), 'k', (0,10), 'Débit de la pompe', 'debiti', 'temps') plt.show()
et le résultat :
Pièce jointe 504402