如何求解具有时滞的Henon的lyapunov指数，能否发个程序参考一下

galaxy0311li

求解具有时滞的Henon的lyapunov指数， 例如
%   x(k+1)=-p*x(k-3)*x(k)+y(k-2)+1
%   y(k+1)=q*x(k-1)
如何在程序中描述为相应的迭代系统，写出线性化jacob矩阵，请指教，最好发个相关求解程序

galaxy0311li

x(n+1)=a*x(n-1)+p*tanh(x(n-4))-sin(y(n-2))
y(n+1)=a*y(n-1)+b21*tanh(x(n-4))+P*sin(y(n-2))

a=0.2; p=-0.4; b21=0:0.01:3;

这个系统怎么求lyapunov指数？？

liliangbiao

回复 1 # galaxy0311li 的帖子

从时滞入手，
y(k)=x(k-3);
y(k+1)=x(k-2);
y(k+2)=x(k-1);
....
一次类推即可。得到的是一个多维方程组，而不是简单的二维方程。