0 个投票

I was wondering if there was a simpler/faster way of calculating the area between the two curves at each point (both of which are normailsed so the total area is 1). I am currently using a simple for loop to get the area between each set of points for the two curves, then subtracting it to get the difference i.e.
for i = 1:length(normpeak)-1
resN(i) = trapz(time(i:i+1),peak1(i:i+1));
resF(i) = trapz(time(i:i+1),peak2(i:i+1));
end
res = resN-resF;
This works, but unfortunately it is very slow. I was reading up on the trapz function and from what I understood it calculates the total integral by performing something very similar to what i'm doing now (see image) but sums the data at the end.
Is there a way to stop trapz from performing this summing of the data, so I can speed up my code?

请先登录，再回答此问题。

 采纳的回答

Star Strider
Star Strider 2015-4-23

0 个投票

Use the cumtrapz function to get the cumulative integral at every point.
I don’t have your data so I can’t test it, but this approach could work:
res = cumtrapz(time, peak1) - cumtrapz(time, peak2);

更多回答（1 个）

Binbin Qi
Binbin Qi 2015-4-23

0 个投票

if you can get the points of every curve, you can use polyarea function to get the area, and then diff

1 个评论
显示 -1更早的评论 隐藏 -1更早的评论

Christian Reece
Christian Reece 2015-4-23
Thanks for the answer, but this doesn't seem to work. I think I was too vague before but what I need is the difference in area at each point on the curve.

请先登录，再进行评论。

类别

在 帮助中心 和 File Exchange 中查找有关 Numerical Integration and Differentiation 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by