時系列データの積分

x = t （一次関数）の積分を考える。

累積和 z=cumsum(x)
累積台形数値積分 z2=cumtrapz(x)

clear; close all;
t = 0:20;
x= t;
y  = 0.5*x.^2;
z = cumsum(x);
z2 = cumtrapz(x);
plot(t,y);
hold on;
plot(t,z);
plot(t,z2);
legend('analytical','Cumsum','Cumtrapz')

cumtrapzで理論解に一致