How to integrate dataset over a certain domain ?

6 次查看(过去 30 天)
clear all
close all
clc
load data.txt
S_v = data(:,2);
T = data (:,1);
I_H = trapz(S_v); %Spectral_Intensity
Here it integrates over the whole domain. But how can I integrate within the specified range ?

采纳的回答

Star Strider
Star Strider 2021-4-22
Perhaps:
load data.txt
S_v = data(:,2);
T = data (:,1);
Trange = (T >= 0.1) & (T <= 2.5);
I_H = trapz(T(Trange), S_v(Trange)); %Spectral_Intensity
will do what you want. (Not able to test it without ‘data.txt’.)
Do not use clear or clc. It is simply inefficient.
  2 个评论
Sumit Saha
Sumit Saha 2021-4-23
@Star Strider Thank you so much. I just want to ask that for such integration using trapz command is completely okay? Because I'm calculating all these things for earthquake engineering.
Star Strider
Star Strider 2021-4-23
As always, my pleasure!
With vectors or arrays, the trapz funciton is the only option. The other integration functions require function handle arguments, and integrate the functions themselves between specific limits.

请先登录,再进行评论。

更多回答(0 个)

类别

Help CenterFile Exchange 中查找有关 Numerical Integration and Differentiation 的更多信息

产品


版本

R2017b

Community Treasure Hunt

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

Start Hunting!

Translated by