跳轉到內容

分形/島嶼尾流

來自華夏公益教科書,開放的書籍,開放的世界

p/q 尾流中最大的島嶼

  • 對於 p/q <= 1/3
    • 位於第一分支(從左到右閱讀)
    • 週期 = q+1
    • 外部射線 ?
  • 對於 p/q > 1/3
    • 位於最後的分支(從左到右閱讀)
    • 週期 = q+2

尾流 1/2

[編輯 | 編輯原始碼]

週期 3 島嶼是尾流 1/2 中最大的島嶼


找到週期 3 分量(島嶼)的一些子球泡的角度,這些子球泡位於主天線,具有外部射線 (3/7,4/7)

平面描述 :[1]

-1.76733 +0.00002 i @ 0.05

可以使用 Wolf Jung 的 Mandel 程式 進行檢查

The angle  3/7  or  p011 has  preperiod = 0  and  period = 3.
The conjugate angle is  4/7  or  p100 .
The kneading sequence is  AB*  and the internal address is  1-2-3 .
The corresponding parameter rays are landing at the root of a primitive component of period 3.

另請參閱


週期 4 島嶼

  • 中心 X -0.15710375803
  • 中心 Y +1.03258348530
  • 畫素步長 +0.0000375

尾流 1/7

[編輯 | 編輯原始碼]

尾流 12/25

  • 根 c = -0.738203140939397 +0.124839088573366 i
  • 主心形的主要尾流 12/25 由角度為引數射線包圍
    • 11184809/33554431 或 p0101010101010101010101001 和
    • 11184810/33554431 或 p0101010101010101010101010。
  • 衛星分量 c = -0.739829393511579 +0.125072144080321 i 週期 = 25 的中心

尾流 12/25 中最大的島嶼

  • 心形
    • 中心 c = -0.744245042107463 +0.127908444364520 i
    • 週期 = 27
    • 尖點


Haskell 程式

size formula degree p a b = do
  -- z = x + i y = 0
  x <- real
  y <- real
  x .= 0
  y .= 0
  -- matrix L
  lxa <- real
  lxb <- real
  lya <- real
  lyb <- real
  -- L = identity
  lxa .= 1
  lxb .= 0
  lya .= 0
  lyb .= 1
  -- matrix B
  bxa <- real
  bxb <- real
  bya <- real
  byb <- real
  -- B = identity
  bxa .= 1
  bxb .= 0
  bya .= 0
  byb .= 1
  -- loop one period
  j <- int
  j .= 1
  while_ (j .< p) $ do
    -- allocate next z, l
    xn <- real
    yn <- real
    lxan <- real
    lxbn <- real
    lyan <- real
    lybn <- real
    -- calculate next z, l
    let R2 fx fy = formula (R2 a b) (R2 x y)
        -- calculate derivatives
        fdxa = ddx d x fx
        fdxb = ddx d y fx
        fdya = ddx d x fy
        fdyb = ddx d y fy
        d v u
              | u == x && v == x = lxa
              | u == x && v == y = lxb
              | u == y && v == x = lya
              | u == y && v == y = lyb
              | u == v = 1
              | otherwise = 0
    -- z = f(z, c)
    xn .= fx
    yn .= fy
    x .= xn
    y .= yn
    -- L = J_f(z, c)
    lxan .= fdxa
    lxbn .= fdxb
    lyan .= fdya
    lybn .= fdyb
    lxa .= lxan
    lxb .= lxbn
    lya .= lyan
    lyb .= lybn
    -- B = B + 1/L
    det <- real
    det .= lxa * lyb - lxb * lya
    bxa .= bxa + lyb / det
    bxb .= bxb - lxb / det
    bya .= bya - lya / det
    byb .= byb + lxa / det
    -- loop counter
    j .= j + 1
  -- l = sqrt (abs (det L))
  l <- real
  l .= sqrt (abs (lxa * lyb - lxb * lya))
  -- beta = sqrt (abs (det B))
  beta <- real
  beta .= sqrt (abs (bxa * byb - bxb * bya))
  -- compute l^d b
  d <- float
  d .= degree / (degree - 1)
  llb <- real
  llb .= l ** d * beta
  return_ (1 / llb)

參考文獻

[編輯 | 編輯原始碼]
  1. Robert P. Munafo 於 2008 年 2 月 28 日的 R2F(1/2B1)S。
華夏公益教科書