Функция trapz
Функция trapz выполняет численное интегрирование
одномерного массива данных методом трапеций.
Первым параметром функция принимает массив значений функции.
Вторым параметром можно передать координаты точек по оси x.
Третьим параметром задается шаг между точками, если координаты не указаны.
Синтаксис
np.trapz(y, [x], [dx])
Параметры
y - массив значений функции в точках интегрирования.
x (необязательный) - массив координат точек по оси x.
Должен иметь ту же длину, что и массив y.
dx (необязательный) - шаг между точками по оси x.
Используется, если параметр x не указан. По умолчанию равен 1.
Пример с шагом по умолчанию
Вычислим интеграл функции y = x на отрезке от 0 до 4 с шагом 1:
import numpy as np
y = np.array([0, 1, 2, 3, 4])
res = np.trapz(y)
print(res)
Результат выполнения кода:
8.0
Пример с заданием координат x
Вычислим интеграл функции y = x^2 на отрезке от 0 до 2 с шагом 0.5:
import numpy as np
x = np.array([0, 0.5, 1, 1.5, 2])
y = np.array([0, 0.25, 1, 2.25, 4])
res = np.trapz(y, x)
print(res)
Результат выполнения кода:
2.75
Пример с заданием шага dx
Вычислим интеграл функции y = sin(x) на отрезке от 0 до π с шагом 0.1:
import numpy as np
dx = 0.1
x = np.arange(0, np.pi + dx, dx)
y = np.sin(x)
res = np.trapz(y, dx=dx)
print(res)
Результат выполнения кода:
1.9998355038874436
