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

标题: (z^2-Q^2)(z^2-c^2)系列分形 [打印本页]

作者: lnszdzg 时间: 2015-1-28 20:49 标题: (z^2-Q^2)(z^2-c^2)系列分形

前几天，老师们完了(z-1)(z+1)(z-b+a)(z-b-a) 的牛迭M集征解(越玩越有趣)
今天，在Mathcad中也实现了这个，而且发现：这一系列的结构是相似的，可以统一在一个分形之中。

图片附件: 1.jpg (2015-1-28 21:49, 14.48 KB) / 下载次数 2390
http://develop.forums.netpad.net.cn/attachment.php?aid=23631&k=2ea208542a71543e5c764b1937a9a94c&t=1749930301&sid=bLX5X5

图片附件: 2.jpg (2015-1-28 21:49, 12.52 KB) / 下载次数 2405
http://develop.forums.netpad.net.cn/attachment.php?aid=23632&k=c841fa711570b13d939a70500753a3f1&t=1749930301&sid=bLX5X5

图片附件: 3.jpg (2015-1-28 21:49, 17.06 KB) / 下载次数 2426
http://develop.forums.netpad.net.cn/attachment.php?aid=23633&k=2fa2666b697a284a4a9e27a352dd94c6&t=1749930301&sid=bLX5X5

图片附件: 4.jpg (2015-1-28 21:49, 19.34 KB) / 下载次数 2396
http://develop.forums.netpad.net.cn/attachment.php?aid=23634&k=80df6de39844ebaef82bd6286f0f62d5&t=1749930301&sid=bLX5X5

图片附件: 5.jpg (2015-1-28 21:49, 18.75 KB) / 下载次数 2349
http://develop.forums.netpad.net.cn/attachment.php?aid=23635&k=2d9350600f8282d77a9f3078040a0a73&t=1749930301&sid=bLX5X5

图片附件: 6.jpg (2015-1-28 21:49, 17.87 KB) / 下载次数 2453
http://develop.forums.netpad.net.cn/attachment.php?aid=23636&k=3a4b9d126f5d472093a7eedd3db58475&t=1749930301&sid=bLX5X5

图片附件: 7.jpg (2015-1-28 21:49, 17.72 KB) / 下载次数 2353
http://develop.forums.netpad.net.cn/attachment.php?aid=23637&k=d95f502c47b8396adf31fb844fb8d842&t=1749930301&sid=bLX5X5

图片附件: a.jpg (2015-1-28 21:50, 17.75 KB) / 下载次数 2415
http://develop.forums.netpad.net.cn/attachment.php?aid=23638&k=ac970e9c6c89a80de42b93735eb15e13&t=1749930301&sid=bLX5X5

图片附件: 8.jpg (2015-1-28 21:49, 18.83 KB) / 下载次数 2400
http://develop.forums.netpad.net.cn/attachment.php?aid=23639&k=fb7141cc4314ff61090919ee3e10e197&t=1749930301&sid=bLX5X5

图片附件: 9.jpg (2015-1-28 21:49, 17.78 KB) / 下载次数 2420
http://develop.forums.netpad.net.cn/attachment.php?aid=23640&k=bbe4be74c7210f019ca55fee34586186&t=1749930301&sid=bLX5X5

作者: lnszdzg 时间: 2015-1-28 20:53

而且，还有一些很有意思的变种

应剑客之约，让我邀请大家到Mathcad吧做客。
希望几位老师能到Mathcad吧看看
http://tieba.baidu.com/f?kw=mathcad&fr=ala0

图片附件: b.jpg (2015-1-28 22:01, 15.04 KB) / 下载次数 2355
http://develop.forums.netpad.net.cn/attachment.php?aid=23641&k=0a0ba17ce9f5ffda11deff9140cd2486&t=1749930301&sid=bLX5X5

图片附件: c.jpg (2015-1-28 22:01, 18.69 KB) / 下载次数 2426
http://develop.forums.netpad.net.cn/attachment.php?aid=23642&k=ebe196d9c76d1176191e27d35edcc8ac&t=1749930301&sid=bLX5X5

