两个混沌分叉图
前提:
import matplotlib.pyplot as plt
import numpy as np
singer map
python代码
def singermap10(x0, a, b, iters):
x = x0
u = np.arange(a, b, 0.001)
for i in range(iters):
x = u * (7.86 * x - 23.31 * x ** 2 + 28.75 * x ** 3 - 13.302875 * x ** 4)
if i >= iters - 100:
plt.plot(u, x, ',k')
plt.xlabel('u')
plt.ylabel('x')
plt.title('Singer map')
plt.show()
singermap10(0.2, 0.9, 1.08, 500)
运行结果:
logistic map
python代码
def logisticmap7(a, b, x0, iters):
"""
logistic映射绘图
:param a: 控制参数下边界
:param b: 控制参数上边界
:param x0: x的初值
:param iter: 迭代次数
:return: 图像
"""
x = x0
u = np.arange(a, b, 0.001)
for i in range(iters):
x = u * x * (1 - x)
if i >= iters - 100:
plt.plot(x, ',k')
plt.xlabel('u')
plt.ylabel('x')
plt.title('Logistic map')
plt.show()
logisticmap7(2, 4, 0.2, 500)
运行结果: