xuefeiyang

色彩的分离：

下载 (85.31 KB)

2010-10-31 16:57

下载 (10.22 KB)

2010-10-31 17:00

xuefeiyang

这个课题大家动动脑筋，众人划桨开大船！

xuefeiyang

在画板中实现这段代码：
init:
  pixel2=-abs(real(pixel))+flip(imag(pixel))
  d=log(1/pixel2-conj(cos(1/pixel2)))
  c=pixel2^3-conj(.1/pixel2-flip(.01/pixel2))
  z=abs(pixel2^12)-flip(.01/pixel2)-conj(.1/pixel2)+.3
  zmod = (-.7456,-.132) + flip(|d/80*pixel2|)
  cmod = |c/7.5| + zmod
  c1=1.5*z^1.2,c2=2.25*conj(z)
  c3=3.375*flip(z),c4=5.0625*flip(z)
  l1=round(real(p1)),l2=round(imag(p1))
  l3=round(real(p2)),l4=round(imag(p2))
  float bailout=16,int iter=0
loop:
  IF (iter==l1)
    cmod = |c1/7.5| + zmod, z = 0
  ELSEIF (iter==l2)
    cmod = |c2/7.5| + zmod, z = 0
  ELSEIF (iter==l3)
    cmod = |c3/7.5| + zmod, z = 0
  ELSEIF (iter==l4)
    cmod = |c4/7.5| + zmod, z = 0
  ENDIF
  z = z^2 + cmod
  iter=iter+1
bailout:
  |z|<=bailout
default:
  title = "Carr 3202"
  periodicity = 0
  maxiter = 700
  magn = 1.4
  center = (0,0)
  method = multipass
  param p1
    caption = "Iter Limits 1,2"
    default = (200,300)
    hint = "0< Re(P1) < Im(P1) < Re(P2) < Im(P2) <maxiter \
            for best results"
  endparam
  param p2
    caption = "Iter Limits 3,4"
    default = (400,550)
    hint = "0< Re(P1) < Im(P1) < Re(P2) < Im(P2) <maxiter \
            for best results"
  endparam
}

xuefeiyang

IF (iter==l1)
    cmod = |c1/7.5| + zmod, z = 0
  ELSEIF (iter==l2)
    cmod = |c2/7.5| + zmod, z = 0
  ELSEIF (iter==l3)
    cmod = |c3/7.5| + zmod, z = 0
  ELSEIF (iter==l4)
    cmod = |c4/7.5| + zmod, z = 0

这段代码什么意思？“＝＝”表示什么意思？

xuefeiyang

复变分形的计算问题很难解决，暂时先放下，奉行“舍得”现实做法。到目前为止，分形绘制的扫描框架与着色方法基本可以了。可以玩一把分形的结构了。就先从UF上的carr904开始作起。
Carr-904 {
  z=((1-pixel)/3)^0.5, c=(1-1/pixel)^0.5:
  z=z*z+c
  |real(z)| <=100
  ;SOURCE: carr.frm
}

xuefeiyang

有哪位朋友知道函数ident的表达式是什么？

xuefeiyang

谢谢二位。UF上有不少术语还不懂。在哪里可以找到UF上的函数的意思呢？

xuefeiyang

255# 柳烟


这是一个M集与J集的混合集，左边是M集，右边是J集，肯定不会有经典的小M集出来了。

xuefeiyang

引入参数K,m。

下载 (20.66 KB)

2010-11-14 10:29

xuefeiyang

Carr1998 (XAXIS) {; Updated for UF2 by Erik Reckase, March 2000
                  ; Modified Sylvie Gallet frm. [101324,3444],1996
init:
  z=c=pixel, c1=1.5*z, c2=2.25*z, c3=3.375*z, c4=5.0625*z
  l1=round(real(p1)), l2=round(imag(p1))
  l3=round(real(p2)), l4=round(imag(p2))
  float bailout=16, int iter=0, pp2 = pixel/imag(p2)
loop:
  IF (iter==l1)
    z = 0, c = c1
  ELSEIF (iter==l2)
    z = 0, c = c2
  ELSEIF (iter==l3)
    z = 0, c = c3
  ELSEIF (iter==l4)
    z = 0, c = c4
  ENDIF
  c = c + pp2
  z=z*z+c
  iter=iter+1
bailout:
  |real(z)| <= bailout
default:
  title = "Carr 1998"
  periodicity = 0
  maxiter = 500
  magn = 1
  center = (0,0)
  method = multipass
  param p1
    caption = "Iter Limits 1,2"
    default = (100,150)
    hint = "0< Re(P1) < Im(P1) < Re(P2) < Im(P2) <maxiter \
            for best results"
  endparam
  param p2
    caption = "Iter Limits 3,4"
    default = (200,250)
    hint = "0< Re(P1) < Im(P1) < Re(P2) < Im(P2) <maxiter \
            for best results"
  endparam
}
作这个分形，不知怎么了，作出这样一个怪胎来！

下载 (11.57 KB)

2010-11-14 14:43

原图这样的：

下载 (49.23 KB)

2010-11-14 14:43