如何在python中模拟单向Henon映射？

Question

我试着用下面的公式来模拟一个单向的Henon图：

代码如下（e在0和1之间变化（间隔为0.5），a=1.4，bx=by=0.3）：

x = np.zeros(200000)
y = np.zeros(200000)
xsel = []
ysel = []
for i in range(2, 200000):
    x[i] = a - (x[i-1]**2) + (bx*x[i-2])
    y[i] = a - ( ( (e*x[i-1]) + ((1-e)*y[i-1]) )*y[i-1] ) + (by*y[i-2])
    if i >=100000:
xsel.append(x[i])
ysel.append(y[i])

datax.append(xsel)
datay.append(ysel)

它生成200000个数据点，最后100000个数据点用于进一步计算。使用这些数据计算的传递熵显示了意想不到的结果，我怀疑我现在实现了Henon映射模拟。有人能验证我申请模拟的逻辑吗