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

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

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

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

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

图片附件: 捕获.PNG (2014-4-27 20:07, 13.8 KB) / 下载次数 2289
http://www.inrm3d.cn/attachment.php?aid=21446&k=2413fec2c5616e409d355a1497aab7bd&t=1715817113&sid=g44kE9

图片附件: New.gif (2014-7-9 09:47, 6.05 KB) / 下载次数 2287
http://www.inrm3d.cn/attachment.php?aid=21447&k=67ee3e86eebfe032fd8324093f195258&t=1715817113&sid=g44kE9

作者: 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 (2014-4-27 20:39, 17.3 KB) / 下载次数 2309
http://www.inrm3d.cn/attachment.php?aid=21448&k=bc92909d189711b66baa4396ecd15f0e&t=1715817113&sid=g44kE9

图片附件: New.gif (2014-7-9 09:48, 11.05 KB) / 下载次数 2329
http://www.inrm3d.cn/attachment.php?aid=21449&k=f94d1bea326464afd9d2afecf82465c7&t=1715817113&sid=g44kE9

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

115下载太慢了，急死人的！
2#的M集横向拉长了，是制图问题还是纵横像素比设置问题？

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

与程序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) / 下载次数 1841
http://www.inrm3d.cn/attachment.php?aid=21450&k=ae260a50ebc2a559485a436863828d78&t=1715817113&sid=g44kE9

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

与程序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) / 下载次数 1878
http://www.inrm3d.cn/attachment.php?aid=21451&k=907ea8614d7317cab5c270e57c0ded6b&t=1715817113&sid=g44kE9

作者: 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 (2014-4-27 21:01, 18.64 KB) / 下载次数 1742
http://www.inrm3d.cn/attachment.php?aid=21452&k=39b0268b16805a522cf84bd4143a9d53&t=1715817113&sid=g44kE9

图片附件: 捕获6.jpg (2014-4-28 18:35, 26.42 KB) / 下载次数 1739
http://www.inrm3d.cn/attachment.php?aid=21453&k=4b7045c7ca36674d7015f60032383d4b&t=1715817113&sid=g44kE9

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

伪3D效果

图片附件: 捕获7.jpg (2014-4-27 21:11, 26.31 KB) / 下载次数 1773
http://www.inrm3d.cn/attachment.php?aid=21454&k=621ab6c9c5ceeab6090e2acdc6cf05a0&t=1715817113&sid=g44kE9

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

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

图片附件: 未命名.JPG (2014-4-27 23:28, 31.23 KB) / 下载次数 1745
http://www.inrm3d.cn/attachment.php?aid=21459&k=4b80725d311bd3473654527d7e7f50eb&t=1715817113&sid=g44kE9

作者: 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），试用起来也不错
将文件解压后，
这样

再这样

这是网址http://pan.baidu.com/s/1pJHKXdD

图片附件: New.gif (2014-7-9 09:50, 5.68 KB) / 下载次数 1521
http://www.inrm3d.cn/attachment.php?aid=21465&k=af6e1fb7a103f893913271c5ecff0d1f&t=1715817113&sid=g44kE9

图片附件: 捕获77.PNG (2014-4-29 08:55, 8.6 KB) / 下载次数 1490
http://www.inrm3d.cn/attachment.php?aid=21466&k=c0ecb51414c872604a43ccebe64c24b5&t=1715817113&sid=g44kE9

作者: 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 (2014-4-29 09:22, 21.83 KB) / 下载次数 1516
http://www.inrm3d.cn/attachment.php?aid=21468&k=2d3e68430bd52eb139a2fda39d68384f&t=1715817113&sid=g44kE9

图片附件: New.gif (2014-7-9 09:51, 15.86 KB) / 下载次数 1510
http://www.inrm3d.cn/attachment.php?aid=21469&k=a7c9fa4556c908a3e4b80a1929881933&t=1715817113&sid=g44kE9

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

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

