inRm3D: 画板论坛 - Powered by Discuz! Board

标题: 恼人的陷阱问题 [打印本页]

作者: 柳烟 时间: 2013-12-29 16:15 标题: 恼人的陷阱问题

未命名.jpg

这是我新造的z^4+1Ｎ集，变换为z----z+1/z.

z^n-2变换Ｎ集.gsp (20.49 KB) .我这两天遇到的怪事，听我慢慢道来。下面我将原来的旧作谢尔宾斯基三角形发上来。见：疑难分形征解

Sierpinski Triangle II (旧文件).gsp (17.5 KB)

图片附件: 未命名.jpg (2013-12-29 16:15, 170.8 KB) / 下载次数 1567
http://inrm3d.cn/attachment.php?aid=20861&k=b41229debdbf437a852b50ee0ebe83a9&t=1732400592&sid=AsH8KM

附件: z^n-2变换Ｎ集.gsp (2013-12-29 16:15, 20.49 KB) / 下载次数 3400
http://inrm3d.cn/attachment.php?aid=20862&k=b638d619bf24a8ae62eb8b815ac15699&t=1732400592&sid=AsH8KM

附件: Sierpinski Triangle II (旧文件).gsp (2013-12-29 16:18, 17.5 KB) / 下载次数 3390
http://inrm3d.cn/attachment.php?aid=20863&k=3e41d643ba7f02ca89f1e16b776d4220&t=1732400592&sid=AsH8KM

作者: 柳烟 时间: 2013-12-29 16:22

下面将上面两个ＧＳＰ文件的数据加载到一起，将谢尔宾斯基迭代终点Zx,Zy加载到牛集的迭代起点，结果比较正确，但是遇到谢氏三角形外当迭代次数大时，外面数据溢出，还好，谢氏三角形内还算正常。

整合.gsp (30.19 KB)
未命名.jpg

外面白色全是数据溢出。

附件: 整合.gsp (2013-12-29 17:22, 30.19 KB) / 下载次数 3218
http://inrm3d.cn/attachment.php?aid=20865&k=0253cf14c4b7ce08990de5a5fb906a17&t=1732400592&sid=AsH8KM

图片附件: 未命名.jpg (2013-12-29 17:22, 137.65 KB) / 下载次数 1640
http://inrm3d.cn/attachment.php?aid=20866&k=ed2023aa2675676090f97a71a99405ae&t=1732400592&sid=AsH8KM

作者: 柳烟 时间: 2013-12-29 17:30

新造谢氏三角形扫图：
未命名.jpg

新造谢氏三角形.gsp (17.74 KB)
文件是没有问题的。下面我将前面的Ｎ集当陷阱，与此文件整合的结果，怪事出现了，问题究竟出在那里呢？
未命名.jpg

问题文件.gsp (30.05 KB)
不知是那个细节出了问题，导致三角形内部Ｎ集中的et为常数，故如斯。大家有空时，可协助找找，这两天我差点被这分形累死，谢氏三角形造了好几遍，始终未找到原因。

图片附件: 未命名.jpg (2013-12-29 17:30, 159.48 KB) / 下载次数 1593
http://inrm3d.cn/attachment.php?aid=20867&k=a104290b7e0bcc696967702a7079fcdb&t=1732400592&sid=AsH8KM

附件: 新造谢氏三角形.gsp (2013-12-29 17:30, 17.74 KB) / 下载次数 3148
http://inrm3d.cn/attachment.php?aid=20868&k=e07924ae5724cb891d0fd540a6d9eee5&t=1732400592&sid=AsH8KM

图片附件: 未命名.jpg (2013-12-29 17:46, 18.87 KB) / 下载次数 1614
http://inrm3d.cn/attachment.php?aid=20869&k=c2f97567ea07e0117c1e2aec463e2425&t=1732400592&sid=AsH8KM

