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

标题: DEM法作J集大练兵 [打印本页]

作者: 柳烟 时间: 2011-10-28 00:11 标题: DEM法作J集大练兵

http://www.ijon.de/mathe/julia/some_julia_sets_3.html
f(z)=1/(z3+dz+c) mit c=0 und d=-3(1+i), dargestellt auf [-3;3]x[-3;3].

图片附件: b005.png (2011-10-28 00:11, 9.87 KB) / 下载次数 2197
http://forums.netpad.net.cn/attachment.php?aid=14872&k=07dd844523b5c52ee631ae53d44c5045&t=1759985112&sid=j0Y4jx

作者: 柳烟 时间: 2011-10-28 00:25

我试着弄，结果扫不出图来。我的问题来了：
我按xyj200909 板友的方法作 J集，D0我认为仍是1.导数求出应是：

帮我看看这导数是否有问题？

图片附件: 未命名.JPG (2011-10-28 00:25, 2.18 KB) / 下载次数 2166
http://forums.netpad.net.cn/attachment.php?aid=14873&k=92bd82dd09606fd88fd1a5f8f02990eb&t=1759985112&sid=j0Y4jx

作者: dyk 时间: 2011-10-28 10:04

我认为Z'=1

作者: dyk 时间: 2011-10-28 10:08

DEM法作黑白分形速度快一点,但难调整,用et或模逼近容易些,就是慢点.下面是用模逼近扫出的,用时2个半.

图片附件: J58.GIF (2011-10-29 11:14, 24.33 KB) / 下载次数 2165
http://forums.netpad.net.cn/attachment.php?aid=14874&k=46ed52efe2605d14a305e678f299f4f8&t=1759985112&sid=j0Y4jx

作者: 柳烟 时间: 2011-10-28 13:38

4# dyk
请教一下模逼近法大致思路，谢谢。

作者: 柳烟 时间: 2011-10-28 13:41

2楼的导数没有问题，我昨晚不知如何搞的，也是按这导数做，结果扫不出图，现在重头到尾再仔细做一遍，居然成了，昨晚问题出在何处，也懒得去检查了，满屏算式，繁。按xyj200909的作色法，扫一金色J集。

图片附件: (z^3+dz)的倒数的J集.jpg (2011-10-29 11:14, 32.05 KB) / 下载次数 2154
http://forums.netpad.net.cn/attachment.php?aid=14875&k=a42d8b34b3ca902c440f16b5aee9f00d&t=1759985112&sid=j0Y4jx

作者: dyk 时间: 2011-10-28 14:35

5# 柳烟

按照常规作法得到终点Z，定义由点z到Z的自定义变换,再构造一个点z'，使其与点z非常接近，对z'实施自定义的变换后得到Z',计算这两终点的模的差，其差的绝对值在某个范围内的点着白色（或黑色），其它的点着黑色（或白色）即可。思路就这么简单。把源文件附上，一看便知。

作者: dyk 时间: 2011-10-28 14:39

无奈，才250k,怎么超出限制，源文件传不上去，压缩一下试一试，还是不行

作者: 柳烟 时间: 2011-10-28 15:18

8# dyk
好象超过200K就不能存入论坛。用分卷压缩后，再存入论坛。

作者: 柳烟 时间: 2011-10-28 15:20

7# dyk
谢谢指点迷途，懂了。

作者: 柳烟 时间: 2011-10-28 17:35

f(z)=1/((0,15+0,15i)z5+z3+(-3+3i)z), dargestellt auf [-3;3]x[-3;3].

我整出的效果图是：

大家整整看，如何整出原图。

图片附件: b036.png (2011-10-28 17:35, 13.21 KB) / 下载次数 1423
http://forums.netpad.net.cn/attachment.php?aid=14878&k=c6aa76d3d6332d7646630d85b5a0f6e6&t=1759985112&sid=j0Y4jx

图片附件: 未命名.JPG (2011-10-28 17:35, 47.41 KB) / 下载次数 1436
http://forums.netpad.net.cn/attachment.php?aid=14879&k=ff3dede8fa265d45b449300ec1e40046&t=1759985112&sid=j0Y4jx

作者: 柳烟 时间: 2011-10-28 19:22

再调一个反色的，并附上文件。

(z^3+dz)的倒数的J集.gsp (21.52 KB)

图片附件: 未命名.JPG (2011-10-28 19:22, 30.06 KB) / 下载次数 1466
http://forums.netpad.net.cn/attachment.php?aid=14880&k=e422f0d28e02ee8278ad44bdcb1bb42f&t=1759985112&sid=j0Y4jx

