Fonctions par partie

Bonjour,

J'aimerais créer un graphique où je représente une courbe U dont l'équation est soit une fonction f soit une fonction g selon la valeur prise par la courbe U
(la fonction f correspond à une équation de charge d'un condensateur de type E*(1-exp(-t/tau)) alors que g correspond à une décharge (de type E*exp(-t/tau))

J'ai testé ceci (dans un premier pour tester l'afficahge selon la valeur du temps) mais j'obtiens une erreur (The truth value of an array with more than one element is ambiguous. Use a.any() or a.all())

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 31 32

import math from numpy import arange import numpy as np import matplotlib.pyplot as plt   def Ucharge(t,tau): return E*(1-np.exp(-t/tau))   def Udecharge(t,tau): return E*np.exp(-t/tau)   t0, tf=0.,200 N=1000 E=10   t=np.linspace(t0,tf,N+1)   R1=10e3 C=1e-3   tau1=R1*C print(tau1)   R2=10 tau2=R2*C print(tau2)   if t<100: Uc=Ucharge(t,tau1)   plt.figure plt.plot(t,Uc)

---

Si quelqu'un pouvait m'aider, ce serait sympa

L'erreur "The truth value of an array with more than one element is ambiguous. Use a.any() or a.all()" se produit car la condition t < 100 est appliquée à l'ensemble de l'array t, qui contient plus d'un élément.

Il n'est pas possible de traiter l'ensemble de l'array t comme une seule valeur booléenne, ce qui crée une ambiguïté. Vous pouvez utiliser la fonction np.where() pour créer un nouvel array en fonction de la condition t<100 :

Cela crée un nouvel array Uc avec la même taille que t, où chaque élément est égal à Ucharge(t, tau1) si la condition t<100 est vraie, sinon l'élément est égal à 0.

Ensuite, vous pouvez tracer Uc en utilisant plt.plot(t, Uc).

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 31 32

import math from numpy import arange import numpy as np import matplotlib.pyplot as plt   def Ucharge(t,tau): return E*(1-np.exp(-t/tau))   def Udecharge(t,tau): return E*np.exp(-t/tau)   t0, tf = 0., 200 N = 1000 E = 10   t = np.linspace(t0, tf, N+1)   R1 = 10e3 C = 1e-3   tau1 = R1*C print(tau1)   R2 = 10 tau2 = R2*C print(tau2)   Uc = np.where(t<100, Ucharge(t, tau1), 0) plt.figure() plt.plot(t, Uc) plt.show()

---

Notez que j'ai également ajouté plt.show() à la fin pour afficher le graphique tracé.