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)