我正在尝试提出一个代码,该代码将使我能够绘制一个图表,以使周期翻倍分叉。
我正在使用等式x = rx − 1(1 − x)
,并试图用0.5到4的R值对其进行建模。这是我正在使用
startr = 0.5
finalr = 4
max_time = 200
x = [0.1]
r= np.linspace(.5,4,200)
for n in range(0,200):
x = np.append(r * x[n] * (1-x[n]))
plt.plot(x, label='x');
plt.xlabel('t');
这不断被踢出去
TypeError: append() missing 1 required positional argument: 'values'
这是numpy.append()
的两个绝对必要的参数,从numpy参考中获取。
arr :将array_like值附加到此数组的副本上。
值: 像array_一样,这些值将附加到ARR的副本上。它一定是 正确的形状(与ARR相同的形状,不包括轴)。如果轴为 未指定,值可以是任何形状,并且将在 使用。
因此,尝试使用
np.append(x, r * x[n] * (1-x[n]))
在您的循环内。
逻辑地图 保存文件并运行,图形的png图像文件将保存在同一文件夹中
import numpy as np
import matplotlib.pyplot as plt
Many =50000
x = np.random.rand(Many)
r = np.linspace(0,4.0, num= Many)
for i in range(1, 54):
x_a = 1-x
Data= np.multiply(x,r)
Data= np.multiply(Data, x_a)
x = Data
plt.title(r'Logistic map: $x_{n+1} = r x_{n} (1-x_{n}).$ n = '+ str(i) )
plt.ylabel('x-Random number')
plt.xlabel('r-Rate')
plt.scatter(r, Data, s=0.1, c='k')
plt.show()
plt.savefig(str(i) + " Logistic Map.png", dpi = 300)
plt.clf()