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

标题: Uf中的分形 [打印本页]

作者: changxde 时间: 2010-11-29 16:52 标题: Uf中的分形

图片附件: Uf-aho.jpg (2010-11-30 17:59, 27.03 KB) / 下载次数 5758
http://forums.netpad.net.cn/attachment.php?aid=8578&k=3579420fed5fa984076f67a4ae0e1b60&t=1779870543&sid=d8ctJd

附件: Uf-aho-CLJ.gsp (2010-11-29 16:54, 52.41 KB) / 下载次数 5814
http://forums.netpad.net.cn/attachment.php?aid=8579&k=f42987c88ff9f26a467bb89601603feb&t=1779870543&sid=d8ctJd

作者: xiaongxp 时间: 2010-11-29 21:44

这色调真抢眼，舒服！

作者: changxde 时间: 2010-11-30 17:34

两个简短的Uf小程序：
cha001{
init:
  z=pixel  ,c=0
loop:
  c=(1-c)*real(sqrt(c))/5+z*.75
  z=z*z+c
bailout:
  |z| <= 4
default:
  periodicity = 0
  maxiter = 100
  magn = 1
  center = (-.4,0)
  method = multipass
}

cha002{
init:
  z=pixel  ,c=0
loop:
  z=z*z+c
  c=(1-c)*real(sqrt(c))/5+z*.75
bailout:
  |z| <= 4
default:
periodicity = 0
}

图片附件: cha001.jpg (2010-11-30 18:00, 17.34 KB) / 下载次数 5782
http://forums.netpad.net.cn/attachment.php?aid=8598&k=0e466048349006df85250f23f27662e7&t=1779870543&sid=d8ctJd

图片附件: cha002.jpg (2010-11-30 18:00, 17.19 KB) / 下载次数 5726
http://forums.netpad.net.cn/attachment.php?aid=8599&k=2c58bb5b49a81f82e1de3ede2ab71afa&t=1779870543&sid=d8ctJd

作者: 榕坚 时间: 2010-11-30 20:01

虽然简单，但却是UF循环体的好范例。练一个，顺便问一下，怎么我找遍了UF也没找着这两个分形的名称：

图片附件: cha002-3.JPG (2010-11-30 20:01, 24.38 KB) / 下载次数 5561
http://forums.netpad.net.cn/attachment.php?aid=8601&k=12e853da716d349e20d1853593afce77&t=1779870543&sid=d8ctJd

作者: changxde 时间: 2010-12-1 14:36

4# 榕坚
这两个是carr1403的变形，名字是随便摁的。

作者: changxde 时间: 2010-12-1 14:54

3#程序的画板工具

图片附件: cha0010.GIF (2010-12-7 12:45, 13.53 KB) / 下载次数 5804
http://forums.netpad.net.cn/attachment.php?aid=8616&k=63d0e9d3371c00741c455e760b217e92&t=1779870543&sid=d8ctJd

图片附件: cha0011.GIF (2010-12-7 12:45, 17.43 KB) / 下载次数 5757
http://forums.netpad.net.cn/attachment.php?aid=8617&k=07f397ff8aafebcbad46aa61e5f24d94&t=1779870543&sid=d8ctJd

图片附件: cha0020.GIF (2010-12-7 12:45, 14.18 KB) / 下载次数 5822
http://forums.netpad.net.cn/attachment.php?aid=8618&k=58097d0ed5d5e830670f16eaca8eb4a6&t=1779870543&sid=d8ctJd

图片附件: cha0021.GIF (2010-12-7 12:48, 31.52 KB) / 下载次数 5752
http://forums.netpad.net.cn/attachment.php?aid=8619&k=747db9e4146921f1519ad223217dce1e&t=1779870543&sid=d8ctJd

附件: 3#.gsp (2010-12-1 14:55, 13.83 KB) / 下载次数 5698
http://forums.netpad.net.cn/attachment.php?aid=8620&k=a71cd5d4df9680bf191a59253b12cbdc&t=1779870543&sid=d8ctJd

作者: changxde 时间: 2010-12-2 10:43

根据carr1434，1435改编的两个Uf小程序

Cha009 {
init:
  z=(1-pixel)^3, c=(-0.7456,0.212)
loop:
  z=z*z+c
  z=z+c/4
bailout:
  |z| <=4
default:
  title = "Cha009"
  periodicity = 0
  maxiter = 500
  magn = 1
  center = (1,0)
  method = multipass
}

