Python代码绘制Gumowski-Mira分形图

以下是使用Python代码绘制Gumowski-Mira公式分形图的示例。

import numpy as np
import matplotlib.pyplot as plt

def gumowski_mira(x, a, b, c, d):
    return (b * x + a * x**3) / (1 + c * x**2) + d * x

def generate_gumowski_mira(a, b, c, d, x0, n):
    x = np.zeros(n)
    x[0] = x0
    for i in range(1, n):
        x[i] = gumowski_mira(x[i-1], a, b, c, d)
    return x

def plot_gumowski_mira(a, b, c, d, x0, n):
    x = generate_gumowski_mira(a, b, c, d, x0, n)
    plt.plot(x[:-1], x[1:], ',')
    plt.title('Gumowski-Mira Fractal')
    plt.xlabel('x')
    plt.ylabel('y')
    plt.show()

plot_gumowski_mira(0.99, 0.002, 0.1, 0.1, 0.1, 100000)

在这个示例中，我们使用Gumowski-Mira公式生成分形图。其中，'a'、'b'、'c' 和 'd' 是参数，'x0' 是初始值，'n' 是生成的点数。函数generate_gumowski_mira 用于生成数据，而 plot_gumowski_mira 函数用于绘制分形图。

本示例使用了以下参数：'a'=0.99，'b'=0.002，'c'=0.1，'d'=0.1，'x0'=0.1，'n'=100000。您可以尝试不同的参数值来生成不同的分形图。