python - fork 图python - 没有图像？

Question

我正在尝试为以下迭代 map 制作分岔图:x_n+1 = x_n * e^(r(1-x_n)) .

首先我定义了 map :

def newmap(x,r):
    return x*math.exp(r*(1-x))

然后我尝试了这个:

def bifurcation_diagram(rmin=0, rmax=4, r_N=2000, N_min=4000, N = 1000):
    rspace = np.linspace(rmin, rmax, r_N)
    x = 0
    rset = []
    xset = []
    for r in rspace:
        for i in range(N_min + N):
            x = newmap(x,r)
            if i > N_min:
                rset.append(r)
                xset.append(x)
    plt.figure(figsize=(16,7))
    plt.xlim((rmin,rmax))
    plt.ylim((0,5))
    plt.scatter(rset,xset,s=0.3,c='C0',linewidth=0)
    plt.xlabel(r'r', fontsize=20)
    plt.ylabel(r'$x_{end}$', fontsize=29, rotation=0)
    plt.show()

当我尝试时 bifurcation_diagram()我得到一个空白的情节。

我不确定我哪里出错了。

Answer 1

问题是x=0是不动点，x=1 也是不动点.如果您切换 x=0至 x=0.1但除此之外，将其留在原处，第一个 r值(value)驱动x到(这些值是什么)吸引不动点 1 .您需要输入 x=0.1主循环内部:

for r in rspace:
    x = 0.01

(和以前一样)。