MATLAB Answers

0

Finite Integral with periodic values of x

Asked by Manish Kumar on 12 Mar 2019
Latest activity Commented on by Torsten
on 12 Mar 2019
Dear All,
Range of my x axis is 0.3 to 4.5. I want to integrate my x and y data with x range from 0.3 to 4.5 with the interval of 0.01. First column of the excel sheet is x and 2nd column is y. I write a code as below:
for i = 0.3:4.5:0.01
idx = (x>=i & x<=4.5);
xp = x(idx);
yp = y(idx);
f(i,:)=trapz(xp,yp)
end
I need improvement on this script. Please help.

  0 Comments

Sign in to comment.

1 Answer

Answer by Torsten
on 12 Mar 2019
Edited by Torsten
on 12 Mar 2019
 Accepted Answer

idx = (x>=0.3 & x<=4.5);
xp = x(idx);
yp = y(idx);
F = cumtrapz(xp,yp);
F = F(end) - F;
plot(xp,F)

  2 Comments

i want to calculate the integral by varying the limit. Lower limit is 0.3, 0.31,0.32,0.33....4.5 and upper limit is 4.5
If you insist varying the lower limit by a fixed amount, you can use
index = 0;
for iter = 0.3:0.01:4.5
index = index + 1;
idx = (x>=iter & x<=4.5);
xp = x(idx);
yp = y(idx);
xiter(index) = iter;
Fiter(index) = trapz(xp,yp);
end
plot(xiter,Fiter)

Sign in to comment.