【混沌动力学引论】Logistics映射（抛物线映射）的分岔图绘制

步骤

1. 定义Logistics映射的函数
2. 分别对参量$\mu$和变量$x$进行迭代，得出两个维度的数值列表
3. 将两个维度的数值列表输入plt.plot()函数中

代码

# 画Logistics混沌分岔图
import numpy as np
import matplotlib.pyplot as plt

# Logistics方程（抛物线映射）
def logistics(x, mu):
    return 1 - mu * x ** 2

# mu参数范围
muArr = np.arange(0, 2, 0.01)
# x迭代值列表
x2dArr = [];
# 对mu, x进行迭代
for mu in muArr:
    x = 1
    xArr = [];
    for j in range(0, 500):
        # 去掉过渡期的迭代值
        if j > 200:
            x = logistics(x, mu)
            xArr.append(x)
    x2dArr.append(xArr)

# 绘图
plt.plot(muArr, x2dArr, ',k', alpha=0.25)
plt.show()

运行结果：