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) / 下载次数 3356
http://inrm3d.cn/attachment.php?aid=15908&k=3c8092b5f8ae674f3a2a367e3ab7f39c&t=1732399773&sid=78hZo7

作者: 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) / 下载次数 3266
http://inrm3d.cn/attachment.php?aid=15912&k=44f171c4bbae517e2d6978d7de2cc04b&t=1732399773&sid=78hZo7

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

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

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

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

图片附件: 未命名.jpg (2012-1-2 00:02, 57.55 KB) / 下载次数 3303
http://inrm3d.cn/attachment.php?aid=15913&k=ad939c4e483c440fc19f488deee4bd1e&t=1732399773&sid=78hZo7

附件: M集TruchetTiles(BOTH)特效.gsp (2012-1-2 12:41, 28.41 KB) / 下载次数 4123
http://inrm3d.cn/attachment.php?aid=15924&k=af582c4d38a7e63929209ffd12d617fd&t=1732399773&sid=78hZo7

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

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

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

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

图片附件: 未命名.jpg (2012-1-2 14:01, 50.58 KB) / 下载次数 3263
http://inrm3d.cn/attachment.php?aid=15923&k=39cc9c5224ac556a4a74d8eb3d828daf&t=1732399773&sid=78hZo7

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

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

图片附件: Fractal1.jpg (2012-1-2 13:33, 64.65 KB) / 下载次数 2849
http://inrm3d.cn/attachment.php?aid=15925&k=58f21fe3609e9e6d5c34316d4b38ef55&t=1732399773&sid=78hZo7

图片附件: Fractal2.jpg (2012-1-2 13:57, 33.07 KB) / 下载次数 2581
http://inrm3d.cn/attachment.php?aid=15926&k=15cbe311c132275f43b9bbc05f019677&t=1732399773&sid=78hZo7

作者: 柳烟 时间: 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) / 下载次数 2919
http://inrm3d.cn/attachment.php?aid=15927&k=1efcd458d5ae914633cc43ed5c218c07&t=1732399773&sid=78hZo7

图片附件: 未命名.jpg (2012-3-12 20:16, 185.57 KB) / 下载次数 2920
http://inrm3d.cn/attachment.php?aid=15928&k=ae52b4900a927df352bc2423e82f11d8&t=1732399773&sid=78hZo7

附件: 国际象棋棋盘扫法(更新).gsp (2012-3-12 20:16, 16.33 KB) / 下载次数 4102
http://inrm3d.cn/attachment.php?aid=15929&k=8d0e2ab287ef53c3cdf257eeef21f67b&t=1732399773&sid=78hZo7

作者: 柳烟 时间: 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) / 下载次数 1853
http://inrm3d.cn/attachment.php?aid=15930&k=788c13ef5086e09127230d150488ae21&t=1732399773&sid=78hZo7

附件: M集CheckerBoard特效.gsp (2012-1-2 20:18, 23.38 KB) / 下载次数 3136
http://inrm3d.cn/attachment.php?aid=15931&k=ce8d64c11dc3773ab70760d6b96ffcc5&t=1732399773&sid=78hZo7

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

图片附件: FamilyGradient(BOTH)-1.JPG (2012-1-2 21:02, 41.08 KB) / 下载次数 1566
http://inrm3d.cn/attachment.php?aid=15932&k=f3abca10a74c9b565f2f44ad6f0c645e&t=1732399773&sid=78hZo7

图片附件: FamilyGradient(BOTH)-3.JPG (2012-1-2 21:22, 30.24 KB) / 下载次数 1553
http://inrm3d.cn/attachment.php?aid=15933&k=c40f8905c5cb698752b70eac78385877&t=1732399773&sid=78hZo7

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

图片附件: 未命名.jpg (2012-1-4 08:30, 46.82 KB) / 下载次数 1916
http://inrm3d.cn/attachment.php?aid=15934&k=91950724b70fb75d072c0bc1669f9c1a&t=1732399773&sid=78hZo7

图片附件: 未命名.jpg (2012-1-4 08:31, 25.32 KB) / 下载次数 1861
http://inrm3d.cn/attachment.php?aid=15935&k=d5d2fc4e319781e9b2d4e770e40eb783&t=1732399773&sid=78hZo7

图片附件: 未命名.jpg (2012-1-4 08:33, 87.76 KB) / 下载次数 2044
http://inrm3d.cn/attachment.php?aid=15937&k=c223e36175dd4ead3782e14357c1363b&t=1732399773&sid=78hZo7

