plot integral
4 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
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 commentaires
Réponses (2)
Sarah Wait Zaranek
le 26 Mai 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 commentaires
Charles Cunningham
le 1 Juil 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 commentaires
Voir également
Catégories
En savoir plus sur Numerical Integration and Differentiation dans Help Center et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!