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

标题: newton集根作色法与分形阵列（征解） [打印本页]

作者: 柳烟 时间: 2012-7-23 17:02 标题: newton集根作色法与分形阵列（征解）

胡乱整一幅用根作色法调色的牛集分形，很不满意。不知有没有一本万利的分块作色法，从FSK软件来看，好象有。

z3-1N集.gsp (23.55 KB)

图片附件: 1.jpg (2015-2-15 08:19, 13.73 KB) / 下载次数 2712
http://forums.netpad.net.cn/attachment.php?aid=18035&k=c870855e92ae7786955008060b14534c&t=1765213663&sid=Ll9cjL

附件: z3-1N集.gsp (2012-7-25 15:07, 23.55 KB) / 下载次数 5146
http://forums.netpad.net.cn/attachment.php?aid=18041&k=dd7c5e8a7b06e4437263400bf4542be1&t=1765213663&sid=Ll9cjL

作者: 柳烟 时间: 2012-7-23 20:38

图片附件: 2.jpg (2015-2-15 08:20, 16.77 KB) / 下载次数 2753
http://forums.netpad.net.cn/attachment.php?aid=18038&k=3393e8a06897e6f0699c22410d862177&t=1765213663&sid=Ll9cjL

图片附件: 3.jpg (2015-2-15 08:20, 20.3 KB) / 下载次数 2674
http://forums.netpad.net.cn/attachment.php?aid=18040&k=527582382756b3aa2e8d29ef1d984619&t=1765213663&sid=Ll9cjL

作者: xiaongxp 时间: 2012-7-23 21:00

2# 柳烟
着色非常鲜亮。

作者: 柳烟 时间: 2012-7-23 21:35

3# xiaongxp
问好向老师。

如何让等势圈抹去，并出现高亮效果，目前我还没有好计谋。

图片附件: 4.jpg (2015-2-15 08:21, 21.34 KB) / 下载次数 2733
http://forums.netpad.net.cn/attachment.php?aid=18043&k=fef8ce1d11cb423636102ecbfb3d43a9&t=1765213663&sid=Ll9cjL

作者: changxde 时间: 2012-7-24 10:21

不用et可去等势圈，
用彩带可去等势圈，FSK好像也是彩带。
UF中有没有相关着色，我没发现。

作者: xiaongxp 时间: 2012-7-24 10:31

4# 柳烟
彩色N集的等势线抹去后就没有层次感了，只能设法使其变得色彩过渡更平滑。若只凸显项链结构，用浮雕法弱化外部等势线并灰度着色为佳。

图片附件: 11.jpg (2012-7-24 10:31, 47.41 KB) / 下载次数 2847
http://forums.netpad.net.cn/attachment.php?aid=18048&k=03f6bcba7af28a8766faea0cd0d270da&t=1765213663&sid=Ll9cjL

作者: 柳烟 时间: 2012-7-24 13:31

一并问好以上二位老友。以前不可思议的分形，通过大家的不懈努力并求索，不断获得解决。用根作色法，扫一张网上的分形，此分形原来作过，着色老是不尽人意，今用根作色，并用调色板颜色染之，自觉可以乱真了。

http://www.inrm3d.cn/viewthread. ... F7%D6%AE&page=1

z3-2z+2N集.gsp (25.31 KB)

图片附件: 5.jpg (2015-2-15 08:22, 18.82 KB) / 下载次数 2731
http://forums.netpad.net.cn/attachment.php?aid=18049&k=95628514d8b09983fcf14c976727960e&t=1765213663&sid=Ll9cjL

附件: z3-2z+2N集.gsp (2012-12-17 19:40, 25.31 KB) / 下载次数 4971
http://forums.netpad.net.cn/attachment.php?aid=19134&k=e29805d8491d950e95dbcb40c2ab7147&t=1765213663&sid=Ll9cjL

作者: 柳烟 时间: 2012-7-24 16:08

再扫一幅：

图片附件: 6.jpg (2015-2-15 08:23, 17.02 KB) / 下载次数 2709
http://forums.netpad.net.cn/attachment.php?aid=18050&k=e025482c29d3197e105adab25026bdd7&t=1765213663&sid=Ll9cjL

作者: 柳烟 时间: 2012-7-24 16:49

来自网上图片，函数式f(z)=z^5+2*z^4+3*z^3-4*z^2-5*z-6，如何求出这个方程的近似根呢？

图片附件: newtonbasin1.jpg (2012-7-25 07:22, 29.33 KB) / 下载次数 2747
http://forums.netpad.net.cn/attachment.php?aid=18052&k=a8f3587df109f52ea9a97b6414ddb18e&t=1765213663&sid=Ll9cjL

