返回列表 回复 发帖

恼人的陷阱问题

未命名.jpg
这是我新造的z^4+1N集,变换为z----z+1/z.
z^n-2变换N集.gsp (20.49 KB) .我这两天遇到的怪事,听我慢慢道来。下面我将原来的旧作谢尔宾斯基三角形发上来。见:疑难分形征解
Sierpinski Triangle II (旧文件).gsp (17.5 KB)
下面将上面两个GSP文件的数据加载到一起,将谢尔宾斯基迭代终点Zx,Zy加载到牛集的迭代起点,结果比较正确,但是遇到谢氏三角形外当迭代次数大时,外面数据溢出,还好,谢氏三角形内还算正常。
整合.gsp (30.19 KB)
未命名.jpg
外面白色全是数据溢出。
新造谢氏三角形扫图:
未命名.jpg
新造谢氏三角形.gsp (17.74 KB)
文件是没有问题的。下面我将前面的N集当陷阱,与此文件整合的结果,怪事出现了,问题究竟出在那里呢?
未命名.jpg
问题文件.gsp (30.05 KB)
不知是那个细节出了问题,导致三角形内部N集中的et为常数,故如斯。大家有空时,可协助找找,这两天我差点被这分形累死,谢氏三角形造了好几遍,始终未找到原因。
又造了一遍,居然内部逸出,活见鬼了。
2# 柳烟
谢氏三角形外的迭代终点太远,画板计算器已不能处理了。把外部迭代终点的坐标缩小10^-50倍,可见外部的环绕带
整合.gsp (30.98 KB)
5# xiaongxp
谢过向老师,帮忙找问题。新文件与旧谢尔三角形造法一样,不知何做,后面做的都不象样。
刚才我重新按旧文件亦步亦趋作,再重新造牛集,我没有用棣美弗公式的工具作,完全死算,结果不但没有数据溢出,而且十分正常,真搞不懂。整了那么多分形,这个分形算是最捣蛋的。
7# 柳烟


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

1 (2).JPG (149.32 KB)

1 (2).JPG

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