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

标题: UF与复分形特效 [打印本页]

作者: 柳烟 时间: 2012-1-1 20:19 标题: UF与复分形特效

M集TruchetTiles(BOTH)特效：
未命名.jpg

附上UF中该特效作色的程序代码：
TruchetTiles(BOTH) {
; By Samuel Monnier, 31.10.00
init:
  z = 0
  zc = 0
  zz = 0
  cr = 0
  float d1 = 0
  float d2 = 0
loop:

final:
  z = #z/@size*exp(1i*pi/180*@rot)
  zc = round(z)
  zz = z - zc

  cr = 2*(@seed/(zc+124) - round(@seed/(zc+124)))

  if real(cr) > 0
if @mode == 0
   d1 = abs(cabs(zz+(.5,.5))-.5)
   d2 = abs(cabs(zz-(.5,.5))-.5)
else
   d1 = abs(real(zz) - imag(zz) -.5)
   d2 = abs(real(zz) - imag(zz) +.5)
endif
if d2 < d1
   d1 = d2
endif
  else
if @mode == 0
   d1 = abs(cabs(zz+(.5,-.5))-.5)
   d2 = abs(cabs(zz-(.5,-.5))-.5)
else
   d1 = abs(real(zz) + imag(zz) -.5)
   d2 = abs(real(zz) + imag(zz) +.5)
endif
if d2 < d1
   d1 = d2
endif
  endif

  #index = d1^@thick

default:
  title = "Truchet Tiles"
  helpfile = "sam-help/variouscol2.htm"
  helptopic = "truchet"

  param mode
caption = "Mode"
default = 0
enum = "Roundy" "Squarry"
  endparam

  param size
caption = "Pattern Size"
default = 1.0
  endparam

  param rot
caption = "Pattern Rotation"
default = 0.0
  endparam

  param thick
caption = "Thickness"
default = 0.1
  endparam

  param seed
caption = "Seed"
default = 1234567890
  endparam
}

图片附件: 未命名.jpg (2012-1-2 00:00, 71.43 KB) / 下载次数 2970
http://www.inrm3d.cn/attachment.php?aid=15908&k=b5188d477dd9779e71ea3e06308a5f12&t=1715845475&sid=vIKCiE

作者: xiaongxp 时间: 2012-1-1 22:05

1# 柳烟
此图太妙了！是有机一体的还是二次扫描图？

作者: 柳烟 时间: 2012-1-1 22:10

2# xiaongxp
一次扫成的。

作者: 柳烟 时间: 2012-1-1 22:34

改变参数，再扫：

图片附件: 未命名.jpg (2012-1-2 00:01, 63.13 KB) / 下载次数 2891
http://www.inrm3d.cn/attachment.php?aid=15912&k=c52039d698abad5c9a88e3ca8c7a7c74&t=1715845475&sid=vIKCiE

作者: 柳烟 时间: 2012-1-1 23:43

UF中的这个特效带两个开关项，要在一个画板文件中，整合这两个开关项，不难，无非多几个判断而已。分形着色，用了本坛榕坚兄发在该坛的一种特效中的着色法。一楼图是其中一个开关项之一，今再扫其姐妹篇：
未命名.jpg

一楼的画板源文件，待我整理好后，再发在此。大家一同努力，将分形推向高潮。
附上1楼的源文件，以方便交流：

M集TruchetTiles(BOTH)特效.gsp (28.41 KB)

图片附件: 未命名.jpg (2012-1-2 00:02, 57.55 KB) / 下载次数 2911
http://www.inrm3d.cn/attachment.php?aid=15913&k=453d6bbfde817776d97ab0cb61a0dad0&t=1715845475&sid=vIKCiE

附件: M集TruchetTiles(BOTH)特效.gsp (2012-1-2 12:41, 28.41 KB) / 下载次数 3766
http://www.inrm3d.cn/attachment.php?aid=15924&k=3311b94f8b5614e660952a8186548990&t=1715845475&sid=vIKCiE

作者: dyk 时间: 2012-1-2 11:29

太漂亮了，赞。背景着色的好方法。

作者: 柳烟 时间: 2012-1-2 12:34

将程序代码中的某些句子动一动，扫得一图：

图片附件: 未命名.jpg (2012-1-2 14:01, 50.58 KB) / 下载次数 2903
http://www.inrm3d.cn/attachment.php?aid=15923&k=4dda0e11192d481ba70f5bac8f40d9bc&t=1715845475&sid=vIKCiE

作者: 榕坚 时间: 2012-1-2 13:33

这UF中的宝贝可真多，怎么也学不完：

图片附件: Fractal1.jpg (2012-1-2 13:33, 64.65 KB) / 下载次数 2492
http://www.inrm3d.cn/attachment.php?aid=15925&k=0120b19720df1785678e79bda50e4188&t=1715845475&sid=vIKCiE

图片附件: Fractal2.jpg (2012-1-2 13:57, 33.07 KB) / 下载次数 2267
http://www.inrm3d.cn/attachment.php?aid=15926&k=4c7d32e1d3c99b566d8267fb05948286&t=1715845475&sid=vIKCiE

作者: 柳烟 时间: 2012-1-2 14:47

8# 榕坚
说得不错，UF确实是学习复分形的绝好教材，是分形艺术的“敦煌莫高窟”。

作者: 柳烟 时间: 2012-1-2 19:26

这个复分形没整成功，程序解读起来费神，弄来不对头，等势圈老是不能消亡了。摆出来，算是征解。不过，倒学会了扫描法扫出国际象棋棋盘：

国际象棋棋盘：

国际象棋棋盘扫法(更新).gsp (16.33 KB)

图片附件: Fractal2.jpg (2012-1-2 20:50, 30.99 KB) / 下载次数 2594
http://www.inrm3d.cn/attachment.php?aid=15927&k=fd8643a81551ca7611447ef5d33477e0&t=1715845475&sid=vIKCiE

图片附件: 未命名.jpg (2012-3-12 20:16, 185.57 KB) / 下载次数 2611
http://www.inrm3d.cn/attachment.php?aid=15928&k=dfdd3b9229727d2036dd0cdd0861aea2&t=1715845475&sid=vIKCiE

附件: 国际象棋棋盘扫法(更新).gsp (2012-3-12 20:16, 16.33 KB) / 下载次数 3706
http://www.inrm3d.cn/attachment.php?aid=15929&k=a567061d2b01c3d5421deea4aa22f6fc&t=1715845475&sid=vIKCiE

作者: 柳烟 时间: 2012-1-2 19:58

我刚才胡整海整，歪打正着，居然作出了，怪哉！我今天就干这个特效，总是带来整去不对劲，感谢上天。
未命名.jpg