作者: 柳烟 时间: 2012-7-24 23:16

根作色法先要求出方程的根，这要用到常老师的三四次方程求根工具。见该论坛。勉强弄一幅：

图片附件: 7.jpg (2015-2-15 08:24, 15.86 KB) / 下载次数 2761
http://forums.netpad.net.cn/attachment.php?aid=18055&k=eea9deb150cb79ccaf554d57abc05f66&t=1765213663&sid=Ll9cjL

图片附件: 8.jpg (2015-2-15 08:24, 15.71 KB) / 下载次数 2713
http://forums.netpad.net.cn/attachment.php?aid=18079&k=3ec62fa4e70598cd6b978acb2edde46e&t=1765213663&sid=Ll9cjL

作者: 柳烟 时间: 2012-7-25 14:23

将根作色法用于牛M集，还行。

忍不住，再扫一幅：

以上各分形着色算法，均采用#1楼源文件给出的算法进行。
这是(z^2-1)(z^2-c^2)的牛M集

图片附件: 9.jpg (2015-2-15 08:25, 12.04 KB) / 下载次数 1406
http://forums.netpad.net.cn/attachment.php?aid=18057&k=ce5a3d39944bbf1b24dcde08a88fd4d6&t=1765213663&sid=Ll9cjL

图片附件: a.jpg (2015-2-15 08:25, 20.43 KB) / 下载次数 1406
http://forums.netpad.net.cn/attachment.php?aid=18063&k=13e721bd5bac350d460b4f5bc67ef1a7&t=1765213663&sid=Ll9cjL

作者: 柳烟 时间: 2012-7-26 12:10

柳烟作品：z^6-3iz^3-2的牛集：

图片附件: b.jpg (2015-2-15 08:25, 18.31 KB) / 下载次数 1367
http://forums.netpad.net.cn/attachment.php?aid=18062&k=bc6ef34a82b74917579a16c774a75ab4&t=1765213663&sid=Ll9cjL

作者: 柳烟 时间: 2012-7-26 18:05

用根作色法演绎UF中的Ikenaga Roots-Mandel：

图片附件: c.jpg (2015-2-15 08:26, 14.63 KB) / 下载次数 1386
http://forums.netpad.net.cn/attachment.php?aid=18071&k=641c8a3e998d88fe43035d98610ca95d&t=1765213663&sid=Ll9cjL

图片附件: d.jpg (2015-2-15 08:26, 14.22 KB) / 下载次数 1417
http://forums.netpad.net.cn/attachment.php?aid=18074&k=8027c93f6988b5d92a18136ad2470fdb&t=1765213663&sid=Ll9cjL

作者: 柳烟 时间: 2012-8-1 02:43

一边看伦登奥运举重，一边关注扫图。
将根作色法与光滑牛集作色法结合起来，造一幅：

图片附件: e.jpg (2015-2-15 08:27, 24.73 KB) / 下载次数 1319
http://forums.netpad.net.cn/attachment.php?aid=18174&k=05e4d9ed1869e1dd9d183e8ea3ac3d75&t=1765213663&sid=Ll9cjL

作者: 柳烟 时间: 2015-2-15 00:00

用HSV作色：

图片附件: f.jpg (2015-2-15 08:28, 16.4 KB) / 下载次数 1332
http://forums.netpad.net.cn/attachment.php?aid=23740&k=3b65effc9486764b01dd0e7ae41cf530&t=1765213663&sid=Ll9cjL

作者: 柳烟 时间: 2015-2-15 17:45

又是别样风味：

图片附件: g.jpg (2015-2-15 18:46, 22.59 KB) / 下载次数 1308
http://forums.netpad.net.cn/attachment.php?aid=23741&k=d3f3457d4c1ec435a9702013211f53e0&t=1765213663&sid=Ll9cjL

作者: 柳烟 时间: 2015-2-16 09:01

f(z)=z^5+2*z^4+3*z^3-4*z^2-5*z-6的N集，根作色。此多项式的五个根用magle软件求出。

f(z)=z5+2z4+3z3-4z2-5z-6.gsp (32.73 KB)

图片附件: g.jpg (2015-2-16 13:39, 20.8 KB) / 下载次数 1399
http://forums.netpad.net.cn/attachment.php?aid=23743&k=6d7d6a763a135df9ddbec39d6c46d70c&t=1765213663&sid=Ll9cjL