附件: 问题文件.gsp (2013-12-29 17:49, 30.05 KB) / 下载次数 3290
http://inrm3d.cn/attachment.php?aid=20870&k=afd6d2fecfa6b11d4315b32b8172bfa3&t=1732400592&sid=AsH8KM

作者: 柳烟 时间: 2013-12-29 18:10

又造了一遍，居然内部逸出，活见鬼了。

作者: xiaongxp 时间: 2013-12-29 18:47

2# 柳烟
谢氏三角形外的迭代终点太远，画板计算器已不能处理了。把外部迭代终点的坐标缩小10^-50倍，可见外部的环绕带

整合.gsp (30.98 KB)

附件: 整合.gsp (2013-12-29 18:47, 30.98 KB) / 下载次数 3123
http://inrm3d.cn/attachment.php?aid=20871&k=fc96670059de722d8b45c36599016098&t=1732400592&sid=AsH8KM

作者: 柳烟 时间: 2013-12-29 19:18

5# xiaongxp
谢过向老师，帮忙找问题。新文件与旧谢尔三角形造法一样，不知何做，后面做的都不象样。

作者: 柳烟 时间: 2013-12-29 22:08

刚才我重新按旧文件亦步亦趋作，再重新造牛集，我没有用棣美弗公式的工具作，完全死算，结果不但没有数据溢出，而且十分正常，真搞不懂。整了那么多分形，这个分形算是最捣蛋的。

作者: 榕坚 时间: 2013-12-30 07:47

7# 柳烟

棣莫弗公式没错，但错在那个指数公式，会出现数据无意义（即溢出）。如果M集用那个公式做，一般情况下看不出有什么问题，但在扫描针芒处的放大图中只要放大到10^12左右就出现混沌了。因此，该公式不但会影响图形质量而且还会影响扫描速度。
如：M集，中心（-1.75889131716031173631，0.0191049125280461795458）
放大倍数：4*5.2946279E11。迭代次数：10000

图片附件: 1 (2).JPG (2013-12-30 08:28, 149.32 KB) / 下载次数 1620
http://inrm3d.cn/attachment.php?aid=20875&k=6617d33b4d4aa1c23fce0352447783a1&t=1732400592&sid=AsH8KM

作者: 柳烟 时间: 2013-12-30 09:15

#3楼问题文件已经找到了。一般陷阱下的复分形，将坠入陷阱的复分形初始值由复函数公式算出的值的迭代终点判断是否落入陷阱进行处理。此谢尔宾斯基三角形由于不象复分形那样有个逃逸半径，所以应该用缩放后的值的迭代终点作为判断是否落入陷阱的依据。结果正常了。

作者: 柳烟 时间: 2013-12-30 12:11

图片附件: 未命名.jpg (2013-12-30 17:46, 27.24 KB) / 下载次数 1568
http://inrm3d.cn/attachment.php?aid=20876&k=37d475b66fa15c314d06089d48baddb4&t=1732400592&sid=AsH8KM

作者: 柳烟 时间: 2013-12-30 15:46

图片附件: 未命名.jpg (2013-12-30 17:47, 39.02 KB) / 下载次数 993
http://inrm3d.cn/attachment.php?aid=20877&k=eeb7e28475fdf4df3a9c61f2ff85c924&t=1732400592&sid=AsH8KM

作者: 柳烟 时间: 2013-12-31 09:55

一般来讲，陷阱的开状，决定了掉入陷阱的复分形的构件。如果是球陷阱，则复分形将由大小不同的球排列，挤出复分形的形状。用ＧＳＰ造陷阱下的复分形，为何对有些分形，总是不能达到这样的效果呢？而ＵＦ中的陷阱下的复分形，均遵守这条不变的定律呢？大家探究探究。

作者: 柳烟 时间: 2014-1-1 13:08

ＧＳＰ整不出，ＵＦ中整一个：

图片附件: Fractal2.jpg (2014-1-1 13:08, 188.97 KB) / 下载次数 0
http://inrm3d.cn/attachment.php?aid=20882&k=f528eb464d08018061889c2ba63e1c27&t=1732400592&sid=AsH8KM

