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

标题: 复变分形陷阱征解 [打印本页]

作者: 柳烟 时间: 2013-3-19 20:02 标题: 复变分形陷阱征解

光滑N集.gsp (19.44 KB)

图片附件: 未命名.jpg (2013-3-20 08:19, 20.02 KB) / 下载次数 1524
http://develop.forums.netpad.net.cn/attachment.php?aid=19630&k=6b98f155f650538238c5fc2c77a5939a&t=1750241456&sid=HzJBH6

附件: 光滑N集.gsp (2013-3-19 20:04, 19.44 KB) / 下载次数 3272
http://develop.forums.netpad.net.cn/attachment.php?aid=19631&k=c394bb5569ab9a8f4c690f0c51c2b18e&t=1750241456&sid=HzJBH6

作者: 柳烟 时间: 2013-3-19 21:01

将就上楼文件，略作修改，立马扫得旧日扫过的旧图，别有一番滋味：

图片附件: 未命名.jpg (2013-3-20 08:20, 33.46 KB) / 下载次数 1536
http://develop.forums.netpad.net.cn/attachment.php?aid=19632&k=cf3381ae15f14ed51eb5899c61b130d4&t=1750241456&sid=HzJBH6

图片附件: 未命名.jpg (2013-3-20 08:21, 25.51 KB) / 下载次数 1538
http://develop.forums.netpad.net.cn/attachment.php?aid=19633&k=d12ca8c72cc7a71d9eb0c7737e73ee39&t=1750241456&sid=HzJBH6

作者: 柳烟 时间: 2013-3-19 22:30

我想光滑牛顿M集的等势圈，结果惨败。上二楼光滑的牛顿J集。

作者: 柳烟 时间: 2013-3-20 20:30

研究复分形作陷阱，扫图一幅：

图片附件: 未命名.jpg (2013-3-21 00:31, 31.02 KB) / 下载次数 1533
http://develop.forums.netpad.net.cn/attachment.php?aid=19634&k=eef1756307be9a573186d9cf829bef62&t=1750241456&sid=HzJBH6

作者: 柳烟 时间: 2013-3-20 23:13