附件: f(z)=z5+2z4+3z3-4z2-5z-6.gsp (2015-2-16 09:32, 32.73 KB) / 下载次数 3063
http://forums.netpad.net.cn/attachment.php?aid=23744&k=b553e12278d78d57d524cf8bf6f9ae15&t=1765213663&sid=Ll9cjL

作者: 柳烟 时间: 2015-2-17 19:31

(z2-1)(z2-c2)的N集.gsp (25.4 KB)

图片附件: 01.jpg (2015-2-22 15:19, 17.5 KB) / 下载次数 1298
http://forums.netpad.net.cn/attachment.php?aid=23750&k=932615bdeb1e1527a1e04fd221cb1a6d&t=1765213663&sid=Ll9cjL

图片附件: 02.jpg (2015-2-22 15:19, 27.47 KB) / 下载次数 1327
http://forums.netpad.net.cn/attachment.php?aid=23751&k=1a351e7ccced9136b3f572b1e6f980d0&t=1765213663&sid=Ll9cjL

附件: (z2-1)(z2-c2)的N集.gsp (2015-2-17 19:48, 25.4 KB) / 下载次数 2970
http://forums.netpad.net.cn/attachment.php?aid=23752&k=ea30b66aa4d3a9aeece62b60dd2d2eda&t=1765213663&sid=Ll9cjL

作者: 柳烟 时间: 2015-2-19 00:06

将M集加点对称平移变换：

M集20150119.gsp (13.78 KB)

图片附件: 1.gif (2015-2-22 15:20, 62.55 KB) / 下载次数 1400
http://forums.netpad.net.cn/attachment.php?aid=23756&k=0e779e45b57c077d7306db66fcd8bdeb&t=1765213663&sid=Ll9cjL

附件: M集20150119.gsp (2015-2-19 00:15, 13.78 KB) / 下载次数 2958
http://forums.netpad.net.cn/attachment.php?aid=23757&k=2c764aa32af844f5c68cad55809b5144&t=1765213663&sid=Ll9cjL

作者: xiaongxp 时间: 2015-2-19 00:20

柳老师新年好！
柳老师好用功哦，大年年关扫图，令人感动。

作者: 柳烟 时间: 2015-2-19 01:04

20# xiaongxp
向老师新年快乐！晚上无事，又无睡意，对打麻将等又没兴趣，看看电视，玩玩电脑，混时间呗。

作者: 柳烟 时间: 2015-2-19 01:30

Mandelbrot {
;
; Generic Mandelbrot set.
;
init:
  z = @start
  c= #pixel
  x=real(c)
  y=imag(c)
  int k=0
  while (k<@TH)
  x=abs(x)-2.5
  y=abs(y)-2
  k = k+1
  endwhile
  x=abs(x)-2
  y=abs(y)-1
  c=x+flip(y)
loop:
  z = z^@power + c
bailout:
  |z| <= @bailout
default:
  title = "Mandelbrot"
  center = (-0.5, 0)
  helpfile = "Uf*.chm"
  helptopic = "Html\formulas\standard\mandelbrot.html"
$IFDEF VER50
  rating = recommended
$ENDIF
  param start
caption = "Starting point"
default = (0,0)
hint = "The starting point parameter can be used to distort the Mandelbrot \
         set. Use (0, 0) for the standard Mandelbrot set."
  endparam
  param power
caption = "Power"
default = (2,0)

  endparam
  param TH
caption = "th"
default = 5

  endparam
  float param bailout
caption = "Bailout value"
default = 4.0
min = 1.0
$IFDEF VER40
exponential = true
$ENDIF
hint = "This parameter defines how soon an orbit bails out while \
         iterating. Larger values give smoother outlines; values around 4 \
         give more interesting shapes around the set. Values less than 4 \
         will distort the fractal."
  endparam
}

图片附件: 2.gif (2015-2-22 15:22, 34.04 KB) / 下载次数 1212
http://forums.netpad.net.cn/attachment.php?aid=23758&k=c2621fd9256c875f057ece6dee4cc24c&t=1765213663&sid=Ll9cjL

作者: 柳烟 时间: 2015-2-19 01:32

如何将一个乌龟通过变换，变成指针均往左或右的乌龟阵势？有了点想法，但不知能否实现。

作者: zhongba 时间: 2015-2-19 10:00

23# 柳烟

附件: 未命名2.gsp (2015-2-19 10:00, 3.75 KB) / 下载次数 2059
http://forums.netpad.net.cn/attachment.php?aid=23759&k=e3e4f1337e6198c9b55222a5d30575e6&t=1765213663&sid=Ll9cjL

