Area under the Peak

Hello,

I am using the following code to find the width, height and area of the peak,

load('data.csv')
X = data(:,1);
Y= data(:,2);
>> plot(X,Y)
findpeaks(Y,X,'MinPeakProminence',4,'Annotate','extents')
[pks,loc,w,p] =findpeaks(Y,X);
x = 1:length(data);
x_pks = x(loc);
area = trapz(x_pks, pks)

I am getting an error to find the area,as mentioned below

Array indices must be positive integers or logical values.

please suggest me the any solution.

Answer 1

[Sulaymon Eshkabilov]

Here is the corrected code;

clearvars; clc

load('data.csv')
X = data(:,1);
Y = data(:,2);
plot(X,Y)
findpeaks(Y,X,'MinPeakProminence',4,'Annotate','extents')
[pks,Loc,width,psrom] =findpeaks(Y,X);
findpeaks(Y,X)
x = X;
% if there are more than 1 peak, then it is necessary to specify which one
x1=Loc(1)-width(1)/2;
x2=Loc(1)+width(1)/2;
area = trapz([x1,x2], [pks(1), pks(1)])

[iqra kiran]

Thank you.

[Sulaymon Eshkabilov]

Most welcome