Solved

Integration of time series data

Posted on 2009-07-16
3
1,593 Views
Last Modified: 2016-02-10
Helo- I have columns of time history data as shown below. How do I do integration of these data and write in a separate column, i. e., in"Integrated values" column? Do I have to use trapizoidal rule, or simpson rule? Any example of MATLAB command or code will help. Thanks.

time Value   Integrated values
0        0.1
0.1     2
0.2     3
0.4     1
0
Comment
Question by:okd
[X]
Welcome to Experts Exchange

Add your voice to the tech community where 5M+ people just like you are talking about what matters.

  • Help others & share knowledge
  • Earn cash & points
  • Learn & ask questions
  • 2
3 Comments
 
LVL 15

Expert Comment

by:yuk99
ID: 24873434
What method to use for integration is up to you. Trapezoidal integration is build into Matlab. If you imported time as x and Value as y, you can do integration as z=trapz(x,y);
If you also need help to import/output data, please tell more details on your original data format and what you want as output. How large is the dataset? You can also attach an example here.
0
 

Author Comment

by:okd
ID: 24896016
Thank you. Could you also send me how to import and output data in MATLAB. The database is like I have sent you, but it has 500 rows. I want to plot the time vs. z=trapz(x,y) values. Matlab commands for plotting and import/output data will help the most. Thanks again, - okd
0
 
LVL 15

Accepted Solution

by:
yuk99 earned 125 total points
ID: 24897374
Here is the code. It read the data from tab-delimited file tsdata.txt and returns 2-column matrix data. To make it clear I put time into time vector and value into value vector (rows 2 and 3). TRAPZ function will return only one value - the total area under the curve for the whole vector. To get integrated value for every time point you can use CUMTRAPZ.
% data import
data=dlmread('tsdata.txt','\t',1,0);
% let's split variables for clarity
time=data(:,1);
value=data(:,2);
% calculating integrated values
intvalue=cumtrapz(time,value);
% plotting
plot(time,intvalue)

Open in new window

0

Featured Post

Independent Software Vendors: We Want Your Opinion

We value your feedback.

Take our survey and automatically be enter to win anyone of the following:
Yeti Cooler, Amazon eGift Card, and Movie eGift Card!

Question has a verified solution.

If you are experiencing a similar issue, please ask a related question

Suggested Solutions

Title # Comments Views Activity
Choosing the right language for new project 8 66
numbers ascending pyramid 101 240
Way to decrease size of apk file 9 109
troubleshoot a python script 8 74
There is an easy way, in .NET, to centralize the treatment of all unexpected errors. First of all, instead of launching the application directly in a Form, you need first to write a Sub called Main, in a module. Then, set the Startup Object to th…
This article is meant to give a basic understanding of how to use R Sweave as a way to merge LaTeX and R code seamlessly into one presentable document.
This tutorial explains how to use the VisualVM tool for the Java platform application. This video goes into detail on the Threads, Sampler, and Profiler tabs.
The viewer will learn how to clear a vector as well as how to detect empty vectors in C++.

749 members asked questions and received personalized solutions in the past 7 days.

Join the community of 500,000 technology professionals and ask your questions.

Join & Ask a Question