Если вы знаете горизонтальный интервал вашего вектора, вы можете использовать trapz , чтобы интегрировать его без функции. Например, чтобы интегрировать y=sin(x)
от 0 до pi с 358 секциями,
x=0:pi/357:pi;
y=sin(x);
area=trapz(x,y);
Если вы просто используете trapz(y)
, вы получите намного большее число, поскольку предполагается, что расстояние между точками по умолчанию равно 1. Эту проблему можно исправить, умножив на расстояние между точками x:
area=pi/357*trapz(y);