# "Subscript indices must either be real positive integers or logicals" error while using trapz function

1 vue (au cours des 30 derniers jours)
ishita agrawal le 19 Mai 2018
Commenté : Image Analyst le 19 Mai 2018
Hello, I am using following code to calculate area under the peak.
D = amplitude(:);
l = 59;
chg(1) = 1;
chg(2) = l;
for k1 = 1:numel(chg)-1
segment_area_1k1(a,j) = trapz(chg(k1:k1+1), D(chg(k1:k1+1)));
end
end
The code upon run shows error "subscript indices must either be real positive integers or logicals". I have attached amplitude data for reference. Can anyone solve this problem?
##### 0 commentairesAfficher -2 commentaires plus anciensMasquer -2 commentaires plus anciens

Connectez-vous pour commenter.

### Réponses (2)

Ameer Hamza le 19 Mai 2018
Modifié(e) : Ameer Hamza le 19 Mai 2018
What are the values of a and j? They must be a positive integer. Also if you just want to calculate the area under the curve, then pass an entire vector to trapz()
area = trapz(D);
If you want to get segmented area, you can avoid for loop altogether. A better approach is to do it like this
segmentedArea = diff(cumtrapz(D))
##### 4 commentairesAfficher 2 commentaires plus anciensMasquer 2 commentaires plus anciens
ishita agrawal le 19 Mai 2018
Actually, this code worked perfectly for other datasets. This time, it is working for some values in the dataset and then this error occurred.
There is something I want to ask. When I tried your code (trapz(D)), the value of calculated area was different than what I calculated using my code. I couldn't understand the reason.
Image Analyst le 19 Mai 2018
Give your complete code, including the code where you read in the data file, and your complete error message, including line numbers and ALL the red text.

Connectez-vous pour commenter.

Image Analyst le 19 Mai 2018
And see this very short tutorial on how to format your code: format your code
##### 0 commentairesAfficher -2 commentaires plus anciensMasquer -2 commentaires plus anciens

Connectez-vous pour commenter.

### Catégories

En savoir plus sur Numerical Integration and Differentiation dans Help Center et File Exchange

### Community Treasure Hunt

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

Start Hunting!

Translated by