作者: 柳烟 时间: 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) / 下载次数 1896
http://inrm3d.cn/attachment.php?aid=15940&k=ce901390181ead9a9e3912614eed0ba6&t=1732399773&sid=78hZo7

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

迎新年，贴上一个窗花

图片附件: 捕获.JPG (2012-1-3 12:57, 28 KB) / 下载次数 1574
http://inrm3d.cn/attachment.php?aid=15942&k=68ab36e7d5add170ae7f9c1612196942&t=1732399773&sid=78hZo7

作者: 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) / 下载次数 1514
http://inrm3d.cn/attachment.php?aid=15967&k=cb5c71d221b9c17d9b07785c21634d5b&t=1732399773&sid=78hZo7

作者: 柳烟 时间: 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) / 下载次数 1795
http://inrm3d.cn/attachment.php?aid=15970&k=b8d8f80332a06c48568c6a6098aeaaee&t=1732399773&sid=78hZo7

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

weave特效.gsp (24.91 KB)

图片附件: 未命名.jpg (2012-1-6 17:59, 75.22 KB) / 下载次数 1808
http://inrm3d.cn/attachment.php?aid=15971&k=dd4b3045c057e23b8242b4764eb3a1c6&t=1732399773&sid=78hZo7

附件: weave特效.gsp (2012-1-4 23:08, 24.91 KB) / 下载次数 2758
http://inrm3d.cn/attachment.php?aid=15974&k=c7df8c1c710ba607008947b72c985faf&t=1732399773&sid=78hZo7

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

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

图片附件: 未命名.jpg (2012-1-6 18:00, 63.2 KB) / 下载次数 1873
http://inrm3d.cn/attachment.php?aid=15975&k=f979c322ee351598206f439c3e23de15&t=1732399773&sid=78hZo7

作者: 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) / 下载次数 1803
http://inrm3d.cn/attachment.php?aid=15997&k=f5de09ed034fdf1c97062f85f00f4dfa&t=1732399773&sid=78hZo7

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

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

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

图片附件: 未命名.jpg (2012-1-6 17:57, 20.5 KB) / 下载次数 1814
http://inrm3d.cn/attachment.php?aid=15998&k=e0563859253118b192239ab464b0d102&t=1732399773&sid=78hZo7

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

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

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

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

图片附件: 捕获.JPG (2012-1-6 15:54, 28.78 KB) / 下载次数 1556
http://inrm3d.cn/attachment.php?aid=15999&k=c8bd62323ab353f0238ccf9484872ad6&t=1732399773&sid=78hZo7

作者: 柳烟 时间: 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) / 下载次数 1717
http://inrm3d.cn/attachment.php?aid=16000&k=c82c146a4f2b3a3d8968529181e632a2&t=1732399773&sid=78hZo7

作者: 柳烟 时间: 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) / 下载次数 1934
http://inrm3d.cn/attachment.php?aid=16003&k=107591211b36bc123fea79acd829b05d&t=1732399773&sid=78hZo7

图片附件: 未命名.jpg (2012-1-6 23:35, 187.21 KB) / 下载次数 1954
http://inrm3d.cn/attachment.php?aid=16004&k=c9565dbe9bc22e9d5591184d47291789&t=1732399773&sid=78hZo7

图片附件: 未命名.jpg (2012-1-7 00:26, 76.53 KB) / 下载次数 1963
http://inrm3d.cn/attachment.php?aid=16006&k=d832556f85e43a97e55caa5d4c0e1878&t=1732399773&sid=78hZo7

附件: ManyJulia.gsp (2012-1-7 00:26, 19.39 KB) / 下载次数 2765
http://inrm3d.cn/attachment.php?aid=16007&k=e85162abe45bdaaca487b880e65fe462&t=1732399773&sid=78hZo7

作者: 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) / 下载次数 2538
http://inrm3d.cn/attachment.php?aid=16013&k=a429d14adf6ff3b002087b6773153b83&t=1732399773&sid=78hZo7

图片附件: 未命名.jpg (2012-1-7 19:52, 178.76 KB) / 下载次数 2518
http://inrm3d.cn/attachment.php?aid=16014&k=fc30506bc462ce15edbb700ffd628130&t=1732399773&sid=78hZo7

作者: 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) / 下载次数 2221
http://inrm3d.cn/attachment.php?aid=16025&k=217c0f31647ee20ae8a5f2ae19128b2c&t=1732399773&sid=78hZo7

作者: 柳烟 时间: 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) / 下载次数 2535
http://inrm3d.cn/attachment.php?aid=16035&k=32b9fdc3c427645c976f16d0328fe8d7&t=1732399773&sid=78hZo7