图片附件: d.jpg (2015-1-28 22:01, 20.18 KB) / 下载次数 2408
http://develop.forums.netpad.net.cn/attachment.php?aid=23643&k=087d5bc0c4412a7361f29f14961922e0&t=1749930301&sid=bLX5X5

作者: lnszdzg 时间: 2015-1-28 21:20

UF变种

图片附件: e.jpg (2015-1-28 22:03, 14.9 KB) / 下载次数 2382
http://develop.forums.netpad.net.cn/attachment.php?aid=23644&k=40c388f47f3135eb5b9734a1d0f4b815&t=1749930301&sid=bLX5X5

作者: lnszdzg 时间: 2015-1-28 21:25

图片附件: f.jpg (2015-1-28 22:04, 18.12 KB) / 下载次数 2365
http://develop.forums.netpad.net.cn/attachment.php?aid=23645&k=3fcfdc889523b2adc82e06d253adca24&t=1749930301&sid=bLX5X5

图片附件: g.jpg (2015-1-28 22:04, 17.3 KB) / 下载次数 2423
http://develop.forums.netpad.net.cn/attachment.php?aid=23646&k=42602bfff1ee275e3b78ca76c33100ef&t=1749930301&sid=bLX5X5

图片附件: h.jpg (2015-1-28 22:04, 25.46 KB) / 下载次数 2361
http://develop.forums.netpad.net.cn/attachment.php?aid=23647&k=dd4a125812cd1bce6e8bb894081f7cd9&t=1749930301&sid=bLX5X5

作者: lnszdzg 时间: 2015-1-28 22:35

图片附件: 3.jpg (2015-1-28 22:44, 30.65 KB) / 下载次数 2418
http://develop.forums.netpad.net.cn/attachment.php?aid=23650&k=6267859362bfb5b7a643da21ce14d2ee&t=1749930301&sid=bLX5X5

作者: 柳烟 时间: 2015-1-28 22:45

杜老兄越玩越高超，快入化境了，祝贺。

作者: lnszdzg 时间: 2015-1-29 15:52

也是变种

图片附件: 捕获201501291.PNG (2015-1-29 15:52, 51.17 KB) / 下载次数 2393
http://develop.forums.netpad.net.cn/attachment.php?aid=23654&k=e4a359af56269d4ab9a85f6c924f97d8&t=1749930301&sid=bLX5X5

作者: lnszdzg 时间: 2015-1-29 15:54

UF代码，很简单，却很有效哈
Newton_Mandelbrot1 {
init:
  complex Q = @Q
  int p = @power
  c = #pixel
  z = sqrt((Q^2+c^p)/6 )
oop:
  err = (z^4-(c^p+Q^2)*z^2+c^p)/(4*z^3-2*(c^p+Q^2)*z)
  z = z - err
bailout:
  |err| > @bailout
default:
  title = "Newton_Mandelbrot1"
  maxiter = 2000
  periodicity = 0
  center = (0,0)
  magn = 1.0
param power
caption = "power"
default = 1
  endparam
  param Q
caption = "Q"
default = (1.0,0.0)
  endparam
param bailout
caption = "Bailout"
default = 0.000001
endparam
}

作者: lnszdzg 时间: 2015-1-29 15:57

“蝴蝶结”的个数有p确定，如果增加变换，得到的图形也很美。

作者: lnszdzg 时间: 2015-1-29 16:03

邀请大家到Mathcad吧做客。
http://tieba.baidu.com/f?kw=mathcad&fr=ala0

作者: lnszdzg 时间: 2015-1-29 23:19

最近画板论坛老是打不开,郁闷啊！

作者: lnszdzg 时间: 2015-2-1 21:53

不会平滑啊

图片附件: 3.jpg (2015-2-1 22:16, 11.37 KB) / 下载次数 2315
http://develop.forums.netpad.net.cn/attachment.php?aid=23685&k=ade3ce93073ba7a1629ad2477592c573&t=1749930301&sid=bLX5X5

作者: xiaongxp 时间: 2015-2-1 22:02

12# lnszdzg
好图！
杜老师，N集不平滑才美呢，其等et圈的密度反映了根的收敛速度。

