c_2331.m

来自「分形与图形设计」· M 代码 · 共 21 行

a=[0 0 0 0 0.18 0 0 0 0 0 0 0 0.1
    0.83 0 0 0 0.86 0.1 0 -0.12 0.84 0 1.62 0 0.85
    0.22 -0.23 0 0.24 0.22 0 0 0 0.32 0 0.82 0 0.07
    -0.22 0.23 0 0.24 0.22 0 0 0 0.32 0 0.82 0  0.07];
x0=0;y0=0;z0=0;
axis tight
for i=1:5000
    r=rand;
    total=a(1,13);
    k=1;
    while (total<r)
        k=k+1;
        total=total+a(k,13);
    end
    x1=a(k,1)*x0+a(k,2)*y0+a(k,3)*z0+a(k,10);
    y1=a(k,4)*x0+a(k,5)*y0+a(k,6)*z0+a(k,11);
    z1=a(k,7)*x0+a(k,8)*y0+a(k,9)*z0+a(k,12);
    x0=x1;y0=y1;z0=z1;
    plot3(x1,y1,z1);
    hold on;
end