如何實現分段數值積分？

Question

在Matlab中是否有任何方法來縮短這個函數或者是否有任何內置函數可以用來代替這個？下面的函數適用於我，但我想從我的腳本中排除它。

function p = nintegrate(x, fx) 
    acc = 0; 
    for i = 1:length(x)-1 
      p(i) = acc; 
      delta_x = x(i+1) - x(i); 
      acc = acc + delta_x * fx(i); 
    end 
end 

Answer 1

看看cumtrapz，它做了類似的事情。它使用了trapezoidal integration，這是一個稍微好一些的數值積分方法。在你的情況下的使用應該是（未經測試）：

acc = cumtrapz(x, fx);