图片附件: 捕获10.PNG (2014-4-29 10:09, 28.33 KB) / 下载次数 1532
http://www.inrm3d.cn/attachment.php?aid=21470&k=35698909e4da5656b2ac79bd954a9dde&t=1715817113&sid=g44kE9

图片附件: New.gif (2014-7-9 09:52, 18.29 KB) / 下载次数 1528
http://www.inrm3d.cn/attachment.php?aid=21471&k=58f8aa29469fffd050fefe8fc2c0c23e&t=1715817113&sid=g44kE9

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

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

图片附件: 捕获12.PNG (2014-4-29 10:21, 32.29 KB) / 下载次数 1487
http://www.inrm3d.cn/attachment.php?aid=21472&k=e1c710632f3615e1e0b740befe564133&t=1715817113&sid=g44kE9

图片附件: New.gif (2014-7-9 09:52, 16.25 KB) / 下载次数 1486
http://www.inrm3d.cn/attachment.php?aid=21473&k=570e17a5fff484cfc879476932ca96d5&t=1715817113&sid=g44kE9

作者: 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 (2014-7-9 09:53, 8.89 KB) / 下载次数 1441
http://www.inrm3d.cn/attachment.php?aid=21474&k=16dc093ce5fac168fbdd1ea60fe78e0f&t=1715817113&sid=g44kE9

图片附件: New.jpg (2014-7-9 09:53, 14.41 KB) / 下载次数 1425
http://www.inrm3d.cn/attachment.php?aid=21475&k=005287bbef1f30df0a78786724ada052&t=1715817113&sid=g44kE9

作者: 柳烟 时间: 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) / 下载次数 1527
http://www.inrm3d.cn/attachment.php?aid=21477&k=ee9364d39d66c98f064e71ef94abc141&t=1715817113&sid=g44kE9

作者: 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 (2014-7-9 09:54, 31.35 KB) / 下载次数 1633
http://www.inrm3d.cn/attachment.php?aid=21492&k=6a04d92760898463671893c86e8e8ad6&t=1715817113&sid=g44kE9

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

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

图片附件: New.jpg (2014-7-9 09:55, 10.17 KB) / 下载次数 1653
http://www.inrm3d.cn/attachment.php?aid=21493&k=9849c6b1df367ec13a6b97f9ba087838&t=1715817113&sid=g44kE9

作者: 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 (2014-7-9 09:55, 24.51 KB) / 下载次数 1674
http://www.inrm3d.cn/attachment.php?aid=21496&k=ea95b5ed5d5210df2a6f7e7eac61e6f5&t=1715817113&sid=g44kE9

图片附件: New.jpg (2014-7-9 09:56, 20.86 KB) / 下载次数 1654
http://www.inrm3d.cn/attachment.php?aid=21497&k=5ec8da9d6836c65ce466be8f0e05d7aa&t=1715817113&sid=g44kE9

图片附件: New.jpg (2014-7-9 09:56, 18.42 KB) / 下载次数 1702
http://www.inrm3d.cn/attachment.php?aid=21498&k=03049d213a568245e95684f782085bb3&t=1715817113&sid=g44kE9

图片附件: New.jpg (2014-7-9 09:57, 22.11 KB) / 下载次数 1677
http://www.inrm3d.cn/attachment.php?aid=21499&k=4e7e8cb434bf6dee3b9b21173c93e5d7&t=1715817113&sid=g44kE9

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

35# lnszdzg
这个分形有特色。

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

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

图片附件: New.jpg (2014-7-9 09:57, 32.54 KB) / 下载次数 1666
http://www.inrm3d.cn/attachment.php?aid=21504&k=c7ebcc5e81957079bb0c73f146edb878&t=1715817113&sid=g44kE9

