Board logo

标题: 几何画板分形的mathcad之实现 [打印本页]

作者: lnszdzg    时间: 2014-4-27 19:58     标题: 几何画板分形的mathcad之实现

应mjj_ljh老师和柳烟老师之约,开设这个话题。
其实,我做分形是从巷老师的《画板分形39例》入门的,后来在画板论坛里的学习中逐渐开始做的。
感谢论坛里各位老师的无私奉献和热情指导!
既然开了话题,咱就不绕弯子了,下从M集说起吧。
作者: lnszdzg    时间: 2014-4-27 20:07

捕获.PNG
New.gif

图片附件: 捕获.PNG (2014-4-27 20:07, 13.8 KB) / 下载次数 2701
http://develop.forums.netpad.net.cn/attachment.php?aid=21446&k=5f0ea2b454aa2bb53cd111d854c29511&t=1750250889&sid=InA5CQ



图片附件: New.gif (2014-7-9 09:47, 6.05 KB) / 下载次数 2699
http://develop.forums.netpad.net.cn/attachment.php?aid=21447&k=03454a6bcc92084f381de4b8a18c3025&t=1750250889&sid=InA5CQ


作者: xiaongxp    时间: 2014-4-27 20:12

翘首以待,期盼已久,一定跟你好好学习,谢谢。
请推荐一个免费的mathcad下载链接,我自己下的就是不能正常运行。
作者: lnszdzg    时间: 2014-4-27 20:19

对程序的解释:
m=400,这是像素的个数,即总像素为400*400个,其中x轴方向和y轴方向个为400个点(画板中不予考虑)
stpx和stpy分别为x和y轴方向上的步长(画板中不予考虑);
下面的两个for循环是对c的取值实部从-2.51到1.5,虚部从-2到+2按照步长取值(扫描);
第三个for循环就是迭代了(相当于UF中的loop句)
break if |z|>=4,当z的模大于4跳出循环,z=z^2+c是迭代
Nj,k=iter,就迭代次数赋值给矩阵N,(直接以iter作图,当然后面会介绍其他作图方法)
最后几句输出矩阵N
作图如上(黑白的)
作者: lnszdzg    时间: 2014-4-27 20:23

3# xiaongxp

http://pan.baidu.com/s/1mgI6VAo
这是最新版的3.0(我没安装,不知是否好用)
过会我传15给您
作者: lnszdzg    时间: 2014-4-27 20:29

翘首以待,期盼已久,一定跟你好好学习,谢谢。
请推荐一个免费的mathcad下载链接,我自己下的就是不能正常运行。
xiaongxp 发表于 2014-4-27 20:12
PTC_MathCAD_15_M010_Win32...体中文版-精诚网WWW.CAXIT.COM.iso
115网盘礼包码:5lber9510pgi
http://115.com/lb/5lber9510pgi
作者: xiaongxp    时间: 2014-4-27 20:35

6# lnszdzg
谢谢。
作者: 柳烟    时间: 2014-4-27 20:38

非常感谢lnszdzg给我们这些菜鸟带来的福音,一定抓紧时间好好学习,谢谢.
作者: lnszdzg    时间: 2014-4-27 20:39

好,下面开始简单的上色(从画板分形39例学的)
捕获2.PNG
New.gif

图片附件: 捕获2.PNG (2014-4-27 20:39, 17.3 KB) / 下载次数 2706
http://develop.forums.netpad.net.cn/attachment.php?aid=21448&k=852184c6812891df782f3981adbd81fe&t=1750250889&sid=InA5CQ



图片附件: New.gif (2014-7-9 09:48, 11.05 KB) / 下载次数 2730
http://develop.forums.netpad.net.cn/attachment.php?aid=21449&k=33725a60f7619a9f2e8862d34db32008&t=1750250889&sid=InA5CQ


作者: xiaongxp    时间: 2014-4-27 20:43

115下载太慢了,急死人的!
2#的M集横向拉长了,是制图问题还是纵横像素比设置问题?
作者: lnszdzg    时间: 2014-4-27 20:45

捕4获.PNG
与程序1比较,主要区别在此。
s=0.05(iter-log(0.1+|ln|z||))着色参数(跟巷老师学的,其实我是“拿来主义”啊)
Rk,j Gk,j Bk,j三个着色矩阵对应于R G B

图片附件: 捕4获.PNG (2014-4-27 20:45, 5.17 KB) / 下载次数 2149
http://develop.forums.netpad.net.cn/attachment.php?aid=21450&k=da26d531fe525a3d438c8f0744eaa371&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-4-27 20:46

捕4获.PNG


与程序1比较,主要区别在此。
s=0.05(iter-log(0.1+|ln|z||))着色参数(跟巷老师学的,其实我是“拿来主义”啊)
Rk,j Gk,j Bk,j三个着色矩阵对应于R G B

图片附件: 捕4获.PNG (2014-4-27 20:46, 5.17 KB) / 下载次数 2147
http://develop.forums.netpad.net.cn/attachment.php?aid=21451&k=47cfac31e51cc4f8fcea11916dbd2c2e&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-4-27 20:48

10# xiaongxp

在Mathcad中图形可以随意拉长和压缩,这是我在截图是没有注意图形的比例造成的。
作者: lnszdzg    时间: 2014-4-27 20:51

115下载太慢了,急死人的!
2#的M集横向拉长了,是制图问题还是纵横像素比设置问题?
xiaongxp 发表于 2014-4-27 20:43
http://vdisk.weibo.com/s/dwcr_yUmd1U7E
作者: lnszdzg    时间: 2014-4-27 21:01

好了,我们对M集再来点别的着色方法
捕获5.PNG

