What do these steps do?

Question

Vineet Patel 2018-12-2

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/433271-what-do-these-steps-do

回答： Walter Roberson 2018-12-2

采纳的回答： Walter Roberson

在 MATLAB Online 中打开

Hi All,

Can someone tell me what these steps mean/do? And if theres a way of simplifing the code.

xmax = [];
[t,x] = ode45(f,tspan,x0);
n = find(t>60);
x = x(n,:);
n = length(x(:,1));
k = 1;
k = k+1;
j = 1;
for i=2 : n-1
        if (x(i-1,1)) < x(i,1) && x(i,1) > (x(i+1,1))
            xmax(k,j)=x(i,1);
            j=j+1;
        end
end

Any help would be great,

Thank you

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Walter Roberson 2018-12-2

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/433271-what-do-these-steps-do#answer_350049

Run an ode45 over tspan. Some but not all of the timesteps correspond to times greater than 60; extract the ode values that were generated for those times.

That done, look for local peaks in the first column.

There are other ways to look for peaks, including vectorized ways. findpeaks() is more flexible about what a peak looks like.