附件: 未命名3.gsp (2015-2-19 10:41, 14.68 KB) / 下载次数 2087
http://forums.netpad.net.cn/attachment.php?aid=23760&k=f45c9d963f118a484c4610ad4c8ab888&t=1765213663&sid=Ll9cjL

作者: 柳烟 时间: 2015-2-19 16:25

24# zhongba
谢谢精妙文件，设计十分巧妙。

作者: zhongba 时间: 2015-2-19 21:46

24# zhongba
精简了一下

附件: [20日凌晨再精简] M集阵列.gsp (2015-2-20 01:08, 12.35 KB) / 下载次数 1886
http://forums.netpad.net.cn/attachment.php?aid=23761&k=704a9219a7a4ec9d3031a20063b9485c&t=1765213663&sid=Ll9cjL

作者: 柳烟 时间: 2015-2-19 21:51

受24楼文件启发：
Mandelbrot同向阵列 {

init:
  z = @start
  c= #pixel
  x=real(c)
  y=imag(c)
  x=x+13-trunc(x+13)
  y=y+13-trunc(y+13)
  c=4*(x+flip(y)+(-0.7,-0.5))
loop:
z = z^@power + c
bailout:
  |z| <= @bailout
default:
  title = "Mandelbrot同向阵列"
  center = (0, 0)
  helpfile = "Uf*.chm"
  helptopic = "Html\formulas\standard\mandelbrot.html"
$IFDEF VER50
  rating = recommended
$ENDIF
  param start
caption = "Starting point"
default = (0,0)
hint = "The starting point parameter can be used to distort the Mandelbrot \
         set. Use (0, 0) for the standard Mandelbrot set."
  endparam
  param power
caption = "Power"
default = (2,0)

  endparam
  float param bailout
caption = "Bailout value"
default = 4.0
min = 1.0
$IFDEF VER40
exponential = true
$ENDIF
hint = "This parameter defines how soon an orbit bails out while \
         iterating. Larger values give smoother outlines; values around 4 \
         give more interesting shapes around the set. Values less than 4 \
         will distort the fractal."
  endparam

}

图片附件: Fractal2.png (2015-2-19 21:51, 30.71 KB) / 下载次数 1205
http://forums.netpad.net.cn/attachment.php?aid=23762&k=471bb6728d26fc1b850ec7bcd6d3d6ac&t=1765213663&sid=Ll9cjL

作者: 柳烟 时间: 2015-2-19 22:50

将就原M集，加点变换：

M集（20150219同向乌龟阵）.gsp (15.94 KB)

图片附件: 3.gif (2015-2-22 15:24, 38.29 KB) / 下载次数 1219
http://forums.netpad.net.cn/attachment.php?aid=23763&k=eb74d03613608d88ba32af998de87dd3&t=1765213663&sid=Ll9cjL

附件: M集（20150219同向乌龟阵）.gsp (2015-2-19 23:05, 15.94 KB) / 下载次数 2578
http://forums.netpad.net.cn/attachment.php?aid=23764&k=3efeb021630620c3959f6c78a8d1fa54&t=1765213663&sid=Ll9cjL

作者: 柳烟 时间: 2015-2-19 23:24

M集反向乌龟阵UF法：
Mandelbrot阵列 {
;
; Generic Mandelbrot set.
;
init:
  z = @start
  c= #pixel
  x=tan(cos(real(c)))-1
  y=tan(cos(imag(c)))
  c=x+flip(y)
loop:
  z = z^@power + c
bailout:
  |z| <= @bailout
default:
  title = "Mandelbrot反向乌龟阵"
  center = (-0.5, 0)
  helpfile = "Uf*.chm"
  helptopic = "Html\formulas\standard\mandelbrot.html"
$IFDEF VER50
  rating = recommended
$ENDIF
  param start
caption = "Starting point"
default = (0,0)
hint = "The starting point parameter can be used to distort the Mandelbrot \
         set. Use (0, 0) for the standard Mandelbrot set."
  endparam
  param power
caption = "Power"
default = (2,0)
hint = "This parameter sets the exponent for the Mandelbrot formula. \
         Increasing the real part to 3, 4, and so on, will add discs to \
         the Mandelbrot figure. Non-integer real values and non-zero \
         imaginary values will create distorted Mandelbrot sets. Use (2, 0) \
         for the standard Mandelbrot set."
  endparam
  float param bailout
caption = "Bailout value"
default = 4.0
min = 1.0
$IFDEF VER40
exponential = true
$ENDIF
hint = "This parameter defines how soon an orbit bails out while \
         iterating. Larger values give smoother outlines; values around 4 \
         give more interesting shapes around the set. Values less than 4 \
         will distort the fractal."
  endparam
switch:
  type = "Julia"
  seed = #pixel
  power = power
  bailout = bailout
}