捕获6.jpg

图片附件: 捕获5.PNG (2014-4-27 21:01, 18.64 KB) / 下载次数 2025
http://develop.forums.netpad.net.cn/attachment.php?aid=21452&k=856532b38e06400912c7418ebbfb8703&t=1750250889&sid=InA5CQ



图片附件: 捕获6.jpg (2014-4-28 18:35, 26.42 KB) / 下载次数 2012
http://develop.forums.netpad.net.cn/attachment.php?aid=21453&k=b772c76a545aec894f66ebb8ce1537ea&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-4-27 21:11

捕获7.jpg

伪3D效果

图片附件: 捕获7.jpg (2014-4-27 21:11, 26.31 KB) / 下载次数 2041
http://develop.forums.netpad.net.cn/attachment.php?aid=21454&k=7c7b7cc25ce8d1f6a75fa531236069c5&t=1750250889&sid=InA5CQ


作者: 柳烟    时间: 2014-4-27 23:28

未命名.JPG
从前面115盘下了软件,点setup,没反应,安不上软件,不知怎回事.

图片附件: 未命名.JPG (2014-4-27 23:28, 31.23 KB) / 下载次数 2004
http://develop.forums.netpad.net.cn/attachment.php?aid=21459&k=8f48b66059792ac08b7332788b3fab49&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-4-28 08:36

21459
从前面115盘下了软件,点setup,没反应,安不上软件,不知怎回事.
柳烟 发表于 2014-4-27 23:28
看看这个:
http://www.cr173.com/soft/23176.html
作者: 柳烟    时间: 2014-4-28 21:43

这个软件不好安,安装途中,安这个软件安那个软件,怎么这么费事?昨晚试安,今晚又安,还没摆平.
作者: lnszdzg    时间: 2014-4-28 21:47

这个软件不好安,安装途中,安这个软件安那个软件,怎么这么费事?昨晚试安,今晚又安,还没摆平.
柳烟 发表于 2014-4-28 21:43
http://tieba.baidu.com/p/2945249354
这里有介绍,您看看
作者: 柳烟    时间: 2014-4-28 22:11

20# lnszdzg
非常感谢.
作者: mjj_ljh    时间: 2014-4-28 22:35

这个软件是够难装的,明天有时间接着装!好事多麿啊!
Inszdzg老师贵姓,你的英文名太难记也太难打了,一不留神就打错了。
作者: lnszdzg    时间: 2014-4-29 08:55

我有一个面安装版的,就是版本有点低(2001),试用起来也不错
将文件解压后,
这样 New.gif
再这样 捕获77.PNG
这是网址http://pan.baidu.com/s/1pJHKXdD

图片附件: New.gif (2014-7-9 09:50, 5.68 KB) / 下载次数 1733
http://develop.forums.netpad.net.cn/attachment.php?aid=21465&k=125586ca5b2d67a0a011ff502c95b8c8&t=1750250889&sid=InA5CQ



图片附件: 捕获77.PNG (2014-4-29 08:55, 8.6 KB) / 下载次数 1693
http://develop.forums.netpad.net.cn/attachment.php?aid=21466&k=14550db898c871ae2208b8e0945874b6&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-4-29 09:00

将mathcad1.exe改名为mathcad.exe(必须的)
找到regtool.exe单击,再单击注册组件,然后关闭。运行mcad或mathcad即可。
虽然版本低,但运行挺不错的。尤其是符号运算的核心和maple是一样的。
作者: lnszdzg    时间: 2014-4-29 09:02

这个软件是够难装的,明天有时间接着装!好事多麿啊!
Inszdzg老师贵姓,你的英文名太难记也太难打了,一不留神就打错了。
mjj_ljh 发表于 2014-4-28 22:35
哈哈,我们的英文名字彼此彼此啦!
我姓杜,你们可以叫我“老杜”啊!
作者: lnszdzg    时间: 2014-4-29 09:22

继续M集。
下面我们对M集进行装饰和美化。
从几位老师的作品和网上现有的资料,以及UF中的着色来看,对于分形的装饰和美化,无非有这样几种:
1、平滑过渡:即在对分形着色时,使各种颜色之间的过渡柔和、淡化,不出现明显的梯度;
2、陷阱点缀;通过各种陷阱的试用,使分形表现出另外的一种美;
3、夸张渲染:即通过颜色的明显对比,产生强烈的视觉冲击。
这里我们先说平滑。
对于平滑算法,几位老师的不尽相同,我摘录几个。
1)UF中的平滑:index=0.05*(#numiter+il*ip-il*log(log(|z|)))
捕获8.PNG

New.gif

图片附件: 捕获8.PNG (2014-4-29 09:22, 21.83 KB) / 下载次数 1718
http://develop.forums.netpad.net.cn/attachment.php?aid=21468&k=effcb0946703ed44c4b65f8c70271b3b&t=1750250889&sid=InA5CQ



图片附件: New.gif (2014-7-9 09:51, 15.86 KB) / 下载次数 1708
http://develop.forums.netpad.net.cn/attachment.php?aid=21469&k=71c85a086e0bedfab7aed2395dbb8472&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-4-29 10:09

2)、侯老师的插值平滑法http://blog.csdn.net/housisong/article/details/6174894
捕获10.PNG

New.gif

图片附件: 捕获10.PNG (2014-4-29 10:09, 28.33 KB) / 下载次数 1726
http://develop.forums.netpad.net.cn/attachment.php?aid=21470&k=79ba2dbff2f6ac01c794bb273310b83b&t=1750250889&sid=InA5CQ



