# How to integrate the area below the peak when I dont have function of the peak?

9 Ansichten (letzte 30 Tage)
sanjay Krishnamurthy am 23 Mär. 2022
I wanted to integrate the area below the shaded peak as shown in the image. The plot is taken from an experimental data so I dont have a function of my curve. How can I specify the limit and get the area below the shaded peak? The simple "trapz" function output is showing "NaN" and I understand I am not specifying the limits here. Can someone help me here. thank you.
% Integrate area below the peak
plot(AR_T,AR_mwmg) %AR_T and AR_mwmg are data stored in AR_data.mat file
int_AR = trapz(AR_T,AR_mwmg)
##### 3 Kommentare1 älteren Kommentar anzeigen1 älteren Kommentar ausblenden
sanjay Krishnamurthy am 23 Mär. 2022
So, How can I proceed to obtain my objective to find the area below the peak?
AndresVar am 23 Mär. 2022

Melden Sie sich an, um zu kommentieren.

### Akzeptierte Antwort

AndresVar am 23 Mär. 2022
Bearbeitet: AndresVar am 23 Mär. 2022
Use conditional indexing to define the limits of extent for the peak.
clear
% remove nan values (why is there a nan at the end?)
idxNotNan = ~isnan(AR_T) & ~isnan(AR_mwmg);
x = AR_T(idxNotNan);
y = AR_mwmg(idxNotNan);
%findpeaks(y,x,'Annotate','extents','WidthReference','halfheight','NPeaks',1,'SortStr','descend');
%[pkY, pkX, pkWidth, pkProm]=findpeaks(y,x,'NPeaks',1,'SortStr','descend')
idxPk = x>210 & x<325; % more or less...
xPk = x(idxPk);
yPk = y(idxPk);
figure;
plot(x,y,LineWidth=2,DisplayName='data')
hold on;
plot(xPk,yPk,'--',LineWidth=2,DisplayName='peak')
legend('show');
AreaFull=trapz(x,y)
AreaFull = 0.1168
AreaPk=trapz(xPk,yPk)
AreaPk = 0.1166
##### 1 Kommentar-1 ältere Kommentare anzeigen-1 ältere Kommentare ausblenden
sanjay Krishnamurthy am 23 Mär. 2022
Thank you so much for the detailed solution. it works well

Melden Sie sich an, um zu kommentieren.

### Kategorien

Mehr zu Numerical Integration and Differentiation finden Sie in Help Center und File Exchange

### Community Treasure Hunt

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

Start Hunting!

Translated by