图片附件: New.jpg (2014-7-9 09:58, 34.63 KB) / 下载次数 1636
http://www.inrm3d.cn/attachment.php?aid=21505&k=bd2ef777eb2f59e57f9a3898dcca4dbf&t=1715817113&sid=g44kE9

作者: 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 (2014-7-9 09:59, 41.5 KB) / 下载次数 1329
http://www.inrm3d.cn/attachment.php?aid=21506&k=dc7f965715981f3075ed913f6d68c993&t=1715817113&sid=g44kE9

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

图片附件: New.jpg (2014-7-9 09:59, 34.91 KB) / 下载次数 1346
http://www.inrm3d.cn/attachment.php?aid=21507&k=03981b8e4dc88f45988bc18c23abaf01&t=1715817113&sid=g44kE9

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

还是蓝白调好看

图片附件: New.jpg (2014-7-9 10:00, 35.08 KB) / 下载次数 1400
http://www.inrm3d.cn/attachment.php?aid=21509&k=a43d9c8ca159192e8ce9530b11038a2c&t=1715817113&sid=g44kE9

作者: 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 (2014-7-9 10:00, 18.52 KB) / 下载次数 1405
http://www.inrm3d.cn/attachment.php?aid=21517&k=b75a8f62b20c5bbb6f9230a8bdd8da6f&t=1715817113&sid=g44kE9

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

图片附件: New.jpg (2014-7-9 10:01, 9.59 KB) / 下载次数 1393
http://www.inrm3d.cn/attachment.php?aid=21518&k=06e2145af4543dcc1c0faa9b258a5ef9&t=1715817113&sid=g44kE9

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

图片附件: New.jpg (2014-7-9 10:01, 36.13 KB) / 下载次数 1395
http://www.inrm3d.cn/attachment.php?aid=21519&k=ac7344a9355090a770713f1a460ef985&t=1715817113&sid=g44kE9

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

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) / 下载次数 1155
http://www.inrm3d.cn/attachment.php?aid=21553&k=63fb4d3c920a0cf92e3aecbb262ad730&t=1715817113&sid=g44kE9

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

图片附件: New.jpg (2014-7-9 10:03, 18.29 KB) / 下载次数 1167
http://www.inrm3d.cn/attachment.php?aid=21554&k=239488f59b592099c65454d083d61704&t=1715817113&sid=g44kE9

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

52# lnszdzg

点陷阱，用了HLS颜色

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

图片附件: New.jpg (2014-7-9 10:03, 25.35 KB) / 下载次数 1149
http://www.inrm3d.cn/attachment.php?aid=21555&k=d5f42c2f7e83413917c62ae28e6ae61d&t=1715817113&sid=g44kE9

图片附件: New.jpg (2014-7-9 10:04, 25.77 KB) / 下载次数 1161
http://www.inrm3d.cn/attachment.php?aid=21556&k=bc3fadd856cf49d619dd956027fac817&t=1715817113&sid=g44kE9

作者: 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 (2014-7-9 10:04, 31 KB) / 下载次数 1078
http://www.inrm3d.cn/attachment.php?aid=21563&k=c849c2ec9340256be00bbaeaf1a78607&t=1715817113&sid=g44kE9

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

同臂同色

图片附件: New.jpg (2014-7-9 10:05, 24.38 KB) / 下载次数 1089
http://www.inrm3d.cn/attachment.php?aid=21564&k=ecf2a6f08dc4be4b5f57610c8d99a290&t=1715817113&sid=g44kE9

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

图片附件: New.jpg (2014-7-9 10:05, 22.95 KB) / 下载次数 1066
http://www.inrm3d.cn/attachment.php?aid=21565&k=784c14106cc6f26e0432fbc5d0b648ce&t=1715817113&sid=g44kE9

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

图片附件: New.jpg (2014-7-9 10:05, 30.27 KB) / 下载次数 1095
http://www.inrm3d.cn/attachment.php?aid=21567&k=3c48a6be41b5b1dda900a75e3219f7f5&t=1715817113&sid=g44kE9

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

