返回列表 回复 发帖
以上几张图片十分漂亮!
不知道在mathcad中扫描放大图的速度如何,我在mahematica中放大一个迭代次数只有500的图把电脑的内存一下就上升到90%,电脑就像酒喝醉了一样。
Untitled-1.jpg
换了好几个mathcad安装文件,都不能完成安装。
换了好几个mathcad安装文件,都不能完成安装。
xiaongxp 发表于 2014-4-29 20:12
出现什么样的提示?
将z的幅角加进来,可以修饰M集的边界

New.jpg
2014-7-9 09:54
将三角不等式着色法中的指数加大为
sum = sum + ((cabs(#z) - lowbound) / (az2+ac - lowbound))^@apower,
就有这种效果
New.jpg
2014-7-9 09:55
在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

New.jpg
2014-7-9 09:56

New.jpg
2014-7-9 09:56

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

New.jpg
2014-7-9 09:57



New.jpg
2014-7-9 09:58
上图是利用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代码只能用图片,太长了上传很不便啊
返回列表