Newton分形着色再探讨

我们以下图为例来逐步实现

这是这个分形的UF代码
Newton8 {
;LNSZDZG于2016年3月15日 修改原分形代码
;原代码的初值是 c = cos(pi/4) + flip(sin(pi/4))
; 迭代函数是 z = z - (c^2*z^4-(c^4+1)*z^2+c^2)/(4*c^2*z^3-2*(c^4+1)*z)
init:
complex zold = (10,0)
complex R = @R
z = #pixel
z = z + 1/z
loop:
zold = z
z = z - (z^R +1)/(R*z^(R-1))
bailout:
|z - zold| > @bailout
default:
title = "Newton8"
maxiter = 1000
periodicity = 0
center = (0,0)
magn = 1
param R
caption = "z的指数"
default = 4.0
endparam
param bailout
caption = "Bailout"
default = 0.00001
endparam
switch:
type = "Newton8 "
power = @power
bailout = @bailout
}

先用MC将这个分形“做出来”
这是程序和效果图

1、平滑eT环
将原参数s=0.025*iter修改为
s=0.025*(iter - ln(|log(|z-zold|)|)*ln(|log(|bailout)|))
其余参数不变，eT环比较平滑了

2、分区域着色
增加参数AR = arg(z)/2π，修改
R[j,k=AR*250,
B[j,k=(1-AR)*250
其余参数不变

程序是这样的

将着色参数修改为

并将图片属性修改为HSV，效果就是这样

这是灰度的效果

用点“旁门左道”——注意下面的
R , B*0.7，G

忘了说了，上图的属性是HLS！！！

收工了！欢迎朋友们一起谈论。