小编Joh*_*ohn的帖子

我有一组周期性的数据(但不是正弦的).我在一个向量中有一组时间值,在第二个向量中有一组幅度.我想快速估算一下这个函数的周期.有什么建议？

具体来说,这是我目前的代码.我想近似矢量x(:,2)对矢量t的周期.最后,我想在很多初始条件下执行此操作并计算每个条件的周期并绘制结果.

function xdot = f (x,t)
         xdot(1) =x(2);
         xdot(2) =-sin(x(1));
endfunction

x0=[1;1.75];     #eventually, I'd like to try lots of values for x0(2)
t = linspace (0, 50, 200);


x = lsode ("f", x0, t)

plot(x(:,1),x(:,2));

谢谢!

约翰