图片附件: Fractal2.jpg (2012-1-9 00:10, 149.01 KB) / 下载次数 2514
http://inrm3d.cn/attachment.php?aid=16036&k=f4fee9d3f323c52b8756c0b1f19758ec&t=1732399773&sid=78hZo7

作者: 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) / 下载次数 1484
http://inrm3d.cn/attachment.php?aid=16045&k=0b5943426a33fd2ef95a719aa76048a3&t=1732399773&sid=78hZo7

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

作者: 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) / 下载次数 1436
http://inrm3d.cn/attachment.php?aid=16050&k=9117b7749cc67195fe8304e8a678996b&t=1732399773&sid=78hZo7

作者: 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) / 下载次数 1454
http://inrm3d.cn/attachment.php?aid=16051&k=2658d2895c728a8122d2b248f242cc98&t=1732399773&sid=78hZo7

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

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

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

图片附件: 未命名.jpg (2012-3-1 23:20, 31.04 KB) / 下载次数 2791
http://inrm3d.cn/attachment.php?aid=16064&k=dc2228731c9077285d3b6d0cd419eb6d&t=1732399773&sid=78hZo7

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

图片附件: 未命名.jpg (2012-3-2 16:25, 31.85 KB) / 下载次数 2775
http://inrm3d.cn/attachment.php?aid=16117&k=35eb98301bb511519c76a5186e391adb&t=1732399773&sid=78hZo7

图片附件: 未命名1.jpg (2012-3-2 16:25, 41.71 KB) / 下载次数 2781
http://inrm3d.cn/attachment.php?aid=16118&k=7a82cef41922846e28cf0e208b07cfd9&t=1732399773&sid=78hZo7

图片附件: 未命名.jpg (2012-3-2 16:26, 39.41 KB) / 下载次数 2784
http://inrm3d.cn/attachment.php?aid=16119&k=027b8565f7c88e409627380c9bbba84d&t=1732399773&sid=78hZo7

作者: 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) / 下载次数 3129
http://inrm3d.cn/attachment.php?aid=16265&k=d856830e29ba907d96df815c4a443b92&t=1732399773&sid=78hZo7

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

图片附件: 未命名.jpg (2012-3-2 16:30, 34.78 KB) / 下载次数 2749
http://inrm3d.cn/attachment.php?aid=16267&k=e236ab921163e0e66ef2dc2904121d70&t=1732399773&sid=78hZo7

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

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

图片附件: 1.jpg (2012-3-2 16:31, 27.72 KB) / 下载次数 2772
http://inrm3d.cn/attachment.php?aid=16455&k=a71866e55fbd376dcf3c308084b13183&t=1732399773&sid=78hZo7

附件: Attractor Julia修正文件.gsp (2012-3-1 21:23, 17.87 KB) / 下载次数 3693
http://inrm3d.cn/attachment.php?aid=16456&k=710318fe8ff013b1a6939845b0eec630&t=1732399773&sid=78hZo7

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

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

图片附件: 未命名.jpg (2012-3-2 16:31, 35.42 KB) / 下载次数 2741
http://inrm3d.cn/attachment.php?aid=16457&k=274af87cb5b9e0eaeff37d8ce138fce0&t=1732399773&sid=78hZo7

作者: 柳烟 时间: 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) / 下载次数 2000
http://inrm3d.cn/attachment.php?aid=16470&k=7dc793a8b87b4551efdac91322e93d7e&t=1732399773&sid=78hZo7

图片附件: 未命名.JPG (2012-3-8 12:52, 15.41 KB) / 下载次数 1954
http://inrm3d.cn/attachment.php?aid=16477&k=622c70f041062458a65b24a7b25aa05c&t=1732399773&sid=78hZo7

附件: Barnsley 2 (Julia)1.gsp (2012-3-3 21:21, 18.35 KB) / 下载次数 2885
http://inrm3d.cn/attachment.php?aid=16478&k=2911e5a57dd54748a6cd598604bcb253&t=1732399773&sid=78hZo7

图片附件: 未命名.jpg (2012-3-5 22:06, 126.94 KB) / 下载次数 2291
http://inrm3d.cn/attachment.php?aid=16507&k=5bc25a6fc9794004ff175a5c0f59ecdc&t=1732399773&sid=78hZo7

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

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

图片附件: 未命名.jpg (2012-3-3 22:13, 60.69 KB) / 下载次数 2244
http://inrm3d.cn/attachment.php?aid=16480&k=d51d8047e6be026219dd49d025d6d65b&t=1732399773&sid=78hZo7

