Intégrale d'une fonction définie par un tableau de valeurs
Si le tableau contenant les valeurs de l'abscisse est connu
## Définition de la fonction intégration (
méthode des trapèzes, méthode d'ordre 2)
def integration(x, y, iy0):
# x et y des listes ou des tableaux de valeurs, de même dimension
# iy0 la valeur initiale de l'intégrale, iy le tableau de valeurs de l’intégrale
iy=[iy0]
for i in range(len(x)-1):
aire = (y[i] + y[i + 1]) /2 * (x[i+1] - x[i]) # aire du trapèze élémentaire
iy.append(iy[-1] + aire) # ajoute la valeur "iy[-1] + aire" au tableau iy
return iy
Exemple d'utilisation de cette fonction :
# t, v sont les tableaux des valeurs brutes de même dimension
# x0 est la valeur initiale
# x est le tableau des valeurs intégrées avec la méthode des trapèzes d'ordre 2
x=integration(t,v,x0)
Si le tableau contenant les valeurs de l'abscisse n'est pas connu
## Définition de la fonction intégration (
méthode des trapèzes, méthode d'ordre 2)
def integration(dx, y, iy0):
# y une liste ou un tableau de valeurs,
# dx l'intervalle d'abscisse entre deux valeurs de y
# iy0 la valeur initiale de l'intégrale, iy le tableau de valeurs de l'intégrale
iy=[iy0]
for i in range(len(y)-1):
aire = (y[i] + y[i + 1]) /2 * dx # aire du trapèze élémentaire
iy.append(iy[-1] + aire) # ajoute la valeur "iy[-1] + aire" au tableau iy
return iy
Exemple d'utilisation de cette fonction :
# v est le tableau des valeurs brutes
# dt l'intervalle de temps entre deux valeurs de v
# x0 est la valeur initiale
# x est le tableau des valeurs intégrées avec la méthode des trapèzes d'ordre 2
x=integration(dt,v,x0)