How can I calculate three different shaded areas?

조회 수: 1 (최근 30일)
wd w
wd w 2024년 9월 20일
댓글: Star Strider 2024년 9월 20일
I have three groups of sample points as attached, and intend to calculate three different shaded areas under curve from x = -40 to x = 30 by using for loop.

채택된 답변

Star Strider
Star Strider 2024년 9월 20일
Try this —
T1 = readtable('example.xlsx')
T1 = 6171x6 table
x1 y1 x2 y2 x3 y3 _______ ______ _______ ______ _______ ______ -69.644 10.562 -69.643 10.583 -69.631 10.605 -69.641 10.574 -69.641 10.596 -69.629 10.618 -69.639 10.586 -69.639 10.608 -69.626 10.631 -69.636 10.599 -69.636 10.621 -69.624 10.644 -69.634 10.611 -69.634 10.633 -69.622 10.657 -69.632 10.623 -69.632 10.645 -69.619 10.67 -69.629 10.636 -69.629 10.658 -69.617 10.683 -69.627 10.648 -69.627 10.67 -69.614 10.696 -69.625 10.66 -69.625 10.682 -69.612 10.709 -69.622 10.673 -69.622 10.694 -69.609 10.721 -69.62 10.685 -69.62 10.707 -69.607 10.734 -69.615 10.71 -69.615 10.732 -69.602 10.758 -69.61 10.734 -69.61 10.758 -69.597 10.783 -69.606 10.759 -69.605 10.784 -69.593 10.807 -69.601 10.784 -69.6 10.81 -69.588 10.831 -69.596 10.809 -69.595 10.836 -69.583 10.856
for k = 1:2:size(T1,2)
xv = T1{:,k};
yv = T1{:,k+1};
idxrng = (xv >= -40) & (xv <= 30);
AUC(ceil(k/2)) = trapz(xv(idxrng), yv(idxrng));
end
AUC
AUC = 1×3
98.2744 90.7472 80.9077
<mw-icon class=""></mw-icon>
<mw-icon class=""></mw-icon>
figure
tiledlayout(3,1)
for k = 1:2:size(T1,2)
% k
xv = T1{:,k};
yv = T1{:,k+1};
nexttile
plot(xv, yv)
hold on
idxrng = (xv >= -40) & (xv <= 30);
patch([xv(idxrng); flip(xv(idxrng))], [zeros(nnz(idxrng),1); flip(yv(idxrng))], 'r', 'FaceAlpha',0.5)
hold off
grid
end
.
  댓글 수: 2
wd w
wd w 2024년 9월 20일
Star Strider, thank you very much!
Star Strider
Star Strider 2024년 9월 20일
As always, my pleasure!

댓글을 달려면 로그인하십시오.

추가 답변 (0개)

카테고리

Help CenterFile Exchange에서 Numerical Integration and Differentiation에 대해 자세히 알아보기

태그

제품


릴리스

R2024b

Community Treasure Hunt

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

Start Hunting!

Translated by