w-m三维分形函数matlab代码

w-m三维分形函数是一种可以通过迭代生成各种形态的函数，并在三维空间中展示的分形函数。下面是一个使用Matlab编写的w-m三维分形函数的简单代码：

% 设置参数
N = 1000000;   % 迭代次数
scale = 0.05;  % 缩放比例
x0 = 0;        % 初始x坐标
y0 = 0;        % 初始y坐标
z0 = 30;       % 初始z坐标

% 初始化数组
X = zeros(1, N);
Y = zeros(1, N);
Z = zeros(1, N);

% 迭代计算
for i = 2:N
    r = rand();
    if r <= 0.01
        x = 0;
        y = 0.16 * y0;
        z = 0;
    elseif r <= 0.86
        x = 0.85 * x0 + 0.04 * y0;
        y = -0.04 * x0 + 0.85 * y0 + 1.6;
        z = 0.0007 * z0;
    elseif r <= 0.93
        x = -0.15 * x0 + 0.28 * y0;
        y = 0.26 * x0 + 0.24 * y0 + 0.44;
        z = 0.0007 * z0 + 1.6;
    else
        x = 0.2 * x0 - 0.26 * y0;
        y = 0.23 * x0 + 0.22 * y0 + 1.6;
        z = 0.0007 * z0 + 0.44;
    end
    X(i) = x;
    Y(i) = y;
    Z(i) = z;
    x0 = x;
    y0 = y;
    z0 = z;
end

% 绘制图形
scatter3(X, Y, Z, 1, 'black');

这段代码首先设置迭代次数N、缩放比例scale和初始坐标x0、y0、z0。初始化数组X、Y和Z用于存储每次迭代后的坐标值。然后，通过迭代计算得到每次迭代后的新坐标，并将其存储在对应的数组中。最后，使用scatter3函数绘制三维散点图，以展示w-m三维分形函数的形态。