plot integral
3 views (last 30 days)
Show older comments
good morning
I have this equationy(t)=integral of((C(t)*dt)) (0-->t) i haven't the equation of C(t) but i have the values of C(t) in each minute i need to plot y(t) how i can do it???
thanks in advance
2 Comments
Answers (2)
Sarah Wait Zaranek
on 26 May 2011
This can be done using cumtrapz --
Doc here:
Code example below using points taken from a sine wave:
X = linspace(0,2*pi,100);
Y = sin(X);
Z = cumtrapz(X,Y);
plot(X , [Y; Z])
0 Comments
Charles Cunningham
on 1 Jul 2013
Sarah, thank you for pointing out the use of cumtrapz. That answer also helped me make progress with a scrpit I'm wokring on. I did notice a discrepancy with the result for Z and it took a little bit to figure out why. So I thought I'd share what I found.
The integral of sin(x)dx is -cos(x).
After executing your code I expected to see the plot of Z to match a plot of -cos(x). However, it did not. To illustrate this point plot the following (based on Sarah's code):
plot(X , [Z; -cos(X)])
Notice that the result Z is -cos(x) offset by a constant.
I was then reminded of 2 things:
- The integral of f(x)dx is f(x)/dx + constant.
- That cumtrapz is an esimation to find an integral.
The estimation includes the constant. So you can correct the estimation of this integral by subtracting that constant. Then the proper result will be plotted. Update Z to the following equation:
Z = cumtrapz(X,Y) - 1;
0 Comments
See Also
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!