图片附件: Fractal4.png (2015-2-19 23:24, 51.01 KB) / 下载次数 1315
http://forums.netpad.net.cn/attachment.php?aid=23765&k=c5bda8ed98a816bd53842503f4a02c00&t=1765213663&sid=Ll9cjL

作者: xiaongxp 时间: 2015-2-20 22:52

M集阵列2.gsp (12.44 KB)

图片附件: M集阵列2.jpg (2015-2-20 22:52, 37.63 KB) / 下载次数 1477
http://forums.netpad.net.cn/attachment.php?aid=23766&k=43c8085440f58f67b330f823a47ba5e4&t=1765213663&sid=Ll9cjL

附件: M集阵列2.gsp (2015-2-20 22:52, 12.44 KB) / 下载次数 3152
http://forums.netpad.net.cn/attachment.php?aid=23767&k=5a7fc2bf7377b2cd5744a5b581ac4821&t=1765213663&sid=Ll9cjL

作者: 柳烟 时间: 2015-2-21 08:16

M集阵列（20140220）.gsp (16.96 KB)

图片附件: 4.gif (2015-2-22 15:26, 40.49 KB) / 下载次数 1237
http://forums.netpad.net.cn/attachment.php?aid=23768&k=a21d4453f4e710b050fb8958a723621c&t=1765213663&sid=Ll9cjL

附件: M集阵列（20140220）.gsp (2015-2-21 08:16, 16.96 KB) / 下载次数 2961
http://forums.netpad.net.cn/attachment.php?aid=23769&k=68a97c33fd50e76b31d2eb0c1c0e2ba9&t=1765213663&sid=Ll9cjL

作者: 柳烟 时间: 2015-2-21 11:27

M集行阵（20150221(1)）.gsp (18.42 KB)

图片附件: 5.gif (2015-2-22 15:27, 52.43 KB) / 下载次数 1370
http://forums.netpad.net.cn/attachment.php?aid=23770&k=dbec4c27e6a4fa901ac8a0e32461d92a&t=1765213663&sid=Ll9cjL

附件: M集行阵（20150221(1)）.gsp (2015-2-21 11:28, 18.42 KB) / 下载次数 2599
http://forums.netpad.net.cn/attachment.php?aid=23771&k=b622115ce725653f2b13d00cbb82ba4a&t=1765213663&sid=Ll9cjL

作者: 柳烟 时间: 2015-2-21 13:38

图片附件: 6.gif (2015-2-22 15:29, 35.12 KB) / 下载次数 1255
http://forums.netpad.net.cn/attachment.php?aid=23772&k=700a0b2dbb629719216cd47784e08cb4&t=1765213663&sid=Ll9cjL

作者: 柳烟 时间: 2015-2-21 19:47

Mandelbrot阵列20150221(1号) {

init:
  z = @start
  c= 4*#pixel
  x=real(c)
  y=imag(c)
  int  k=0
  while (k<@hp)
  x=abs(x)-6
  y=abs(y)-6
  k=k+1
  endwhile
  x=abs(x)-3
  y=abs(y)-3
  m=sqrt(3)*x-y
  n=sqrt(3)*x+y
  IF n>0&&y>0
  s=2*pi/3

  elseif m>0&&y<0
  s=-2*pi/3
  endif
  c=x+flip(y)
  c=c*exp(1i*s)+1
loop:
  z = z^@power + c
bailout:
  |z| <= @bailout
default:
  title = "Mandelbrot乌龟阵20150221(1号)"
  center = (0, 0)
  helpfile = "Uf*.chm"
  helptopic = "Html\formulas\standard\mandelbrot.html"
$IFDEF VER50
  rating = recommended
$ENDIF
  param start
caption = "Starting point"
default = (0,0)
hint = "The starting point parameter can be used to distort the Mandelbrot \
         set. Use (0, 0) for the standard Mandelbrot set."
  endparam
  param power
caption = "Power"
default = (2,0)

  endparam
param hp
caption = "hp"
default = 1

  endparam
  float param bailout
caption = "Bailout value"
default = 4.0
min = 1.0
$IFDEF VER40
exponential = true
$ENDIF
hint = "This parameter defines how soon an orbit bails out while \
         iterating. Larger values give smoother outlines; values around 4 \
         give more interesting shapes around the set. Values less than 4 \
         will distort the fractal."
  endparam

}

图片附件: 7.gif (2015-2-22 15:30, 32.68 KB) / 下载次数 1244
http://forums.netpad.net.cn/attachment.php?aid=23773&k=c9d7cbe672251cbad896e6be3ad59133&t=1765213663&sid=Ll9cjL