图片附件: New.gif (2014-7-9 09:52, 18.29 KB) / 下载次数 1730
http://develop.forums.netpad.net.cn/attachment.php?aid=21471&k=b0ca5fa14dd5a9223ae586901f0d2840&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-4-29 10:21

侯老师的新插值公式http://blog.csdn.net/housisong/article/details/6174894
捕获12.PNG

New.gif

图片附件: 捕获12.PNG (2014-4-29 10:21, 32.29 KB) / 下载次数 1696
http://develop.forums.netpad.net.cn/attachment.php?aid=21472&k=2b941026857d5faa3ee93411190f0e2f&t=1750250889&sid=InA5CQ



图片附件: New.gif (2014-7-9 09:52, 16.25 KB) / 下载次数 1694
http://develop.forums.netpad.net.cn/attachment.php?aid=21473&k=bc8d1df87309122b38ed8f202946810f&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-4-29 10:29

画板分形中的平滑方法我就不再班门抡斧了,几位老师都是我的导师。
对于新手,可以参考http://www.inrm3d.cn/viewthread.php?tid=4738&page=1#pid43153
作者: lnszdzg    时间: 2014-4-29 10:44

4、三角不等式平滑法(3、就是画板中的方法)
New.jpg

New.jpg

图片附件: New.jpg (2014-7-9 09:53, 8.89 KB) / 下载次数 1625
http://develop.forums.netpad.net.cn/attachment.php?aid=21474&k=cf0094d54d12ddd2b8355e6df240238e&t=1750250889&sid=InA5CQ



图片附件: New.jpg (2014-7-9 09:53, 14.41 KB) / 下载次数 1616
http://develop.forums.netpad.net.cn/attachment.php?aid=21475&k=4d2df64edc2b10a563751335cc4bfc6e&t=1750250889&sid=InA5CQ


作者: 柳烟    时间: 2014-4-29 12:19

以上几张图片十分漂亮!
作者: 榕坚    时间: 2014-4-29 19:20

不知道在mathcad中扫描放大图的速度如何,我在mahematica中放大一个迭代次数只有500的图把电脑的内存一下就上升到90%,电脑就像酒喝醉了一样。

图片附件: Untitled-1.jpg (2014-4-29 21:47, 36.83 KB) / 下载次数 1714
http://develop.forums.netpad.net.cn/attachment.php?aid=21477&k=349901683e1a0ede59ac920de8b82aeb&t=1750250889&sid=InA5CQ


作者: xiaongxp    时间: 2014-4-29 20:12

换了好几个mathcad安装文件,都不能完成安装。
作者: lnszdzg    时间: 2014-4-29 21:54

换了好几个mathcad安装文件,都不能完成安装。
xiaongxp 发表于 2014-4-29 20:12
出现什么样的提示?
作者: lnszdzg    时间: 2014-5-1 13:01

将z的幅角加进来,可以修饰M集的边界

New.jpg

图片附件: New.jpg (2014-7-9 09:54, 31.35 KB) / 下载次数 1839
http://develop.forums.netpad.net.cn/attachment.php?aid=21492&k=55dd60a313352dedd7b3a80da2bba104&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-5-1 13:48

