how to intergrate my acceleration vs time curve to velocity vs time curve with simpson's rule, instead of using cumtrapz?
8 views (last 30 days)
Show older comments
i have a data set of acceleration vs time, i need to intergrate it into velocity vs time and also displacement vs time. anyone can help me with numerical simpson's rule instead of using cumtrapz?
0 Comments
Answers (1)
Thiago Henrique Gomes Lobato
on 12 Apr 2020
You need to implement the simpson's rule on your own. A naive implementation would be something like that
f = @(t)t.^2;
N = 1000;
ti = 0;
tend = 5;
h=(tend-ti)/N;
t=linspace(ti,tend,N);
F = 0;
for i = 1:N-1
F = F + (h/6)*(f(t(i))+(4*f((t(i)+t(i+1))/2))+f(t(i+1)));
end
% To verify the result
AnalyticalResult = tend^3/3;
Error = (AnalyticalResult-F)/AnalyticalResult
Error =
1.0000e-03
0 Comments
See Also
Products
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!