Matlab 代码实现 Gumowski-Mira 分形图

以下是实现 Gumowski-Mira 公式的 Matlab 代码：

% Gumowski-Mira 公式的 Matlab 代码

% 参数设置
a = 0.999;  % 控制分形结构的参数
b = 0.2;    % 控制非线性特性的参数
n = 10000;  % 迭代次数

% 初始化
x = zeros(1, n);
y = zeros(1, n);
x(1) = 0.1;
y(1) = 0;

% 迭代
for i = 2:n
    x(i) = a * x(i-1) + 2 * b * y(i-1) * (1 - x(i-1)^2);
    y(i) = -b * x(i-1) + a * y(i-1) + 0.1;
end

% 绘图
plot(x, y, '.');
title('Gumowski-Mira 分形图');
xlabel('x');
ylabel('y');

运行上述代码即可得到 Gumowski-Mira 分形图。您可以根据需要调整参数 a、b 和迭代次数 n，以生成不同的分形结构。