Julia 分形

贡献者： addis

本文处于草稿阶段。
本文缺少预备知识，初学者可能会遇到困难。

　　¹ 这里给出一个 Julia 分形的演示程序，运行结果如图 1 。

未完成：说明未完成

图 1：运行结果

代码 1：julia_set.m

function julia_set
close all;
N = 1000;
scale = 1.5;
x = linspace(-1, 1, N); y = x;
[X, Y] = ndgrid(x, y);
I = julia2(X*scale, Y*scale);
I(I > 1000) = 1000;
surfCart(x, x, I);
caxis([0, 500]);
end

function I = julia2(X, Y)
    c = -0.8 + 0.15745i;
    Z = X + 1i*Y;
    for i=1:80
        Z = Z .* Z + c;
    end
    I = abs(Z); % I = I/max(I(:));
end

1. ^ 参考 Wikipedia 相关页面

致读者：小时百科一直以来坚持所有内容免费，这导致我们处于严重的亏损状态。长此以往很可能会最终导致我们不得不选择大量广告以及内容付费等。因此，我们请求广大读者热心打赏 ，使网站得以健康发展。如果看到这条信息的每位读者能慷慨打赏 10 元，我们一个星期内就能脱离亏损，并保证在接下来的一整年里向所有读者继续免费提供优质内容。但遗憾的是只有不到 1% 的读者愿意捐款，他们的付出帮助了 99% 的读者免费获取知识，我们在此表示感谢。