作者: xklppp 时间: 2015-3-2 18:37

(z^m-q^a)(z^m-c^b)
z0=[m!m!(q^a+c^b)/(2m)!]^(1/m)

z^m[c^b(z^m-c^b)-q^a]+c^b
z0={m!m![q^a+c^(2b)]/(2m)!}^(1/m)

参数平面：c=d*c^f+e/c^g
动力平面：z=d*z^f+e/z^g

作者: xiaongxp 时间: 2015-3-9 23:39

邀请大家到Mathcad吧做客。
http://tieba.baidu.com/f?kw=mathcad&fr=ala0
lnszdzg 发表于 2015-1-29 16:03

MC吧的学习氛围令人羡慕，无私、热情、有深度、……

作者: lnszdzg100 时间: 2015-3-20 19:30

15# xiaongxp

我上不去画板论坛原来是被学校屏蔽了

作者: lnszdzg100 时间: 2015-3-21 20:27

图片附件: 捕获12101.jpg (2015-3-21 20:27, 27.81 KB) / 下载次数 2058
http://develop.forums.netpad.net.cn/attachment.php?aid=23990&k=c12ffabd913a0cea2601c29112d50920&t=1749930301&sid=bLX5X5

作者: zhongba 时间: 2015-3-21 20:43

17# lnszdzg100
赞！
这个几何画板办不到啰。

作者: 榕坚 时间: 2015-3-22 09:58

18# zhongba

可以尝试：

图片附件: 捕获.JPG (2015-3-22 09:58, 40.82 KB) / 下载次数 2174
http://develop.forums.netpad.net.cn/attachment.php?aid=23992&k=e3e1fd2774a9f2ff5ac6404685c72f0a&t=1749930301&sid=bLX5X5

作者: lnszdzg100 时间: 2015-3-22 20:21

19# 榕坚

这个不错，很好啊！

作者: 榕坚 时间: 2015-3-22 21:41

头像的大小可以随着迭代过程变化了，但都是正放的，怎么让它随着旋转呢？

作者: xiaongxp 时间: 2015-3-22 22:50

21# 榕坚
用三点粘贴图片，可保图片跟随旋转。
圆形图片太难找了

图片的逆迭代J集.gsp (84.64 KB)

图片附件: 图片的逆迭代J集.jpg (2015-3-22 23:00, 36.6 KB) / 下载次数 2590
http://develop.forums.netpad.net.cn/attachment.php?aid=23993&k=f35db1f3a34ec8396c1dd897f0419cfb&t=1749930301&sid=bLX5X5

附件: 图片的逆迭代J集.gsp (2015-3-22 23:00, 84.64 KB) / 下载次数 4157
http://develop.forums.netpad.net.cn/attachment.php?aid=23994&k=ddf7b0a8cda2ca097caca0d39ba86a99&t=1749930301&sid=bLX5X5

作者: 榕坚 时间: 2015-3-24 16:24

可以转了：

图片附件: 2.JPG (2015-3-24 16:24, 40.12 KB) / 下载次数 2398
http://develop.forums.netpad.net.cn/attachment.php?aid=24001&k=cae58892e5835d0be9561ec0c46811e1&t=1749930301&sid=bLX5X5

作者: 榕坚 时间: 2015-3-24 20:28

有点像了：

图片附件: 1.jpg (2015-3-24 20:40, 27.05 KB) / 下载次数 2026
http://develop.forums.netpad.net.cn/attachment.php?aid=24002&k=6a43204d1ee85601132beb6996c01da3&t=1749930301&sid=bLX5X5

作者: lnszdzg100 时间: 2015-3-27 10:32

22# xiaongxp

图片是经过加工才呈圆形的。

作者: lnszdzg100 时间: 2015-3-27 10:34

24# 榕坚

榕老师厉害啊！

作者: lnsz_dzg 时间: 2021-5-20 16:41

已经忘完了

作者: lnsz_dzg 时间: 2021-5-20 16:43