Cha010 {
init:
  z=(1-pixel)^4, c=(-0.7456,0.212)
loop:
  z=z+c/5
  z=z*z+c

bailout:
  |z| <=4
default:
  title = "Cha010"
  periodicity = 0
  maxiter = 500
  magn = 1
  center = (1,0)
  method = multipass
}

图片附件: cha009.jpg (2010-12-7 12:50, 31.55 KB) / 下载次数 5735
http://forums.netpad.net.cn/attachment.php?aid=8639&k=f8473bbc1af5f0e82222e66285d97838&t=1779870543&sid=d8ctJd

图片附件: cha010.jpg (2010-12-7 12:50, 32.13 KB) / 下载次数 5768
http://forums.netpad.net.cn/attachment.php?aid=8640&k=48ff0a88f1ec697c60787c9bc770eb73&t=1779870543&sid=d8ctJd

作者: 榕坚 时间: 2010-12-2 11:19

试试这个UF的分形，我怎么做了几遍都与原图不大一致：
Carr1968 {
; Updated for UF2 by Erik Reckase, March 2000
         ; Modified Sylvie Gallet frm.
init:
  z=c=pixel
  int compt=0, int limit=round(real(p1)), float bailout=4
loop:
  IF (compt>=limit)
c = (-.7259,.205999)
  ENDIF
  z=z*z+c
  compt=compt+1
bailout:
  |z|<=bailout
default:
  title = "Carr 1968"
  periodicity = 0
  maxiter = 500
  magn = 1.4
  center = (0,0)
  method = multipass
  param p1
caption = "Iter Limit"
default = (100,0)
hint = "The real part of this parameter signals a change \
         in the value of c. The imaginary part is not used."
  endparam
}

图片附件: Fractal1.jpg (2010-12-2 11:19, 59.56 KB) / 下载次数 6201
http://forums.netpad.net.cn/attachment.php?aid=8641&k=e44ababe97b7cacea52a4c1ef389135f&t=1779870543&sid=d8ctJd

作者: 榕坚 时间: 2010-12-2 11:33

这是我做的结果，莫名其妙的多了一些东西：

图片附件: 捕获.JPG (2010-12-2 11:33, 31.16 KB) / 下载次数 5553
http://forums.netpad.net.cn/attachment.php?aid=8642&k=1f5cf921cef58b5d083650dc91bf5713&t=1779870543&sid=d8ctJd

作者: 柳烟 时间: 2010-12-2 16:17

榕兄可增大阈值看看，我通过改变阈值后，变正常了。搞了这么久，我发现UF中的阈值直接用于画板制作分形，有时不大对劲。

图片附件: 未命名.JPG (2010-12-2 19:09, 27.35 KB) / 下载次数 5771
http://forums.netpad.net.cn/attachment.php?aid=8647&k=b5eef792d7862ebcb14ef462c1fc63e0&t=1779870543&sid=d8ctJd

作者: changxde 时间: 2010-12-2 16:32

8# 榕坚
看看源文件，它的前100次是M集，然后在M集内做400次的J集。
下图不是Uf的方法，而是采用叠加的方法，

图片附件: carr1968.jpg (2010-12-2 16:59, 48.94 KB) / 下载次数 2426
http://forums.netpad.net.cn/attachment.php?aid=8648&k=ca2d047a08ba70f0a1b0aa52f977101c&t=1779870543&sid=d8ctJd

附件: MJ结合.gsp (2010-12-2 16:59, 50.09 KB) / 下载次数 2888
http://forums.netpad.net.cn/attachment.php?aid=8649&k=d998aae8e73108bbea811d6eaba4be94&t=1779870543&sid=d8ctJd

作者: changxde 时间: 2010-12-2 17:04

图片附件: cha009g.jpg (2010-12-7 12:39, 32.13 KB) / 下载次数 2358
http://forums.netpad.net.cn/attachment.php?aid=8650&k=7463b929e1e546c3a58d51b9156ba53f&t=1779870543&sid=d8ctJd

附件: cha009.gsp (2010-12-2 17:05, 49.43 KB) / 下载次数 3006
http://forums.netpad.net.cn/attachment.php?aid=8651&k=201a42df39a8887bd3e8c19e94101cfa&t=1779870543&sid=d8ctJd

作者: 榕坚 时间: 2010-12-2 20:15

