changxde

请教一个问题，看你们多次提到em，我不清楚em特征值是什么，请指教。

changxde

这个图的参数是多少，我用UF试了很长时间也没找到满意的参数。

下载 (10.56 KB)

2010-6-27 20:21

changxde

可是中间那个平行四边形就是出不来

changxde

下载 (17.25 KB)

2010-6-29 21:15

学习做一个IFS

changxde

61# 分形几何

我是按解方程求出12个根进行迭代得到的那个图，为什么会有缺口呢

changxde

迭代10万次才到这效果

下载 (56.72 KB)

2010-7-5 17:30

http://www.jamesh.id.au/fractals/orbit/threeply.html

changxde

给定的值-55，-1，-42

changxde

109# math
这是一个IFS分形,用到茂比乌斯变换,帮助中的e文看不懂,另外还有梅老师推荐的那本书http://www.inrm3d.cn/viewthread.php?tid=1171&extra=page%3D1%26amp%3Bfilter%3Ddigest中也有这个内容.可惜也是e文.

changxde

global:
  
  Mobius a0
  Mobius b0
  root = IIf(AlternateRoot, 0, 1)
  
  switch (Recipe) {
    case Recipes.Grandma
      Mobius.GrandmaRecipe(TraceA, TraceB, root, a0, b0)
      specialWordP = GrandmaP
      specialWordQ = GrandmaQ

    case Recipes.Jergensen
      Mobius.JergensenRecipe(TraceA, TraceB, root, a0, b0)
      specialWordP = GrandmaP
      specialWordQ = GrandmaQ

    case Recipes.Riley
      Mobius.RileyRecipe(RileyParameter, a0, b0)
      specialWordP = RileyP
      specialWordQ = RileyQ

    case Recipes.Maskit
      Mobius.MaskitRecipe(MaskitParameter, a0, b0)
      specialWordP = MaskitP
      specialWordQ = MaskitQ
  }
  Mobius R0
  Mobius R1
  '
  ' Transform/Conjugate the base generators a0 and b0.
  '
  if (ApplyUnitCircleGroupTransform) {
    R0 = Mobius.UnitCircleGroup(DegreeToRadian(UCGAngle), UCGPoint)
    R1 = Mobius.Inverse(R0)
    a0 = Mobius.Multiply3(R0, a0, R1)
    b0 = Mobius.Multiply3(R0, b0, R1)
  }
  if (ConjugateBy <> MobiusSpecialMatrixTypes.Identity) {
    R0 = Mobius.SpecialMatrix(ConjugateBy)
    R1 = Mobius.Inverse(R0)
    a0 = Mobius.Multiply3(R0, a0, R1)
    b0 = Mobius.Multiply3(R0, b0, R1)
  }
  Mobius.Normalize(a0)
  Mobius.Normalize(b0)
  Mobius A1 = Mobius.Inverse(a0)
  Mobius B1 = Mobius.Inverse(b0)
  '
  ' Define the commutators.
  '
  Mobius abAB = Mobius.Multiply4(a0, b0, A1, B1)
  Mobius baBA = Mobius.Inverse(abAB)
  '
  ' Fill w[] with the base set and the commutators.
  ' Include the special word if required.
  '
  if (specialWordP <> 0 && specialWordQ <> 0) {
    Mobius sw0 = Mobius.PQWord(specialWordP, specialWordQ, a0, B1)
    Mobius sw1 = Mobius.Inverse(sw0)
    Mobius w[] = a0, b0, A1, B1, abAB, baBA, sw0, sw1
  } else {
    Mobius w[] = a0, b0, A1, B1, abAB, baBA
  }
  const Complex base = Array.Dim1(w[])
  const Complex count = base * Power
  const Mobius generator[count]
  const Complex indexLUT[10000]
  Complex p[base]
  
  for (i = 0, i < base, i += 1) {
    p = IIf(i < 4, 1, CommutatorWeight)
  }
  Math.NormalizeWeights(p[], base)
  Math.GenerateIndexLookupTable(p[], base, indexLUT[])

  for (i = 0, i < Power, i += 1) {
    Array.Copy(w[], 0, generator[], i*base, base)
   
    for (j = 0, j < base, j += 1) {
      w[j] = Mobius.Multiply(w[j], generator[j])
    }
  }
  
initialize:
  
  z = Mobius.Sink(abAB)
  
iterate:
  
  attractorIndex = Math.GenerateIndex(indexLUT[])

  if (Random.Number() < Weight) {
    index = attractorIndex
  } else {
    index = attractorIndex * Random.Integer(Power)
  }
  z = Mobius.TransformPoint(generator[index], z)
  
properties:.....

在http://www.fractalsciencekit.com/program/orbeqn.htm中的帮助中有些没有表达式.

changxde

改变FSK中的参数得到下列图形