Matlab绘制三维Mandelbrot集代码示例

以下是绘制三维Mandelbrot集的Matlab代码：

% 设置参数
n = 200; % 网格尺寸
maxiter = 100; % 最大迭代次数
R = 2; % 区域半径

% 初始化网格
x = linspace(-R, R, n);
y = linspace(-R, R, n);
z = zeros(n,n);
[X, Y] = meshgrid(x, y);

% 迭代计算
for i = 1:n
    for j = 1:n
        c = complex(X(i,j), Y(i,j));
        z(i,j) = mandelbrot(c, maxiter);
    end
end

% 绘制图像
surf(X,Y,z), shading flat
colormap(jet)
view(3)

其中，mandelbrot 函数的实现如下：

function n = mandelbrot(c, maxiter)
z = 0;
n = 0;
while abs(z) <= 2 && n < maxiter
    z = z^2 + c;
    n = n + 1;
end
end

该函数接受一个复数 'c' 和最大迭代次数 'maxiter'，返回在迭代次数不超过 'maxiter' 的情况下，从初值 'z=0' 开始迭代 'z^2 + c' 的次数 'n'。如果 'n' 达到了最大迭代次数但仍然满足 '|z|<=2'，则认为 'c' 属于 Mandelbrot 集合中，否则不属于。