作者: 柳烟 时间: 2014-1-5 20:10

GSP实现复分形陷阱，复杂且困难，即使弄出，扫描速度成问题。就是ＵＦ中造出，速度也受影响。在ＵＦ中搞两张图片。

图片附件: Fractal2.jpg (2014-1-5 20:10, 161.12 KB) / 下载次数 993
http://inrm3d.cn/attachment.php?aid=20903&k=6f15ee9cd7ea2a8c903393b2aef06a0b&t=1732400592&sid=AsH8KM

图片附件: Fractal2.jpg (2014-1-5 20:42, 158.28 KB) / 下载次数 978
http://inrm3d.cn/attachment.php?aid=20904&k=10f22021914d4ebf4ba7072519b3f3a7&t=1732400592&sid=AsH8KM

作者: 榕坚 时间: 2014-1-5 20:14

14# 柳烟

让谢氏陷阱对称就更好看了。这构图有趣。

作者: 柳烟 时间: 2014-1-5 23:17

15# 榕坚
我尝试在ＵＦ中调成轴对称，搞了老半天，结果不行，发现随便如何调，竟是个中心对称。我将ucl文件发上来，你试试，如何调成轴对称。
柳烟谢氏三角陷阱文件：
柳烟谢氏三角文件1 {

init:
float stdz = 0
float dz = 0.0
loop:
  trz = exp(flip(-pi/180*@rot))*(#z-@center)/@size
  x = 1/sqrt(@r)*real(trz)
  y = sqrt(@r)*imag(trz)
  if @freq != 0
x = 2*sin(2*pi*@freq*x)
y = 2*sin(2*pi*@freq*y)
  endif
  trz = x + flip(y)
  z1 = trz
  float i = 0
  iter = @niter

  while i < @niter
i = i + 1
  IF ((imag(z1)>=.575*real(z1)) && (-.575*real(z1)<=imag(z1)))
z1 = 2*z1-1i
  ELSEIF (real(z1)<=0)
z1 = 2*z1+(.8660254, 0.5)
  ELSEIF (real(z1)>0)
z1 = 2*z1+(-.8660254, 0.5)
  ELSE
z1 = 2*z1
  ENDIF
   z=z1
if |z| > 1e20
   iter = i
   i = @niter
endif
endwhile
  float logp = 1/log(2)
  float logb = log(log(1e20))
   dz = real(iter + logp*logb - logp*log(cabs(z1)))
if dz >stdz
   stdz = dz
endif
final:
  #index = .1*stdz
default:
  title = "柳烟谢氏三角陷阱"
  helpfile = "sam-help/juliatrap.htm"

param center
caption = "Center"
default = (1,1)
  endparam
  param rot
caption = "Rotation"
default = -30.0
  endparam
  param size
caption = "Size"
default = 0.38750
  endparam
  param r
caption = "Ratio Width/Heigh"
default = 1.0
  endparam
  param @niter
caption = "Julia Iterations"
default = 100
  endparam
  param freq
caption = "Trap Frequency"
default = 0.0
  endparam
}

作者: 柳烟 时间: 2014-1-5 23:59

图片附件: Fractal2.jpg (2014-1-5 23:59, 187.47 KB) / 下载次数 0
http://inrm3d.cn/attachment.php?aid=20905&k=50a031ebe67cfa65b1cdc30b0f3b4842&t=1732400592&sid=AsH8KM

作者: 榕坚 时间: 2014-1-6 07:57

16# 柳烟

是要这效果吗？

图片附件: Fractal1.jpg (2014-1-6 07:57, 116.35 KB) / 下载次数 1032
http://inrm3d.cn/attachment.php?aid=20906&k=34fc607c8eac61dd4750cc7ddad08f3d&t=1732400592&sid=AsH8KM

作者: 柳烟 时间: 2014-1-6 08:59

是这个效果。

欢迎光临 inRm3D: 画板论坛 (http://inrm3d.cn/)