Maximizing value for an array with integral function.

2 次查看(过去 30 天)
I have 'time(t)' and 'acceleration(A)' values. t1 and t2 are the lower and upper limits, where (t2-t1 = 10).
f = max{((integration(A,t,t1,t2))^2.5) * (t2-t1)}
Please help me to put this function in my matlab code.
for t1=1:length(t)
for t2=t1:t1+10
f = (((int(A,t,t1,t2)))^2.5)*(t2-t1);
end
end
I am using this for loop, but getting error for integration.

采纳的回答

Torsten
Torsten 2018-2-22
编辑:Torsten 2018-2-22
fmax = -Inf;
for t1=1:numel(t)-10
fmax = max((trapz(t(t1:t1+10),A(t1:t1+10)))^2.5*(t(t1+10)-t(t1)),fmax);
end
fmax
Best wishes
Torsten.
  8 个评论
Lal
Lal 2018-2-27
Thank you Torsten. I understand now where I am doing mistake. Your suggestion helped a lot.
And one more thing in the result t(i) and t(i+15) is fixed to "t(i+15) - t(i) = 15".
But It can be any point between t(i+15) and t(i).
For example:
t_low = 7
t_up = 10
for f_max = 100
(because of this reason I introduced "i" and "j". As you said its two-dimensional array )
Thanks
Torsten
Torsten 2018-2-27
f_max = -Inf;
for i = 1:85
for j = 1:15
f = (trapz(t(i:i+j),A(i:i+j)))^2.5*(t(i+j)-t(i))
if f > f_max
tlow_max = t(i);
tup_max = t(i+j);
f_max = f;
end
end
end
f_max
tlow_max
tup_max

请先登录,再进行评论。

更多回答(0 个)

类别

Help CenterFile Exchange 中查找有关 Numerical Integration and Differentiation 的更多信息

Community Treasure Hunt

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

Start Hunting!

Translated by