图片附件: 未命名.JPG (2011-10-28 21:46, 31.77 KB) / 下载次数 1447
http://forums.netpad.net.cn/attachment.php?aid=14882&k=24f3e1bf8631723135ccc4024be1ae2d&t=1759985112&sid=j0Y4jx

附件: (z^3+dz)的倒数的J集.gsp (2011-10-28 21:46, 21.52 KB) / 下载次数 2824
http://forums.netpad.net.cn/attachment.php?aid=14883&k=52f5d1eea80c5263460fad580a1c16b2&t=1759985112&sid=j0Y4jx

作者: 榕坚 时间: 2011-10-28 20:54

11# 柳烟

逃逸值太大了，调小一些可能就会把逃逸部分去掉。

作者: 柳烟 时间: 2011-10-28 21:47

13# 榕坚
我试了，不成。

作者: dyk 时间: 2011-10-29 08:36

用et逼近，很好调整。

图片附件: J60.GIF (2011-10-29 11:16, 26.82 KB) / 下载次数 1434
http://forums.netpad.net.cn/attachment.php?aid=14885&k=734c86b4997614c4bb022c0876b60c36&t=1759985112&sid=j0Y4jx

作者: xuefeiyang 时间: 2011-10-29 09:44

15# dyk

请问迭代了多少次？绘制这样一幅图需要多长时间？

作者: 榕坚 时间: 2011-10-29 09:47

11# 柳烟
你再检查一下算式，这个图形的做法没有什么特殊性：

图片附件: 捕获.JPG (2011-10-29 09:47, 38.77 KB) / 下载次数 1300
http://forums.netpad.net.cn/attachment.php?aid=14887&k=9206ff0a861cc9d2f8b90323fe794874&t=1759985112&sid=j0Y4jx

作者: 柳烟 时间: 2011-10-29 10:23

15# dyk
作得太好了，附上源文件吧，让我学学这et逼近法。谢谢。

作者: dyk 时间: 2011-10-29 10:24

16# xuefeiyang

n=500,采样数量1000,速度：0.15，用时5

作者: 柳烟 时间: 2011-10-29 10:25

17# 榕坚
看来我的计算或者其它出问题了，我再推倒重来一遍看看。

作者: 柳烟 时间: 2011-10-29 10:26

19# dyk
用时5个小时？我真服了你的耐心。

作者: dyk 时间: 2011-10-29 10:33

我的机器太旧运行太慢，晚上调好，睡觉。早上起来图完成了。其实n=300即可。

作者: 榕坚 时间: 2011-10-29 10:37

22# dyk

我扫的那幅n=100

作者: 柳烟 时间: 2011-10-29 13:46

我推倒做了二次，每次都不一样，怪怪的。我再干一次，若不成功，我将病态文件发上来，大家看看，我完全按xy老师的制作流程进行。

作者: dyk 时间: 2011-10-29 14:13

18# 柳烟

我是跟向老师学的，传源文件，供参考。

附件: 黑白分形J.gsp (2011-10-29 14:13, 184.08 KB) / 下载次数 1879
http://forums.netpad.net.cn/attachment.php?aid=14892&k=9ad57a9fe7fd3ccfbe35d9278e3c2c90&t=1759985112&sid=j0Y4jx

作者: dyk 时间: 2011-10-29 14:19

23# 榕坚

100时，不够清晰，我就调大迭代次数和采样数量睡觉去了。一般我机扫黑白图都在两个小时左右，没有耐心，让机器运转，人休息。

作者: xuefeiyang 时间: 2011-10-29 14:25

这个分形需要的最大迭代次数为374次左右。太大浪费时间，太小了，细节有些是显示不出来的。

作者: 柳烟 时间: 2011-10-29 14:59

推倒第三次，可能成功了，我正在扫图，老是把图扫不清晰。dxy老师张帖在此的那张，清晰度太好了。
扫好了，费时2个钟头左右，按dxy老师给出的参数进行，然扫出的图与他发表的图一比，差了一个挡次了。

((0,15+0,15i)z5+z3+(-3+3i)z)的倒数的dem.gsp (25.45 KB)
再扫一个女的：

图片附件: 未命名.JPG (2011-10-29 16:44, 37.17 KB) / 下载次数 1560
http://forums.netpad.net.cn/attachment.php?aid=14895&k=aacaa81e3cb5663a413845a6530a9fc3&t=1759985112&sid=j0Y4jx

附件: ((0,15+0,15i)z5+z3+(-3+3i)z)的倒数的dem.gsp (2011-10-29 16:51, 25.45 KB) / 下载次数 2749
http://forums.netpad.net.cn/attachment.php?aid=14896&k=cc9e429ee8d971fdde093fc17a1e281b&t=1759985112&sid=j0Y4jx

