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) / 下载次数 2890
http://inrm3d.cn/attachment.php?aid=15908&k=f85c26238e24211aa2f610a968027993&t=1714685954&sid=H5L7kq

作者: 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) / 下载次数 2800
http://inrm3d.cn/attachment.php?aid=15912&k=33a2382b6af0d15373190cefc51e2782&t=1714685954&sid=H5L7kq

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

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

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

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

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

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

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

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

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

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

图片附件: 未命名.jpg (2012-1-2 14:01, 50.58 KB) / 下载次数 2812
http://inrm3d.cn/attachment.php?aid=15923&k=d8310bf262424314dc37ba5b756650ee&t=1714685954&sid=H5L7kq

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

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

图片附件: Fractal1.jpg (2012-1-2 13:33, 64.65 KB) / 下载次数 2429
http://inrm3d.cn/attachment.php?aid=15925&k=aa6d4e14fd1b21be272448cbb349fbc9&t=1714685954&sid=H5L7kq

图片附件: Fractal2.jpg (2012-1-2 13:57, 33.07 KB) / 下载次数 2195
http://inrm3d.cn/attachment.php?aid=15926&k=20d7e3c15d6dfa7054b9b34397ad2086&t=1714685954&sid=H5L7kq

作者: 柳烟 时间: 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) / 下载次数 2512
http://inrm3d.cn/attachment.php?aid=15927&k=fa9b0363f27fae9de94ec25e3d5af08c&t=1714685954&sid=H5L7kq

图片附件: 未命名.jpg (2012-3-12 20:16, 185.57 KB) / 下载次数 2527
http://inrm3d.cn/attachment.php?aid=15928&k=6e5ef2820a555d0a717a561f93f40095&t=1714685954&sid=H5L7kq

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

作者: 柳烟 时间: 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) / 下载次数 1698
http://inrm3d.cn/attachment.php?aid=15930&k=c7223bdf1ff40f7b78b84ae3a475c4c9&t=1714685954&sid=H5L7kq

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

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

图片附件: FamilyGradient(BOTH)-1.JPG (2012-1-2 21:02, 41.08 KB) / 下载次数 1437
http://inrm3d.cn/attachment.php?aid=15932&k=31d9fb56aaddc245d08e57ca67e551a0&t=1714685954&sid=H5L7kq

图片附件: FamilyGradient(BOTH)-3.JPG (2012-1-2 21:22, 30.24 KB) / 下载次数 1428
http://inrm3d.cn/attachment.php?aid=15933&k=8fec7ea7922a4f043c37e44a2e202730&t=1714685954&sid=H5L7kq

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

图片附件: 未命名.jpg (2012-1-4 08:30, 46.82 KB) / 下载次数 1757
http://inrm3d.cn/attachment.php?aid=15934&k=d4f317686640cc9a81b46e62b53c9658&t=1714685954&sid=H5L7kq

图片附件: 未命名.jpg (2012-1-4 08:31, 25.32 KB) / 下载次数 1694
http://inrm3d.cn/attachment.php?aid=15935&k=490f22c8bf992c1001e46003b0cd163f&t=1714685954&sid=H5L7kq

图片附件: 未命名.jpg (2012-1-4 08:33, 87.76 KB) / 下载次数 1874
http://inrm3d.cn/attachment.php?aid=15937&k=0bdbd855933e7cce746d4adcaee66cd4&t=1714685954&sid=H5L7kq

作者: 柳烟 时间: 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) / 下载次数 1736
http://inrm3d.cn/attachment.php?aid=15940&k=f9b5cec3552919dbccebef717ef353d2&t=1714685954&sid=H5L7kq

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

迎新年，贴上一个窗花

图片附件: 捕获.JPG (2012-1-3 12:57, 28 KB) / 下载次数 1450
http://inrm3d.cn/attachment.php?aid=15942&k=d204606e16c6b4d2cf82d0eb3c89ae3f&t=1714685954&sid=H5L7kq

作者: 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) / 下载次数 1389
http://inrm3d.cn/attachment.php?aid=15967&k=9a96fe46c453ca7559d7773a2cc1e9a3&t=1714685954&sid=H5L7kq