作者: 柳烟 时间: 2015-2-21 20:09

Mandelbrot阵列20150221(2号) {

init:
  z = @start
  c= 7*#pixel
  x=real(c)
  y=imag(c)
  int  k=0
  while (k<@hp)
  x=abs(x)-6
  y=abs(y)-6
  k=k+1
  endwhile
  x=abs(x)-3
  y=abs(y)-3
  m=x-y
  n=x+y
  IF m<0&&n<0
  s=0
elseif m<0&&n>0
  s=pi/2
elseif m>0&&n>0
  s=-pi
elseif m>0&&n<0
  s=3*pi/2
  endif
  c=x+flip(y)
  c=c*exp(1i*s)+1
loop:
  z = z^@power + c
bailout:
  |z| <= @bailout
default:
  title = "Mandelbrot乌龟阵20150221(2号)"
  center = (0, 0)
  helpfile = "Uf*.chm"
  helptopic = "Html\formulas\standard\mandelbrot.html"
$IFDEF VER50
  rating = recommended
$ENDIF
  param start
caption = "Starting point"
default = (0,0)
hint = "The starting point parameter can be used to distort the Mandelbrot \
         set. Use (0, 0) for the standard Mandelbrot set."
  endparam
  param power
caption = "Power"
default = (2,0)

  endparam
param hp
caption = "hp"
default = 1

  endparam
  float param bailout
caption = "Bailout value"
default = 4.0
min = 1.0
$IFDEF VER40
exponential = true
$ENDIF
hint = "This parameter defines how soon an orbit bails out while \
         iterating. Larger values give smoother outlines; values around 4 \
         give more interesting shapes around the set. Values less than 4 \
         will distort the fractal."
  endparam
}

图片附件: Fractal4.png (2015-2-21 20:09, 65.17 KB) / 下载次数 1339
http://forums.netpad.net.cn/attachment.php?aid=23774&k=f40e3b46bcd3334d43ce3e6153198a9b&t=1765213663&sid=Ll9cjL

作者: 柳烟 时间: 2015-2-21 23:36

造阵列，受启发，终于完成了五角星科赫雪花的UF代码：
五角星科赫雪花20150221 {
init:
z = pixel
s=0.5
arg=atan2(z)
if abs(arg+pi/10)>4*pi/5
z=z*exp(1i*4*pi/5)
endif
loop:
arg=atan2(z)
x = real(z)
y = imag(z)
L13=-y+0.30902
L35=1.11803*x+1.53884*y+0.587785
L25=1.80902*x+0.587785*y-0.587785
L24=1.80902*x-0.587785*y+0.587785
L14=1.11803*x-1.53884*y-0.587785
L1=0.32492*x+y
L2=1.37638*x-y
m1=x-(1+s)*0.95106
n1=y-(1+s)*0.30902
m2=x-(1+s)*0.58779
n2=y+(1+s)*0.80902
n3=x+(1+s)*0.58779
if  abs(arg-pi/10) <pi/5
v=m1+flip(n1)
b=-pi/5
elseif abs(arg-13*pi/20) <7*pi/20
v=x+flip(y-(1+s))
b=-3*pi/5
elseif abs(arg+3*pi/10) <pi/5
v=m2+flip(n2)
b=pi/5
elseif abs(arg+pi/4)<3*pi/4
v=n3+flip(n2)
b=3*pi/5
endif
z=1/s*v*exp(1i*b)
bailout:
(L13<=0&&L24>=0&&L25>=0)||(L13<=0&&L24<0)||(L14>=0&&L13>0&&L25>=0&&L35>=0)||(L14>=0&&L35<0)||(L13>0&&L35<0&&L24<0&&L14<0)
default:
  title = "五角星科赫雪花20150221"
  helpfile = "sam-help/various.htm"
  helptopic = "sierp2"
  center = (0,0)

}

可惜的是，原始陷阱搞掉了，在UF中不知如何用代码实现。

图片附件: Fractal5.png (2015-2-21 23:37, 46.54 KB) / 下载次数 1233
http://forums.netpad.net.cn/attachment.php?aid=23775&k=9c2bd041129d519a4cc2ef9216230cf7&t=1765213663&sid=Ll9cjL

图片附件: 1.gif (2015-2-22 15:38, 47.41 KB) / 下载次数 1264
http://forums.netpad.net.cn/attachment.php?aid=23776&k=bd056a90f8ce8d1831a5e457d2ebdda3&t=1765213663&sid=Ll9cjL