将三角不等式着色法中的指数加大为
sum = sum + ((cabs(#z) - lowbound) / (az2+ac - lowbound))^@apower,
就有这种效果
New.jpg

图片附件: New.jpg (2014-7-9 09:55, 10.17 KB) / 下载次数 1853
http://develop.forums.netpad.net.cn/attachment.php?aid=21493&k=127f59322a95e175999f87fcecfbeda9&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-5-1 17:07

在26#程序中增加以下语句:
;type分形类型(type=0M集,type=1Julia集)
lx=-2.501+0.5type  ;中心平移
pixcel=lx+(j-1)stpx+(-1.5+(k-1)stpy)
z=type*pixcel ;z的初值
c=(1-type)*pixcel+type*seed ;c的初值
可以将M和Julia集合并在同一程序中。

将Trapped=0修改为Trapped*s=0;(s=0为先下后上,s=1为后来居上)
这样一个程序就可以包含M和J集,并能将先下后上和后来居上也合并了。

New.jpg
New.jpg
New.jpg
New.jpg

图片附件: New.jpg (2014-7-9 09:55, 24.51 KB) / 下载次数 1877
http://develop.forums.netpad.net.cn/attachment.php?aid=21496&k=9299da5a05313322f199a9c769fcb469&t=1750250889&sid=InA5CQ



图片附件: New.jpg (2014-7-9 09:56, 20.86 KB) / 下载次数 1850
http://develop.forums.netpad.net.cn/attachment.php?aid=21497&k=c3b057471ad912a85719ab817ee55dbc&t=1750250889&sid=InA5CQ



图片附件: New.jpg (2014-7-9 09:56, 18.42 KB) / 下载次数 1907
http://develop.forums.netpad.net.cn/attachment.php?aid=21498&k=798ff1c1ced4b908a943eb4ccb39bec8&t=1750250889&sid=InA5CQ



图片附件: New.jpg (2014-7-9 09:57, 22.11 KB) / 下载次数 1879
http://develop.forums.netpad.net.cn/attachment.php?aid=21499&k=3ed6d4526d3c892be79308f471b83d0a&t=1750250889&sid=InA5CQ


作者: 柳烟    时间: 2014-5-1 19:11

35# lnszdzg
这个分形有特色。
作者: lnszdzg    时间: 2014-5-2 11:12

2、以下继续聊一聊陷阱
几位老师分别给出了点陷阱、圆陷阱、环陷阱等等,其实只要制造陷阱的思路清楚了,陷阱问题也就顺理成章了。
虽然UF中给出了各种优美的陷阱,但我还是非常欣赏几位老师在画板中制造陷阱的方法(尤其是IF,if的构造特别奇妙)
这是UF中的一个陷阱

New.jpg


New.jpg

图片附件: New.jpg (2014-7-9 09:57, 32.54 KB) / 下载次数 1864
http://develop.forums.netpad.net.cn/attachment.php?aid=21504&k=61b41ff91c70774a66928a545f6a4502&t=1750250889&sid=InA5CQ



图片附件: New.jpg (2014-7-9 09:58, 34.63 KB) / 下载次数 1841
http://develop.forums.netpad.net.cn/attachment.php?aid=21505&k=e437dca998e3610b112da879b75d2489&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-5-2 11:15

上图是利用Mathcad做的,
UF代码
mt-beads { ; Mark Townsend, 13 July 1999
init:
  trapped_z = 0
  int iter = 0
  int trapped_iter = 0
  float distance = 0
  bool trapped = false
  
  if @no_of_iters == 0
    int max_iter = #maxiter - @skip
  else
    int max_iter = @skip + @no_of_iters
  endif  
  
loop:
  if !trapped || !@trapping && iter >= @skip && iter < max_iter
    if (abs(real(#z)) < abs(cos(imag(#z))))
      trapped = true
      trapped_iter = iter
      trapped_z = #z
      distance = real(abs((abs(real(#z)) - abs(cos(imag(#z))))))
    elseif (abs(imag(#z)) < abs(sin(real(#z))))
      trapped = true
      trapped_iter = iter
      trapped_z = #z
      distance = real(abs((abs(imag(#z)) - abs(sin(real(#z))))))
    endif
  endif
  iter = iter + 1
final:
  if trapped
    if @color == 0
      #index = distance
    elseif @color == 1
      #index = cabs(trapped_z)
    elseif @color == 2
      #index = abs(real(trapped_z))
    elseif @color == 3
      #index = abs(imag(trapped_z))
    elseif @color == 4
      float angle = atan2(trapped_z)
      if angle < 0
        angle = angle + 2 * #pi
      endif
      angle = 1 / (2 * #pi) * angle
      #index = angle
    elseif @color == 5
      #index = 0.01 * trapped_iter
    endif  
  else
    #solid = true
  endif   
default:
  title = "Beads"
  param @color
    caption = "Color mode"
    enum = "Distance" "Magnitude" "Real" "Imag" "Angle" "Iteration"
  endparam  
  
  param skip
    caption = "Iterations to skip"
    default = 0
  endparam  
   
  param no_of_iters
    caption = "No. of iterations"
    hint = "This is the number of iterations to monitor the orbit."
    default = 0
    min = 0
  endparam  
  
  param trapping
   caption = "Bailout on first"
   default = false
  endparam  

}

Mathcad代码只能用图片,太长了上传很不便啊
作者: lnszdzg    时间: 2014-5-2 11:28

揭开前三层,就很有特色了
New.jpg

图片附件: New.jpg (2014-7-9 09:59, 41.5 KB) / 下载次数 1551
http://develop.forums.netpad.net.cn/attachment.php?aid=21506&k=2a531d6b3bab534d6952afd7b2a1e9b0&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-5-2 11:30

New.jpg

图片附件: New.jpg (2014-7-9 09:59, 34.91 KB) / 下载次数 1568
http://develop.forums.netpad.net.cn/attachment.php?aid=21507&k=1f02eb6978309fe816f5f18fe16bb69d&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-5-2 11:56

还是蓝白调好看
New.jpg

图片附件: New.jpg (2014-7-9 10:00, 35.08 KB) / 下载次数 1654
http://develop.forums.netpad.net.cn/attachment.php?aid=21509&k=6e1313aafdf6eeadbe0112bc448c1a2f&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-5-2 19:21

不知道在mathcad中扫描放大图的速度如何,我在mahematica中放大一个迭代次数只有500的图把电脑的内存一下就上升到90%,电脑就像酒喝醉了一样。
榕坚 发表于 2014-4-29 19:20
跟新作图没有区别,因为作图机理不同。
Mathcad作图都必须转化为矩阵,所以只要矩阵的大小一样,就无所谓了。
作者: 榕坚    时间: 2014-5-2 21:08

44# lnszdzg


主要是放大图需要的迭代次数增加了,mathematica也是计算后一次性成图。我做了那个图后就不敢再试了。
作者: 柳烟    时间: 2014-5-2 22:03

maple 15这软件,比起mathematica好装,好象也能造作分形,可是整不明白。杜老师讲解mathematica造作分形,晓畅明白,辛苦了。
作者: lnszdzg    时间: 2014-5-3 10:45

几位老师,Mathcad15的安装确实比较麻烦,如果有兴趣可以试试免安装版的

这是网址http://pan.baidu.com/s/1pJHKXdD
作者: lnszdzg    时间: 2014-5-3 11:52

New.jpg

三龟竞水

图片附件: New.jpg (2014-7-9 10:00, 18.52 KB) / 下载次数 1655
http://develop.forums.netpad.net.cn/attachment.php?aid=21517&k=da87f0a8a8c8ea39550485efa996d1db&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-5-3 12:09

New.jpg

图片附件: New.jpg (2014-7-9 10:01, 9.59 KB) / 下载次数 1657
http://develop.forums.netpad.net.cn/attachment.php?aid=21518&k=abf30f8e7c789ead30b10023a9039e55&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-5-3 12:25

New.jpg

图片附件: New.jpg (2014-7-9 10:01, 36.13 KB) / 下载次数 1643
http://develop.forums.netpad.net.cn/attachment.php?aid=21519&k=9b657f54a804f441cf67f4060c9ce677&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-5-6 10:04

New.jpg
UF中的一个陷阱,Mathcad实现了
UF代码
HypocycloidTrap {
  ; Paul Carlson
  ; Coloring Method extraced from "HNspirJ2".
  ; By Ken Childress July 1999.  
init:
  bool trapped = FALSE
  float range_num = 0
  float num_ranges = real(@RangesColors)
  float colors_in_range = imag(@RangesColors)
  float index_factor = (colors_in_range - 1) / @width
  float color_index = 0
  int iter = 0
  complex i = (0,1)
  complex k = 0
  complex angle = 0
  complex ka = 0
  complex rz = 0
  complex iz = 0
  complex v = 0
  float dist = 0

loop:
  IF !trapped && (iter >= @skip)
    angle = atan(imag(#z) / real(#z))
    k = 1.2 - @hfactor
    ka = k * angle / @hfactor
    rz = k * cos(angle) + @hfactor * cos(ka)
    iz = k * sin(angle) + @hfactor * sin(ka)
    v = rz + i * iz
    dist = abs(|#z| - |v|)
   
    IF (dist < @width && iter >= @skip)
      trapped = TRUE
      range_num = iter % num_ranges
      color_index = index_factor * dist +  range_num * colors_in_range
    ENDIF
   
  ENDIF

  iter = iter + 1

final:
  IF trapped
    #index = ((real(color_index)+1) % 256) / 256 + @randomness * real(#random)
  ELSE
    IF @solid
      #solid = TRUE
    ELSE
      #index = @back_color / 256
    ENDIF
  ENDIF
   
default:
  title = "Hypocycloid Trap"

  param width
    caption = "Stalk Width Factor"
    default = 0.03
    hint = "Controls stalk width."
  endparam

  param hfactor
    caption = "Hypocycloid Factor"
    default = 0.4
    hint = "Hypocycloid Factor."
  endparam
   
  param skip
    caption = "Iterations to skip"
    default = 0
    hint = "Iterations to skip."
  endparam

  param RangesColors
    caption = "Ranges and Colors"
    default = (2, 125)
    hint = "Number of Color Ranges and Number of Colors in each Range."
  endparam

  param back_color
    caption = "Background Color"
    hint = "This is the color used for orbits that aren't trapped. It has \
           no effect if the background is solid."
    default = 254
    min = 0
    max = 255
  endparam

  param solid
    caption = "Solid Background"
    hint = "If this is enabled orbits that aren't trapped become solid."
    default = false
  endparam
  
  param randomness
    caption = "Random Texture"
    default = 0.0
    hint = "This adds a random texture to the coloring."
  endparam
}

图片附件: New.jpg (2014-7-9 10:02, 21.1 KB) / 下载次数 1367
http://develop.forums.netpad.net.cn/attachment.php?aid=21553&k=eebe875a52903ac49cc07988729a1c2b&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-5-6 10:26

New.jpg

图片附件: New.jpg (2014-7-9 10:03, 18.29 KB) / 下载次数 1373
http://develop.forums.netpad.net.cn/attachment.php?aid=21554&k=7f92c8b29075730eccc9057a59f4328c&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-5-6 10:52

52# lnszdzg


点陷阱,用了HLS颜色
作者: lnszdzg    时间: 2014-5-6 14:16

New.jpg
New.jpg

图片附件: New.jpg (2014-7-9 10:03, 25.35 KB) / 下载次数 1369
http://develop.forums.netpad.net.cn/attachment.php?aid=21555&k=ff69358e5eb35a993067937f7720c2d1&t=1750250889&sid=InA5CQ



图片附件: New.jpg (2014-7-9 10:04, 25.77 KB) / 下载次数 1370
http://develop.forums.netpad.net.cn/attachment.php?aid=21556&k=bf94e96cfd825cef1efe3210bd548b14&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-5-6 14:19

UF代码

Rings_Coloring {
  ; Paul W. Carlson
  ; Variant #1 extracted from "RsNewtMset1".
  ; Variant #2 extracted from "R4Julia".
  ; By Ken Childress June 1999.
init:
  bool first = TRUE
  bool trapped = FALSE
  int iter = 0
  float HalfT = 0.5 * @T
  float Phi = #pi * 0.125
  float Ro = @Rm + HalfT
  float Py = @Rm * sin(Phi)
  float Px = @Rm * cos(Phi)
  float Dsqd = @Rm * @Rm + Ro * Ro - 2 * Ro * Px
  float color_index = 252
  float index_factor = @colors_in_range - 1
   
  ; Initialization for Variant #1
  complex k = (0.5, 0.5)
  float index_factor1 = (@colors_in_range - 1) / @Rm

loop:
  IF @Variant == 0
  
    IF (!trapped && (abs(cabs(#z) - @Rm) < HalfT) && \
      first == FALSE && iter >= @skip)
      
      trapped = TRUE
      float X = real(#z)
      float Y = imag(#z)
      float Xabs = abs(X)
      float Yabs = abs(Y)
      IF (Xabs >= Yabs)
        float WtoPsqd = (Xabs-Px)*(Xabs-Px) + (Yabs-Py)*(Yabs-Py)
      ELSE
        WtoPsqd = (Xabs-Py)*(Xabs-Py) + (Yabs-Px)*(Yabs-Px)
      ENDIF
      IF (X >= 0 && Y >= 0)
        IF (Xabs >= Yabs)
          float Segment = 0
        ELSE
          Segment = 1
        ENDIF
      ELSEIF (X < 0 && Y >= 0)
        IF (Xabs < Yabs)
          Segment = 2
        ELSE
          Segment = 3
        ENDIF
      ELSEIF (X < 0 && Y < 0)
        IF (Xabs >= Yabs)
          Segment = 4
        ELSE
          Segment = 5
        ENDIF
      ELSE
        IF (Xabs < Yabs)
          Segment = 6
        ELSE
          Segment = 7
        ENDIF
      ENDIF
      float Ratio = sqrt(WtoPsqd/Dsqd)
      ; Modulate Segment for less than eight ranges.
      float range_num = Segment % @num_ranges
      color_index = index_factor * Ratio + range_num * @colors_in_range
    ENDIF
    first = FALSE
   
  ELSEIF @Variant == 1
  
    IF !trapped && iter >= @skip
      float dist = abs(|#z - k| - 0.01)
      IF dist < @Rm
        trapped = TRUE
        float range_num = iter % @num_ranges
        color_index = index_factor1 * dist + range_num * @colors_in_range
      ENDIF
    ENDIF

  ENDIF
  
  iter = iter + 1

final:
  if trapped
    ; Adjust color index to gradient range.
    #index = (color_index + 1) % 256 / 256 + @randomness * real(#random)
  else
    if @solid
      #solid = true
    else
    ; For the purists we don't have to be solid.
      #index = @back_color / 256
    endif
  endif

default:
  title = "Ring Coloring"
  
  param Variant
    caption = "Ring Variant"
    enum = "Variant #1" "Variant #2"
    default = 0
    hint = "The Ring variant to use."
  endparam

  param Rm
    caption = "Ring radius"
    default = 1.0
    hint = "The radius to the midline of the ring."
  endparam

  param T
    caption = "Ring thickness"
    default = 0.2
    hint = "The thickness of the ring. This is not used for Variant #2."
  endparam
  
  param skip
    caption = "Iterations to Skip"
    default = 0
    hint = "The number of iterations to skip.  This is not used for Normal."
  endparam

  param num_ranges
    Caption = "Color ranges"
    hint = "This is the number of separate color ranges."
    default = 8
  endparam

  param colors_in_range
    Caption = "Colors in range"
    hint = "This is the number of colors in each color range."
    default = 30
  endparam

  param back_color
    caption = "Background Color"
    hint = "This is the color used for orbits that aren't trapped. It has \
           no effect if the background is solid."
    default = 254
  endparam

  param solid
    caption = "Solid Background"
    hint = "If this is enabled orbits that aren't trapped become solid."
    default = false
  endparam
  
  param randomness
    caption = "Random Texture"
    default = 0.0
    hint = "This adds a random texture to the coloring."
  endparam
}
作者: lnszdzg    时间: 2014-5-8 21:14

New.jpg

一园一环陷阱

图片附件: New.jpg (2014-7-9 10:04, 31 KB) / 下载次数 1266
http://develop.forums.netpad.net.cn/attachment.php?aid=21563&k=e8d0506bd5cf4866308f8a7161edb868&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-5-8 21:22

New.jpg

同臂同色

图片附件: New.jpg (2014-7-9 10:05, 24.38 KB) / 下载次数 1289
http://develop.forums.netpad.net.cn/attachment.php?aid=21564&k=2ce11abe5898a11ee76a963926497f07&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-5-8 21:40

New.jpg

图片附件: New.jpg (2014-7-9 10:05, 22.95 KB) / 下载次数 1260
http://develop.forums.netpad.net.cn/attachment.php?aid=21565&k=54e994fb557c5ce31db939a905a21a80&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-5-8 22:02

New.jpg

图片附件: New.jpg (2014-7-9 10:05, 30.27 KB) / 下载次数 1289
http://develop.forums.netpad.net.cn/attachment.php?aid=21567&k=dd1d41a8d23be8f7fd7baf7c7dc4ce52&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-5-11 18:11

New.jpg
New.jpg
New.jpg
New.jpg

图片附件: New.jpg (2014-7-9 10:07, 28.82 KB) / 下载次数 1259
http://develop.forums.netpad.net.cn/attachment.php?aid=21585&k=550df23bd3cf42e326e3827dfe3da640&t=1750250889&sid=InA5CQ



图片附件: New.jpg (2014-7-9 10:07, 31.51 KB) / 下载次数 1281
http://develop.forums.netpad.net.cn/attachment.php?aid=21586&k=75ed5be86023f91286efb06608aec56d&t=1750250889&sid=InA5CQ



图片附件: New.jpg (2014-7-9 10:08, 44.14 KB) / 下载次数 1343
http://develop.forums.netpad.net.cn/attachment.php?aid=21587&k=575503cbbbf199bb24bd8ed82b1d7969&t=1750250889&sid=InA5CQ



图片附件: New.jpg (2014-7-9 10:08, 28.16 KB) / 下载次数 1290
http://develop.forums.netpad.net.cn/attachment.php?aid=21588&k=9f5b39c7db1f0dfbeb3902284a83e77e&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-5-13 13:03

New.jpg

图片附件: New.jpg (2014-7-9 10:09, 43.13 KB) / 下载次数 1239
http://develop.forums.netpad.net.cn/attachment.php?aid=21598&k=20aa9c283f90948386f6c11cb7d30391&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-5-20 16:45

New.jpg

图片附件: New.jpg (2014-7-9 10:09, 6.59 KB) / 下载次数 1189
http://develop.forums.netpad.net.cn/attachment.php?aid=21647&k=95a0b39298fa67d666539e3d5fcb7e9c&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-5-20 17:02

New.jpg

图片附件: New.jpg (2014-7-9 10:10, 23.68 KB) / 下载次数 1212
http://develop.forums.netpad.net.cn/attachment.php?aid=21648&k=b629eca47f4f59dc8865ca04f98826b1&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-5-20 17:32

New.jpg

图片附件: New.jpg (2014-7-9 10:10, 23.88 KB) / 下载次数 1176
http://develop.forums.netpad.net.cn/attachment.php?aid=21649&k=03a03750d6d287e1871d57b226cf9894&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-5-20 17:34

New.jpg

图片附件: New.jpg (2014-7-9 10:10, 36.49 KB) / 下载次数 1282
http://develop.forums.netpad.net.cn/attachment.php?aid=21650&k=3ca6efa736c5a798d8a3706146a13b5d&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-5-21 18:07

New.jpg

图片附件: New.jpg (2014-7-9 10:11, 25.83 KB) / 下载次数 1244
http://develop.forums.netpad.net.cn/attachment.php?aid=21657&k=59cbf7198be90f2759302c78fa273956&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-5-21 18:09

New.jpg

蝴蝶结陷阱

图片附件: New.jpg (2014-7-9 10:11, 30.96 KB) / 下载次数 1252
http://develop.forums.netpad.net.cn/attachment.php?aid=21658&k=ae8a950daaae343b79012632e7aa9e69&t=1750250889&sid=InA5CQ


作者: guxiaoping    时间: 2014-5-21 18:31

67# lnszdzg
用什么曲线作陷阱?
作者: lnszdzg    时间: 2014-5-21 18:47

New.jpg

图片附件: New.jpg (2014-7-9 10:12, 57.46 KB) / 下载次数 1439
http://develop.forums.netpad.net.cn/attachment.php?aid=21659&k=b78f1000dec35caa5ab386647e2296ce&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-5-21 19:20

68# guxiaoping

像蝴蝶结哈,方程式||z|-|sin(Re(z))||
作者: lnszdzg    时间: 2014-5-21 19:39

捕获81.jpg

曲线:|cos(z)-sin(Re(z))|

图片附件: 捕获81.jpg (2014-5-22 13:36, 44.32 KB) / 下载次数 1251
http://develop.forums.netpad.net.cn/attachment.php?aid=21660&k=62a764ee016b895c5ad0afa5bd2b55d3&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-5-21 19:45

曲线:
捕获83.jpg

图片附件: 捕获83.jpg (2014-5-22 13:39, 42.64 KB) / 下载次数 1215
http://develop.forums.netpad.net.cn/attachment.php?aid=21661&k=4ead164a71c65d37e9f554dfea6af83d&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-5-21 20:00

捕获84.jpg

图片附件: 捕获84.jpg (2014-5-22 13:40, 52.95 KB) / 下载次数 1364
http://develop.forums.netpad.net.cn/attachment.php?aid=21663&k=12bbfbbd404ed42d6c39d82385ae3f9d&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-5-21 20:18

捕获85.jpg

图片附件: 捕获85.jpg (2014-5-22 13:41, 42.67 KB) / 下载次数 1241
http://develop.forums.netpad.net.cn/attachment.php?aid=21664&k=9ca5484bcc767bc563a17bd9a6c34883&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-5-21 20:23

捕获86.jpg

图片附件: 捕获86.jpg (2014-5-22 13:42, 42.65 KB) / 下载次数 1235
http://develop.forums.netpad.net.cn/attachment.php?aid=21665&k=54c71b79902f76325859823aa69dc504&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-5-22 11:05

捕获91.jpg

香蕉啊:Re(z)^(2/3)+Im(z)^2

图片附件: 捕获91.jpg (2014-5-22 13:46, 27.64 KB) / 下载次数 1056
http://develop.forums.netpad.net.cn/attachment.php?aid=21666&k=5542670ff4d0428806098c5b3da27f45&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-5-22 11:13

New.jpg

图片附件: New.jpg (2014-7-9 10:12, 24.81 KB) / 下载次数 1068
http://develop.forums.netpad.net.cn/attachment.php?aid=21667&k=4556b37698cb2f0103f18d6f4f7d553a&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-5-22 14:28

DEM算法下的陷阱
New.jpg

图片附件: New.jpg (2014-7-9 10:13, 33.08 KB) / 下载次数 1079
http://develop.forums.netpad.net.cn/attachment.php?aid=21668&k=568c0947bbff2d1144aebbc2aade9482&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-5-23 17:00

New.jpg

图片附件: New.jpg (2014-7-9 10:13, 49.82 KB) / 下载次数 1080
http://develop.forums.netpad.net.cn/attachment.php?aid=21671&k=d0600c0862f6dcc965604ef5fc2001e4&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-5-23 17:02

New.jpg

图片附件: New.jpg (2014-7-9 10:13, 48.37 KB) / 下载次数 1073
http://develop.forums.netpad.net.cn/attachment.php?aid=21672&k=f4eebf7514a4ff54bd69eaeac953a625&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-5-23 17:03

New.jpg

图片附件: New.jpg (2014-7-9 10:14, 29.07 KB) / 下载次数 1219
http://develop.forums.netpad.net.cn/attachment.php?aid=21673&k=a01f386fac9267e1c14c98a647fecee5&t=1750250889&sid=InA5CQ


作者: 柳烟    时间: 2014-5-23 20:54

81# lnszdzg
这幅不错,啥陷阱?
作者: xiaongxp    时间: 2014-5-23 21:27

81# lnszdzg
双扭线?
作者: lnszdzg    时间: 2014-5-23 22:27

四叶玫瑰线
作者: lnszdzg    时间: 2014-5-24 11:21

捕获98.jpg
五叶玫瑰

图片附件: 捕获98.jpg (2014-5-24 11:30, 31.19 KB) / 下载次数 1189
http://develop.forums.netpad.net.cn/attachment.php?aid=21679&k=39af2d8ac9c2cddfda736de60c6a0191&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-5-24 11:24

方程与四叶玫瑰不同,这次采用了玫瑰线的统一表达式,简单多了
argz=arg(z)
rn=R*|sin(n/2*argz)|
|z|=rn^2
这里n是玫瑰花瓣数量,R是画板的大小
作者: lnszdzg    时间: 2014-5-24 16:15

New.jpg
桃花朵朵开

图片附件: New.jpg (2014-7-9 10:14, 25.93 KB) / 下载次数 1200
http://develop.forums.netpad.net.cn/attachment.php?aid=21684&k=422ae2586f9fd5c46327cb8b70d02193&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-5-24 17:36

New.jpg

图片附件: New.jpg (2014-7-9 10:15, 27.63 KB) / 下载次数 1212
http://develop.forums.netpad.net.cn/attachment.php?aid=21685&k=f1edff04bb780ae4dae73e656eea9752&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-5-25 21:22

喇叭花与桃花
捕获119.PNG
New.jpg

图片附件: 捕获119.PNG (2014-5-25 21:22, 25.08 KB) / 下载次数 1237
http://develop.forums.netpad.net.cn/attachment.php?aid=21697&k=0c75972570688ceb915a11484af8778a&t=1750250889&sid=InA5CQ



图片附件: New.jpg (2014-7-9 10:15, 21.16 KB) / 下载次数 1197
http://develop.forums.netpad.net.cn/attachment.php?aid=21698&k=d997b19c2f0f3992a403a9fed6bf0780&t=1750250889&sid=InA5CQ


作者: xiaongxp    时间: 2014-5-25 22:00

这些真花漂亮!
作者: lnszdzg    时间: 2014-5-26 16:06

New.jpg

五星

图片附件: New.jpg (2014-7-9 10:16, 16.76 KB) / 下载次数 830
http://develop.forums.netpad.net.cn/attachment.php?aid=21715&k=11170f227d5d20a0d9b272a86912932f&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-5-26 16:39

New.jpg

图片附件: New.jpg (2014-7-9 10:16, 19.17 KB) / 下载次数 835
http://develop.forums.netpad.net.cn/attachment.php?aid=21716&k=8b6aadba649efd2acbc0020dffa74e25&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-5-26 16:44

New.jpg

图片附件: New.jpg (2014-7-9 10:16, 19.93 KB) / 下载次数 822
http://develop.forums.netpad.net.cn/attachment.php?aid=21717&k=6232acc0307fe15b3cd1ccc77770f66a&t=1750250889&sid=InA5CQ


作者: 柳烟    时间: 2014-5-26 20:28

New.jpg

图片附件: New.jpg (2014-7-9 10:17, 23.31 KB) / 下载次数 989
http://develop.forums.netpad.net.cn/attachment.php?aid=21721&k=547a1558a12367e785348a12c058c9be&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-5-26 22:03

94# 柳烟

柳老师这个很漂亮。
来个带花心的

New.jpg

图片附件: New.jpg (2014-7-9 10:17, 18.08 KB) / 下载次数 854
http://develop.forums.netpad.net.cn/attachment.php?aid=21722&k=bee3afbf7e2793100839cea2b0a99b36&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-5-27 09:58

New.jpg

图片附件: New.jpg (2014-7-9 10:18, 31.98 KB) / 下载次数 844
http://develop.forums.netpad.net.cn/attachment.php?aid=21724&k=139ef57b5297df96d8da0f7e9232d5f4&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-5-27 10:36

心形线下的M集
New.jpg

掏空肚子

New.jpg

New.jpg

图片附件: New.jpg (2014-7-9 10:18, 21.3 KB) / 下载次数 815
http://develop.forums.netpad.net.cn/attachment.php?aid=21725&k=5a7b19ab6929600a0067456fc0d9719c&t=1750250889&sid=InA5CQ



图片附件: New.jpg (2014-7-9 10:19, 23.62 KB) / 下载次数 842
http://develop.forums.netpad.net.cn/attachment.php?aid=21726&k=44b63bb7177d89221f36f0d1d2187981&t=1750250889&sid=InA5CQ



图片附件: New.jpg (2014-7-9 10:19, 26.54 KB) / 下载次数 820
http://develop.forums.netpad.net.cn/attachment.php?aid=21727&k=5a771426d973c95492334e2401841592&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-5-27 10:44

New.jpg
双纽线
New.jpg
三叶玫瑰

图片附件: New.jpg (2014-7-9 10:20, 27.18 KB) / 下载次数 849
http://develop.forums.netpad.net.cn/attachment.php?aid=21728&k=63c2446283e9d515044cabc0a9dac24e&t=1750250889&sid=InA5CQ



图片附件: New.jpg (2014-7-9 10:20, 28.88 KB) / 下载次数 831
http://develop.forums.netpad.net.cn/attachment.php?aid=21729&k=445e8b0f621a2cb27784b361c2a38ddc&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-5-27 21:29

棋盘背景
捕获138.PNG

图片附件: 捕获138.PNG (2014-5-27 21:29, 49.69 KB) / 下载次数 824
http://develop.forums.netpad.net.cn/attachment.php?aid=21739&k=b09622d0651951b4abe1d0f6d9c74531&t=1750250889&sid=InA5CQ


作者: lnszdzg    时间: 2014-5-28 16:18

DEM算法下的J集,对像素的要求太高了,速度很慢啊。
捕获1392.jpg

图片附件: 捕获1392.jpg (2014-5-28 16:18, 27.78 KB) / 下载次数 982
http://develop.forums.netpad.net.cn/attachment.php?aid=21742&k=eda9bf26dedeceb0f2be323cf6df295a&t=1750250889&sid=InA5CQ






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