Discussion :

[segolsed]

Bonjour à tous,
j'ai tracé 6 courbes correspondants aux 6 débits de piston d'une pompe hydraulique

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


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.

[wiztricks]

Salut,

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

[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.

[wiztricks]

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

[BufferBob]

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 :

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

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()

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é...

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 :

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

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 :