作者: zhongba 时间: 2015-2-22 00:21

36# 柳烟
五角星科赫雪花中的相邻两级迭代五角星的半径比约为1-0.618

作者: 柳烟 时间: 2015-2-22 11:46

图片附件: 8.gif (2015-2-22 15:31, 32.14 KB) / 下载次数 1252
http://forums.netpad.net.cn/attachment.php?aid=23778&k=d2bb94039d66e434bc6ef6f6e810e6e9&t=1765213663&sid=Ll9cjL

作者: 柳烟 时间: 2015-2-22 23:35

MJ阵列20150222 {

init:
  z1 = @start
  c1= #pixel
  x=real(c1)
  c2=6*(c1+(23,23)-trunc(c1+(23,23)))-3
  x=real(c2)
if x>0
  z=z1
  c=1.2*(c2+(-2,-3))
  else
  z=1.2*(c2+(1.5,-3))
  c=@JJ
  endif
loop:
  z = z^@power + c
bailout:
  |z| <= @bailout
default:
  title = "MJ集阵列20150222"
  center = (0, 0)
  helpfile = "Uf*.chm"
  helptopic = "Html\formulas\standard\mandelbrot.html"
$IFDEF VER50
  rating = recommended
$ENDIF
  param start
caption = "Starting point"
default = (0,0)

  endparam
param JJ
caption = "J集定位"
default = (0.356180,0.07463)

  endparam
  param power
caption = "Power"
default = (2,0)

  endparam
param hp
caption = "hp"
default = 1

  endparam
  float param bailout
caption = "Bailout value"
default = 4.0
min = 1.0
$IFDEF VER40
exponential = true
$ENDIF
hint = "This parameter defines how soon an orbit bails out while \
         iterating. Larger values give smoother outlines; values around 4 \
         give more interesting shapes around the set. Values less than 4 \
         will distort the fractal."
  endparam

}
我想将行距压缩一下，结果搞了许久，没成功。用了向老师前面文件同向乌龟中的函数trunc.

图片附件: Fractal1.png (2015-2-22 23:35, 19.58 KB) / 下载次数 1293
http://forums.netpad.net.cn/attachment.php?aid=23779&k=4070677163bcc85b339b5e1486badf72&t=1765213663&sid=Ll9cjL

作者: 柳烟 时间: 2015-2-22 23:45

MJ集平移阵列20150223.gsp (16.04 KB)

如何在GSP中将行距减少。当我减少行距时，MJ集沿纵轴方向压扁了，不好看。

附件: MJ集平移阵列20150223.gsp (2015-2-22 23:45, 16.04 KB) / 下载次数 3113
http://forums.netpad.net.cn/attachment.php?aid=23780&k=0ea9f5b1b49e3e66d249824f5a723b50&t=1765213663&sid=Ll9cjL

图片附件: New.gif (2015-2-22 23:57, 32.68 KB) / 下载次数 1290
http://forums.netpad.net.cn/attachment.php?aid=23781&k=068a1df4e60ec42ab79bdd81782eca58&t=1765213663&sid=Ll9cjL

作者: 柳烟 时间: 2015-2-23 00:03

图片附件: 3.gif (2015-2-24 21:29, 34.46 KB) / 下载次数 1394
http://forums.netpad.net.cn/attachment.php?aid=23782&k=407fe26f4bcac814d2939cf0b4c5a691&t=1765213663&sid=Ll9cjL

作者: 柳烟 时间: 2015-2-23 17:47

M集（20150223四乌龟阵列1号）.gsp (17.58 KB)

图片附件: 2.gif (2015-2-24 21:28, 27.05 KB) / 下载次数 1400
http://forums.netpad.net.cn/attachment.php?aid=23786&k=7ce2a267bd365b86e55e25062c9069ab&t=1765213663&sid=Ll9cjL

附件: M集（20150223四乌龟阵列1号）.gsp (2015-2-23 17:52, 17.58 KB) / 下载次数 3169
http://forums.netpad.net.cn/attachment.php?aid=23787&k=0081f2945b4d72ff2a4b5fd6e629622f&t=1765213663&sid=Ll9cjL

作者: 柳烟 时间: 2015-2-23 18:57

