How can I calculate trapz value based on an event counter and then plot it?

1 次查看（过去 30 天）
Dennis
Dennis 2022-9-13
评论： Dennis 2022-9-16
Hey there,
the following lines of code should implement a calculation of an integral for certain events in the data.
The very first colum is the time, the second the corresponding value and the third column the event counter.
The values are saved in the integ_mtrx. It worked with cumtrapz, but i am having my problems with the trapz function.
Here the events 0 should be calculated. Events with 0 in the second colum shouldn't be calculated ( like event 1,2,3,4,..)
for every event the integral should be calculated:
for i= 0:(vzwcounter)
clear temp;
temp = Integ_MTRX(Integ_MTRX(:,3)==i,:);
B(i,1)= trapz(temp(:,1),temp(:,2));
end
  4 个评论
显示 2更早的评论
Jan
Jan 2022-9-14
I'm lost. Should I assume that the first image contains the contents of Integ_MTRX? What is the value of vzwcounter? What is the contents of the second image and why is this an integral?
Which are the "counter variables in the 1st image"? What do you expect trapz(temp(:,1),temp(:,2)) to reply, if temp is a [1x2] vector?
Can you post some code which runs by copy&paste?
Dennis
Dennis 2022-9-14
I attached you the Integ_MTRX and vzwcounter data.
Both images are from Integ_MTRX and show a snapshot from the data.
Based on the data I need to calculate the integral.
trapz(temp(:,1),temp(:,2)) should be skipped whenever the value in the third column, which are the values of vzwcounter, is changing for each entry. vzwcounter is just counting up events in a different function. for each event the integral should be calculated (e.g. event 175 in the second image, first column is the time (x) and second column are the values (y)).

请先登录，再进行评论。

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

采纳的回答

Jan
Jan 2022-9-15
As far as I understand your code was almost working. I added an if condition to skip scalar data and shifted the index for the output B by 1 because 0 is no valid index in Matlab:
B = nan(1, vzwcounter + 1);
for k = 0:vzwcounter
temp = Integ_MTRX(Integ_MTRX(:, 3) == k, :);
if size(temp, 1) > 1
B(k + 1) = trapz(temp(:, 1), temp(:, 2), 1);
end
end

更多回答（0 个）

类别

在 Help Center 和 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