M集CheckerBoard特效.gsp (23.38 KB)
CheckerBoard {
; Marcelo Anelli
  ; Caution!!:
  ; this method does not work with fast Mandel and fast Julia.
  init:
float p = 0.0
float w = 320 / #width * 0.03 * sqrt(@p1)
int tx = 0
int ty = 0
float rx = real( #screenpixel )
float ry = imag( #screenpixel )
tx = floor( rx * w )
if @ratio
   ty = floor( ry * w * #width / #height )
else
   ty = floor( ry * w )
endif
int t = ( tx + ty ) % 2
if (t == 0)
   p = @pmin
else
   p = @pmax
endif
  final:
#index = p
default:
  title = "CheckerBoard"
  param p1
caption = "Width"
default = 1.0
hint = "Not depending of magnification."
  endparam
  param ratio
default = false
caption = "Aspect ratio"
hint = "Enable to match the aspect ratio of the square and the image"
  endparam
  param pmin
caption = "White Square Color"
min = 0.0
max = 1.0
default = 0.1
  endparam
  param pmax
caption = "Black Square Color"
min = 0.0
max = 1.0
default = 0.7
  endparam
}

图片附件: 未命名.jpg (2012-1-2 20:48, 33.56 KB) / 下载次数 1720
http://www.inrm3d.cn/attachment.php?aid=15930&k=46e99baee3ea358ef8ef13d01d11642e&t=1715845475&sid=vIKCiE

附件: M集CheckerBoard特效.gsp (2012-1-2 20:18, 23.38 KB) / 下载次数 2940
http://www.inrm3d.cn/attachment.php?aid=15931&k=448b7fae462e246a0e40706caafd4b64&t=1715845475&sid=vIKCiE

作者: 榕坚 时间: 2012-1-2 21:02

图片附件: FamilyGradient(BOTH)-1.JPG (2012-1-2 21:02, 41.08 KB) / 下载次数 1453
http://www.inrm3d.cn/attachment.php?aid=15932&k=3b95eca814e4015c1c40598de33f2dc5&t=1715845475&sid=vIKCiE

图片附件: FamilyGradient(BOTH)-3.JPG (2012-1-2 21:22, 30.24 KB) / 下载次数 1441
http://www.inrm3d.cn/attachment.php?aid=15933&k=b5c4dbe42eca18c3e14b44ae44a3df4d&t=1715845475&sid=vIKCiE

作者: 柳烟 时间: 2012-1-2 23:00

图片附件: 未命名.jpg (2012-1-4 08:30, 46.82 KB) / 下载次数 1769
http://www.inrm3d.cn/attachment.php?aid=15934&k=692eaab913142dd650ceabda731026d9&t=1715845475&sid=vIKCiE

图片附件: 未命名.jpg (2012-1-4 08:31, 25.32 KB) / 下载次数 1711
http://www.inrm3d.cn/attachment.php?aid=15935&k=66e4cb76ba446662cab64eedcd62d04b&t=1715845475&sid=vIKCiE

图片附件: 未命名.jpg (2012-1-4 08:33, 87.76 KB) / 下载次数 1889
http://www.inrm3d.cn/attachment.php?aid=15937&k=3e5c9ced6e42297e652db00e0aa9b201&t=1715845475&sid=vIKCiE

作者: 柳烟 时间: 2012-1-3 00:09

12# 榕坚
不同的着色模式，引发不同的视觉冲击。再扫一图：
未命名.jpg

特效代码如下：
FamilyGradient(BOTH) {
; By Samuel Monnier, 20.3.00
; Designed for the "Familly formulas"
; (Julifamilly and Pcfamilly in my file
; and Pinskyfamilly in kpk.ufm.
init:
  complex z = pixel
  complex zc = 0
loop:

final:
  z = z - @tilem
  z = z*exp(1i*pi/180*@rottile)
  zc = round(z/@sizesq)*@sizesq
  z = z*exp(-1i*pi/180*@rottile)
  zc = zc*exp(-1i*pi/180*@rottile)
  z = z + @tilem
  zc = zc + @tilem
  z = 2*(z - zc)/@sizesq
  z = z*exp(1i*pi/180*@rottile)
  if @mode == 0
#index = cabs(z)
  elseif @mode == 1
float x = abs(real(z))
float y = abs(imag(z))
if x > y
   #index = x
else
   #index = y
endif
  endif
default:
  title = "Gradient for <<Familly>>"
  helpfile = "sam-help/variouscol.htm"
  helptopic = "gradforfam"

  param mode
caption = "Gradient Mode"
default = 1
enum = "Circles" "Squares"
  endparam

  param tilem
caption = "Tiling displacement"
default = (0,0)
  endparam

  param sizesq
caption = "Size of the tiles"
default = 0.6
  endparam

  param rottile
caption = "Tiling rotation"
default = 0.0
  endparam
}

图片附件: 未命名.jpg (2012-1-4 08:35, 35.36 KB) / 下载次数 1757
http://www.inrm3d.cn/attachment.php?aid=15940&k=67be5472d529fb9414f3f8efb3f8f1fd&t=1715845475&sid=vIKCiE

作者: 榕坚 时间: 2012-1-3 12:57

迎新年，贴上一个窗花

图片附件: 捕获.JPG (2012-1-3 12:57, 28 KB) / 下载次数 1462
http://www.inrm3d.cn/attachment.php?aid=15942&k=8a808b501e93711395ab7420c1299e43&t=1715845475&sid=vIKCiE

作者: yandongtai 时间: 2012-1-3 13:58

各位对几何画板与分形的热爱与痴迷真让人佩服。学习了！

作者: 柳烟 时间: 2012-1-3 20:59

UF中的有些特效太迷人，但破译起来太费神，判断多，令人眼花了乱，只好放弃了。

作者: 榕坚 时间: 2012-1-3 21:02

UF中肯定也有边界构造的特效，能把它找出来吗？

作者: 榕坚 时间: 2012-1-3 21:07

17# 柳烟
有些特效确实已经超出了几何画板的能力范围，象一些3D的特效连UF本身运行起来也非常慢，只能选作几何画板能为的。

作者: 柳烟 时间: 2012-1-3 22:17

19# 榕坚
完全正确。

作者: dyk 时间: 2012-1-4 10:19

老师们的刻苦钻研精神是画板分形向高峰攀登的源动力，致敬！向你们学习、学习、再学习。

作者: dyk 时间: 2012-1-4 10:32

我对UF可以说是一点不懂，只是简单地基本操作，不知道E文，采用乱试看效果，更谈不上读代码。从试验的角度看，分形主要有公式，内部着色和外部着色，公式决定了分形的结构，通过内外着色来实现不同的效果。是这样的吗？

作者: 柳烟 时间: 2012-1-4 12:54

22# dyk
应该是这样的。我们也对Ｅ文一窍不通，只是按代码程序，大概加估计进行的。

作者: dyk 时间: 2012-1-4 13:30

23# 柳烟

谢谢柳老师，我非常赞赏你的毅力和智慧。

作者: 榕坚 时间: 2012-1-4 16:18

画板演绎UF的乐趣：先是被分形奇妙的结构所吸引，利用UF已完成的范例充分欣赏一幅幅漂亮的图形，接下来就有了用几何画板创作的冲动。能一次成功那当然是最理想的，但对于一些复杂些的分形往往要承受一次次失败的打击，放弃……继续……放弃……继续，感受直到某一次扫出正确图形时的喜悦，当中当然离不开失败后一次次仔细琢磨原代码（原代码其实就是计算机的高级语言或类似，只要学过BASIC的一般是没有问题的，当中当然有一些英文注释，看得懂最好，看不懂也没多大关系，因为数学语言是世界通用的）。

图片附件: 2.JPG (2012-1-4 16:21, 84.24 KB) / 下载次数 1396
http://www.inrm3d.cn/attachment.php?aid=15967&k=97a04315ddc34db313ee3d012432210c&t=1715845475&sid=vIKCiE

作者: 柳烟 时间: 2012-1-4 21:13

25# 榕坚
说得太好了！今天搞了几个特效，以失败告终。今晚搞另一个，失败了好几次，终于迎来希望的曙光。
mt-weave {
; Mark Townsend, 15 November 2000
final:
  float this = 0
  float that = 0
  float x = 4096 + real(#z) * @size
  float y = 4096 + imag(#z) * @size
  int xi = floor(x)
  int yi = floor(y)
  float xr = x - xi
  float yr = y - yi
  float s = (1 - @width) / 2
  if @shading == 0
this = yr
that = xr
  else
this = xr
that = yr
  endif
  if (xi + yi) % 2 == 0
if (yr > s)&& (yr < 1 - s)
   #index = this
else
   if (xr > s)&& (xr < 1 - s)
      #index = that
   else
      #solid = true
   endif
endif
  else
if (xr > s) && (xr < 1 - s)
   #index = that
else
   if (yr > s) && (yr < 1 - s)
      #index = this
   else
      #solid = true
   endif
endif
  endif
default:
  title = "Weave"
  param shading
caption = "Shading"
enum = "Thread" "Ribbon"
  endparam
  param width
caption = "Width"
default = 0.9
max = 1.0
  endparam
  param size
caption = "Scale"
default = 1.0
  endparam
}

图片附件: 未命名.jpg (2012-1-4 22:44, 58.12 KB) / 下载次数 1628
http://www.inrm3d.cn/attachment.php?aid=15970&k=32cf8cd292a6e6f4461b3de2b3f72457&t=1715845475&sid=vIKCiE

作者: 柳烟 时间: 2012-1-4 22:22

weave特效.gsp (24.91 KB)

图片附件: 未命名.jpg (2012-1-6 17:59, 75.22 KB) / 下载次数 1644
http://www.inrm3d.cn/attachment.php?aid=15971&k=6ff1fc7e58598a2d92c95ffb73540901&t=1715845475&sid=vIKCiE

附件: weave特效.gsp (2012-1-4 23:08, 24.91 KB) / 下载次数 2574
http://www.inrm3d.cn/attachment.php?aid=15974&k=301d98f5999d8702cebe83053e33fcfd&t=1715845475&sid=vIKCiE

作者: 柳烟 时间: 2012-1-5 00:09

将就UF胚胎，擅自变动程序，扫UF没有的：

图片附件: 未命名.jpg (2012-1-6 18:00, 63.2 KB) / 下载次数 1702
http://www.inrm3d.cn/attachment.php?aid=15975&k=bad69af349f5b8eaae6462fa3c7d447a&t=1715845475&sid=vIKCiE

作者: dyk 时间: 2012-1-5 09:48

失败乃成功之母，柳老师在失败中成功了，祝贺。

作者: 柳烟 时间: 2012-1-6 13:40

bezier-curve {
; Kerry Mitchell 08apr00
;
; Colors by the orbit's closest approach to a user-defined Bezier curve.
; The curve is determined by specifying beginning and ending anchor points,
; through which the curve passes, and 2 control points, which influence the
; shape of the curve.
;
init:
  float x0=real(@z0)
  float y0=imag(@z0)
  float x1=real(@z1)
  float y1=imag(@z1)
  float x2=real(@z2)
  float y2=imag(@z2)
  float x3=real(@z3)
  float y3=imag(@z3)
  float cx=3*(x1-x0)
  float bx=3*(x2-x1)-cx
  float ax=x3-x0-cx-bx
  float cy=3*(y1-y0)
  float by=3*(y2-y1)-cy
  float ay=y3-y0-cy-by
  float t=0.0
  float r=0.0
  float x=0.0
  float y=0.0
  float u=0.0
  float v=0.0
  int iter=0
  float rmin=1.0e20
  int itermin=0
  zmin=(0.0,0.0)
loop:
  iter=iter+1
  u=real(#z)
  v=imag(#z)
;
; The curve is parameterized with x(t) and y(t).  Step through several t
; values to find the nearest approach of the orbit to the curve.
;
  t=0.0
  while(t<=1.0)
x=((ax*t+bx)*t+cx)*t+x0
y=((ay*t+by)*t+cy)*t+y0
r=(x-u)*(x-u)+(y-v)*(y-v)
if(r<rmin)
   rmin=r
   itermin=iter
   zmin=#z
endif
t=t+@dt
  endwhile
final:
  if(@colorby==1)          ; iteration @ min
#index=0.01*itermin
  elseif(@colorby==2)       ; angle @ min
t=atan2(zmin)
t=t/pi
if(t<0.0)
   t=t+2.0
endif
#index=0.5*t
  elseif(@colorby==3)       ; draw section
u=real(#pixel)
v=imag(#pixel)
t=0.0
rmin=1e20
while(t<=1.0)
   x=((ax*t+bx)*t+cx)*t+x0
   y=((ay*t+by)*t+cy)*t+y0
   r=(x-u)*(x-u)+(y-v)*(y-v)
   if(r<rmin)
      rmin=r
   endif
   t=t+@dt
endwhile
#index=rmin^@nexp
  else                      ; minimum distance
#index=rmin^@nexp
  endif
default:
  title="Bezier Curve"
  helpfile="lkm-help\lkm-bezier.html"
  param z0
caption="1st anchor point"
default=(1.0,0.0)
hint="Curve starts at this point."
  endparam
  param z1
caption="1st control point"
default=(1.0,1.0)
hint="Influences the shape of the curve."
  endparam
  param z2
caption="2nd control point"
default=(0.0,0.0)
hint="Influences the shape of the curve."
  endparam
  param z3
caption="2nd anchor point"
default=(0.0,1.0)
hint="Curve ends at this point."
  endparam
  param dt
caption="step size"
default=0.1
hint="Decrease for smoother line, increase \
   to see dots.  Should be between 0 & 1."
min=0.0
max=1.0
  endparam
  param nexp
caption="power"
default=0.1
min=0.0
hint="Decrease to make thinner lines. Use \
   with 'minimum distance' coloring."
  endparam
  param colorby
caption="color by"
default=0
enum="minimum distance" "iteration @ min" \
   "angle @ min" "show curve"
  endparam
}
ＵＦ中效果图：

我干了两天，劳而无功。

图片附件: Fractal2.jpg (2012-1-6 18:01, 143.78 KB) / 下载次数 1650
http://www.inrm3d.cn/attachment.php?aid=15997&k=bb5d06e112047e0992193dd89fe394f2&t=1715845475&sid=vIKCiE

作者: 柳烟 时间: 2012-1-6 13:53

这是我扫出的较为接近的病态文件：

好象是等势圈将金链子盖住了。

图片附件: 未命名.jpg (2012-1-6 17:57, 20.5 KB) / 下载次数 1679
http://www.inrm3d.cn/attachment.php?aid=15998&k=6584ff0e1899602c9952ff25957b98c7&t=1715845475&sid=vIKCiE

作者: 榕坚 时间: 2012-1-6 15:41

你不要把et做为着色参数试试。

作者: 榕坚 时间: 2012-1-6 15:54

我也弄了一个不象样的东西：

图片附件: 捕获.JPG (2012-1-6 15:54, 28.78 KB) / 下载次数 1456
http://www.inrm3d.cn/attachment.php?aid=15999&k=ec565ccc679f33cae33554fc05adfcf7&t=1715845475&sid=vIKCiE

作者: 柳烟 时间: 2012-1-6 16:28

我是先按照着色程序中的t到t+dt，由判断t<=1,迭代9次，得到rmin后，将所有的z都试过，替换u\v,不成。又将9次得到的rmin,再将z^2+c或按p缩放后的点替换程序中的u\v,再二次迭代，得到终点rmin的终点，再得到index,没有引进et,仍有等势圈。
观上图，比我扫出的，更接近原图，请问迭代规则是如何弄的？

作者: 榕坚 时间: 2012-1-6 18:21

34# 柳烟

应该是只要将Z(没经过缩放）赋值给u,v。再迭代。

作者: 榕坚 时间: 2012-1-6 18:39

不过这个特效很占内存的，只迭代30次（加小循环10次）就非常慢了。别指望它扫放大图。

作者: 榕坚 时间: 2012-1-6 18:45

这个特效更绝：

图片附件: Fractal1.jpg (2012-1-6 18:47, 36.81 KB) / 下载次数 1595
http://www.inrm3d.cn/attachment.php?aid=16000&k=63b64f3d514cc4fb818af83e583a8a3a&t=1715845475&sid=vIKCiE

作者: 柳烟 时间: 2012-1-6 20:07

36# 榕坚
应该非常慢，我弄出的魇品，扫起来都慢如蜗牛。

作者: 柳烟 时间: 2012-1-6 20:08

37# 榕坚
里面有好几个绝的，还有个在等势圈上扫出希尔伯特曲线，我想只能欣赏，用画板实现起来，太遥远了，甚至不可能。

作者: 柳烟 时间: 2012-1-6 23:26

UF中的算法真是诡谲莫测，令人叹服。同时也很纳闷，当时跟着UF学复分形，竟然让这个不太难的奇妙分形成了漏网之鱼！
未命名.jpg

改变参数，居然扫出了M集，神奇！

ManyJulia.gsp (19.39 KB)

图片附件: 未命名.jpg (2012-1-6 23:26, 105.16 KB) / 下载次数 1774
http://www.inrm3d.cn/attachment.php?aid=16003&k=a33c6ee19412640bd0aadc73d48fa6b2&t=1715845475&sid=vIKCiE

图片附件: 未命名.jpg (2012-1-6 23:35, 187.21 KB) / 下载次数 1814
http://www.inrm3d.cn/attachment.php?aid=16004&k=5b0414146d8ad7f5fc470711c26c4f56&t=1715845475&sid=vIKCiE

图片附件: 未命名.jpg (2012-1-7 00:26, 76.53 KB) / 下载次数 1830
http://www.inrm3d.cn/attachment.php?aid=16006&k=03b436af54574e4a805bd2c6bbdc50f1&t=1715845475&sid=vIKCiE

附件: ManyJulia.gsp (2012-1-7 00:26, 19.39 KB) / 下载次数 2579
http://www.inrm3d.cn/attachment.php?aid=16007&k=5065b00b8aad3cc1ddb2810dd9ea7bce&t=1715845475&sid=vIKCiE

作者: xiaongxp 时间: 2012-1-6 23:48

40# 柳烟
妙！真妙！一图扫出多个位置的J集，简直妙不可言！要是以M集为背景，把J集与M集的定位关系反应出来，那简直就妙极了。

作者: 柳烟 时间: 2012-1-7 19:26

旁边的ManyNova :

图片附件: 未命名.jpg (2012-1-7 19:26, 148.86 KB) / 下载次数 2363
http://www.inrm3d.cn/attachment.php?aid=16013&k=6882a71d8fedbf25fcf64e351713ce38&t=1715845475&sid=vIKCiE

图片附件: 未命名.jpg (2012-1-7 19:52, 178.76 KB) / 下载次数 2353
http://www.inrm3d.cn/attachment.php?aid=16014&k=86e592d3d6584553ebe68d369bfafba0&t=1715845475&sid=vIKCiE

作者: dyk 时间: 2012-1-8 09:29

42# 柳烟
柳老师对UF吃的透，对画板玩的精。这种执着的精神值得发扬学习。

作者: 榕坚 时间: 2012-1-8 11:22

41# xiaongxp

可以用几何画板的图层合并来完成：

图片附件: 捕获2.JPG (2012-1-8 11:22, 39.55 KB) / 下载次数 2075
http://www.inrm3d.cn/attachment.php?aid=16025&k=fdb615deeaf0cf21f25cbf611cec9617&t=1715845475&sid=vIKCiE

作者: 柳烟 时间: 2012-1-8 11:29

43# dyk
各有千秋，比如dem法，就是我向此坛板友学的，“精”字实不敢当。

作者: 柳烟 时间: 2012-1-8 11:31

44# 榕坚
很是漂亮，学习。

作者: xiaongxp 时间: 2012-1-8 16:05

44# 榕坚
经典之作！要的就是M集的这种“字典”效果。

作者: 柳烟 时间: 2012-1-9 00:10

我整出的效果图是：
未命名.jpg

而UF的效果图是：

等势圈的差别原因何在，半天不明所以。
String-ArtColoring(BOTH) {
; By Samuel Monnier, 30.9.01
init:
  z = 0
  float r = 0
  float rtest = 0
  float arg = 0
  float arg2 = 0
  float d = 0

final:
  z = #z
  r = cabs(z)
  arg = atan2(z)
  arg2 = arg + @phi*pi/180
  if @shape == 0
rtest = (@w1/(@w1+@w2)*cos(@order1*arg)+@w2/(@w1+@w2)*cos(@order2*arg+@phi*pi/180))*(@obound-@ibound)/2+(@obound+@ibound)/2
  elseif @shape == 1
rtest = (@w1/(@w1+@w2)*(2*(cos(@order1*arg)/2+.5)^@power-1)+@w2/(@w1+@w2)*(2*(cos(@order2*arg2)/2+.5)^@power-1))*(@obound-@ibound)/2+(@obound+@ibound)/2
  elseif @shape == 2
rtest = (@w1/(@w1+@w2)*(2*(2*abs(arg/(2*pi)*@order1-round(arg/(2*pi)*@order1)))^@power-1)+@w2/(@w1+@w2)*(2*(2*abs(arg2/(2*pi)*@order2-round(arg2/(2*pi)*@order2)))^@power-1))*(@obound-@ibound)/2+(@obound+@ibound)/2
  endif
  d = abs(r-rtest)
  #index = d^@thick

default:
  title = "String-Art Coloring"
  helpfile = "sam-help/stringart.htm"
  helptopic = "strartcol"

  param shape
caption = "Folium Type"
default = 0
enum = "Old" "Cosinus" "Sawtooth"
  endparam

  param power
caption = "Folium Power"
default = 1.0
  endparam

  param order1
caption = "Folium 1 Order"
default = 5.0
  endparam

  param w1
caption = "Folium 1 Weight"
default = 1.0
  endparam

  param order2
caption = "Folium 2 Order"
default = 3.0
  endparam

  param w2
caption = "Folium 2 Weight"
default = 0.0
  endparam

  param phi
caption = "Dephasage"
default = 90.0
  endparam

  param obound
caption = "Folium Outer Bound"
default = 2.0
  endparam

  param ibound
caption = "Folium Inner Bound"
default = .4
  endparam

  param thick
caption = "Thickness"
default = 1.0
  endparam
}

图片附件: 未命名.jpg (2012-1-9 00:10, 176.97 KB) / 下载次数 2349
http://www.inrm3d.cn/attachment.php?aid=16035&k=8d71cf76d345503f81c124d6667b29ab&t=1715845475&sid=vIKCiE

图片附件: Fractal2.jpg (2012-1-9 00:10, 149.01 KB) / 下载次数 2342
http://www.inrm3d.cn/attachment.php?aid=16036&k=ff7124573bef11dc7932e89dd002b608&t=1715845475&sid=vIKCiE

作者: dyk 时间: 2012-1-9 09:32

44# 榕坚

请问榕老师画板的图层合并是同一条线着两次吗？

作者: 榕坚 时间: 2012-1-9 10:25

49# dyk

与UF一样，就是你之前的多张图片重叠的方法。底层设置不透明，重在上面的设置不同的透明程度。一条扫描线扫两幅图也可以但效果不佳。

作者: dyk 时间: 2012-1-9 11:21

50# 榕坚
谢谢榕老师。

作者: 柳烟 时间: 2012-1-9 12:23

50# 榕坚
这是分形作色的重要经验之一，学习。

作者: 柳烟 时间: 2012-1-9 13:41

按dxy老师的处理Ｍ集内外着色方法，已圆满解决#48楼问题，借此向dxy老师表示感谢。扫图一幅：

M集String-ArtColoring(BOTH)特效.gsp (24.09 KB)

图片附件: 未命名.jpg (2012-3-2 16:32, 35.66 KB) / 下载次数 1320
http://www.inrm3d.cn/attachment.php?aid=16045&k=7d12929839040bd8d9304f126692dce6&t=1715845475&sid=vIKCiE

附件: M集String-ArtColoring(BOTH)特效.gsp (2012-1-9 13:47, 24.09 KB) / 下载次数 2575
http://www.inrm3d.cn/attachment.php?aid=16046&k=d82f71a895b900b65ccb3e4b12a01112&t=1715845475&sid=vIKCiE

作者: xiaongxp 时间: 2012-1-9 14:29

53# 柳烟
漂亮

作者: 榕坚 时间: 2012-1-9 14:35

53# 柳烟

放大图的小M有变化吗？

作者: 柳烟 时间: 2012-1-9 16:00

将１楼文件，按dxy老师的内外赋色法作一下色，扫得一图：

图片附件: 未命名.jpg (2012-3-2 16:33, 35.95 KB) / 下载次数 1290
http://www.inrm3d.cn/attachment.php?aid=16050&k=93c49f60cdfd8381b961881433dfa69c&t=1715845475&sid=vIKCiE

作者: dyk 时间: 2012-1-9 16:04

56# 柳烟

太漂亮了。

作者: 柳烟 时间: 2012-1-9 16:12

55# 榕坚
放大后，好象小Ｍ肚中没有花色，可能是我的作色算法有些问题，我再查查。谢过榕兄指正。

作者: 榕坚 时间: 2012-1-9 16:17

58# 柳烟

不要查了，UF中也没有，我就是觉得奇怪。

作者: 柳烟 时间: 2012-1-9 16:23

59# 榕坚
我查了，我那图与ＵＦ一致，扫出的小Ｍ集有些有，有些没有。

图片附件: 未命名.jpg (2012-3-2 16:33, 28.72 KB) / 下载次数 1295
http://www.inrm3d.cn/attachment.php?aid=16051&k=c52a724acd619f9bd430126149b534fe&t=1715845475&sid=vIKCiE

作者: 柳烟 时间: 2012-1-9 20:34

UF中有一个与上楼类似的，放大后每个小M集肚子上的花纹一样，但代码简直搞不懂，只好望洋兴叹。

作者: 柳烟 时间: 2012-1-10 00:16

图片附件: 未命名.jpg (2012-3-1 23:20, 31.04 KB) / 下载次数 2443
http://www.inrm3d.cn/attachment.php?aid=16064&k=030b3f39615e1db7316d9f8c6c2a4b7c&t=1715845475&sid=vIKCiE

作者: 柳烟 时间: 2012-1-13 01:15

图片附件: 未命名.jpg (2012-3-2 16:25, 31.85 KB) / 下载次数 2440
http://www.inrm3d.cn/attachment.php?aid=16117&k=bcda3daf1736e059572ab2f775684e36&t=1715845475&sid=vIKCiE

图片附件: 未命名1.jpg (2012-3-2 16:25, 41.71 KB) / 下载次数 2441
http://www.inrm3d.cn/attachment.php?aid=16118&k=2e9c9632973612456ece24493395a57b&t=1715845475&sid=vIKCiE

图片附件: 未命名.jpg (2012-3-2 16:26, 39.41 KB) / 下载次数 2447
http://www.inrm3d.cn/attachment.php?aid=16119&k=3a77b2c909775add96735db144232a56&t=1715845475&sid=vIKCiE

作者: xiaongxp 时间: 2012-1-13 14:14

请问这扫描版的希尔伯特曲线的迭代公式是什么？

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

64# xiaongxp

HilbertCurve(BOTH) {
; By Samuel Monnier, 2.9.00
init:
  z = 0
  int i = 0
  int ttype = 0
  float d = -1e20
loop:

final:

  z = #z/2
  z = abs(z) - (.5,.5)
  while i < @niter
i = i + 1
if ttype == 0
   if real(z) < 0 && imag(z) < 0
      z = 2*z + (.5,.5)
      z = -conj(z)
   elseif real(z) > 0 && imag(z) < 0
      z = 2*z + (-.5,.5)
      z = 1i*z
   elseif real(z) > 0 && imag(z) > 0
      z = 2*z + (-.5,-.5)
      z = 1i*conj(z)
   elseif real(z) < 0 && imag(z) > 0
      z = 2*z + (.5,-.5)
      ttype = 1
   endif
else
   if real(z) < 0 && imag(z) < 0
      z = 2*z + (.5,.5)
      z = 1i*z
      ttype = 0
   elseif real(z) > 0 && imag(z) < 0
      z = 2*z + (-.5,.5)
      z = -1i*conj(z)
      ttype = 0
   elseif real(z) > 0 && imag(z) > 0
      z = 2*z + (-.5,-.5)
      ttype = 0
   elseif real(z) < 0 && imag(z) > 0
      z = 2*z + (.5,-.5)
      z = -conj(z)
      ttype = 0
   endif
endif
  endwhile

  if @style == 0
z = z + (1,1)
if ttype == 0
   if abs(real(z))-1 > d
      d = abs(real(z))-1
   endif
   if abs(imag(z))-1 > d
      d = abs(imag(z))-1
   endif
else
   d = imag(z)-1
endif
  elseif @style == 1 || @style == 2
if ttype == 0
   d = cabs(z+(.5,.5))-.5
else
   if @style == 1
      d = imag(z)
   else
      d = imag(z) - (real(z)^2-.25)^2*3
   endif
endif
  elseif @style == 3

if ttype == 0
   z = z + (.5,.5)
   d = abs(real(z)) + abs(imag(z)) - .5
else
   ;d = imag(z)-.5
   d = abs(real(z)) - imag(z) - .5
endif
  endif

  #index = abs(d)^@power

default:
  title = "Hilbert Curve"
  helpfile = "sam-help/hilbert.htm"

  param style
caption = "Style"
default = 0
enum = "Square" "Round I" "Round II" "Diagonal"
  endparam

  param power
caption = "Thickness"
default = .1
  endparam

  param niter
caption = "Number of Iterations"
default = 4
  endparam
}

作者: 柳烟 时间: 2012-1-13 23:09

扫描版的希尔伯特曲线与扫描版的谢尔兵斯基三角形一样，数学机理还不是很清楚。如何将迭代函数系统IFS，用扫描法整出，是一个有趣而意义深远的课题。象UF中扫描版的三圆的极限集，如何用画板实现，这些课题都难如哥德巴赫猜想，但又充满魅力。
上楼程序具体实现，详见榕坚兄大帖：分形奇图共评赏里，他修改柳烟的病态文件后的未命名文件里。

作者: 柳烟 时间: 2012-2-7 09:17

近一个月出外，没沾过电脑，如今归来，好多东西又搞忘了，又重新起摊。

图片附件: 未命名.jpg (2012-3-2 16:29, 66.86 KB) / 下载次数 2719
http://www.inrm3d.cn/attachment.php?aid=16265&k=dcc6c3bfc272f2bcc49cfe7bcbef93a9&t=1715845475&sid=vIKCiE

作者: 柳烟 时间: 2012-2-7 12:53

图片附件: 未命名.jpg (2012-3-2 16:30, 34.78 KB) / 下载次数 2407
http://www.inrm3d.cn/attachment.php?aid=16267&k=7ce666838f99c53b5eee4e8c71b40d04&t=1715845475&sid=vIKCiE

作者: 柳烟 时间: 2012-3-1 19:42

Attractor Julia修正文件.gsp (17.87 KB)

图片附件: 1.jpg (2012-3-2 16:31, 27.72 KB) / 下载次数 2415
http://www.inrm3d.cn/attachment.php?aid=16455&k=bbeafe6248b77a92e375b8fe1a4dcb6e&t=1715845475&sid=vIKCiE

附件: Attractor Julia修正文件.gsp (2012-3-1 21:23, 17.87 KB) / 下载次数 3362
http://www.inrm3d.cn/attachment.php?aid=16456&k=aeb31b1bb03edd4347d3186fca1813a6&t=1715845475&sid=vIKCiE

作者: 柳烟 时间: 2012-3-1 22:09

上楼分形,很久以前作过,今天重新作,发现原来对UF的理解不透彻,有些问题.现将就UF中的参数,把上楼文件加点反演变换,扫图:

图片附件: 未命名.jpg (2012-3-2 16:31, 35.42 KB) / 下载次数 2405
http://www.inrm3d.cn/attachment.php?aid=16457&k=dce2c8c5edcc07939aa09ca083ebd2f4&t=1715845475&sid=vIKCiE

作者: 柳烟 时间: 2012-3-3 20:02

Barnsley 2 (Julia)1.gsp (18.35 KB)
分形楼房:

将上面图片分形楼的局部再放大一张:

再扫:
[attach]16475[/attach]

图片附件: 未命名.JPG (2012-3-8 12:49, 15.79 KB) / 下载次数 1863
http://www.inrm3d.cn/attachment.php?aid=16470&k=3ba346a5954bea6557440a121bce39ee&t=1715845475&sid=vIKCiE

图片附件: 未命名.JPG (2012-3-8 12:52, 15.41 KB) / 下载次数 1819
http://www.inrm3d.cn/attachment.php?aid=16477&k=7359651a1d280a0b460203dc2dc602df&t=1715845475&sid=vIKCiE

附件: Barnsley 2 (Julia)1.gsp (2012-3-3 21:21, 18.35 KB) / 下载次数 2720
http://www.inrm3d.cn/attachment.php?aid=16478&k=50c93f7e1573bd65919b3f155fad131b&t=1715845475&sid=vIKCiE

图片附件: 未命名.jpg (2012-3-5 22:06, 126.94 KB) / 下载次数 2135
http://www.inrm3d.cn/attachment.php?aid=16507&k=b33b3eb48af034646c3b88aac87c48cb&t=1715845475&sid=vIKCiE

作者: 柳烟 时间: 2012-3-3 22:13

由上楼文件,改变参数后扫出的某部分.

图片附件: 未命名.jpg (2012-3-3 22:13, 60.69 KB) / 下载次数 2082
http://www.inrm3d.cn/attachment.php?aid=16480&k=ff1d9fc9c45f588a038464b70c7d7025&t=1715845475&sid=vIKCiE

作者: 柳烟 时间: 2012-3-4 23:29

UF真是学复变分形的绝好教材,算法千奇百怪,让人神往.以前用画板玩UF中的分形,太粗糙,这几天细玩,觉得很有意思,今天又研究以前没做过的分形,让人精神振奋,完全可以将各种算法综合一下,造出新的分形.这UF的确是取之不尽,用之不完的宝库,在我们手头缺失复分形材料的情况下,这UF是个好老师.
扫一幅J集展览墙,以前面发过的那个复分形算法不一样.在UF中改变各参数,可看到各参数的功能,比画板中玩效率高,成图快,画板扫图实在不敢恭维.

Julifamily.gsp (26.73 KB)

图片附件: 未命名11.jpg (2012-3-4 23:29, 81.24 KB) / 下载次数 2088
http://www.inrm3d.cn/attachment.php?aid=16481&k=662f5ece6dbd9c986dad71228370de92&t=1715845475&sid=vIKCiE

图片附件: 未命名.JPG (2012-3-8 12:54, 20.03 KB) / 下载次数 1876
http://www.inrm3d.cn/attachment.php?aid=16482&k=d94874321bff4531a6e76b32d1c01187&t=1715845475&sid=vIKCiE

附件: Julifamily.gsp (2012-3-4 23:33, 26.73 KB) / 下载次数 2935
http://www.inrm3d.cn/attachment.php?aid=16483&k=d014cc447b569fc46729227aa289a26c&t=1715845475&sid=vIKCiE

作者: mjj_ljh 时间: 2012-3-5 09:42

柳老师uf玩的好啊！向您学习。这是哪个系列？

作者: 柳烟 时间: 2012-3-5 15:24

74# mjj_ljh
梅老师好。此分形在sam.ufm系列中。

作者: 柳烟 时间: 2012-3-5 21:40

Barnsley 1 (Julia).gsp (18.08 KB)

附件: Barnsley 1 (Julia).gsp (2012-3-5 21:40, 18.08 KB) / 下载次数 2945
http://www.inrm3d.cn/attachment.php?aid=16501&k=46cdb31c622931ec50e9d33a2e6b1262&t=1715845475&sid=vIKCiE

图片附件: 4.jpg (2012-3-12 12:55, 35.78 KB) / 下载次数 1864
http://www.inrm3d.cn/attachment.php?aid=16502&k=d8fdfd3ae6b0a37e4bffd1b1a257eb05&t=1715845475&sid=vIKCiE

作者: 柳烟 时间: 2012-3-5 21:43

将前面的J集墙算法与上楼的J集算法糅合在一起,出生一新的J集墙:

放大一局部一小块J集:

柳烟J集墙作品画板论坛.gsp (25.82 KB)

图片附件: 1.jpg (2012-3-12 12:56, 44.06 KB) / 下载次数 1842
http://www.inrm3d.cn/attachment.php?aid=16503&k=d97088fdbe8506088471bea673445830&t=1715845475&sid=vIKCiE

图片附件: 2.jpg (2012-3-12 12:56, 35.02 KB) / 下载次数 1856
http://www.inrm3d.cn/attachment.php?aid=16504&k=c0a02edf9957409b8c9603be6c80e0ac&t=1715845475&sid=vIKCiE

图片附件: 3.jpg (2012-3-12 12:56, 37.11 KB) / 下载次数 1825
http://www.inrm3d.cn/attachment.php?aid=16505&k=b0e43cce26e2ed02f98ec3dc3a703fa5&t=1715845475&sid=vIKCiE

附件: 柳烟J集墙作品画板论坛.gsp (2012-3-5 21:51, 25.82 KB) / 下载次数 2683
http://www.inrm3d.cn/attachment.php?aid=16506&k=9a032e34a5d24e9ff1e91e4d9128d8de&t=1715845475&sid=vIKCiE

作者: 柳烟 时间: 2012-3-6 21:53

耍几个J集:
`1.jpg

图片附件: `1.jpg (2012-3-12 12:57, 30.45 KB) / 下载次数 1805
http://www.inrm3d.cn/attachment.php?aid=16513&k=2a4f48bbed8d8b834330f7257febdea0&t=1715845475&sid=vIKCiE

图片附件: 5.jpg (2012-3-12 12:57, 38.12 KB) / 下载次数 1857
http://www.inrm3d.cn/attachment.php?aid=16514&k=c832ad7a174b061af48ed4885579c888&t=1715845475&sid=vIKCiE

图片附件: 未命名.jpg (2012-3-12 12:57, 25.26 KB) / 下载次数 1817
http://www.inrm3d.cn/attachment.php?aid=16522&k=36cb3c1dde1869019a1fe137becee62e&t=1715845475&sid=vIKCiE

作者: 柳烟 时间: 2012-3-7 15:34

这个J集来自网上，下面我给出网址连接，由于网站上是E文，本人不认识洋文，所以这个J弄来不对劲，大家帮忙看看。上楼的J集也是来自此网址。
http://eldar.mathstat.uoguelph.ca/dashlock/ftax/QJG.html
这是第九张图片。

图片附件: QJG009.JPG (2012-3-7 15:34, 135.52 KB) / 下载次数 2087
http://www.inrm3d.cn/attachment.php?aid=16521&k=ec4bbed21f321aeaa8d7ab71db86dd39&t=1715845475&sid=vIKCiE

作者: 榕坚 时间: 2012-3-7 16:06

79# 柳烟

第9图的那个参数好象有误。

作者: 柳烟 时间: 2012-3-7 22:38

以前研究过的J集,今重新调色,觉得亲切.

再扫另一个味道的图:

图片附件: 未命名.JPG (2012-3-8 12:46, 25.87 KB) / 下载次数 1487
http://www.inrm3d.cn/attachment.php?aid=16530&k=e94e2c8eb1b507662699fe5572a56d22&t=1715845475&sid=vIKCiE

图片附件: 未命名.JPG (2012-3-8 12:47, 34.06 KB) / 下载次数 1511
http://www.inrm3d.cn/attachment.php?aid=16531&k=2ed06d4c589f60a3109116954e5faefa&t=1715845475&sid=vIKCiE

作者: 柳烟 时间: 2012-3-8 12:29

图片附件: 未命名.JPG (2012-3-8 12:47, 31.7 KB) / 下载次数 1513
http://www.inrm3d.cn/attachment.php?aid=16536&k=df9afd6b9b602eed5a33d83ad57390fb&t=1715845475&sid=vIKCiE

作者: 柳烟 时间: 2012-3-8 23:36

RMExpiLin.gsp (44.32 KB)

图片附件: 未命名.jpg (2012-3-12 12:50, 15.05 KB) / 下载次数 1519
http://www.inrm3d.cn/attachment.php?aid=16541&k=ca83f17267d0df4c66cdc56f024f8795&t=1715845475&sid=vIKCiE

图片附件: 未命名.jpg (2012-3-12 12:50, 24.74 KB) / 下载次数 1482
http://www.inrm3d.cn/attachment.php?aid=16557&k=0d15c0530e69eeb8ec9afff98c5faf96&t=1715845475&sid=vIKCiE

附件: RMExpiLin.gsp (2012-3-11 00:24, 44.32 KB) / 下载次数 2751
http://www.inrm3d.cn/attachment.php?aid=16558&k=a7e1dcc727c258953f6e19ebdf039324&t=1715845475&sid=vIKCiE

作者: 柳烟 时间: 2012-3-11 12:28

RMExpiNMp-O.gsp (45.37 KB)

图片附件: 未命名.jpg (2012-3-11 12:28, 112.48 KB) / 下载次数 1614
http://www.inrm3d.cn/attachment.php?aid=16561&k=f33c826185481c88af8cddf8d6ba8fa4&t=1715845475&sid=vIKCiE

附件: RMExpiNMp-O.gsp (2012-3-11 12:34, 45.37 KB) / 下载次数 2513
http://www.inrm3d.cn/attachment.php?aid=16563&k=7f0afe59a52569cec73fdb56fcb57dde&t=1715845475&sid=vIKCiE

作者: 柳烟 时间: 2012-3-11 14:52

柳烟作品:J集:

絮我调色不行,但这些J集结构颇不错的.

图片附件: 未命名.jpg (2012-3-12 12:51, 25.82 KB) / 下载次数 1403
http://www.inrm3d.cn/attachment.php?aid=16566&k=7cace5631d840b81e401ba925ede6dab&t=1715845475&sid=vIKCiE

图片附件: 未命名.jpg (2012-3-12 12:52, 28.77 KB) / 下载次数 1412
http://www.inrm3d.cn/attachment.php?aid=16567&k=2da29578eff1744c323936820bdfc77d&t=1715845475&sid=vIKCiE

图片附件: 未命名.jpg (2012-3-12 12:52, 32.46 KB) / 下载次数 1412
http://www.inrm3d.cn/attachment.php?aid=16569&k=40788b16d47d97e195570d05ae52972b&t=1715845475&sid=vIKCiE

图片附件: 未命名.jpg (2012-3-12 12:53, 29.94 KB) / 下载次数 1410
http://www.inrm3d.cn/attachment.php?aid=16581&k=6cc415a954f6e790b2f07a7e94e01a48&t=1715845475&sid=vIKCiE

作者: 柳烟 时间: 2012-3-12 00:03

柳烟作品:N集:

图片附件: 未命名.jpg (2012-3-12 12:54, 25.06 KB) / 下载次数 1417
http://www.inrm3d.cn/attachment.php?aid=16587&k=cf8c7e6c9758925d99882e0bc65ae1a3&t=1715845475&sid=vIKCiE

作者: 柳烟 时间: 2012-3-17 13:51

征解:

着色代码:
dmj-Lyapunov {
;
; This algorithm computes the Lyapunov exponent
; for Mandelbrot types. This exponent is usually
; negative for divergent orbits (outside points),
; and positive for convergent orbits (inside
; points).
;
; This might have interesting results when used
; with other fractal types, although the results
; would not be mathematically accurate.
;
; Optimizations suggested by Charles Vassallo.
;
init:
  float oldsum = 0
  float sum = 1
  float v = 0
  float il = 1/log(real(@power))
  float lp = log(log(@bailout)/2.0)
  float f = 0.0

loop:
  IF (@trackvariable == 0) ; |z|
v = cabs(#z)
  ELSEIF (@trackvariable == 1) ; real(z)
v = real(#z)
  ELSEIF (@trackvariable == 2) ; imag(z)
v = imag(#z)
  ELSEIF (@trackvariable == 3) ; real(z)/imag(z)
v = real(#z)/imag(#z)
  ELSEIF (@trackvariable == 4) ; imag(z)/real(z)
v = imag(#z)/real(#z)
  ELSEIF (@trackvariable == 5) ; arg(z)
v = atan2(#z)
  ELSEIF (@trackvariable == 6) ; 1/real(z)
v = 1.0/real(#z)
  ELSEIF (@trackvariable == 7) ; 1/imag(z)
v = 1.0/imag(#z)
  ENDIF

  oldsum = sum
;  sum = sum + log(abs(2*v)) ; sum the Lyapunov exponent (slow method)
  sum = sum * (abs(2*v)) ; sum the Lyapunov exponent

final:
  oldsum = log(oldsum)
  sum = log(sum)
  IF (@negative == 1)
sum = -sum/#numiter
oldsum = -oldsum/(#numiter-1)
  ELSEIF (@negative == 2)
sum = abs(sum/#numiter)
oldsum = abs(oldsum/(#numiter-1))
  ELSE
sum = sum/#numiter
oldsum = oldsum/(#numiter-1)
  ENDIF
  IF (@smooth)
f = il*lp - il*log(log(cabs(#z)))
#index = oldsum + (sum-oldsum) * (f+1)
  ELSE
#index = sum
  ENDIF

default:
  title = "Lyapunov"
  helpfile = "dmj-pub\dmj-pub-uf-lyapunov.htm"

  param trackvariable
caption = "Variable to Track"
default = 0
enum = "magnitude of z" "real part of z" "imaginary part of z" \
         "real / imag" "imag / real" "angle of z" "1 / real(z)" "1 / imag(z)"
hint = "Indicates which variable to measure the Lyapunov exponent for."
  endparam
  param negative
caption = "Sign"
default = 2
enum = "positive" "negative" "absolute value"
hint = "Affects the sign of the exponent.  'Negative' and 'absolute \
         value' are useful for inside coloring."
  endparam
  param power
caption = "Exponent"
default = 2.0
hint = "This should be set to match the exponent of the \
         formula you are using.  For Mandelbrot, this is 2."
  endparam
  param bailout
caption = "Bailout"
default = 1e20
min = 1
hint = "This should be set to match the bailout value in \
         the Formula tab.  Use a very high bailout!"
  endparam
  param smooth
caption = "Smooth Coloring"
default = false
hint = "If set, results will be 'smoothed' to hide iteration bands."
  endparam
}

图片附件: Fractal2.png (2012-3-17 13:51, 95.74 KB) / 下载次数 1605
http://www.inrm3d.cn/attachment.php?aid=16635&k=6c1c30e387c5a74dc322bb2ef104162d&t=1715845475&sid=vIKCiE

作者: 柳烟 时间: 2012-3-17 13:58

这是我几天前扫出的魇品,有个别的芽苞上乌龟肚腹上没有圆圈,而没达到UF中的每个芽苞上有圈圈.怪哉的是,我重新弄,反而连各别乌龟上的圆圈圈都弄不出来了,真是活见鬼.

图片附件: 未命名.JPG (2012-3-17 13:58, 34.09 KB) / 下载次数 1386
http://www.inrm3d.cn/attachment.php?aid=16636&k=9521b33982f091d3bf039a5f38e98fd5&t=1715845475&sid=vIKCiE

作者: 柳烟 时间: 2012-3-19 21:42

图片附件: 未命名.jpg (2012-3-19 22:41, 34.02 KB) / 下载次数 1428
http://www.inrm3d.cn/attachment.php?aid=16660&k=c935abce7a5589fbb9c599c34175ecd0&t=1715845475&sid=vIKCiE

作者: 柳烟 时间: 2012-3-20 22:52

这是柳某演绎Uf的作品.

图片附件: 未命名.jpg (2012-3-20 22:52, 122.28 KB) / 下载次数 1613
http://www.inrm3d.cn/attachment.php?aid=16672&k=3c73d867336732c6819808800cfd698b&t=1715845475&sid=vIKCiE

作者: 柳烟 时间: 2012-3-21 21:18

将UF中的分形玩完,我就收山了,画板这软件,毕竟是个小软件.再一个原因是,我有更好的书要读,人生活得累没意思,本人性喜逍遥与自在,祝个位画板技艺日新又日新,日日新.

图片附件: 未命名.jpg (2012-3-21 21:18, 120.11 KB) / 下载次数 3345
http://www.inrm3d.cn/attachment.php?aid=16682&k=b7850c7fc81ae31ef06cefc52923aed6&t=1715845475&sid=vIKCiE

作者: 柳烟 时间: 2012-3-22 20:58

柳烟作品:

图片附件: 未命名.jpg (2012-3-22 20:58, 142.19 KB) / 下载次数 3423
http://www.inrm3d.cn/attachment.php?aid=16695&k=e208fbf9b643365d525a0bd571969070&t=1715845475&sid=vIKCiE

作者: 榕坚 时间: 2012-3-22 21:30

92# 柳烟

这个逃逸区有趣。

作者: 榕坚 时间: 2012-3-22 22:10

87# 柳烟
我最怕的就是着色中有对数函数，几何画板有个缺陷：if n=A else n=B我们通常是用sgn来处理的，它必须A与B都有意义才能计算，即它不是真正的判断函数，这与编程软件是有区别的。

图片附件: dmj-Lyapunov-2.JPG (2012-3-22 22:10, 59.05 KB) / 下载次数 3462
http://www.inrm3d.cn/attachment.php?aid=16696&k=f2205cbc7f622468d23f8c6862543441&t=1715845475&sid=vIKCiE

作者: inRm 时间: 2012-3-22 23:04

92# 柳烟

精彩！

作者: 柳烟 时间: 2012-3-22 23:43

87# 柳烟
我最怕的就是着色中有对数函数，几何画板有个缺陷：if n=A else n=B我们通常是用sgn来处理的，它必须A与B都有意义才能计算，即它不是真正的判断函数，这与编程软件是有区别的。
榕坚发表于 2012-3-22 22:10

请榕坚兄说说画板作此的造法大略,谢谢.

作者: 榕坚 时间: 2012-3-23 08:03

96# 柳烟

这个着色没有特殊的方法啊，你做的结果应该是正确的，把颜色调淡些试试。

作者: 柳烟 时间: 2012-3-24 23:47

今天继续造87#那个特效，仍然不能确保每个牙苞上均内含圈圈。我将文件发在此，该代码不复杂，也不长，就是弄来不能让人圆满。榕兄有空帮忙看看，我老是在原地踏步，跳不出来，越弄越糊涂了。
这是我那个不成器的文件：

222.gsp (23.82 KB)

附件: 222.gsp (2012-3-24 23:47, 23.82 KB) / 下载次数 4347
http://www.inrm3d.cn/attachment.php?aid=16709&k=1be2118dc1719f3415d8aa114f9ae6d0&t=1715845475&sid=vIKCiE

作者: 榕坚 时间: 2012-3-27 19:20

98# 柳烟

这两天忙着学校的一级达标检查。看了一下你的文件，按代码的意思sum是需要进行迭代的，而f是不要迭代的，你再试试。

作者: 柳烟 时间: 2012-4-2 23:43

请问榕坚兄,代码中的#numiter相当于画板造分形的何参数?

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