作者: 柳烟 时间: 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) / 下载次数 1617
http://inrm3d.cn/attachment.php?aid=15970&k=ab7a3722b2453e2e796336ef8552eb04&t=1714685954&sid=H5L7kq

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

weave特效.gsp (24.91 KB)

图片附件: 未命名.jpg (2012-1-6 17:59, 75.22 KB) / 下载次数 1633
http://inrm3d.cn/attachment.php?aid=15971&k=902116a786ac2cd6e6e1bcbd11f84f19&t=1714685954&sid=H5L7kq

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

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

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

图片附件: 未命名.jpg (2012-1-6 18:00, 63.2 KB) / 下载次数 1690
http://inrm3d.cn/attachment.php?aid=15975&k=9136d7dc688fda096057148d51b49d30&t=1714685954&sid=H5L7kq

作者: 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) / 下载次数 1642
http://inrm3d.cn/attachment.php?aid=15997&k=14533e8685ccccd2f743581943c77b01&t=1714685954&sid=H5L7kq

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

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

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

图片附件: 未命名.jpg (2012-1-6 17:57, 20.5 KB) / 下载次数 1652
http://inrm3d.cn/attachment.php?aid=15998&k=578bad5c9ab4023fc6779caa15ba47ba&t=1714685954&sid=H5L7kq

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

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

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

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

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

作者: 柳烟 时间: 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) / 下载次数 1572
http://inrm3d.cn/attachment.php?aid=16000&k=a2ad7fcbfe522239fc885465112680cc&t=1714685954&sid=H5L7kq

作者: 柳烟 时间: 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) / 下载次数 1754
http://inrm3d.cn/attachment.php?aid=16003&k=00e8f8582884654a262c3a8177edf663&t=1714685954&sid=H5L7kq

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

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

附件: ManyJulia.gsp (2012-1-7 00:26, 19.39 KB) / 下载次数 2559
http://inrm3d.cn/attachment.php?aid=16007&k=aa1e76d8c8b5592e367ccc762ca39e57&t=1714685954&sid=H5L7kq

作者: 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) / 下载次数 2337
http://inrm3d.cn/attachment.php?aid=16013&k=70d6624fafad6135ec7a36d834c92a9b&t=1714685954&sid=H5L7kq

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

作者: 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) / 下载次数 2043
http://inrm3d.cn/attachment.php?aid=16025&k=decbf8bf04c5813cae0a339b73585d71&t=1714685954&sid=H5L7kq

作者: 柳烟 时间: 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) / 下载次数 2320
http://inrm3d.cn/attachment.php?aid=16035&k=a1ff9edbc6b3661f3fc43a29d6cbfb69&t=1714685954&sid=H5L7kq

图片附件: Fractal2.jpg (2012-1-9 00:10, 149.01 KB) / 下载次数 2311
http://inrm3d.cn/attachment.php?aid=16036&k=68fcd3809f2030bab29caf661c5857c7&t=1714685954&sid=H5L7kq

作者: 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) / 下载次数 1281
http://inrm3d.cn/attachment.php?aid=16045&k=893ff3ee43c93f0614dbd217918fa964&t=1714685954&sid=H5L7kq

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

作者: 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) / 下载次数 1252
http://inrm3d.cn/attachment.php?aid=16050&k=ac3f72d40ab497d4ef683d9e85f045a7&t=1714685954&sid=H5L7kq

作者: 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) / 下载次数 1259
http://inrm3d.cn/attachment.php?aid=16051&k=d38e4776458d755debb0819b917b1096&t=1714685954&sid=H5L7kq

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

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

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

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

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

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

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

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

作者: 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) / 下载次数 2644
http://inrm3d.cn/attachment.php?aid=16265&k=23bd249ce2f03bb35180d6d02c51d70b&t=1714685954&sid=H5L7kq

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

图片附件: 未命名.jpg (2012-3-2 16:30, 34.78 KB) / 下载次数 2333
http://inrm3d.cn/attachment.php?aid=16267&k=0cf775979d8c9d0cfbf1da849918e334&t=1714685954&sid=H5L7kq

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

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

图片附件: 1.jpg (2012-3-2 16:31, 27.72 KB) / 下载次数 2338
http://inrm3d.cn/attachment.php?aid=16455&k=81c33e5558bb62c367fb17042027c699&t=1714685954&sid=H5L7kq