图片附件: New.jpg (2014-7-9 10:07, 28.82 KB) / 下载次数 1064
http://www.inrm3d.cn/attachment.php?aid=21585&k=559d6e4f38939728c96f15b19ad2c2e1&t=1715817113&sid=g44kE9

图片附件: New.jpg (2014-7-9 10:07, 31.51 KB) / 下载次数 1097
http://www.inrm3d.cn/attachment.php?aid=21586&k=7a5b002e6c3a53ab4ad1d14fb4d91f14&t=1715817113&sid=g44kE9

图片附件: New.jpg (2014-7-9 10:08, 44.14 KB) / 下载次数 1157
http://www.inrm3d.cn/attachment.php?aid=21587&k=7870d844bd5f63392da6f7cdf8ceffad&t=1715817113&sid=g44kE9

图片附件: New.jpg (2014-7-9 10:08, 28.16 KB) / 下载次数 1103
http://www.inrm3d.cn/attachment.php?aid=21588&k=db5eb66a66b07db99b93cce4970f5547&t=1715817113&sid=g44kE9

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

图片附件: New.jpg (2014-7-9 10:09, 43.13 KB) / 下载次数 1075
http://www.inrm3d.cn/attachment.php?aid=21598&k=4c4152d03a71fb571cfbfbca1e34d76b&t=1715817113&sid=g44kE9

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

图片附件: New.jpg (2014-7-9 10:09, 6.59 KB) / 下载次数 1043
http://www.inrm3d.cn/attachment.php?aid=21647&k=58887a6843578868dcb1d3ea6207defb&t=1715817113&sid=g44kE9

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

图片附件: New.jpg (2014-7-9 10:10, 23.68 KB) / 下载次数 1064
http://www.inrm3d.cn/attachment.php?aid=21648&k=e70d9ddddc4e7144a0882efc4b583af9&t=1715817113&sid=g44kE9

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

图片附件: New.jpg (2014-7-9 10:10, 23.88 KB) / 下载次数 1023
http://www.inrm3d.cn/attachment.php?aid=21649&k=c71d00ad14878d0e8ea652984306e98e&t=1715817113&sid=g44kE9

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

图片附件: New.jpg (2014-7-9 10:10, 36.49 KB) / 下载次数 1107
http://www.inrm3d.cn/attachment.php?aid=21650&k=040fde83f642e26f8f16ef1cdd820bc5&t=1715817113&sid=g44kE9

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

图片附件: New.jpg (2014-7-9 10:11, 25.83 KB) / 下载次数 1067
http://www.inrm3d.cn/attachment.php?aid=21657&k=15583ef4b4491169679b702a4f01e2cb&t=1715817113&sid=g44kE9

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

蝴蝶结陷阱

图片附件: New.jpg (2014-7-9 10:11, 30.96 KB) / 下载次数 1071
http://www.inrm3d.cn/attachment.php?aid=21658&k=b516b368ae12aa9878fd364121a63a70&t=1715817113&sid=g44kE9

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

67# lnszdzg
用什么曲线作陷阱?

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

图片附件: New.jpg (2014-7-9 10:12, 57.46 KB) / 下载次数 1246
http://www.inrm3d.cn/attachment.php?aid=21659&k=7c54bf17c6b638681617cc36b25edaaf&t=1715817113&sid=g44kE9

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

68# guxiaoping

像蝴蝶结哈，方程式||z|-|sin(Re(z))||

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

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

图片附件: 捕获81.jpg (2014-5-22 13:36, 44.32 KB) / 下载次数 1066
http://www.inrm3d.cn/attachment.php?aid=21660&k=d725dec1e28cb0acc1d0a2a43c34ed54&t=1715817113&sid=g44kE9

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

曲线：