Mandelbrot 20150223四乌龟阵1号{
  z = @start
  c=#pixel
c=3*(1.5*c+(12,12)-trunc(1.5*c+(12,12)))+(-1.5,-1.5)
x=real(c)
y=imag(c)
m=x-y
n=x+y
if m<0&&n<0
  s=pi/2
elseif m<0&&n>0
s=pi
elseif m>0&&n>0
s=-pi/2
endif
c=3*c*exp(1i*s)+(-0.5,1.8)
loop:
  z = z^@power + c
bailout:
  |z| <= @bailout
default:
  title = "Mandelbrot四乌龟阵列20150223一号"
  center = (0, 0)
  helpfile = "Uf*.chm"
  helptopic = "Html\formulas\standard\mandelbrot.html"
$IFDEF VER50
  rating = recommended
$ENDIF
  param start
caption = "Starting point"
default = (0,0)
hint = "The starting point parameter can be used to distort the Mandelbrot \
         set. Use (0, 0) for the standard Mandelbrot set."
  endparam
  param power
caption = "Power"
default = (2,0)
hint = "This parameter sets the exponent for the Mandelbrot formula. \
         Increasing the real part to 3, 4, and so on, will add discs to \
         the Mandelbrot figure. Non-integer real values and non-zero \
         imaginary values will create distorted Mandelbrot sets. Use (2, 0) \
         for the standard Mandelbrot set."
  endparam
  float param bailout
caption = "Bailout value"
default = 1000000.0
min = 1.0
$IFDEF VER40
exponential = true
$ENDIF
hint = "This parameter defines how soon an orbit bails out while \
         iterating. Larger values give smoother outlines; values around 4 \
         give more interesting shapes around the set. Values less than 4 \
         will distort the fractal."
  endparam

}

图片附件: Fractal2.png (2015-2-23 19:01, 30.61 KB) / 下载次数 1396
http://forums.netpad.net.cn/attachment.php?aid=23788&k=caa9cbd53680b1470ce269c226aaa56d&t=1765213663&sid=Ll9cjL

作者: 柳烟 时间: 2015-2-23 23:03

M集对称变换4乌龟（20150223四乌龟阵列2号）.gsp (18.61 KB)

附件: M集对称变换4乌龟（20150223四乌龟阵列2号）.gsp (2015-2-23 23:03, 18.61 KB) / 下载次数 3406
http://forums.netpad.net.cn/attachment.php?aid=23789&k=f1d65612377271f5ecfeaaf59a61fd35&t=1765213663&sid=Ll9cjL

图片附件: 1.gif (2015-2-24 21:27, 29.27 KB) / 下载次数 1387
http://forums.netpad.net.cn/attachment.php?aid=23790&k=2ed3834e5fde3b7433872ab4356ce368&t=1765213663&sid=Ll9cjL

作者: 柳烟 时间: 2015-2-23 23:21

Mandelbrot 20150223四乌龟阵2号{
  z = @start
  c=4.5*#pixel
  K=0
while (k<@hp)
  c=abs(c)-(3,3)
  k=k+1
  endwhile
c=abs(c)-(1.5,1.5)
m=real(c)-imag(c)
n= real(c)+imag(c)
if m<0&&n<0
  s=pi/2
elseif m<0&&n>0
s=pi
elseif m>0&&n>0
s=-pi/2
endif
c=3*c*exp(1i*s)+(-0.5,1.8)
loop:
  z = z^@power + c
bailout:
  |z| <= @bailout
default:
  title = "Mandelbrot四乌龟阵列20150223二号"
  center = (0, 0)
  helpfile = "Uf*.chm"
  helptopic = "Html\formulas\standard\mandelbrot.html"
$IFDEF VER50
  rating = recommended
$ENDIF
  param start
caption = "Starting point"
default = (0,0)
hint = "The starting point parameter can be used to distort the Mandelbrot \
         set. Use (0, 0) for the standard Mandelbrot set."
  endparam
param hp
caption = "hp"
default = 1
  endparam
  param power
caption = "Power"
default = (2,0)
hint = "This parameter sets the exponent for the Mandelbrot formula. \
         Increasing the real part to 3, 4, and so on, will add discs to \
         the Mandelbrot figure. Non-integer real values and non-zero \
         imaginary values will create distorted Mandelbrot sets. Use (2, 0) \
         for the standard Mandelbrot set."
  endparam
  float param bailout
caption = "Bailout value"
default = 1000000.0
min = 1.0
$IFDEF VER40
exponential = true
$ENDIF
hint = "This parameter defines how soon an orbit bails out while \
         iterating. Larger values give smoother outlines; values around 4 \
         give more interesting shapes around the set. Values less than 4 \
         will distort the fractal."
  endparam

}

图片附件: Fractal2.png (2015-2-23 23:25, 46.41 KB) / 下载次数 1393
http://forums.netpad.net.cn/attachment.php?aid=23791&k=374bd346dadd4ffcaa6c884f39c29b9c&t=1765213663&sid=Ll9cjL

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