附件: Attractor Julia修正文件.gsp (2012-3-1 21:23, 17.87 KB) / 下载次数 3322
http://inrm3d.cn/attachment.php?aid=16456&k=198857ec6836f0e4864fe9068ea4bbd8&t=1714685954&sid=H5L7kq

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

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

图片附件: 未命名.jpg (2012-3-2 16:31, 35.42 KB) / 下载次数 2324
http://inrm3d.cn/attachment.php?aid=16457&k=89a45f9629b0205ed2a6951181db6e8a&t=1714685954&sid=H5L7kq

作者: 柳烟 时间: 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) / 下载次数 1814
http://inrm3d.cn/attachment.php?aid=16470&k=77ff4bd0b1c53ee0ab247102d4d3f58a&t=1714685954&sid=H5L7kq

图片附件: 未命名.JPG (2012-3-8 12:52, 15.41 KB) / 下载次数 1766
http://inrm3d.cn/attachment.php?aid=16477&k=b688d129e812879e94f33267716f33b5&t=1714685954&sid=H5L7kq

附件: Barnsley 2 (Julia)1.gsp (2012-3-3 21:21, 18.35 KB) / 下载次数 2695
http://inrm3d.cn/attachment.php?aid=16478&k=04c176352e05f58eb1c9df53e7ee0ee3&t=1714685954&sid=H5L7kq

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

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

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

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

作者: 柳烟 时间: 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) / 下载次数 2038
http://inrm3d.cn/attachment.php?aid=16481&k=c34fb180111bf8da11cafbace010ebac&t=1714685954&sid=H5L7kq

图片附件: 未命名.JPG (2012-3-8 12:54, 20.03 KB) / 下载次数 1825
http://inrm3d.cn/attachment.php?aid=16482&k=44da05c6c46406aaaf6cf29465892119&t=1714685954&sid=H5L7kq

附件: Julifamily.gsp (2012-3-4 23:33, 26.73 KB) / 下载次数 2910
http://inrm3d.cn/attachment.php?aid=16483&k=db7536b1065d58608367de800f239b51&t=1714685954&sid=H5L7kq

作者: 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) / 下载次数 2918
http://inrm3d.cn/attachment.php?aid=16501&k=532c43b6fc4671fdb93dbe2b66324df0&t=1714685954&sid=H5L7kq

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

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

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

放大一局部一小块J集:

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

图片附件: 1.jpg (2012-3-12 12:56, 44.06 KB) / 下载次数 1797
http://inrm3d.cn/attachment.php?aid=16503&k=9dfcebebfcf8cbe0e4d0045994ccdb05&t=1714685954&sid=H5L7kq

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

图片附件: 3.jpg (2012-3-12 12:56, 37.11 KB) / 下载次数 1781
http://inrm3d.cn/attachment.php?aid=16505&k=16a5d147b812228f5823a655ab1a0918&t=1714685954&sid=H5L7kq

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

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

耍几个J集:
`1.jpg

图片附件: `1.jpg (2012-3-12 12:57, 30.45 KB) / 下载次数 1762
http://inrm3d.cn/attachment.php?aid=16513&k=ffa8d8d193078b0e44ca6e3c4cbe4d74&t=1714685954&sid=H5L7kq

图片附件: 5.jpg (2012-3-12 12:57, 38.12 KB) / 下载次数 1824
http://inrm3d.cn/attachment.php?aid=16514&k=6735136e2c37bd54910493edd1be8e9b&t=1714685954&sid=H5L7kq

图片附件: 未命名.jpg (2012-3-12 12:57, 25.26 KB) / 下载次数 1780
http://inrm3d.cn/attachment.php?aid=16522&k=0e25e560f31dc3cc70f10010862d12eb&t=1714685954&sid=H5L7kq

作者: 柳烟 时间: 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) / 下载次数 2037
http://inrm3d.cn/attachment.php?aid=16521&k=88031ed512d1820f71b1cce18d3765a5&t=1714685954&sid=H5L7kq

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

79# 柳烟

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

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

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

再扫另一个味道的图:

图片附件: 未命名.JPG (2012-3-8 12:46, 25.87 KB) / 下载次数 1472
http://inrm3d.cn/attachment.php?aid=16530&k=597d400d94ca0031a0b19b86d10d588a&t=1714685954&sid=H5L7kq

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

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

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

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