10# 柳烟

加大#bailout的值确实可以起到作用：

图片附件: carr 1968'-1.JPG (2010-12-2 20:15, 51.33 KB) / 下载次数 2491
http://forums.netpad.net.cn/attachment.php?aid=8655&k=872503caffe564e72109a6513e625632&t=1779870543&sid=d8ctJd

作者: 柳烟 时间: 2010-12-2 20:57

13# 榕坚
如果取P1=500,则扫得混合集；如果取P1=0，则扫得纯粹J集；如果取0-500之间，则扫得混合集，在UF中放大挺有意思的，图片美。

作者: 榕坚 时间: 2010-12-2 21:07

14# 柳烟

取p1=500应为纯M集。

作者: 柳烟 时间: 2010-12-2 22:27

15# 榕坚
对。

作者: 柳烟 时间: 2010-12-2 23:14

再扫一图

图片附件: 未命名.JPG (2010-12-2 23:14, 28.38 KB) / 下载次数 2373
http://forums.netpad.net.cn/attachment.php?aid=8660&k=248475b7445d320df1ac82bc7c2042cc&t=1779870543&sid=d8ctJd

作者: changxde 时间: 2010-12-4 11:38

9# 榕坚
10# 柳烟
现象原因分析：
由于Uf中有循环跳出功能，在M集前100次迭代中，某点如果坚持不到最后，中间就退出循环，接下来的J集迭代它就不参与了，所以Uf图形外围就干净。但在Gsp中，没有循环跳出功能，对每一点都要坚持，实在坚持不住，记p=0再坚持直到循环结束。在M集前100次迭代中，某点如果坚持不到最后就是p=0，当M迭代结束，开始J集迭代时，算法改变了，原来p=0的点有些就复活，变成p=1，对et又有贡献了，所以GSP图形外围就出现了一些图形。
解决方法：
1）增大阈值使其改变算法时也不能复活；
2）改进算法使其一旦死亡就不能复活。

作者: 柳烟 时间: 2010-12-4 11:48

18# changxde
分析得很正确。

作者: changxde 时间: 2010-12-4 16:59

求助大家，
我做Carr1968在系统坐标系下做好后，扫描一切正常，但做成工具再扫描里面的J集发生了变化，请大家找找原因。
正常图形:

不正常图形:

图片附件: carr1.jpg (2010-12-4 16:59, 46.92 KB) / 下载次数 2362
http://forums.netpad.net.cn/attachment.php?aid=8713&k=6b893e1fde5ebccf208fb65f51b8698c&t=1779870543&sid=d8ctJd

图片附件: carr2.jpg (2010-12-7 12:40, 28.4 KB) / 下载次数 2393
http://forums.netpad.net.cn/attachment.php?aid=8714&k=6a63bc3f515b5bd1af0133db62554b7c&t=1779870543&sid=d8ctJd

附件: carr1968.gsp (2010-12-4 16:59, 59.46 KB) / 下载次数 2800
http://forums.netpad.net.cn/attachment.php?aid=8715&k=08a60c83dca257f8fd949476af0cbff2&t=1779870543&sid=d8ctJd

作者: 柳烟 时间: 2010-12-4 18:17

我把你制成的工具，加载后扫出的图与你的正常图一致。说明你制成的工具没有问题，请看：

加载carr1968.gsp (23.98 KB)

图片附件: 未命名.JPG (2010-12-4 18:33, 35.03 KB) / 下载次数 2377
http://forums.netpad.net.cn/attachment.php?aid=8716&k=277c33171326cc4e4e4fe2d6cddc61d5&t=1779870543&sid=d8ctJd

附件: 加载carr1968.gsp (2010-12-4 18:33, 23.98 KB) / 下载次数 3511
http://forums.netpad.net.cn/attachment.php?aid=8717&k=948e427371a699f38a477d2bfac2aefe&t=1779870543&sid=d8ctJd

作者: 柳烟 时间: 2010-12-4 20:16

changxde老师，在不改变阈值情形下，你这文件外围处理得干净。给我们很大启示。这种方法是最好的。我试用你的方法作作，如有问题，再向你请教。

作者: changxde 时间: 2010-12-4 20:41

我又试了试，工具还是不行

作者: changxde 时间: 2010-12-4 20:53

22# 柳烟
这方法可以使死灰不再复燃。

作者: 柳烟 时间: 2010-12-5 00:15