我将UF中的复分形陷阱文件，抽出M集陷阱，代码如下：
juliatrap1(BOTH) {
; By Samuel Monnier, 1999
init:
  float dz = 0.0
  float stdz = 0.0
  start = @seed

loop:
  trz = exp(flip(-pi/180*@rot))*(#z-@center)/@size
  x = 1/sqrt(@r)*real(trz)
  y = sqrt(@r)*imag(trz)
  if @freq != 0
x = 2*sin(2*pi*@freq*x)
y = 2*sin(2*pi*@freq*y)
  endif

  start = x + flip(y)

  iter = @niter
  zz = @seed

  float i = 0
  while i < @niter
i = i + 1
zz = zz^2 + start
if |zz| > 1e20
   iter = i
   i = @niter
endif

  endwhile
  float logp = 1/log(2)
  float logb = log(log(1e20))

dz = real(iter + logp*logb - logp*log(log(cabs(zz))))

if dz >stdz
   stdz = dz
endif

final:
  #index = .1*stdz
default:
  title = "Copy of Julia Trap"
  helpfile = "sam-help/juliatrap.htm"

  param seed
caption = "Seed (for 'custom')"
default = (0,0)
  endparam
  param center
caption = "Center"
default = (1,1)
  endparam
  param rot
caption = "Rotation"
default = -30.0
  endparam
  param size
caption = "Size"
default = 0.38750
  endparam
  param r
caption = "Ratio Width/Heigh"
default = 1.0
  endparam
  param @niter
caption = "Julia Iterations"
default = 100
  endparam
  param freq
caption = "Trap Frequency"
default = 0.0
  endparam
}
此代码，我用于谢氏三角形，得到上楼效果，但我将其用于牛集时，连无数小乌龟影子都没看到，怪哉。

作者: 柳烟 时间: 2013-3-20 23:20

UF中的效果图是：

图片非常有趣。我干了两天了，没招。

图片附件: Fractal2.JPG (2013-3-21 00:32, 37.75 KB) / 下载次数 1543
http://develop.forums.netpad.net.cn/attachment.php?aid=19635&k=81f83476d80e1d923ca729d1bb95d8c3&t=1750241456&sid=HzJBH6

作者: 柳烟 时间: 2013-3-21 21:17

M集陷阱，按代码作出后，参数总是与面板上不一致，试试M集坠入M集陷阱，作出的总是不尽人意，将M集用于N集，更相差甚远……等会我将文件发上来，兄台们看看怎么回事。

M集陷阱尝试（感觉有问题）.gsp (30.98 KB)

图片附件: 未命名.jpg (2013-3-22 09:51, 28.6 KB) / 下载次数 1540
http://develop.forums.netpad.net.cn/attachment.php?aid=19640&k=b18721b1e419b4a2fab8ff25d92bb89e&t=1750241456&sid=HzJBH6

附件: M集陷阱尝试（感觉有问题）.gsp (2013-3-21 22:07, 30.98 KB) / 下载次数 3379
http://develop.forums.netpad.net.cn/attachment.php?aid=19641&k=6e1d69a63181f0b7f1911f108a9433e2&t=1750241456&sid=HzJBH6

作者: 柳烟 时间: 2013-3-22 21:05

最近在GSP中玩的晦气，在UF中整点代码，造了个谢尔兵司机陷阱，扫一个M集，颇不错的。估计在秀气画板中造出后，运行不起走了。

图片附件: Fractal1.JPG (2013-3-23 07:51, 45.26 KB) / 下载次数 1537
http://develop.forums.netpad.net.cn/attachment.php?aid=19645&k=49956ade7c17ea612ef4178a199c9d2b&t=1750241456&sid=HzJBH6

作者: 柳烟 时间: 2013-3-23 00:17

另一种谢尔宾斯基陷阱：

图片附件: Fractal1.JPG (2013-3-23 07:52, 40.07 KB) / 下载次数 1569
http://develop.forums.netpad.net.cn/attachment.php?aid=19646&k=867d2582f630cf7987f4e8b20658e111&t=1750241456&sid=HzJBH6

作者: 柳烟 时间: 2013-3-23 00:24

太美了：

图片附件: Fractal1.JPG (2013-3-23 07:53, 48.17 KB) / 下载次数 1566
http://develop.forums.netpad.net.cn/attachment.php?aid=19647&k=9149763c92df7ef5f213f3b79da41b43&t=1750241456&sid=HzJBH6

图片附件: Fractal1.jpg (2013-3-23 00:30, 179.15 KB) / 下载次数 1769
http://develop.forums.netpad.net.cn/attachment.php?aid=19648&k=8efde0ca53e2b76fe0e0aab0258e3a37&t=1750241456&sid=HzJBH6

作者: 柳烟 时间: 2013-3-23 21:32

N集陷阱：

图片附件: Fractal2.png (2013-3-23 21:32, 173.27 KB) / 下载次数 1697
http://develop.forums.netpad.net.cn/attachment.php?aid=19649&k=9d07d43d1db70ec5eb298e928679a395&t=1750241456&sid=HzJBH6

作者: 柳烟 时间: 2013-3-23 23:53

牛集落入自身陷阱：

图片附件: Fractal1.png (2013-3-23 23:53, 152.69 KB) / 下载次数 1721
http://develop.forums.netpad.net.cn/attachment.php?aid=19650&k=0a5bd924b2ac99c6d68654488a9568e0&t=1750241456&sid=HzJBH6

作者: xiaongxp 时间: 2013-3-24 16:17

4# 柳烟
用相似变换也可以达到这一效果

图片附件: 分形叠加(Mandelbrot & Sierpinski).jpg (2013-3-24 16:17, 48.13 KB) / 下载次数 1262
http://develop.forums.netpad.net.cn/attachment.php?aid=19651&k=af1c44976310824e5a76e09d80150dc1&t=1750241456&sid=HzJBH6

附件: 分形叠加(Mandelbrot & Sierpinski).gsp (2013-3-24 16:17, 21.8 KB) / 下载次数 1710
http://develop.forums.netpad.net.cn/attachment.php?aid=19652&k=894f0327189b382a9886732c00d4cb40&t=1750241456&sid=HzJBH6

作者: 柳烟 时间: 2013-3-24 16:25

13# xiaongxp
很是漂亮！学习。

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