作者: 柳烟 时间: 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) / 下载次数 2239
http://inrm3d.cn/attachment.php?aid=16481&k=12b9f28a457a4594227b949cd6d8d569&t=1732399773&sid=78hZo7

图片附件: 未命名.JPG (2012-3-8 12:54, 20.03 KB) / 下载次数 2014
http://inrm3d.cn/attachment.php?aid=16482&k=2fa82204cbef71234693b9cb89388232&t=1732399773&sid=78hZo7

附件: Julifamily.gsp (2012-3-4 23:33, 26.73 KB) / 下载次数 3118
http://inrm3d.cn/attachment.php?aid=16483&k=3e46144b8d6935ff54e1193405c0dd51&t=1732399773&sid=78hZo7

作者: 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) / 下载次数 3138
http://inrm3d.cn/attachment.php?aid=16501&k=2733b063ab077ddf13e9f7e2d407d347&t=1732399773&sid=78hZo7

图片附件: 4.jpg (2012-3-12 12:55, 35.78 KB) / 下载次数 2004
http://inrm3d.cn/attachment.php?aid=16502&k=6c3cd80a6d224dd6afc565fff63680b9&t=1732399773&sid=78hZo7

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

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

放大一局部一小块J集:

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

图片附件: 1.jpg (2012-3-12 12:56, 44.06 KB) / 下载次数 1983
http://inrm3d.cn/attachment.php?aid=16503&k=74b331cfd4788e89f203292c81a0098a&t=1732399773&sid=78hZo7

图片附件: 2.jpg (2012-3-12 12:56, 35.02 KB) / 下载次数 1994
http://inrm3d.cn/attachment.php?aid=16504&k=17e9148394f950018259da081ff15ef7&t=1732399773&sid=78hZo7

图片附件: 3.jpg (2012-3-12 12:56, 37.11 KB) / 下载次数 1990
http://inrm3d.cn/attachment.php?aid=16505&k=c5893d25d4a0c82847e50f24a502bfaa&t=1732399773&sid=78hZo7

附件: 柳烟J集墙作品画板论坛.gsp (2012-3-5 21:51, 25.82 KB) / 下载次数 2867
http://inrm3d.cn/attachment.php?aid=16506&k=614b5b9934923b816eab4279314151f4&t=1732399773&sid=78hZo7

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

耍几个J集:
`1.jpg

图片附件: `1.jpg (2012-3-12 12:57, 30.45 KB) / 下载次数 1975
http://inrm3d.cn/attachment.php?aid=16513&k=6102c2219cd855a52d164ac4247355f2&t=1732399773&sid=78hZo7

图片附件: 5.jpg (2012-3-12 12:57, 38.12 KB) / 下载次数 1999
http://inrm3d.cn/attachment.php?aid=16514&k=c5d35981cc9db215d9fa22a21704fad4&t=1732399773&sid=78hZo7

图片附件: 未命名.jpg (2012-3-12 12:57, 25.26 KB) / 下载次数 1975
http://inrm3d.cn/attachment.php?aid=16522&k=561b984c4555af90c3b6525620f8a988&t=1732399773&sid=78hZo7

作者: 柳烟 时间: 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) / 下载次数 2255
http://inrm3d.cn/attachment.php?aid=16521&k=71cb923014f4d9c551d38bd9722592eb&t=1732399773&sid=78hZo7

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

79# 柳烟

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

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

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

再扫另一个味道的图:

图片附件: 未命名.JPG (2012-3-8 12:46, 25.87 KB) / 下载次数 1611
http://inrm3d.cn/attachment.php?aid=16530&k=17bf5ac65aec2b8db7c35d7b7bc8452b&t=1732399773&sid=78hZo7

图片附件: 未命名.JPG (2012-3-8 12:47, 34.06 KB) / 下载次数 1634
http://inrm3d.cn/attachment.php?aid=16531&k=ab5b6eb7f49546124c7b511324eb430b&t=1732399773&sid=78hZo7

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

图片附件: 未命名.JPG (2012-3-8 12:47, 31.7 KB) / 下载次数 1638
http://inrm3d.cn/attachment.php?aid=16536&k=666e3f13412000d8c463514ad2391369&t=1732399773&sid=78hZo7

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

RMExpiLin.gsp (44.32 KB)

图片附件: 未命名.jpg (2012-3-12 12:50, 15.05 KB) / 下载次数 1634
http://inrm3d.cn/attachment.php?aid=16541&k=fd68c7453920935ccec78ba3c4e1a47a&t=1732399773&sid=78hZo7