24# changxde
终于按你的方法弄成功了，你帮我们解决了困扰多时的问题，并明白了外围出现多余部分的成因。谢谢。我也该休息了。

作者: xuefeiyang 时间: 2010-12-5 16:00

[attach]8727[/attach]

图片附件: Snap195.jpg (2010-12-7 12:34, 21.77 KB) / 下载次数 2276
http://forums.netpad.net.cn/attachment.php?aid=8726&k=247ad9a0d7b52a6d1084510a2d30d753&t=1779870543&sid=d8ctJd

作者: 柳烟 时间: 2010-12-5 21:26

carr2431

用了changxde老师纯净外围环境方法，造carr2431，不错，扫两张图片。

图片附件: 未命名.JPG (2010-12-5 21:26, 24.19 KB) / 下载次数 2325
http://forums.netpad.net.cn/attachment.php?aid=8732&k=415d14830bb337d5836b9b5ac6595d32&t=1779870543&sid=d8ctJd

图片附件: 未命名2.JPG (2010-12-5 21:26, 32.11 KB) / 下载次数 2314
http://forums.netpad.net.cn/attachment.php?aid=8733&k=06313362dc6bbd9aaa572908c3a08e1f&t=1779870543&sid=d8ctJd

作者: 柳烟 时间: 2010-12-6 23:14

carr2379

图片附件: 未命名.JPG (2010-12-6 23:14, 24.67 KB) / 下载次数 2298
http://forums.netpad.net.cn/attachment.php?aid=8765&k=0dd28009de5a54e44e80a41c7b3c78a4&t=1779870543&sid=d8ctJd

图片附件: 未命名.JPG (2010-12-6 23:57, 29.7 KB) / 下载次数 2316
http://forums.netpad.net.cn/attachment.php?aid=8767&k=a4b03749b86f823a15395bb8f87981ee&t=1779870543&sid=d8ctJd

作者: xiaongxp 时间: 2010-12-6 23:27

柳老师的复分形越来越好了，干净利落！

作者: 柳烟 时间: 2010-12-6 23:33

29# xiaongxp
向老师过奖了。不知怎的，扫出图看起来还清晰，但发在坛子上，就有点乌。

作者: xiaongxp 时间: 2010-12-6 23:40

用iSee制图，保存质量设为“最佳”或“很好”，可能效果好些。你上面的图色彩这么丰富，文件才24k，当然画质不理想。

作者: 柳烟 时间: 2010-12-7 00:36

用向老师说的，处理图片，要好看多了。又学习了一招，哎，要学的太多了。

图片附件: 未命名2.jpg (2010-12-7 12:32, 37.19 KB) / 下载次数 1999
http://forums.netpad.net.cn/attachment.php?aid=8768&k=9cc35588c6ad189c11f0f1bd655eb595&t=1779870543&sid=d8ctJd

图片附件: 未命名2(1).jpg (2010-12-7 12:32, 33.01 KB) / 下载次数 2051
http://forums.netpad.net.cn/attachment.php?aid=8769&k=05fedd0cacba0885b96975df40f06c3d&t=1779870543&sid=d8ctJd

作者: xuefeiyang 时间: 2010-12-11 19:01

carr3202

图片附件: carr 3202.jpg (2010-12-16 14:04, 23.73 KB) / 下载次数 2030
http://forums.netpad.net.cn/attachment.php?aid=8855&k=146a009282bfc503108e02f78e8c84ad&t=1779870543&sid=d8ctJd

作者: 榕坚 时间: 2010-12-14 16:06

33# xuefeiyang

这个carr 3202我做的结果与你一样，但与UF的结果有差别。估计要改一下算法：

图片附件: Carr3202(YAXIS)'.JPG (2010-12-14 16:06, 67.5 KB) / 下载次数 2164
http://forums.netpad.net.cn/attachment.php?aid=8911&k=217cf22f0a48f13aad83d6aeef311d34&t=1779870543&sid=d8ctJd

图片附件: Fractal1.jpg (2010-12-14 16:06, 53.03 KB) / 下载次数 2133
http://forums.netpad.net.cn/attachment.php?aid=8912&k=756ae0580dd5a8c52a2b4451c5e0ad6a&t=1779870543&sid=d8ctJd

作者: changxde 时间: 2010-12-15 15:18

在UF中修改L1的值为300就得到与上面类似的图形

