12345678
返回列表 回复 发帖

博导

Rank: 4Rank: 4

UID
188 
帖子
660 
精华
7 
积分
405 
41#
分形几何发表于 2010-6-13 17:27 | 只看该作者
39# 柳烟


10个8个函数没问题的。
回复 引用

TOP

版主

Rank: 7Rank: 7Rank: 7

UID
539 
帖子
87 
精华
0 
积分
338 
42#
柳烟发表于 2010-6-13 18:51 | 只看该作者
35# 柳烟

复变函数中的对数函数的定义你知道吧？按那里的规定就是了！log(z)=(ln(abs(z)),arg(z))
分形几何 发表于 2010-6-13 17:26
胡老师这一说，想起来了，对的。
回复 引用

TOP

版主

Rank: 7Rank: 7Rank: 7

UID
539 
帖子
87 
精华
0 
积分
338 
43#
柳烟发表于 2010-6-13 18:53 | 只看该作者
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图片：
Fractal1.jpg
未命名1.JPG
此图不对称。
z((a3z + a2)z + a1)z五.gsp (29.5 KB)
回复 引用

TOP

版主

Rank: 7Rank: 7Rank: 7

UID
539 
帖子
87 
精华
0 
积分
338 
44#
柳烟发表于 2010-6-13 19:01 | 只看该作者
上楼文件可能有点问题，与UF当中比较，GSP中的触须怎么没接上呢？可能是我在那个地方出了乱子。
回复 引用

TOP

版主

Rank: 7Rank: 7Rank: 7

UID
539 
帖子
87 
精华
0 
积分
338 
45#
柳烟发表于 2010-6-13 20:26 | 只看该作者
重新做了一遍，最右边接上了，其它地方又没接上，与UF效果相比，最右边GSP中间，还多了一个小M集，怪。不知是不是阈值设置原因，我试试看。改变阈值，仍不行。
但此图对称。
未命名.JPG
回复 引用

TOP

院士

Rank: 5Rank: 5Rank: 5

UID
7 
帖子
2888 
精华
10 
积分
2836 
46#
榕坚发表于 2010-6-13 22:01 | 只看该作者
45# 柳烟


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

捕获343.JPG (40.33 KB)

捕获343.JPG

回复 引用

TOP

版主

Rank: 7Rank: 7Rank: 7

UID
539 
帖子
87 
精华
0 
积分
338 
47#
柳烟发表于 2010-6-13 23:37 | 只看该作者
将45楼的文件造出J集后，扫出一张图片，还勉强。
未命名.JPG
回复 引用

TOP

版主

Rank: 7Rank: 7Rank: 7

UID
539 
帖子
87 
精华
0 
积分
338 
48#
柳烟发表于 2010-6-14 00:26 | 只看该作者
我又重新搞了一遍，这是第三遍了，应该说计算没有问题，怎么还是45楼模样，这是什么原因呢？
未命名.JPG
回复 引用

TOP

版主

Rank: 7Rank: 7Rank: 7

UID
539 
帖子
87 
精华
0 
积分
338 
49#
柳烟发表于 2010-6-14 01:46 | 只看该作者
终于成功了，还是刚才造的文件，若对C先作色，再合并Zq,则触须不交；我先合并z与q，再对C着色，结果完全正常，真是一桩怪事。太迟了，我该睡觉了。
未命名.JPG
614((a3z+a2)z+a1)zM集五.gsp (44.35 KB)
回复 引用

TOP

版主

Rank: 7Rank: 7Rank: 7

UID
539 
帖子
87 
精华
0 
积分
338 
50#
柳烟发表于 2010-6-14 22:29 | 只看该作者
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图片：
Fractal1.jpg
GSP效果：
未命名.JPG
z (((0.25z + a3)z + a2)z + a1)z六.gsp (31.89 KB)
回复 引用

TOP

12345678
返回列表