Integration of one array variable wrt other array variable
4 views (last 30 days)
Show older comments
Hi
I am trying to integrate the current with respect to time to get polarization value. i have both time and current array can anyone help me with this to get proper integration
I am currently using following code:
Thickness = 15*10^-4;
Area = 350
data = readtable('Good PE Loop.txt')
time = data{:,1}
current= data{:,4}
polarization =data{:,5}
field = data{:,2}/Thickness
plot(time,current)
p1 = cumtrapz(current)
p = (10^5/3.5)*p1
plot(time,p)
plot(field,p)
but the value of polarization which im getteing is not correct. and im unable detect the mess up!
thanks in advance!!
- datafile is attached. (my calculated values should be as in the array polarization)
0 Comments
Answers (1)
Amanpreetsingh Arora
on 1 Feb 2021
If the time spacing for integration is non uniform, then you will also need to pass the time vector to "cumtrapz" or else it assumes unit spacing for integration. Try passing time vector in the following fashion and see if it yields expected results.
p1 = cumtrapz(time,current)
For more information on function "cumtrapz", refer to the following documentation.
5 Comments
Amanpreetsingh Arora
on 1 Feb 2021
Does the integral have an initial condition? If it has then you will have to add that initial value to "p" becuase "cumtrapz" computes the definate integral considering zero as initial condition.
See Also
Categories
Find more on Calculus in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!