图片附件: carr3202''.jpg (2010-12-15 15:18, 72.64 KB) / 下载次数 2245
http://forums.netpad.net.cn/attachment.php?aid=8929&k=5d3eb2db610a1772f25e0c30855425f3&t=1779870543&sid=d8ctJd

作者: changxde 时间: 2010-12-15 15:34

把源程序修改一如下
Carr3202'(YAXIS) {
; Updated for UF2 by Erik Reckase, February 2000
               ; Modified Sylvie Gallet frm.
init:
  d=log(1/pixel-conj(cos(1/pixel)))
  c=(pixel^3-conj(.1/pixel-flip(.01/pixel)))/7.5
  z=abs(pixel^12)-flip(.01/pixel)-conj(.1/pixel)+.3
  zmod = (-.7456,-.132) + flip(|d/80*pixel|)
  cmod = |c| + zmod
  c1=0.2*z^1.2,c2=0.3*conj(z)
  c3=0.45*flip(z),c4=0.675*flip(z)
  l1=round(real(p1)),l2=round(imag(p1))
  l3=round(real(p2)),l4=round(imag(p2))
  float bailout=16,int iter=0
loop:
  IF (iter==l1)
cmod = |c1| + zmod, z = 0
  ELSEIF (iter==l2)
cmod = |c2| + zmod, z = 0
  ELSEIF (iter==l3)
cmod = |c3| + zmod, z = 0
  ELSEIF (iter==l4)
cmod = |c4| + zmod, z = 0
  ENDIF
  z = z^2 + cmod
  iter=iter+1
bailout:
  |z|<=bailout
default:
  title = "Carr 3202'"
  periodicity = 0
  maxiter = 700
  magn = 1.4
  center = (0,0)
  method = multipass
  param p1
caption = "Iter Limits 1,2"
default = (200,300)
hint = "0< Re(P1) < Im(P1) < Re(P2) < Im(P2) <maxiter \
         for best results"
  endparam
  param p2
caption = "Iter Limits 3,4"
default = (400,550)
hint = "0< Re(P1) < Im(P1) < Re(P2) < Im(P2) <maxiter \
         for best results"
  endparam
}
加载图像和原图没有差别，然后放大图像就发生了变化，缩小是下图

图片附件: carr3202.jpg (2010-12-15 15:34, 71.76 KB) / 下载次数 2091
http://forums.netpad.net.cn/attachment.php?aid=8930&k=9502002cb7be03dd7f86ef3870959dec&t=1779870543&sid=d8ctJd

作者: 柳烟 时间: 2010-12-16 22:09

我干了几天carr3202色一片，无图象，今天扫出中间一点，毫无美感可言，大家帮忙看看。我这文件，我查了好几遍，我按程序的先后顺序，已将各数据计算好了，容易阅读，麻烦大家看看是何原因，多谢了。问题文件中的文档1，我已将各数据计算好，还没缩放与迭代，文档2是我造好的，扫出图不对劲。

问题文件.gsp (66.04 KB)

图片附件: 未命名.jpg (2010-12-17 11:57, 12.14 KB) / 下载次数 2013
http://forums.netpad.net.cn/attachment.php?aid=8942&k=cbf8112a1c61c980450164bcd696ce5c&t=1779870543&sid=d8ctJd

附件: 问题文件.gsp (2010-12-16 22:09, 66.04 KB) / 下载次数 3449
http://forums.netpad.net.cn/attachment.php?aid=8943&k=dd456c3ea18a78404273aa444132d2bd&t=1779870543&sid=d8ctJd

作者: 榕坚 时间: 2010-12-17 08:56

37# 柳烟

修改了一下：

图片附件: 捕获.JPG (2010-12-17 08:56, 35.85 KB) / 下载次数 1918
http://forums.netpad.net.cn/attachment.php?aid=8945&k=b71a2d2fd371cdef755898ddf71de699&t=1779870543&sid=d8ctJd

附件: 问题文件.gsp (2010-12-17 08:56, 66.75 KB) / 下载次数 2878
http://forums.netpad.net.cn/attachment.php?aid=8946&k=9a243b50dd7b04d20d38bbe615c4faac&t=1779870543&sid=d8ctJd

作者: 柳烟 时间: 2010-12-17 11:02

谢过榕坚兄帮忙，我刚才看了你修改的文件，再查了UF中的帮助，原来我把/---/错当成了模的符号，原来是表示实部与虚部的平方和。