图片附件: 捕获83.jpg (2014-5-22 13:39, 42.64 KB) / 下载次数 1024
http://www.inrm3d.cn/attachment.php?aid=21661&k=37565c705fb7f13ab7d243f4ab9dc20a&t=1715817113&sid=g44kE9

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

图片附件: 捕获84.jpg (2014-5-22 13:40, 52.95 KB) / 下载次数 1153
http://www.inrm3d.cn/attachment.php?aid=21663&k=2909ec5e052fbd6491896b32279e90b8&t=1715817113&sid=g44kE9

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

图片附件: 捕获85.jpg (2014-5-22 13:41, 42.67 KB) / 下载次数 1046
http://www.inrm3d.cn/attachment.php?aid=21664&k=a378f09653499a12577961d4de46a7c5&t=1715817113&sid=g44kE9

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

图片附件: 捕获86.jpg (2014-5-22 13:42, 42.65 KB) / 下载次数 1047
http://www.inrm3d.cn/attachment.php?aid=21665&k=c2c5a04db414137571b7de687deec448&t=1715817113&sid=g44kE9

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

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

图片附件: 捕获91.jpg (2014-5-22 13:46, 27.64 KB) / 下载次数 890
http://www.inrm3d.cn/attachment.php?aid=21666&k=2b5892c5d70960a1ba842cfba9d634bb&t=1715817113&sid=g44kE9

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

图片附件: New.jpg (2014-7-9 10:12, 24.81 KB) / 下载次数 907
http://www.inrm3d.cn/attachment.php?aid=21667&k=c108718bfa6d854da6f7d8b1c10b7a9e&t=1715817113&sid=g44kE9

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

DEM算法下的陷阱

图片附件: New.jpg (2014-7-9 10:13, 33.08 KB) / 下载次数 907
http://www.inrm3d.cn/attachment.php?aid=21668&k=d1439cf22e66049f6079813b395721f9&t=1715817113&sid=g44kE9

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

图片附件: New.jpg (2014-7-9 10:13, 49.82 KB) / 下载次数 923
http://www.inrm3d.cn/attachment.php?aid=21671&k=09471ed54831436fb352b66640cbad88&t=1715817113&sid=g44kE9

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

图片附件: New.jpg (2014-7-9 10:13, 48.37 KB) / 下载次数 907
http://www.inrm3d.cn/attachment.php?aid=21672&k=bd14e36d68c864c5b01285fc137c6fed&t=1715817113&sid=g44kE9

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

图片附件: New.jpg (2014-7-9 10:14, 29.07 KB) / 下载次数 1047
http://www.inrm3d.cn/attachment.php?aid=21673&k=332d2dcf83c58ce271bf576139dbe24c&t=1715817113&sid=g44kE9

作者: 柳烟 时间: 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 (2014-5-24 11:30, 31.19 KB) / 下载次数 1035
http://www.inrm3d.cn/attachment.php?aid=21679&k=9d5f6a83905e3e1a7a89780a3e079cc2&t=1715817113&sid=g44kE9

作者: 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 (2014-7-9 10:14, 25.93 KB) / 下载次数 1039
http://www.inrm3d.cn/attachment.php?aid=21684&k=b7d72335b6ea849107c876b44ee409bc&t=1715817113&sid=g44kE9

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

图片附件: New.jpg (2014-7-9 10:15, 27.63 KB) / 下载次数 1059
http://www.inrm3d.cn/attachment.php?aid=21685&k=c2284957fc5d0de2037cb0af8b830574&t=1715817113&sid=g44kE9

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

喇叭花与桃花

图片附件: 捕获119.PNG (2014-5-25 21:22, 25.08 KB) / 下载次数 1079
http://www.inrm3d.cn/attachment.php?aid=21697&k=55b74950e479050569e9c49f66367020&t=1715817113&sid=g44kE9