图片附件: 未命名.JPG (2011-10-29 18:08, 78.83 KB) / 下载次数 1731
http://forums.netpad.net.cn/attachment.php?aid=14897&k=3643c78d2364bdc5d48afe4bc12fca62&t=1759985112&sid=j0Y4jx

作者: 榕坚 时间: 2011-10-29 15:54

19# dyk
迭代500次，采样数1000.那应该还是不错的电脑了。这是迭代100次，采样数500的彩色图：

图片附件: 捕获.JPG (2011-10-29 15:57, 49.55 KB) / 下载次数 1311
http://forums.netpad.net.cn/attachment.php?aid=14894&k=aa1e2dacdfdde5bfd19b4d490019b00b&t=1759985112&sid=j0Y4jx

作者: dyk 时间: 2011-10-29 16:36

29# 榕坚

漂亮，做彩图比做黑白的要快一些。估计这幅用不了半小时。

作者: 柳烟 时间: 2011-10-29 19:32

将就dyk老师的源文件，扫一个女的：

图片附件: 未命名.JPG (2011-10-29 19:32, 83.31 KB) / 下载次数 2065
http://forums.netpad.net.cn/attachment.php?aid=14900&k=96111ca1df76a08e7c18e3f92811d0e9&t=1759985112&sid=j0Y4jx

作者: dyk 时间: 2011-10-29 22:06

柳老师用的慢速吧，黑色背景有条纹。

作者: 柳烟 时间: 2011-10-29 22:36

32# dyk
够慢的了，几乎看不到扫描线动，扫出的，好象在原图中看不到条纹。

作者: 柳烟 时间: 2011-10-29 23:20

问题：如何将dem法用于N集制造？能否说说大致的思维框架？我茫无目的的尝试后，弄出来空白一片，有时扫出来几个点。

作者: dyk 时间: 2011-10-29 23:49

33# 柳烟

那就是存成JPG时质量有损失。

作者: dyk 时间: 2011-10-30 11:48

来个简单的f(z)=z^3+z+.05 n=80

图片附件: J61.jpg (2011-10-30 11:48, 57.47 KB) / 下载次数 2018
http://forums.netpad.net.cn/attachment.php?aid=14904&k=b9b460e6d5c5b1e9990a6cb37b03aa7e&t=1759985112&sid=j0Y4jx

作者: dyk 时间: 2011-10-30 13:51

上面分形的两个变形

图片附件: J61_1.jpg (2011-10-30 13:51, 48.91 KB) / 下载次数 1836
http://forums.netpad.net.cn/attachment.php?aid=14908&k=065245047aeed0e6656f6880036e066d&t=1759985112&sid=j0Y4jx

图片附件: J61_2.jpg (2011-10-30 13:51, 48.84 KB) / 下载次数 1791
http://forums.netpad.net.cn/attachment.php?aid=14909&k=a5e8b8eaa392e5a433d685121386b535&t=1759985112&sid=j0Y4jx

作者: 柳烟 时间: 2011-10-30 16:02

未命名.jpg

素材来自一楼所给链接，原分形公式是：f(z)=z6+1,01*z, dargestellt auf [-1,2;1,2]x[-1,2;1,2].我按dem法造出后，发现原网站所给的图片各乳沟很深，而我的图片各乳沟很浅，不美，所以将复函数公式中的第二项1,01*z,改为1,1z,扫的图片。

f(z)=z6+1,1z的J集DEM.gsp (20 KB)

图片附件: 未命名.jpg (2011-10-30 16:02, 50.58 KB) / 下载次数 2090
http://forums.netpad.net.cn/attachment.php?aid=14912&k=613971141252f7eb640672a7a9f127c7&t=1759985112&sid=j0Y4jx

附件: f(z)=z6+1,1z的J集DEM.gsp (2011-10-30 16:05, 20 KB) / 下载次数 2943
http://forums.netpad.net.cn/attachment.php?aid=14913&k=fa14fe3c30051b1672b28a7bbd4c25b8&t=1759985112&sid=j0Y4jx

作者: 榕坚 时间: 2011-10-30 16:34

38# 柳烟

这应该是迭代次数没办法上去的原因。

作者: 柳烟 时间: 2011-10-30 19:45

一般彩色分形，迭代1000次，用这dem法，或微差法，机子的迭代次数一下就迭下来了，大约缩一半。扫一张好图，费时费力，老外的这软件太秀气了。

作者: mjj_ljh 时间: 2011-10-30 20:03

柳烟老师的用词还真让人忍俊不禁

。

作者: 柳烟 时间: 2011-10-30 22:15

41# mjj_ljh
呵呵，谢谢夸奖，胡乱说说，问好梅老师。

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