作者: 榕坚 时间: 2010-12-17 11:09

39# 柳烟
中间还有一个大红的圆环不知道原因出在哪。

作者: 柳烟 时间: 2010-12-17 18:31

不知何原因，我将|----|的式子全部改过来后，迭代竟然除了一个圆环及环中间的纯色，什么都没有。

作者: changxde 时间: 2010-12-19 16:41

图片附件: cos.jpg (2010-12-19 16:41, 44.04 KB) / 下载次数 1801
http://forums.netpad.net.cn/attachment.php?aid=8974&k=c1ea68d47fac077273431a0a0f073657&t=1779870543&sid=d8ctJd

作者: 榕坚 时间: 2010-12-19 17:06

42# changxde

真漂亮！

作者: changxde 时间: 2010-12-19 17:17

放大一下

图片附件: cos1.jpg (2010-12-19 17:17, 46.11 KB) / 下载次数 1766
http://forums.netpad.net.cn/attachment.php?aid=8975&k=d942cb058d9c989ec4aeddd0c6f10240&t=1779870543&sid=d8ctJd

图片附件: cos2.jpg (2010-12-19 17:29, 41.17 KB) / 下载次数 1784
http://forums.netpad.net.cn/attachment.php?aid=8976&k=09565e9059a3fcab84a15e680ff72f35&t=1779870543&sid=d8ctJd

作者: 柳烟 时间: 2010-12-21 19:50

40# 榕坚
原因找着了，原来是工具问题。我换了一个工具z^k，就对了。得到了与你们一样的结果。

作者: 榕坚 时间: 2010-12-21 20:41

45# 柳烟

是哪一个工具造成的？计算函数的吗？

作者: 榕坚 时间: 2010-12-21 21:01

图片附件: Newton-MSet-Mix1-5.JPG (2010-12-21 21:01, 96.3 KB) / 下载次数 1951
http://forums.netpad.net.cn/attachment.php?aid=9007&k=6aef821cdb0b2a9e6c3bdb56ddfd201d&t=1779870543&sid=d8ctJd

作者: 柳烟 时间: 2010-12-21 23:44

46# 榕坚
向老师工具箱中的z^K，象前面说的分形中莫名其妙的横线，还有前面我一直造不出的分形。我另外造作工具后，结果对了。向老师来短信说，他已更新了那工具，谢谢向老师及朋友们。

作者: 柳烟 时间: 2010-12-21 23:45

纳闷：carr3500:

图片附件: 未命名.jpg (2010-12-21 23:45, 149.82 KB) / 下载次数 2004
http://forums.netpad.net.cn/attachment.php?aid=9016&k=44f8c9204a57836c9c4d6fb44f49073e&t=1779870543&sid=d8ctJd

作者: 榕坚 时间: 2010-12-22 07:38

49# 柳烟

比我第一次做的时候好多了：

图片附件: Carr3500(YAXIS).JPG (2010-12-22 07:38, 48.36 KB) / 下载次数 1725
http://forums.netpad.net.cn/attachment.php?aid=9017&k=66e2dea741cad3569bc60ecfe0a64224&t=1779870543&sid=d8ctJd

作者: 柳烟 时间: 2010-12-22 09:54

50# 榕坚
感觉到运算的工作量有点让人喘不过气，出问题了，查起来相当麻烦，不如再推倒重来还省事，哎！

作者: 榕坚 时间: 2010-12-22 10:10

是啊，做这一系列我都是这样的。

作者: changxde 时间: 2010-12-22 20:51

图片附件: cos3.jpg (2010-12-22 20:51, 35.86 KB) / 下载次数 1259
http://forums.netpad.net.cn/attachment.php?aid=9036&k=1ca2fada6baada5b4966bc6aff1620b1&t=1779870543&sid=d8ctJd

作者: changxde 时间: 2010-12-22 21:05

图片附件: Newton-MSet-Mix1.jpg (2010-12-22 21:05, 80.04 KB) / 下载次数 1318
http://forums.netpad.net.cn/attachment.php?aid=9037&k=51afca2238018cfc3c0c49c369bd7afb&t=1779870543&sid=d8ctJd

作者: changxde 时间: 2010-12-22 21:07

图片附件: cos4.jpg (2010-12-22 21:07, 36.43 KB) / 下载次数 1214
http://forums.netpad.net.cn/attachment.php?aid=9038&k=67a3236f175b9b4a1e721148a2a3bbf6&t=1779870543&sid=d8ctJd