图片附件: New.jpg (2014-7-9 10:15, 21.16 KB) / 下载次数 1046
http://www.inrm3d.cn/attachment.php?aid=21698&k=d41d9665edb61f93e6d4781f14acfabd&t=1715817113&sid=g44kE9

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

这些真花漂亮！

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

五星

图片附件: New.jpg (2014-7-9 10:16, 16.76 KB) / 下载次数 572
http://www.inrm3d.cn/attachment.php?aid=21715&k=1b92259abc69574dd358f62ee7bad8a8&t=1715817113&sid=g44kE9

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

图片附件: New.jpg (2014-7-9 10:16, 19.17 KB) / 下载次数 569
http://www.inrm3d.cn/attachment.php?aid=21716&k=b13f21ed56203ceebd16e002a091ceaf&t=1715817113&sid=g44kE9

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

图片附件: New.jpg (2014-7-9 10:16, 19.93 KB) / 下载次数 569
http://www.inrm3d.cn/attachment.php?aid=21717&k=83003e47d77414559c52c7c1437c58d6&t=1715817113&sid=g44kE9

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

图片附件: New.jpg (2014-7-9 10:17, 23.31 KB) / 下载次数 715
http://www.inrm3d.cn/attachment.php?aid=21721&k=22c9eee272bd2ceb70804fc78c4e2759&t=1715817113&sid=g44kE9

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

94# 柳烟

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

图片附件: New.jpg (2014-7-9 10:17, 18.08 KB) / 下载次数 583
http://www.inrm3d.cn/attachment.php?aid=21722&k=f7671264f1b361bc0b23de7ba327c790&t=1715817113&sid=g44kE9

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

图片附件: New.jpg (2014-7-9 10:18, 31.98 KB) / 下载次数 586
http://www.inrm3d.cn/attachment.php?aid=21724&k=240bdaf79bbcdc4b797b1b9c654c7ee5&t=1715817113&sid=g44kE9

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

心形线下的M集

掏空肚子

图片附件: New.jpg (2014-7-9 10:18, 21.3 KB) / 下载次数 555
http://www.inrm3d.cn/attachment.php?aid=21725&k=6c8e6b751203ebd2b493fa03d8604cca&t=1715817113&sid=g44kE9

图片附件: New.jpg (2014-7-9 10:19, 23.62 KB) / 下载次数 569
http://www.inrm3d.cn/attachment.php?aid=21726&k=f05b18744bb2e7ffbc5418272edeb6cd&t=1715817113&sid=g44kE9

图片附件: New.jpg (2014-7-9 10:19, 26.54 KB) / 下载次数 570
http://www.inrm3d.cn/attachment.php?aid=21727&k=ef1a60c05e6b99bd46691a5d2dafdaf6&t=1715817113&sid=g44kE9

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

双纽线

三叶玫瑰

图片附件: New.jpg (2014-7-9 10:20, 27.18 KB) / 下载次数 588
http://www.inrm3d.cn/attachment.php?aid=21728&k=cc2e0da82d7a735f351cbc66f5a61cd2&t=1715817113&sid=g44kE9

图片附件: New.jpg (2014-7-9 10:20, 28.88 KB) / 下载次数 579
http://www.inrm3d.cn/attachment.php?aid=21729&k=8961c50c7590113a54c2637b88b30441&t=1715817113&sid=g44kE9

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

棋盘背景

图片附件: 捕获138.PNG (2014-5-27 21:29, 49.69 KB) / 下载次数 577
http://www.inrm3d.cn/attachment.php?aid=21739&k=9e69b54340ac4f37eada261a54a00916&t=1715817113&sid=g44kE9

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

DEM算法下的J集，对像素的要求太高了，速度很慢啊。

图片附件: 捕获1392.jpg (2014-5-28 16:18, 27.78 KB) / 下载次数 711
http://www.inrm3d.cn/attachment.php?aid=21742&k=3206a2c44a422ba00ed03e004474a6de&t=1715817113&sid=g44kE9

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