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

标题: UF分形中的画板实现第3集 [打印本页]

作者: 柳烟 时间: 2010-8-17 07:14 标题: UF分形中的画板实现第3集

General_Sterling_18.gsp (110 KB)

图片附件: 未命名1.JPG (2010-8-17 07:14, 25.21 KB) / 下载次数 1794
http://inrm3d.cn/attachment.php?aid=6570&k=407d297d085f8b7dce22c23da80f7171&t=1732398679&sid=462I2T

附件: General_Sterling_18.gsp (2010-8-17 07:16, 110 KB) / 下载次数 3384
http://inrm3d.cn/attachment.php?aid=6571&k=f72dcf87f26e22b8260e8e46331aef15&t=1732398679&sid=462I2T

图片附件: 未命名2.JPG (2010-8-17 07:23, 27.24 KB) / 下载次数 1859
http://inrm3d.cn/attachment.php?aid=6572&k=65e0da89aef989d4ed5a56cfa5438494&t=1732398679&sid=462I2T

图片附件: 未命名3.JPG (2010-8-17 07:32, 46.66 KB) / 下载次数 1804
http://inrm3d.cn/attachment.php?aid=6573&k=13a632c59c2691fd15ee7dc23e186ca5&t=1732398679&sid=462I2T

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

附UF中的代码：
General_Sterling_18 {

; Original: z=z-(z*z2+z2*c-z+c-zsin)/(3*z2+2*c*z-1-zcos)

init:
  float zds
  complex z1
  complex z2

if @usepixel
z = #pixel*@seed
  else
;Original From Sterling
z = @seed
  endif

loop:
  z1 = z
  z2 = z^2

  z = z - (@p1*z*z2+@p2*z2*#pixel-z+@p4*#pixel-@p3*@fn1(z))/(@p5*z2+@p6*#pixel*z-1-@p7*@fn2(z))
  zds = |z1-z|

bailout:
  (zds > 1e-25) && (zds < 1e20)

default:
  title = "General_Sterling_18"
  periodicity = 0
  maxiter = 190

  float param Seed
caption = "Seed"
default = -1/3
  endparam

  complex func fn1
caption = "Fn1"
hint = "This is the first function."
default = sin()
  endfunc

  complex func fn2
caption = "Fn2"
hint = "This is the Second function."
default = cos()
  endfunc

  complex param p1
caption = "p1"
default = 1
  endparam

  complex param p2
caption = "p2"
default = 1
  endparam

  complex param p3
caption = "p3"
default = 1
  endparam

  complex param p4
caption = "p4"
default = 1
  endparam

  complex param p5
caption = "p5"
default = 3
  endparam

  complex param p6
caption = "p6"
default = 2
  endparam

  complex param p7
caption = "p7"
default = 1
  endparam

  Bool param usepixel
caption = "Use Pixels (C*Seed)"
default = True
  endparam
}

作者: 柳烟 时间: 2010-8-17 11:13

将就一楼的文件，改变参数，再扫一图：

图片附件: 未命名1.JPG (2010-8-17 11:13, 17.98 KB) / 下载次数 1777
http://inrm3d.cn/attachment.php?aid=6576&k=bdb7d4414263c1dc254af1bb1b0130e5&t=1732398679&sid=462I2T

作者: xiaongxp 时间: 2010-8-17 12:41

看不懂源代码，只有欣赏的份儿哟。

作者: 柳烟 时间: 2010-8-17 19:29

看不懂源代码，只有欣赏的份儿哟。
xiaongxp 发表于 2010-8-17 12:41

我也不太懂代码,因为天天看代码,对合面板上的东西,能悟之三四.

作者: 柳烟 时间: 2010-8-17 19:48

代码如下:
bwp4 {
init:
  z = @InitValue
  t = #pixel
loop:
  z = sin(z) - t
  t = 1 / (z * @Value)
bailout:
  |z| < @Bailout
default:
  title = "BWP 04"
  param InitValue
caption = "Initial Value"
default = (0, 0)
  endparam
  param Value
caption = "Value"
default = (50, 0)
  endparam
  param Bailout
caption = "Bailout Value"
default = 4.0
  endparam
}

BWP 04.gsp (34.31 KB)

图片附件: 未命名1.JPG (2010-8-17 19:48, 41.04 KB) / 下载次数 1800
http://inrm3d.cn/attachment.php?aid=6587&k=2f256468ca5c8022c3c60538f45ba565&t=1732398679&sid=462I2T

图片附件: 未命名2.JPG (2010-8-17 19:59, 39.68 KB) / 下载次数 1828
http://inrm3d.cn/attachment.php?aid=6588&k=d87a1f8abfc6dd5978eac407a51ba6aa&t=1732398679&sid=462I2T

附件: BWP 04.gsp (2010-8-17 19:59, 34.31 KB) / 下载次数 3444
http://inrm3d.cn/attachment.php?aid=6589&k=865f30a0ad8adf2b4611f6ab6f088f78&t=1732398679&sid=462I2T

作者: 柳烟 时间: 2010-8-18 00:37

将上楼文件放大后，扫出一图：

图片附件: 未命名3.JPG (2010-8-18 00:37, 40.9 KB) / 下载次数 1800
http://inrm3d.cn/attachment.php?aid=6607&k=d977d154a4a89d4edd863c58a2cd40e4&t=1732398679&sid=462I2T

作者: 榕坚 时间: 2010-8-18 11:34

3# 柳烟

怎么会有那小宝贝出现？改了哪些参数呢？

作者: 柳烟 时间: 2010-8-18 15:15

8# 榕坚
将源文件中的有关sinz、cosz的系数，好象是p3与p7设置为0

作者: 柳烟 时间: 2010-8-18 15:17

代码如下：
bwp7 {
init:
  z = @InitValue
  t = #pixel
loop:
  z = tanh(z) - t
  t = 1 / (z * @Value)
bailout:
  |z| < @Bailout
default:
  title = "BWP 07"
  param InitValue
caption = "Initial Value"
default = (0, 0)
  endparam
  param Value
caption = "Value"
default = (50, 0)
  endparam
  param Bailout
caption = "Bailout Value"
default = 4.0
  endparam
}

BWP07.gsp (11.47 KB)

图片附件: 未命名1.JPG (2010-8-18 15:17, 41.62 KB) / 下载次数 1783
http://inrm3d.cn/attachment.php?aid=6631&k=1ec2e83ffba3d37e0e3a57563835fed0&t=1732398679&sid=462I2T

图片附件: 未命名1.JPG (2010-8-18 15:34, 29.35 KB) / 下载次数 1826
http://inrm3d.cn/attachment.php?aid=6632&k=370e0bf0d6cc2fd45105524047a1b20f&t=1732398679&sid=462I2T

附件: BWP07.gsp (2010-8-18 15:39, 11.47 KB) / 下载次数 3455
http://inrm3d.cn/attachment.php?aid=6633&k=07b7bd29d831d7f6e46a34ca26a6c6a9&t=1732398679&sid=462I2T

作者: 柳烟 时间: 2010-8-18 15:27

我扫出的图片有时挺不错的，但是发在坛子上，颜色变来不好看，我是用windows中的附件处理的图片，请教板友们，有什么更好的处理gsp中的图片法子？

作者: 柳烟 时间: 2010-8-18 16:45

再扫一图：

图片附件: 未命名3.JPG (2010-8-18 16:45, 45.45 KB) / 下载次数 1246
http://inrm3d.cn/attachment.php?aid=6637&k=b0addca230cb3df146be1f0fd173b548&t=1732398679&sid=462I2T

图片附件: 未命名1.JPG (2010-8-18 22:51, 45.6 KB) / 下载次数 1225
http://inrm3d.cn/attachment.php?aid=6641&k=13db47c62ee251eef216f84a5735f0fe&t=1732398679&sid=462I2T

作者: 柳烟 时间: 2010-8-18 22:04

代码：
bwp3 {
init:
  z1 = @InitValue
  t = #pixel
loop:
  z2 = z1
  z1 = sin(z1) - t
  t = 1 / (z1 * @Value)
  z = z1 - z2
bailout:
  |z| < @Bailout
default:
  title = "BWP 03"
  param InitValue
caption = "Initial Value"
default = (0, 0)
  endparam
  param Value
caption = "Value"
default = (50, 0)
  endparam
  param Bailout
caption = "Bailout Value"
default = 4.0
  endparam
}

BWP03.gsp (17.65 KB)

图片附件: 未命名1.JPG (2010-8-18 22:04, 37.54 KB) / 下载次数 1383
http://inrm3d.cn/attachment.php?aid=6638&k=aaee340174ed89fbf8f720cc6736b112&t=1732398679&sid=462I2T

附件: BWP03.gsp (2010-8-18 22:04, 17.65 KB) / 下载次数 2723
http://inrm3d.cn/attachment.php?aid=6639&k=c3129743cd3bad871b702857d8447073&t=1732398679&sid=462I2T

作者: 柳烟 时间: 2010-8-18 23:12

BMP03的changxde 变换图片

BWP03的变换图片.gsp (19.27 KB)

图片附件: 未命名1.JPG (2010-8-18 23:12, 34.4 KB) / 下载次数 1380
http://inrm3d.cn/attachment.php?aid=6643&k=9e562dcf2c4d15f2f98c041b6a81595d&t=1732398679&sid=462I2T

附件: BWP03的变换图片.gsp (2010-8-18 23:17, 19.27 KB) / 下载次数 2692
http://inrm3d.cn/attachment.php?aid=6644&k=0dbc0d0c1bb4e3502827004dcb22cfaf&t=1732398679&sid=462I2T

作者: 柳烟 时间: 2010-8-18 23:51

图片附件: 未命名2.JPG (2010-8-18 23:51, 42.68 KB) / 下载次数 1371
http://inrm3d.cn/attachment.php?aid=6649&k=8f082d58e59cd887d639c63eefd90a20&t=1732398679&sid=462I2T

作者: xiaongxp 时间: 2010-8-19 18:06

12# 柳烟
这组图若用mjj的起泡法作色，将所有圆作成有光泽的球，一定更加漂亮。

作者: 柳烟 时间: 2010-8-19 23:44

16# xiaongxp
向老师建议好，由于近来忙，对梅老师的起泡法还未曾钻研，等以后研究了，再来作。

作者: 柳烟 时间: 2010-8-19 23:47

代码：
bwp16 {
init:
  z = @InitValue
loop:
  z = @Mult3 * ((@Mult1 * (z ^ @Exp1)) + \
  (@Mult2 * (z ^ @Exp2))) + @Fun(#pixel)
bailout:
  |z| < @Bailout
default:
  title = "BWP 16"
  param InitValue
caption = "Initial Value"
default = (0, 0)
  endparam
  param Mult1
caption = "Multiplier 1"
default = (0.5, 0)
  endparam
  param Mult2
caption = "Multiplier 2"
default = (0.5, 0)
  endparam
  param Mult3
caption = "Multiplier 3"
default = (1, 0)
  endparam
  param Exp1
caption = "Exponent 1"
default = (2, 0)
  endparam
  param Exp2
caption = "Exponent 2"
default = (2, 0)
  endparam
  param Bailout
caption = "Bailout Value"
default = 4.0
  endparam
  func Fun
caption = "Function"
default = ident()
  endfunc
switch:
  type = "bwp17"
  Seed = #pixel
  Mult1 = @Mult1
  Mult2 = @Mult2
  Mult3 = @Mult3
  Exp1 = @Exp1
  Exp2 = @Exp2
  Fun = @Fun
  Bailout = @Bailout
}
扫出的图片多如繁星，我这文件，改变参数代表的开关顶，就可扫出不同的图片，有些是UF中能扫出的，有些是不能扫出的，先扫几张UF中有的。连我自已都扫不赢，大家自已扫去。

BWP16一号.gsp (33.05 KB)

未命名2.JPG

附件: BWP16一号.gsp (2010-8-19 23:47, 33.05 KB) / 下载次数 2632
http://inrm3d.cn/attachment.php?aid=6671&k=06b2e22b060dd4ed5627166cfd39dfe8&t=1732398679&sid=462I2T

图片附件: 未命名1.JPG (2010-8-19 23:47, 18.59 KB) / 下载次数 1106
http://inrm3d.cn/attachment.php?aid=6672&k=f0422b15682b826b09b22467b5902a82&t=1732398679&sid=462I2T

图片附件: 未命名2.JPG (2010-8-19 23:47, 36.16 KB) / 下载次数 1231
http://inrm3d.cn/attachment.php?aid=6673&k=02cd92b665519fb478693f654db33578&t=1732398679&sid=462I2T

作者: 柳烟 时间: 2010-8-20 01:21

下略扫扫几张UF中没有，而楼上文件中，改变参数就能得到的图片：
未命名6.JPG

图片附件: 未命名5.JPG (2010-8-20 01:21, 12.49 KB) / 下载次数 1149
http://inrm3d.cn/attachment.php?aid=6675&k=e67de917b4883d85f688040fe4ca9c46&t=1732398679&sid=462I2T

图片附件: 未命名6.JPG (2010-8-20 01:31, 9.78 KB) / 下载次数 1207
http://inrm3d.cn/attachment.php?aid=6676&k=cc35cf5e3fa2bb35adbe5cd6551d4f9a&t=1732398679&sid=462I2T

图片附件: 未命名7.JPG (2010-8-20 02:55, 31.55 KB) / 下载次数 1176
http://inrm3d.cn/attachment.php?aid=6677&k=84ac64af14ac653af1b5dff5f0fd19f9&t=1732398679&sid=462I2T

作者: 柳烟 时间: 2010-8-20 10:13

BMP16 二号.gsp (51.93 KB)
未命名3.JPG

图片附件: 未命名1.JPG (2010-8-20 10:13, 14.11 KB) / 下载次数 1200
http://inrm3d.cn/attachment.php?aid=6679&k=ee5328bca63d209bfec3e2e07145bd43&t=1732398679&sid=462I2T

图片附件: 未命名3.JPG (2010-8-20 10:45, 15.78 KB) / 下载次数 1208
http://inrm3d.cn/attachment.php?aid=6681&k=6c14d5fbbb78d7a2a42bf548a84d7a3a&t=1732398679&sid=462I2T

附件: BMP16 二号.gsp (2010-8-20 10:45, 51.93 KB) / 下载次数 2565
http://inrm3d.cn/attachment.php?aid=6682&k=6a56e376969e4e9ef1dacb3f557bf791&t=1732398679&sid=462I2T

图片附件: 未命名3.JPG (2010-8-20 11:09, 21.56 KB) / 下载次数 1177
http://inrm3d.cn/attachment.php?aid=6684&k=9dd36b07b794b7fb98fd9b3920c4a39c&t=1732398679&sid=462I2T

作者: 柳烟 时间: 2010-8-20 11:48

bwp24 {
init:
  z = @InitValue
loop:
  z = (z ^ @Exponent + #pixel) / (@Value + #pixel)
bailout:
  |z| < @Bailout
default:
  title = "BWP 24"
  param InitValue
caption = "Initial Value"
default = (0, 0)
  endparam
  param Exponent
caption = "Exponent"
default = (2, 0)
  endparam
  param Value
caption = "Value"
default = (1, 0)
  endparam
  param Bailout
caption = "Bailout Value"
default = 4.0
  endparam
switch:
  type = "bwp25"
  Seed = #pixel
  Exponent = @Exponent
  Value = @Value
  Bailout = @Bailout
}

BWP24.gsp (17.88 KB)

附件: BWP24.gsp (2010-8-20 11:48, 17.88 KB) / 下载次数 2239
http://inrm3d.cn/attachment.php?aid=6685&k=727fc4f5c2f5c91ea3e7a6ae92d2da3d&t=1732398679&sid=462I2T

图片附件: 未命名5.JPG (2010-8-20 11:48, 19.86 KB) / 下载次数 1021
http://inrm3d.cn/attachment.php?aid=6686&k=dbe224a4521d2fe062e4fc8b6a1af0ea&t=1732398679&sid=462I2T

作者: 柳烟 时间: 2010-8-26 00:58

a115 {
;Natural symmetric fractals
init:
  z = c = #pixel
  if @modeJ == "Mandel"
;z = @seed ;optional settings for
;c = #pixel  ;mandel mode
  endif
  if @modeM == "Julia"
z = pixel ;Now julia mode is active
c = @seed
  endif
loop:
  z=@fn1(@fn2(z*c)^@p2)^@p1 / @fn3(@fn4(c*z)^@p4)^@p3
  z=z^@p5*c
bailout:
|z| < @bailo
default:
  param modeM  ;Mode holder
caption="switchmode"
enum = "Mandel""Julia"
default = 0
enabled = false
  endparam
  param modeJ  ;Mode holder
caption="switchto"
enum = "Mandel""Julia"
default = 1
enabled = false
  endparam

float param bailo
caption="bailout"
  default = 50
endparam
func fn1
default = sin()
endfunc
float param p1
  caption="powerFn1"
  default=0
endparam
func fn2
default = asin()
endfunc
float param p2
  caption="powerFn2"
  default=1
endparam
func fn3
default = sin()
endfunc
float param p3
  caption="powerFn3"
  default=1
endparam
func fn4
default = sqr()
endfunc
float param p4
  caption="powerFn4"
  default=1
endparam
float param p5
  caption="powerOfZ"
  default=1
endparam
param seed
  caption = "Seed"
  visible=(@modeM == "Julia")
  default = (1,0)
endparam
switch:
  type= "a115"
  seed= #pixel
  bailo=bailo
  modeJ=@modeM ;Everytime when switched occured
  modeM=@modeJ ;values of modeM and modeJ will be
  fn1=fn1    ;exchanged.
  p1=p1
  fn2=fn2
  p2=p2
  fn3=fn3
  p3=p3
  fn4=fn4
  p4=p4
  p5=p5

}

由于误读了程序中的一个细微地方，没造出UF中的那图，但造出的图仍然美不胜收。

a115分形的变种.gsp (17.57 KB)

图片附件: mset1.jpg (2010-8-26 00:58, 54.92 KB) / 下载次数 1159
http://inrm3d.cn/attachment.php?aid=6803&k=47cc49dc3433fd9a274a1aeec4e5cb0f&t=1732398679&sid=462I2T

附件: a115分形的变种.gsp (2010-8-26 01:02, 17.57 KB) / 下载次数 2239
http://inrm3d.cn/attachment.php?aid=6804&k=bebff261a2a7d7f14e1d137629133ee6&t=1732398679&sid=462I2T

作者: 柳烟 时间: 2010-8-26 01:27

接着造，仍未达目的，但扫出这张，也美，那里出问题，还不知。

图片附件: mset1.jpg (2010-8-26 01:27, 42.69 KB) / 下载次数 992
http://inrm3d.cn/attachment.php?aid=6805&k=c4f93f88879b03a26665c2ebc855c9ac&t=1732398679&sid=462I2T

作者: 榕坚 时间: 2010-8-29 21:43

23# 柳烟

问题是否出在那个asin函数上，复数的反三角是怎么定义的呢？

作者: 柳烟 时间: 2010-8-29 23:12

24# 榕坚
我曾用此站上的老师造的此工具，用来作其他的分形，没问题呀。

作者: 柳烟 时间: 2010-8-30 00:46

UF的效果图是这样的：

图片附件: 未命名.JPG (2010-8-30 00:46, 40.46 KB) / 下载次数 985
http://inrm3d.cn/attachment.php?aid=6896&k=9a0201d242c46fa8fd36efae2f10b033&t=1732398679&sid=462I2T

作者: 榕坚 时间: 2010-8-31 09:46

这个代码很是奇怪，那个@seed根本没有作用。面板上的开关项有些是灰色不能选择的，读代码可不能按步就班。

作者: 柳烟 时间: 2010-8-31 12:46

我感觉到代码中给出的一长串函数表达式有问题.

作者: 榕坚 时间: 2010-8-31 14:38

28# 柳烟

表达式应该没有问题，我试过了。就是Z与C的问题了，还有就是visible=(@modeM == "Julia")这个代码可能就是造成面板上开关项的灰色。

作者: 榕坚 时间: 2010-8-31 17:51

上午试了一下出来的结果还有点象，只是有些空白而且中间的圆圈没有连完整。接下来想把空白弄掉，就把隐藏的点显示出来想继续修改，结果跳出一个几何画板内部文件出错，眼睁睁地看着弄了一个上午的文件泡汤了。下午继续从头开始弄了两遍都是与楼主类似的结果。感到莫名其妙。
谁有什么好办法可以解决几何画板跳出出错对画框后类似office那样可以恢复的，我已经遇到好多次了。都是在显示隐藏的点时出现的，特别是当要显示合并后的隐藏点时经常出问题。我用的是英文5.01版

图片附件: 捕获a115.JPG (2010-8-31 17:51, 60.07 KB) / 下载次数 1142
http://inrm3d.cn/attachment.php?aid=6936&k=c8515fd94a269e7885cf1d7c79735974&t=1732398679&sid=462I2T

图片附件: 捕获1.JPG (2010-8-31 17:51, 58.88 KB) / 下载次数 1149
http://inrm3d.cn/attachment.php?aid=6937&k=24bb8c5cf0545af6bb03aae5a084bbc7&t=1732398679&sid=462I2T

图片附件: 捕获2.JPG (2010-8-31 17:51, 40.93 KB) / 下载次数 1031
http://inrm3d.cn/attachment.php?aid=6938&k=3bdc8600d44126f1a7d661a8f542f93c&t=1732398679&sid=462I2T

作者: 柳烟 时间: 2010-8-31 19:03

老兄所说的问题,我还没遇到过.上楼第一张图片对头了,按面板上的参数,结合代吗,你用的是不是函数:
c/sin(c^2.z^2)呢？我试了好几次，都不能达到UF中的效果。

作者: 榕坚 时间: 2010-8-31 20:10

那些函数我是照着做的，如果要化简的话应该化为：(Z*C)^0/sin((Z*C)^2)，注意：(Z*C)^0不能用1代替。我在UF中试过了。

作者: xiaongxp 时间: 2010-8-31 20:18

非常佩服二位的智慧和毅力，老巷望尘莫及。

作者: 柳烟 时间: 2010-8-31 22:02

32# 榕坚
我就是用1代替的,可能错误正在于此,为1应该有条件，我再重新作作看。

作者: 柳烟 时间: 2010-8-31 22:04

33# xiaongxp
向老师是我入门导师,柳某终生难忘,分形确实魅力引人,总是觉得时间不够用,要研究的问题太多了.

作者: xiaongxp 时间: 2010-8-31 22:17

柳老师过谦了，大家是一起切磋板技，相互成就，共同提高的。你们对问题的矜持、对UF源文件的解读能力，真是令我钦佩。这一点，我是自叹弗如，永远也学不会的。

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