作者: 榕坚 时间: 2010-12-25 14:36

54# changxde

这个颜色看起来太舒服了。

作者: changxde 时间: 2011-1-8 17:00

图片附件: screenshot2.JPG (2011-1-8 17:00, 12.53 KB) / 下载次数 1228
http://forums.netpad.net.cn/attachment.php?aid=9522&k=d62f742b36440d83d1f31a1084c53bd7&t=1779870543&sid=d8ctJd

作者: changxde 时间: 2011-1-10 16:56

图片附件: screenshot3.JPG (2011-1-10 16:56, 22.27 KB) / 下载次数 1254
http://forums.netpad.net.cn/attachment.php?aid=9570&k=f1f98ad34bf202c88148664de3a19fb1&t=1779870543&sid=d8ctJd

附件: Carr1963.gsp (2011-1-10 16:56, 13.43 KB) / 下载次数 1929
http://forums.netpad.net.cn/attachment.php?aid=9571&k=782624fc7bf340fac13369c9308e55e1&t=1779870543&sid=d8ctJd

作者: changxde 时间: 2011-1-10 20:35

carr1964
仔细观察还是和UF不一样的。

图片附件: screenshot4.JPG (2011-1-10 20:35, 22.82 KB) / 下载次数 1457
http://forums.netpad.net.cn/attachment.php?aid=9573&k=47a82754869d15083527b7f30be35723&t=1779870543&sid=d8ctJd

作者: 榕坚 时间: 2011-1-11 08:04

如果不让迭代起死回生就应该一样了：

图片附件: Carr1964.JPG (2011-1-11 08:04, 35.2 KB) / 下载次数 1206
http://forums.netpad.net.cn/attachment.php?aid=9576&k=2a6d7b5e548196c5ed19d590e50e708f&t=1779870543&sid=d8ctJd

作者: changxde 时间: 2011-1-12 14:38

还是要考虑终止循环。

作者: changxde 时间: 2011-1-13 10:25

UF中的carr1958缩小到全景后，再放大缩小，可得到下面几种不同的全景图，原因何在。不知画板如何。

图片附件: carr19583.JPG (2011-1-13 10:16, 24.48 KB) / 下载次数 1584
http://forums.netpad.net.cn/attachment.php?aid=9609&k=87b355fb52ced2ae93b6d95476ba7e79&t=1779870543&sid=d8ctJd

图片附件: carr19581.JPG (2011-1-13 10:16, 35.66 KB) / 下载次数 1567
http://forums.netpad.net.cn/attachment.php?aid=9610&k=084a3e828bb3d840dd7d4adef808f2dc&t=1779870543&sid=d8ctJd

图片附件: carr19582.JPG (2011-1-13 10:16, 15.27 KB) / 下载次数 1589
http://forums.netpad.net.cn/attachment.php?aid=9611&k=81d016f810e4f15db2131e3088dcaacd&t=1779870543&sid=d8ctJd

作者: 榕坚 时间: 2011-1-13 10:37

62# changxde

不对啊，在我这里1958得到的全图始终是这一幅图：

图片附件: Fractal1.jpg (2011-1-13 10:37, 44.22 KB) / 下载次数 1432
http://forums.netpad.net.cn/attachment.php?aid=9612&k=cee638eaf7c81ddfb0aa613122904bfb&t=1779870543&sid=d8ctJd

作者: changxde 时间: 2011-1-13 10:46

又一副不同的结果（还是1958）。

图片附件: Fractal1.jpg (2011-1-13 10:46, 63.51 KB) / 下载次数 1706
http://forums.netpad.net.cn/attachment.php?aid=9613&k=8228b67ae436ac95b756ced4d50aa701&t=1779870543&sid=d8ctJd

作者: changxde 时间: 2011-1-13 10:48

1958的结构是对称的

作者: 榕坚 时间: 2011-1-13 10:54

奇怪了，从UF的缩略图看好象是对称的。可是我打开以后怎么不一样呢？

图片附件: 捕获.JPG (2011-1-13 10:54, 4.78 KB) / 下载次数 1440
http://forums.netpad.net.cn/attachment.php?aid=9615&k=ee197e46337fdbda4fc0d01922f98472&t=1779870543&sid=d8ctJd

作者: changxde 时间: 2011-1-13 11:00

