Python中的分叉和Lyapunov指数

Question

在涉及python方面，我是一个相对入门的人，目前我正在尝试找出一些python解决我遇到的问题。 我正在尝试计算应该创建的分歧图的lyapunov指数。

方程为X_（n + 1）= a sin（pi x_（n））， 其中a = 0.9（用于计算指数）

这是目前我为创建值数组而设置的代码。

import numpy as np
np.set_printoptions(threshold=np.nan)
import matplotlib.pyplot as plt

a = np.linspace(0,1) 
xn = np.array([.001], dtype = float)

for i in range(0,10000):
    y = a*np.sin(np.pi*xn[i])
    xn = np.append(xn,y)

plt.plot(a,xn[-1])

但是，很显然，当我绘制xn时，我只是出现了点混乱的点，而不是分叉图。我希望我可以得到一些指导，以朝着正确的图前进，我希望可以用它来接近最终目标。

感谢您的帮助，谢谢！

Answer 1

我不确定您要完成的工作，而且我对分叉的了解还不足以独自解决这个问题，但我能够得到一些可行的方法。主要警告是，如果alpha起始于小于0.158，则不会产生正确的输出。

sub-delim

产生图：