Definite integral over discrete data

Question

Zoro 2015-11-25

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/257511-definite-integral-over-discrete-data

评论： Star Strider 2015-11-26

Hi, I have a set of discrete data "A" over frequency "f". I want to integrate "A" over "f" from "f1" to "f2". I don't want to integrate over all frequency range. Using trapz(f,A) it gives me the integral over all frequencies... How can I do that? -Thank you

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Star Strider 2015-11-25

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/257511-definite-integral-over-discrete-data#answer_201165

编辑：Star Strider 2015-11-25

在 MATLAB Online 中打开

Without your data I can only outline the approach:

ix1 = find(f <= f1, 1, 'last');
ix2 = find(f >= f2, 1, 'first');
Alim = A(ix1:ix2);
flim = f(ix1:ix2);
int_f1_f1 = trapz(flim, Alim)

This finds the last index of ‘f’ less than or equal to ‘f1’ and the first index greater than or equal to ‘f2’. It then uses those indices to create ‘flim’ and ‘Alim’.

This is UNTESTED CODE but should work.

2 个评论
显示无隐藏无

Zoro 2015-11-25

It worked... Thanks!

Star Strider 2015-11-25

My pleasure!

请先登录，再进行评论。

Answer 2

Adam 2015-11-25

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/257511-definite-integral-over-discrete-data#answer_201162

编辑：Adam 2015-11-25

trapz( f( startIdx:endIdx ), A( startIdx:endIdx ) )

where startIdx and endIdx are the indices of the frequencies f1 and f2 which you can find by e.g.

startIdx = find( f >= f1, 1 );
endIdx = find( f >= f2, 1 );

if you don't already have them. I use >= there because exact equality tests are not stable with doubles, but you can change the test there to whatever you prefer.

4 个评论
显示 2更早的评论隐藏 2更早的评论

Adam 2015-11-26

Thanks! I've never come across it before.

Star Strider 2015-11-26

My pleasure! You will if you post often enough, although the spam trap has gotten progressively more sophisticated since the spam inundations in the Spring of 2014.

There should be a link on the spam notice you see that tells you to contact MathWorks if you believe the spam trap flagged it erroneously. I always click on it and include the URL of the post in my message, although I suspect several people with privileges to do so see the false positive and remove it, as I did with yours.

请先登录，再进行评论。