Integration of one array variable wrt other array variable

2 次查看（过去 30 天）
Somnath Kale
Somnath Kale 2021-1-29
编辑： Somnath Kale 2021-2-1
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)

请先登录，再回答此问题。

回答（1 个）

Amanpreetsingh Arora
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.
https://in.mathworks.com/help/matlab/ref/cumtrapz.html
  5 个评论
显示 3更早的评论
Amanpreetsingh Arora
Amanpreetsingh Arora 2021-2-1
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.
Somnath Kale
Somnath Kale 2021-2-1
编辑：Somnath Kale 2021-2-1
can u just add condition from -40 to +40 and wirte small code?

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Numerical Integration and Differential Equations 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by