图片附件: 未命名.jpg (2012-3-12 12:50, 24.74 KB) / 下载次数 1604
http://inrm3d.cn/attachment.php?aid=16557&k=d68a986b37fe37c32d7941e708eab15a&t=1732399773&sid=78hZo7

附件: RMExpiLin.gsp (2012-3-11 00:24, 44.32 KB) / 下载次数 2931
http://inrm3d.cn/attachment.php?aid=16558&k=e0d9b273747c56c9583265ccdaeeb579&t=1732399773&sid=78hZo7

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

RMExpiNMp-O.gsp (45.37 KB)

图片附件: 未命名.jpg (2012-3-11 12:28, 112.48 KB) / 下载次数 1737
http://inrm3d.cn/attachment.php?aid=16561&k=ebe20ca00309443cbc6c379263a028c0&t=1732399773&sid=78hZo7

附件: RMExpiNMp-O.gsp (2012-3-11 12:34, 45.37 KB) / 下载次数 2667
http://inrm3d.cn/attachment.php?aid=16563&k=356828b2b5f7b05a81bb0e4f971ae9f7&t=1732399773&sid=78hZo7

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

柳烟作品:J集:

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

图片附件: 未命名.jpg (2012-3-12 12:51, 25.82 KB) / 下载次数 1500
http://inrm3d.cn/attachment.php?aid=16566&k=3d3b911efd8615f4bfb3cb33d9c2aa7e&t=1732399773&sid=78hZo7

图片附件: 未命名.jpg (2012-3-12 12:52, 28.77 KB) / 下载次数 1526
http://inrm3d.cn/attachment.php?aid=16567&k=557c61a40be406b16a95ead243c6edb3&t=1732399773&sid=78hZo7

图片附件: 未命名.jpg (2012-3-12 12:52, 32.46 KB) / 下载次数 1500
http://inrm3d.cn/attachment.php?aid=16569&k=bd401d25e04e7e70960b7f77cf568bdb&t=1732399773&sid=78hZo7

图片附件: 未命名.jpg (2012-3-12 12:53, 29.94 KB) / 下载次数 1511
http://inrm3d.cn/attachment.php?aid=16581&k=a35dc56df7b59fc4f81cb00a94a7e375&t=1732399773&sid=78hZo7

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

柳烟作品:N集:

图片附件: 未命名.jpg (2012-3-12 12:54, 25.06 KB) / 下载次数 1528
http://inrm3d.cn/attachment.php?aid=16587&k=2dc872f7baf9e6bd607b4409c300dddc&t=1732399773&sid=78hZo7

作者: 柳烟 时间: 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) / 下载次数 1719
http://inrm3d.cn/attachment.php?aid=16635&k=4c2d5c3a66c4987fe0e26b62bc25223a&t=1732399773&sid=78hZo7

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

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

图片附件: 未命名.JPG (2012-3-17 13:58, 34.09 KB) / 下载次数 1501
http://inrm3d.cn/attachment.php?aid=16636&k=6cbb8b501fa7c14d86e5440bc3aeff8b&t=1732399773&sid=78hZo7

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

图片附件: 未命名.jpg (2012-3-19 22:41, 34.02 KB) / 下载次数 1544
http://inrm3d.cn/attachment.php?aid=16660&k=9f808ec3eccd334a9b70ac9dba3f6206&t=1732399773&sid=78hZo7

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

这是柳某演绎Uf的作品.

图片附件: 未命名.jpg (2012-3-20 22:52, 122.28 KB) / 下载次数 1721
http://inrm3d.cn/attachment.php?aid=16672&k=72625003076c7fec6706a1097b7073b0&t=1732399773&sid=78hZo7

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

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

图片附件: 未命名.jpg (2012-3-21 21:18, 120.11 KB) / 下载次数 3879
http://inrm3d.cn/attachment.php?aid=16682&k=398eb0265d30a0a87030426bd1783f7d&t=1732399773&sid=78hZo7

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

柳烟作品:

图片附件: 未命名.jpg (2012-3-22 20:58, 142.19 KB) / 下载次数 3963
http://inrm3d.cn/attachment.php?aid=16695&k=f965d2ab51627be8f69eb89b90767503&t=1732399773&sid=78hZo7

作者: 榕坚 时间: 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) / 下载次数 4075
http://inrm3d.cn/attachment.php?aid=16696&k=5f8e9ba3e73f26d0127676487930cf08&t=1732399773&sid=78hZo7

作者: 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) / 下载次数 4825
http://inrm3d.cn/attachment.php?aid=16709&k=8a0399c54f3ec2d38c3937b30fce4f73&t=1732399773&sid=78hZo7

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

98# 柳烟

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

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

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

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