xiaongxp

190# 柳烟
迭代次数的奇偶不同，分别对应着割法与补法Koch雪花

下载 (9.86 KB)

2015-1-7 19:02

Koch雪花[两个相似变换·外迭代].gsp (11.17 KB)

下载次数: 1856

2015-1-7 19:02

柳烟

向老师作得相当不错，最近向老师文件，对画板扫描IFS分形作出了重要贡献，容我慢慢消化。

柳烟

下载 (141.67 KB)

2015-1-10 08:54

UF中独缺这谢氏地毯，今用前面外迭代法，造一个。
代码如下：
Sierpinski正方形地毯** {

init:
  z = #pixel
  bool bail = false
loop:
    x = real(z)
    y = imag(z)
    if abs(x) <0.4&&abs(y) <0.4
      bail = true
    elseif y>0.4&& abs(x) <0.4
      z = 3*z  - 1i*2.4
   elseif y<-0.4&& abs(x) <0.4
      z = 3*z  + 1i*2.4
  elseif x>0.4&& abs(y) <0.4
      z = 3*z  -2.4
  elseif x<-0.4&& abs(y) <0.4
      z = 3*z  +2.4
  elseif x<-0.4&& y>0.4
      z = 3*z  +2.4-1i*2.4
elseif x>0.4&& y>0.4
      z = 3*z -2.4-1i*2.4
elseif x<-0.4&& y<-0.4
      z = 3*z +2.4+1i*2.4
elseif x>0.4&& y<-0.4
      z = 3*z -2.4+1i*2.4

    endif
bailout:
  bail == false
default:
  title = "Sierpinski正方形地毯"
  helpfile = "sam-help/various.htm"
  helptopic = "sierp2"
  center = (0,0.0)
}

柳烟

UF中的勾股树：

下载 (31.06 KB)

2015-1-9 21:56

附我编写的代码：
勾股数 {

init:
  z = #pixel
  float x = real(z)
  float y = imag(z)
  float a =1*pi/2
  float b=1*pi/2-1/2*a
  bool bail = false

loop:
    x = real(z)
    y = imag(z)
    if x>-0.3&& x<0.3&& y>-0.6&&y<0
      bail = true
    elseif x> cos(a)
      z = z-0.3
      z=z*exp(1i*b)+0.3
    z=1/cos(b)*(z-0.3)+0.3-1i*0.6
    elseif x<cos(a)
      z = z+0.3
      z=z*exp(-1i*(pi*0.5-b))-0.3
     z=1/sin(b)*(z+0.3)-0.3-1i*0.6

    endif
bailout:
  bail == false
default:
  title = "勾股树"
  helpfile = "sam-help/various.htm"
  helptopic = "sierp2"
  center = (0,0.5)
}

柳烟

下载 (39.39 KB)

2015-1-9 21:57

UF的产品，有些地方很毛躁，不知能不能想办法弄光滑些，如果有UF玩家知，烦告知，不胜感激。

xiaongxp

用画板作点陷阱效果还可以

下载 (28.04 KB)

2015-1-10 21:03

下载 (25.78 KB)

2015-1-9 22:51

xiaongxp

下载 (14.72 KB)

2015-1-10 23:15

下载 (13.4 KB)

2015-1-10 23:15

勾股树3[外迭代].gsp (11 KB)

下载次数: 2047

2015-1-10 23:15

柳烟

将UF科赫雪花开关项1代码提取出来，再取代码中的一段落，看看究竟有何几何意义，想了一天，有了成果，于是得到开关项1UF中科赫雪花补法的又一种外迭代法：
KochCurve {
init:
  z = #pixel
  x = real(z)
  y = imag(z)
  sq3 = sqrt(3)
  bail2 =false
loop:
    z = 3*z
    x = real(z)
    y = imag(z)
     if (y > 0) && (sq3*x - y + sq3 > 0) \
       && (sq3*x + y - sq3 < 0)
      bail2 = true
    endif
    z = z/3
    x = real(z)

    if x < -1/3
      z = 3*z + 2
    elseif x > 1/3
      z = 3*z - 2
    else
if x < 0
        z = z + 1/3
        z = z*exp(-1i*pi/3)
        z = 3*z - 1
      elseif x> 0
        z = z - 1/3
        z = z*exp(1i*pi/3)
        z = 3*z + 1
      endif
    endif

bailout:
bail2 == false
default:
  title = "Koch Curve1"
  helpfile = "sam-help/kochcurves.htm"
  helptopic = "kcurve"
  magn = 1.5
  center = (0.0002,0)
  maxiter = 50
}
此代码已经GSP化了，与前面向老师的外迭代法原理惊人一致。这抽取的代码只是科赫雪花三支中的一支。
UF科赫雪花开关项1中补法代码1支（外迭代法）.gsp (16.25 KB)

下载次数: 1984

2015-1-10 21:12

柳烟

上楼文件的扫图：

下载 (2.87 KB)

2015-1-11 09:21

代码意义研究草图

下载 (10.54 KB)

2015-1-10 21:00

柳烟

下载 (7.46 KB)

2015-1-11 09:19

科赫雪花开关1中1支的旋转变换（外迭代法).gsp (18.75 KB)

下载次数: 1909

2015-1-11 00:23

附文件草图。

下载 (21.63 KB)

2015-1-11 00:23