Circlesbelt {
init:
  int iter = 0
  float num = @cnum
  float Xabs = 0
  float Yabs = 0
  float Dsgd0 = 0

  float Circle = 0
  float ZtoPsqd = 0
  float Rc = @r;/#magn

  float Rm = @Rr
  float RcSqd = Rc
  bool Trapped = false
loop:
if  (@tran==0 &&!Trapped) || @tran==1
  iter = iter + 1
  if iter > @skip
Xabs = real(#z)
Yabs = imag(#z)
Dsgd0 = abs(sqrt((Xabs+0.0)^2 + (Yabs -0)^2) - Rm)
if Dsgd0 < RcSqd
   Trapped = true
   ZtoPsqd= Dsgd0
      Circle = iter
      float argz = atan2(#z)
         if argz<0
            argz = argz + 2*#pi
         endif
         argz = argz/(2*#pi)
endif
  endif
endif
final:
  if !trapped
#solid = true
  else
Circle = ((Circle -1 + @off) % num)/num
float Ratio = sqrt(ZtoPsqd/Rcsqd)
; float ColorIndex = 110 * Ratio + Circle * 140
;#index = (ColorIndex + 1) % 256 /256
#color = hsl(argz*6.0 ,Ratio*0.5 +0.5,(0.90-0.6*Ratio)*0.9)
  endif
default:
  title = "1 belt"
param tran
caption="type Trap"
default=1
enum="first" "last"
endparam
param cnum
caption = "color num"
default = 6.0
  endparam

  param r
caption = "Circle radius"
default = 0.3
  endparam
  param Rr
caption = "Out radius"
default = 1.5
  endparam

  param skip
caption = "Iters to skip"
default = 0
hint = "Iterations to skip."
  endparam
  param off
caption = "Hue cycle"
default = 0
min = 0
max = 7
hint = "This rotates the coloring order of the balls."
  endparam
}

；忘了，忘完了，贴出代码，为了纪念

作者: lnsz_dzg 时间: 2021-5-20 16:57

图片附件: 001.png (2021-5-20 16:57, 68.86 KB) / 下载次数 2555
http://develop.forums.netpad.net.cn/attachment.php?aid=28916&k=4b1e32dbce104857af7494340acbe8c0&t=1749930301&sid=bLX5X5

图片附件: 002.png (2021-5-20 16:57, 70.88 KB) / 下载次数 2610
http://develop.forums.netpad.net.cn/attachment.php?aid=28917&k=a8b8fb8967d89280a79dafd960b23593&t=1749930301&sid=bLX5X5

图片附件: 003.png (2021-5-20 16:57, 64.57 KB) / 下载次数 2544
http://develop.forums.netpad.net.cn/attachment.php?aid=28918&k=cdbdeb433d66a7de1f70af4bb46dcd5d&t=1749930301&sid=bLX5X5

图片附件: 004.png (2021-5-20 16:57, 89.47 KB) / 下载次数 2558
http://develop.forums.netpad.net.cn/attachment.php?aid=28919&k=5b5ff33b5dea51ac6ebdad5554101668&t=1749930301&sid=bLX5X5

图片附件: 005.png (2021-5-20 16:57, 55.17 KB) / 下载次数 2556
http://develop.forums.netpad.net.cn/attachment.php?aid=28920&k=1ff7b41a6b3dd523c45b228778dc9048&t=1749930301&sid=bLX5X5

图片附件: 006.png (2021-5-20 16:57, 58.54 KB) / 下载次数 2623
http://develop.forums.netpad.net.cn/attachment.php?aid=28921&k=3694594ee625d815110732f94f7db1fb&t=1749930301&sid=bLX5X5

图片附件: 007.png (2021-5-20 16:57, 79.82 KB) / 下载次数 2630
http://develop.forums.netpad.net.cn/attachment.php?aid=28922&k=ad6d4143340532d6c0a41daca3080da9&t=1749930301&sid=bLX5X5

图片附件: 008.png (2021-5-20 16:57, 94.52 KB) / 下载次数 2665
http://develop.forums.netpad.net.cn/attachment.php?aid=28923&k=e3eec970f3b0f08873e35ce548a31831&t=1749930301&sid=bLX5X5

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