分形几何

39# 柳烟


10个8个函数没问题的。

柳烟

35# 柳烟

复变函数中的对数函数的定义你知道吧？按那里的规定就是了！log(z)=(ln(abs(z)),arg(z))
分形几何 发表于 2010-6-13 17:26

胡老师这一说，想起来了，对的。

柳烟

UF中代码如下：
tepfix {
; Generic Mandelbrot set

init:
  complex a3 = 0.3333333333333333333
  complex opix = 3/(4*#pixel)
  complex xm1 = (#pixel + opix)/sqrt(3)

  complex p = (#pixel - opix)/xm1
  complex q = 2*(xm1 + 1)/p + a3*p
  complex a2 = -0.5*(p + q)
  complex a1 = p*q
  z = q
loop:
  z = ((a3*z + a2)*z + a1)*z
UF图片：

下载 (17.64 KB)

2010-6-13 18:53

下载 (32.49 KB)

2010-6-13 18:56

此图不对称。
z((a3z + a2)z + a1)z五.gsp (29.5 KB)

下载次数: 2549

2010-6-13 18:55

柳烟

上楼文件可能有点问题，与UF当中比较，GSP中的触须怎么没接上呢？可能是我在那个地方出了乱子。

柳烟

重新做了一遍，最右边接上了，其它地方又没接上，与UF效果相比，最右边GSP中间，还多了一个小M集，怪。不知是不是阈值设置原因，我试试看。改变阈值，仍不行。
但此图对称。

下载 (38.04 KB)

2010-6-13 20:26

榕坚

45# 柳烟


是否把sqrt当成平方来算了？还有迭代次数不够。也可能是着色问题。我之前好象做过，右边不可能有一个小M集的。

柳烟

将45楼的文件造出J集后，扫出一张图片，还勉强。

下载 (64.84 KB)

2010-6-13 23:37

柳烟

我又重新搞了一遍，这是第三遍了，应该说计算没有问题，怎么还是45楼模样，这是什么原因呢？

下载 (42.12 KB)

2010-6-14 00:26

柳烟

终于成功了，还是刚才造的文件，若对C先作色，再合并Zq,则触须不交；我先合并z与q，再对C着色，结果完全正常，真是一桩怪事。太迟了，我该睡觉了。

下载 (28.62 KB)

2010-6-14 01:46

614((a3z+a2)z+a1)zM集五.gsp (44.35 KB)

下载次数: 2267

2010-6-14 01:46

柳烟

UF代码：
Trailfix {
; Generic Mandelbrot set
init:
  complex s = (#pixel*#pixel - 2)^(1/3)
  complex p = sqrt(2)*s/#pixel
  complex r = p - s
  complex q = r*(r - 2*p)/(2*r - 6*p)
  complex a3 = -0.3333333333333333333*(p + q + r)
  complex a2 = 0.5*(p*q + p*r + q*r)
  complex a1 = -p*q*r
  z = q
loop:
  z = (((0.25*z + a3)*z + a2)*z + a1)*z
UF图片：

下载 (32.39 KB)

2010-6-14 22:29

GSP效果：

下载 (45.35 KB)

2010-6-14 22:29

z (((0.25z + a3)z + a2)z + a1)z六.gsp (31.89 KB)

下载次数: 2245

2010-6-14 22:29