与64#同...

图片附件: carr19585.JPG (2011-1-13 11:00, 63.6 KB) / 下载次数 1799
http://forums.netpad.net.cn/attachment.php?aid=9616&k=131cbdcf31f4fd6a876fd6ecf252723b&t=1779870543&sid=d8ctJd

作者: changxde 时间: 2011-1-13 11:04

UF程序是
Carr1958 (XYAXIS) { ; Updated for UF2 by Erik Reckase, March 2000
                  ; Modified Sylvie Gallet frm. REVISED 4/12/96
                  ; iter <= limit ==> Mandelbrot set of  z^7+c
                  ; iter >  limit ==> Newton's method applied to  z^6-1 = 0
init:
  op = 1/pixel, z = 1/(pixel-sin(op)-tan(op))
  c = 1/(pixel-conj(op)-flip(op))
  float b1=64, float b2=.000001, zn=pixel
  float limit = real(p1), float iter = 1
  bool test0 = TRUE, bool test3 = FALSE
loop:
  test0 = !(test0)
  IF (iter==limit)
z = zn
  ENDIF
  n = z^6-1
  IF (iter<limit)
z = z*n + z + c + (-.7456,.2)
IF (!test3)
   test3 = (|z|>b1)
ENDIF
IF (test3 && test0)
   z = 0
ENDIF
float zfinal = |z|
float bfinal = b1
  ELSE
IF (z==0)
   d = 0
ELSE
   d = 4*z^5
   z = z - n/d
ENDIF
float zfinal = b2
float bfinal = |n|
  ENDIF
  iter = iter+2
bailout:
  zfinal <= bfinal
default:
  title = "Carr 1958"
  periodicity = 0
  maxiter = 500
  magn = .8
  center = (0,0)
  method = multipass
  param p1
caption = "Iteration Limit"
default = (100,0)
hint = "The real part of this parameter signifies the point \
         at which the iterating formula changes.  The imag \
         part is not used, but remains for Fractint compat."
  endparam
}
检查一下为什么会是这样

作者: 柳烟 时间: 2011-1-13 13:41

我这里的全景图，始终是62楼中间的那付。

作者: 榕坚 时间: 2011-1-13 14:46

明白为什么会这样了：柳老师在outside中选的是BASIC模式，changxde老师选的是NONE模式，而我选的是SMOOTH模式。因此三个人得到不同的结果。在做其它分形时选这三个模式结构应该是一样的只是着色效果不同，而在这个分形中结构为什么会发生变化呢？

作者: 榕坚 时间: 2011-1-13 20:06

几何画板不但放大受到限制而且缩小的倍数受限更利害，carr 1958为例按代码做出的原图与缩小到极限的图形如下，继续缩小就没有反应了：

图片附件: Carr1958(XYAXIS)-1.JPG (2011-1-13 20:06, 62.68 KB) / 下载次数 3613
http://forums.netpad.net.cn/attachment.php?aid=9626&k=89e54b23830da609455b43d1ada4948e&t=1779870543&sid=d8ctJd

图片附件: carr 1958(XYAXIS)-2.JPG (2011-1-13 20:06, 87.36 KB) / 下载次数 3621
http://forums.netpad.net.cn/attachment.php?aid=9627&k=e5471462f35e8280154a786f0b20ba71&t=1779870543&sid=d8ctJd

作者: 榕坚 时间: 2011-1-14 07:37

想了一个办法让几何画板也有完整图：

图片附件: carr 1958(XYAXIS)-3.JPG (2011-1-14 07:37, 36.62 KB) / 下载次数 3226
http://forums.netpad.net.cn/attachment.php?aid=9632&k=2e86888521f3bfb7c1db235dd8ca8cd2&t=1779870543&sid=d8ctJd

作者: inRm 时间: 2011-1-14 08:34

以上三幅也是GSP做的？神奇。

作者: 榕坚 时间: 2011-1-14 09:03

73# inRm

是GSP完成的。

作者: 柳烟 时间: 2011-1-14 09:04

72# 榕坚
颜色调得太好了。

作者: changxde 时间: 2011-1-14 20:07

72# 榕坚
能把其他几个图案都做出来么？

作者: 榕坚 时间: 2011-1-14 21:50

76# changxde

还没，我在想那几种着色效果的具体原理如果弄明白了可能就知道该怎么做了。也可能是UF的计算原理造成的那就没办法了。

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