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) / 下载次数 2849
http://forums.netpad.net.cn/attachment.php?aid=14872&k=5fa53be0591bf04b989ef219c4ffa021&t=1792500730&sid=A0m2r2

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

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

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

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

作者: 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) / 下载次数 2821
http://forums.netpad.net.cn/attachment.php?aid=14874&k=a2cb0424ac04ae08ab6e05e82c06ac62&t=1792500730&sid=A0m2r2

作者: 柳烟 时间: 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) / 下载次数 2803
http://forums.netpad.net.cn/attachment.php?aid=14875&k=cfd62abb32bd298f6e0f11b20f26a56b&t=1792500730&sid=A0m2r2

作者: 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) / 下载次数 1812
http://forums.netpad.net.cn/attachment.php?aid=14878&k=37df7a8df3cc17c9d96841a0c268b372&t=1792500730&sid=A0m2r2

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

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

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

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

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

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

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

作者: 榕坚 时间: 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) / 下载次数 1813
http://forums.netpad.net.cn/attachment.php?aid=14885&k=e02d083f2c19159d74c4ee074e03c412&t=1792500730&sid=A0m2r2

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

15# dyk

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

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

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

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

作者: 柳烟 时间: 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) / 下载次数 2234
http://forums.netpad.net.cn/attachment.php?aid=14892&k=ae3415a47acad8841c00f2a5e863a098&t=1792500730&sid=A0m2r2

作者: 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) / 下载次数 1969
http://forums.netpad.net.cn/attachment.php?aid=14895&k=f8e450edb74bd14ea840d59f061a1591&t=1792500730&sid=A0m2r2

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

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

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

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

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

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

29# 榕坚

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

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

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

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

作者: 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) / 下载次数 2402
http://forums.netpad.net.cn/attachment.php?aid=14904&k=1e2fe7931850b5ec1008462fbc88bbe6&t=1792500730&sid=A0m2r2

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

上面分形的两个变形

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

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

作者: 柳烟 时间: 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) / 下载次数 2527
http://forums.netpad.net.cn/attachment.php?aid=14912&k=405a376e21985a32bbe97cd7bebe968f&t=1792500730&sid=A0m2r2

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

作者: 榕坚 时间: 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/)