MATLAB Answers

trapz whe x is a matrix

3 ビュー (過去 30 日間)
Jose Costa
Jose Costa 2019 年 11 月 19 日
コメント済み: dpb 2019 年 11 月 22 日
I have a Y matrix size 2500x11 were each line corresponds to the 11 vales of curve points for trapezoidal integration.
the x values are not uniform and they vary so i have a second matrix for X size 2500x11.
I' tring to integrate using trapz (X,Y,2) but I get error "X must be a vector.".
How ca one use trapz with a matrix for x values? In octave it works.

  0 件のコメント

Sign in to comment.

採用された回答

dpb
dpb 2019 年 11 月 19 日
TMW only implemented variable X spacing for the full array, not by column (or row). Have to loop to do it...worthy of an enhancement request for the added functionality it would seem.
N=size(Y,1);
for i=1:N
Q(i)=trapz(X(:,i),Y(:,i));
end

  3 件のコメント

Jose Costa
Jose Costa 2019 年 11 月 22 日
Thanks, But loops are always very slow. Prefere in this case to explicitly code the trapezoidal method and use matricial calculation.
In the 50x11 example, loop time was 0.0456 and matrix time was 8.6840e^-4
An enhancement would be good. In Octave it exists.
Jose Costa
Jose Costa 2019 年 11 月 22 日
sorry, 2500x11.
dpb
dpb 2019 年 11 月 22 日
MATLAB is not Octave! Converse is also true... :)
Feel free to submit the enhancement request...I'd suggest referencing this thread if do so.

Sign in to comment.

その他の回答 (0 件)

タグ

製品


リリース

R2019b

Translated by