Specific integral of a vector
3 visualizaciones (últimos 30 días)
Mostrar comentarios más antiguos
I have a graph that is plotted from two vectors, val_x and x. it is plotted below as plot(x,val_x). I want to integrate this plot from a to b and get the shaded area under the graph, even when a and b are different from the values in x (example: x is increasing from 0 to 1 with steps of 0.1 and b=0.15 a=0.25). I have tried to represent the vector as a function through polynomial regression, but the function could not catch the accurate shape of the plot.
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/148296/image.png)
0 comentarios
Respuestas (1)
Christiaan
el 25 de Mzo. de 2015
Editada: Christiaan
el 25 de Mzo. de 2015
Dear Ser/Madame,
To create more points in your function, you could use an interpolation function like pchip or spline to create more points. A example how it works is:
clc; clear all; close all;
a = 0.5;
b = 0.7;
%%PART 1: not enough points
h1 = subplot(2,1,1);
x = linspace(0,1,20);
y = 1/2*sin(x*pi).*exp(x*pi);
xpoints=x(find(x>=a & x<=b));
ypoints=y(find(x>=a & x<=b));
plot(h1,x,y);axes(h1); hold on;
area(xpoints,ypoints); hold on;
line([a a],get(h1,'YLim'),'Color',[1 0 0]);
line([b b],get(h1,'YLim'),'Color',[1 0 1]);
title(h1,'title 1');
%%PART 2: more points
h2 = subplot(2,1,2);
xnew = linspace(min(x),max(x),length(x)*50);
ynew = pchip(x,y,xnew);
xpoints_new=xnew(find(xnew>=a & xnew<=b));
ypoints_new=ynew(find(xnew>=a & xnew<=b));
plot(h2,xnew,ynew);axes(h2); hold on;
area(xpoints_new,ypoints_new); hold on;
line([a a],get(h1,'YLim'),'Color',[1 0 0]);
line([b b],get(h1,'YLim'),'Color',[1 0 1]);
title(h2,'title 2');
If you like you can increase the number 50 to any number you like.
Good luck! Christiaan
1 comentario
graadiggrei
el 26 de Mzo. de 2015
Thank you very much! I chose the xnew values to go from "a" to "b" with a certain step, and then integrating it normally with trapz to get the answer I was looking for.
Somehow Matlab central has disconnected the question from my user, so I am unable to accept this as the answer to the question, but it solved the problem!
Ver también
Categorías
Más información sobre Numerical Integration and Differentiation en Help Center y File Exchange.
Productos
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!