RMExpiLin.gsp (44.32 KB)

图片附件: 未命名.jpg (2012-3-12 12:50, 15.05 KB) / 下载次数 1511
http://inrm3d.cn/attachment.php?aid=16541&k=48817dacd5b64919b1247bbc1af74b2d&t=1714685954&sid=H5L7kq

图片附件: 未命名.jpg (2012-3-12 12:50, 24.74 KB) / 下载次数 1473
http://inrm3d.cn/attachment.php?aid=16557&k=1a898d593eb3b6653663db3f8e0766f0&t=1714685954&sid=H5L7kq

附件: RMExpiLin.gsp (2012-3-11 00:24, 44.32 KB) / 下载次数 2729
http://inrm3d.cn/attachment.php?aid=16558&k=b420afb2c91f6864a24f753be00ac577&t=1714685954&sid=H5L7kq

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

RMExpiNMp-O.gsp (45.37 KB)

图片附件: 未命名.jpg (2012-3-11 12:28, 112.48 KB) / 下载次数 1604
http://inrm3d.cn/attachment.php?aid=16561&k=2aaa50165edb4fd4851c6fcdbde9b953&t=1714685954&sid=H5L7kq

附件: RMExpiNMp-O.gsp (2012-3-11 12:34, 45.37 KB) / 下载次数 2492
http://inrm3d.cn/attachment.php?aid=16563&k=6ae86ddb3dea05b5464fdb8a2dab8042&t=1714685954&sid=H5L7kq

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

柳烟作品:J集:

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

图片附件: 未命名.jpg (2012-3-12 12:51, 25.82 KB) / 下载次数 1399
http://inrm3d.cn/attachment.php?aid=16566&k=6f4a177685905f58f2142254b4d1ac78&t=1714685954&sid=H5L7kq

图片附件: 未命名.jpg (2012-3-12 12:52, 28.77 KB) / 下载次数 1405
http://inrm3d.cn/attachment.php?aid=16567&k=54996172f871258f7ee8018c2529f3df&t=1714685954&sid=H5L7kq

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

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

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

柳烟作品:N集:

图片附件: 未命名.jpg (2012-3-12 12:54, 25.06 KB) / 下载次数 1407
http://inrm3d.cn/attachment.php?aid=16587&k=965d78ac4a03e54eb81783746ca7eb53&t=1714685954&sid=H5L7kq

作者: 柳烟 时间: 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) / 下载次数 1599
http://inrm3d.cn/attachment.php?aid=16635&k=ae6d47e7a01bb7ba5057d024b355f62b&t=1714685954&sid=H5L7kq

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

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

图片附件: 未命名.JPG (2012-3-17 13:58, 34.09 KB) / 下载次数 1379
http://inrm3d.cn/attachment.php?aid=16636&k=fad98970ed3ab09b386bdfa16a220184&t=1714685954&sid=H5L7kq

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

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

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

这是柳某演绎Uf的作品.

图片附件: 未命名.jpg (2012-3-20 22:52, 122.28 KB) / 下载次数 1604
http://inrm3d.cn/attachment.php?aid=16672&k=67fec892c891a72af965a890ec69907a&t=1714685954&sid=H5L7kq

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

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

图片附件: 未命名.jpg (2012-3-21 21:18, 120.11 KB) / 下载次数 3266
http://inrm3d.cn/attachment.php?aid=16682&k=03cd5eb546ffeb37611d27840087543e&t=1714685954&sid=H5L7kq

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

柳烟作品:

图片附件: 未命名.jpg (2012-3-22 20:58, 142.19 KB) / 下载次数 3352
http://inrm3d.cn/attachment.php?aid=16695&k=2143b7906911b89ebcfd34f51bd79025&t=1714685954&sid=H5L7kq

作者: 榕坚 时间: 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) / 下载次数 3371
http://inrm3d.cn/attachment.php?aid=16696&k=db831bfbf4e074c43003eb1492229c4e&t=1714685954&sid=H5L7kq

作者: 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) / 下载次数 4298
http://inrm3d.cn/attachment.php?aid=16709&k=6f40ce2cffbb8bf01665984941c0f4ca&t=1714685954&sid=H5L7kq

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

98# 柳烟

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

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

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

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